A systems biology approach uncovers the core gene regulatory network governing iridophore fate choice from the neural crest
Authors: Petratou K. , Subkhankulova T. , Lister J. A. , Rocco A. , Schwetlick H. , Kelsh R. N. 
Published: May 2018

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. 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 novel, 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 asox10-dependent positive feedback loop betweentfecandltkdriving iridophore specification; the molecular basis ofsox10maintenance throughout iridophore development; and the cooperation betweensox10andtfecin driving expression ofpnp4a, a key differentiation gene. We also assess a candidate repressor ofmitfa, a melanocyte-specific target ofsox10. Surprisingly, our data challenge the reported role of Foxd3, an establishedmitfarepressor, 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.

Author Summary

Multipotent neural crest (NC) progenitors generate an astonishing array of derivatives, including neuronal, skeletal components and pigment cells, 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 present expression data, as well as loss and gain of function results, guiding the derivation of a core iridophore specification GRN. Moreover, we use a process of mathematical modelling and rigorous computational exploration of the GRN to predict gene expression dynamics, assessing them by criteria suited to the qualitative nature of our current understanding of iridophore development. 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. The core iridophore GRN defined here is a key stepping stone towards exploring how chromatophores fate decisions are made in multipotent NC progenitors.