Genetics, development and evolution of adaptive pigmentation in vertebrates

Skin transcriptomic and selection signature analyses identify ASIP as a key gene in cattle coat color determination

Objective

Coat color is a complex trait and plays an important role in breed identification. However, information regarding genes associated with coat color in cattle is limited, especially at the skin transcriptome level.

Methods

We investigated the differential expressed genes (DEGs) and genomic selection signal underlying the coat color variation between black and brown cattle breeds. A total of 19 cattle (Brangus, Angus, Simmental, and Guanling) were performed skin transcriptome analysis and 262 cattle (Angus and Simmental) were performed whole genome analysis.

Results

Angus cattle (black coat) had a significantly higher melanin content in both their hair and skin compared to that of Simmental and Guanling cattle (brown coat). Transcriptomic analysis identified 14,118 expressed genes, with principal component analysis and hierarchical clustering revealing clear differences between black and brown cattle. DEGs analysis across four pairwise breed comparisons highlighted 343 downregulated and 54 upregulated genes common to all comparisons, with the ASIP gene (agouti signaling protein) emerging as a key gene linked to melanogenesis. The ASIP expression was several dozen-fold higher in brown cattle than in black cattle, suggesting a crucial role in coat color determination. Path-way enrichment and gene set enrichment analysis (GSEA) identified the "Melanogenesis" pathway as significantly enriched and central to coat color variation. Genes such as FZD10, WNT6, and ASIP showed differential expression patterns that correlated with coat color. Genomic analysis revealed strong selection signals in the ASIP gene region, with several SNPs exhibiting high linkage disequilibrium. Notably, the mutation type was predominant in Simmental cattle, while the reference allele was more common in Angus cattle.

Conclusion

Based on the skin transcriptomic and genomic analyses, we found that ASIP was significantly differential expressed between black and brown cattle breeds and under strong positive selection. These findings provide valuable insights into the genetic basis of coat color variation in cattle and highlight the ASIP gene as a critical determinant of this trait.

A novel, rapidly progressive ataxia due to a spontaneous Myo5a mutation in mice impairs transport proteins and alters mitochondria

Spontaneous mouse mutants have helped define genetic contributions to many phenotypes. Here we report a spontaneous Novel Ataxic Phenotype in mice. Ataxia findings were evident at post-natal day 11 in NAP mice and rapidly worsened, resulting in preweaning lethality. Using genome sequencing and genome-wide mapping, we identified a 3' donor splice variant in exon 14 of Myo5a, encoding an actin-based motor protein. The variant in Myo5a (c.1752g>a) excises exon 14 and ablates MYO5A protein expression, which is implicated in intracellular transport and Griscelli syndrome type I in humans. NAP mice displayed expansion of PAX6-positive cells in the external granule layer of the cerebellum, and mass spectrometry analysis of cerebellar extracts uncovered differentially abundant proteins involved in short-range organelle transport, and specifically proteins implicated with early endosomes. Using cerebellar lysates and primary neurons, we provide evidence for an interaction between MYO5A and ANKFY1, a known effector for the endosomal protein, RAB5A. We also found neurons from NAP mice had elongated mitochondria, linking MYO5A to mitochondrial homeostasis. This allele provides new insight into Myo5a function in developmental neuropathology.

Repeatable Selection on Large Ancestry Blocks in an Avian Hybrid Zone

Hybrid zones create natural tests of genetic incompatibilities by combining loci from 2 species in the same genetic background in the wild, making them useful for identifying loci involved in both intrinsic and ecological (extrinsic) isolation. Two Swainson's thrush subspecies form a hybrid zone in western North America. These coastal and inland subspecies exhibit dramatic differences in migration routes; their hybrids exhibit poor migratory survival, suggesting that ecological incompatibilities maintain this zone. We used a panel of ancestry informative markers to identify repeated patterns of selection and introgression across 4 hybrid populations that span the entire length of the Swainson's thrush hybrid zone. Two repeatable patterns consistent with selection against incompatibilities-steep genomic clines and few transitions between ancestry states-were found in large genetic blocks on chromosomes 1 and 5. The block on chromosome 1 showed evidence for inland subspecies introgression while the block on chromosome 5 exhibited coastal subspecies introgression. Some regions previously associated with migratory phenotypes, including migratory orientation, or exhibiting misexpression between the subspecies exhibited signatures of selection in the hybrid zone. Both selection and introgression across the genome were shaped by genomic structural features and evolutionary history, with stronger selection and reduced introgression in regions of low recombination, high subspecies differentiation, positive selection within the subspecies, and on macrochromosomes. Cumulatively, these results suggest that linkage among loci interacts with divergent selection and past divergent evolution between species to strengthen barriers to gene flow within hybrid zones.

(Limited) Predictability of thermal adaptation in invertebrates

Evolutionary genomic approaches provide powerful tools to understand variation in and evolution of physiological processes. Untargeted genomic or transcriptomic screens can identify functionally annotated candidate genes linked to specific physiological processes, in turn suggesting evolutionary roles for these processes. Such studies often aim to inform modeling of the potential of natural populations to adapt to climate change, but these models are most accurate when evolutionary responses are repeatable, and thus predictable. Here, we synthesize the evolutionary genetic and comparative transcriptomic literature on terrestrial and marine invertebrates to assess whether evolutionary responses to temperature are repeatable within populations, across populations and across species. There is compelling evidence for repeatability, sometimes even across species. However, responses to laboratory selection and geographic variation across thermal gradients appear to be highly idiosyncratic. We also survey whether genetic/transcriptomic studies repeatedly identify candidate genes in three functional groups previously associated with the response to thermal stress: heat shock protein (Hsp) genes, proteolysis genes and immunity genes. Multiple studies across terrestrial and marine species identify candidates included in these gene sets. Yet, each of the gene sets are identified in only a minority of studies. Together, these patterns suggest that there is limited predictability of evolutionary responses to natural selection, including across studies within species. We discuss specific patterns for the candidate gene sets, implications for predictive modeling, and other potential applications of evolutionary genetics in elucidating physiology and gene function. Finally, we discuss limitations of inferences from available evolutionary genetic studies and directions for future research.

Development of in vitro hair pigmentation model using hair follicle organoids

Hair color is formed through a series of processes such as melanin synthesis and storage in melanosomes, transfer from melanocytes, and reception by hair matrix cells in the hair bulb. Because gray hair is caused by the deterioration of a single or multiple of these processes, understanding the mechanisms responsible for these processes is crucial for developing therapeutic strategies. Recently, a robust approach for preparing hair follicle organoids (HFOs) was reported, in which hair follicle morphogenesis, including hair shaft elongation, was tracked in vitro. Here, we investigated whether HFOs could be used to assess genes involved in hair pigmentation. HFOs generated hair follicles and pigmented shafts during the in vitro culturing process. The knockdown of genes associated with melanosome production (Bcl2 and Mitf) and transport (MyoX, PAR2, and Rab11b) significantly increased the number of gray hairs in HFOs. This organoid model may be a promising platform for better understanding hair pigmentation and screening drugs for gray hair.

Functions of thyroid hormone signaling in regulating melanophore, iridophore, erythrophore, and pigment pattern formation in spotted scat (Scatophagus argus)

BACKGROUND

Spotted scat, a marine aquaculture fish, has variable body color development stages during their ontogenesis. However, the regulatory mechanism of body color patterns formation was poorly understood. Thyroid hormones (TH) function as an important endocrine factor in regulating metamorphosis. In this study, exogenous thyroid hormones 3,5,3'-L-triiodothyronine (T3) and its inhibitor thiourea (TU) were used to treat spotted scat juveniles during the metamorphosis stage (from 60 to 90 dpf). The function and molecular mechanism of thyroid hormone signaling in regulating body color patterns formation was revealed, using the micro-observation of pigments cells distribution, colorimetric evaluation and carotenoids concentration measurement by spectrophotometry, and comparative transcriptome analysis.

RESULTS

Spotted scat body color patterns consisted of whole body black color, black bar, black and red spots, and its final pattern was formed through the metamorphosis. When spotted scat were treated with the inhibitor TU to disrupt thyroid hormone signaling, the levels of T3 and T4 were significantly decreased, the melanophores numbers were significantly increased, as well as the expression of genes involved in melanin synthesis and melanophore differentiation (tyr, tyrp1, dct, mitf, pmel, oca2, slc24a5, and erbb3) was significantly increased. Besides, the expression of genes associated with carotenoids and pteridine metabolism (apod, pnpla2, rdh12, stard10, xdh, abca1, retsat, scarb1, rgs2, and gch1) and carotenoids accumulation were stimulated, when thyroid hormone signaling was disrupted by TU. On the contrary, the levels of T3 and T4 were significantly elevated in spotted scat treated with T3, which could weaken the skin redness and reduce the number of black spots and melanophores, as well as the number and diameter of larval erythrophores. Notably, unlike melanophores and erythrophores, the differentiation of iridophore was promoted by thyroid hormones, gene related to iridophore differentiation (fhl2-l, fhl2, ltk, id2a, alx4) and guanine metabolism (gmps, hprt1, ppat, impdh1b) were up-regulated after T3 treatment, but they were down-regulated after TU treatment.

CONCLUSIONS

Above results showed that thyroid hormone signaling might play critical roles in regulation pigments synthesis and deposition, thereby affecting pigment cells (melanophores, iridophores and erythrophores) formation and body color patterns. The mechanisms of hyperthyroid and hypothyroid on different pigment cells development were different. Excess thyroid hormone might impact the rearrangement of melanophore by regulating cell cycle, resulting in the abnormalities of black spots in spotted scat. Meanwhile, the excessed thyroid hormone could reduce the number and diameter of larval erythrophores, as well as weaken the skin redness of juvenile erythrophores, but they were enhanced by the disruption of thyroid hormone. However, the formation of iridophore differentiation and guanine synthesis genes expression were stimulated by thyroid hormones. These findings provide new insights for exploring the formation of body color patterns in fish, and help to elucidate the molecular mechanism of thyroid hormone in regulating pigment cell development and body coloration, and may also contribute to selective breeding of ornamental fish.

Structure, function and formation of the amniote skin pattern

From feather and hair dotted arrays to pigmented stripes and spots, the spatial distribution of skin appendages and colouration often forms visible ornaments crucial for fitness in the coat of birds and mammals. These geometrical motifs are extremely diverse in nature. Yet, phenotypic surveys evidenced common themes in variation: the orientation, appendage-specificity or pigmentation of a given region may be conserved across groups or species. Here, we review naturalist observations of natural variation in the anatomy and ecological function of the skin pattern in amniotes. We then describe several decades of genetics, mathematical modelling and experimental embryology work aiming at understanding the molecular and morphogenetic mechanisms responsible for pattern formation. We discuss how these studies provided evidence that the morphological trends and differences representative of the phenotypic landscape of skin patterns in wild amniote species is rooted in the mechanisms controlling the production of distinct compartments in the embryonic skin.

A sequence of SVA retrotransposon insertions in ASIP shaped human pigmentation

Retrotransposons comprise about 45% of the human genome1, but their contributions to human trait variation and evolution are only beginning to be explored2,3. Here, we find that a sequence of SVA retrotransposon insertions in an early intron of the ASIP (agouti signaling protein) gene has probably shaped human pigmentation several times. In the UK Biobank (n = 169,641), a recent 3.3-kb SVA insertion polymorphism associated strongly with lighter skin pigmentation (0.22 [0.21-0.23] s.d.; P = 2.8 × 10-351) and increased skin cancer risk (odds ratio = 1.23 [1.18-1.27]; P = 1.3 × 10-28), appearing to underlie one of the strongest common genetic influences on these phenotypes within European populations4-6. ASIP expression in skin displayed the same association pattern, with the SVA insertion allele exhibiting 2.2-fold (1.9-2.6) increased expression. This effect had an unusual apparent mechanism: an earlier, nonpolymorphic, human-specific SVA retrotransposon 3.9 kb upstream appeared to have caused ASIP hypofunction by nonproductive splicing, which the new (polymorphic) SVA insertion largely eliminated. Extended haplotype homozygosity indicated that the insertion allele has risen to allele frequencies up to 11% in European populations over the past several thousand years. These results indicate that a sequence of retrotransposon insertions contributed to a species-wide increase, then a local decrease, of human pigmentation.

Distinct genetic origins of eumelanin levels and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. To understand the genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

Genomic sequencing reveals convergent adaptation during experimental evolution in two budding yeast species

Convergent evolution is central in the origins of multicellularity. Identifying the basis for convergent multicellular evolution is challenging because of the diverse evolutionary origins and environments involved. Haploid Kluyveromyces lactis populations evolve multicellularity during selection for increased settling in liquid media. Strong genomic and phenotypic convergence is observed between K. lactis and previously selected S. cerevisiae populations under similar selection, despite their >100-million-year divergence. We find K. lactis multicellularity is conferred by mutations in genes ACE2 or AIM44, with ACE2 being predominant. They are a subset of the six genes involved in the S. cerevisiae multicellularity. Both ACE2 and AIM44 regulate cell division, indicating that the genetic convergence is likely due to conserved cellular replication mechanisms. Complex population dynamics involving multiple ACE2/AIM44 genotypes are found in most K. lactis lineages. The results show common ancestry and natural selection shape convergence while chance and contingency determine the degree of divergence.

Expression patterns of melanin-related genes are linked to crypsis and conspicuousness in a pumpkin toadlet

Colour signals play pivotal roles in different communication systems, and the evolution of these characters has been associated with behavioural ecology, integumentary production processes and perceptual mechanisms of the species involved. Here, we present the first insight into the molecular and histological basis of skin colour polymorphism within a miniaturized species of pumpkin toadlet, potentially representing the lowest size threshold for colour polytypism in tetrapods. Brachycephalus actaeus exhibits a coloration ranging from cryptic green to conspicuous orange skin, and our findings suggest that colour morphs differ in their capability to be detected by potential predators. We also found that the distribution and abundance of chromatophores are variable in the different colour morphs. The expression pattern of coloration related genes was predominantly associated with melanin synthesis (including dct, edn1, mlana, oca2, pmel, slc24a5, tyrp1 and wnt9a). Up-regulation of melanin genes in grey, green and brown skin was associated with higher melanophore abundance than in orange skin, where xanthophores predominate. Our findings provide a significant foundation for comparing and understanding the diverse pathways that contribute to the evolution of pigment production in the skin of amphibians.

The genomics of mimicry: Gene expression throughout development provides insights into convergent and divergent phenotypes in a Müllerian mimicry system

A common goal in evolutionary biology is to discern the mechanisms that produce the astounding diversity of morphologies seen across the tree of life. Aposematic species, those with a conspicuous phenotype coupled with some form of defence, are excellent models to understand the link between vivid colour pattern variations, the natural selection shaping it, and the underlying genetic mechanisms underpinning this variation. Mimicry systems in which species share a conspicuous phenotype can provide an even better model for understanding the mechanisms of colour production in aposematic species, especially if comimics have divergent evolutionary histories. Here we investigate the genetic mechanisms by which mimicry is produced in poison frogs. We assembled a 6.02-Gbp genome with a contig N50 of 310 Kbp, a scaffold N50 of 390 Kbp and 85% of expected tetrapod genes. We leveraged this genome to conduct gene expression analyses throughout development of four colour morphs of Ranitomeya imitator and two colour morphs from both R. fantastica and R. variabilis which R. imitator mimics. We identified a large number of pigmentation and patterning genes differentially expressed throughout development, many of them related to melanophores/melanin, iridophore development and guanine synthesis. We also identify the pteridine synthesis pathway (including genes such as qdpr and xdh) as a key driver of the variation in colour between morphs of these species, and identify several plausible candidates for colouration in vertebrates (e.g. cd36, ep-cadherin and perlwapin). Finally, we hypothesise that keratin genes (e.g. krt8) are important for producing different structural colours within these frogs.

Colonial-driven extinction of the blue antelope despite genomic adaptation to low population size

Low genomic diversity is generally indicative of small population size and is considered detrimental by decreasing long-term adaptability.1,2,3,4,5,6 Moreover, small population size may promote gene flow with congeners and outbreeding depression.7,8,9,10,11,12,13 Here, we examine the connection between habitat availability, effective population size (Ne), and extinction by generating a 40× nuclear genome from the extinct blue antelope (Hippotragus leucophaeus). Historically endemic to the relatively small Cape Floristic Region in southernmost Africa,14,15 populations were thought to have expanded and contracted across glacial-interglacial cycles, tracking suitable habitat.16,17,18 However, we found long-term low Ne, unaffected by glacial cycles, suggesting persistence with low genomic diversity for many millennia prior to extinction in ∼AD 1800. A lack of inbreeding, alongside high levels of genetic purging, suggests adaptation to this long-term low Ne and that human impacts during the colonial era (e.g., hunting and landscape transformation), rather than longer-term ecological processes, were central to its extinction. Phylogenomic analyses uncovered gene flow between roan (H. equinus) and blue antelope, as well as between roan and sable antelope (H. niger), approximately at the time of divergence of blue and sable antelope (∼1.9 Ma). Finally, we identified the LYST and ASIP genes as candidates for the eponymous bluish pelt color of the blue antelope. Our results revise numerous aspects of our understanding of the interplay between genomic diversity and evolutionary history and provide the resources for uncovering the genetic basis of this extinct species' unique traits.

Evolution and genetic architecture of sex-limited polymorphism in cuckoos

Sex-limited polymorphism has evolved in many species including our own. Yet, we lack a detailed understanding of the underlying genetic variation and evolutionary processes at work. The brood parasitic common cuckoo (Cuculus canorus) is a prime example of female-limited color polymorphism, where adult males are monochromatic gray and females exhibit either gray or rufous plumage. This polymorphism has been hypothesized to be governed by negative frequency-dependent selection whereby the rarer female morph is protected against harassment by males or from mobbing by parasitized host species. Here, we show that female plumage dichromatism maps to the female-restricted genome. We further demonstrate that, consistent with balancing selection, ancestry of the rufous phenotype is shared with the likewise female dichromatic sister species, the oriental cuckoo (Cuculus optatus). This study shows that sex-specific polymorphism in trait variation can be resolved by genetic variation residing on a sex-limited chromosome and be maintained across species boundaries.

Exploring the origin of a unique mutant allele in twin-tail goldfish using CRISPR/Cas9 mutants

Artificial selection has been widely applied to genetically fix rare phenotypic features in ornamental domesticated animals. For many of these animals, the mutated loci and alleles underlying rare phenotypes are known. However, few studies have explored whether these rare genetic mutations might have been fixed due to competition among related mutated alleles or if the fixation occurred due to contingent stochastic events. Here, we performed genetic crossing with twin-tail ornamental goldfish and CRISPR/Cas9-mutated goldfish to investigate why only a single mutated allele-chdS with a E127X stop codon (also called chdAE127X)-gives rise to the twin-tail phenotype in the modern domesticated goldfish population. Two closely related chdS mutants were generated with CRISPR/Cas9 and compared with the E127X allele in F2 and F3 generations. Both of the CRISPR/Cas9-generated alleles were equivalent to the E127X allele in terms of penetrance/expressivity of the twin-tail phenotype and viability of carriers. These findings indicate that multiple truncating mutations could have produced viable twin-tail goldfish. Therefore, the absence of polymorphic alleles for the twin-tail phenotype in modern goldfish likely stems from stochastic elimination or a lack of competing alleles in the common ancestor. Our study is the first experimental comparison of a singular domestication-derived allele with CRISPR/Cas9-generated alleles to understand how genetic fixation of a unique genotype and phenotype may have occurred. Thus, our work may provide a conceptual framework for future investigations of rare evolutionary events in domesticated animals.

Yosemite toad (Anaxyrus canorus) transcriptome reveals interplay between speciation genes and adaptive introgression

Genomes are heterogeneous during the early stages of speciation, with small 'islands' of DNA appearing to reflect strong adaptive differences, surrounded by vast seas of relative homogeneity. As species diverge, secondary contact zones between them can act as an interface and selectively filter through advantageous alleles of hybrid origin. Such introgression is another important adaptive process, one that allows beneficial mosaics of recombinant DNA ('rivers') to flow from one species into another. Although genomic islands of divergence appear to be associated with reproductive isolation, and genomic rivers form by adaptive introgression, it is unknown whether islands and rivers tend to be the same or different loci. We examined three replicate secondary contact zones for the Yosemite toad (Anaxyrus canorus) using two genomic data sets and a morphometric data set to answer the questions: (1) How predictably different are islands and rivers, both in terms of genomic location and gene function? (2) Are the adaptive genetic trait loci underlying tadpole growth and development reliably islands, rivers or neither? We found that island and river loci have significant overlap within a contact zone, suggesting that some loci are first islands, and later are predictably converted into rivers. However, gene ontology enrichment analysis showed strong overlap in gene function unique to all island loci, suggesting predictability in overall gene pathways for islands. Genome-wide association study outliers for tadpole development included LPIN3, a lipid metabolism gene potentially involved in climate change adaptation, that is island-like for all three contact zones, but also appears to be introgressing (as a river) across one zone. Taken together, our results suggest that adaptive divergence and introgression may be more complementary forces than currently appreciated.

Differential expression and alternative splicing analyses of multiple tissues reveal albinism-associated genes in the Wels catfish (Silurus glanis)

Albinism is a widespread departure from a typical body colouration due to altered melanin production. The Wels catfish (Silurus glanis) is among the largest freshwater fish species in the world, and albino individuals occur both in the wild and in aquaculture. Here, we performed transcriptome-wide analysis of albino and normally pigmented S. glanis using four tissues (skin, dorsal fin, whole eye and liver) to identify genes associated with albinism by exploring patterns of differential expression (DE) and differential alternative splicing (DAS). Multi-tissue analyses revealed a large number of genes in skin (n = 1355) and fin (n = 614) tissue associated with the albino phenotype in S. glanis, while the number of DE genes in eye and liver tissues was lower (n = 188, n = 189, respectively). Several DE genes across multiple tissues were detected as the most promising candidates (e.g., hsp4, hsp90b1, raph1, uqcrfs1, adcy-family and wnt-family) potentially causally linked to the albino phenotype in Wels catfish. Moreover, our findings supported earlier observations of physiological differences between albino and normally pigmented individuals, particularly in energy metabolism and immune response. In contrast, there were only a few pigmentation-related genes observed among DAS genes (4 in skin, 2 in fin), the overlap between DAS and DE genes was low (n = 25) and did not include known pigmentation-related genes. This suggests that DAS and DE in Wels catfish are, to a large extent, independent processes, and the observed alternative splicing cases are probably not causally linked with albinism in S. glanis. This work provides the first transcriptome-wide multi-tissue insights into the albinism of Wels catfish and serves as a valuable resource for further understanding the genetic mechanisms of pigmentation in fish.

Differential expression of ASIP transcripts reveals genetic mechanism underpinning black-tail independence from body plumage in yellow-bodied chickens

The genetic foundation of chicken body plumage color has been extensively studied. However, little attention has been paid to the inheritance patterns and molecular mechanisms underlying the formation of distal feather colors (tail and wingtip). Differences in these colors are common; for example, the Chinese Huiyang Beard chicken has black tail feathers, but yellow body plumage. Here, the hybrid offspring of Huiyang Beard and White Leghorn chickens were used to study the inheritance patterns of tail-feather color. The expression levels of pigment genes in differently colored feather follicles were analyzed using quantitative real-time PCR. The results showed that genetic regulation of tail-feather color was independent of body-plumage color. The Dominant White locus inhibited eumelanin synthesis in tail feathers without affecting the formation of yellow body plumage, whereas the Silver locus had the opposite effect. The expression of agouti signaling protein (ASIP) gene class 1 transcripts was significantly lower in black tail-feather follicles than in yellow body follicles, whereas tyrosinase-related protein 1 (TYRP1) gene expression was significantly higher in black tail feathers. These differentially expressed genes were confirmed to exert an effect on eumelanin and pheomelanin formation in feathers, thus influencing the regulation of chicken tail-feather color. In conclusion, this study lays the foundation for further research on the genetic mechanisms of regional differences in feather color, contributing to a better understanding of plumage pigmentation in chickens.

Taking a color photo: A homozygous 25-bp deletion in Bace2 may cause brown-and-white coat color in giant pandas

Brown-and-white giant pandas (hereafter brown pandas) are distinct coat color mutants found exclusively in the Qinling Mountains, Shaanxi, China. However, its genetic mechanism has remained unclear since their discovery in 1985. Here, we identified the genetic basis for this coat color variation using a combination of field ecological data, population genomic data, and a CRISPR-Cas9 knockout mouse model. We de novo assembled a long-read-based giant panda genome and resequenced the genomes of 35 giant pandas, including two brown pandas and two family trios associated with a brown panda. We identified a homozygous 25-bp deletion in the first exon of Bace2, a gene encoding amyloid precursor protein cleaving enzyme, as the most likely genetic basis for brown-and-white coat color. This deletion was further validated using PCR and Sanger sequencing of another 192 black giant pandas and CRISPR-Cas9 edited knockout mice. Our investigation revealed that this mutation reduced the number and size of melanosomes of the hairs in knockout mice and possibly in the brown panda, further leading to the hypopigmentation. These findings provide unique insights into the genetic basis of coat color variation in wild animals.

Single-cell profiling of MC1R-inhibited melanocytes

The human red hair color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to define a MC1R-inhibited Gene Signature (MiGS) comprising a large set of previously unidentified genes which may be implicated in melanogenesis and oncogenic transformation. We show that one of the candidate MiGS genes, TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into further mechanisms by which melanocytes with reduced MC1R signaling may regulate pigmentation and offer new candidates of study toward understanding how individuals with the RHC phenotype are predisposed to melanoma.

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 Genetic Basis of Melanism in Abert's Squirrel (Sciurus aberti)

Melanism is widespread in different taxa and has been hypothesized to provide adaptive benefits in certain environments. Melanism is typically caused by mutations in one of two regulatory genes: the Melanocortin 1 Receptor (MC1R) or the Agouti Signaling Protein (ASIP). Melanism has repeatedly evolved among tree squirrels and their relatives (tribe Sciurini) in at least 12 different species based on our review of the literature. The causal mutations for melanism have been characterized in two species so far. This study examines Abert's Squirrel (Sciurus aberti), which has a melanistic morph whose genetic basis has not yet been established. We sequenced the MC1R and ASIP genes for five wild-type and seven melanistic S. aberti individuals to search for melanism-associated mutations. A novel single base pair mutation in the ASIP gene, unique to S. aberti, was found to be associated with melanism in the species, indicating that melanism in S. aberti evolved independently from other tree squirrels and thus represents an example of convergent evolution. The independent evolution of melanism in this species suggests that there is an adaptive advantage to the melanistic phenotype. The geographic range and habitat of S. aberti suggest possible benefits associated with thermoregulation, post-forest-fire camouflage, or other untested hypotheses.

Diet influences resource allocation in chemical defence but not melanin synthesis in an aposematic moth

For animals that synthesise their chemical compounds de novo, resources, particularly proteins, can influence investment in chemical defences and nitrogen-based wing colouration such as melanin. Competing for the same resources often leads to trade-offs in resource allocation. We manipulated protein availability in the larval diet of the wood tiger moth, Arctia plantaginis, to test how early life resource availability influences relevant life history traits, melanin production and chemical defences. We expected higher dietary protein to result in more effective chemical defences in adult moths and a higher amount of melanin in the wings. According to the resource allocation hypothesis, we also expected individuals with less melanin to have more resources to allocate to chemical defences. We found that protein-deprived moths had a slower larval development, and their chemical defences were less unpalatable for bird predators, but the expression of melanin in their wings did not differ from that of moths raised on a high-protein diet. The amount of melanin in the wings, however, unexpectedly correlated positively with chemical defences. Our findings demonstrate that the resources available in early life have an important role in the efficacy of chemical defences, but melanin-based warning colours are less sensitive to resource variability than other fitness-related traits.

Genetic and behavioural factors affecting interpopulation colour pattern variation in two congeneric chameleon species

We conducted a study on interpopulation variation of colour patterns in two congeneric chameleon species, which have an analogous life history. Both species are able to rapidly change colour pattern, and their context-dependent colour patterns often vary across a wide geographical range. Specifically, we tested four hypotheses that can explain the observed interpopulation variation of colour patterns by a series of behavioural field trials where the colour patterns of individuals were recorded and later analysed by a deep neural network algorithm. We used redundancy analysis to relate genetic, spectral and behavioural predictors to interpopulation colour pattern distance. Our results showed that both isolation by distance (IBD) and alternative mating tactics were significant predictors for interpopulation colour pattern variation in Chamaeleo chamaeleon males. By contrast, in Chamaeleo dilepis, the interpopulation colour pattern variation was largely explained by IBD, and evidence for alternative mating tactics was absent. In both chameleon species, the environmental colours showed no evidence of influencing chameleon interpopulation colour pattern variation, regardless of sex or behavioural context. This contrasting finding suggests that interpopulation context-dependent colour pattern variations in each species are maintained under a different set of selective pressures or circumstances.

Involvement of miR-495 in the skin pigmentation of rainbow trout (Oncorhynchus mykiss) through the regulation of mc1r

MicroRNAs (miRNAs) play crucial roles in skin pigmentation in animals. Rainbow trout (Oncorhynchus mykiss) is a key economic fish species worldwide, and skin color directly affects its economic value. However, the functions of miRNAs in rainbow trout skin pigmentation remain largely unknown. Herein, we overexpressed and silenced miR-495 in vitro and in vivo to investigate its functions. The analysis of spatial and temporal expression patterns suggested that miR-495 is a potential regulator during the process of skin pigmentation. In vitro, mc1r was validated as a direct target for miR-495 by dual-luciferase reporter assay, and overexpression of miR-495 significantly inhibited mc1r expression; in contrast, mc1r and its downstream gene mitf levels were markedly upregulated by decreased miR-495. In vivo, overexpressed miR-495 by injecting agomiR-495 led to a substantial decrease in the expression of mc1r and mitf in dorsal skin and liver, while the opposite results were obtained after miR-495 silencing by antagomiR-495. These findings suggested that miR-495 can target mc1r to regulate rainbow trout skin pigmentation, which provide a potential basis for using miRNAs as target drugs to treat pigmentation disorders and melanoma.

Ultrastructural analysis of throat dermal tissue and chromatophore components in the threespine stickleback (Gasterosteus aculeatus)

The threespine stickleback (Gasterosteus aculeatus) is an important model for studying the evolution of nuptial coloration, but histological analyses of color are largely lacking. Previous analyses of one nuptial coloration trait, orange-red coloration along the body, have indicated carotenoids are the main pigment producing this color. In addition, recent gene expression studies found variation in the correlates of throat coloration between the sexes and between populations, raising the possibility of variation in the mechanisms underlying superficially similar coloration. We used transmission electron microscopy (TEM) to investigate the histological correlates of color in the throat dermal tissue of threespine stickleback from Western North America, within and between sexes, populations, and ecotypes. Ultrastructural analysis revealed carotenoid-containing erythrophores to be the main chromatophore component associated with orange-red coloration in both males and females across populations. In individuals where some darkening of the throat tissue was present, with no obvious orange-red coloration, erythrophores were not detected. Melanophore presence was more population-specific in expression, including being the only chromatophore component detected in a population of darker fish. We found no dermal chromatophore units within colorless throat tissue. This work confirms the importance of carotenoids and the erythrophore in producing orange-red coloration across sexes, as well as melanin within the melanophore in producing darkened coloration, but does not reveal broad histological differences among populations with similar coloration.

Two High-Quality Cygnus Genome Assemblies Reveal Genomic Variations Associated with Plumage Color

As an exemplary model for examining molecular mechanisms responsible for extreme phenotypic variations, plumage color has garnered significant interest. The Cygnus genus features two species, Cygnus olor and Cygnus atratus, that exhibit striking disparities in plumage color. However, the molecular foundation for this differentiation has remained elusive. Herein, we present two high-quality genomes for C. olor and C. atratus, procured using the Illumina and Nanopore technologies. The assembled genome of C. olor was 1.12 Gb in size with a contig N50 of 26.82 Mb, while its counterpart was 1.13 Gb in size with a contig N50 of 21.91 Mb. A comparative analysis unveiled three genes (TYR, SLC45A2, and SLC7A11) with structural variants in the melanogenic pathway. Notably, we also identified a novel gene, PWWP domain containing 2A (PWWP2A), that is related to plumage color, for the first time. Using targeted gene modification analysis, we demonstrated the potential genetic effect of the PWWP2A variant on pigment gene expression and melanin production. Finally, our findings offer insight into the intricate pattern of pigmentation and the role of polygenes in birds. Furthermore, these two high-quality genome references provide a comprehensive resource and perspective for comparative functional and genetic studies of evolution within the Cygnus genus.

Colour polymorphism associated with a gene duplication in male wood tiger moths

Colour is often used as an aposematic warning signal, with predator learning expected to lead to a single colour pattern within a population. However, there are many puzzling cases where aposematic signals are also polymorphic. The wood tiger moth, Arctia plantaginis, displays bright hindwing colours associated with unpalatability, and males have discrete colour morphs which vary in frequency between localities. In Finland, both white and yellow morphs can be found, and these colour morphs also differ in behavioural and life-history traits. Here, we show that male colour is linked to an extra copy of a yellow family gene that is only present in the white morphs. This white-specific duplication, which we name valkea, is highly upregulated during wing development. CRISPR targeting valkea resulted in editing of both valkea and its paralog, yellow-e, and led to the production of yellow wings. We also characterise the pigments responsible for yellow, white, and black colouration, showing that yellow is partly produced by pheomelanins, while black is dopamine-derived eumelanin. Our results add to a growing number of studies on the genetic architecture of complex and seemingly paradoxical polymorphisms, and the role of gene duplications and structural variation in adaptive evolution.

Pigment Identification and Gene Expression Analysis during Erythrophore Development in Spotted Scat (Scatophagus argus) Larvae

Red coloration is considered an economically important trait in some fish species, including spotted scat, a marine aquaculture fish. Erythrophores are gradually covered by melanophores from the embryonic stage. Despite studies of black spot formation and melanophore coloration in the species, little is known about erythrophore development, which is responsible for red coloration. 1-phenyl 2-thiourea (PTU) is a tyrosinase inhibitor commonly used to inhibit melanogenesis and contribute to the visualization of embryonic development. In this study, spotted scat embryos were treated with 0.003% PTU from 0 to 72 h post fertilization (hpf) to inhibit melanin. Erythrophores were clearly observed during the embryonic stage from 14 to 72 hpf, showing an initial increase (14 to 36 hpf), followed by a gradual decrease (36 to 72 hpf). The number and size of erythrophores at 36 hpf were larger than those at 24 and 72 hpf. At 36 hpf, LC-MS and absorbance spectrophotometry revealed that the carotenoid content was eight times higher than the pteridine content, and β-carotene and lutein were the main pigments related to red coloration in spotted scat larvae. Compared with their expression in the normal hatching group, rlbp1b, rbp1.1, and rpe65a related to retinol metabolism and soat2 and apoa1 related to steroid hormone biosynthesis and steroid biosynthesis were significantly up-regulated in the PTU group, and rh2 associated with phototransduction was significantly down-regulated. By qRT-PCR, the expression levels of genes involved in carotenoid metabolism (scarb1, plin6, plin2, apoda, bco1, and rep65a), pteridine synthesis (gch2), and chromatophore differentiation (slc2a15b and csf1ra) were significantly higher at 36 hpf than at 24 hpf and 72 hpf, except for bco1. These gene expression profiles were consistent with the developmental changes of erythrophores. These findings provide insights into pigment cell differentiation and gene function in the regulation of red coloration and contribute to selective breeding programs for ornamental aquatic animals.

Uncovering the genetic architecture of parallel evolution

Identifying the genetic architecture underlying adaptive traits is exceptionally challenging in natural populations. This is because associations between traits not only mask the targets of selection but also create correlated patterns of genomic divergence that hinder our ability to isolate causal genetic effects. Here, we examine the repeated evolution of components of the auxin pathway that have contributed to the replicated loss of gravitropism (i.e. the ability of a plant to bend in response to gravity) in multiple populations of the Senecio lautus species complex in Australia. We use a powerful approach which combines parallel population genomics with association mapping in a Multiparent Advanced Generation Inter-Cross (MAGIC) population to break down genetic and trait correlations to reveal how adaptive traits evolve during replicated evolution. We sequenced auxin and shoot gravitropism-related gene regions in 80 individuals from six natural populations (three parallel divergence events) and 133 individuals from a MAGIC population derived from two of the recently diverged natural populations. We show that artificial tail selection on gravitropism in the MAGIC population recreates patterns of parallel divergence in the auxin pathway in the natural populations. We reveal a set of 55 auxin gene regions that have evolved repeatedly during the evolution of the species, of which 50 are directly associated with gravitropism divergence in the MAGIC population. Our work creates a strong link between patterns of genomic divergence and trait variation contributing to replicated evolution by natural selection, paving the way to understand the origin and maintenance of adaptations in natural populations.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

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.

Genealogical characterization of regional populations and dorsal coat color variation in the house mouse Mus musculus from Asia based on haplotype structure analysis of a gene-rich region harboring Mc1r

We analyzed 196 haplotype sequences from a gene-rich region (250 kb) that includes Mc1r, a gene involved in coat color regulation, to gain insight into the evolution of coat color variation in subspecies of the house mouse Mus musculus. Phylogenetic networks revealed haplotype groups from the major subspecies of M. m. castaneus (CAS), M. m. domesticus (DOM), and M. m. musculus (MUS). Using haplotype sequences assigned to each of CAS and MUS through phylogenetic analysis, we proposed migration routes associated with prehistoric humans from west to east across Eurasia. Comparing nucleotide diversity among subspecies-specific haplotypes in different geographic areas showed a marked reduction during migration, particularly in MUS-derived haplotypes from Korea and Japan, suggesting intensive population bottlenecks during migration. We found that a C>T polymorphism at site 302 (c.302C>T) in the Mc1r coding region correlated with a darkening of dorsal fur color in both CAS and MUS. However, C/C homozygous mice in MUS showed marked variation in lightness, indicating the possibility of another genetic determinant that affects the lightness of dorsal fur color. Detailed sequence comparisons of haplotypes revealed that short fragments assigned to DOM were embedded in CAS-assigned fragments, indicating ancient introgression between subspecies. The estimated age of c.302C>T also supports the hypothesis that genetic interaction between subspecies occurred in ancient times. This suggests that the genome of M. musculus evolved through gene flow between subspecies over an extended period before the movement of the species in conjunction with prehistoric humans.

ebony underpins Batesian mimicry in melanic stoneflies

The evolution of Batesian mimicry - whereby harmless species avoid predation through their resemblance to harmful species - has long intrigued biologists. In rare cases, Batesian mimicry is linked to intraspecific colour variation, in which only some individuals within a population resemble a noxious 'model'. Here, we assess intraspecific colour variation within a widespread New Zealand stonefly, wherein highly melanized individuals of Zelandoperla closely resemble a chemically defended aposematic stonefly, Austroperla cyrene. We assess convergence in the colour pattern of these two species, compare their relative palatability to predators, and use genome-wide association mapping to assess the genetic basis of this resemblance. Our analysis reveals that melanized Zelandoperla overlap significantly with Austroperla in colour space but are significantly more palatable to predators, implying that they are indeed Batesian mimics. Analysis of 194,773 genome-wide SNPs reveals an outlier locus (ebony) strongly differentiating melanic versus non-melanic Zelandoperla. Genotyping of 338 specimens from a single Zelandoperla population indicates that ebony explains nearly 70% of the observed variance in melanism. As ebony has a well-documented role in insect melanin biosynthesis, our findings indicate this locus has a conserved function across deeply divergent hexapod lineages. Distributional records suggest a link between the occurrence of melanic Zelandoperla and the forested ecosystems where the model Austroperla is abundant, suggesting the potential for adaptive shifts in this system underpinned by environmental change.

Effects of Dietary Cystine and Tyrosine Supplementation on Melanin Synthesis in the Pacific Oyster (Crassostrea gigas)

Melanogenesis is a multistep process to produce melanin for dark pigmentation in skin coloration. Previous studies in vertebrates demonstrated that cystine and tyrosine amino acids are involved in the melanin synthesis. However, very little is known about the melanogenesis in bivalve. In this study, cystine supplementation for 30 days significantly upregulated the expression of CgB-aat1, CgCbs and CgTyr and pheomelanin content in the Pacific oyster Crassostrea gigas. Transmission electron microscope (TEM) results revealed more melanosomes in the connective tissue and melanin granules were secreted in epithelium of mantle. In contrast, tyrosine supplementation had no clear effect on melanogenesis except the gene expression changes of CgB-aat1 and CgCbs. In addition, prolonged supplementation of cystine or tyrosine for 60 days had a negative impact on melanogenesis. Indeed, after 60 days, expression of most of the melanin synthesis-related genes under study was decreased, and melanin content was significantly reduced, indicating that cystine and tyrosine might inhibit production of eumelanin and pheomelanin, respectively. In addition, in vitro analysis using primary cell culture from mantle tissue indicated that incubation with cystine, tyrosine, or B-AAT1 polypeptide, CBS/TYR recombinant proteins induced the increase of CgB-aat1 and CgCbs expression in a dose-dependent manner, suggesting the presence of a regulatory network in response to cystine and tyrosine amino acids intakes in pheomelanin synthesis-related gene expression. Taken together, these data indicate that cystine-CgB-aat1-CgCbs-CgTyr axis is a potential regulator of the pheomelanin biosynthesis pathway, and thus plays an important role in the mantle pigmentation in C. gigas. This work provides a new clue for selective cultivation of oyster strains with specific shell colors in bivalve breeding.

Analysis of the genetic loci of pigment pattern evolution in vertebrates

Vertebrate pigmentation patterns are amongst the best characterised model systems for studying the genetic basis of adaptive evolution. The wealth of available data on the genetic basis for pigmentation evolution allows for analysis of trends and quantitative testing of evolutionary hypotheses. We employed Gephebase, a database of genetic variants associated with natural and domesticated trait variation, to examine trends in how cis-regulatory and coding mutations contribute to vertebrate pigmentation phenotypes, as well as factors that favour one mutation type over the other. We found that studies with lower ascertainment bias identified higher proportions of cis-regulatory mutations, and that cis-regulatory mutations were more common amongst animals harbouring a higher number of pigment cell classes. We classified pigmentation traits firstly according to their physiological basis and secondly according to whether they affect colour or pattern, and identified that carotenoid-based pigmentation and variation in pattern boundaries are preferentially associated with cis-regulatory change. We also classified genes according to their developmental, cellular, and molecular functions. We found a greater proportion of cis-regulatory mutations in genes implicated in upstream developmental processes compared to those involved in downstream cellular functions, and that ligands were associated with a higher proportion of cis-regulatory mutations than their respective receptors. Based on these trends, we discuss future directions for research in vertebrate pigmentation evolution.

Single-Cell Transcriptome and Pigment Biochemistry Analysis Reveals the Potential for the High Nutritional and Medicinal Value of Purple Sea Cucumbers

The sea cucumber Apostichopus japonicus has important nutritional and medicinal value. Unfortunately, we know little of the source of active chemicals in this animal, but the plentiful pigments of these animals are thought to function in intriguing ways for translation into clinical and food chemistry usage. Here, we found key cell groups with the gene activity predicted for the color morphology of sea cucumber body using single-cell RNA-seq. We refer to these cell populations as melanocytes and quinocytes, which are responsible for the synthesis of melanin and quinone pigments, respectively. We integrated analysis of pigment biochemistry with the transcript profiles to illuminate the molecular mechanisms regulating distinct pigment formation in echinoderms. In concert with the correlated pigment analysis from each color morph, this study expands our understanding of medically important pigment production, as well as the genetic mechanisms for color morphs, and provides deep datasets for exploring advancements in the fields of bioactives and nutraceuticals.

An allelic series at the EDNRB2 locus controls diverse piebalding patterns in the domestic pigeon

Variation in pigment patterns within and among vertebrate species reflects underlying changes in cell migration and function that can impact health, reproductive success, and survival. The domestic pigeon (Columba livia) is an exceptional model for understanding the genetic changes that give rise to diverse pigment patterns, as selective breeding has given rise to hundreds of breeds with extensive variation in plumage color and pattern. Here, we map the genetic architecture of a suite of pigmentation phenotypes known as piebalding. Piebalding is characterized by patches of pigmented and non-pigmented feathers, and these plumage patterns are often breed-specific and stable across generations. Using a combination of quantitative trait locus mapping in F2 laboratory crosses and genome-wide association analysis, we identify a locus associated with piebalding across many pigeon breeds. This shared locus harbors a candidate gene, EDNRB2, that is a known regulator of pigment cell migration, proliferation, and survival. We discover multiple distinct haplotypes at the EDNRB2 locus in piebald pigeons, which include a mix of protein-coding, noncoding, and structural variants that are associated with depigmentation in specific plumage regions. These results identify a role for EDNRB2 in pigment patterning in the domestic pigeon, and highlight how repeated selection at a single locus can generate a diverse array of stable and heritable pigment patterns.

A unique single nucleotide polymorphism in Agouti Signalling Protein (ASIP) gene changes coat colour of Sri Lankan leopard (Panthera pardus kotiya) to dark black

The Sri Lankan leopard (Panthera pardus kotiya) is an endangered subspecies restricted to isolated and fragmented populations in Sri Lanka. Among them, melanistic leopards have been recorded on a few occasions. Literature suggests the evolution of melanism several times in the Felidae family, with three species having distinct mutations. Nevertheless, the mutations or other variations in the remaining species, including Sri Lankan melanistic leopard, are unknown. We used reference-based assembled nuclear genomes of Sri Lankan wild type and melanistic leopards and de novo assembled mitogenomes of the same to investigate the genetic basis, adaptive significance, and evolutionary history of the Sri Lankan melanistic leopard. Interestingly, we identified a single nucleotide polymorphism in exon-4 Sri Lankan melanistic leopard, which may completely ablate Agouti Signalling Protein (ASIP) function. The wild type leopards in Sri Lanka did not carry this mutation, suggesting the cause for the occurrence of melanistic leopords in the population. Comparative analysis of existing genomic data in the literature suggests it as a P. p. kotiya specific mutation and a novel mutation in the ASIP-gene of the Felidae family, contributing to naturally occurring colour polymorphism. Our data suggested the coalescence time of Sri Lankan leopards at ~0.5 million years, sisters to the Panthera pardus lineage. The genetic diversity was low in Sri Lankan leopards. Further, the P. p. kotiya melanistic leopard is a different morphotype of the P. p. kotiya wildtype leopard resulting from the mutation in the ASIP-gene. The ability of black leopards to camouflage, along with the likelihood of recurrence and transfer to future generations, suggests that this rare mutation could be environment-adaptable.

Distinct genetic origins of eumelanin intensity and barring patterns in cichlid fishes

Pigment patterns are incredibly diverse across vertebrates and are shaped by multiple selective pressures from predator avoidance to mate choice. A common pattern across fishes, but for which we know little about the underlying mechanisms, is repeated melanic vertical bars. In order to understand genetic factors that modify the level or pattern of vertical barring, we generated a genetic cross of 322 F2 hybrids between two cichlid species with distinct barring patterns, Aulonocara koningsi and Metriaclima mbenjii. We identify 48 significant quantitative trait loci that underlie a series of seven phenotypes related to the relative pigmentation intensity, and four traits related to patterning of the vertical bars. We find that genomic regions that generate variation in the level of eumelanin produced are largely independent of those that control the spacing of vertical bars. Candidate genes within these intervals include novel genes and those newly-associated with vertical bars, which could affect melanophore survival, fate decisions, pigment biosynthesis, and pigment distribution. Together, this work provides insights into the regulation of pigment diversity, with direct implications for an animal's fitness and the speciation process.

Identification of a new Corin atrial natriuretic peptide-converting enzyme substrate: Agouti-signaling protein (ASIP)

Corin is a transmembrane tethered enzyme best known for processing the hormone atrial natriuretic peptide (ANP) in cardiomyocytes to control electrolyte balance and blood pressure. Loss of function mutations in Corin prevent ANP processing and lead to hypertension. Curiously, Corin loss of function variants also result in lighter coat color pigmentation in multiple species. Corin pigmentation effects are dependent on a functional Agouti locus encoding the agouti-signaling protein (ASIP) based on a genetic interaction. However, the nature of this conserved role of Corin has not been defined. Here we report that ASIP is a direct proteolytic substrate of the Corin enzyme.

Bear essentials

Strangely, American black bears come in many colours. New work by Puckett et al. shows that a missense alteration in the gene encoding tyrosinase-related protein 1 (TYRP1) likely interferes with melanin synthesis and is responsible for the cinnamon colour variant in the southwest USA. However, the adaptive significance of colour polymorphisms in this large carnivore remains opaque.

Loss of MC1R signaling implicates TBX3 in pheomelanogenesis and melanoma predisposition

The human Red Hair Color (RHC) trait is caused by increased pheomelanin (red-yellow) and reduced eumelanin (black-brown) pigment in skin and hair due to diminished melanocortin 1 receptor (MC1R) function. In addition, individuals harboring the RHC trait are predisposed to melanoma development. While MC1R variants have been established as causative of RHC and are a well-defined risk factor for melanoma, it remains unclear mechanistically why decreased MC1R signaling alters pigmentation and increases melanoma susceptibility. Here, we use single-cell RNA-sequencing (scRNA-seq) of melanocytes isolated from RHC mouse models to reveal a Pheomelanin Gene Signature (PGS) comprising genes implicated in melanogenesis and oncogenic transformation. We show that TBX3, a well-known anti-senescence transcription factor implicated in melanoma progression, is part of the PGS and binds both E-box and T-box elements to regulate genes associated with melanogenesis and senescence bypass. Our results provide key insights into mechanisms by which MC1R signaling regulates pigmentation and how individuals with the RHC phenotype are predisposed to melanoma.

Molecular mechanisms of adaptive evolution in wild animals and plants

Wild animals and plants have developed a variety of adaptive traits driven by adaptive evolution, an important strategy for species survival and persistence. Uncovering the molecular mechanisms of adaptive evolution is the key to understanding species diversification, phenotypic convergence, and inter-species interaction. As the genome sequences of more and more non-model organisms are becoming available, the focus of studies on molecular mechanisms of adaptive evolution has shifted from the candidate gene method to genetic mapping based on genome-wide scanning. In this study, we reviewed the latest research advances in wild animals and plants, focusing on adaptive traits, convergent evolution, and coevolution. Firstly, we focused on the adaptive evolution of morphological, behavioral, and physiological traits. Secondly, we reviewed the phenotypic convergences of life history traits and responding to environmental pressures, and the underlying molecular convergence mechanisms. Thirdly, we summarized the advances of coevolution, including the four main types: mutualism, parasitism, predation and competition. Overall, these latest advances greatly increase our understanding of the underlying molecular mechanisms for diverse adaptive traits and species interaction, demonstrating that the development of evolutionary biology has been greatly accelerated by multi-omics technologies. Finally, we highlighted the emerging trends and future prospects around the above three aspects of adaptive evolution.

Theme and variation in the evolution of insect sex determination

The development of dimorphic adult sexes is a critical process for most animals, one that is subject to intense selection. Work in vertebrate and insect model species has revealed that sex determination mechanisms vary widely among animal groups. However, this variation is not uniform, with a limited number of conserved factors. Therefore, sex determination offers an excellent context to consider themes and variations in gene network evolution. Here we review the literature describing sex determination in diverse insects. We have screened public genomic sequence databases for orthologs and duplicates of 25 genes involved in insect sex determination, identifying patterns of presence and absence. These genes and a 3.5 reference set of 43 others were used to infer phylogenies and compared to accepted organismal relationships to examine patterns of congruence and divergence. The function of candidate genes for roles in sex determination (virilizer, female-lethal-2-d, transformer-2) and sex chromosome dosage compensation (male specific lethal-1, msl-2, msl-3) were tested using RNA interference in the milkweed bug, Oncopeltus fasciatus. None of these candidate genes exhibited conserved roles in these processes. Amidst this variation we wish to highlight the following themes for the evolution of sex determination: (1) Unique features within taxa influence network evolution. (2) Their position in the network influences a component's evolution. Our analyses also suggest an inverse association of protein sequence conservation with functional conservation.

Melanism and color saturation of butterfly assemblages: A comparison between a tropical rainforest and a xeric white forest

The study of butterfly coloration has helped to identify the ecological pressures involved in the evolution of animal coloration. However, almost all studies that addressed this issue have focused on species that inhabit more temperate environments, leaving the species and ecological factors of tropical regions mostly understudied. Here, our purpose was to evaluate whether butterfly assemblages from two distinct Neotropical biomes (i.e., tropical rainforest and xeric white forest) differ regarding their melanism and/or color saturation. Our hypotheses were that (1) tropical rainforest butterflies should be more melanic and color saturated, and that (2) butterflies from more open/arid tropical environments should be more melanic on their dorsal wing surfaces than on their ventral wings. Therefore, we quantified melanism and color saturation from dorsal and ventral surfaces of 121 different butterfly species. Comparisons show that rainforest butterflies, when contrasted to white forest butterflies, have more melanic dorsal wing surfaces, which might be seen as a form of protection against parasites. Our data also show that rainforest butterflies, but not white forest species, have darker dorsal wing surfaces, when compared to their own ventral surfaces, a trend that was also found for species inhabiting both biomes, which might be associated to thermoregulatory advantages. At last, our results also point that butterflies' dorsal wing sides present a higher variance between species (regardless of Biome), when compared to their own ventral wing side, an indication that some ecological factor (e.g., predation avoidance) might be exerting a strong homogenizing force on ventral wing coloration.

First record of the servaline morph in a serval (Leptailurus serval Schreber, 1776) in Akagera National Park, Rwanda

The serval (Leptailurus serval) is a medium-sized felid that is widespread throughout Sub-Saharan Africa. We report on the first scientific record of the servaline morph in a serval in Akagera National Park, Rwanda. Though this morph has been documented in the forested ecosystems of West and Central Africa, the servaline morph has not been recorded in a predominantly savannah ecosystem. Further research into the prevalence of this morph outside of forest ecosystems and the adaptive advantage of the servaline coat pattern in servals is recommended. Our findings highlight the importance of conserving savannah ecosystems to preserve morphological diversity in a species.

Involvement of B-aat1 and Cbs in regulating mantle pigmentation in the Pacific oyster (Crassostrea gigas)

BACKGROUND

Shell color formation is an important physiological process in bivalves, the molecular genetic basis has potential application in bivalve aquaculture, but there is still remaining unclear about this issue. The cystine/glutamate transporter (Slc7a11) and cystathionine beta-synthase (Cbs) are integral genes in pheomelanin synthesis pathway, which is vital to skin pigmentation.

METHODS AND RESULTS

Here, the sequences of b (0, +) -type amino acid transporter 1 (B-aat1) and Cbs in Pacific oyster (Crassostrea gigas) (CgB-aat1, CgCbs) were characterized. Phylogenetically, the deduced amino acid sequences of CgB-aat1 and CgCbs both possessed conserved features. Genes were both ubiquitously expressed in six tested tissues with more abundant expression level in central mantle. Besides, the polyclonal antibodies of CgB-aat1, CgCbs, CgTyr, and CgTyrp2 were successfully prepared. Immunofluorescence analysis revealed that CgB-aat1 and CgCbs proteins were both expressed in gill rudiments of eyed-larvae and concentrated mainly in cytoplasm of epithelial cell and nerve axons in mantle. Additionally, after CgB-aat1 or CgCbs silencing, expressions at mRNA and protein levels of CgB-aat1 and CgCbs involved in pheomelanin synthesis were significantly suppressed, and CgTyr, CgTyrp1 and CgTyrp2 related to eumelanin synthesis were also down-regulated but no apparent differences, respectively. Moreover, micrographic examination found less brown-granules at mantle edge in CgB-aat1 interference group.

CONCLUSION

These results implied that pheomelanin synthesis was possible induced by CgB-aat1-CgTyr-CgCbs axis, and it played an essential role on mantle pigmentation in the oysters. These findings provide the useful genetic knowledge and enrich the physiological information for the shell color formation in bivalve aquaculture.