Athanasia C. Tzika- Laboratory of Artificial and Natural Evolution (LANE), Department of Genetics and Evolution, University of Geneva, Geneva, Switzerland
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.
1 Introduction
Although the zebrafish remains the main reference in the animal coloration field (Patterson and Parichy, 2019), reptilian species are gaining ground as new models (Milinkovitch and Tzika, 2007) 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 (Kuriyama et al., 2020), 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 (Andrade et al., 2019; Ullate-Agote et al., 2020; Ullate-Agote and Tzika, 2021; Garcia-Elfring et al., 2023) and the establishment of gene-editing protocols using CRISPR-Cas9 in the same lineages (Rasys et al., 2019; Tzika et al., 2023) 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.
Garcia-Elfing et al. (Garcia-Elfring et al., 2023) recently investigated the piebald ball python morph (Python regius), 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 TFEC 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 TFEC in the brown anole lizard (Anolis sergei), resulting in the loss of iridophores and reduced body coloration, particularly in the snout, forelimbs, and hindlimbs. But as Garcia-Elfing et al. (Garcia-Elfring et al., 2023) conclude, it remains “unresolved whether piebald ball pythons have iridophores in either pigmented or white skin”, 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 et al. and the results presented here suggest that the same molecule can affect reptilian and vertebrate coloration in different ways.
2 Materials and methods
2.1 Animal experimentation
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.
2.2 Histology and imaging
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 cm × 2 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 mm2 were fixed, sectioned, and imaged as previously described (Ullate-Agote et al., 2020). Sample processing and imaging were performed at the Electron Microscopy Facility, University of Lausanne (Switzerland).
2.3 TFEC amplification
For all animals (Table 1), genomic DNA was extracted from whole blood using the DNeasy Blood & Tissue Kit (Qiagen, 51104) and genotyping by Sanger sequencing, targeting the TFEC SNP identified in the original publication (Garcia-Elfring et al., 2023), was performed with the following primers: TFECg_103193F (CAGTGCAACTCAAAGGGAACA) and TFECg_103880R (GCAGACCCATGAAATCAATGGA).
TABLE 1. Individuals genotyped for the presence of the SNP suggested to cause the piebald ball python phenotype. The SNP is highlighted in bold in the ‘sequence’ column. The ‘genotype’ is deduced by the appearance and the pedigree of the animals. ‘Stripe’ is another ball python morph unrelated to the piebald morph.
3 Results
Garcia-Elfing et al. (Garcia-Elfring et al., 2023) 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 TFEC, and the transcription factor produced is truncated. A single copy of TFEC 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 TFEC with the piebald allele (Table 1) in individuals that were not included in the original study.
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 (Figure 1A). The same result is obtained when illuminating sections of white skin patches (Figure 1B) from a piebald python. Similarly, no iridophores were detected in the iris of the wild-type ball python (Figure 1C), where iridophores have previously been reported in other reptiles, such as the Texas rat snake (Ullate-Agote and Tzika, 2021). This finding must be contrasted to the result obtained when performing the same experiment with skin and eye samples of the corn snake (Pantherophis guttatus): the bright reflection of light by iridophores is visible both in the wild-type skin (Figure 1C) and in the skin of the recessive Palmetto phenotype (Figure 1D), as well the eyes of the wild-type corn (Figure 1F). 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.
FIGURE 1. Photography of animals (left) and microscopy image of a deparaffinized dorsal skin section (right) of (A) a wild-type adult ball python (black skin), (B) a piebald ball python (white patch), (D) a wild-type corn snake (dorsal side), and (E) 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 (C) and a wild-type corn snake (F). The black bar and arrows indicate the melanophores in (A). Cyan bars highlight the layer of iridophores in D, E, and F. Scale bars: 100 μm.
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 (Figure 2A), ii) epidermal melanophores and xanthophores can be found in the dorsal brown skin (Figure 2B), and iii) the lateral light brown skin (Figure 2C) has a similar composition but the epidermal melanophores seem more scarce. The dermal and epidermal melanophores present similar subcellular morphology (Figure 2D). The presence of concentric lamellae in the xanthosomes (Figure 2E) 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 (Figure 2F). The subcellular structure of the ball python melanophores and xanthophores strongly resembles that of other snakes, such as the corn snake [Figure 4 in Ullate-Agote et al. (2020)] and the Texas rat snake [Figure 2 in Ullate-Agote and Tzika (2021)].
FIGURE 2. TEM imaging of the dorsal black (A), dorsal brown (B) and lateral light brown (C) skin of an adult wild-type ball python. Magnification of the dermal and epidermal melanophores from A (D), a xanthophore from B (E), and the collagen fibers from C (F). TEM imaging of the dorsal brown (G), lateral white (H), and ventral white (I) 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 μm (A,B,C,G,H, I), 1 μm (D, F), 0.5 μm (E).
The dorsal brown skin of a piebald individual (Figure 1G) contains epidermal melanophores and xanthophores. In the white lateral patches (Figure 1H) and the white ventral skin (Figure 1I), 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 TFEC brown anole mutants (Garcia-Elfring et al., 2023) 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 (i.e., the spatial distribution of colored motifs) is affected in these gene-edited lizards. In conclusion, the data presented here show that the phenotype of TFEC mutants of the brown anole does not recapitulate the effect of the TFEC mutation in piebald ball pythons, advocating that extrapolations of results between different species can be illuminating but must be made with caution.
4 Discussion
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 Tfec knockout presents abnormal hair pigmentation (Groza et al., 2023), whereas mice with a partially truncated Tfec, as in the piebald pythons, are normally pigmented (Steingrimsson et al., 2002). Mitf mutations in mice result in decreased or absent pigmentation and, occasionally, in white spotting (Steingrimsson et al., 2004). In zebrafish, mitfa mutants are characterized by the loss of melanophores and an increased number of iridophores (Lister et al., 1999), and tfec mutations mainly impact the differentiation of iridophores (Petratou et al., 2021). Reduced MITF expression in the Texas Rat snake results in the loss of melanophores and xanthophores (Ullate-Agote and Tzika, 2021). 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 (Kelsh et al., 2021)], 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 Figure 3. Depending on the chromatophore types present in a species, these transcription factors are likely to take up different functions in the chromatophore fate determination.
FIGURE 3. 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.
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 (Fofonjka and Milinkovitch, 2021; Jahanbakhsh and Milinkovitch, 2022), 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 (Garcia-Elfring et al., 2023), would also illuminate how differentiation, maturation, migration, and survival of different chromatophores—as well as their interactions—have evolved to produce the remarkable diversity of colors and patterns in Squamate reptiles.
Data availability statement
The original contributions presented in the study are included in the article, further inquiries can be directed to the corresponding author.
Ethics statement
The animal study was approved by the Geneva veterinary cantonal authorities. The study was conducted in accordance with the local legislation and institutional requirements.
Author contributions
ACT: Conceptualization, Formal Analysis, Funding acquisition, Investigation, Methodology, Project administration, Validation, Visualization, Writing–original draft, Writing–review and editing.
Funding
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énéral de l’Université de Genève (23_28), and the Emile Plantamour Fund (2024/12 - Société Académique de Genève) to ACT.
Acknowledgments
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.
Conflict of interest
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.
Publisher’s note
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.
References
Andrade, P., Pinho, C., Perez, I. D. L. G., Afonso, S., Brejcha, J., Rubin, C. J., et al. (2019). Regulatory changes in pterin and carotenoid genes underlie balanced color polymorphisms in the wall lizard. Proc. Natl. Acad. Sci. U. S. A. 116, 5633–5642. doi:10.1073/pnas.1820320116
Fofonjka, A., and Milinkovitch, M. C. (2021). Reaction-diffusion in a growing 3D domain of skin scales generates a discrete cellular automaton. Nat. Commun. 12, 2433. doi:10.1038/s41467-021-22525-1
Garcia-Elfring, A., Sabin, C. E., Iouchmanov, A. L., Roffey, H. L., Samudra, S. P., Alcala, A. J., et al. (2023). Piebaldism and chromatophore development in reptiles are linked to the tfec gene. Curr. Biol. 33, 755–763 e3. doi:10.1016/j.cub.2023.01.004
Groza, T., Gomez, F. L., Mashhadi, H. H., Munoz-Fuentes, V., Gunes, O., Wilson, R., et al. (2023). The International Mouse Phenotyping Consortium: comprehensive knockout phenotyping underpinning the study of human disease. Nucleic Acids Res. 51, D1038–D1045. doi:10.1093/nar/gkac972
Jahanbakhsh, E., and Milinkovitch, M. C. (2022). Modeling convergent scale-by-scale skin color patterning in multiple species of lizards. Curr. Biol. 32, 5069–5082. doi:10.1016/j.cub.2022.10.044
Kelsh, R. N., Sosa, K. C., Farjami, S., Makeev, V., Dawes, J. H. P., and Rocco, A. (2021). Cyclical fate restriction: a new view of neural crest cell fate specification. Development 148, dev176057. doi:10.1242/dev.176057
Kuriyama, T., Murakami, A., Brandley, M., and Hasegawa, M. (2020). Blue, black, and stripes: evolution and development of color production and pattern formation in lizards and snakes. Front. Ecol. Evol. 8. doi:10.3389/fevo.2020.00232
Lister, J. A., Robertson, C. P., Lepage, T., Johnson, S. L., and Raible, D. W. (1999). Nacre encodes a zebrafish microphthalmia-related protein that regulates neural-crest-derived pigment cell fate. Development 126, 3757–3767. doi:10.1242/dev.126.17.3757
Milinkovitch, M. C., and Tzika, A. (2007). Escaping the mouse trap: the selection of new evo-devo model species. J. Exp. Zoology Part B-Molecular Dev. Evol. 308b, 337–346. doi:10.1002/jez.b.21180
Patterson, L. B., and Parichy, D. M. (2019). Zebrafish pigment pattern formation: insights into the development and evolution of adult form. Annu. Rev. Genet. 53 (53), 505–530. doi:10.1146/annurev-genet-112618-043741
Petratou, K., Spencer, S. A., Kelsh, R. N., and Lister, J. A. (2021). The MITF paralog tfec is required in neural crest development for fate specification of the iridophore lineage from a multipotent pigment cell progenitor. Plos One 16, e0244794. doi:10.1371/journal.pone.0244794
Rasys, A. M., Park, S., Ball, R. E., Alcala, A. J., Lauderdale, J. D., and Menke, D. B. 2019. CRISPR-Cas9 gene editing in lizards through microinjection of unfertilized oocytes. Cell Rep., 28, 2288–2292. doi:10.1016/j.celrep.2019.07.089
Steingrimsson, E., Copeland, N. G., and Jenkins, N. A. (2004). Melanocytes and the Microphthalmia transcription factor network. Annu. Rev. Genet. 38, 365–411. doi:10.1146/annurev.genet.38.072902.092717
Steingrimsson, E., Tessarollo, L., Pathak, B., Hou, L., Arnheiter, H., Copeland, N. G., et al. (2002). Mitf and Tfe3, two members of the Mitf-Tfe family of bHLH-Zip transcription factors, have important but functionally redundant roles in osteoclast development. Proc. Natl. Acad. Sci. U. S. A. 99, 4477–4482. doi:10.1073/pnas.072071099
Tzika, A. C., Ullate-Agote, A., Zakany, S., Kummrow, M., and Milinkovitch, M. C. (2023). Somitic positional information guides self-organized patterning of snake scales. Sci. Adv. 9, eadf8834. doi:10.1126/sciadv.adf8834
Ullate-Agote, A., Burgelin, I., Debry, A., Langrez, C., Montange, F., Peraldi, R., et al. (2020). Genome mapping of a LYST mutation in corn snakes indicates that vertebrate chromatophore vesicles are lysosome-related organelles. Proc. Natl. Acad. Sci. U. S. A. 117, 26307–26317. doi:10.1073/pnas.2003724117
Keywords: reptiles, coloration, skin, chromatophores, ball python, snakes, TFEC, anole lizard
Citation: Tzika AC (2024) On the role of TFEC in reptilian coloration. Front. Cell Dev. Biol. 12:1358828. doi: 10.3389/fcell.2024.1358828
Received: 20 December 2023; Accepted: 22 January 2024;
Published: 07 February 2024.
Edited by:
Alejandra García-Gasca, National Council of Science and Technology (CONACYT), MexicoReviewed by:
Robert Neil Kelsh, University of Bath, United KingdomDouglas Menke, University of Georgia, United States
Copyright © 2024 Tzika. 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.
*Correspondence: Athanasia C. Tzika, athanasia.tzika@unige.ch