Skin Transcriptome Profiles Associated with Skin Color in Chickens

Nutritional and medicinal benefits have been attributed to the consumption of tissues from the black-boned chickens in oriental countries. Lueyang black-boned chicken is one of the native chicken breeds. However, some birds may instead have white or lighter skin, which directly causes economic losses every year. Previous studies of pigmentation have focused on a number of genes that may play important roles in coat color regulation. Illumina2000 sequencing technology was used to catalog the global gene expression profiles in the skin of the Lueyang chicken with white versus black skin. A total of 18,608 unigenes were assembled from the reads obtained from the skin of the white and black chickens. A total of 649 known genes were differentially expressed in the black versus white chickens, with 314 genes that were up regulated and 335 genes that were down-regulated, and a total of 162 novel genes were differentially expressed in the black versus white chickens, consisting of 73 genes that were up-regulated (including 4 highly expressed genes that were expressed exclusively in the skin of the black chickens) and 89 genes that were down-regulated. There were also a total of 8 known coat-color genes expressed in previous studies (ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R). In this study, 4 of which showed greater expression in the black chickens, and several were up-regulated, such as KIT, ASIP, TYR and OCA2. To our surprise, KITLG, MITF and MC1R showed no significant difference in expression between the black- and white-skinned chickens, and the expression of TYRP1 was not detected in either skin color. The expression of ASIP, TYR, KIT, TYRP1, OCA2, KITLG, MITF and MC1R was validated by real-time quantitative polymerase chain reaction (qPCR), and the results of the qPCR were consistent with the RNA-seq. This study provides several candidate genes that may be associated with the development of black versus white skin. More importantly, the fact that the MC1R gene showed no significant difference in expression between the black and white chickens is of particular interest for future studies that aim to elucidate its functional role in the regulation of skin color.     

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.     

A comparison of approaches to estimate the inbreeding coefficient and pairwise relatedness using genomic and pedigree data in a sheep population

Genome-wide SNP data provide a powerful tool to estimate pairwise relatedness among individuals and individual inbreeding coefficient. The aim of this study was to compare methods for estimating the two parameters in a Finnsheep population based on genome-wide SNPs and genealogies, separately. This study included ninety-nine Finnsheep in Finland that differed in coat colours (white, black, brown, grey, and black/white spotted) and were from a large pedigree comprising 319 119 animals. All the individuals were genotyped with the Illumina Ovine SNP50K BeadChip by the International Sheep Genomics Consortium. We identified three genetic subpopulations that corresponded approximately with the coat colours (grey, white, and black and brown) of the sheep. We detected a significant subdivision among the colour types (F _ST = 5.4%, P<0.05). We applied robust algorithms for the genomic estimation of individual inbreeding (F _SNP) and pairwise relatedness (Φ _SNP) as implemented in the programs KING and PLINK, respectively. Estimates of the two parameters from pedigrees (F _PED and Φ _PED) were computed using the RelaX2 program. Values of the two parameters estimated from genomic and genealogical data were mostly consistent, in particular for the highly inbred animals (e.g. inbreeding coefficient F>0.0625) and pairs of closely related animals (e.g. the full- or half-sibs). Nevertheless, we also detected differences in the two parameters between the approaches, particularly with respect to the grey Finnsheep. This could be due to the smaller sample size and relative incompleteness of the pedigree for them.We conclude that the genome-wide genomic data will provide useful information on a per sample or pairwise-samples basis in cases of complex genealogies or in the absence of genealogical data.     

Identification of single nucleotide polymorphisms in the agouti signaling protein (ASIP) gene in some goat breeds in tropical and temperate climates

The agouti-signaling protein (ASIP) plays a major role in mammalian pigmentation as an antagonist to melanocortin-1 receptor gene to stimulate pheomelanin synthesis, a major pigment conferring mammalian coat color. We sequenced a 352 bp fragment of ASIP gene spanning part of exon 2 and part of intron 2 in 215 animals representing six goat breeds from Nigeria and the United States: West African Dwarf, predominantly black; Red Sokoto, mostly red; and Sahel, mostly white from Nigeria; black and white Alpine, brown and white Spanish and white Saanen from the US. Twenty haplotypes from nine mutations representing three intronic, one silent and five missense (p.S19R, p.N35K, p.L36V, p.M42L and p.L45W) mutations were identified in Nigerian goats. Approximately 89 % of Nigerian goats carry haplotype 1 (TGCCATCCG) which seems to be the wild type configuration of mutations in this region of the gene. Although we found no association between these polymorphisms in the ASIP gene and coat color in Nigerian goats, in-silico functional analysis predicts putative deleterious functional impact of the p.L45W mutation on the basic amino-terminal domain of ASIP. In the American goats, two intronic mutations, g.293G>A and g.327C>A, were identified in the Alpine breed, although the g.293G>A mutation is common to American and Nigerian goat populations. All Sannen and Sahel goats in this study belong to haplotypes 1 of both populations which seem to be the wild-type composite ASIP haplotype. Overall, there was no clear association of this portion of the ASIP gene interrogated in this study with coat color variation. Therefore, additional genomic analyses of promoter sequence, the entire coding and non-coding regions of the ASIP gene will be required to obtain a definite conclusion.     

QTL and association analysis for skin and fibre pigmentation in sheep provides evidence of a major causative mutation and epistatic effects

The pursuits of white features and white fleeces free of pigmented fibre have been important selection objectives for many sheep breeds. The cause and inheritance of non‐white colour patterns in sheep has been studied since the early 19th century. Discovery of genetic causes, especially those which predispose pigmentation in white sheep, may lead to more accurate selection tools for improved apparel wool. This article describes an extended QTL study for 13 skin and fibre pigmentation traits in sheep. A total of 19 highly significant, 10 significant and seven suggestive QTL were identified in a QTL mapping experiment using an Awassi × Merino × Merino backcross sheep population. All QTL on chromosome 2 exceeded a LOD score of greater than 4 (range 4.4–30.1), giving very strong support for a major gene for pigmentation on this chromosome. Evidence of epistatic interactions was found for QTL for four traits on chromosomes 2 and 19. The ovine TYRP1 gene on OAR 2 was sequenced as a strong positional candidate gene. A highly significant association (P < 0.01) of grandparental haplotypes across nine segregating SNP/microsatellite markers including one non‐synonymous SNP with pigmentation traits could be shown. Up to 47% of the observed variation in pigmentation was accounted for by models using TYRP1 haplotypes and 83% for models with interactions between two QTL probabilities, offering scope for marker‐assisted selection for these traits.     

The brown coat colour of Coppernecked goats is associated with a non‐synonymous variant at the TYRP1 locus on chromosome 8

The recent development of a goat SNP genotyping microarray enables genome‐wide association studies in this important livestock species. We investigated the genetic basis of the black and brown coat colour in Valais Blacknecked and Coppernecked goats. A genome‐wide association analysis using goat SNP50 BeadChip genotypes of 22 cases and 23 controls allowed us to map the locus for the brown coat colour to goat chromosome 8. The TYRP1 gene is located within the associated chromosomal region, and TYRP1 variants cause similar coat colour phenotypes in different species. We thus considered TYRP1 as a strong positional and functional candidate. We resequenced the caprine TYRP1 gene by Sanger and Illumina sequencing and identified two non‐synonymous variants, p.Ile478Thr and p.Gly496Asp, that might have a functional impact on the TYRP1 protein. However, based on the obtained pedigree and genotype data, the brown coat colour in these goats is not due to a single recessive loss‐of‐function allele. Surprisingly, the genotype distribution and the pedigree data suggest that the ⁴⁹⁶Asp allele might possibly act in a dominant manner. The ⁴⁹⁶Asp allele was present in 77 of 81 investigated Coppernecked goats and did not occur in black goats. This strongly suggests heterogeneity underlying the brown coat colour in Coppernecked goats. Functional experiments or targeted matings will be required to verify the unexpected preliminary findings.     

A premature stop codon in the TYRP1 gene is associated with brown coat colour in the European rabbit (Oryctolagus cuniculus)

Classical genetic studies in European rabbits (Oryctolagus cuniculus) suggested the presence of two alleles at the brown coat colour locus: a wild‐type B allele that gives dense black pigment throughout the coat and a recessive b allele that in the homozygous condition (b/b genotype) produces brown rabbits that are unable to develop black pigmentation. In several other species, this locus is determined by mutations in the tyrosinase‐related protein 1 (TYRP1) gene, encoding a melanocyte enzyme needed for the production of dark eumelanin. In this study, we investigated the rabbit TYRP1 gene as a strong candidate for the rabbit brown coat colour locus. A total of 3846 bp of the TYRP1 gene were sequenced in eight rabbits of different breeds and identified 23 single nucleotide polymorphisms (SNPs; 12 in intronic regions, five in exons and six in the 3′‐untranslated region) and an insertion/deletion of 13 bp, in the 3′‐untranslated region, organised in a few haplotypes. A mutation in exon 2 (g.41360196G>A) leads to a premature stop codon at position 190 of the deduced amino acid sequence (p.Trp190ter). Therefore, translation predicts a truncated TYRP1 protein lacking almost completely the tyrosinase domain. Genotyping 203 rabbits of 32 different breeds identified this mutation only in brown Havana rabbits. Its potential functional relevance in disrupting the TYRP1 protein and its presence only in brown animals strongly argue for this non‐sense mutation being a causative mutation for the recessive b allele at the brown locus in Oryctolagus cuniculus.     

A missense mutation in the agouti signaling protein gene (ASIP) is associated with the no light points coat phenotype in donkeys

Background

Seven donkey breeds are recognized by the French studbook and are characterized by a black, bay or grey coat colour including light cream-to-white points (LP). Occasionally, Normand bay donkeys give birth to dark foals that lack LP and display the no light points (NLP) pattern. This pattern is more frequent and officially recognized in American miniature donkeys. The LP (or pangare) phenotype resembles that of the light bellied agouti pattern in mouse, while the NLP pattern resembles that of the mammalian recessive black phenotype; both phenotypes are associated with the agouti signaling protein gene (ASIP).

Findings

We used a panel of 127 donkeys to identify a recessive missense c.349 T > C variant in ASIP that was shown to be in complete association with the NLP phenotype. This variant results in a cysteine to arginine substitution at position 117 in the ASIP protein. This cysteine is highly-conserved among vertebrate ASIP proteins and was previously shown by mutagenesis experiments to lie within a functional site. Altogether, our results strongly support that the identified mutation is causative of the NLP phenotype.

Conclusions

Thus, we propose to name the c.[349 T > C] allele in donkeys, the a^nlp allele, which enlarges the panel of coat colour alleles in donkeys and ASIP recessive loss-of-function alleles in animals.

Electronic supplementary material

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

Coat Color Variation and Pigmentation Gene Expression in Rhesus Macaques (Macaca mulatta)

Light-dark coat color variation is a common aspect of color diversity within and across mammalian taxa. This variation in pelage brightness is associated with aspects of evolutionary ecology, particularly for primates, but little is known about the genetic mechanisms underlying light-dark differences in pelage pigmentation. Previous work, focusing particularly on macaques (Genus Macaca), has found no clear relationship between color variation and coding sequences of key pigmentation genes. This suggests that other loci and/or gene regulatory differences underlie this variation and raises the question of how patterns of gene expression differ in light verses dark hair follicles. Here, we examine relative expression levels of pigmentation genes in hair follicles from free-ranging rhesus macaques (Macaca mulatta) showing stark light-dark coat color variation. We quantified the brightness (reflectance) of plucked hair tufts using a spectrophotometer. We extracted RNA from the follicles and used quantitative RT-PCR to measure the relative amounts of gene product (mRNA) for seven candidate pigmentation genes (MITF, MC1R, MGRN1, ATRN, SLC24A5, TYRP1, and DCT). Expression values were normalized with the house-keeping gene ACTB. All candidate genes were expressed at similar levels in dark, intermediate, and light hair, and thus, light-dark variation in macaque coat color is unlikely to be due to differences in the expression of these key pigmentation genes. This study represents the first examination of gene expression and natural color variation in a non-human primate population. Our results indicate that even in a system, like pigmentation, where a candidate-gene approach is promising, identifying important intra-specific gene regulatory differences remains challenging.     

Sequence characterization of the melanocortin 1 receptor (MC1R) gene in sheep with different coat colours and identification of the putative e allele at the ovine Extension locus

Sequence of the melanocortin 1 receptor (MC1R) gene (the Extension locus) was obtained from a panel of 73 animals belonging to 9 Italian sheep breeds or populations (Appenninica, Bergamasca, Comisana, Cornigliese-like, Delle Langhe, Massese, Merinizzata Italiana, Sarda and Valle del Belice) with different coat colours. Evaluation of the identified polymorphisms on this phenotype was reported with in silico predictions and comparative approaches within and across breeds and across species. Five novel single nucleotide polymorphisms (SNPs), organized in three haplotypes, were detected. Another haplotype, including the two missense mutations already described for the E^D allele, was identified in few Massese sheep. One SNP (c.199C > T) caused a predicted amino acid substitution (p.R67C) in a highly conserved position of the first intracellular loop of the MC1R protein. The same substitution causes recessive pheomelanism in other species. We propose that the p.67C allele represents the recessive e allele at the ovine Extension series that was, so far, not completely recognized in sheep by classical genetic studies. This polymorphism was analysed in a total of 388 sheep of the 9 investigated breeds. The p.67C allele was identified only in the Valle del Belice breed (allele frequency of 21.3% in 176 analysed animals of this breed) in which the presence of epistatic white-determining loci might mask, at least in part, its effects. Confirming the effect of this novel allele on coat colour will lead to new perspectives on the composition of specialized coloured sheep lines.     

Compelling evidence that a single nucleotide substitution in TYRP1 is responsible for coat-colour polymorphism in a free-living population of Soay sheep

Identifying the genes that underlie phenotypic variation in natural populations is a central objective of evolutionary genetics. Here, we report the identification of the gene and causal mutation underlying coat colour variation in a free-living population of Soay sheep (Ovis aries). We targeted tyrosinase-related protein 1 (TYRP1), a positional candidate gene based on previous work that mapped the Coat colour locus to an approximately 15cM window on chromosome 2. We identified a non-synonymous substitution in exon IV that was perfectly associated with coat colour. This polymorphism is predicted to cause the loss of a cysteine residue that is highly evolutionarily conserved and likely to be of functional significance. We eliminated the possibility that this association is due to the presence of strong linkage disequilibrium with an unknown regulatory mutation by demonstrating that there is no difference in relative TYRP1 expression between colour morphs. Analysis of this putative causal mutation in a complex pedigree of more than 500 sheep revealed almost perfect co-segregation with coat colour (chi2-test, p<0.0001, LOD=110.20), and very tight linkage between Coat colour and TYRP1 (LOD=29.50).     

A deletion spanning the promoter and first exon of the hair cycle-specific ASIP transcript isoform in black and tan rabbits

Black and tan animals have tan-coloured ventral body surfaces separated by sharp boundaries from black-coloured dorsal body surfaces. In the a^t mouse mutant, a retroviral 6 kb insertion located in the hair cycle-specific promoter of the murine Asip gene encoding agouti signalling protein causes the black and tan phenotype. In rabbits, three ASIP alleles are thought to exist, including an a^t allele causing a black and tan coat colour that closely resembles the mouse black and tan phenotype. The goal of our study was to identify the functional genetic variant causing the rabbit a^t allele. We performed a WGS-based comparative analysis of the ASIP gene in one black and tan and three wt agouti-coloured rabbits. The analysis identified 75 a^t-associated variants including an 11 kb deletion. The deletion is located in the region of the hair cycle-specific ASIP promoter and thus in a region homologous to the site of the retroviral insertion causing the a^t allele in mice. We observed perfect association of the genotypes at this deletion with the coat colour phenotype in 49 rabbits. The comparative analysis and the previous knowledge about the regulation of ASIP expression suggest that the 11 kb deletion is the most likely causative variant for the black and tan phenotype in rabbits.     

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^Soay and two novel loss‐of‐function variants. We propose to designate the frame‐shift variant TYRP1:c.86_87delGA as b^VS1 and the nonsense variant TYRP1:c.1066C>T as b^VS2. Interestingly, the b^VS1 allele occurs only in local breeds of Switzerland whereas the b^VS2 allele seems to be more widespread across Europe.     

Pigmentation in Black-boned sheep (<Emphasis Type="Italic">Ovis aries</Emphasis>): association with polymorphism of the <Emphasis Type="Italic">Tyrosinase</Emphasis> gene

Measurements were made in Black-boned (n = 40) and normal (n = 23) sheep (Ovis aries) from a flock in Nanping County of Yunnan Province, China, as well as a group (n = 21) of Romney Marsh sheep (O. aries) with the view to explaining the basis of the dark pigmentation occurring in the Black-boned animals. Plasma colour was significantly darker (P < 0.01) in Black-boned sheep than in their normal flock mates, which in turn had significantly darker plasma (P < 0.01) than the Romney Marsh sheep. Similar significant (P < 0.01) differences were measured for plasma tyrosinase activity and both groups of sheep from Nanping County had similar plasma concentrations of glutathione which were significantly smaller (P < 0.01) than for the Romney Marsh sheep.A partial fragment of 750 bp of exon 1 of the gene encoding tyrosinase was constructed and found to contain two silent mutation sites (G192C and C462T) but there was no effect on amino acid sequences of tyrosinase. Using restriction fragment length polymorphism analyses two allelic variants of site G192C were identified giving rise to the genotypes GG, GC and CC; the frequencies of allele G being 0.914, 0.824 and 0.286 in the Black-boned sheep, their flock mates and the Romney Marsh sheep respectively. Plasma tyrosinase activity was similar for genotypes GG and GC and for both genotypes significantly higher (P < 0.05) than for genotype CC. The sheep from Nanping County displayed only the GG and GC genotypes and had predominantly black or black and white coat colour whereas the Romney Marsh sheep were of either genotype GC or CC and exhibited only white coat colouration. It is not appears that the dark pigmentation of the Black-boned sheep arises because of polymorphisms in the exon 1 of tyrosinase gene. However, this result could explain the differences between Black-boned and Romney Marsh sheep but not for differences between Black-boned and Nanping Normal sheep. Moreover, this result has provided evidence of genetic markers in the form of polymorphisms of the tyrosinase gene which may help to find the black traits causing mutations. There would be merit in further studies using histochemical and molecular techniques to elucidate the causes of the dark pigmentation in these Black-boned sheep.