Atoms to Phenotypes: Molecular Design Principles of Cellular Energy Metabolism.

We report a 100-million atom-scale model of an entire cell organelle, a photosynthetic chromatophore vesicle from a purple bacterium, that reveals the cascade of energy conversion steps culminating in the generation of ATP from sunlight. Molecular dynamics simulations of this vesicle elucidate how the integral membrane complexes influence local curvature to tune photoexcitation of pigments. Brownian dynamics of small molecules within the chromatophore probe the mechanisms of directional charge transport under various pH and salinity conditions. Reproducing phenotypic properties from atomistic details, a kinetic model evinces that low-light adaptations of the bacterium emerge as a spontaneous outcome of optimizing the balance between the chromatophore's structural integrity and robust energy conversion. Parallels are drawn with the more universal mitochondrial bioenergetic machinery, from whence molecular-scale insights into the mechanism of cellular aging are inferred. Together, our integrative method and spectroscopic experiments pave the way to first-principles modeling of whole living cells.

The dynamic behavior of chromatophores marks the transition from bands to spots in leopard geckos.

Reptilian skin coloration is spectacular and diverse, yet little is known about the ontogenetic processes that govern its establishment and the molecular signaling pathways that determine it. Here, we focus on the development of the banded pattern of leopard gecko hatchlings and the transition to black spots in the adult. With our histological analyses, we show that iridophores are present in the white and yellow bands of the hatchling and they gradually perish in the adult skin. Furthermore, we demonstrate that melanophores can autonomously form spots in the absence of the other chromatophores both on the regenerated skin of the tail and on the dorsal skin of the Mack Super Snow (MSS) leopard geckos. This color morph is characterized by uniform black coloration in hatchlings and black spots in adulthood; we establish that their skin is devoid of xanthophores and iridophores at both stages. Our genetic analyses identified a 13-nucleotide deletion in the PAX7 transcription factor of MSS geckos, affecting its protein coding sequence. With our single-cell transcriptomics analysis of embryonic skin, we confirm that PAX7 is expressed in iridophores and xanthophores, suggesting that it plays a key role in the differentiation of both chromatophores. Our in situ hybridizations on whole-mount embryos document the dynamics of the skin pattern formation and how it is impacted in the PAX7 mutants. We hypothesize that the melanophores-iridophores interactions give rise to the banded pattern of the hatchlings and black spot formation is an intrinsic capacity of melanophores in the postembryonic skin.

DNA-binding and protein structure of nuclear factors likely acting in genetic information processing in the Paulinella chromatophore.

The chromatophores in Paulinella are evolutionary-early-stage photosynthetic organelles. Biological processes in chromatophores depend on a combination of chromatophore and nucleus-encoded proteins. Interestingly, besides proteins carrying chromatophore-targeting signals, a large arsenal of short chromatophore-targeted proteins (sCTPs; <90 amino acids) without recognizable targeting signals were found in chromatophores. This situation resembles endosymbionts in plants and insects that are manipulated by host-derived antimicrobial peptides. Previously, we identified an expanded family of sCTPs of unknown function, named here "DNA-binding (DB)-sCTPs". DB-sCTPs contain a ~45 amino acid motif that is conserved in some bacterial proteins with predicted functions in DNA processing. Here, we explored antimicrobial activity, DNA-binding capacity, and structures of three purified recombinant DB-sCTPs. All three proteins exhibited antimicrobial activity against bacteria involving membrane permeabilization, and bound to bacterial lipids in vitro. A combination of in vitro assays demonstrated binding of recombinant DB-sCTPs to chromatophore-derived genomic DNA sequences with an affinity in the low nM range. Additionally, we report the 1.2 Å crystal structure of one DB-sCTP. In silico docking studies suggest that helix α2 inserts into the DNA major grove and the exposed residues, that are highly variable between different DB-sCTPs, confer interaction with the DNA bases. Identification of photosystem II subunit CP43 as a potential interaction partner of one DB-sCTP, suggests DB-sCTPs to be involved in more complex regulatory mechanisms. We hypothesize that membrane binding of DB-sCTPs is related to their import into chromatophores. Once inside, they interact with the chromatophore genome potentially providing nuclear control over genetic information processing.

Lizards exploit the changing optics of developing chromatophore cells to switch defensive colors during ontogeny.

Many animals undergo changes in functional colors during development, requiring the replacement of integument or pigment cells. A classic example of defensive color switching is found in hatchling lizards, which use conspicuous tail colors to deflect predator attacks away from vital organs. These tail colors usually fade to concealing colors during ontogeny. Here, we show that the ontogenetic blue-to-brown tail color change in Acanthodactylus beershebensis lizards results from the changing optical properties of single types of developing chromatophore cells. The blue tail colors of hatchlings are produced by incoherent scattering from premature guanine crystals in underdeveloped iridophore cells. Cryptic tail colors emerge during chromatophore maturation upon reorganization of the guanine crystals into a multilayer reflector concomitantly with pigment deposition in the xanthophores. Ontogenetic changes in adaptive colors can thus arise not via the exchange of different optical systems, but by harnessing the timing of natural chromatophore development. The incoherent scattering blue color here differs from the multilayer interference mechanism used in other blue-tailed lizards, indicating that a similar trait can be generated in at least two ways. This supports a phylogenetic analysis showing that conspicuous tail colors are prevalent in lizards and that they evolved convergently. Our results provide an explanation for why certain lizards lose their defensive colors during ontogeny and yield a hypothesis for the evolution of transiently functional adaptive colors.

Elucidating the control and development of skin patterning in cuttlefish.

Few animals provide a readout that is as objective of their perceptual state as camouflaging cephalopods. Their skin display system includes an extensive array of pigment cells (chromatophores), each expandable by radial muscles controlled by motor neurons. If one could track the individual expansion states of the chromatophores, one would obtain a quantitative description-and potentially even a neural description by proxy-of the perceptual state of the animal in real time. Here we present the use of computational and analytical methods to achieve this in behaving animals, quantifying the states of tens of thousands of chromatophores at sixty frames per second, at single-cell resolution, and over weeks. We infer a statistical hierarchy of motor control, reveal an underlying low-dimensional structure to pattern dynamics and uncover rules that govern the development of skin patterns. This approach provides an objective description of complex perceptual behaviour, and a powerful means to uncover the organizational principles that underlie the function, dynamics and morphogenesis of neural systems.

Molecular parallelisms between pigmentation in the avian iris and the integument of ectothermic vertebrates.

Birds exhibit striking variation in eye color that arises from interactions between specialized pigment cells named chromatophores. The types of chromatophores present in the avian iris are lacking from the integument of birds or mammals, but are remarkably similar to those found in the skin of ectothermic vertebrates. To investigate molecular mechanisms associated with eye coloration in birds, we took advantage of a Mendelian mutation found in domestic pigeons that alters the deposition of yellow pterin pigments in the iris. Using a combination of genome-wide association analysis and linkage information in pedigrees, we mapped variation in eye coloration in pigeons to a small genomic region of ~8.5kb. This interval contained a single gene, SLC2A11B, which has been previously implicated in skin pigmentation and chromatophore differentiation in fish. Loss of yellow pigmentation is likely caused by a point mutation that introduces a premature STOP codon and leads to lower expression of SLC2A11B through nonsense-mediated mRNA decay. There were no substantial changes in overall gene expression profiles between both iris types as well as in genes directly associated with pterin metabolism and/or chromatophore differentiation. Our findings demonstrate that SLC2A11B is required for the expression of pterin-based pigmentation in the avian iris. They further highlight common molecular mechanisms underlying the production of coloration in the iris of birds and skin of ectothermic vertebrates.

A bipartite chromatophore transit peptide and N-terminal protein processing in the Paulinella chromatophore.

The amoeba Paulinella chromatophora contains photosynthetic organelles, termed chromatophores, which evolved independently from plastids in plants and algae. At least one-third of the chromatophore proteome consists of nucleus-encoded (NE) proteins that are imported across the chromatophore double envelope membranes. Chromatophore-targeted proteins exceeding 250 amino acids (aa) carry a conserved N-terminal extension presumably involved in protein targeting, termed the chromatophore transit peptide (crTP). Short imported proteins do not carry discernable targeting signals. To explore whether the import of proteins is accompanied by their N-terminal processing, here we identified N-termini of 208 chromatophore-localized proteins by a mass spectrometry-based approach. Our study revealed extensive N-terminal acetylation and proteolytic processing in both NE and chromatophore-encoded (CE) fractions of the chromatophore proteome. Mature N-termini of 37 crTP-carrying proteins were identified, of which 30 were cleaved in a common processing region. Surprisingly, only the N-terminal ∼50 aa (part 1) become cleaved upon import. This part contains a conserved adaptor protein-1 complex-binding motif known to mediate protein sorting at the trans-Golgi network followed by a predicted transmembrane helix, implying that part 1 anchors the protein co-translationally in the endoplasmic reticulum and mediates trafficking to the chromatophore via the Golgi. The C-terminal part 2 contains conserved secondary structural elements, remains attached to the mature proteins, and might mediate translocation across the chromatophore inner membrane. Short imported proteins remain largely unprocessed. Finally, this work illuminates N-terminal processing of proteins encoded in an evolutionary-early-stage organelle and suggests host-derived posttranslationally acting factors involved in regulation of the CE chromatophore proteome.

Generation of Electric Potential Difference by Chromatophores from Photosynthetic Bacteria in the Presence of Trehalose under Continuous Illumination.

Measurement of electrical potential difference (Δψ) in membrane vesicles (chromatophores) from the purple bacterium Rhodobacter sphaeroides associated with the surface of a nitrocellulose membrane filter (MF) impregnated with a phospholipid solution in decane or immersed into it in the presence of exogenous mediators and disaccharide trehalose demonstrated an increase in the amplitude and stabilization of the signal under continuous illumination. The mediators were the ascorbate/N,N,N'N'-tetramethyl-p-phenylenediamine pair and ubiquinone-0 (electron donor and acceptor, respectively). Although stabilization of photoelectric responses upon long-term continuous illumination was observed for both variants of chromatophore immobilization, only the samples immersed into the MF retained the functional activity of reaction centers (RCs) for a month when stored in the dark at room temperature, which might be due to the preservation of integrity of chromatophore proteins inside the MF pores. The stabilizing effect of the bioprotector trehalose could be related to its effect on both the RC proteins and the phospholipid bilayer membrane. The results obtained will expand current ideas on the use of semi-synthetic structures based on various intact photosynthetic systems capable of converting solar energy into its electrochemical form.

Bioluminescent backlighting illuminates the complex visual signals of a social squid in the deep sea.

Visual signals rapidly relay information, facilitating behaviors and ecological interactions that shape ecosystems. However, most known signaling systems can be restricted by low light levels-a pervasive condition in the deep ocean, the largest inhabitable space on the planet. Resident visually cued animals have therefore been hypothesized to have simple signals with limited information-carrying capacity. We used cameras mounted on remotely operated vehicles to study the behavior of the Humboldt squid, Dosidicus gigas, in its natural deep-sea habitat. We show that specific pigmentation patterns from its diverse repertoire are selectively displayed during foraging and in social scenarios, and we investigate how these behaviors may be used syntactically for communication. We additionally identify the probable mechanism by which D. gigas, and related squids, illuminate these patterns to create visual signals that can be readily perceived in the deep, dark ocean. Numerous small subcutaneous (s.c.) photophores (bioluminescent organs) embedded throughout the muscle tissue make the entire body glow, thereby backlighting the pigmentation patterns. Equipped with a mechanism by which complex information can be rapidly relayed through a visual pathway under low-light conditions, our data suggest that the visual signals displayed by D. gigas could share design features with advanced forms of animal communication. Visual signaling by deep-living cephalopods will likely be critical in understanding how, and how much, information can be shared in one of the planet's most challenging environments for visual communication.

Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles.

Reptiles exhibit a spectacular diversity of skin colors and patterns brought about by the interactions among three chromatophore types: black melanophores with melanin-packed melanosomes, red and yellow xanthophores with pteridine- and/or carotenoid-containing vesicles, and iridophores filled with light-reflecting platelets generating structural colors. Whereas the melanosome, the only color-producing endosome in mammals and birds, has been documented as a lysosome-related organelle, the maturation paths of xanthosomes and iridosomes are unknown. Here, we first use 10x Genomics linked-reads and optical mapping to assemble and annotate a nearly chromosome-quality genome of the corn snake Pantherophis guttatus The assembly is 1.71 Gb long, with an N50 of 16.8 Mb and L50 of 24. Second, we perform mapping-by-sequencing analyses and identify a 3.9-Mb genomic interval where the lavender variant resides. The lavender color morph in corn snakes is characterized by gray, rather than red, blotches on a pink, instead of orange, background. Third, our sequencing analyses reveal a single nucleotide polymorphism introducing a premature stop codon in the lysosomal trafficking regulator gene (LYST) that shortens the corresponding protein by 603 amino acids and removes evolutionary-conserved domains. Fourth, we use light and transmission electron microscopy comparative analyses of wild type versus lavender corn snakes and show that the color-producing endosomes of all chromatophores are substantially affected in the LYST mutant. Our work provides evidence characterizing xanthosomes in xanthophores and iridosomes in iridophores as lysosome-related organelles.

Exploring the mechanisms and impacts of melatonin on fish colouration.

The pineal hormone melatonin is a multi-functional molecule with a recognized role in pigment aggregation in chromatophores, mediating its actions through binding to subtypes of its specific receptors. Since its discovery, melatonin has been known to be responsible for pigment aggregation towards the cell centre in fishes, including their embryos, as an adaptation to reduced light and thus results in pale body colouration. Diversity exists in the sensitivity of melanophores towards melatonin at interspecies, intraspecific levels, seasons, and amongst chromatophores at different regions of the animal body. In most of the fishes, melatonin leads to their skin paling at night. It is indicated that the melatonin receptors have characteristically maintained to show the same aggregating effects in fishes and other vertebrates in the evolutionary hierarchy. However, besides this aggregatory effect, melatonin is also responsible for pigment dispersion in certain fishes. Here is the demand in our review to explore further the nature of the dispersive behaviour of melatonin through the so-called ß-melatonin receptors. It is clear that the pigment translocations in lower vertebrates under the effect of melatonin are mediated through the melatonin receptors coupled with other hormonal receptors as well. Therefore, being richly supplied with a variety of receptors, chromatophores and melanocytes can be used as in vitro test models for pharmacological applications of known and novel drugs. In this review, we present diverse effects of melatonin on chromatophores of fishes in particular with appropriate implications on most of the recent findings.

Amoeba Genome Reveals Dominant Host Contribution to Plastid Endosymbiosis.

Eukaryotic photosynthetic organelles, plastids, are the powerhouses of many aquatic and terrestrial ecosystems. The canonical plastid in algae and plants originated >1 Ga and therefore offers limited insights into the initial stages of organelle evolution. To address this issue, we focus here on the photosynthetic amoeba Paulinella micropora strain KR01 (hereafter, KR01) that underwent a more recent (∼124 Ma) primary endosymbiosis, resulting in a photosynthetic organelle termed the chromatophore. Analysis of genomic and transcriptomic data resulted in a high-quality draft assembly of size 707 Mb and 32,361 predicted gene models. A total of 291 chromatophore-targeted proteins were predicted in silico, 208 of which comprise the ancestral organelle proteome in photosynthetic Paulinella species with functions, among others, in nucleotide metabolism and oxidative stress response. Gene coexpression analysis identified networks containing known high light stress response genes as well as a variety of genes of unknown function ("dark" genes). We characterized diurnally rhythmic genes in this species and found that over 49% are dark. It was recently hypothesized that large double-stranded DNA viruses may have driven gene transfer to the nucleus in Paulinella and facilitated endosymbiosis. Our analyses do not support this idea, but rather suggest that these viruses in the KR01 and closely related P. micropora MYN1 genomes resulted from a more recent invasion.

Evidence for a robust photosystem II in the photosynthetic amoeba Paulinella.

Paulinella represents the only known case of an independent primary plastid endosymbiosis, outside Archaeplastida, that occurred c. 120 (million years ago) Ma. These photoautotrophs grow very slowly in replete culture medium with a doubling time of 6-7 d at optimal low light, and are highly sensitive to photodamage under moderate light levels. We used genomic and biophysical methods to investigate the extreme slow growth rate and light sensitivity of Paulinella, which are key to photosymbiont integration. All photosystem II (PSII) genes except psb28-2 and all cytochrome b6 f complex genes except petM and petL are present in Paulinella micropora KR01 (hereafter, KR01). Biophysical measurements of the water oxidation complex, variable chlorophyll fluorescence, and photosynthesis-irradiance curves show no obvious evidence of PSII impairment. Analysis of photoacclimation under high-light suggests that although KR01 can perform charge separation, it lacks photoprotection mechanisms present in cyanobacteria. We hypothesize that Paulinella species are restricted to low light environments because they are deficient in mitigating the formation of reactive oxygen species formed within the photosystems under peak solar intensities. The finding that many photoprotection genes have been lost or transferred to the host-genome during endosymbiont genome reduction, and may lack light-regulation, is consistent with this hypothesis.

Disproportionate effect of cationic antiseptics on the quantum yield and fluorescence lifetime of bacteriochlorophyll molecules in the LH1-RC complex of R. rubrum chromatophores.

Photosynthetic membrane complexes of purple bacteria are convenient and informative macromolecular systems for studying the mechanisms of action of various physicochemical factors on the functioning of catalytic proteins both in an isolated state and as part of functional membranes. In this work, we studied the effect of cationic antiseptics (chlorhexidine, picloxydine, miramistin, and octenidine) on the fluorescence intensity and the efficiency of energy transfer from the light-harvesting LH1 complex to the reaction center (RC) of Rhodospirillum rubrum chromatophores. The effect of antiseptics on the fluorescence intensity and the energy transfer increased in the following order: chlorhexidine, picloxydine, miramistin, octenidine. The most pronounced changes in the intensity and lifetime of fluorescence were observed with the addition of miramistin and octenidine. At the same concentration of antiseptics, the increase in fluorescence intensity was 2-3 times higher than the increase in its lifetime. It is concluded that the addition of antiseptics decreases the efficiency of the energy migration LH1 → RC and increases the fluorescence rate constant kfl. We associate the latter with a change in the polarization of the microenvironment of bacteriochlorophyll molecules upon the addition of charged antiseptic molecules. A possible mechanism of antiseptic action on R. rubrum chromatophores is considered. This work is a continuation of the study of the effect of antiseptics on the energy transfer and fluorescence intensity in chromatophores of purple bacteria published earlier in Photosynthesis Research (Strakhovskaya et al. in Photosyn Res 147:197-209, 2021).

The physiological cost of colour change: evidence, implications and mitigations.

Animals benefit from phenotypic plasticity in changing environments, but this can come at a cost. Colour change, used for camouflage, communication, thermoregulation and UV protection, represents one of the most common plastic traits in nature and is categorised as morphological or physiological depending on the mechanism and speed of the change. Colour change has been assumed to carry physiological costs, but current knowledge has not advanced beyond this basic assumption. The costs of changing colour will shape the evolution of colour change in animals, yet no coherent research has been conducted in this area, leaving a gap in our understanding. Therefore, in this Review, we examine the direct and indirect evidence of the physiological cost of colour change from the cellular to the population level, in animals that utilise chromatophores in colour change. Our Review concludes that the physiological costs result from either one or a combination of the processes of (i) production, (ii) translocation and (iii) maintenance of pigments within the colour-containing cells (chromatophores). In addition, both types of colour change (morphological and physiological) pose costs as they require energy for hormone production and neural signalling. Moreover, our Review upholds the hypothesis that, if repetitively used, rapid colour change (i.e. seconds-minutes) is more costly than slow colour change (days-weeks) given that rapidly colour-changing animals show mitigations, such as avoiding colour change when possible. We discuss the potential implications of this cost on colour change, behaviour and evolution of colour-changing animals, generating testable hypotheses and emphasising the need for future work to address this gap.

Effect of Dipyridamole on Membrane Energization and Energy Transfer in Chromatophores of Rba. sphaeroides.

Effect of dipyridamole (DIP) at concentrations up to 1 mM on fluorescent characteristics of light-harvesting complexes LH2 and LH1, as well as on conditions of photosynthetic electron transport chain in the bacterial chromatophores of Rba. sphaeroides was investigated. DIP was found to affect efficiency of energy transfer from the light-harvesting complex LH2 to the LH1-reaction center core complex and to produce the long-wavelength ("red") shift of the absorption band of light-harvesting bacteriochlorophyll molecules in the IR spectral region at 840-900 nm. This shift is associated with the membrane transition to the energized state. It was shown that DIP is able to reduce the photooxidized bacteriochlorophyll of the reaction center, which accelerated electron flow along the electron transport chain, thereby stimulating generation of the transmembrane potential on the chromatophore membrane. The results are important for clarifying possible mechanisms of DIP influence on the activity of membrane-bound functional proteins. In particular, they might be significant for interpreting numerous therapeutic effects of DIP.

Generation of no-yellow-pigment Xenopus tropicalis by slc2a7 gene knockout.

BACKGROUND: Amphibians possess three kinds of dermal chromatophore: melanophores, iridophores, and xanthophores. Knockout Xenopus tropicalis that lack the pigmentation of melanophores and iridophores have been reported. The identification of the causal genes for xanthophore pigmentation or differentiation could lead to the creation of a see-through frog without three chromatophores. The genes causing xanthophore differentiation mutants are slc2a11b and slc2a15b in Japanese medaka (Oryzias latipes). RESULTS: To obtain a heritable line of X tropicalis mutants without yellow pigment, we generated slc2a7 and slc2a15a knockout animals because they have the greatest similarity to the O latipes slc2a11b and slc2a15b genes. The slc2a7 knockout frog had a bluish skin and there were no visible yellow pigments in stereo microscope and skin section observations. Furthermore, no pterinosomes, which are characteristic of xanthophores, were observed via transmission electron microscopy in the skin of knockout animals. CONCLUSIONS: We report the successful generation of a heritable no-yellow-pigment X tropicalis mutant after knock out of the slc2a7 gene. This finding will enable the creation of a see-through frog with no chromatophores.

Visualization of Sox10-positive chromatoblasts by GFP fluorescence in flounder larvae and juveniles using electroporation.

Japanese flounder are left-right asymmetrical, with features, such as dark, ocular-side specific pigmentation. This pigmentation arises during metamorphic stages, along with the asymmetric differentiation of adult-type chromatophores. Additionally, among juveniles, tank-reared specimens commonly show ectopic pigmentation on their blind sides. In both cases, neural crest-derived Sox10-positive progenitor cells at the dorsal fin base are hypothesized to contribute to chromatophore development. Here, we developed a method to visualize Sox10-positive cells via green fluorescent protein (GFP) fluorescence to directly monitor their migration and differentiation into chromatophores in vivo. Electroporation was applied to introduce GFP reporter vectors into the dorsal fin base of larvae and juveniles. Cre-loxP system vectors were also tested to enable cell labeling even after a decrease in sox10 expression levels. In larvae, undifferentiated Sox10-positive progenitor cells were labeled in the dorsal fin base, whereas newly differentiated adult-type chromatophores were seen dispersed on the ocular side. In juveniles, Sox10-positive cells were identified in the connective tissue of the dorsal fin base and observed prominently in areas of ectopic pigmentation, including several labeled melanophores. Thus, it was suggested that during metamorphic stages, Sox10-positive cells at the dorsal fin base contribute to adult-type chromatophore development, whereas in juveniles, they persist as precursors in the connective tissue, which in response to stimuli migrate to generate ectopic pigmentation. These findings contribute to elucidating pigmentation mechanisms, as well as abnormalities seen in hatchery-reared flounders. The electroporation method may be adapted to diverse animals as an accessible gene transfer method in various research fields, including developmental and biomedical studies.

Fish skin pigmentation in aquaculture: The influence of rearing conditions and its neuroendocrine regulation.

Skin pigmentation pattern is a species-specific characteristic that depends on the number and the spatial combination of several types of chromatophores. This feature can change during life, for example in the metamorphosis or reproductive cycle, or as a response to biotic and/or abiotic environmental cues (nutrition, UV incidence, surrounding luminosity, and social interactions). Fish skin pigmentation is one of the most important quality criteria dictating the market value of both aquaculture and ornamental species because it serves as an external signal to infer its welfare and the culture conditions used. For that reason, several studies have been conducted aiming to understand the mechanisms underlying fish pigmentation as well as the influence exerted by rearing conditions. In this context, the present review focuses on the current knowledge on endocrine regulation of fish pigmentation as well as on the aquaculture conditions affecting skin coloration. Available information on Iberoamerican fish species cultured is presented.

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.