Zebrafish pigment cells develop directly from persistent highly multipotent progenitors

Neural crest cells (NCCs) are highly multipotent stem cells. A long-standing controversy exists over the mechanism of NCC fate specification, specifically regarding the presence and potency of intermediate progenitors. The direct fate restriction (DFR) model, based on early in vivo clonal studies, hypothesised that intermediates are absent and that migrating cells maintain full multipotency1-6. However, most authors favour progressive fate restriction (PFR) models, with fully multipotent early NCCs (ENCCs) transitioning to partially- restricted intermediates before committing to individual fates7-12. Here, single cell transcriptional profiling of zebrafish pigment cell development leads to us proposing a Cyclical Fate Restriction mechanism of NCC development that reconciles the DFR and PFR models. Our clustering of single NCC Nanostring transcriptional profiles identifies only broadly multipotent intermediate states between ENCCs and differentiated melanocytes and iridophores. Leukocyte tyrosine kinase (Ltk) marks the multipotent progenitor and iridophores, consistent with biphasic ltk expression13-15. Ltk inhibitor and constitutive activation studies support expression at an early multipotent stage, whilst lineage-tracing of ltk-expressing cells reveals their multipotency extends beyond pigment cell-types to neural fates. We conclude that pigment cell development does not involve a conventional PFR mechanism, but instead occurs directly and more dynamically from a broadly multipotent intermediate state. We investigated the transcriptional profiles of 1227 zebrafish NCCs throughout pigment cell development (18-72 hpf). To obtain sensitive quantitation of key mRNA expression, we used the NanoString platform and a set of 45 genes, focused on those known or suspected to have important roles in melanocyte or iridophore fate choice, but including multiple markers of early NCCs, plus key markers of other major NC-derived fates Processed data are provided only for samples GSM5619061-GSM5619791.