<?xml version="1.0" encoding="UTF-8"?>
<!DOCTYPE article PUBLIC "-//NLM//DTD Journal Publishing DTD v2.3 20070202//EN" "journalpublishing.dtd">
<article article-type="brief-report" dtd-version="2.3" xml:lang="EN" xmlns:mml="http://www.w3.org/1998/Math/MathML" xmlns:xlink="http://www.w3.org/1999/xlink">
<front>
<journal-meta>
<journal-id journal-id-type="publisher-id">Front. Cell Dev. Biol.</journal-id>
<journal-title>Frontiers in Cell and Developmental Biology</journal-title>
<abbrev-journal-title abbrev-type="pubmed">Front. Cell Dev. Biol.</abbrev-journal-title>
<issn pub-type="epub">2296-634X</issn>
<publisher>
<publisher-name>Frontiers Media S.A.</publisher-name>
</publisher>
</journal-meta>
<article-meta>
<article-id pub-id-type="publisher-id">1358828</article-id>
<article-id pub-id-type="doi">10.3389/fcell.2024.1358828</article-id>
<article-categories>
<subj-group subj-group-type="heading">
<subject>Cell and Developmental Biology</subject>
<subj-group>
<subject>Brief Research Report</subject>
</subj-group>
</subj-group>
</article-categories>
<title-group>
<article-title>On the role of TFEC in reptilian coloration</article-title>
<alt-title alt-title-type="left-running-head">Tzika</alt-title>
<alt-title alt-title-type="right-running-head">
<ext-link ext-link-type="uri" xlink:href="https://doi.org/10.3389/fcell.2024.1358828">10.3389/fcell.2024.1358828</ext-link>
</alt-title>
</title-group>
<contrib-group>
<contrib contrib-type="author" corresp="yes">
<name>
<surname>Tzika</surname>
<given-names>Athanasia C.</given-names>
</name>
<xref ref-type="corresp" rid="c001">&#x2a;</xref>
<uri xlink:href="https://loop.frontiersin.org/people/1030464/overview"/>
<role content-type="https://credit.niso.org/contributor-roles/conceptualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/formal-analysis/"/>
<role content-type="https://credit.niso.org/contributor-roles/funding-acquisition/"/>
<role content-type="https://credit.niso.org/contributor-roles/investigation/"/>
<role content-type="https://credit.niso.org/contributor-roles/methodology/"/>
<role content-type="https://credit.niso.org/contributor-roles/project-administration/"/>
<role content-type="https://credit.niso.org/contributor-roles/validation/"/>
<role content-type="https://credit.niso.org/contributor-roles/visualization/"/>
<role content-type="https://credit.niso.org/contributor-roles/writing-original-draft/"/>
<role content-type="https://credit.niso.org/contributor-roles/Writing - review &#x26; editing/"/>
</contrib>
</contrib-group>
<aff>
<institution>Laboratory of Artificial and Natural Evolution (LANE)</institution>, <institution>Department of Genetics and Evolution</institution>, <institution>University of Geneva</institution>, <addr-line>Geneva</addr-line>, <country>Switzerland</country>
</aff>
<author-notes>
<fn fn-type="edited-by">
<p>
<bold>Edited by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/1593992/overview">Alejandra Garc&#xed;a-Gasca</ext-link>, National Council of Science and Technology (CONACYT), Mexico</p>
</fn>
<fn fn-type="edited-by">
<p>
<bold>Reviewed by:</bold> <ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/313118/overview">Robert Neil Kelsh</ext-link>, University of Bath, United Kingdom</p>
<p>
<ext-link ext-link-type="uri" xlink:href="https://loop.frontiersin.org/people/2624661/overview">Douglas Menke</ext-link>, University of Georgia, United States</p>
</fn>
<corresp id="c001">&#x2a;Correspondence: Athanasia C. Tzika, <email>athanasia.tzika@unige.ch</email>
</corresp>
</author-notes>
<pub-date pub-type="epub">
<day>07</day>
<month>02</month>
<year>2024</year>
</pub-date>
<pub-date pub-type="collection">
<year>2024</year>
</pub-date>
<volume>12</volume>
<elocation-id>1358828</elocation-id>
<history>
<date date-type="received">
<day>20</day>
<month>12</month>
<year>2023</year>
</date>
<date date-type="accepted">
<day>22</day>
<month>01</month>
<year>2024</year>
</date>
</history>
<permissions>
<copyright-statement>Copyright &#xa9; 2024 Tzika.</copyright-statement>
<copyright-year>2024</copyright-year>
<copyright-holder>Tzika</copyright-holder>
<license xlink:href="http://creativecommons.org/licenses/by/4.0/">
<p>This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.</p>
</license>
</permissions>
<abstract>
<p>Reptilian species, particularly snakes and lizards, are emerging models of animal coloration. Here, I focus on the role of the TFEC transcription factor in snake and lizard coloration based on a study on wild-type and piebald ball pythons. Genomic mapping previously identified a TFEC mutation linked to the piebald ball python phenotype. The association of TFEC with skin coloration was further supported by gene-editing experiments in the brown anole lizard. However, novel histological analyses presented here reveal discrepancies between the ball python and the anole TFEC mutants phenotype, cautioning against broad generalizations. Indeed, both wild-type and piebald ball pythons completely lack iridophores, whereas the TFEC anole lizard mutants lose their iridophores compared to the wild-type anole. Based on these findings, I discuss the potential role of the MiT/TFE family in skin pigmentation across vertebrate lineages and advocate the need for developmental analyses and additional gene-editing experiments to explore the reptilian coloration diversity.</p>
</abstract>
<kwd-group>
<kwd>reptiles</kwd>
<kwd>coloration</kwd>
<kwd>skin</kwd>
<kwd>chromatophores</kwd>
<kwd>ball python</kwd>
<kwd>snakes</kwd>
<kwd>TFEC</kwd>
<kwd>anole lizard</kwd>
</kwd-group>
<custom-meta-wrap>
<custom-meta>
<meta-name>section-at-acceptance</meta-name>
<meta-value>Evolutionary Developmental Biology</meta-value>
</custom-meta>
</custom-meta-wrap>
</article-meta>
</front>
<body>
<sec sec-type="intro" id="s1">
<title>1 Introduction</title>
<p>Although the zebrafish remains the main reference in the animal coloration field (<xref ref-type="bibr" rid="B10">Patterson and Parichy, 2019</xref>), reptilian species are gaining ground as new models (<xref ref-type="bibr" rid="B9">Milinkovitch and Tzika, 2007</xref>) thanks to the great diversity of color and color patterns they exhibit. Multiple studies have investigated the distribution of chromatophores in their skin at late stages of development or in adults to understand this colorful variety (<xref ref-type="bibr" rid="B7">Kuriyama et al., 2020</xref>), but little is known about the developmental processes involved in the differentiation, migration and self-organised patterning of these cells. The characterization of spontaneously-occurring mutations that affect the coloration of snakes and lizards (<xref ref-type="bibr" rid="B1">Andrade et al., 2019</xref>; <xref ref-type="bibr" rid="B16">Ullate-Agote et al., 2020</xref>; <xref ref-type="bibr" rid="B17">Ullate-Agote and Tzika, 2021</xref>; <xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>) and the establishment of gene-editing protocols using CRISPR-Cas9 in the same lineages (<xref ref-type="bibr" rid="B12">Rasys et al., 2019</xref>; <xref ref-type="bibr" rid="B15">Tzika et al., 2023</xref>) will greatly advance our understanding of these processes, but it remains to be seen how transferrable this information is within reptiles, and Squamates (snakes and lizards) in particular.</p>
<p>Garcia-Elfing <italic>et al.</italic> (<xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>) recently investigated the piebald ball python morph (<italic>Python regius</italic>), a recessive mutant phenotype characterized by the presence of white patches on its otherwise black and brown dorsal skin pattern. Their genome mapping analyses, and candidate-gene approach provide evidence for a mutation in the gene <italic>TFEC</italic> as the genetic determinant of the piebald phenotype in ball pythons. The authors also generated, through CRISPR-Cas9 gene-editing, a targeted knock-out mutation of <italic>TFEC</italic> in the brown anole lizard (<italic>Anolis sergei</italic>), resulting in the loss of iridophores and reduced body coloration, particularly in the snout, forelimbs, and hindlimbs. But as Garcia-Elfing <italic>et al.</italic> (<xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>) conclude, it remains &#x201c;unresolved whether piebald ball pythons have iridophores in either pigmented or white skin&#x201d;, and they strongly recommend histological analyses to resolve this matter. Here, I present such histological data demonstrating that both wild-type and piebald ball pythons lack iridophores. Thus, the brown anole phenotype does not recapitulate the observations in ball pythons. Nevertheless, these results significantly advance our understanding of reptilian coloration. The work of Garcia-Elfing <italic>et al.</italic> and the results presented here suggest that the same molecule can affect reptilian and vertebrate coloration in different ways.</p>
</sec>
<sec sec-type="materials|methods" id="s2">
<title>2 Materials and methods</title>
<sec id="s2-1">
<title>2.1 Animal experimentation</title>
<p>Ball pythons and corn snakes were housed and bred at the LANE animal facility running under veterinary cantonal permit no. 1008. The individuals were sampled following Swiss law regulations and under the experimentation permit GE24/33145.</p>
</sec>
<sec id="s2-2">
<title>2.2 Histology and imaging</title>
<p>Three wild-type (histology, TEM, eye) and one piebald (histology, TEM, eye) pythons and two wild-type (histology, eye) and one Palmetto (histology) corn snakes were sampled. Skin (roughly 1&#xa0;cm &#xd7; 2&#xa0;cm) and eye samples were fixed in 4% paraformaldehyde and dehydrated in ethanol before embedding in paraffin blocks. Seven-micrometer microtome sections were deparaffinized and directly imaged with the VHX-7000 (Keyence). For transmission electron microscopy, skin pieces of 1&#xa0;mm<sup>2</sup> were fixed, sectioned, and imaged as previously described (<xref ref-type="bibr" rid="B16">Ullate-Agote et al., 2020</xref>). Sample processing and imaging were performed at the Electron Microscopy Facility, University of Lausanne (Switzerland).</p>
</sec>
<sec id="s2-3">
<title>2.3 <italic>TFEC</italic> amplification</title>
<p>For all animals (<xref ref-type="table" rid="T1">Table 1</xref>), genomic DNA was extracted from whole blood using the DNeasy Blood &#x26; Tissue Kit (Qiagen, 51104) and genotyping by Sanger sequencing, targeting the <italic>TFEC</italic> SNP identified in the original publication (<xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>), was performed with the following primers: TFECg_103193F (CAG&#x200b;TGC&#x200b;AAC&#x200b;TCA&#x200b;AAG&#x200b;GGA&#x200b;ACA) and TFECg_103880R (GCA&#x200b;GAC&#x200b;CCA&#x200b;TGA&#x200b;AAT&#x200b;CAA&#x200b;TGG&#x200b;A).</p>
<table-wrap id="T1" position="float">
<label>TABLE 1</label>
<caption>
<p>Individuals genotyped for the presence of the SNP suggested to cause the piebald ball python phenotype. The SNP is highlighted in bold in the &#x2018;sequence&#x2019; column. The &#x2018;genotype&#x2019; is deduced by the appearance and the pedigree of the animals. &#x2018;Stripe&#x2019; is another ball python morph unrelated to the piebald morph.</p>
</caption>
<table>
<thead valign="top">
<tr>
<th align="center">Individual ID</th>
<th align="center">Phenotype</th>
<th align="center">Genotype</th>
<th align="center">Sequence</th>
</tr>
</thead>
<tbody valign="top">
<tr>
<td align="center">PREG003</td>
<td align="center">wild-type</td>
<td align="center">&#x2b;/&#x2b;</td>
<td align="center">CAC&#x200b;AGA&#x200b;TAC&#x200b;A<bold>C</bold>GAG&#x200b;CAA&#x200b;TGG&#x200b;C</td>
</tr>
<tr>
<td align="center">PREG006</td>
<td align="center">wild-type</td>
<td align="center">piebald/&#x2b;</td>
<td align="center">CACAGATACA<bold>Y</bold>GAGCAATGGC</td>
</tr>
<tr>
<td align="center">PREG008</td>
<td align="center">piebald</td>
<td align="center">piebald/piebald</td>
<td align="center">CAC&#x200b;AGA&#x200b;TAC&#x200b;A<bold>T</bold>GAG&#x200b;CAA&#x200b;TGG&#x200b;C</td>
</tr>
<tr>
<td align="center">PREG025</td>
<td align="center">stripe</td>
<td align="center">&#x2b;/&#x2b;</td>
<td align="center">CAC&#x200b;AGA&#x200b;TAC&#x200b;A<bold>C</bold>GAG&#x200b;CAA&#x200b;TGG&#x200b;C</td>
</tr>
<tr>
<td align="center">PREG037.24</td>
<td align="center">wild-type</td>
<td align="center">piebald/&#x2b;</td>
<td align="center">CACAGATACA<bold>Y</bold>GAGCAATGGC</td>
</tr>
<tr>
<td align="center">PREG043</td>
<td align="center">wild-type</td>
<td align="center">&#x2b;/&#x2b;</td>
<td align="center">CAC&#x200b;AGA&#x200b;TAC&#x200b;A<bold>C</bold>GAG&#x200b;CAA&#x200b;TGG&#x200b;C</td>
</tr>
<tr>
<td align="center">PREG052</td>
<td align="center">piebald</td>
<td align="center">piebald/piebald</td>
<td align="center">CAC&#x200b;AGA&#x200b;TAC&#x200b;A<bold>T</bold>GAG&#x200b;CAA&#x200b;TGG&#x200b;C</td>
</tr>
</tbody>
</table>
</table-wrap>
</sec>
</sec>
<sec sec-type="results" id="s3">
<title>3 Results</title>
<p>Garcia-Elfing <italic>et al.</italic> (<xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>) provide convincing evidence, from their genomic mapping analyses, that the mutation responsible for the piebald phenotype in ball pythons resides in an 8-Mb interval on the Chromosome 7 of the Burmese python genome. Among the variants detected within the protein-coding sequence of 32 genes in the interval, only one is expected to affect the structure of the corresponding protein. Indeed, this variant introduces a STOP codon in exon 5 of <italic>TFEC,</italic> and the transcription factor produced is truncated. A single copy of <italic>TFEC</italic> is present in reptilian genomes. Based on Sanger sequencing of exon 5 from 7 animals (3 wild-type, 2 homozygous piebald and 2 heterozygous piebald), I was able to confirm the co-segregation of the C-to-T transition in <italic>TFEC</italic> with the piebald allele (<xref ref-type="table" rid="T1">Table 1</xref>) in individuals that were not included in the original study.</p>
<p>I proceeded to characterize histologically the chromatophore composition of the ball python skin both in the wild-type and in piebald mutants. When coaxial light illuminates deparaffinised sections of the black skin of a wild-type ball python, the characteristic reflective shine of iridophores is absent (<xref ref-type="fig" rid="F1">Figure 1A</xref>). The same result is obtained when illuminating sections of white skin patches (<xref ref-type="fig" rid="F1">Figure 1B</xref>) from a piebald python. Similarly, no iridophores were detected in the iris of the wild-type ball python (<xref ref-type="fig" rid="F1">Figure 1C</xref>), where iridophores have previously been reported in other reptiles, such as the Texas rat snake (<xref ref-type="bibr" rid="B17">Ullate-Agote and Tzika, 2021</xref>). This finding must be contrasted to the result obtained when performing the same experiment with skin and eye samples of the corn snake (<italic>Pantherophis guttatus</italic>): the bright reflection of light by iridophores is visible both in the wild-type skin (<xref ref-type="fig" rid="F1">Figure 1C</xref>) and in the skin of the recessive Palmetto phenotype (<xref ref-type="fig" rid="F1">Figure 1D</xref>), as well the eyes of the wild-type corn (<xref ref-type="fig" rid="F1">Figure 1F</xref>). The Palmetto corn snake morph is characterized by fully white dorsal and ventral skin with scarce patches of coloration; the causative mutation for this phenotype remains unknown. Note that the full-depth skin pieces for paraffin sections were large, so any iridophores present could not have been missed.</p>
<fig id="F1" position="float">
<label>FIGURE 1</label>
<caption>
<p>Photography of animals (left) and microscopy image of a deparaffinized dorsal skin section (right) of <bold>(A)</bold> a wild-type adult ball python (black skin), <bold>(B)</bold> a piebald ball python (white patch), <bold>(D)</bold> a wild-type corn snake (dorsal side), and <bold>(E)</bold> a Palmetto corn snake (white patch). Photography of an eye (left) and microscopy image of a deparaffinized eye section from a wild-type ball python <bold>(C)</bold> and a wild-type corn snake <bold>(F)</bold>. The black bar and arrows indicate the melanophores in <bold>(A)</bold>. Cyan bars highlight the layer of iridophores in D, E, and F. Scale bars: 100&#xa0;&#x3bc;m.</p>
</caption>
<graphic xlink:href="fcell-12-1358828-g001.tif"/>
</fig>
<p>Transmission electron microscopy (TEM) imaging of the skin of a wild-type adult ball python shows that i) epidermal and dermal melanophores are present in the dorsal black skin (<xref ref-type="fig" rid="F2">Figure 2A</xref>), ii) epidermal melanophores and xanthophores can be found in the dorsal brown skin (<xref ref-type="fig" rid="F2">Figure 2B</xref>), and iii) the lateral light brown skin (<xref ref-type="fig" rid="F2">Figure 2C</xref>) has a similar composition but the epidermal melanophores seem more scarce. The dermal and epidermal melanophores present similar subcellular morphology (<xref ref-type="fig" rid="F2">Figure 2D</xref>). The presence of concentric lamellae in the xanthosomes (<xref ref-type="fig" rid="F2">Figure 2E</xref>) suggests that they contain pteridines, and they could thus be pterinosomes, but chemical analyses are necessary to confirm this. Here, I use the terms xanthophores and xanthosomes as more generic, rather than to specify the content (pteridines vs carotenoids) of these chromatophores. All pieces of skin contain large amounts of collagen fibers in various orientations (<xref ref-type="fig" rid="F2">Figure 2F</xref>). The subcellular structure of the ball python melanophores and xanthophores strongly resembles that of other snakes, such as the corn snake [Figure 4 in <xref ref-type="bibr" rid="B16">Ullate-Agote et al. (2020)</xref>] and the Texas rat snake [Figure 2 in <xref ref-type="bibr" rid="B17">Ullate-Agote and Tzika (2021)</xref>].</p>
<fig id="F2" position="float">
<label>FIGURE 2</label>
<caption>
<p>TEM imaging of the dorsal black <bold>(A)</bold>, dorsal brown <bold>(B)</bold> and lateral light brown <bold>(C)</bold> skin of an adult wild-type ball python. Magnification of the dermal and epidermal melanophores from A <bold>(D)</bold>, a xanthophore from B <bold>(E)</bold>, and the collagen fibers from C <bold>(F)</bold>. TEM imaging of the dorsal brown <bold>(G)</bold>, lateral white <bold>(H)</bold>, and ventral white <bold>(I)</bold> skin of a piebald adult ball python. The dashed lines mark the dermis/epidermis boundary. eM: epidermal melanophore, dM: dermal melanophore, X: xanthophore, cf: collagen fibers. Scale bars: 5&#xa0;&#x3bc;m <bold>(A,B,C,G,H, I)</bold>, 1&#xa0;&#x3bc;m <bold>(D, F)</bold>, 0.5&#xa0;&#x3bc;m <bold>(E)</bold>.</p>
</caption>
<graphic xlink:href="fcell-12-1358828-g002.tif"/>
</fig>
<p>The dorsal brown skin of a piebald individual (<xref ref-type="fig" rid="F1">Figure 1G</xref>) contains epidermal melanophores and xanthophores. In the white lateral patches (<xref ref-type="fig" rid="F1">Figure 1H</xref>) and the white ventral skin (<xref ref-type="fig" rid="F1">Figure 1I</xref>), only dense collagen fibers can be seen. Thus, the white color we perceive in ball pythons, both on the ventral skin of wild-type animals and on the white dorso-lateral patches of piebald mutants, is likely caused by light scattering from the dense network of collagen fibers, rather than by disorganized lattices of guanine nano-crystals in iridophores. Conversely, TEM imaging of the skin on <italic>TFEC</italic> brown anole mutants (<xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>) showed that they maintain both melanophores and xanthophores, but lose iridophores, explaining why the entire skin of these animals is translucent. Note that, although their coloration is reduced, it is unclear if the pattern itself (<italic>i.e.</italic>, the spatial distribution of colored motifs) is affected in these gene-edited lizards. In conclusion, the data presented here show that the phenotype of <italic>TFEC</italic> mutants of the brown anole does not recapitulate the effect of the <italic>TFEC</italic> mutation in piebald ball pythons, advocating that extrapolations of results between different species can be illuminating but must be made with caution.</p>
</sec>
<sec sec-type="discussion" id="s4">
<title>4 Discussion</title>
<p>TFEC is part of the microphthalmia/transcription factor E (MiT/TFE) family of transcription factors, which additionally includes MITF (Melanocyte Inducing Transcription Factor), TFEB (Transcription Factor EB), and TFE3 (Transcription Factor Binding to IGHM Enhancer 3). So far, TFEB and TFE3 have not been associated with the development of animal skin coloration. In mice, a <italic>Tfec</italic> knockout presents abnormal hair pigmentation (<xref ref-type="bibr" rid="B4">Groza et al., 2023</xref>), whereas mice with a partially truncated <italic>Tfec,</italic> as in the piebald pythons, are normally pigmented (<xref ref-type="bibr" rid="B14">Steingrimsson et al., 2002</xref>). <italic>Mitf</italic> mutations in mice result in decreased or absent pigmentation and, occasionally, in white spotting (<xref ref-type="bibr" rid="B13">Steingrimsson et al., 2004</xref>). In zebrafish, <italic>mitfa</italic> mutants are characterized by the loss of melanophores and an increased number of iridophores (<xref ref-type="bibr" rid="B8">Lister et al., 1999</xref>), and <italic>tfec</italic> mutations mainly impact the differentiation of iridophores (<xref ref-type="bibr" rid="B11">Petratou et al., 2021</xref>). Reduced <italic>MITF</italic> expression in the Texas Rat snake results in the loss of melanophores and xanthophores (<xref ref-type="bibr" rid="B17">Ullate-Agote and Tzika, 2021</xref>). If we focus on the phenotype of gene knockout animals (induced or spontaneously-occurring), rather than on the proposed models describing the fate determination processes of chromatophores [progressive, direct, and cyclical (<xref ref-type="bibr" rid="B6">Kelsh et al., 2021</xref>)], we can speculate that the role of MITF and TFEC varies in different vertebrate lineages. Possible associations based on the scarce existing data on reptiles are presented in the simplified scheme of <xref ref-type="fig" rid="F3">Figure 3</xref>. Depending on the chromatophore types present in a species, these transcription factors are likely to take up different functions in the chromatophore fate determination.</p>
<fig id="F3" position="float">
<label>FIGURE 3</label>
<caption>
<p>Schematic representation of the possible association of MITF and TFEC in the differentiation of the three chromatophore types from neural crest precursor cells in zebrafish, snakes and lizards. The fate determination processes are disregarded in this simplified scheme.</p>
</caption>
<graphic xlink:href="fcell-12-1358828-g003.tif"/>
</fig>
<p>To elucidate the actual role of these transcription factors in reptilian coloration, it is thus necessary to investigate further their evolutionary history across lineages, for example, by comparing their protein structure and expression levels during development. In this study, I only sampled adult individuals, so I cannot exclude the possibility that iridophores are present in ball pythons during embryogenesis and disappear as the animals grow. Indeed, there are reptiles, like the ocellated and other lizards (<xref ref-type="bibr" rid="B2">Fofonjka and Milinkovitch, 2021</xref>; <xref ref-type="bibr" rid="B5">Jahanbakhsh and Milinkovitch, 2022</xref>), whose skin coloration and pattern changes continuously. This is not the case though for ball pythons; they maintain the same skin coloration and pattern throughout their lives. Nevertheless, it is necessary to investigate the differentiation of their chromatophores during development. Transcriptomic analyses can help us identify gene markers of reptilian chromatophores, such that we track them during embryogenesis. Undoubtedly, gene-editing experiments in multiple species, as performed in the brown anole (<xref ref-type="bibr" rid="B3">Garcia-Elfring et al., 2023</xref>), would also illuminate how differentiation, maturation, migration, and survival of different chromatophores&#x2014;as well as their interactions&#x2014;have evolved to produce the remarkable diversity of colors and patterns in Squamate reptiles.</p>
</sec>
</body>
<back>
<sec sec-type="data-availability" id="s5">
<title>Data availability statement</title>
<p>The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.</p>
</sec>
<sec id="s6">
<title>Ethics statement</title>
<p>The animal study was approved by the Geneva veterinary cantonal authorities. The study was conducted in accordance with the local legislation and institutional requirements.</p>
</sec>
<sec id="s7">
<title>Author contributions</title>
<p>ACT: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Validation, Visualization, Writing&#x2013;original draft, Writing&#x2013;review and editing.</p>
</sec>
<sec sec-type="funding-information" id="s8">
<title>Funding</title>
<p>The author declares financial support was received for the research, authorship, and/or publication of this article. Funding was provided by an SNSF grant (310030_204466), the Ernst and Lucie Schmidheiny Foundation (10_2023), the Fonds G&#xe9;n&#xe9;ral de l&#x2019;Universit&#xe9; de Gen&#xe8;ve (23_28), and the Emile Plantamour Fund (2024/12 - Soci&#xe9;t&#xe9; Acad&#xe9;mique de Gen&#xe8;ve) to ACT.</p>
</sec>
<ack>
<p>I would like to thank Ingrid Burgelin, Carine Langrez and Valerie Haechler for technical assistance, as well as Adrien Debry and Florent Montange for animal maintenance and imaging.</p>
</ack>
<sec sec-type="COI-statement" id="s9">
<title>Conflict of interest</title>
<p>The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.</p>
</sec>
<sec sec-type="disclaimer" id="s10">
<title>Publisher&#x2019;s note</title>
<p>All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.</p>
</sec>
<ref-list>
<title>References</title>
<ref id="B1">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Andrade</surname>
<given-names>P.</given-names>
</name>
<name>
<surname>Pinho</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Perez</surname>
<given-names>I. D. L. G.</given-names>
</name>
<name>
<surname>Afonso</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Brejcha</surname>
<given-names>J.</given-names>
</name>
<name>
<surname>Rubin</surname>
<given-names>C. J.</given-names>
</name>
<etal/>
</person-group> (<year>2019</year>). <article-title>Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>116</volume>, <fpage>5633</fpage>&#x2013;<lpage>5642</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.1820320116</pub-id>
</citation>
</ref>
<ref id="B2">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Fofonjka</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Milinkovitch</surname>
<given-names>M. C.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Reaction-diffusion in a growing 3D domain of skin scales generates a discrete cellular automaton</article-title>. <source>Nat. Commun.</source> <volume>12</volume>, <fpage>2433</fpage>. <pub-id pub-id-type="doi">10.1038/s41467-021-22525-1</pub-id>
</citation>
</ref>
<ref id="B3">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Garcia-Elfring</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Sabin</surname>
<given-names>C. E.</given-names>
</name>
<name>
<surname>Iouchmanov</surname>
<given-names>A. L.</given-names>
</name>
<name>
<surname>Roffey</surname>
<given-names>H. L.</given-names>
</name>
<name>
<surname>Samudra</surname>
<given-names>S. P.</given-names>
</name>
<name>
<surname>Alcala</surname>
<given-names>A. J.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>Piebaldism and chromatophore development in reptiles are linked to the tfec gene</article-title>. <source>Curr. Biol.</source> <volume>33</volume>, <fpage>755</fpage>&#x2013;<lpage>763 e3</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2023.01.004</pub-id>
</citation>
</ref>
<ref id="B4">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Groza</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Gomez</surname>
<given-names>F. L.</given-names>
</name>
<name>
<surname>Mashhadi</surname>
<given-names>H. H.</given-names>
</name>
<name>
<surname>Munoz-Fuentes</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Gunes</surname>
<given-names>O.</given-names>
</name>
<name>
<surname>Wilson</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>2023</year>). <article-title>The International Mouse Phenotyping Consortium: comprehensive knockout phenotyping underpinning the study of human disease</article-title>. <source>Nucleic Acids Res.</source> <volume>51</volume>, <fpage>D1038</fpage>&#x2013;<lpage>D1045</lpage>. <pub-id pub-id-type="doi">10.1093/nar/gkac972</pub-id>
</citation>
</ref>
<ref id="B5">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Jahanbakhsh</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Milinkovitch</surname>
<given-names>M. C.</given-names>
</name>
</person-group> (<year>2022</year>). <article-title>Modeling convergent scale-by-scale skin color patterning in multiple species of lizards</article-title>. <source>Curr. Biol.</source> <volume>32</volume>, <fpage>5069</fpage>&#x2013;<lpage>5082</lpage>. <pub-id pub-id-type="doi">10.1016/j.cub.2022.10.044</pub-id>
</citation>
</ref>
<ref id="B6">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kelsh</surname>
<given-names>R. N.</given-names>
</name>
<name>
<surname>Sosa</surname>
<given-names>K. C.</given-names>
</name>
<name>
<surname>Farjami</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Makeev</surname>
<given-names>V.</given-names>
</name>
<name>
<surname>Dawes</surname>
<given-names>J. H. P.</given-names>
</name>
<name>
<surname>Rocco</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Cyclical fate restriction: a new view of neural crest cell fate specification</article-title>. <source>Development</source> <volume>148</volume>, <comment>dev176057</comment>. <pub-id pub-id-type="doi">10.1242/dev.176057</pub-id>
</citation>
</ref>
<ref id="B7">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Kuriyama</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Murakami</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Brandley</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Hasegawa</surname>
<given-names>M.</given-names>
</name>
</person-group> (<year>2020</year>). <article-title>Blue, black, and stripes: evolution and development of color production and pattern formation in lizards and snakes</article-title>. <source>Front. Ecol. Evol.</source> <volume>8</volume>. <pub-id pub-id-type="doi">10.3389/fevo.2020.00232</pub-id>
</citation>
</ref>
<ref id="B8">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Lister</surname>
<given-names>J. A.</given-names>
</name>
<name>
<surname>Robertson</surname>
<given-names>C. P.</given-names>
</name>
<name>
<surname>Lepage</surname>
<given-names>T.</given-names>
</name>
<name>
<surname>Johnson</surname>
<given-names>S. L.</given-names>
</name>
<name>
<surname>Raible</surname>
<given-names>D. W.</given-names>
</name>
</person-group> (<year>1999</year>). <article-title>Nacre encodes a zebrafish microphthalmia-related protein that regulates neural-crest-derived pigment cell fate</article-title>. <source>Development</source> <volume>126</volume>, <fpage>3757</fpage>&#x2013;<lpage>3767</lpage>. <pub-id pub-id-type="doi">10.1242/dev.126.17.3757</pub-id>
</citation>
</ref>
<ref id="B9">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Milinkovitch</surname>
<given-names>M. C.</given-names>
</name>
<name>
<surname>Tzika</surname>
<given-names>A.</given-names>
</name>
</person-group> (<year>2007</year>). <article-title>Escaping the mouse trap: the selection of new evo-devo model species</article-title>. <source>J. Exp. Zoology Part B-Molecular Dev. Evol.</source> <volume>308b</volume>, <fpage>337</fpage>&#x2013;<lpage>346</lpage>. <pub-id pub-id-type="doi">10.1002/jez.b.21180</pub-id>
</citation>
</ref>
<ref id="B10">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Patterson</surname>
<given-names>L. B.</given-names>
</name>
<name>
<surname>Parichy</surname>
<given-names>D. M.</given-names>
</name>
</person-group> (<year>2019</year>). <article-title>Zebrafish pigment pattern formation: insights into the development and evolution of adult form</article-title>. <source>Annu. Rev. Genet.</source> <volume>53</volume> (<issue>53</issue>), <fpage>505</fpage>&#x2013;<lpage>530</lpage>. <pub-id pub-id-type="doi">10.1146/annurev-genet-112618-043741</pub-id>
</citation>
</ref>
<ref id="B11">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Petratou</surname>
<given-names>K.</given-names>
</name>
<name>
<surname>Spencer</surname>
<given-names>S. A.</given-names>
</name>
<name>
<surname>Kelsh</surname>
<given-names>R. N.</given-names>
</name>
<name>
<surname>Lister</surname>
<given-names>J. A.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>The MITF paralog tfec is required in neural crest development for fate specification of the iridophore lineage from a multipotent pigment cell progenitor</article-title>. <source>Plos One</source> <volume>16</volume>, <comment>e0244794</comment>. <pub-id pub-id-type="doi">10.1371/journal.pone.0244794</pub-id>
</citation>
</ref>
<ref id="B12">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Rasys</surname>
<given-names>A. M.</given-names>
</name>
<name>
<surname>Park</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Ball</surname>
<given-names>R. E.</given-names>
</name>
<name>
<surname>Alcala</surname>
<given-names>A. J.</given-names>
</name>
<name>
<surname>Lauderdale</surname>
<given-names>J. D.</given-names>
</name>
<name>
<surname>Menke</surname>
<given-names>D. B.</given-names>
</name>
</person-group> <year>2019</year>. <article-title>CRISPR-Cas9 gene editing in lizards through microinjection of unfertilized oocytes</article-title>. <source>Cell Rep.</source>, <volume>28</volume>, <fpage>2288</fpage>&#x2013;<lpage>2292</lpage>. <pub-id pub-id-type="doi">10.1016/j.celrep.2019.07.089</pub-id>
</citation>
</ref>
<ref id="B13">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Steingrimsson</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Copeland</surname>
<given-names>N. G.</given-names>
</name>
<name>
<surname>Jenkins</surname>
<given-names>N. A.</given-names>
</name>
</person-group> (<year>2004</year>). <article-title>Melanocytes and the Microphthalmia transcription factor network</article-title>. <source>Annu. Rev. Genet.</source> <volume>38</volume>, <fpage>365</fpage>&#x2013;<lpage>411</lpage>. <pub-id pub-id-type="doi">10.1146/annurev.genet.38.072902.092717</pub-id>
</citation>
</ref>
<ref id="B14">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Steingrimsson</surname>
<given-names>E.</given-names>
</name>
<name>
<surname>Tessarollo</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Pathak</surname>
<given-names>B.</given-names>
</name>
<name>
<surname>Hou</surname>
<given-names>L.</given-names>
</name>
<name>
<surname>Arnheiter</surname>
<given-names>H.</given-names>
</name>
<name>
<surname>Copeland</surname>
<given-names>N. G.</given-names>
</name>
<etal/>
</person-group> (<year>2002</year>). <article-title>Mitf and Tfe3, two members of the Mitf-Tfe family of bHLH-Zip transcription factors, have important but functionally redundant roles in osteoclast development</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>99</volume>, <fpage>4477</fpage>&#x2013;<lpage>4482</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.072071099</pub-id>
</citation>
</ref>
<ref id="B15">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Tzika</surname>
<given-names>A. C.</given-names>
</name>
<name>
<surname>Ullate-Agote</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Zakany</surname>
<given-names>S.</given-names>
</name>
<name>
<surname>Kummrow</surname>
<given-names>M.</given-names>
</name>
<name>
<surname>Milinkovitch</surname>
<given-names>M. C.</given-names>
</name>
</person-group> (<year>2023</year>). <article-title>Somitic positional information guides self-organized patterning of snake scales</article-title>. <source>Sci. Adv.</source> <volume>9</volume>, <comment>eadf8834</comment>. <pub-id pub-id-type="doi">10.1126/sciadv.adf8834</pub-id>
</citation>
</ref>
<ref id="B16">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ullate-Agote</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Burgelin</surname>
<given-names>I.</given-names>
</name>
<name>
<surname>Debry</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Langrez</surname>
<given-names>C.</given-names>
</name>
<name>
<surname>Montange</surname>
<given-names>F.</given-names>
</name>
<name>
<surname>Peraldi</surname>
<given-names>R.</given-names>
</name>
<etal/>
</person-group> (<year>2020</year>). <article-title>Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles</article-title>. <source>Proc. Natl. Acad. Sci. U. S. A.</source> <volume>117</volume>, <fpage>26307</fpage>&#x2013;<lpage>26317</lpage>. <pub-id pub-id-type="doi">10.1073/pnas.2003724117</pub-id>
</citation>
</ref>
<ref id="B17">
<citation citation-type="journal">
<person-group person-group-type="author">
<name>
<surname>Ullate-Agote</surname>
<given-names>A.</given-names>
</name>
<name>
<surname>Tzika</surname>
<given-names>A. C.</given-names>
</name>
</person-group> (<year>2021</year>). <article-title>Characterization of the leucistic Texas rat snake pantherophis obsoletus</article-title>. <source>Front. Ecol. Evol.</source> <volume>9</volume>. <pub-id pub-id-type="doi">10.3389/fevo.2021.583136</pub-id>
</citation>
</ref>
</ref-list>
</back>
</article>