Oscillatory differentiation dynamics fundamentally restricts the resolution of pseudotime reconstruction algorithms.

The challenge to understand differentiation and cell lineages in development has resulted in many bioinformatics software tools, notably those working with gene expression data obtained via single-cell RNA sequencing obtained at snapshots in time. Reconstruction methods for trajectories often proceed by dimension reduction, data clustering and then computation of a tree graph in which edges indicate closely related clusters. Cell lineages can then be deduced by following paths through the tree. In the case of multi-potent cells undergoing differentiation, this trajectory reconstruction involves the reconstruction of multiple distinct lineages corresponding to commitment to each of a set of distinct fates. Recent work suggests that there may be cases in which the cell differentiation process involves trajectories that explore, in a dynamic and oscillatory fashion, propensity to differentiate into a number of possible cell fates before commitment finally occurs. Here, we show theoretically that the presence of such oscillations provides intrinsic constraints on the quality and resolution of the trajectory reconstruction process, even for idealized noise-free data. These constraints point to inherent common limitations of current methodologies and serve both to provide additional challenge in the development of software tools and also may help to understand features observed in recent experiments.

Myron Gordon Award Lecture 2023: Painting the neural crest: How studying pigment cells illuminates neural crest cell biology.

It has been 30 (!!) years since I began working on zebrafish pigment cells, as a postdoc in the laboratory of Prof. Christiane Nüsslein-Volhard. There, I participated in the first large-scale mutagenesis screen in zebrafish, focusing on pigment cell mutant phenotypes. The isolation of colourless, shady, parade and choker mutants allowed us (as a postdoc in Prof. Judith Eisen's laboratory, and then in my own laboratory at the University of Bath since 1997) to pursue my ambition to address long-standing problems in the neural crest field. Thus, we have studied how neural crest cells choose individual fates, resulting in our recent proposal of a new, and potentially unifying, model which we call Cyclical Fate Restriction, as well as addressing how pigment cell patterns are generated. A key feature of our work in the last 10 years has been the use of mathematical modelling approaches to clarify our biological models and to refine our interpretations. None of this would have been possible without a hugely talented group of laboratory members and other collaborators from around the world-it has been, and I am sure will continue to be, a pleasure and privilege to work with you all!

A gene regulatory network combining Pax3/7, Sox10 and Mitf generates diverse pigment cell types in medaka and zebrafish.

Neural crest cells generate numerous derivatives, including pigment cells, and are a model for studying how fate specification from multipotent progenitors is controlled. In mammals, the core gene regulatory network for melanocytes (their only pigment cell type) contains three transcription factors, Sox10, Pax3 and Mitf, with the latter considered a master regulator of melanocyte development. In teleosts, which have three to four pigment cell types (melanophores, iridophores and xanthophores, plus leucophores e.g. in medaka), gene regulatory networks governing fate specification are poorly understood, although Mitf function is considered conserved. Here, we show that the regulatory relationships between Sox10, Pax3 and Mitf are conserved in zebrafish, but the role for Mitf is more complex than previously emphasized, affecting xanthophore development too. Similarly, medaka Mitf is necessary for melanophore, xanthophore and leucophore formation. Furthermore, expression patterns and mutant phenotypes of pax3 and pax7 suggest that Pax3 and Pax7 act sequentially, activating mitf expression. Pax7 modulates Mitf function, driving co-expressing cells to differentiate as xanthophores and leucophores rather than melanophores. We propose that pigment cell fate specification should be considered to result from the combinatorial activity of Mitf with other transcription factors.

Zebrafish pigment cells develop directly from persistent highly multipotent progenitors.

Neural crest cells are highly multipotent stem cells, but it remains unclear how their fate restriction to specific fates occurs. The direct fate restriction model hypothesises that migrating cells maintain full multipotency, whilst progressive fate restriction envisages fully multipotent cells transitioning to partially-restricted intermediates before committing to individual fates. Using zebrafish pigment cell development as a model, we show applying NanoString hybridization single cell transcriptional profiling and RNAscope in situ hybridization that neural crest cells retain broad multipotency throughout migration and even in post-migratory cells in vivo, with no evidence for partially-restricted intermediates. We find that leukocyte tyrosine kinase early expression marks a multipotent stage, with signalling driving iridophore differentiation through repression of fate-specific transcription factors for other fates. We reconcile the direct and progressive fate restriction models by proposing that pigment cell development occurs directly, but dynamically, from a highly multipotent state, consistent with our recently-proposed Cyclical Fate Restriction model.

New advances in CRISPR/Cas-mediated precise gene-editing techniques.

Over the past decade, CRISPR/Cas-based gene editing has become a powerful tool for generating mutations in a variety of model organisms, from Escherichia coli to zebrafish, rodents and large mammals. CRISPR/Cas-based gene editing effectively generates insertions or deletions (indels), which allow for rapid gene disruption. However, a large proportion of human genetic diseases are caused by single-base-pair substitutions, which result in more subtle alterations to protein function, and which require more complex and precise editing to recreate in model systems. Precise genome editing (PGE) methods, however, typically have efficiencies of less than a tenth of those that generate less-specific indels, and so there has been a great deal of effort to improve PGE efficiency. Such optimisations include optimal guide RNA and mutation-bearing donor DNA template design, modulation of DNA repair pathways that underpin how edits result from Cas-induced cuts, and the development of Cas9 fusion proteins that introduce edits via alternative mechanisms. In this Review, we provide an overview of the recent progress in optimising PGE methods and their potential for generating models of human genetic disease.

Trunk Neural Crest Migratory Position and Asymmetric Division Predict Terminal Differentiation.

The generation of complex structures during embryogenesis requires the controlled migration and differentiation of cells from distant origins. How these processes are coordinated and impact each other to form functional structures is not fully understood. Neural crest cells migrate extensively giving rise to many cell types. In the trunk, neural crest cells migrate collectively forming chains comprised of cells with distinct migratory identities: one leader cell at the front of the group directs migration, while followers track the leader forming the body of the chain. Herein we analysed the relationship between trunk neural crest migratory identity and terminal differentiation. We found that trunk neural crest migration and fate allocation is coherent. Leader cells that initiate movement give rise to the most distal derivativities. Interestingly, the asymmetric division of leaders separates migratory identity and fate. The distal daughter cell retains the leader identity and clonally forms the Sympathetic Ganglia. The proximal sibling migrates as a follower and gives rise to Schwann cells. The sympathetic neuron transcription factor phox2bb is strongly expressed by leaders from early stages of migration, suggesting that specification and migration occur concomitantly and in coordination. Followers divide symmetrically and their fate correlates with their position in the chain.

Notch controls the cell cycle to define leader versus follower identities during collective cell migration.

Coordination of cell proliferation and migration is fundamental for life, and its dysregulation has catastrophic consequences, such as cancer. How cell cycle progression affects migration, and vice versa, remains largely unknown. We address these questions by combining in silico modelling and in vivo experimentation in the zebrafish trunk neural crest (TNC). TNC migrate collectively, forming chains with a leader cell directing the movement of trailing followers. We show that the acquisition of migratory identity is autonomously controlled by Notch signalling in TNC. High Notch activity defines leaders, while low Notch determines followers. Moreover, cell cycle progression is required for TNC migration and is regulated by Notch. Cells with low Notch activity stay longer in G1 and become followers, while leaders with high Notch activity quickly undergo G1/S transition and remain in S-phase longer. In conclusion, TNC migratory identities are defined through the interaction of Notch signalling and cell cycle progression.

Review: The Role of Wnt/ß-Catenin Signalling in Neural Crest Development in Zebrafish.

The neural crest (NC) is a multipotent cell population in vertebrate embryos with extraordinary migratory capacity. The NC is crucial for vertebrate development and forms a myriad of cell derivatives throughout the body, including pigment cells, neuronal cells of the peripheral nervous system, cardiomyocytes and skeletogenic cells in craniofacial tissue. NC induction occurs at the end of gastrulation when the multipotent population of NC progenitors emerges in the ectodermal germ layer in the neural plate border region. In the process of NC fate specification, fate-specific markers are expressed in multipotent progenitors, which subsequently adopt a specific fate. Thus, NC cells delaminate from the neural plate border and migrate extensively throughout the embryo until they differentiate into various cell derivatives. Multiple signalling pathways regulate the processes of NC induction and specification. This review explores the ongoing role of the Wnt/ß-catenin signalling pathway during NC development, focusing on research undertaken in the Teleost model organism, zebrafish (Danio rerio). We discuss the function of the Wnt/ß-catenin signalling pathway in inducing the NC within the neural plate border and the specification of melanocytes from the NC. The current understanding of NC development suggests a continual role of Wnt/ß-catenin signalling in activating and maintaining the gene regulatory network during NC induction and pigment cell specification. We relate this to emerging models and hypotheses on NC fate restriction. Finally, we highlight the ongoing challenges facing NC research, current gaps in knowledge, and this field's potential future directions.

Cell Fate Decisions in the Neural Crest, from Pigment Cell to Neural Development.

The neural crest shows an astonishing multipotency, generating multiple neural derivatives, but also pigment cells, skeletogenic and other cell types. The question of how this process is controlled has been the subject of an ongoing debate for more than 35 years. Based upon new observations of zebrafish pigment cell development, we have recently proposed a novel, dynamic model that we believe goes some way to resolving the controversy. Here, we will firstly summarize the traditional models and the conflicts between them, before outlining our novel model. We will also examine our recent dynamic modelling studies, looking at how these reveal behaviors compatible with the biology proposed. We will then outline some of the implications of our model, looking at how it might modify our views of the processes of fate specification, differentiation, and commitment.

Contribution of sox9b to pigment cell formation in medaka fish.

SoxE-type transcription factors, Sox10 and Sox9, are key regulators of the development of neural crest cells. Sox10 specifies pigment cell, glial, and neuronal lineages, whereas Sox9 is reportedly closely associated with skeletogenic lineages in the head, but its involvement in pigment cell formation has not been investigated genetically. Thus, it is not fully understood whether or how distinctly these genes as well as their paralogs in teleosts are subfunctionalized. We have previously shown using the medaka fish Oryzias latipes that pigment cell formation is severely affected by the loss of sox10a, yet unaffected by the loss of sox10b. Here we aimed to determine whether Sox9 is involved in the specification of pigment cell lineage. The sox9b homozygous mutation did not affect pigment cell formation, despite lethality at the early larval stages. By using sox10a, sox10b, and sox9b mutations, compound mutants were established for the sox9b and sox10 genes and pigment cell phenotypes were analyzed. Simultaneous loss of sox9b and sox10a resulted in the complete absence of melanophores and xanthophores from hatchlings and severely defective iridophore formation, as has been previously shown for sox10a-/- ; sox10b-/- double mutants, indicating that Sox9b as well as Sox10b functions redundantly with Sox10a in pigment cell development. Notably, leucophores were present in sox9b-/- ; sox10a-/- and sox10a-/- ; sox10b-/- double mutants, but their numbers were significantly reduced in the sox9b-/- ; sox10a-/- mutants. These findings highlight that Sox9b is involved in pigment cell formation, and plays a more critical role in leucophore development than Sox10b.

Novel generic models for differentiating stem cells reveal oscillatory mechanisms.

Understanding cell fate selection remains a central challenge in developmental biology. We present a class of simple yet biologically motivated mathematical models for cell differentiation that generically generate oscillations and hence suggest alternatives to the standard framework based on Waddington's epigenetic landscape. The models allow us to suggest two generic dynamical scenarios that describe the differentiation process. In the first scenario, gradual variation of a single control parameter is responsible for both entering and exiting the oscillatory regime. In the second scenario, two control parameters vary: one responsible for entering, and the other for exiting the oscillatory regime. We analyse the standard repressilator and four variants of it and show the dynamical behaviours associated with each scenario. We present a thorough analysis of the associated bifurcations and argue that gene regulatory networks with these repressilator-like characteristics are promising candidates to describe cell fate selection through an oscillatory process.

A suitable anaesthetic protocol for metamorphic zebrafish.

Zebrafish are frequently used as a means to investigate development. These studies increasingly require repeated anaesthesia of zebrafish during juvenile (i.e. metamorphic) stages. The effects of anaesthesia during this time remain poorly studied. The aim of this study was to develop a reliable method that can be used for frequently repeated anaesthesia during juvenile stages. Initially, we assessed different concentrations of MS-222, the most commonly used fish anaesthetic, for 30 minute anaesthesia with recovery. We showed that suitable MS-222 doses could be identified for the smallest (7mm) and largest (20mm) fish. However, we found that juvenile fish within a specific metamorphic window (sized between 8-16 mm) were vulnerable to MS-222 and no standard concentration of MS-222 provided reliable anaesthesia under these conditions. Hence we focussed our efforts on identifying a protocol for these stages. We tested six different published anaesthesia protocols P1-P6 where P1, P2 corresponds to 0.01% MS-222, P3, P4: 0.085% 2-phenoxyethanol and P5, P6: 0.00025%/0.0050% Propofol/Lidocaine. In protocols P1, P3, P5 fish were maintained by immersion, whilst in P2, P4 and P6: fish were maintained on an anaesthetic-doused cotton-pad. We assessed reliable anaesthesia using 10 fish for 10 minutes, with full recovery. Our data allowed us to eliminate two of these protocols as unsuitable for short term anaesthesia with recovery of juvenile fish. Extending these studies to explore repeated anaesthesia at 4 day intervals for 20 days under the remaining four protocols, we showed that P1 and P4 were both suitable for repeated anaesthesia, and that P4 was most suitable for imaging. We confirmed that P4 remained suitable when the frequency of anaesthesia was increased to every 2 days. We conclude that this protocol provides a refinement to the current protocol for repeated anaesthesia with recovery of juvenile zebrafish in the vulnerable metamorphic window.

The MITF paralog tfec is required in neural crest development for fate specification of the iridophore lineage from a multipotent pigment cell progenitor.

Understanding how fate specification of distinct cell-types from multipotent progenitors occurs is a fundamental question in embryology. Neural crest stem cells (NCSCs) generate extraordinarily diverse derivatives, including multiple neural, skeletogenic and pigment cell fates. Key transcription factors and extracellular signals specifying NCSC lineages remain to be identified, and we have only a little idea of how and when they function together to control fate. Zebrafish have three neural crest-derived pigment cell types, black melanocytes, light-reflecting iridophores and yellow xanthophores, which offer a powerful model for studying the molecular and cellular mechanisms of fate segregation. Mitfa has been identified as the master regulator of melanocyte fate. Here, we show that an Mitf-related transcription factor, Tfec, functions as master regulator of the iridophore fate. Surprisingly, our phenotypic analysis of tfec mutants demonstrates that Tfec also functions in the initial specification of all three pigment cell-types, although the melanocyte and xanthophore lineages recover later. We show that Mitfa represses tfec expression, revealing a likely mechanism contributing to the decision between melanocyte and iridophore fate. Our data are consistent with the long-standing proposal of a tripotent progenitor restricted to pigment cell fates. Moreover, we investigate activation, maintenance and function of tfec in multipotent NCSCs, demonstrating for the first time its role in the gene regulatory network forming and maintaining early neural crest cells. In summary, we build on our previous work to characterise the gene regulatory network governing iridophore development, establishing Tfec as the master regulator driving iridophore specification from multipotent progenitors, while shedding light on possible cellular mechanisms of progressive fate restriction.

Cyclical fate restriction: a new view of neural crest cell fate specification.

Neural crest cells are crucial in development, not least because of their remarkable multipotency. Early findings stimulated two hypotheses for how fate specification and commitment from fully multipotent neural crest cells might occur, progressive fate restriction (PFR) and direct fate restriction, differing in whether partially restricted intermediates were involved. Initially hotly debated, they remain unreconciled, although PFR has become favoured. However, testing of a PFR hypothesis of zebrafish pigment cell development refutes this view. We propose a novel 'cyclical fate restriction' hypothesis, based upon a more dynamic view of transcriptional states, reconciling the experimental evidence underpinning the traditional hypotheses.

Taste buds are not derived from neural crest in mouse, chicken, and zebrafish.

Our lineage tracing studies using multiple Cre mouse lines showed a concurrent labeling of abundant taste bud cells and the underlying connective tissue with a neural crest (NC) origin, warranting a further examination on the issue of whether there is an NC derivation of taste bud cells. In this study, we mapped NC cell lineages in three different models, Sox10-iCreERT2/tdT mouse, GFP+ neural fold transplantation to GFP- chickens, and Sox10-Cre/GFP-RFP zebrafish model. We found that in mice, Sox10-iCreERT2 specifically labels NC cell lineages with a single dose of tamoxifen at E7.5 and that the labeled cells were widely distributed in the connective tissue of the tongue. No labeled cells were found in taste buds or the surrounding epithelium in the postnatal mice. In the GFP+/GFP- chicken chimera model, GFP+ cells migrated extensively to the cranial region of chicken embryos ipsilateral to the surgery side but were absent in taste buds in the base of oral cavity and palate. In zebrafish, Sox10-Cre/GFP-RFP faithfully labeled known NC-derived tissues but did not label taste buds in lower jaw or the barbel. Our data, together with previous findings in axolotl, indicate that taste buds are not derived from NC cells in rodents, birds, amphibians or teleost fish.

Enteric glia as a source of neural progenitors in adult zebrafish.

The presence and identity of neural progenitors in the enteric nervous system (ENS) of vertebrates is a matter of intense debate. Here, we demonstrate that the non-neuronal ENS cell compartment of teleosts shares molecular and morphological characteristics with mammalian enteric glia but cannot be identified by the expression of canonical glial markers. However, unlike their mammalian counterparts, which are generally quiescent and do not undergo neuronal differentiation during homeostasis, we show that a relatively high proportion of zebrafish enteric glia proliferate under physiological conditions giving rise to progeny that differentiate into enteric neurons. We also provide evidence that, similar to brain neural stem cells, the activation and neuronal differentiation of enteric glia are regulated by Notch signalling. Our experiments reveal remarkable similarities between enteric glia and brain neural stem cells in teleosts and open new possibilities for use of mammalian enteric glia as a potential source of neurons to restore the activity of intestinal neural circuits compromised by injury or disease.

A quantitative modelling approach to zebrafish pigment pattern formation.

Pattern formation is a key aspect of development. Adult zebrafish exhibit a striking striped pattern generated through the self-organisation of three different chromatophores. Numerous investigations have revealed a multitude of individual cell-cell interactions important for this self-organisation, but it has remained unclear whether these known biological rules were sufficient to explain pattern formation. To test this, we present an individual-based mathematical model incorporating all the important cell-types and known interactions. The model qualitatively and quantitatively reproduces wild type and mutant pigment pattern development. We use it to resolve a number of outstanding biological uncertainties, including the roles of domain growth and the initial iridophore stripe, and to generate hypotheses about the functions of leopard. We conclude that our rule-set is sufficient to recapitulate wild-type and mutant patterns. Our work now leads the way for further in silico exploration of the developmental and evolutionary implications of this pigment patterning system.

Dicer1 is required for pigment cell and craniofacial development in zebrafish.

The multidomain RNase III endoribonuclease DICER is required for the generation of most functional microRNAs (miRNAs). Loss of Dicer affects developmental processes at different levels. Here, we characterized the zebrafish Dicer1 mutant, dicer1sa9205, which has a single point mutation induced by N-ethyl-N-nitrosourea mutagenesis. Heterozygous dicer1sa9205 developed normally, being phenotypically indistinguishable from wild-type siblings. Homozygous dicer1sa9205 mutants display smaller eyes, abnormal craniofacial development and aberrant pigmentation. Reduced numbers of both iridophores and melanocytes were observed in the head and ventral trunk of dicer1sa9205 homozygotes; the effect on melanocytes was stronger and detectable earlier in development. The expression of microphthalmia-associated transcription factor a (mitfa), the master gene for melanocytes differentiation, was enhanced in dicer1-depleted fish. Similarly, the expression of SRY-box containing gene 10 (sox10), required for mitfa activation, was higher in mutants than in wild types. In silico and in vivo analyses of either sox10 or mitfa 3'UTRs revealed conserved potential miRNA binding sites likely involved in the post-transcriptional regulation of both genes. Based on these findings, we propose that dicer1 participates in the gene regulatory network governing zebrafish melanocyte differentiation by controlling the expression of mitfa and sox10.

Endothelin neurotransmitter signalling controls zebrafish social behaviour.

The formation of social groups is an adaptive behaviour that can provide protection from predators, improve foraging and facilitate social learning. However, the costs of proximity can include competition for resources, aggression and kleptoparasitism meaning that the decision whether to interact represents a trade-off. Here we show that zebrafish harbouring a mutation in endothelin receptor aa (ednraa) form less cohesive shoals than wild-types. ednraa-/- mutants exhibit heightened aggression and decreased whole-body cortisol levels suggesting that they are dominant. These behavioural changes correlate with a reduction of parvocellular arginine vasopressin (AVP)-positive neurons in the preoptic area, an increase in the size of magnocellular AVP neurons and a higher concentration of 5-HT and dopamine in the brain. Manipulation of AVP or 5-HT signalling can rescue the shoaling phenotype of ednraa-/- providing an insight into how the brain controls social interactions.

Endothelin receptor Aa regulates proliferation and differentiation of Erb-dependent pigment progenitors in zebrafish.

Skin pigment patterns are important, being under strong selection for multiple roles including camouflage and UV protection. Pigment cells underlying these patterns form from adult pigment stem cells (APSCs). In zebrafish, APSCs derive from embryonic neural crest cells, but sit dormant until activated to produce pigment cells during metamorphosis. The APSCs are set-aside in an ErbB signaling dependent manner, but the mechanism maintaining quiescence until metamorphosis remains unknown. Mutants for a pigment pattern gene, parade, exhibit ectopic pigment cells localised to the ventral trunk, but also supernumerary cells restricted to the Ventral Stripe. Contrary to expectations, these melanocytes and iridophores are discrete cells, but closely apposed. We show that parade encodes Endothelin receptor Aa, expressed in the blood vessels, most prominently in the medial blood vessels, consistent with the ventral trunk phenotype. We provide evidence that neuronal fates are not affected in parade mutants, arguing against transdifferentiation of sympathetic neurons to pigment cells. We show that inhibition of BMP signaling prevents specification of sympathetic neurons, indicating conservation of this molecular mechanism with chick and mouse. However, inhibition of sympathetic neuron differentiation does not enhance the parade phenotype. Instead, we pinpoint ventral trunk-restricted proliferation of neural crest cells as an early feature of the parade phenotype. Importantly, using a chemical genetic screen for rescue of the ectopic pigment cell phenotype of parade mutants (whilst leaving the embryonic pattern untouched), we identify ErbB inhibitors as a key hit. The time-window of sensitivity to these inhibitors mirrors precisely the window defined previously as crucial for the setting aside of APSCs in the embryo, strongly implicating adult pigment stem cells as the source of the ectopic pigment cells. We propose that a novel population of APSCs exists in association with medial blood vessels, and that their quiescence is dependent upon Endothelin-dependent factors expressed by the blood vessels.