Hidden invasiveness of non-native Schlegel's Japanese gecko (Reptilia: Squamata: Gekkonidae) and three-way competition among natives and non-natives in Japan.

In this study, we investigated the invasiveness of Gekko japonicus, a prevalent gecko species in Japan and an ancient non-native species, focusing on its competition with both the undescribed endemic Gekko species (referred to as Nishiyamori in Japanese) and G. hokouensis. These species are co-distributed with G. japonicus, leading us to hypothesize that G. japonicus was invasive upon its initial introduction. We employed niche analysis and population genetics through ddRAD-seq to assess the historical invasiveness of G. japonicus by comparing regions with and without interspecies competition. Our niche analysis across the Goto Islands, Hiradojima Island (colonized by G. japonicus) and the Koshikishima Islands (not colonized by G. japonicus) indicated that endemic Gekko sp. alter their microhabitat usage in response to invasions by other gecko species, despite having similar suitable habitats and microhabitat preferences. Population genetic analysis revealed significant population declines in Gekko sp. within areas of introduced competition, in contrast to stable populations in areas without such competition. These findings suggest a tripartite competitive relationship among the gecko species, with G. japonicus and G. hokouensis invasions restricting the distribution of the endemic Gekko sp. Consequently, G. japonicus may have historically acted as an invasive species. Acknowledging the historical dynamics of current biodiversity is crucial for addressing complex ecological issues and making informed conservation decisions.

Post-meiotic mechanism of facultative parthenogenesis in gonochoristic whiptail lizard species.

Facultative parthenogenesis (FP) has historically been regarded as rare in vertebrates, but in recent years incidences have been reported in a growing list of fish, reptile, and bird species. Despite the increasing interest in the phenomenon, the underlying mechanism and evolutionary implications have remained unclear. A common finding across many incidences of FP is either a high degree of homozygosity at microsatellite loci or low levels of heterozygosity detected in next-generation sequencing data. This has led to the proposal that second polar body fusion following the meiotic divisions restores diploidy and thereby mimics fertilization. Here, we show that FP occurring in the gonochoristic Aspidoscelis species A. marmoratus and A. arizonae results in genome-wide homozygosity, an observation inconsistent with polar body fusion as the underlying mechanism of restoration. Instead, a high-quality reference genome for A. marmoratus and analysis of whole-genome sequencing from multiple FP and control animals reveals that a post-meiotic mechanism gives rise to homozygous animals from haploid, unfertilized oocytes. Contrary to the widely held belief that females need to be isolated from males to undergo FP, females housed with conspecific and heterospecific males produced unfertilized eggs that underwent spontaneous development. In addition, offspring arising from both fertilized eggs and parthenogenetic development were observed to arise from a single clutch. Strikingly, our data support a mechanism for facultative parthenogenesis that removes all heterozygosity in a single generation. Complete homozygosity exposes the genetic load and explains the high rate of congenital malformations and embryonic mortality associated with FP in many species. Conversely, for animals that develop normally, FP could potentially exert strong purifying selection as all lethal recessive alleles are purged in a single generation.

Substrate specialisation drives an unexpectedly diverse radiation in barking geckos (Ptenopus: Gekkonidae).

Barking geckos (genus Ptenopus) are terrestrial, burrowing lizards endemic to southern Africa, currently with three recognised species. Two species are range-restricted (P. kochi and P. carpi) and display clear differences in substrate preference (soft sand vs. hard gravel). The third and most widespread species, P. garrulus, occurs on a variety of substrates of differing hardness, across potential geographic barriers, and over a steep climatic gradient. Variations in morphology and advertisement calls indicates that P. garrulus may be a species complex. Two subspecies of P. garrulus are currently recognised: P. g. maculatus and P. g. garrulus. To investigate species boundaries, we produced the first comprehensive phylogeny for the genus. We used a novel application of multiple regression on matrices models to assess multiple environmental drivers of diversification, as contrasted to isolation by distance. We show that P. kochi, P. carpi, and P. g. garrulus are valid species, but that P. g. maculatus is a paraphyletic complex of five previously unrecognised taxa. Specialisation onto different substrates was likely the main driver of divergence, with parapatric occurrence of two to four clades occurring at each of the three substrate transition zones identified a priori. The region encompasses diverse bioclimatic regions and potential geographic barriers, and these likely played a role in some divergence events.

Cryptic diversity and phylogeographic patterns of Mediodactylus species in the Eastern Mediterranean region.

Cryptic diversity poses a great obstacle in our attempts to assess the current biodiversity crisis and may hamper conservation efforts. The gekkonid genus Mediodactylus, a well-known case of hidden species and genetic diversity, has been taxonomically reclassified several times during the last decade. Focusing on the Mediterranean populations, a recent study within the M. kotschyi species complex using classic mtDNA/nuDNA markers suggested the existence of five distinct species, some being endemic and some possibly threatened, yet their relationships have not been fully resolved. Here, we generated genome-wide SNPs (using ddRADseq) and applied molecular species delimitation approaches and population genomic analyses to further disentangle these relationships. Τhe most extensive nuclear dataset, so far, encompassing 2,360 loci and âˆ¼ 699,000 bp from across the genome of Mediodactylus gecko, enabled us to resolve previously obscure phylogenetic relationships among the five, recently elevated, Mediodactylus species and to support the hypothesis that the taxon includes several new, undescribed species. Population genomic analyses within each of the proposed species showed strong genetic structure and high levels of genetic differentiation among populations.

Dragons in the tropics - Phylogeography and speciation in Diporiphora lizards and common geographic breaks in co-distributed taxa.

Co-distributed taxa can respond both similarly or differently to the same climatic and geological events, resulting in a range of phylogeographic patterns across the region. Using a nested approach on a taxonomically diverse yet morphologically conservative group of agamid lizards, we first aimed to evaluate more precisely the extent of phylogeographic structuring within the genus. Then, focusing on four lineages within the more widespread species, we assessed the impact of biogeographic barriers on phylogeographic structuring and demographic history of species, comparing to patterns previously observed in co-distributed taxa. These species occur in the Australian Monsoonal Tropics, a vast tropical savanna system with high richness and endemism associated with environmental heterogeneity and past climate fluctuations. The employment of genomic data helped to determine the relationships between specific taxa that were previously difficult to place. We found a local influence of biogeographic and climatic breaks on population dynamics, analogous to other species. We detected high levels of population structure in the West Kimberley and Arnhem Plateau, which are already known for high endemism. However, we also highlighted unique lineages in areas that have been overlooked until recently, in the South Kimberley and West Top End. Climatic and geographical features in the Arnhem Plateau act as a soft barrier between populations in the east and west regions of the Top End. These observations reflect patterns observed for other vertebrates across this rich biome, indicating how climatic variation, species' ecology, and landscape features interact to shape regional diversity and endemism.

Spatio-temporal expression patterns of glycine-rich beta proteins and cysteine-rich beta proteins in setae development of Gekko japonicus.

BACKGROUND: Setae on the pad lamellae of the Japanese gecko Gekko japonicus (Schlegel, 1836), a vital epidermal derivative, are primarily composed of cornified beta-proteins (CBPs) and play a pivotal role in adhesion and climbing. The amino acid composition of CBPs might be a determining factor influencing their functional properties. However, the molecular mechanisms governed by CBP genes with diverse amino acid compositions in setae development remain unexplored. RESULTS: Based on RNA-seq analyses, this study confirmed that all G. japonicus CBPs (GjCBPs) are involved in setae formation. Cysteine-rich CBPs encoding genes (ge-cprp-17 to ge-cprp-26) and glycine-rich CBPs encoding genes (ge-gprp-17 to ge-gprp-22) were haphazardly selected, with quantitative real-time PCR revealing their expression patterns in embryonic pad lamellae and dorsal epidermis. It is inferred that glycine-rich CBPs are integral to the formation of both dorsal scales and lamellar setae, cysteine-rich CBPs are primarily associated with setae development. Additionally, fluorescence in situ hybridization revealed spatiotemporal differences in the expression of a glycine-rich CBP encoding gene (ge-gprp-19) and a cysteine-rich CBP encoding gene (ge-cprp-17) during dorsal scales and/or lamellar development. CONCLUSIONS: All 66 CBPs are involved in the formation of setae. Glycine-rich CBPs hold a significant role in the development of dorsal scales and lamellar setae, whereas most cysteine-rich CBPs appear to be essential components of G. japonicus setae. Even GjCBPs with similar amino acid compositions may play diverse functions. The clear spatio-temporal expression differences between the glycine-rich and cysteine-rich CBP encoding genes during epidermal scale and/or setae formation were observed. Embryonic developmental stages 39 to 42 emerged as crucial phases for setae development. These findings lay the groundwork for deeper investigation into the function of GjCBPs in the development of G. japonicus setae.

Plasticity of Gene Expression and Thermal Tolerance: Implications for Climate Change Vulnerability in a Tropical Forest Lizard.

AbstractTropical ectotherms are thought to be especially vulnerable to climate change because they have evolved in temporally stable thermal environments and therefore have decreased tolerance for thermal variability. Thus, they are expected to have narrow thermal tolerance ranges, live close to their upper thermal tolerance limits, and have decreased thermal acclimation capacity. Although models often predict that tropical forest ectotherms are especially vulnerable to rapid environmental shifts, these models rarely include the potential for plasticity of relevant traits. We measured phenotypic plasticity of thermal tolerance and thermal preference as well as multitissue transcriptome plasticity in response to warmer temperatures in a species that previous work has suggested is highly vulnerable to climate warming, the Panamanian slender anole lizard (Anolis apletophallus). We found that many genes, including heat shock proteins, were differentially expressed across tissues in response to short-term warming. Under long-term warming, the voluntary thermal maxima of lizards also increased, although thermal preference exhibited only limited plasticity. Using these data, we modeled changes in the activity time of slender anoles through the end of the century under climate change and found that plasticity should delay declines in activity time by at least two decades. Our results suggest that slender anoles, and possibly other tropical ectotherms, can alter the expression of genes and phenotypes when responding to shifting environmental temperatures and that plasticity should be considered when predicting the future of organisms under a changing climate.

The efficacy of single mitochondrial genes at reconciling the complete mitogenome phylogeny-a case study on dwarf chameleons.

Although genome-scale data generation is becoming more tractable for phylogenetics, there are large quantities of single gene fragment data in public repositories and such data are still being generated. We therefore investigated whether single mitochondrial genes are suitable proxies for phylogenetic reconstruction as compared to the application of full mitogenomes. With near complete taxon sampling for the southern African dwarf chameleons (Bradypodion), we estimated and compared phylogenies for the complete mitogenome with topologies generated from individual mitochondrial genes and various combinations of these genes. Our results show that the topologies produced by single genes (ND2, ND4, ND5, COI, and COIII) were analogous to the complete mitogenome, suggesting that these genes may be reliable markers for generating mitochondrial phylogenies in lieu of generating entire mitogenomes. In contrast, the short fragment of 16S commonly used in herpetological systematics, produced a topology quite dissimilar to the complete mitogenome and its concatenation with ND2 weakened the resolution of ND2. We therefore recommend the avoidance of this 16S fragment in future phylogenetic work.

Fossil-informed biogeographic analysis suggests Eurasian regionalization in crown Squamata during the early Jurassic.

Background: Squamata (lizards, snakes, and amphisbaenians) is a Triassic lineage with an extensive and complex biogeographic history, yet no large-scale study has reconstructed the ancestral range of early squamate lineages. The fossil record indicates a broadly Pangaean distribution by the end- Cretaceous, though many lineages (e.g., Paramacellodidae, Mosasauria, Polyglyphanodontia) subsequently went extinct. Thus, the origin and occupancy of extant radiations is unclear and may have been localized within Pangaea to specific plates, with potential regionalization to distinct Laurasian and Gondwanan landmasses during the Mesozoic in some groups. Methods: We used recent tectonic models to code extant and fossil squamate distributions occurring on nine discrete plates for 9,755 species, with Jurassic and Cretaceous fossil constraints from three extinct lineages. We modeled ancestral ranges for crown Squamata from an extant-only molecular phylogeny using a suite of biogeographic models accommodating different evolutionary processes and fossil-based node constraints from known Jurassic and Cretaceous localities. We hypothesized that the best-fit models would not support a full Pangaean distribution (i.e., including all areas) for the origin of crown Squamata, but would instead show regionalization to specific areas within the fragmenting supercontinent, likely in the Northern Hemisphere where most early squamate fossils have been found. Results: Incorporating fossil data reconstructs a localized origin within Pangaea, with early regionalization of extant lineages to Eurasia and Laurasia, while Gondwanan regionalization did not occur until the middle Cretaceous for Alethinophidia, Scolecophidia, and some crown Gekkotan lineages. While the Mesozoic history of extant squamate biogeography can be summarized as a Eurasian origin with dispersal out of Laurasia into Gondwana, their Cenozoic history is complex with multiple events (including secondary and tertiary recolonizations) in several directions. As noted by previous authors, squamates have likely utilized over-land range expansion, land-bridge colonization, and trans-oceanic dispersal. Tropical Gondwana and Eurasia hold more ancient lineages than the Holarctic (Rhineuridae being a major exception), and some asymmetries in colonization (e.g., to North America from Eurasia during the Cenozoic through Beringia) deserve additional study. Future studies that incorporate fossil branches, rather than as node constraints, into the reconstruction can be used to explore this history further.

High-Quality, Chromosome-Level Reference Genomes of the Viviparous Caribbean Skinks Spondylurus nitidus and S. culebrae.

New World mabuyine skinks are a diverse radiation of morphologically cryptic lizards with unique reproductive biologies. Recent studies examining population-level data (morphological, ecological, and genomic) have uncovered novel biodiversity and phenotypes, including the description of dozens of new species and insights into the evolution of their highly complex placental structures. Beyond the potential for this diverse group to serve as a model for the evolution of viviparity in lizards, much of the taxonomic diversity is concentrated in regions experiencing increasing environmental instability from climate and anthropogenic change. Consequently, a better understanding of genome structure and diversity will be an important tool in the adaptive management and conservation of this group. Skinks endemic to Caribbean islands are particularly vulnerable to global change with several species already considered likely extinct and several remaining species either endangered or threatened. Combining PacBio long-read sequencing, Hi-C, and RNAseq data, here we present the first genomic resources for this group by describing new chromosome-level reference genomes for the Puerto Rican Skink Spondylurus nitidus and the Culebra Skink S. culebrae. Results indicate two high quality genomes, both ∼1.4 Gb, assembled nearly telomere to telomere with complete mitochondrion assembly and annotation.

Morphometric dataset of Varanus salvator for non-invasive sex identification using machine learning.

Reliable sex identification in Varanus salvator traditionally relied on invasive methods like genetic analysis or dissection, as less invasive techniques such as hemipenes inversion are unreliable. Given the ecological importance of this species and skewed sex ratios in disturbed habitats, a dataset that allows ecologists or zoologists to study the sex determination of the lizard is crucial. We present a new dataset containing morphometric measurements of V. salvator individuals from the skin trade, with sex confirmed by dissection post- measurement. The dataset consists of a mixture of primary and secondary data such as weight, skull size, tail length, condition etc. and can be used in modelling studies for ecological and conservation research to monitor the sex ratio of this species. Validity was demonstrated by training and testing six machine learning models. This dataset has the potential to streamline sex determination, offering a non-invasive alternative to complement existing methods in V. salvator research, mitigating the need for invasive procedures.

Adaptive introgression reveals the genetic basis of a sexually selected syndrome in wall lizards.

The joint expression of particular colors, morphologies, and behaviors is a common feature of adaptation, but the genetic basis for such "phenotypic syndromes" remains poorly understood. Here, we identified a complex genetic architecture associated with a sexually selected syndrome in common wall lizards, by capitalizing on the adaptive introgression of coloration and morphology into a distantly related lineage. Consistent with the hypothesis that the evolution of phenotypic syndromes in vertebrates is facilitated by developmental linkage through neural crest cells, most of the genes associated with the syndrome are involved in neural crest cell regulation. A major locus was a ~400-kb region, characterized by standing structural genetic variation and previously implied in the evolutionary innovation of coloration and beak size in birds. We conclude that features of the developmental and genetic architecture contribute to maintaining trait integration, facilitating the extensive and rapid introgressive spread of suites of sexually selected characters.

Cytogenetic Analysis of Satellitome of Madagascar Leaf-Tailed Geckos.

Satellite DNA (satDNA) consists of sequences of DNA that form tandem repetitions across the genome, and it is notorious for its diversity and fast evolutionary rate. Despite its importance, satDNA has been only sporadically studied in reptile lineages. Here, we sequenced genomic DNA and PCR-amplified microdissected W chromosomes on the Illumina platform in order to characterize the monomers of satDNA from the Henkel's leaf-tailed gecko U. henkeli and to compare their topology by in situ hybridization in the karyotypes of the closely related Günther's flat-tail gecko U. guentheri and gold dust day gecko P. laticauda. We identified seventeen different satDNAs; twelve of them seem to accumulate in centromeres, telomeres and/or the W chromosome. Notably, centromeric and telomeric regions seem to share similar types of satDNAs, and we found two that seem to accumulate at both edges of all chromosomes in all three species. We speculate that the long-term stability of all-acrocentric karyotypes in geckos might be explained from the presence of specific satDNAs at the centromeric regions that are strong meiotic drivers, a hypothesis that should be further tested.

Genetic structure and landscape effects on gene flow in the Neotropical lizard Norops brasiliensis (Squamata: Dactyloidae).

One key research goal of evolutionary biology is to understand the origin and maintenance of genetic variation. In the Cerrado, the South American savanna located primarily in the Central Brazilian Plateau, many hypotheses have been proposed to explain how landscape features (e.g., geographic distance, river barriers, topographic compartmentalization, and historical climatic fluctuations) have promoted genetic structure by mediating gene flow. Here, we asked whether these landscape features have influenced the genetic structure and differentiation in the lizard species Norops brasiliensis (Squamata: Dactyloidae). To achieve our goal, we used a genetic clustering analysis and estimate an effective migration surface to assess genetic structure in the focal species. Optimized isolation-by-resistance models and a simulation-based approach combined with machine learning (convolutional neural network; CNN) were then used to infer current and historical effects on population genetic structure through 12 unique landscape models. We recovered five geographically distributed populations that are separated by regions of lower-than-expected gene flow. The results of the CNN showed that geographic distance is the sole predictor of genetic variation in N. brasiliensis, and that slope, rivers, and historical climate had no discernible influence on gene flow. Our novel CNN approach was accurate (89.5%) in differentiating each landscape model. CNN and other machine learning approaches are still largely unexplored in landscape genetics studies, representing promising avenues for future research with increasingly accessible genomic datasets.

Species divergence in valleys: the phylogeny of Phrynocephalus forsythii complex and description of a new species.

Background: Geographic isolation caused by high-altitude valleys promotes the formation of geographic segregation of species, leading to species differentiation. The subgenus Oreosaura contains viviparous species from the Tibetan Plateau and the vicinity of the Tarim Basin, which can be divided into three species complexes according to their geographical distribution: Phrynocephalus vlangalii, Phrynocephalus theobaldi, and Phrynocephalus forsythii. However, molecular data for the P. forsythii complex are limited and the diversity of this species complex has been greatly underestimated. Therefore, this study aimed to explore the species diversity of Oreosaura and species differentiation within the P. forsythii complex. Methods: We analysed the species diversity of Oreosaura by combining previous data, constructed a phylogenetic tree of the subgenus based on cytochrome c oxidase subunit I and 16S sequences, and estimated the divergence time. Results: The results suggest significant genetic differences between the Tarim Basin populations and adjacent mountain valley populations of the P. forsythii complex and that the combination of deep valley landscapes in the high mountains and ice-age events have contributed to the differentiation of the viviparous toad-headed agama lizard, which is a key factor in the phylogenetics of the P. forsythii complex. Furthermore, we identified a population collected from Wuqia County, Xinjiang, as a new species, Phrynocephalus kangsuensis sp. nov. The results will provide data for phylogenetic studies following the P. forsythii complex and help demonstrate that valleys promote the formation of Phrynocephalus species.

The role of historical biogeography in shaping colour morph diversity in the common wall lizard.

The maintenance of polymorphisms often depends on multiple selective forces, but less is known on the role of stochastic or historical processes in maintaining variation. The common wall lizard (Podarcis muralis) is a colour polymorphic species in which local colour morph frequencies are thought to be modulated by natural and sexual selection. Here, we used genome-wide single-nucleotide polymorphism data to investigate the relationships between morph composition and population biogeography at a regional scale, by comparing morph composition with patterns of genetic variation of 54 populations sampled across the Pyrenees. We found that genetic divergence was explained by geographic distance but not by environmental features. Differences in morph composition were associated with genetic and environmental differentiation, as well as differences in sex ratio. Thus, variation in colour morph frequencies could have arisen via historical events and/or differences in the permeability to gene flow, possibly shaped by the complex topography and environment. In agreement with this hypothesis, colour morph diversity was positively correlated with genetic diversity and rates of gene flow and inversely correlated with the likelihood of the occurrence of bottlenecks. Concurrently, we did not find conclusive evidence for selection in the two colour loci. As an illustration of these effects, we observed that populations with higher proportions of the rarer yellow and yellow-orange morphs had higher genetic diversity. Our results suggest that processes involving a decay in overall genetic diversity, such as reduced gene flow and/or bottleneck events have an important role in shaping population-specific morph composition via non-selective processes.

Transcriptomic profiling revealed immune-related signaling pathways in response to experimental infection of Leishmania donovani in two desert lizards from Northwest China.

Little is known about the immune response of lizards to Leishmania parasties. In this study, we conducted the first liver transcriptome analysis of two lizards (Phrynocephalus przewalskii and Eremias multiocellata) challenged with L. donovani, endemic to the steppe desert region of northwestern China. Our results revealed that multiple biological processes and immune-related signaling pathways are closely associated with the immune response to experimental L. donovani infection in the two lizards, and that both lizards show similar changes to mammals in terms of immunity to Leishmania. However, the interspecific divergence of the two lizards leads to different transcriptomic changes. In particular, in contrast to P. przewalskii, the challenged E. mutltiocellata was characterized by the induction of down-regulation of most DEGs. These findings will contribute to the scarce resources on lizard immunity and provide a reference for further research on immune mechanisms in reptiles.

Digest: Liolaemus lizards as an emerging model system to study hybridization as a driver of rapid radiations.

Recent perspectives on speciation genomics emphasize the pivotal role of hybridization in driving rapid radiations. The Liolaemus lizard genus displays impressive species richness with around 290 species widely distributed across southern South America. Sánchez et al. conducted a comprehensive study on the 5-million-year-old Liolaemus kingii group, which includes 14 species. The research provides new key insights to enhance our understanding of this rapid radiation, including its diversification in space and time and the consequences of hybridization in its morphological evolution and taxonomy.

Population genomics of flat-tailed horned lizards (Phrynosoma mcallii) informs conservation and management across a fragmented Colorado Desert landscape.

Phrynosoma mcallii (flat-tailed horned lizards) is a species of conservation concern in the Colorado Desert of the United States and Mexico. We analysed ddRADseq data from 45 lizards to estimate population structure, infer phylogeny, identify migration barriers, map genetic diversity hotspots, and model demography. We identified the Colorado River as the main geographic feature contributing to population structure, with the populations west of this barrier further subdivided by the Salton Sea. Phylogenetic analysis confirms that northwestern populations are nested within southeastern populations. The best-fit demographic model indicates Pleistocene divergence across the Colorado River, with significant bidirectional gene flow, and a severe Holocene population bottleneck. These patterns suggest that management strategies should focus on maintaining genetic diversity on both sides of the Colorado River and the Salton Sea. We recommend additional lands in the United States and Mexico that should be considered for similar conservation goals as those in the Rangewide Management Strategy. We also recommend periodic rangewide genomic sampling to monitor ongoing attrition of diversity, hybridization, and changing structure due to habitat fragmentation, climate change, and other long-term impacts.

Genetic Basis and Evolutionary Forces of Sexually Dimorphic Color Variation in a Toad-Headed Agamid Lizard.

In the animal kingdom, sexually dimorphic color variation is a widespread phenomenon that significantly influences survival and reproductive success. However, the genetic underpinnings of this variation remain inadequately understood. Our investigation into sexually dimorphic color variation in the desert-dwelling Guinan population of the toad-headed agamid lizard (Phrynocephalus putjatai) utilized a multidisciplinary approach, encompassing phenotypic, ultrastructural, biochemical, genomic analyses, and behavioral experiments. Our findings unveil the association between distinct skin colorations and varying levels of carotenoid and pteridine pigments. The red coloration in males is determined by a genomic region on chromosome 14, housing four pigmentation genes: BCO2 and three 6-pyruvoyltetrahydropterin synthases. A Guinan population-specific nonsynonymous single nucleotide polymorphism in BCO2 is predicted to alter the electrostatic potential within the binding domain of the BCO2-ß-carotene complex, influencing their interaction. Additionally, the gene MAP7 on chromosome 2 emerges as a potential contributor to the blue coloration in subadults and adult females. Sex-specific expression patterns point to steroid hormone-associated genes (SULT2B1 and SRD5A2) as potential upstream regulators influencing sexually dimorphic coloration. Visual modeling and field experiments support the potential selective advantages of vibrant coloration in desert environments. This implies that natural selection, potentially coupled with assortative mating, might have played a role in fixing color alleles, contributing to prevalence in the local desert habitat. This study provides novel insights into the genetic basis of carotenoid and pteridine-based color variation, shedding light on the evolution of sexually dimorphic coloration in animals. Moreover, it advances our understanding of the driving forces behind such intricate coloration patterns.