Reconstruction of the global neural crest gene regulatory network in vivo

Ruth M Williams, Ivan Candido-Ferreira, Emmanouela Repapi, Daria Gavriouchkina, Upeka Senanayake, Jelena Telenius, Stephen Taylor, Jim Hughes, Tatjana Sauka-Spengler

Preprint posted on December 30, 2018

How to build a vertebrate neural crest: reconstructing neural crest gene regulatory networks in vivo

Selected by Hannah Brunsdon



Gene regulatory networks (GRNs) are formed of molecular regulators encoded in DNA, which process information from transcription factors and signalling molecules to regulate many growth and developmental processes. Much of what we know about vertebrate GRNs is inferred from manipulation of candidate genes, and GRNs of simpler organisms – the sheer complexity of vertebrate gene regulation makes it extremely difficult to gain an unbiased view of GRN interactions in vertebrate development. This has not deterred the authors of a recent preprint from the Sauka-Spengler lab, however. Their work shows the output of an enormous effort in reconstructing the chick neural crest (NC) GRN using both transcriptomic and epigenomic profiling, thus creating the first high-resolution, genome-wide reconstruction of a vertebrate GRN in vivo.


Key findings

To gain an overview of the early cranial NC transcriptional landscape, the authors isolated NC cells from chick embryos at both pre-migratory stages (5-6ss) and early migratory stages (8-10ss) for RNA seq. The data clustered into 13 gene expression modules that are associated with the changing nature of the NC as cells undergo EMT and delaminate from the neural tube. To investigate the regulation of these processes at the epigenomic level, the authors performed ATAC-seq in conjunction with Next-Generation Capture-C on stage-matched samples. This revealed differences in chromatin accessibility upstream of core NC genes, revealing the presence of enhancer elements, as well as physical interactions between these and putative target genes.

Diving deeper into the ATAC-seq data, the authors performed k-means clustering analyses based on patterns of chromatin accessibility. Of particular interest was “k-cluster 3”, which contained elements only accessible in NC cells, and “k-cluster 9”-containing elements, enriched in both the NC and neighbouring neuroepithelial ectoderm. Cis-regulatory elements from both these clusters were then assigned to their nearest-neighbour genes, in order to build predictive networks of how different combinations of active enhancers cooperate in NC gene activation.

What I think makes this preprint so special, is that predictions made from their different datasets are tested and visualised in vivo. For example, the locus of Sox10 – a known master regulator of NC cell fate – was investigated by ChIP-seq to find novel putative enhancers. The activity of these enhancers was visualised in vivo by electroporating chick embryos with fluorescently-tagged enhancer reporter constructs. Individual Sox10 enhancer activity was spatially and temporally diverse, showing complex cis-regulation of its expression, as seen below:


From preprint Figure 4. In vivo activity of novel Sox10 enhancers at the 7 and 8 somite stage.


Taking this a step further, CRISPR-mediated genome editing was employed to test the functionality of these Sox10 enhancers. Guide RNAs targeting each enhancer were electroporated into chick embryos alongside the transcriptional repressor dCas9-Krab, thus selectively disrupting their activity. Subsequent RT-qPCR analysis on isolated neural tubes showed that disruption of single enhancers did not affect Sox10 expression for long, but simultaneous repression of certain enhancer combinations caused long-term repression.

So, does this heterogeneity of cis-regulatory elements influence NC lineage decisions? GO term enrichment analysis was performed on genes associated to k-cluster 3 and 9 elements. This confirmed that open k-cluster 3 elements drive NC development and differentiation programmes, whereas open k-cluster 9 elements prime cells for more neural fates. To investigate if these differences produce meaningful transcriptomic changes within individual NC cells, single cell RNA-seq was performed at the point of NC delamination. PCA analysis of the data revealed three gene clusters. Sc-cluster 1 and 3 strongly overlapped with k-cluster 3 genes and were characteristic of the canonical NC, whereas sc-cluster 2 was more characteristic of NC-derived neural progenitors, and strongly correlated with k-cluster 9 genes. Therefore, well before commitment of cells to distinct NC-derived lineages, there is already considerable heterogeneity in regulatory elements, which must be important for establishing diverse NC progenitor lineages.

In the final part of their paper, the authors identified transcription factor co-binding sites within their integrated transcriptomic and epigenomic data, thus building a comprehensive, genome-wide reconstruction of k-cluster 3 and k-cluster 9-dependent GRNs, and identifying the hierarchies of transcription factors that comprised each GRN’s minimal regulatory cores. Through CRISPR-Cas9-mediated disruption of these minimal cores and subsequent RNA-seq, the authors confirmed that the balance of k-cluster 3 and 9 GRNs, with likely cross-repressive interactions, regulates the balance between NC-derived neural and mesenchymal lineages.

If it wasn’t already clear, this research has generated a huge amount of data, tying together transcriptional, epigenetic, single-cell and bulk RNA-seq experiments, as well as in vivo validation. Thankfully for neural crest enthusiasts unfamiliar with manipulating large datasets, the authors have assembled their data into a user-friendly Shiny App (R-based, available from GitHub). Below are screenshots of the App interface visualising the upstream regulatory elements and k-cluster-dependent transcriptional networks of TFAP2B:


From Figure 7 in the preprint. Images from the Shiny App showing the regulatory elements, and transcription factor networks for TFAP2B


Why I chose this preprint

This preprint represents an incredible amount of work and team effort, and I was amazed by the number of Next-Generation seq analyses that were performed throughout the study, but nevertheless effectively pulled together to build a cohesive story about the early NC GRN. Secondly, whereas other large scale transcriptomics/epigenomics studies sometimes conclude with their computational predictive models, this preprint matches their intricate dataset analyses with ambitious in vivo work too.  For people like myself, who are more comfortable with wet lab experiments than manipulating huge datasets, the provision of a friendly app interface to enable the research community to play with the data themselves is a very valuable resource.


Questions for the authors

  • Given the complexity of the paper, this question will probably now sound really simplistic, but what is it about the chick as a model organism that made it ideal for building NC GRNs?
  • Do you have any plans to investigate other k-clusters and developmental stages to compare their GRNs to k-cluster 6 and 9? It would be interesting to know how k-cluster 6, which is involved in early establishment of the NC could differ to, and interact with, k-cluster 3 and 9 for example. Following on from that, by sampling at other stages would it be possible to follow NC GRNs from NC establishment to final lineage differentiation in this way?



Posted on: 28th January 2019

Read preprint (No Ratings Yet)

  • Have your say

    Your email address will not be published. Required fields are marked *

    This site uses Akismet to reduce spam. Learn how your comment data is processed.

    Sign up to customise the site to your preferences and to receive alerts

    Register here

    Also in the developmental biology category:

    The visual system of the genetically tractable crustacean Parhyale hawaiensis: diversification of eyes and visual circuits associated with low-resolution vision

    Ana Patricia Ramos, Ola Gustafsson, Nicolas Labert, et al.

    Selected by Alexa Sadier

    Distinct RhoGEFs activate apical and junctional actomyosin contractility under control of G proteins during epithelial morphogenesis

    Alain Garcia De Las Bayonas, Jean-Marc Philippe, Annemarie C. Lellouch, et al.

    Selected by Ivana Viktorinová


    The Hunchback temporal transcription factor determines motor neuron axon and dendrite targeting in Drosophila

    Austin Q Seroka, Chris Q Doe

    Selected by Abagael Lasseigne


    Preformed Chromatin Topology Assists Transcriptional Robustness of Shh during Limb Development

    Christina Paliou, Philine Guckelberger, Robert Schöpflin, et al.

    Selected by Rafael Galupa


    Molecular Logic of Spinocerebellar Tract Neuron Diversity and Connectivity

    Myungin Baek, Vilas Menon, Thomas Jessell, et al.

    Selected by Yen-Chung Chen

    In vivo glucose imaging in multiple model organisms with an engineered single-wavelength sensor

    Jacob P. Keller, Jonathan S. Marvin, Haluk Lacin, et al.

    Selected by Stephan Daetwyler


    Active behaviour during early development shapes glucocorticoid reactivity

    Luis A. Castillo-Ramírez, Soojin Ryu, Rodrigo J. De Marco

    Selected by Kathleen Gilmour

    A unicellular relative of animals generates an epithelium-like cell layer by actomyosin-dependent cellularization

    Omaya Dudin, Andrej Ondracka, Xavier Grau-Bové, et al.

    Selected by Paul Gerald L. Sanchez and Stefano Vianello


    Members of the Arabidopsis auxin receptor gene family are essential early in embryogenesis and have broadly overlapping functions

    Michael J Prigge, Nikita Kadakia, Kathleen Greenham, et al.

    Selected by Chandra Shekhar Misra


    The spindle assembly checkpoint functions during early development in non-chordate embryos

    Janet Chenevert, Marianne Roca, Lydia Besnardeau, et al.

    Selected by Maiko Kitaoka

    Over-activation of BMP signaling in neural crest cells precipitates heart outflow tract septation

    Jean Francois Darrigrand, Mariana Valente, Pauline Martinez, et al.

    Selected by Giulia Boezio

    Mutations in the Insulator Protein Suppressor of Hairy Wing Induce Genome Instability

    Shih-Jui Hsu, Emily C. Stow, James R. Simmons, et al.

    Selected by Maiko Kitaoka


    Multi-immersion open-top light-sheet microscope for high-throughput imaging of cleared tissues

    Adam K. Glaser, Nicholas P. Reder, Ye Chen, et al.

    Selected by Tim Fessenden


    Dynamic Erasure of Random X-Chromosome Inactivation during iPSC Reprogramming

    Adrian Janiszewski, Irene Talon, Juan Song, et al.

    Selected by Sergio Menchero

    MicroRNA-mediated control of developmental lymphangiogenesis

    Hyun Min Jung, Ciara Hu, Alexandra M Fister, et al.

    Selected by Rudra Nayan Das

    Endogenous CRISPR arrays for scalable whole organism lineage tracing

    James Cotterell, James Sharpe

    Selected by Irepan Salvador-Martinez

    Also in the genomics category:

    Accurate detection of m6A RNA modifications in native RNA sequences

    Huanle Liu, Oguzhan Begik, Morghan C Lucas, et al.

    Selected by Christian Bates


    Crowdfunded whole-genome sequencing of the celebrity cat Lil BUB identifies causal mutations for her osteopetrosis and polydactyly

    Mike Bridavsky, Heiner Kuhl, Arthur Woodruf, et al.

    Selected by Jesus Victorino, Gabriel Aughey


    Prospective, brain-wide labeling of neuronal subclasses with enhancer-driven AAVs

    Lucas T Graybuck, Adriana Sedeño-Cortés, Thuc Nghi Nguyen, et al.

    Selected by Jesus Victorino

    Self-reporting transposons enable simultaneous readout of gene expression and transcription factor binding in single cells

    Arnav Moudgil, Michael N Wilkinson, Xuhua Chen, et al.

    Selected by James Briscoe


    High-throughput functional analysis of lncRNA core promoters elucidates rules governing tissue-specificity

    Kaia Mattioli, Pieter-Jan Volders, Chiara Gerhardinger, et al.

    Selected by Clarice Hong

    Reconstruction of the global neural crest gene regulatory network in vivo

    Ruth M Williams, Ivan Candido-Ferreira, Emmanouela Repapi, et al.

    Selected by Hannah Brunsdon

    Charting a tissue from single-cell transcriptomes

    Mor Nitzan, Nikos Karaiskos, Nir Friedman, et al.

    Selected by Irepan Salvador-Martinez

    Single cell RNA-Seq reveals distinct stem cell populations that drive sensory hair cell regeneration in response to loss of Fgf and Notch signaling

    Mark E. Lush, Daniel C. Diaz, Nina Koenecke, et al.


    Distinct progenitor populations mediate regeneration in the zebrafish lateral line.

    Eric D Thomas, David Raible

    Selected by Rudra Nayan Das


    Maintenance of spatial gene expression by Polycomb-mediated repression after formation of a vertebrate body plan

    Julien Rougot, Naomi D Chrispijn, Marco Aben, et al.

    Selected by Yen-Chung Chen


    The embryonic transcriptome of Arabidopsis thaliana

    Falko Hofmann, Michael A Schon, Michael D Nodine

    Selected by Chandra Shekhar Misra


    Simultaneous multiplexed amplicon sequencing and transcriptome profiling in single cells

    Mridusmita Saikia, Philip Burnham, Sara H Keshavjee, et al.


    High-throughput targeted long-read single cell sequencing reveals the clonal and transcriptional landscape of lymphocytes

    Mandeep Singh, Ghamdan Al-Eryani, Shaun Carswell, et al.

    Selected by Samantha Seah

    The microbial basis of impaired wound healing: differential roles for pathogens, "bystanders", and strain-level diversification in clinical outcomes

    Lindsay Kalan, Jacquelyn S Meisel, Michael A Loesche, et al.

    Selected by Snehal Kadam

    Comparative analysis of droplet-based ultra-high-throughput single-cell RNA-seq systems

    Xiannian Zhang, Tianqi Li, Feng Liu, et al.

    Selected by Samantha Seah

    PUMILIO hyperactivity drives premature aging of Norad-deficient mice

    Florian Kopp, Mehmet Yalvac, Beibei Chen, et al.

    Selected by Carmen Adriaens

    LCM-seq reveals unique transcriptional adaption mechanisms of resistant neurons in spinal muscular atrophy

    Susanne Nichterwitz, Helena Storvall, Jik Nijssen, et al.


    Axon-seq decodes the motor axon transcriptome and its modulation in response to ALS

    Jik Nijssen, Julio Cesar Aguila Benitez, Rein Hoogstraaten, et al.

    Selected by Yen-Chung Chen

    LADL: Light-activated dynamic looping for endogenous gene expression control

    Mayuri Rege, Ji Hun Kim, Jacqueline Valeri, et al.

    Selected by Ivan Candido-Ferreira