In Silico Analysis of Gene Expression Network Components Underlying Pigmentation Phenotypes in the Python Identified Evolutionarily Conserved Clusters of Transcription Factor Binding Sites

A Ball Python Colour Morph Implicates MC1R in Melanophore-Xanthophore Distribution and Pattern Formation

Reptiles showcase an extensive array of skin colours and patterns, yet little is known about the genetics of reptile colouration. Here, we investigate the genetic basis of the Clown colour morph found in captive-bred ball pythons (Python regius) to study skin pigmentation and patterning in snakes. We obtained samples by crowdsourcing shed skin from commercial breeders and hobbyists. We applied a case-control design, whole-genome pool sequencing, variant annotation, histological analyses, and electron microscopy imaging. We identified a missense mutation in a transmembrane region of the melanocortin-1 receptor (MC1R) associated with the Clown phenotype. In classic avian and mammalian model species, MC1R is known for controlling the type and amount of melanin produced. In contrast, our results suggest that MC1R signalling might play a key role in pattern formation in ball pythons, affecting xanthophore-melanophore distribution. This work highlights the varied functions of MC1R across different vertebrate lineages and promotes a novel model system to study reptile colouration.

Piebaldism and chromatophore development in reptiles are linked to the tfec gene

Reptiles display great diversity in color and pattern, yet much of what we know about vertebrate coloration comes from classic model species such as the mouse and zebrafish.1,2,3,4 Captive-bred ball pythons (Python regius) exhibit a remarkable degree of color and pattern variation. Despite the wide range of Mendelian color phenotypes available in the pet trade, ball pythons remain an overlooked species in pigmentation research. Here, we investigate the genetic basis of the recessive piebald phenotype, a pattern defect characterized by patches of unpigmented skin (leucoderma). We performed whole-genome sequencing and used a case-control approach to discover a nonsense mutation in the gene encoding the transcription factor tfec, implicating this gene in the leucodermic patches in ball pythons. We functionally validated tfec in a lizard model (Anolis sagrei) using the gene editing CRISPR/Cas9 system and TEM imaging of skin. Our findings show that reading frame mutations in tfec affect coloration and lead to a loss of iridophores in Anolis, indicating that tfec is required for chromatophore development. This study highlights the value of captive-bred ball pythons as a model species for accelerating discoveries on the genetic basis of vertebrate coloration.

A community-science approach identifies genetic variants associated with three color morphs in ball pythons (Python regius)

Color morphs in ball pythons (Python regius) provide a unique and largely untapped resource for understanding the genetics of coloration in reptiles. Here we use a community-science approach to investigate the genetics of three color morphs affecting production of the pigment melanin. These morphs-Albino, Lavender Albino, and Ultramel-show a loss of melanin in the skin and eyes, ranging from severe (Albino) to moderate (Lavender Albino) to mild (Ultramel). To identify genetic variants causing each morph, we recruited shed skins of pet ball pythons via social media, extracted DNA from the skins, and searched for putative loss-of-function variants in homologs of genes controlling melanin production in other vertebrates. We report that the Albino morph is associated with missense and non-coding variants in the gene TYR. The Lavender Albino morph is associated with a deletion in the gene OCA2. The Ultramel morph is associated with a missense variant and a putative deletion in the gene TYRP1. Our study is one of the first to identify genetic variants associated with color morphs in ball pythons and shows that pet samples recruited from the community can provide a resource for genetic studies in this species.

Developmental mechanisms of longitudinal stripes in the Japanese four-lined snake

The developmental mechanisms of color patterns formation and its evolution remain unclear in reptilian sauropsids. We, therefore, studied the pigment cell mechanisms of stripe pattern formation during embryonic development of the snake Elaphe quadrivirgata. We identified 10 post-ovipositional embryonic developmental stages based on external morphological characteristics. Examination for the temporal changes in differentiation, distribution, and density of pigment cells during embryonic development revealed that melanophores first appeared in myotome and body cavity but not in skin surface at Stage 5. Epidermal melanophores were first recognized at Stage 7, and dermal melanophores and iridophores appeared in Stage 9. Stripe pattern first appeared to establish at Stage 8 as a spatial density gradient of epidermal melanophores between the regions of future dark brown longitudinal stripes and light colored background. Our study, thus, provides a comprehensive pigment-cell-based understanding of stripe pattern formation during embryonic development. We briefly discuss the importance of the gene expression studies by considering the biologically relevant theoretical models with standard developmental staging for understanding reptilian color pattern evolution.