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.

A systems biology approach uncovers the core gene regulatory network governing iridophore fate choice from the neural crest.

Multipotent neural crest (NC) progenitors generate an astonishing array of derivatives, including neuronal, skeletal components and pigment cells (chromatophores), but the molecular mechanisms allowing balanced selection of each fate remain unknown. In zebrafish, melanocytes, iridophores and xanthophores, the three chromatophore lineages, are thought to share progenitors and so lend themselves to investigating the complex gene regulatory networks (GRNs) underlying fate segregation of NC progenitors. Although the core GRN governing melanocyte specification has been previously established, those guiding iridophore and xanthophore development remain elusive. Here we focus on the iridophore GRN, where mutant phenotypes identify the transcription factors Sox10, Tfec and Mitfa and the receptor tyrosine kinase, Ltk, as key players. Here we present expression data, as well as loss and gain of function results, guiding the derivation of an initial iridophore specification GRN. Moreover, we use an iterative process of mathematical modelling, supplemented with a Monte Carlo screening algorithm suited to the qualitative nature of the experimental data, to allow for rigorous predictive exploration of the GRN dynamics. Predictions were experimentally evaluated and testable hypotheses were derived to construct an improved version of the GRN, which we showed produced outputs consistent with experimentally observed gene expression dynamics. Our study reveals multiple important regulatory features, notably a sox10-dependent positive feedback loop between tfec and ltk driving iridophore specification; the molecular basis of sox10 maintenance throughout iridophore development; and the cooperation between sox10 and tfec in driving expression of pnp4a, a key differentiation gene. We also assess a candidate repressor of mitfa, a melanocyte-specific target of sox10. Surprisingly, our data challenge the reported role of Foxd3, an established mitfa repressor, in iridophore regulation. Our study builds upon our previous systems biology approach, by incorporating physiologically-relevant parameter values and rigorous evaluation of parameter values within a qualitative data framework, to establish for the first time the core GRN guiding specification of the iridophore lineage.

Calculation of long time classical trajectories: algorithmic treatment and applications for molecular systems.

We study the problem of finding a path that joins a given initial state with a final one, where the evolution is governed by classical (Hamiltonian) dynamics. A new algorithm for the computation of long time transition trajectories connecting two configurations is presented. In particular, a strategy for finding transition paths between two stable basins is established. The starting point is the formulation of the equation of motion of classical mechanics in the framework of Jacobi's principle; a shortening procedure inspired by Birkhoff's method is then applied to find geodesic solutions. Numerical examples are given for Muller's potential and the collinear reaction H(2) + H --> H + H(2).