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

  1. Alexandre P Thiery
  2. Ailin Leticia Buzzi
  3. Eva Hamrud
  4. Chris Cheshire
  5. Nicholas M Luscombe
  6. James Briscoe
  7. Andrea Streit  Is a corresponding author
  1. King's College London, United Kingdom
  2. The Francis Crick Institute, United Kingdom

Abstract

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.

Data availability

10x single cell RNAseq was carried out in two batches and is available under two separate accession numbers (ArrayExpress: E-MTAB-10408 and E-MTAB-1144). Our NGS alignments and downstream analysis have been wrapped into custom Nextflow pipelines allowing for full reproducibility. For the code used in this analysis, including links to our Docker containers, see our GitHub repository at https://github.com/alexthiery/10x_neural_plate_border. Finally, we have developed a user friendly ShinyApp to allow public exploration of our single cell RNAseq data at https://shiny.crick.ac.uk/thiery_neural_plate_border/.

The following data sets were generated
The following previously published data sets were used

Article and author information

Author details

  1. Alexandre P Thiery

    Centre for Craniofacial and Regenerative Biology, King's College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  2. Ailin Leticia Buzzi

    Centre for Craniofacial and Regenerative Biology, King's College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  3. Eva Hamrud

    Centre for Craniofacial and Regenerative Biology, King's College London, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  4. Chris Cheshire

    Bioinformatics and Computational Biology Laboratory, The Francis Crick Institute, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  5. Nicholas M Luscombe

    Bioinformatics and Computational Biology Laboratory, The Francis Crick Institute, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
  6. James Briscoe

    Bioinformatics and Computational Biology Laboratory, The Francis Crick Institute, London, United Kingdom
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0002-1020-5240
  7. Andrea Streit

    Centre for Craniofacial and Regenerative Biology, King's College London, London, United Kingdom
    For correspondence
    andrea.streit@kcl.ac.uk
    Competing interests
    The authors declare that no competing interests exist.
    ORCID icon "This ORCID iD identifies the author of this article:" 0000-0001-7664-7917

Funding

Biotechnology and Biological Sciences Research Council (BB/S005536/1)

  • Alexandre P Thiery
  • Ailin Leticia Buzzi
  • Andrea Streit

Biotechnology and Biological Sciences Research Council (BB/R006342/1)

  • Alexandre P Thiery
  • Ailin Leticia Buzzi
  • Andrea Streit

Wellcome Trust (108874/B/15/Z)

  • Nicholas M Luscombe

Wellcome Trust (FC001051)

  • Chris Cheshire
  • Nicholas M Luscombe
  • James Briscoe

Cancer Research UK (FC001051)

  • Chris Cheshire
  • Nicholas M Luscombe
  • James Briscoe

Medical Research Council (FC001051)

  • Chris Cheshire
  • Nicholas M Luscombe
  • James Briscoe

Wellcome Trust (108874/Z/15/Z)

  • Eva Hamrud

The funders had no role in study design, data collection and interpretation, or the decision to submit the work for publication.

Reviewing Editor

  1. Richard A Schneider, University of California, San Francisco, United States

Version history

  1. Preprint posted: February 16, 2022 (view preprint)
  2. Received: August 15, 2022
  3. Accepted: August 1, 2023
  4. Accepted Manuscript published: August 2, 2023 (version 1)
  5. Version of Record published: August 14, 2023 (version 2)

Copyright

© 2023, Thiery et al.

This article is distributed under the terms of the Creative Commons Attribution License permitting unrestricted use and redistribution provided that the original author and source are credited.

Metrics

  • 1,234
    views
  • 261
    downloads
  • 8
    citations

Views, downloads and citations are aggregated across all versions of this paper published by eLife.

Download links

A two-part list of links to download the article, or parts of the article, in various formats.

Downloads (link to download the article as PDF)

Open citations (links to open the citations from this article in various online reference manager services)

Cite this article (links to download the citations from this article in formats compatible with various reference manager tools)

  1. Alexandre P Thiery
  2. Ailin Leticia Buzzi
  3. Eva Hamrud
  4. Chris Cheshire
  5. Nicholas M Luscombe
  6. James Briscoe
  7. Andrea Streit
(2023)
scRNA-sequencing in chick suggests a probabilistic model for cell fate allocation at the neural plate border
eLife 12:e82717.
https://doi.org/10.7554/eLife.82717

Share this article

https://doi.org/10.7554/eLife.82717

Further reading

    1. Computational and Systems Biology
    2. Developmental Biology
    Arya Y Nakhe, Prasanna K Dadi ... David A Jacobson
    Research Article

    The gain-of-function mutation in the TALK-1 K+ channel (p.L114P) is associated with maturity-onset diabetes of the young (MODY). TALK-1 is a key regulator of β-cell electrical activity and glucose-stimulated insulin secretion. The KCNK16 gene encoding TALK-1 is the most abundant and β-cell-restricted K+ channel transcript. To investigate the impact of KCNK16 L114P on glucose homeostasis and confirm its association with MODY, a mouse model containing the Kcnk16 L114P mutation was generated. Heterozygous and homozygous Kcnk16 L114P mice exhibit increased neonatal lethality in the C57BL/6J and the CD-1 (ICR) genetic background, respectively. Lethality is likely a result of severe hyperglycemia observed in the homozygous Kcnk16 L114P neonates due to lack of glucose-stimulated insulin secretion and can be reduced with insulin treatment. Kcnk16 L114P increased whole-cell β-cell K+ currents resulting in blunted glucose-stimulated Ca2+ entry and loss of glucose-induced Ca2+ oscillations. Thus, adult Kcnk16 L114P mice have reduced glucose-stimulated insulin secretion and plasma insulin levels, which significantly impairs glucose homeostasis. Taken together, this study shows that the MODY-associated Kcnk16 L114P mutation disrupts glucose homeostasis in adult mice resembling a MODY phenotype and causes neonatal lethality by inhibiting islet insulin secretion during development. These data suggest that TALK-1 is an islet-restricted target for the treatment for diabetes.

    1. Developmental Biology
    2. Structural Biology and Molecular Biophysics
    Arne Elofsson, Ling Han ... Luca Jovine
    Research Article

    A crucial event in sexual reproduction is when haploid sperm and egg fuse to form a new diploid organism at fertilization. In mammals, direct interaction between egg JUNO and sperm IZUMO1 mediates gamete membrane adhesion, yet their role in fusion remains enigmatic. We used AlphaFold to predict the structure of other extracellular proteins essential for fertilization to determine if they could form a complex that may mediate fusion. We first identified TMEM81, whose gene is expressed by mouse and human spermatids, as a protein having structural homologies with both IZUMO1 and another sperm molecule essential for gamete fusion, SPACA6. Using a set of proteins known to be important for fertilization and TMEM81, we then systematically searched for predicted binary interactions using an unguided approach and identified a pentameric complex involving sperm IZUMO1, SPACA6, TMEM81 and egg JUNO, CD9. This complex is structurally consistent with both the expected topology on opposing gamete membranes and the location of predicted N-glycans not modeled by AlphaFold-Multimer, suggesting that its components could organize into a synapse-like assembly at the point of fusion. Finally, the structural modeling approach described here could be more generally useful to gain insights into transient protein complexes difficult to detect experimentally.