Genomics of coloration in natural animal populations

Genetic basis and expression of ventral colour in polymorphic common lizards

Colour is an important visual cue that can correlate with sex, behaviour, life history or ecological strategies, and has evolved divergently and convergently across animal lineages. Its genetic basis in non-model organisms is rarely known, but such information is vital for determining the drivers and mechanisms of colour evolution. Leveraging genetic admixture in a rare contact zone between oviparous and viviparous common lizards (Zootoca vivipara), we show that females (N = 558) of the two otherwise morphologically indistinguishable reproductive modes differ in their ventral colouration (from pale to vibrant yellow) and intensity of melanic patterning. We find no association between female colouration and reproductive investment, and no evidence for selection on colour. Using a combination of genetic mapping and transcriptomic evidence, we identified two candidate genes associated with ventral colour differentiation, DGAT2 and PMEL. These are genes known to be involved in carotenoid metabolism and melanin synthesis respectively. Ventral melanic spots were associated with two genomic regions, including a SNP close to protein tyrosine phosphatase (PTP) genes. Using genome re-sequencing data, our results show that fixed coding mutations in the candidate genes cannot account for differences in colouration. Taken together, our findings show that the evolution of ventral colouration and its associations across common lizard lineages is variable. A potential genetic mechanism explaining the flexibility of ventral colouration may be that colouration in common lizards, but also across squamates, is predominantly driven by regulatory genetic variation.

The genomic architecture of continuous plumage colour variation in the European barn owl (Tyto alba)

The maintenance of colour variation in wild populations has long fascinated evolutionary biologists, although most studies have focused on discrete traits exhibiting rather simple inheritance patterns and genetic architectures. However, the study of continuous colour traits and their potentially oligo- or polygenic genetic bases remains rare in wild populations. We studied the genetics of the continuously varying white-to-rufous plumage coloration of the European barn owl (Tyto alba) using a genome-wide association approach on the whole-genome data of 75 individuals. We confirmed a mutation at the melanocortin-1-receptor gene (MC1R) is involved in the coloration and identified two new regions, located in super-scaffolds 9 and 42. The combination of the three regions explains most of the colour variation (80.37%, 95% credible interval 58.45-100%). One discovered region, located in the sex chromosome, differs between the most extreme colorations in owls sharing a specific MC1R genotype. This region may play a role in the colour sex dimorphism of this species, possibly in interaction with the autosomal MC1R. We thus provide insights into the genetic architecture of continuous colour variation, pointing to an oligogenic basis with potential epistatic effects among loci that should aid future studies understanding how continuous colour variation is maintained in nature.

Genetic variants in melanogenesis proteins TYRP1 and TYR are associated with the golden rhesus macaque phenotype

Nonhuman primates (NHPs) are vital translational research models due to their high genetic, physiological, and anatomical homology with humans. The "golden" rhesus macaque (Macaca mulatta) phenotype is a naturally occurring, inherited trait with a visually distinct pigmentation pattern resulting in light blonde colored fur. Retinal imaging also reveals consistent hypopigmentation and occasional foveal hypoplasia. Here, we describe the use of genome-wide association in 2 distinct NHP populations to identify candidate variants in genes linked to the golden phenotype. Two missense variants were identified in the Tyrosinase-related protein 1 gene (Asp343Gly and Leu415Pro) that segregate with the phenotype. An additional and distinct association was also found with a Tyrosinase variant (His256Gln), indicating the light-colored fur phenotype can result from multiple genetic mechanisms. The implicated genes are related through their contribution to the melanogenesis pathway. Variants in these 2 genes are known to cause pigmentation phenotypes in other species and to be associated with oculocutaneous albinism in humans. The novel associations presented in this study will permit further investigations into the role these proteins and variants play in the melanogenesis pathway and model the effects of genetic hypopigmentation and altered melanogenesis in a naturally occurring nonhuman primate model.

Genetics of yellow-orange color variation in a pair of sympatric sulphur butterflies

Continuous color polymorphisms can serve as a tractable model for the genetic and developmental architecture of traits. Here we investigated continuous color variation in Colias eurytheme and Colias philodice, two species of sulphur butterflies that hybridize in sympatry. Using quantitative trait locus (QTL) analysis and high-throughput color quantification, we found two interacting large-effect loci affecting orange-to-yellow chromaticity. Knockouts of red Malpighian tubules (red), likely involved in endosomal maturation, result in depigmented wing scales. Additionally, the transcription factor bric-a-brac can act as a modulator of orange pigmentation. We also describe the QTL architecture of other continuously varying traits, together supporting a large-X effect model where the genetic control of species-defining traits is enriched on sex chromosomes. This study sheds light on the range of possible genetic architectures that can underpin a continuously varying trait and illustrates the power of using automated measurement to score phenotypes that are not always conspicuous to the human eye.

Redundancy analysis, genome-wide association studies and the pigmentation of brown trout (Salmo trutta L.)

The association of molecular variants with phenotypic variation is a main issue in biology, often tackled with genome-wide association studies (GWAS). GWAS are challenging, with increasing, but still limited, use in evolutionary biology. We used redundancy analysis (RDA) as a complimentary ordination approach to single- and multitrait GWAS to explore the molecular basis of pigmentation variation in brown trout (Salmo trutta) belonging to wild populations impacted by hatchery fish. Based on 75,684 single nucleotide polymorphic (SNP) markers, RDA, single- and multitrait GWAS allowed the extraction of 337 independent colour patterning loci (CPLs) associated with trout pigmentation traits, such as the number of red and black spots on flanks. Collectively, these CPLs (i) mapped onto 35 out of 40 brown trout linkage groups indicating a polygenic genomic architecture of pigmentation, (ii) were found to be associated with 218 candidate genes, including 197 genes formerly mentioned in the literature associated to skin pigmentation, skin patterning, differentiation or structure notably in a close relative, the rainbow trout (Onchorhynchus mykiss), and (iii) related to functions relevant to pigmentation variation (e.g., calcium- and ion-binding, cell adhesion). Annotated CPLs include genes with well-known pigmentation effects (e.g., PMEL, SLC45A2, SOX10), but also markers associated with genes formerly found expressed in rainbow or brown trout skins. RDA was also shown to be useful to investigate management issues, especially the dynamics of trout pigmentation submitted to several generations of hatchery introgression.

Genetic and environmental drivers of colour and pattern in the Australian jacky dragon (Amphibolurus muricatus)

The underlying drivers of variation in the colouration (colour and pattern) of animals can be genetic, non‐genetic, or more likely, a combination of both. Understanding the role of heritable genetic elements, as well as non‐genetic factors such as age, habitat or temperature, in shaping colouration can provide insight into the evolution and function of these traits, as well as the speed of response to changing environments. This project examined the genetic and non‐genetic drivers of continuous variation in colouration in a lizard, the jacky dragon (Amphibolurus muricatus). We leveraged a large captive experiment that manipulated parental and offspring thermal environment to simultaneously estimate the genetic and non‐genetic drivers of variation in colouration. We found that the overall brightness, the elongation of the longitudinal stripes on the dorsum and the contrast between light and dark patches of the pattern were all heritable. Colouration varied according to the age of the hatchling; however, the thermal environment of neither the parents nor offspring contributed significantly to colouration. It appears that developmental plasticity and maternal effects associated with temperature are not important drivers of variation in our measures of colouration.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

The genetic basis of a novel reproductive strategy in Sulawesi ricefishes: How modularity and a low number of loci shape pelvic brooding

The evolution of complex phenotypes like reproductive strategies is challenging to understand, as they often depend on multiple adaptations that only jointly result in a specific functionality. Sulawesi ricefishes (Adrianichthyidae) evolved a reproductive strategy termed as pelvic brooding. In contrast to the more common transfer brooding, female pelvic brooders carry an egg bundle connected to their body for weeks until the fry hatches. To examine the genetic architecture of pelvic brooding, we crossed the pelvic brooding Oryzias eversi and the transfer brooding Oryzias nigrimas (species divergence time: ∼3.6 my). We hypothesize, that a low number of loci and modularity have facilitated the rapid evolution of pelvic brooding. Traits associated to pelvic brooding, like rib length, pelvic fin length, and morphology of the genital papilla, were correlated in the parental species but correlations were reduced or lost in their F1 and F2 hybrids. Using the Castle-Wright estimator, we found that generally few loci underlie the studied traits. Further, both parental species showed modularity in their body plans. In conclusion, morphological traits related to pelvic brooding were based on a few loci and the mid-body region likely could evolve independently from the remaining body parts. Both factors presumably facilitated the evolution of pelvic brooding.

CamoEvo: An open access toolbox for artificial camouflage evolution experiments

Camouflage research has long shaped our understanding of evolution by natural selection, and elucidating the mechanisms by which camouflage operates remains a key question in visual ecology. However, the vast diversity of color patterns found in animals and their backgrounds, combined with the scope for complex interactions with receiver vision, presents a fundamental challenge for investigating optimal camouflage strategies. Genetic algorithms (GAs) have provided a potential method for accounting for these interactions, but with limited accessibility. Here, we present CamoEvo, an open-access toolbox for investigating camouflage pattern optimization by using tailored GAs, animal and egg maculation theory, and artificial predation experiments. This system allows for camouflage evolution within the span of just 10-30 generations (∼1-2 min per generation), producing patterns that are both significantly harder to detect and that are optimized to their background. CamoEvo was built in ImageJ to allow for integration with an array of existing open access camouflage analysis tools. We provide guides for editing and adjusting the predation experiment and GA as well as an example experiment. The speed and flexibility of this toolbox makes it adaptable for a wide range of computer-based phenotype optimization experiments.

Fur colour in the Arctic fox: genetic architecture and consequences for fitness

Genome-wide association studies provide good opportunities for studying the genetic basis of adaptive traits in wild populations. Yet, previous studies often failed to identify major effect genes. In this study, we used high-density single nucleotide polymorphism and individual fitness data from a wild non-model species. Using a whole-genome approach, we identified the MC1R gene as the sole causal gene underlying Arctic fox Vulpes lagopus fur colour. Further, we showed the adaptive importance of fur colour genotypes through measures of fitness that link ecological and evolutionary processes. We found a tendency for blue foxes that are heterozygous at the fur colour locus to have higher fitness than homozygous white foxes. The effect of genotype on fitness was independent of winter duration but varied with prey availability, with the strongest effect in years of increasing rodent populations. MC1R is located in a genomic region with high gene density, and we discuss the potential for indirect selection through linkage and pleiotropy. Our study shows that whole-genome analyses can be successfully applied to wild species and identify major effect genes underlying adaptive traits. Furthermore, we show how this approach can be used to identify knowledge gaps in our understanding of interactions between ecology and evolution.

Avian Coloration Genetics: Recent Advances and Emerging Questions

The colorful phenotypes of birds have long provided rich source material for evolutionary biologists. Avian plumage, beaks, skin, and eggs-which exhibit a stunning range of cryptic and conspicuous forms-inspired early work on adaptive coloration. More recently, avian color has fueled discoveries on the physiological, developmental, and-increasingly-genetic mechanisms responsible for phenotypic variation. The relative ease with which avian color traits can be quantified has made birds an attractive system for uncovering links between phenotype and genotype. Accordingly, the field of avian coloration genetics is burgeoning. In this review, we highlight recent advances and emerging questions associated with the genetic underpinnings of bird color. We start by describing breakthroughs related to 2 pigment classes: carotenoids that produce red, yellow, and orange in most birds and psittacofulvins that produce similar colors in parrots. We then discuss structural colors, which are produced by the interaction of light with nanoscale materials and greatly extend the plumage palette. Structural color genetics remain understudied-but this paradigm is changing. We next explore how colors that arise from interactions among pigmentary and structural mechanisms may be controlled by genes that are co-expressed or co-regulated. We also identify opportunities to investigate genes mediating within-feather micropatterning and the coloration of bare parts and eggs. We conclude by spotlighting 2 research areas-mechanistic links between color vision and color production, and speciation-that have been invigorated by genetic insights, a trend likely to continue as new genomic approaches are applied to non-model species.

Ecological basis and genetic architecture of crypsis polymorphism in the desert clicker grasshopper (Ligurotettix coquilletti)

Color polymorphic species can offer exceptional insight into the ecology and genetics of adaptation. Although the genetic architecture of animal coloration is diverse, many color polymorphisms are associated with large structural variants and maintained by biotic interactions. Grasshoppers are notably polymorphic in both color and karyotype, which makes them excellent models for understanding the ecological drivers and genetic underpinnings of color variation. Banded and uniform morphs of the desert clicker grasshopper (Ligurotettix coquilletti) are found across the western deserts of North America. To address the hypothesis that predation maintains local color polymorphism and shapes regional crypsis variation, we surveyed morph frequencies and tested for covariation with two predation environments. Morphs coexisted at intermediate frequencies at most sites, consistent with local balancing selection. Morph frequencies covaried with the appearance of desert substrate-an environment used only by females-suggesting that ground-foraging predators are major agents of selection on crypsis. We next addressed the hypothesized link between morph variation and genome structure. To do so, we designed an approach for detecting inversions and indels using only RADseq data. The banded morph was perfectly correlated with a large putative indel. Remarkably, indel dominance differed among populations, a rare example of dominance evolution in nature.

Variation in Genetic Mechanisms for Plumage Polymorphism in Skuas (Stercorarius)

Coloration is evolutionarily labile and so provides an excellent trait for examining the repeatability of evolution. Here, we investigate the repeatability of the evolution of polymorphic variation in ventral plumage coloration in skuas (Stercorarius: Stercorariidae). In 2 species, arctic (S. parasiticus) and pomarine skuas (S. pomarinus), plumage polymorphism was previously shown to be associated with coding changes at the melanocortin-1 receptor (MC1R) locus. Here, we show that polymorphism in a third species, the south polar skua (S. maccormicki), is not associated with coding variation at MC1R or with variation at a Z-linked second candidate locus, tyrosinase-related protein 1 (TYRP1). Hence, convergent evolution of plumage polymorphisms in skuas is only partly repeatable at the level of the genetic locus involved. Interestingly, the pattern of repeatability in skuas is aligned not with phylogeny but with the nature of the phenotypic variation. In particular, south polar skuas show a strong sex bias to coloration that is absent in the other species, and it may be that this has a unique genetic architecture.

Epistatic interactions between pterin and carotenoid genes modulate intra-morph color variation in a lizard

Color polymorphisms have become a major topic in evolutionary biology and substantial efforts have been devoted to the understanding of the mechanisms responsible for originating such colorful systems. Within-morph continuous variation, on the other hand, has been neglected in most of the studies. Here, we combine spectrophotometric/visual modeling and genetic data to study the mechanisms promoting continuous variation within categorical color morphs of Podarcis muralis. Our results suggest that intra-morph variability in the pterin-based orange morph is greater compared to white and yellow morphs. We also show that continuous variation within the orange morph is partially discriminable by conspecifics. Genotyping results indicate that allelic variants at the BCO2 locus (responsible for deposition of yellow carotenoids) contribute to generate continuous variation in orange individuals. However, other intrinsic and/or extrinsic mechanisms, such as body size, might be involved, opening a new avenue for future research on the drivers of continuous variation within-morphs.

A combined RAD-Seq and WGS approach reveals the genomic basis of yellow color variation in bumble bee Bombus terrestris

Bumble bees exhibit exceptional diversity in their segmental body coloration largely as a result of mimicry. In this study we sought to discover genes involved in this variation through studying a lab-generated mutant in bumble bee Bombus terrestris, in which the typical black coloration of the pleuron, scutellum, and first metasomal tergite is replaced by yellow, a color variant also found in sister lineages to B. terrestris. Utilizing a combination of RAD-Seq and whole-genome re-sequencing, we localized the color-generating variant to a single SNP in the protein-coding sequence of transcription factor cut. This mutation generates an amino acid change that modifies the conformation of a coiled-coil structure outside DNA-binding domains. We found that all sequenced Hymenoptera, including sister lineages, possess the non-mutant allele, indicating different mechanisms are involved in the same color transition in nature. Cut is important for multiple facets of development, yet this mutation generated no noticeable external phenotypic effects outside of setal characteristics. Reproductive capacity was reduced, however, as queens were less likely to mate and produce female offspring, exhibiting behavior similar to that of workers. Our research implicates a novel developmental player in pigmentation, and potentially caste, thus contributing to a better understanding of the evolution of diversity in both of these processes.

Fitness and fur colouration: Testing the camouflage and thermoregulation hypotheses in an Arctic mammal

Selection for crypsis has been recognized as an important ecological driver of animal colouration, whereas the relative importance of thermoregulation is more contentious with mixed empirical support. A potential thermal advantage of darker individuals has been observed in a wide range of animal species. Arctic animals that exhibit colour polymorphisms and undergo seasonal colour moults are interesting study subjects for testing the two alternative hypotheses: demographic performance of different colour morphs might be differentially affected by snow cover with a cryptic advantage for lighter morphs, or conversely by winter temperature with a thermal advantage for darker morphs. In this study, we explored whether camouflage and thermoregulation might explain differences in reproduction and survival between the white and blue colour morphs of the Arctic fox Vulpes lagopus under natural conditions. Juvenile and adult survival, breeding propensity and litter size were measured for 798 captive-bred and released or wild-born Arctic foxes monitored during an 11-year period (2007-2017) in two subpopulations in south-central Norway. We investigated the proportion of the two colour morphs and compared their demographic performance in relation to spatial variation in duration of snow cover, onset of snow season and winter temperatures. After population re-establishment, a higher proportion of blue individuals was observed among wild-born Arctic foxes compared to the proportion of blue foxes released from the captive population. Our field study provides the first evidence for an effect of colour morph on the reproductive performance of Arctic foxes under natural conditions, with a higher breeding propensity of the blue morph compared to the white one. Performance of the two colour morphs was not differentially affected by the climatic variables, except for juvenile survival. Blue morph juveniles showed a tendency for higher survival under colder winter temperatures but lower survival under warmer temperatures compared to white morph juveniles. Overall, our findings do not consistently support predictions of the camouflage or the thermoregulation hypotheses. The higher success of blue foxes suggests an advantage of the dark morph not directly related to disruptive selection by crypsis or thermoregulation. Our results rather point to physiological adaptations and behavioural traits not necessarily connected to thermoregulation, such as stress response, immune function, sexual behaviour and aggressiveness. Our findings highlight the need to explore the potential role of genetic linkage or pleiotropy in influencing the fitness of white and blue Arctic foxes as well as other species with colour polymorphisms.

Overview of the betta fish genome regarding species radiation, parental care, behavioral aggression, and pigmentation model relevant to humans

BACKGROUND

The Siamese fighting fish (Betta splendens, also known as the betta) is well known in aquarium markets, and also presents an exciting new research model for studying parental care, aggressive behavior, and cryptically diverse pigmentation. However, concentrated efforts are required, both in the context of conservation biology and in its genetics, to address the problems of ongoing outbreeding depression, loss of biodiversity, and lack of scientific biological information.

OBJECTIVE

The evolutionary dynamics of the betta must be better understood at the genomic scale in order to resolve the phylogenetic status of unrecognized species, develop molecular markers to study variation in traits, and identify interesting sets of genes encoding various bioresource functions.

METHODS

The recent revolution in multi-omics approaches such as genomics, transcriptomics, epigenomics, and proteomics has uncovered genetic diversity and gained insights into many aspects of betta bioresources.

RESULTS

Here, we present current research and future plans in an ongoing megaproject to characterize the betta genome as de novo assemblies, genes and repeat annotations, generating data to study diverse biological phenomena. We highlight key questions that require answers and propose new directions and recommendations to develop bioresource management to protect and enhance the betta genus.

CONCLUSION

Successful accomplishment of these plans will allow the creation of a reference annotated genome and provide valuable information at the molecular level that can be utilized to sustain biodiversity and eco-management of the betta to improve breeding programs for future biomedical research.

Extensive hybridization reveals multiple coloration genes underlying a complex plumage phenotype

Coloration is an important target of both natural and sexual selection. Discovering the genetic basis of colour differences can help us to understand how this visually striking phenotype evolves. Hybridizing taxa with both clear colour differences and shallow genomic divergences are unusually tractable for associating coloration phenotypes with their causal genotypes. Here, we leverage the extensive admixture between two common North American woodpeckers-yellow-shafted and red-shafted flickers-to identify the genomic bases of six distinct plumage patches involving both melanin and carotenoid pigments. Comparisons between flickers across approximately 7.25 million genome-wide SNPs show that these two forms differ at only a small proportion of the genome (mean F_ST = 0.008). Within the few highly differentiated genomic regions, we identify 368 SNPs significantly associated with four of the six plumage patches. These SNPs are linked to multiple genes known to be involved in melanin and carotenoid pigmentation. For example, a gene (CYP2J19) known to cause yellow to red colour transitions in other birds is strongly associated with the yellow versus red differences in the wing and tail feathers of these flickers. Additionally, our analyses suggest novel links between known melanin genes and carotenoid coloration. Our finding of patch-specific control of plumage coloration adds to the growing body of literature suggesting colour diversity in animals could be created through selection acting on novel combinations of coloration genes.

The Biology of Polymorphic Melanic Side-Spotting Patterns in Poeciliid Fishes

Melanin-based color patterns are an emerging model for studying molecular and evolutionary mechanisms driving phenotypic correlations. Extensive literature exists on color patterns and their correlated traits in the family Poeciliidae, indicating that these fishes are tractable models. We review the biology of polymorphic melanic side-spotting patterns characterized by macromelanophores forming irregular spotted patterns across fishes’ flanks. These patterns are present in the generaGambusia, Limia, Phalloceros, Poecilia, andXiphophorus. Their presence is controlled by dominant genes on autosomes or sex chromosomes. Variation in expression is under polygenic control; however, these genes’ identities are still largely unknown. In someGambusia holbrookiandPoecilia latipinna, expression is dependent on low temperature exposure, but underlying molecular mechanisms are unknown. Spotted fish develop melanoma in rare cases and are a well-developed model for melanoma research. Little is known about other physiological correlates except that spottedG. holbrookimales exhibit higher basal cortisol levels than unspotted males and that metabolic rate does not differ between morphs in someXiphophorusspecies. Behavioral differences between morphs are widespread, but specific to population, species, and social context. SpottedG. holbrookimales appear to be more social and more dominant. Juvenile spottedG. holbrookihave lower behavioral flexibility, and spottedX. variatusexhibit greater stress resistance. Findings conflict on whether morphs differ in sexual behavior and in sexual selection by females. Melanic side-spotting patterns are uncommon (<30%) in populations, although extreme high-frequency populations exist. This low frequency is surprising for dominant genes, indicating that a variety of selective pressures influence both these patterns and their correlated traits. Little is known about reproductive life history traits. SpottedG. holbrookiare larger and have higher survival when uncommon, but underlying mechanisms remain unknown. Spotted morphs appear to have a strong selective advantage during predation. Predators prefer to attack and consume unspotted morphs; however, this preference disappears when spottedG. holbrookimales are common, indicating negative frequency-dependent selection. Spotted morphs are preferred socially under turbid conditions, but other environmental factors that shape phenotypic correlations and morph fitness have not been studied. Finally, we present questions for future studies on melanic side-spotting patterns.

Population genetic structure underlying the geographic variation in beetle structural colour with multiple transition zones

We studied the population genetic structure underlying the geographic variation in the structural colour of the geotrupid dung beetle, Phelotrupes auratus, which exhibits metallic body colours of different reflectance wavelengths perceived as red, green and indigo. These forms occur parapatrically in an area of Japan. The colour variation was not related to variation in climatic factors. Using single nucleotide polymorphisms (SNPs) from restriction-site-associated DNA sequences, we discriminated five groups of populations (west/red, south/green, south/indigo, south/red and east/red) by a combination of genetic clusters (west, south and east) and three colour forms. There were three transition zones for the colour forms: two between the red and green forms were hybrid zones with steep genetic clines, which implies the existence of barriers to gene flow between regions with different colours. The remaining transition zone between the green and indigo forms lacked genetic differentiation, despite the evident colour changes, which implies regionally specific selection on the different colours. In a genomewide association study, we identified four SNPs that were associated with the red/green or indigo colour and were not linked with one another, which implies that the coloration was controlled by multiple loci, each affecting the expression of a different colour range. These loci may have controlled the transitions between different combinations of colours. Our study demonstrates that geographic colour variation within a species can be maintained by nonuniform interactions among barriers to gene flow, locally specific selection on different colours, and the effects of different colour loci.

Molecular Plasticity in Animal Pigmentation: Emerging Processes Underlying Color Changes

Animal coloration has been rigorously studied and has provided morphological implications for fitness with influences over social behavior, predator–prey interactions, and sexual selection. In vertebrates, its study has developed our understanding across diverse fields ranging from behavior to molecular biology. In the search for underlying molecular mechanisms, many have taken advantage of pedigree-based and genome-wide association screens to reveal the genetic architecture responsible for pattern variation that occurs in early development. However, genetic differences do not provide a full picture of the dynamic changes in coloration that are most prevalent across vertebrates at the molecular level. Changes in coloration that occur in adulthood via phenotypic plasticity rely on various social, visual, and dietary cues independent of genetic variation. Here, I will review the contributions of pigment cell biology to animal color changes and recent studies describing their molecular underpinnings and function. In this regard, conserved epigenetic processes such as DNA methylation play a role in lending plasticity to gene regulation as it relates to chromatophore function. Lastly, I will present African cichlids as emerging models for the study of pigmentation and molecular plasticity for animal color changes. I posit that these processes, in a dialog with environmental stimuli, are important regulators of variation and the selective advantages that accompany a change in coloration for vertebrate animals.

The evolution of insect body coloration under changing climates

Insects have been influential models in research on color variation, its evolutionary drivers and the mechanistic basis of such variation. More recently, several studies have indicated that insect color is responding to rapid climate change. However, it remains challenging to ascertain drivers of color variation among populations and species, and across space and time, as multiple biotic and abiotic factors can interact and mediate color change. Here, we describe some of the challenges and recent advances made in this field. First, we outline the main alternative hypotheses that exist for insect color variation in relation to climatic factors. Second, we review the existing evidence for contemporary adaptive evolution of insect color in response to climate change and then discuss factors that can promote or hinder the evolution of color in response to climate change. Finally, we propose future directions and highlight gaps in this research field. Pigments and structures producing insect color can vary concurrently or independently, and may evolve at different rates, with poorly understood effects on gene frequencies and fitness. Disentangling multiple competing hypotheses explaining insect coloration should be key to assign color variation as an evolutionary response to climate change.

Being red, blue and green: the genetic basis of coloration differences in the strawberry poison frog (Oophaga pumilio)

Background

Animal coloration is usually an adaptive attribute, under strong local selection pressures and often diversified among species or populations. The strawberry poison frog (Oophaga pumilio) shows an impressive array of color morphs across its distribution in Central America. Here we quantify gene expression and genetic variation to identify candidate genes involved in generating divergence in coloration between populations of red, green and blue O. pumilio from the Bocas del Toro archipelago in Panama.

Results

We generated a high quality non-redundant reference transcriptome by mapping the products of genome-guided and de novo transcriptome assemblies onto a re-scaffolded draft genome of O. pumilio. We then measured gene expression in individuals of the three color phenotypes and identified color-associated candidate genes by comparing differential expression results against a list of a priori gene sets for five different functional categories of coloration – pteridine synthesis, carotenoid synthesis, melanin synthesis, iridophore pathways (structural coloration), and chromatophore development. We found 68 candidate coloration loci with significant expression differences among the color phenotypes. Notable upregulated examples include pteridine synthesis genes spr, xdh and pts (in red and green frogs); carotenoid metabolism genes bco2 (in blue frogs), scarb1 (in red frogs), and guanine metabolism gene psat1 (in blue frogs). We detected significantly higher expression of the pteridine synthesis gene set in red and green frogs versus blue frogs. In addition to gene expression differences, we identified 370 outlier SNPs on 162 annotated genes showing signatures of diversifying selection, including eight pigmentation-associated genes.

Conclusions

Gene expression in the skin of the three populations of frogs with differing coloration is highly divergent. The strong signal of differential expression in pteridine genes is consistent with a major role of these genes in generating the coloration differences among the three morphs. However, the finding of differentially expressed genes across pathways and functional categories suggests that multiple mechanisms are responsible for the coloration differences, likely involving both pigmentary and structural coloration. In addition to regulatory differences, we found potential evidence of differential selection acting at the protein sequence level in several color-associated loci, which could contribute to the color polymorphism.

Ecotypic changes of alpine birds to climate change

In endotherm animals, several traits are related to climate. For example, Bergmann’s rule predicts a decrease in body size within species and across closely related species with increasing temperature, whereas Gloger’s rule states that birds and mammals should be darker in humid and warm environments compared to colder and drier areas. However, it is still not clear whether ecotypic responses to variation in the local environment can also apply to morphological and colouration changes through time in response to climate change. We present a 100-year-long time series on morphological and melanin-based colours of snowfinch (325 Montifringilla, 92 Pyrgilauda and 30 Onychostruthus) museum specimens. Here we show that the tarsus length of the species has decreased and the saturation of the melanin-based colour has increased, which was correlated with the increase of temperature and precipitations. As ecotypic variations are tightly linked to individual behavioural and physiological responses to environmental variations, differently sized and coloured individuals are expected to be differently penalized by global changes. This study opens the pertinent question about whether ecotypic responses can enhance population persistence in the context of global change.

Melanin-based coloration and host-parasite interactions under global change

The role of parasites in shaping melanin-based colour polymorphism, and the consequences of colour polymorphism for disease resistance, remain debated. Here we review recent evidence of the links between melanin-based coloration and the behavioural and immunological defences of vertebrates against their parasites. First we propose that (1) differences between colour morphs can result in variable exposure to parasites, either directly (certain colours might be more or less attractive to parasites) or indirectly (variations in behaviour and encounter probability). Once infected, we propose that (2) immune variation between differently coloured individuals might result in different abilities to cope with parasite infection. We then discuss (3) how these different abilities could translate into variable sexual and natural selection in environments varying in parasite pressure. Finally, we address (4) the potential role of parasites in the maintenance of melanin-based colour polymorphism, especially in the context of global change and multiple stressors in human-altered environments. Because global change will probably affect both coloration and the spread of parasitic diseases in the decades to come, future studies should take into account melanin-based coloration to better predict the evolutionary responses of animals to changing disease risk in human-altered environments.

In situ quantification of individual mRNA transcripts in melanocytes discloses gene regulation of relevance to speciation

Functional validation of candidate genes involved in adaptation and speciation remains challenging. Here, we exemplify the utility of a method quantifying individual mRNA transcripts in revealing the molecular basis of divergence in feather pigment synthesis during early-stage speciation in crows. Using a padlock probe assay combined with rolling circle amplification, we quantified cell-type-specific gene expression in the histological context of growing feather follicles. Expression of Tyrosinase Related Protein 1 (TYRP1), Solute Carrier Family 45 member 2 (SLC45A2) and Hematopoietic Prostaglandin D Synthase (HPGDS) was melanocyte-limited and significantly reduced in follicles from hooded crow, explaining the substantially lower eumelanin content in grey versus black feathers. The central upstream Melanocyte Inducing Transcription Factor (MITF) only showed differential expression specific to melanocytes - a feature not captured by bulk RNA-seq. Overall, this study provides insight into the molecular basis of an evolutionary young transition in pigment synthesis, and demonstrates the power of histologically explicit, statistically substantiated single-cell gene expression quantification for functional genetic inference in natural populations.

Discrete or indiscrete? Redefining the colour polymorphism of the land snail Cepaea nemoralis

Biologists have long tried to describe and name the different phenotypes that make up the shell polymorphism of the land snail Cepaea nemoralis. Traditionally, the view is that the ground colour of the shell is one of a few major colour classes, either yellow, pink or brown, but in practise it is frequently difficult to distinguish the colours, and define different shades of the same colour. To understand whether colour variation is in reality continuous, and to investigate how the variation may be perceived by an avian predator, we applied psychophysical models of colour vision to shell reflectance measures. We found that both achromatic and chromatic variation are indiscrete in Cepaea nemoralis, being continuously distributed over many perceptual units. Nonetheless, clustering analysis based on the density of the distribution did reveal three groups, roughly corresponding to human-perceived yellow, pink and brown shells. We also found large-scale geographic variation in the frequency of these groups across Europe, and some covariance between shell colour and banding patterns. Although further studies are necessary, the observation of continuous variation in colour is intriguing because the traditional theory is that the underlying supergene that determines colour has evolved to prevent phenotypes from “dissolving” into continuous trait distributions. The findings thus have significance for understanding the Cepaea polymorphism, and the nature of the selection that acts upon it, as well as more generally highlighting the need to measure colour objectively in other systems.

Embryonic origin and genetic basis of cave associated phenotypes in the isopod crustacean Asellus aquaticus

Characteristics common to animals living in subterranean environments include the reduction or absence of eyes, lessened pigmentation and enhanced sensory systems. How these characteristics have evolved is poorly understood for the majority of cave dwelling species. In order to understand the evolution of these changes, this study uses an invertebrate model system, the freshwater isopod crustacean, Asellus aquaticus, to examine whether adult differences between cave and surface dwelling individuals first appear during embryonic development. We hypothesized that antennal elaboration, as well as eye reduction and pigment loss, would be apparent during embryonic development. We found that differences in pigmentation, eye formation, and number of segments of antenna II were all present by the end of embryonic development. In addition, we found that cave and surface hatchlings do not significantly differ in the relative size of antenna II and the duration of embryonic development. To investigate whether the regions responsible for eye and pigment differences could be genetically linked to differences in article number, we genotyped F2 hybrids for the four previously mapped genomic regions associated with eye and pigment differences and phenotyped these F2 hybrids for antenna II article number. We found that the region previously known to be responsible for both presence versus absence of pigment and eye size also was significantly associated with article number. Future experiments will address whether pleiotropy and/or genetic linkage play a role in the evolution of cave characteristics in Asellus aquaticus.

Seeing spots: quantifying mother-offspring similarity and assessing fitness consequences of coat pattern traits in a wild population of giraffes (Giraffa camelopardalis)

Polymorphic phenotypes of mammalian coat coloration have been important to the study of genetics and evolution, but less is known about the inheritance and fitness consequences of individual variation in complex coat pattern traits such as spots and stripes. Giraffe coat markings are highly complex and variable and it has been hypothesized that variation in coat patterns most likely affects fitness by camouflaging neonates against visually hunting predators. We quantified complex coat pattern traits of wild Masai giraffes using image analysis software, determined the similarity of spot pattern traits between mother and offspring, and assessed whether variation in spot pattern traits was related to fitness as measured by juvenile survival. The methods we described could comprise a framework for objective quantification of complex mammal coat pattern traits based on photographic coat pattern data. We demonstrated that some characteristics of giraffe coat spot shape were likely to be heritable, as measured by mother-offspring regression. We found significant variation in juvenile survival among phenotypic groups of neonates defined by multivariate clustering based on spot trait measurement variables. We also found significant variation in neonatal survival associated with spot size and shape covariates. Larger spots (smaller number of spots) and irregularly shaped or rounder spots (smaller aspect ratio) were correlated with increased survival. These findings will inform investigations into developmental and genetic architecture of complex mammal coat patterns and their adaptive value.

Toward Understanding the Repeated Occurrence of Associations between Melanin-Based Coloration and Multiple Phenotypes

Melanin is the most widespread pigment in organisms. Melanin-based coloration has been repeatedly observed to be associated with the same traits and in the same direction in different vertebrate and insect species. However, whether any factors that are common to different taxa account for the repeated evolution of melanin-phenotype associations remains unclear. We propose to approach this question from the perspective of convergent and parallel evolution to clarify to what extent different species have evolved the same associations owing to a shared genetic basis and being subjected to similar selective pressures. Our current understanding of the genetic basis of melanin-phenotype associations allows for both convergent and parallel evolution, but this understanding is still limited. Further research is needed to clarify the generality and interdependencies of the different proposed mechanisms (supergenes, pleiotropy based on hormones, or neural crest cells). The general ecological scenarios whereby melanin-based coloration is under selection-protection from ultraviolet radiation, thermoregulation in cold environments, or as a signal of social status-offer a good opportunity to study how melanin-phenotype associations evolve. Reviewing these scenarios shows that some traits associated with melanin-based coloration might be selected together with coloration by also favoring adaptation but that other associated traits might impede adaptation, which may be indicative of genetic constraints. We therefore encourage further research on the relative roles that selection and genetic constraints play in shaping multiple melanin-phenotype associations. Placed into a phylogenetic context, this will help clarify to what extent these associations result from convergent or parallel evolutionary processes and why melanin-phenotype associations are so common across the tree of life.

Admixture mapping in a hybrid zone reveals loci associated with avian feather coloration

Identifying the genetic bases for colour patterns has provided important insights into the control and expression of pigmentation and how these characteristics influence fitness. However, much more is known about the genetic bases for traits based on melanin pigments than for traits based on another major class of pigments, carotenoids. Here, we use natural admixture in a hybrid zone between Audubon's and myrtle warblers (Setophaga coronata auduboni/S. c. coronata) to identify genomic regions associated with both types of pigmentation. Warblers are known for rapid speciation and dramatic differences in plumage. For each of five plumage coloration traits, we found highly significant associations with multiple single-nucleotide polymorphisms (SNPs) across the genome and these were clustered in discrete regions. Regions near significantly associated SNPs were enriched for genes associated with keratin filaments, fibrils that make up feathers. A carotenoid-based trait that differs between the taxa-throat colour-had more than a dozen genomic regions of association. One cluster of SNPs for this trait overlaps the Scavenger Receptor Class F Member 2 (SCARF2) gene. Other scavenger receptors are presumed to be expressed at target tissues and involved in the selective movement of carotenoids into the target cells, making SCARF2 a plausible new candidate for carotenoid processing. In addition, two melanin-based plumage traits-colours of the eye line and eye spot-show very strong associations with a single genomic region mapping to chromosome 20 in the zebra finch. These findings indicate that only a subset of the genomic regions differentiated between these two warblers are associated with the plumage differences between them and demonstrate the utility of reduced-representation genomic scans in hybrid zones.

The current and future state of animal coloration research

Animal colour patterns are a model system for understanding evolution because they are unusually accessible for study and experimental manipulation. This is possible because their functions are readily identifiable. In this final paper of the symposium we provide a diagram of the processes affecting colour patterns and use this to summarize their functions and put the other papers in a broad context. This allows us to identify significant 'holes' in the field that only become obvious when we see the processes affecting colour patterns, and their interactions, as a whole. We make suggestions about new directions of research that will enhance our understanding of both the evolution of colour patterns and visual signalling but also illuminate how the evolution of multiple interacting traits works.This article is part of the themed issue 'Animal coloration: production, perception, function and application'.