Molecular and pharmacological characterization of dominant black coat color in sheep

Whole genome resequencing reveals genomic regions related to red plumage in ducks

Plumage color is a characteristic trait of ducks that originates as a result of natural and artificial selection. As a conspicuous phenotypic feature, it is a breed characteristic. Previous studies have identified some genes associated with the formation of black and white plumage in ducks. However, studies on the genetic basis underlying the red plumage phenotype in ducks are limited. Here, genome-wide association analysis (GWAS) and selection signal detection (Fst, θπ ratio, and cross-population composite likelihood ratio [XP-CLR]) were conducted to identify candidate regions and genes underlying duck plumage color phenotype. Selection signal detection revealed 29 overlapping genes (including ENPP1 and ULK1) significantly associated with red plumage color in Ji'an Red ducks. ENSAPLG00000012679, ESRRG, and SPATA5 were identified as candidate genes associated with red plumage using GWAS. Selection signal detection revealed that 19 overlapping genes (including GMDS, PDIA6, and ODC1) significantly correlated with light brown plumage in Brown Tsaiya ducks. GWAS to narrow down the significant regions further revealed nine candidate genes (AKT1, ATP6V1C2, GMDS, LRP4, MAML3, PDIA6, PLD5, TMEM63B, and TSPAN8). Notably, in Brown Tsaiya ducks, GMDS, ODC1, and PDIA6 exhibit significantly differentiated allele frequencies among other feather-colored ducks, while in Ji'an Red ducks, ENSAPLG00000012679 has different allele frequency distributions compared with that in other feather-colored ducks. This study offers new insights into the variation and selection of the red plumage phenotype using GWAS and selective signals.

MC1R Gene Variants and Their Relationship with Coat Color in South American Camelids

In domestic camelids, fleece color is an essential characteristic because it defines the direction of production. Variants were determined in the MC1R gene that showed a relationship with coat color in alpacas and llamas at the level of the coding region. This report sequenced the MC1R gene from 290 alpacas (142 white, 84 black, 50 brown, and 14 light fawn), five brown llamas, nine vicuñas, and three guanacos to analyze the association between coat color and the MC1R gene among South American camelids. A total of nineteen polymorphisms were identified. Seven polymorphisms were significant; three of them were of nonsynonymous type (c.82A > G, c.376G > A, and c.901C > T), two were of synonymous type (c.126 T > C and c.933G > A), one was in the promoter region (-42C > G), and one was in the 3' UTR (+5T > C). More polymorphisms were found in domestic camelids than in wild camelids. Besides polymorphism, the association of polymorphisms might cause white and dark pigmentation in the fleece of South American camelids. In addition, the MC1R protein would answer the pigmentation in alpacas.

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.

Selective sweeps on different pigmentation genes mediate convergent evolution of island melanism in two incipient bird species

Insular organisms often evolve predictable phenotypes, like flightlessness, extreme body sizes, or increased melanin deposition. The evolutionary forces and molecular targets mediating these patterns remain mostly unknown. Here we study the Chestnut-bellied Monarch (Monarcha castaneiventris) from the Solomon Islands, a complex of closely related subspecies in the early stages of speciation. On the large island of Makira M. c. megarhynchus has a chestnut belly, whereas on the small satellite islands of Ugi, and Santa Ana and Santa Catalina (SA/SC) M. c. ugiensis is entirely iridescent blue-black (i.e., melanic). Melanism has likely evolved twice, as the Ugi and SA/SC populations were established independently. To investigate the genetic basis of melanism on each island we generated whole genome sequence data from all three populations. Non-synonymous mutations at the MC1R pigmentation gene are associated with melanism on SA/SC, while ASIP, an antagonistic ligand of MC1R, is associated with melanism on Ugi. Both genes show evidence of selective sweeps in traditional summary statistics and statistics derived from the ancestral recombination graph (ARG). Using the ARG in combination with machine learning, we inferred selection strength, timing of onset and allele frequency trajectories. MC1R shows evidence of a recent, strong, soft selective sweep. The region including ASIP shows more complex signatures; however, we find evidence for sweeps in mutations near ASIP, which are comparatively older than those on MC1R and have been under relatively strong selection. Overall, our study shows convergent melanism results from selective sweeps at independent molecular targets, evolving in taxa where coloration likely mediates reproductive isolation with the neighboring chestnut-bellied subspecies.

Chestnut-bellied Monarchs (Monarcha castaneiventris ugiensis) from two archipelagos in the Solomon Islands have evolved entirely black plumage from a chestnut ancestor (Monarcha castaneiventris megarhynchus), a phenomenon known as island melanism. We obtain and analyze whole genome sequences using traditional summary statistics and new methods that combine inference of the ancestral recombination graph with machine learning. We find multiple lines of evidence for independent selective sweeps on the MC1R and ASIP genes, a receptor/ligand pair which regulates the production of melanin. Melanism on each archipelago is mediated by mutations in one of these two genes. Mutations in and around MC1R underwent a recent soft sweep experiencing strong selection on the islands of Santa Ana and Santa Catalina, whereas selection was also strong but comparatively older for ASIP on the island of Ugi. We show how melanism originated under positive selection on independent molecular targets, evolving convergently in taxa where coloration mediates reproductive isolation.

Genetics of the phenotypic evolution in sheep: a molecular look at diversity-driving genes

BACKGROUND

After domestication, the evolution of phenotypically-varied sheep breeds has generated rich biodiversity. This wide phenotypic variation arises as a result of hidden genomic changes that range from a single nucleotide to several thousands of nucleotides. Thus, it is of interest and significance to reveal and understand the genomic changes underlying the phenotypic variation of sheep breeds in order to drive selection towards economically important traits.

REVIEW

Various traits contribute to the emergence of variation in sheep phenotypic characteristics, including coat color, horns, tail, wool, ears, udder, vertebrae, among others. The genes that determine most of these phenotypic traits have been investigated, which has generated knowledge regarding the genetic determinism of several agriculturally-relevant traits in sheep. In this review, we discuss the genomic knowledge that has emerged in the past few decades regarding the phenotypic traits in sheep, and our ultimate aim is to encourage its practical application in sheep breeding. In addition, in order to expand the current understanding of the sheep genome, we shed light on research gaps that require further investigation.

CONCLUSIONS

Although significant research efforts have been conducted in the past few decades, several aspects of the sheep genome remain unexplored. For the full utilization of the current knowledge of the sheep genome, a wide practical application is still required in order to boost sheep productive performance and contribute to the generation of improved sheep breeds. The accumulated knowledge on the sheep genome will help advance and strengthen sheep breeding programs to face future challenges in the sector, such as climate change, global human population growth, and the increasing demand for products of animal origin.

Melanocortin-1 receptor mutations and pigmentation: Insights from large animals

The melanocortin-1 receptor (MC1R) is a G protein-coupled receptor expressed in cutaneous and hair follicle melanocytes, and plays a central role in coat color determination in vertebrates. Numerous MC1R variants have been identified in diverse species. Some of these variants have been associated with specific hair and skin color phenotypes in humans as well as coat color in animals. Gain-of-function mutations of the MC1R gene cause dominant or partially dominant black/dark coat color, and loss-of-function mutations of the MC1R gene cause recessive or partially recessive red/yellow/pale coat color phenotypes. These have been well documented in a large number of mammals, including human, dog, cattle, horse, sheep, pig, and fox. Higher similarities between large mammals and humans makes them better models to understand pathogenesis of human diseases caused by MC1R mutations. High identities in MC1Rs and similar variants identified in both humans and large mammals also provide an opportunity for receptor structure and function study. In this review, we aim to summarize the naturally occurring mutations of MC1R in humans and large animals.

Functional divergence of the pigmentation gene melanocortin-1 receptor (MC1R) in six endemic Macaca species on Sulawesi Island

Coat color is often highly variable within and between animal taxa. Among hundreds of pigmentation-related genes, melanocortin-1 receptor (MC1R) plays key roles in regulating the synthesis of the dark eumelanin and the red–yellow pheomelanin. The six species of macaques that inhabit Sulawesi Island diverged rapidly from their common ancestor, M. nemestrina. Unlike most macaques, Sulawesi macaques commonly have a dark coat color, with divergence in shade and color pattern. To clarify the genetic and evolutionary basis for coat color in Sulawesi macaques, we investigated the MC1R sequences and functional properties, including basal cAMP production and α-MSH-induced activity in vitro. We found fixed non-synonymous substitutions in MC1R in each species. Furthermore, we found that six species-specific variants corresponded with variation in agonist-induced and basal activity of MC1R. Inconsistent with the dark coat color, four substitutions independently caused decreases in the basal activity of MC1R in M. hecki, M. nigra, M. tonkeana, and M. ochreata. Selective analysis suggested MC1R of M. nigra and M. nigrescens underwent purifying selection. Overall, our results suggest that fixed differences in MC1R resulted in different functional characteristics and might contribute to divergence in color among the six Sulawesi macaque species.

Identification of Candidate Genes for Pigmentation in Camels Using Genotyping-by-Sequencing

The coat color of dromedary is usually uniform and varies from black to white, although dark- to light-brown colors are the most common phenotypes. This project was designed to gain knowledge on novel color-related variants using genotyping-by-sequencing (GBS). The association between the SNPs and coat color was tested using MLM (mixed linear models) with kinship matrix. Three GWAS models including white color vs. non-white color, black vs. non-black color, and light-brown vs. dark-brown color were performed. There were no distinct genetic clusters detected based on the color phenotypes. However, admixture occurred among all individuals of the four different coat color groups. We identified nine significant SNPs associated with white color after Bonferroni correction, located close to ANKRD26, GNB1, TSPYL4, TEKT5, DEXI, CIITA, TVP23B, CLEC16A, TMPRSS13, FXYD6, MPZL3, ANKRD26, HFM1, CDC7, TGFBR3, and HACE1 genes in neighboring flanking regions. The 13 significant SNPs associated with black color and the candidate genes were: CAPN7, CHRM4, CIITA, CLEC16A, COL4A4, COL6A6, CREB3L1, DEXI, DGKZ, DGKZ, EAF1, HDLBP, INPP5F, MCMBP, MDK, SEC23IP, SNAI1, TBX15, TEKT5, TMEM189, trpS, TSPYL4, TVP23B, and UBE2V1. The SNAI1 gene interacted with MCIR, ASIP and KIT genes. These genes play a key role in the melanin biosynthetic and pigmentation biological process and melanogenesis biological pathway. Further research using a larger sample size and pedigree data will allow confirmation of associated SNPs and the identified candidate genes.

Functional characterization of cAMP signaling of variant porcine MC1R alleles in PK15 cells

The melanocortin 1 receptor (MC1R), encoded by the classical extension (E) coat color locus, is expressed on the surface of melanocytes and plays a critical role in switching melanin synthesis from pheomelanin (red/yellow) to eumelanin (black/brown). Different MC1R alleles associated with various coat color patterns in pigs have been identified over the past decades. However, functional analysis of variant porcine MC1R alleles has not yet been performed. Therefore, in this study, we examined the subcellular localization and cyclic adenosine monophosphate (cAMP) signaling capability of MC1R variants in porcine kidney epithelial cells (PK15) overexpressing different MC1R alleles. Transcriptional slippage may partially restore the reading frame of the E^P allele, possibly accounting for the observed spot phenotype. The A243T substitution in the e allele severely disrupted the membrane localization of the MC1R receptor, resulting in a severely impaired cAMP signaling capability. Both the V95M and L102P substitutions in the E^D1 allele may contribute to the constitutively active function of MC1R, thus accounting for the dominant black phenotype. The D124N substitution in the E^D2 allele severely attenuated the cAMP signaling capability of MC1R; however, whether this mutation contributes to the distinct phenotype of Hampshire pigs requires further investigation. Thus, our results provide new insights into the functional characteristics of MC1R variants and their roles in porcine coat color formation.

Origin, demographics, inbreeding, phylogenetics, and phenogenetics of Karamaniko breed, a major common ancestor of the autochthonous Greek sheep

Greece has a long history in autochthonous sheep, the genetic ancestry of which has been associated with four subtypes known to inhabit Greece at the end of the nineteenth century. Among them, the Karamaniko breed is still surviving, however endangered. This study was designed in order to (a) determine the phylogenetic status, (b) to evaluate the levels of inbreeding, and (c) to assess the genetic basis of coat color of Karamaniko breed. For these purposes, the mitochondrial cyt b gene was sequenced, the AFLP methodology was applied, and the MC1R was genotyped, respectively, in 72 female sheep from the Karamaniko breed. Four different novel cyt b haplotypes were defined and three MC1R genotypes were scored, whereas inbreeding levels estimated using AFLPs by the means of relatedness coefficient (r) were 0.287, with gene diversity at the levels of 0.105. Phylogenetic analysis indicated an eastern Asian tropical and subtropical origin of the Karamaniko breed, close with breeds originating from central Turkey, or a clustering within western European or Mediterranean sheep, mirroring a recent genetic divergence, with a non-random spread towards the formation of lowland breeds. The MC1R genotypes were all associated with the white coat color, in which selective breeding has probably been based on traditional morphological characters. Finally, levels of inbreeding do not constitute an indication for a particular mating plan to prevent unpleasant phenomena such as inbreeding depression, probably because of the special attention paid by the farmers towards the avoidance of relative recurrent mating.

Trends towards revealing the genetic architecture of sheep tail patterning: Promising genes and investigatory pathways

Different sheep breeds have evolved after initial domestication, generating various tail phenotypic patterns. The phenotypic diversity of sheep tail patterns offers ideal materials for comparative analysis of its genetic basis. Evolutionary biologists, animal geneticists, breeders, and producers have been curious to clearly understand the underlying genetics behind phenotypic differences in sheep tails. Understanding the causal gene(s) and mutation(s) underlying these differences will help probe an evolutionary riddle, improve animal production performance, promote animal welfare, and provide lessons that help comprehend human diseases related to fat deposition (i.e., obesity). Historically, fat tails have served as an adaptive response to aridification and climate change. However, the fat tail is currently associated with compromised mating and animal locomotion, fat distribution in the animal body, increased raising costs, reduced consumer preference, and other animal welfare issues such as tail docking. The developing genomic approaches provide unprecedented opportunities to determine causal variants underlying phenotypic differences among populations. In the last decade, researchers have performed several genomic investigations to assess the genomic causality underlying phenotypic variations in sheep tails. Various genes have been suggested with the prominence of several potentially significant causatives, including the BMP2 and PDGFD genes associated with the fat tail phenotype and the TBXT gene linked with the caudal vertebrae number and tail length. Although the potential genes related to sheep tail characteristics have been revealed, the causal variant(s) and mutation(s) of these high-ranking candidate genes are still elusive and need further investigation. The review discusses the potential genes, sheds light on a knowledge gap, and provides possible investigative approaches that could help determine the specific genomic causatives of sheep tail patterns. Besides, characterizing and revealing the genetic determinism of sheep tails will help solve issues compromising sheep breeding and welfare in the future.

Detecting the footprint of selection on the genomes of Murciano-Granadina goats

Artificial selection is one of the major forces modifying the genetic composition of livestock populations. Identifying genes under selection could be useful to elucidate their impact on phenotypic variation. We aimed to identify genomic regions targeted by selection for dairy and pigmentation traits in Murciano-Granadina goats. Performance of a selection scan based on the integrated haplotype score test in a population of 1183 Murciano-Granadina goats resulted in the identification of 77 candidate genomic regions/SNPs. The most significant selective sweeps mapped to chromosomes 1 (69.86 Mb), 4 (41.80-49.95 Mb), 11 (65.74 Mb), 12 (31.24 and 52.51 Mb), 17 (34.76-37.67 Mb), 22 (31.75 Mb), and 26 (26.69-31.05 Mb). By using previously generated RNA-Seq data, we built a catalogue of 6414 genes that are differentially expressed across goat lactation (i.e. 78 days post-partum, early lactation; 216 days post-partum, late lactation; 285 days post-partum, dry period). Interestingly, 183 of these genes mapped to selective sweeps and several of them display functions related with lipid, protein, and carbohydrate metabolism, insulin signaling, cell proliferation, as well as mammary development and involution. Of particular interest are the CSN3 and CSN1S2 genes, which encode two major milk proteins. Additionally, we found three pigmentation genes (GLI3, MC1R, and MITF) co-localizing with selective sweeps. Performance of a genome-wide association study and Sanger sequencing and TaqMan genotyping experiments revealed that the c.801C>G (p.Cys267Trp) polymorphism in the melanocortin 1 receptor (MC1R) gene is the main determinant of the black (GG or GC genotypes) and brown (CC genotypes) colorations of Murciano-Granadina goats.

Museomics dissects the genetic basis for adaptive seasonal colouration in the least weasel

Dissecting the link between genetic variation and adaptive phenotypes provides outstanding opportunities to understand fundamental evolutionary processes. Here, we use a museomics approach to investigate the genetic basis and evolution of winter coat colouration morphs in least weasels (Mustela nivalis), a repeated adaptation for camouflage in mammals with seasonal pelage colour moults across regions with varying winter snow. Whole-genome sequence data was obtained from biological collections and mapped onto a newly assembled reference genome for the species. Sampling represented two replicate transition zones between nivalis and vulgaris colouration morphs in Europe, which typically develop white or brown winter coats, respectively. Population analyses showed that the morph distribution across transition zones is not a by-product of historical structure. Association scans linked a 200 kb genomic region to colouration morph, which was validated by genotyping museum specimens from inter-morph experimental crosses. Genotyping the wild populations narrowed down the association to pigmentation gene MC1R and pinpointed a candidate amino acid change co-segregating with colouration morph. This polymorphism replaces an ancestral leucine residue by lysine at the start of the first extracellular loop of the protein in the vulgaris morph. A selective sweep signature overlapped the association region in vulgaris, suggesting that past adaptation favoured winter-brown morphs and can anchor future adaptive responses to decreasing winter snow. Using biological collections as valuable resources to study natural adaptations, our study showed a new evolutionary route generating winter colour variation in mammals and that seasonal camouflage can be modulated by changes at single key genes.

Specific mutations in the genes of MC1R and TYR have an important influence on the determination of pheomelanin pigmentation in Korean native chickens

Objective:

The TYR (Tyrosinase) and MC1R (Melanocortin 1 receptor) genes are recognized as important genes involved in plumage pigmentation in Korean native chickens. Specifically, most color patterns in chicken result from differential expression of the TYR gene. In this study, the co-segregation of the pigmentation and sequence of the TYR and MC1R genes was investigated through intercrosses between red (R1q1), red with black and black plumage color types of native Korean chickens.

Materials and Methods:

Using DNA, RNA, and tissue by plumage color of each Korean native chickens, the role of major genes in pigmentation of pheomelanin was evaluated. Reverse transcription polymerase chain reaction, sequencing, western blot, and immunohistochemical were performed to determine the effect of TYR and MC1R genes on plumage pigmentation in Korean native chickens.

Results:

The KCO line (Korean chicken Ogol: Black-line) with an EEC _ genotype exhibited black feathers, whereas red and red mixed with black chicken with EeC genotype exhibited white feathers. There were notable differences between the base sequences of MC1R and TYR in three Korean chicken breeds, with the highest variation in TYR. Perhaps this is the key characteristics of Korean chicken. Further, we analyzed the expression patterns of MC1R and TYR genes in each type of Korea native chicken and observed that TYR expression was high in feather follicle (R1q2) of KCO tissue. However, native red (Korean chicken red) and native red with black (Korean chicken red dark) chickens have increased TYR expression in the tissue. However, the expression of MC1R was much different from that of TYR.

Conclusion:

In this study, our results suggest that the differences in position and TYR expression levels exert more influence on plumage pigmentation in native Korean chicken breeds than changes in MC1R expression levels.

Polymorphism of the Melanocortin 1 Receptor (MC1R) Gene and its Role in Determining the Coat Colour of Central European Cattle Breeds

There are many genes responsible for the appearance of different coat colours, among which the melanocortin 1 receptor gene (MC1R) plays an important role. The aim of the study was to characterize genetic variation in Central European cattle breeds based on polymorphism of the MC1R gene and factors determining their coat colour. The study was conducted on 290 individuals of the following breeds: Polish White-Backed (PW), Lithuanian White-Backed (LW), Polish Red (PR), Lithuanian Red (LR), Carpathian Brown (CB), Ukrainian Grey (UG), and Slovak Pinzgau (SP). Polymorphism at the MC1R gene locus was analysed by polymerase chain reaction-restriction fragment length polymorphism (PCR-RFLP) using two restriction enzymes: Cfr10I and SsiI. The proportions of alleles and genotypes in the MC1R locus indicates a strong relationship between polymorphism and the coat colour of cattle: The ED allele proved to be characteristic for the breeds with a white-backed coat (PW and LW), while the dominant allele in the red breeds (PR and LR) was E+. It is noteworthy that coat colour in the SP population was determined only by the recessive e allele, which resulted in the formation of a separate clade in the phylogenetic tree.

Genomic mapping identifies two genetic variants in the MC1R gene for coat colour variation in Chinese Tan sheep

Coat colour is one of the most important economic traits of sheep and is mainly used for breed identification and characterization. This trait is determined by the biochemical function, availability and distribution of phaeomelanin and eumelanin pigments. In our study, we conducted a genome-wide association study to identify candidate genes and genetic variants associated with coat colour in 75 Chinese Tan sheep using the ovine 600K SNP BeadChip. Accordingly, we identified two significant SNPs (rs409651063 at 14.232 Mb and rs408511664 at 14.228 Mb) associated with coat colour in the MC1R gene on chromosome 14 with −log10(P) = 2.47E-14 and 1.00E-13, respectively. The consequence of rs409651063 was a missense variant (g.14231948 G>A) that caused an amino acid change (Asp105Asn); however, the second SNP (rs408511664) was a synonymous substitution and is an upstream variant (g.14228343G>A). Moreover, our PCR analysis revealed that the genotype of white sheep was exclusively homozygous (GG), whereas the genotypes of black-head sheep were mainly heterozygous (GA). Interestingly, allele-specific expression analysis (using the missense variant for the skin cDNA samples from black-head sheep) revealed that only the G allele was expressed in the skin covered with white hair, while both the G and A alleles were expressed in the skin covered with black hair. This finding indicated that the missense mutation that we identified is probably responsible for white coat colour in Tan sheep. Furthermore, qPCR analysis of MC1R mRNA level in the skin samples was significantly higher in black-head than white sheep and very significantly higher in GA than GG individuals. Taken together, these results help to elucidate the genetic mechanism underlying coat colour variation in Chinese indigenous sheep.

A Stop-Gain Mutation within MLPH Is Responsible for the Lilac Dilution Observed in Jacob Sheep

A coat color dilution, called lilac, was observed within the Jacob sheep breed. This dilution results in sheep appearing gray, where black would normally occur. Pedigree analysis suggested an autosomal recessive inheritance. Whole-genome sequencing of a dilute case, a known carrier, and sixteen non-dilute sheep was used to identify the molecular variant responsible for the coat color change. Through investigation of the genes MLPH, MYO5A, and RAB27A, we discovered a nonsynonymous mutation within MLPH, which appeared to match the reported autosomal recessive nature of the lilac dilution. This mutation (NC_019458.2:g.3451931C>A) results in a premature stop codon being introduced early in the protein (NP_001139743.1:p.Glu14*), likely losing its function. Validation testing of additional lilac Jacob sheep and known carriers, unrelated to the original case, showed a complete concordance between the mutation and the dilution. This stop-gain mutation is likely the causative mutation for dilution within Jacob sheep.

Recent advances in amniote palaeocolour reconstruction and a framework for future research

Preserved melanin pigments have been discovered in fossilised integumentary appendages of several amniote lineages (fishes, frogs, snakes, marine reptiles, non-avialan dinosaurs, birds, and mammals) excavated from lagerstätten across the globe. Melanisation is a leading factor in organic integument preservation in these fossils. Melanin in extant vertebrates is typically stored in rod- to sphere-shaped, lysosome-derived, membrane-bound vesicles called melanosomes. Black, dark brown, and grey colours are produced by eumelanin, and reddish-brown colours are produced by phaeomelanin. Specific morphotypes and nanostructural arrangements of melanosomes and their relation to the keratin matrix in integumentary appendages create the so-called 'structural colours'. Reconstruction of colour patterns in ancient animals has opened an exciting new avenue for studying their life, behaviour and ecology. Modern relationships between the shape, arrangement, and size of avian melanosomes, melanin chemistry, and feather colour have been applied to reconstruct the hues and colour patterns of isolated feathers and plumages of the dinosaurs Anchiornis, Sinosauropteryx, and Microraptor in seminal papers that initiated the field of palaeocolour reconstruction. Since then, further research has identified countershading camouflage patterns, and informed subsequent predictions on the ecology and behaviour of these extinct animals. However, palaeocolour reconstruction remains a nascent field, and current approaches have considerable potential for further refinement, standardisation, and expansion. This includes detailed study of non-melanic pigments that might be preserved in fossilised integuments. A common issue among existing palaeocolour studies is the lack of contextualisation of different lines of evidence and the wide variety of techniques currently employed. To that end, this review focused on fossil amniotes: (i) produces an overarching framework that appropriately reconstructs palaeocolour by accounting for the chemical signatures of various pigments, morphology and local arrangement of pigment-bearing vesicles, pigment concentration, macroscopic colour patterns, and taphonomy; (ii) provides background context for the evolution of colour-producing mechanisms; and (iii) encourages future efforts in palaeocolour reconstructions particularly of less-studied groups such as non-dinosaur archosaurs and non-archosaur amniotes.

Discovery of a new nucleotide substitution in the MC1R gene and haplotype distribution in native and non-Japanese chicken breeds

Melanocortin 1-receptor (MC1R) is one of the major genes that controls chicken plumage colour. In this study, we investigated the sequence and haplotype distribution of the MC1R gene in native Japanese chickens, along with non-Japanese chicken breeds. In total, 732 and 155 chickens from 30 Japanese and eight non-Japanese breeds respectively were used. Three synonymous and 11 non-synonymous nucleotide substitutions were detected, resulting in 15 haplotypes (H0-H14). Of these, three were newly found haplotypes (H9, H13 and H14), of which one (H9) was composed of known substitutions C69T, T212C, G274A and G636A. The second one (H13) possessed newly found non-synonymous substitution C919G, apart from the known substitutions C69T, G178A, G274A, G636A and T637C. The third one (H14) comprised a newly discovered substitution C919G in addition to the known C69T, G274A and G409A substitutions. The homozygote for this new haplotype exhibited wt like plumage despite the presence of G274A. In addition to discovering a new nucleotide substitution (C919G) and three new haplotypes, we defined the plumage colour of the bird that was homozygous for the A644C substitution (H5 haplotype) as wheaten-like for the first time; although the substitution has been already reported, its effect was not revealed. Besides detecting the new plumage colour, we also confirmed that the A427G and G274A substitutions contribute in expressing brownish and black plumage colour respectively, as reported by the previous studies. Moreover, we confirmed that the buttercup allele does not express black plumage despite possessing a G274A substitution, under the suppression effect of A644C. In contrast, the birds homozygous for the birchen allele presented solid black plumage, which was contradictory to the previous reports. In conclusion, we revealed a large diversity in the MC1R gene of native Japanese chicken breeds, along with the discovery of a new non-synonymous nucleotide substitution (C919G) and three novel haplotypes (H9, H13 and H14).

Transcriptome profiling analysis reveals key genes of different coat color in sheep skin

Background

To investigate the molecular mechanisms determining the coat color of native breed sheep in Xinjiang.

Methods

Bashibai sheep, Yemule white sheep and Tulufan black sheep were selected. Illumina HiSeq X Ten sequencing technology was used to detect the genes responsible for the white, light brown, black and cyan gray coat colors in sheep. Sequence analysis and functional gene annotation analysis were performed to analyze the results. The signal pathways and differentially expressed genes related to sheep hair color production regulation were screened and finally verified by real-time polymerase chain reaction.

Results

Functional annotation by Kyoto Encyclopedia of Genes and Genomes analysis revealed significant differences in enrichment of immunity-related pathways as well as melanogenesis synthetic and tyrosine metabolism pathways. Our results showed that the DCT, TYR, TYRP1, PMEL, SLC45A2 and MLANA six genes may be associated with the regulation of coat color development and provide a theoretical basis for selecting natural coat colors of sheep.

Skin transcriptome profiles associated with black- and white-coated regions in Boer and Macheng black crossbred goats

To increase the current understanding of the gene-expression profiles in different skin regions associated with different coat colors and identify key genes for the regulation of color patterns in goats, we used the Illumina RNA-Seq method to compare the skin transcriptomes of the black- and white-coated regions containing hair follicles from the Boer and Macheng Black crossbred goat, which has a black head and a white body. Six cDNA libraries derived from skin samples of the white-coated region (n = 3) and black-coated region (n = 3) were constructed from three full-sib goats. On average, we obtained approximately 76.5 and 73.5 million reads for skin samples from black- and white-coated regions, respectively, of which 75.39% and 76.05% were covered in the genome database. A total of 165 differentially expressed genes (DEGs) were detected between these two color regions, among which 110 were upregulated and 55 were downregulated in the skin samples of white- vs. black-coated regions. The results of Gene Ontology and Kyoto Encyclopedia of Genes and Genomes enrichment analyses revealed that some of these DEGs may play an important role in controlling the pigmentation of skin or hair follicles. We identified three key DEGs, i.e., Agouti, DCT, and TYRP1, in the pathway related to melanogenesis in the different skin regions of the crossbred goat. DCT and TYRP1 were downregulated and Agouti was upregulated in the skin of the white-coated region, suggesting a lack of mature melanocytes in this region and that Agouti might play a key developmental role in color-pattern formation. All data sets (Gene Expression Omnibus) are available via public repositories. In addition, MC1R was genotyped in 200 crossbred goats with a black head and neck. Loss-of-function mutations in MC1R as well as homozygosity for the mutant alleles were widely found in this population. The MC1R gene did not seem to play a major role in determining the black head and neck in our crossbred goats. Our study provides insights into the transcriptional regulation of two distinct coat colors, which might serve as a key resource for understanding coat color pigmentation in goats. The region-specific expression of Agouti may be associated with the distribution of pigments across the body in Boer and Macheng Black crossbred goats.

Identification of two TYRP1 loss-of-function alleles in Valais Red sheep

The Valais Red sheep breed is a local breed of the Swiss canton Valais. Although the breed is characterised by its brown colour, black animals occasionally occur and the objective of this study was to identify the causative genetic variants responsible for the obvious difference. A GWAS using high-density SNP data to compare 51 brown and 38 black sheep showed a strong signal on chromosome 2 at the TYRP1 locus. Haplotype analyses revealed three different brown-associated alleles. The WGS of three sheep revealed four protein-changing variants within the TYRP1 gene. Three of these variants were associated with the recessively inherited brown coat colour. This includes the known missense variant TYRP1:c.869G>T designated as b^{S oay} and two novel loss-of-function variants. We propose to designate the frame-shift variant TYRP1:c.86_87delGA as b^{VS 1} and the nonsense variant TYRP1:c.1066C>T as b^{VS 2} . Interestingly, the b^{VS 1} allele occurs only in local breeds of Switzerland whereas the b^{VS 2} allele seems to be more widespread across Europe.

A Novel Pale-Yellow Coat Color of Rabbits Generated via MC1R Mutation With CRISPR/Cas9 System

Coat color is of great importance in animal breed characteristics; it is not only a significant productive trait but also an indispensable economic trait, especially in the rabbit industry. In the present study, the relationship between melanocortin 1 receptor (MC1R) genotypes and coat color phenotypes was observed in five rabbit breeds with popular coat colors that are present in China. These breeds comprised the Lianshan black rabbit (BR), Fujian yellow rabbit (YR), New Zealand white rabbit (WR), Gray Giant rabbit (GR), and Checkered Giant rabbit (CR), which were firstly determined, and the results showed that GR had an E allele; WR, CR, and BR had a 6-bp in-frame deletion (c.281_286del6, E^D allele); and YR had a 30-bp deletion (c.304_333del30 E allele). To explore the feasibility of obtaining a novel rabbit coat color through the mutation of MC1R with the CRISPR/Cas9 system, two single-guide RNAs (sgRNAs) were designed for the MC1R gene, and the editing efficiency was confirmed by injection of rabbits' zygotes. Unlike the donor rabbits whose coat color was originally black, two novel pale-yellow-coated rabbits were generated in the founders. A total of six novel MC1R gene deletions were identified in the two founder rabbits, in which the longest deletion was more than 700 bp. The histological hematoxylin-and-eosin (H&E) staining results indicated that eumelanin amounts were absent in hair follicles of MC1R-knockout (KO) rabbits, when compared with that of donor BR. In addition, the messenger RNA (mRNA) levels of some key downstream genes in the MC1R pathway were all downregulated in MC1R-KO rabbits compared with BR and YR. These results further indicate that loss-of-function MC1R contributed to blocking the synthesis of eumelanin and created a novel pale-yellow coat color in the MC1R-KO rabbits, and gene editing technology may be a useful tool to generate novel phenotypes in rabbit breeding.

Mutations in ASIP and MC1R: dominant black and recessive black alleles segregate in native Swedish sheep populations

By studying genes associated with coat colour, we can understand the role of these genes in pigmentation but also gain insight into selection history. North European short-tailed sheep, including Swedish breeds, have variation in their coat colour, making them good models to expand current knowledge of mutations associated with coat colour in sheep. We studied ASIP and MC1R, two genes with known roles in pigmentation, and their association with black coat colour. We did this by sequencing the coding regions of ASIP in 149 animals and MC1R in 129 animals from seven native Swedish sheep breeds in individuals with black, white or grey fleece. Previously known mutations in ASIP [recessive black allele: g.100_105del (D₅ ) and/or g.5172T>A] were associated with black coat colour in Klövsjö and Roslag sheep breeds and mutations in both ASIP and MC1R (dominant black allele: c.218T>A and/or c.361G>A) were associated with black coat colour in Swedish Finewool. In Gotland, Gute, Värmland and Helsinge sheep breeds, coat colour inheritance was more complex: only 11 of 16 individuals with black fleece had genotypes that could explain their black colour. These breeds have grey individuals in their populations, and grey is believed to be a result of mutations and allelic copy number variation within the ASIP duplication, which could be a possible explanation for the lack of a clear inheritance pattern in these breeds. Finally, we found a novel missense mutation in MC1R (c.452G>A) in Gotland, Gute and Värmland sheep and evidence of a duplication of MC1R in Gotland sheep.

Combined approaches to identify genomic regions involved in phenotypic differentiation between low divergent breeds: Application in Sardinian sheep populations

Selective breeding has led to modifications in the genome of many livestock breeds. In this study, we identified the genomic regions that may explain some of the phenotypic differences between two closely related breeds from Sardinia. A total of 44 animals, 20 Sardinian Ancestral Black (SAB) and 24 Sardinian White (SW), were genotyped using the Illumina Ovine 50K array. A total of 68, 38 and 15 significant markers were identified using the case-control genome-wide association study (GWAS), the Bayesian population differentiation analysis (F_ST ) and the Rsb metric, respectively. Comparisons among the approaches revealed a total of 22 overlapping markers between GWAS and F_ST and one marker between GWAS and Rsb. Three markers detected by Rsb were also located near (<2 Mb) to highly significant regions identified by GWAS and F_ST analyses. Moreover, one candidate marker identified by GWAS and F_ST approaches was located in a run of homozygosity island that was shared by both breeds. We identified several genes involved in many phenotypic differences (such as stature and growth, reproduction, ear size, coat colour, behaviour) between the two analysed breeds. This study shows that combining several genome-wide approaches could improve discovery of regions involved in the variability of breeding traits and responsible for the phenotypic diversity even between closely related breeds. Overall, the combination of such genome-wide methods can be extended to other livestock breeds that share between them a similar genetic background, to understand the process that shapes the patterns of genetic variability between closely related populations.

Exploring differentially expressed genes associated with coat color in goat skin using RNA-seq

Fur color in domestic goats is an important, genetically determined characteristic that is associated with economic value. This study was designed to perform a comprehensive expression profiling of genes expressed in the skin tissues from Laiwu Black goat and Lubei White goat. Comparisons of black and white goat skin transcriptomes revealed 102 differentially expressed genes (DEGs), of which 38 were upregulated and 64 downregulated in black skin compared with white skin. Among the DEGs, we identified six genes involved in pigmentation, including agouti signaling protein (ASIP), CAMP responsive element binding protein 3-like 1 (CREB3L1), dopachrome tautomerase (DCT), premelanosome protein (PMEL), transient receptor potential cation channel subfamily M member 1 (TRPM1), and tyrosinase-related protein 1 (TYRP1). Notably, there were no significant differences in the expression of melanocortin 1 receptor, microphthalmia-associated transcription factor, tyrosinase, and KIT proto-oncogene receptor tyrosine kinase between the black and white skin samples, whereas ASIP expression was detected only in white skin. PMEL, TRPM1, TYRP1, and DCT showed higher expression in black goat skin, but ASIP and CREB3L1 had higher expression in white goat skin. Quantitative polymerase chain reaction results for PMEL, TRPM1, DCT, TYRP1, and CREB3L1 expression were consistent with those for RNA-seq. These results will expand our understanding of the complex molecular mechanisms of skin physiology and melanogenesis in goats, and provide a foundation for future studies.

Genomic signatures of high-altitude adaptation in Ethiopian sheep populations

BACKGROUND

Ethiopian sheep populations such as Arsi-Bale, Horro and Adilo (long fat-tailed, LFT) inhabit mid to high-altitude areas; and Menz sheep (MZ, short fat-tailed) are adapted to cool sub-alpine environments. In contrast, Blackhead Somali sheep (BHS, fat-rumped) thrive well in arid and semi-arid areas characterized by high temperature and low precipitation. The genomic investigation of Ethiopian sheep populations may help to identify genes and biological pathways enable to adapt to the different ecological conditions.

OBJECTIVE

To uncover genomic regions and genes showing evidence of positive selection for altitude adaptation in Ethiopian sheep populations.

METHODS

A total of 72 animals inhabiting high-versus low-altitude environments were genotyped on an Ovine Infinium HD array (~ 600 K). Pairwise genetic differentiation (Fst) was calculated in sliding windows of 20 SNPs and the upper 1% smoothed Fst values were considered to represent positive selection signatures. Genes within < 25 kb of the most differentiated SNPs were considered as selection candidates.

RESULTS

Signatures of selection were detected in genes known to be associated high with altitude adaptation in MZ-BHS pair comparison (PPP1R12A, RELN, PARP2, and DNAH9) and in LFT-BHS pair comparison (VAV3, MSRB3,EIF2AK4, MET, and TACR1). The candidate genes (MITF, FGF5, MTOR, TRHDE, and TUBB3) associated with altitude adaptation and shared between the MZ-BHS and LTF-BHS pair comparisons were also detected as under selection. Further functional analyses reveal that the candidate genes were involved in biological processes and pathways relevant to adaptation under extreme altitudes, including respiratory system development and smoothened signaling pathway.

CONCLUSION

The results of the present study could aid in-depth understanding and exploitation of the underlying genetic mechanisms for sheep and other livestock species adaptation to high-altitude environments.

Polymorphism in MC1R, TYR and ASIP genes in different colored feather chickens

Coat color genetics successfully adapted and applied to different animal species, which provides a good demonstration of the concept of comparative genetics. In this study, we sequenced 945 bp fragments of melanocortin 1 receptor (MC1R) gene, 421 bp fragments of exon 1 of tyrosinase (TYR) gene and 266 bp fragments of exon 3 of agouti signaling protein (ASIP) gene for 250 individuals with five plumage color patterns. We detected a total of three SNPs (T398A, T637C, and G920C) in MC1R and built six haplotypes (H1-H6) based on the three SNPs. H5 and H6 haplotypes were mainly concentrated in white and grey chicken. And diplotypes H2H3 occurred in white feather and black-speckle feather with the same frequency. Moreover, a total of three SNPs (C47G, T120C, and T172C) in TYR were found and built six haplotypes (P1-P6) based on the three SNPs. Among them, haplotype P2, P3 and P6 were not occurred in black chicken, the diplotypes P1P6 and P4P6 were only distributed in white, gray and black-speckled feather. We only detected one SNP (T168C) in ASIP gene and found that genotype TT was advantage genotype in the different plumage color groups of chickens. Collectively, our study suggested an association between plumage color and genetic variation of MC1R, TYR and ASIP in chicken.

Identification and sequence characterization of melanocortin 1 receptor gene (MC1R) in Bos indicus versus (Bos taurus X Bos indicus)

Melanocortin 1 receptor (MC1R) plays a vital role in melanogenesis and determines coat color of mammals. Polymorphic variants in MC1R, causing coat color variation, were described in few mammals; however, such studies were not done in cattle. The objective of the study was to explore the association of MC1R gene polymorphism within Tharparkar (Bos indicus) and Karan Fries (B. indicus X Bos taurus) cattle. Genomic DNA isolated from blood samples of Tharparkar breed by modified Phenol: Chloroform; Isoamyl alcohol method. Using genomic DNA as template for PCR, MC1R gene was amplified and sequenced. The sequences were analyzed and submitted to Genbank with Acc.No MG373615-MG373644. Comparison of sequence alignment with other bovine species using ClustalW revealed 99-96% similarity. MC1R gene phylogenetic analyses were analyzed using MEGA X. The MC1R gene tree, protein domains and genetic variation of cattle were retrieved from Ensemble Asia Cattle Genome Browser. Eight single nucleotide polymorphisms (SNPs) (c.296T > C, c.583T > C, c.663C > T, c.830T > C, c.853G > A, c.880G > A, c.906C > G, c.927C > T) in CDS reveal high genetic variability. Subsequent to amino acid changes p.L99P, p.F195L, p.F277S, p.A285T and p.D293N, p.R302S, respectively found in seven-transmembrane. Mutations appeared in MC1R of B. taurus with white and black coat color as compared to B. indicus with white coat.

Association of a novel SNP in the ASIP gene with skin color in black-bone chicken

The agouti signaling protein gene (ASIP) is a widely studied pigmentation gene that plays an important role in melanin synthesis. To determine the variety of ASIP expression in the Muchuan Black-Bone chicken, we examined genetic variation in the ASIP promoter region. A single nucleotide polymorphism (c.-1826A>T) was found to be associated with the skin color (dorsal and subalar) of black-bone chicken. Individuals with TT and AT genotypes had higher ASIP mRNA levels in the skin than did those with the AA genotype (P < 0.01). In addition, individuals with the TT genotype had higher ASIP mRNA levels than did those with the AT genotype (P < 0.05). Expression of melanogenesis-related genes (melanocortin 1 receptor and tyrosinase genes) was higher in the skin of chickens with the TT and AT genotypes than in those with the AA genotype (P < 0.01). A luciferase assay showed that promoter activity was higher in chickens with the TT genotype than in those with the AA genotype. Putative transcription factor prediction suggested that the c.-1826A>T mutation might shift the promoter binding affinity with differential transcription factors. In summary, we identified a novel mutation in the ASIP gene promoter that may affect chicken skin color by altering ASIP transcriptional activity.

Sequence characterization and polymorphism of melanocortin 1 receptor gene in some goat breeds with different coat color of Mongolia

Objective

Extension and Agouti loci play a key role for proportions of eumelanin and pheomelanin in determining coat color in several species, including goat. Mongolian goats exhibit diverse types of coat color phenotypes. In this study, investigation of the melanocortin 1 receptor (MC1R) coding region in different coat colors in Mongolian goats was performed to ascertain the presence of the extension allele.

Methods

A total of 105 goat samples representing three goat breeds were collected for this study from middle Mongolia. A 938 base pair (bp) long coding region of the MC1R gene was sequenced for three different breeds with different coat colors (Gobi Gurwan Saikhan: complete black, Zalaa Jinstiin Tsagaan: complete white, Mongolian native goat: admixture of different of coat colors). The genotypes of these goats were obtained from analyzing and comparing the sequencing results.

Results

A total of seven haplotypes defined by five substitution were identified. The five single nucleotide polymorphisms included two synonymous mutations (c.183C>T and c.489G>A) and three missense (non-synonymous) mutations (c.676A>G, c.748T>G, and c.770T>A). Comparison of genotypes frequencies of two common missense mutions using chi-sqaure (x²) test revealed significant differences between coat color groups (p<0.001). A logistic regression analysis additionally suggested highly significant association between genotypes and variation of black versus white uniform combination. Alternatively, most investigated goats (60.4%) belonged to H2 (TGAGT) haplotype.

Conclusion

According to the findings obtained in this study on the investigated coat colors, mutations in MC1R gene may have the crucial role for determining eumelanin and pheomelanin phenotypes. Due to the complication of coat color phenotype, more detailed investigation needed.

Genetic Variation in Coat Colour Genes MC1R and ASIP Provides Insights Into Domestication and Management of South American Camelids

The domestication of wild vicuña and guanaco by early pre-Inca cultures is an iconic example of wildlife management and domestication in the Americas. Although domestic llamas and alpacas were clearly selected for key, yet distinct, phenotypic traits, the relative patterns and direction of selection and domestication have not been confirmed using genetic approaches. However, the detailed archaeological records from the region suggest that domestication was a process carried out under significant control and planning, which would have facilitated coordinated and thus extremely effective selective pressure to achieve and maintain desired phenotypic traits. Here we link patterns of sequence variation in two well-characterised genes coding for colour variation in vertebrates and interpret the results in the context of domestication in guanacos and vicuñas. We hypothesise that colour variation in wild populations of guanacos and vicunas were strongly selected against. In contrast, variation in coat colour variation in alpaca was strongly selected for and became rapidly fixed in alpacas. In contrast, coat colour variants in llamas were of less economic value, and thus were under less selective pressure. We report for the first time the full sequence of MC1R and 3 exons of ASIP in 171 wild specimens from throughout their distribution and which represented a range of commonly observed colour patterns. We found a significant difference in the number of non-synonymous substitutions, but not synonymous substitutions among wild and domestics species. The genetic variation in MC1R and ASIP did not differentiate alpaca from llama due to the high degree of reciprocal introgression, but the combination of 11 substitutions are sufficient to distinguish domestic from wild animals. Although there is gene flow among domestic and wild species, most of the non-synonymous variation in MC1R and ASIP was not observed in wild species, presumably because these substitutions and the associated colour phenotypes are not effectively transmitted back into wild populations. Therefore, this set of substitutions unequivocally differentiates wild from domestic animals, which will have important practical application in forensic cases involving the poaching of wild vicuñas and guanacos. These markers will also assist in identifying and studying archaeological remains pre- and post-domestication.

Transcriptomics reveals the molecular processes of light-induced rapid darkening of the non-obligate cave dweller Oreolalax rhodostigmatus (Megophryidae, Anura) and their genetic basis of pigmentation strategy

BACKGROUND

Vertebrates use different pigmentation strategies to adapt to various environments. A large amount of research has been done on disclosing the mechanisms of pigmentation strategies in vertebrates either under light, or, living in constant darkness. However, less attention has been paid to non-obligate, darkness dwellers. Red-spotted toothed toads Oreolalax rhodostigmatus (Megophryidae; Anura) from the karst mountainous region of southwestern China are non-obligate cave dwellers. Most tadpoles of the species possess transparent skin as they inhabit the dark karst caves. But remarkably, the transparent tadpoles can darken just within 15 h once exposed to light. Obviously, it is very significant to reveal molecular mechanisms of the unexpected rapid-darkening phenomenon.

RESULTS

We compared the transcriptomes of O. rhodostigmatus tadpoles with different durations of light exposure to investigate the cellular processes and potential regulation signals for their light-induced rapid darkening. Genes involved in melanogenesis (i.e. TYR, TYRP1 and DCT) and melanocyte proliferation, as well as their transcriptional factor (MITF), showed light-induced transcription, suggesting a dominating role of morphological color change (MCC) in this process. Transcription of genes related to growth factor, MAPK and PI3K-Akt pathways increased with time of light exposure, suggesting that light could induce significant growth signal, which might facilitate the rapid skin darkening. Most importantly, an in-frame deletion of four residues was identified in O. rhodostigmatus melanocortin-1 receptor (MC1R), a critical receptor in MCC. This deletion results in a more negatively charged ligand pocket with three stereo-tandem aspartate residues. Such structural changes likely decrease the constitutive activity of MC1R, but increase its ligands-dependent activity, thus coordinating pigment regression and rapid melanogenesis in the dark and light, respectively.

CONCLUSION

Our study suggested that rapid MCC was responsible for the light-induced rapid darkening of O. rhodostigmatus tadpoles. Genetic mutations of MC1R in them could explain how these non-obligate cave dwellers coordinate pigment regression and robust melanogenesis in darkness and light, respectively. To our knowledge, this is the first study that reports the association between pigmentation phenotype adaptation and MC1R mutations in amphibians and/or in non-obligate cave dwellers.

Genetic background of coat colour in sheep

Abstract. The coat colour of animals is an extremely important trait that affects their behaviour and is decisive for survival in the natural environment. In farm animal breeding, as a result of the selection of a certain coat colour type, animals are characterized by a much greater variety of coat types. This makes them an appropriate model in research in this field. A very important aspect of the coat colour types of farm animals is distinguishing between breeds and varieties based on this trait. Furthermore, for the sheep breeds which are kept for skins and wool, coat/skin colour is an important economic trait. Until now the study of coat colour inheritance in sheep proved the dominance of white colour over pigmented/black coat or skin and of black over brown. Due to the current knowledge of the molecular basis of ovine coat colour inheritance, there is no molecular test to distinguish coat colour types in sheep although some are available for other species, such as cattle, dogs, and horses. Understanding the genetic background of variation in one of the most important phenotypic traits in livestock would help to identify new genes which have a great effect on the coat colour type. Considering that coat colour variation is a crucial trait for discriminating between breeds (including sheep), it is important to broaden our knowledge of the genetic background of pigmentation. The results may be used in the future to determine the genetic pattern of a breed. Until now, identified candidate genes that have a significant impact on colour type in mammals mainly code for factors located in melanocytes. The proposed candidate genes code for the melanocortin 1 receptor (MC1R), agouti signaling protein (ASIP), tyrosinase-related protein 1 (TYRP1), microphthalmia-associated transcription factor MITF, and v-kit Hardy–Zuckerman 4 feline sarcoma viral oncogene homologue (KIT). However, there is still no conclusive evidence of established polymorphisms for specific coat colour types in sheep.

Genome-Wide Association Studies Identify Candidate Genes for Coat Color and Mohair Traits in the Iranian Markhoz Goat

The Markhoz goat provides an opportunity to study the genetics underlying coat color and mohair traits of an Angora type goat using genome-wide association studies (GWAS). This indigenous Iranian breed is valued for its quality mohair used in ceremonial garments and has the distinction of exhibiting an array of coat colors including black, brown, and white. Here, we performed 16 GWAS for different fleece (mohair) traits and coat color in 228 Markhoz goats sampled from the Markhoz Goat Research Station in Sanandaj, Kurdistan province, located in western Iran using the Illumina Caprine 50K beadchip. The Efficient Mixed Model Linear analysis was used to identify genomic regions with potential candidate genes contributing to coat color and mohair characteristics while correcting for population structure. Significant associations to coat color were found within or near the ASIP, ITCH, AHCY, and RALY genes on chromosome 13 for black and brown coat color and the KIT and PDGFRA genes on chromosome 6 for white coat color. Individual mohair traits were analyzed for genetic association along with principal components that allowed for a broader perspective of combined traits reflecting overall mohair quality and volume. A multitude of markers demonstrated significant association to mohair traits highlighting potential candidate genes of POU1F1 on chromosome 1 for mohair quality, MREG on chromosome 2 for mohair volume, DUOX1 on chromosome 10 for yearling fleece weight, and ADGRV1 on chromosome 7 for grease percentage. Variation in allele frequencies and haplotypes were identified for coat color and differentiated common markers associated with both brown and black coat color. This demonstrates the potential for genetic markers to be used in future breeding programs to improve selection for coat color and mohair traits. Putative candidate genes, both novel and previously identified in other species or breeds, require further investigation to confirm phenotypic causality and potential epistatic relationships.

Revealing the selection history of adaptive loci using genome-wide scans for selection: an example from domestic sheep

Background

One of the approaches to detect genetics variants affecting fitness traits is to identify their surrounding genomic signatures of past selection. With established methods for detecting selection signatures and the current and future availability of large datasets, such studies should have the power to not only detect these signatures but also to infer their selective histories. Domesticated animals offer a powerful model for these approaches as they adapted rapidly to environmental and human-mediated constraints in a relatively short time. We investigated this question by studying a large dataset of 542 individuals from 27 domestic sheep populations raised in France, genotyped for more than 500,000 SNPs.

Results

Population structure analysis revealed that this set of populations harbour a large part of European sheep diversity in a small geographical area, offering a powerful model for the study of adaptation. Identification of extreme SNP and haplotype frequency differences between populations listed 126 genomic regions likely affected by selection. These signatures revealed selection at loci commonly identified as selection targets in many species (“selection hotspots”) including ABCG2, LCORL/NCAPG, MSTN, and coat colour genes such as ASIP, MC1R, MITF, and TYRP1. For one of these regions (ABCG2, LCORL/NCAPG), we could propose a historical scenario leading to the introgression of an adaptive allele into a new genetic background. Among selection signatures, we found clear evidence for parallel selection events in different genetic backgrounds, most likely for different mutations. We confirmed this allelic heterogeneity in one case by resequencing the MC1R gene in three black-faced breeds.

Conclusions

Our study illustrates how dense genetic data in multiple populations allows the deciphering of evolutionary history of populations and of their adaptive mutations.

Electronic supplementary material

The online version of this article (10.1186/s12864-018-4447-x) contains supplementary material, which is available to authorized users.

Mutations in different pigmentation genes are associated with parallel melanism in island flycatchers

The independent evolution of similar traits across multiple taxa provides some of the most compelling evidence of natural selection. Little is known, however, about the genetic basis of these convergent or parallel traits: are they mediated by identical or different mutations in the same genes, or unique mutations in different genes? Using a combination of candidate gene and reduced representation genomic sequencing approaches, we explore the genetic basis of and the evolutionary processes that mediate similar plumage colour shared by isolated populations of the Monarcha castaneiventris flycatcher of the Solomon Islands. A genome-wide association study (GWAS) that explicitly controlled for population structure revealed that mutations in known pigmentation genes are the best predictors of parallel plumage colour. That is, entirely black or melanic birds from one small island share an amino acid substitution in the melanocortin-1 receptor (MC1R), whereas similarly melanic birds from another small island over 100 km away share an amino acid substitution in a predicted binding site of agouti signalling protein (ASIP). A third larger island, which separates the two melanic populations, is inhabited by birds with chestnut bellies that lack the melanic MC1R and ASIP allelic variants. Formal FST outlier tests corroborated the results of the GWAS and suggested that strong, directional selection drives the near fixation of the MC1R and ASIP variants across islands. Our results, therefore, suggest that selection acting on different mutations with large phenotypic effects can drive the evolution of parallel melanism, despite the relatively small population size on islands.

Deciphering signature of selection affecting beef quality traits in Angus cattle

Artificial selection towards a desired phenotype/trait has modified the genomes of livestock dramatically that generated breeds that greatly differ in morphology, production and environmental adaptation traits. Angus cattle are among the famous cattle breeds developed for superior beef quality. This paper aimed at exploring genomic regions under selection in Angus cattle that are associated with meat quality traits and other associated phenotypes. The whole genome of 10 Angus cattle was compared with 11 Hanwoo (A-H) and 9 Jersey (A-J) cattle breeds using a cross-population composite likelihood ratio (XP-CLR) statistical method. The top 1% of the empirical distribution was taken as significant and annotated using UMD3.1. As a result, 255 and 210 genes were revealed under selection from A-H and A-J comparisons, respectively. The WebGestalt gene ontology analysis resulted in sixteen (A-H) and five (A-J) significantly enriched KEGG pathways. Several pathways associated with meat quality traits (insulin signaling, type II diabetes mellitus pathway, focal adhesion pathway, and ECM-receptor interaction), and feeding efficiency (olfactory transduction, tight junction, and metabolic pathways) were enriched. Genes affecting beef quality traits (e.g., FABP3, FTO, DGAT2, ACS, ACAA2, CPE, TNNI1), stature and body size (e.g., PLAG1, LYN, CHCHD7, RPS20), fertility and dystocia (e.g., ESR1, RPS20, PPP2R1A, GHRL, PLAG1), feeding efficiency (e.g., PIK3CD, DNAJC28, DNAJC3, GHRL, PLAG1), coat color (e.g., MC1-R) and genetic disorders (e.g., ITGB6, PLAG1) were found to be under positive selection in Angus cattle. The study identified genes and pathways that are related to meat quality traits and other phenotypes of Angus cattle. The findings in this study, after validation using additional or independent dataset, will provide useful information for the study of Angus cattle in particular and beef cattle in general.

Analysis of the genetic variation in mitochondrial DNA, Y-chromosome sequences, and MC1R sheds light on the ancestry of Nigerian indigenous pigs

Background

The history of pig populations in Africa remains controversial due to insufficient evidence from archaeological and genetic data. Previously, a Western ancestry for West African pigs was reported based on loci that are involved in the determination of coat color. We investigated the genetic diversity of Nigerian indigenous pigs (NIP) by simultaneously analyzing variation in mitochondrial DNA (mtDNA), Y-chromosome sequence and the melanocortin receptor 1 (MC1R) gene.

Results

Median-joining network analysis of mtDNA D-loop sequences from 201 NIP and previously characterized loci clustered NIP with populations from the West (Europe/North Africa) and East/Southeast Asia. Analysis of partial sequences of the Y-chromosome in 57 Nigerian boars clustered NIP into lineage HY1. Finally, analysis of MC1R in 90 NIP resulted in seven haplotypes, among which the European wild boar haplotype was carried by one individual and the European dominant black by most of the other individuals (93%). The five remaining unique haplotypes differed by a single synonymous substitution from European wild type, European dominant black and Asian dominant black haplotypes.

Conclusions

Our results demonstrate a European and East/Southeast Asian ancestry for NIP. Analyses of MC1R provide further evidence. Additional genetic analyses and archaeological studies may provide further insights into the history of African pig breeds. Our findings provide a valuable resource for future studies on whole-genome analyses of African pigs.

Electronic supplementary material

The online version of this article (doi:10.1186/s12711-017-0326-1) contains supplementary material, which is available to authorized users.

The Relationship between MC1R Mutation and Plumage Color Variation in Pigeons

The polymorphisms of MC1R gene play a crucial role in coat color variation in mammals; however, the relationship is still unclear in pigeons. In this study, we sequenced 741 bp fragment of the MC1R for 39 individuals with five plumage color patterns (gray plumage, n = 12; black plumage, n = 9; white plumage, n = 3; spotted plumage, n = 12; red plumage, n = 3). A total of three single nucleotide polymorphisms (SNPs) were detected, including G199A, G225A, and A466G, which subsequently determined four haplotypes (H1–H4). Among them, H1 is the predominant haplotype. Association analysis revealed that H1 and H3 were significantly associated with the black plumage trait (P < 0.05), while the H4 was significantly associated with gray plumage trait (P < 0.05). Furthermore, only diplotype H1H1 was significantly associated with black and gray traits of pigeons. Collectively, our study suggested an association between genetic variation of MC1R and plumage color in pigeon.

Identification and pharmacological analyses of eight naturally occurring caprine melanocortin-1 receptor mutations in three different goat breeds

The melanocortin-1 receptor (MC1R) belongs to the family of seven transmembrane G protein-coupled receptors and plays a central role in animal coat color. We have sequenced the full coding region of 954bp of the MC1R gene in 72 goats of three breeds with different coat colors and identified five missense mutations (K226E, F250V, G255D, V265I, and C267W) and one silent mutation (A61A), among which two haplotypes with complete linkage disequilibrium (A61A and F250V, G255D and V265I) were found. We performed detailed functional studies on the six single and two double mutations in transiently transfected HEK293T cells. We found that none of the mutants had decreased cell surface expression. However, all the mutants except A61A had decreased constitutive activities in the cAMP pathway. Five mutations (F250V, G255D, G267W, A61A/F250V, G255D/V265I) exhibited significant defects in ligand binding and consequent agonist-induced cAMP signaling and ERK1/2 activation. Additionally, K226E, with normal ligand binding affinity and cAMP signaling, showed a significant defect in ERK1/2 activation, exhibiting biased signaling. Co-expression studies showed that the five defective mutants did not affect wild-type MC1R signaling, hence they were not dominant negative. In summary, we provided detailed data of these goat MC1R mutations leading to a better understanding of the role of MC1R mutation and coat color in goats.

Allelic variation of melanocortin-1 receptor locus in Saudi indigenous sheep exhibiting different color coats

Objective

This study was designed to characterize the DNA polymorphisms of the melanocortin-1 receptor (MC1R) gene in indigenous Saudi Arabian sheep breeds exhibiting different color coats, along with individuals of the Sawaknee breed, an exotic sheep imported from Sudan.

Methods

The complete coding region of MC1R gene including parts of 3′ and 5′ untranslated regions was amplified and sequenced from three the indigenous Saudi sheep; Najdi (generally black, n = 41), Naeimi (generally white with brown faces, n = 36) and Herri (generally white, n = 18), in addition to 13 Sawaknee sheep.

Results

Five single nucleotide polymorphisms (SNPs) were detected in the MC1R gene: two led to nonsynonymous mutations (c.218 T>A, p.73 Met>Lys and c.361 G>A, p.121 Asp>Asn) and three led to synonymous mutations (c.429 C>T, p.143 Tyr>Tyr; c.600 T>G, p.200 Leu>Leu, and c.735 C>T, p.245 Ile>Ile). Based on these five SNPs, eight haplotypes representing MC1R E^d and E⁺ alleles were identified among the studied sheep breeds. The most common haplotype (H3) of the dominant E^d allele was associated with either black or brown coat color in Najdi and Sawaknee sheep, respectively. Two other haplotypes (H6 and H7) of E^d allele, with only the nonsynonymous mutation A218T, were detected for the first time in Saudi indigenous sheep.

Conclusion

In addition to investigating the MC1R allelic variation in Saudi indigenous sheep populations, the present study supports the assumption that the two independent nonsynonymous Met73Lys and Asp121Asn mutations in MC1R gene are associated with black or red coat colors in sheep breeds.