Agouti-Signaling Protein and Melanocortin-1-Receptor Mutations Associated with Coat Color Phenotypes in Fallow Deer (Dama dama)

Four dominant coat color phenotypes are found in fallow deer (Dama dama). Brown is the most common. Black, menil, and white occur with varying frequencies. In order to gain insights into the molecular genetic background of these phenotypes, 998 fallow animals (772 brown, 62 black, 126 menil, and 38 white) were examined for mutations in the ASIP, MC1R, TYR, and SLC45A2 genes. In ASIP, two mutations (ASIP-M-E2, located at the boundary from exon 2 to intron 2; and ASIP-M-E3, an InDel of five nucleotides) were found, leading to black fallow deer being either homozygous or heterozygous in combination. There were also two mutations found in MC1R. Whereby the mutation MC1R-M1 (leucine to proline, L48P) homozygous leads to a white coat, while the mutation MC1R-M2 (glycine to aspartic acid, G236D) homozygous is associated with the menil phenotype. When both mutations occur together in a heterozygous character state, it results in a menil coat. Since the mutations in the two genes are only present alternatively, 36 genotypes can be identified that form color clusters to which all animals can be assigned. No mutations were found in the TYR and SLC45A2 genes. Our investigations demonstrate that the four dominant coat colors in fallow deer can be explained by ASIP and MC1R mutations only.

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.

Development of a 76k Alpaca (Vicugna pacos) Single Nucleotide Polymorphisms (SNPs) Microarray

Small farm producers’ sustenance depends on their alpaca herds and the production of fiber. Genetic improvement of fiber characteristics would increase their economic benefits and quality of life. The incorporation of molecular marker technology could overcome current limitations for the implementation of genetic improvement programs. Hence, the aim of this project was the generation of an alpaca single nucleotide polymorphism (SNP) microarray. A sample of 150 Huacaya alpacas from four farms, two each in Puno and Cerro de Pasco were used for SNP discovery by genotyping by sequencing (GBS). Reduced representation libraries, two per animal, were produced after DNA digestion with ApeK1 and double digestion with Pst1-Msp1. Ten alpaca genomes, sequenced at depths between 12× to 30×, and the VicPac3.1 reference genome were used for read alignments. Bioinformatics analysis discovered 76,508 SNPs included in the microarray. Candidate genes SNPs (302) for fiber quality and color are also included. The microarray SNPs cover 90.5% of the genome length with a density of about 39 ± 2.51 SNPs/Mb of DNA at an average interval of 26.45 ± 18.57 kbp. The performance was evaluated by genotyping 30 family trios and comparing them to their pedigrees, as well as comparing microarray to GBS genotypes. Concordance values of 0.93 and 0.94 for ApeK1 and Pst1-Msp1 generated SNPs were observed. Similarly, 290 fiber quality and color candidate gene SNPs were validated. Availability of this microarray will facilitate genome-wide association studies, marker-assisted selection and, in time, genomic selection.

Comprehensive genetic testing combined with citizen science reveals a recently characterized ancient MC1R mutation associated with partial recessive red phenotypes in dog

Background

The Melanocortin 1 Receptor (MC1R) plays a central role in regulation of coat color determination in various species and is commonly referred to as the “E (extension) Locus”. Allelic variation of the MC1R gene is associated with coat color phenotypes E^M (melanistic mask), E^G (grizzle/domino) and e^1–3 (recessive red) in dogs. In addition, a previous study of archeological dog specimens over 10,000 years of age identified a variant p.R301C in the MC1R gene that may have influenced coat color of early dogs.

Results

Commercial genotyping of 11,750 dog samples showed the R301C variant of the MC1R gene was present in 35 breeds or breed varieties, at an allele frequency of 1.5% in the tested population. We detected no linkage disequilibrium between R301C and other tested alleles of the E locus. Based on current convention we propose that R301C should be considered a novel allele of the E locus, which we have termed e^A for “e ancient red”. Phenotype analysis of owner-provided dog pictures reveals that the e^A allele has an impact on coat color and is recessive to wild type E and dominant to the e alleles. In dominant black (K^B/*) dogs it can prevent the phenotypic expression of the K locus, and the expressed coat color is solely determined by the A locus. In the absence of dominant black, e^A/e^A and e^A/e genotypes result in the coat color patterns referred to in their respective breed communities as domino in Alaskan Malamute and other Spitz breeds, grizzle in Chihuahua, and pied in Beagle.

Conclusions

This study demonstrates a large genotype screening effort to identify the frequency and distribution of the MC1R R301C variant, one of the earliest mutations captured by canine domestication, and citizen science empowered characterization of its impact on coat color.

Supplementary Information

The online version contains supplementary material available at 10.1186/s40575-020-00095-7.

Utility of genetic variation in coat color genes to distinguish wild, domestic and hybrid South American camelids for forensic and judicial applications

A molecular genetic protocol for distinguishing pure and hybrid South American camelids was developed to provide strong, quantifiable, and unbiased species identification. We detail the application of the approach in the context of a criminal case in the Andes Mountains of central Chile where the defendants were alleged to have illegally hunted three wild guanacos (Lama guanicoe), as opposed to hybrid domestic llama (Lama glama)/wild guanaco crosses, which are unregulated. We describe a workflow that differentiates among wild, domestic and hybrid South American camelids (Lama versus Vicugna) based on mitochondrial cytochrome b genetic variation (to distinguish between Lama and Vicugna), and MC1R and exon 4 variation of the ASIP gene (to differentiate wild from domestic species). Additionally, we infer the population origin and sex of each of the three individuals from a panel of 15 autosomal microsatellite loci and the presence or absence of the SRY gene. Our analyses strongly supported the inference that the confiscated carcasses corresponded with 2 male and 1 female guanacos that were hunted illegally. Statistical power analyses suggested that there was an extremely low probability of misidentifying domestic camelids as wild camelids (an estimated 0 % Type I error rate), or using more conservative approached a 1.17 % chance of misidentification of wild species as domestic camelids (Type II error). Our case report and methodological and analytical protocols demonstrate the power of genetic variation in coat color genes to identify hybrids between wild and domestic camelid species and highlight the utility of the approach to help combat illegal wildlife hunting and trafficking.

Polymorphisms in MC1R and ASIP Genes are Associated with Coat Color Variation in the Arabian Camel

Pigmentation in mammals is primarily determined by the distribution of eumelanin and pheomelanin, the ratio of which is mostly controlled by the activity of melanocortin 1 receptor (MC1R) and agouti signaling protein (ASIP) genes. Using 91 animals from 10 Arabian camel populations, that included the 4 predominant coat color phenotypes observed in the dromedary (light brown, dark brown, black, and white), we investigated the effects of the MC1R and ASIP sequence variants and identified candidate polymorphisms associated with coat color variation. In particular, we identified a single nucleotide polymorphism (SNP), found in the coding region of MC1R (901C/T), linked to the white coat color, whereas a 1-bp deletion (23delT/T) and a SNP (25G/A) in exon 2 of ASIP are associated with both black and dark-brown coat colors. Our results also indicate support that the light-brown coat color is likely the ancestral coat color for the dromedary. These sequence variations at the MC1R and ASIP genes represent the first documented evidence of candidate polymorphisms associated with Mendelian traits in the dromedary.

Chromosome-Level Alpaca Reference Genome VicPac3.1 Improves Genomic Insight Into the Biology of New World Camelids

The development of high-quality chromosomally assigned reference genomes constitutes a key feature for understanding genome architecture of a species and is critical for the discovery of the genetic blueprints of traits of biological significance. South American camelids serve people in extreme environments and are important fiber and companion animals worldwide. Despite this, the alpaca reference genome lags far behind those available for other domestic species. Here we produced a chromosome-level improved reference assembly for the alpaca genome using the DNA of the same female Huacaya alpaca as in previous assemblies. We generated 190X Illumina short-read, 8X Pacific Biosciences long-read and 60X Dovetail Chicago^® chromatin interaction scaffolding data for the assembly, used testis and skin RNAseq data for annotation, and cytogenetic map data for chromosomal assignments. The new assembly VicPac3.1 contains 90% of the alpaca genome in just 103 scaffolds and 76% of all scaffolds are mapped to the 36 pairs of the alpaca autosomes and the X chromosome. Preliminary annotation of the assembly predicted 22,462 coding genes and 29,337 isoforms. Comparative analysis of selected regions of the alpaca genome, such as the major histocompatibility complex (MHC), the region involved in the Minute Chromosome Syndrome (MCS) and candidate genes for high-altitude adaptations, reveal unique features of the alpaca genome. The alpaca reference genome VicPac3.1 presents a significant improvement in completeness, contiguity and accuracy over VicPac2 and is an important tool for the advancement of genomics research in all New World camelids.

TYR Gene in Llamas: Polymorphisms and Expression Study in Different Color Phenotypes

Tyrosinase, encoded by TYR gene, is an enzyme that plays a major role in mammalian pigmentation. It catalyzes the oxidation of L-dihydroxy-phenylalanine (DOPA) to DOPA quinone, a precursor of both types of melanin: eumelanin and pheomelanin. TYR is commonly known as the albino locus since mutations in this gene result in albinism in several species. However, many other TYR mutations have been found to cause diluted phenotypes, like the Himalayan or chinchilla phenotypes in mice. The llama (Lama glama) presents a wide variety of coat colors ranging from non-diluted phenotypes (eumelanic and pheomelanic), through different degrees of dilution, to white. To investigate the possible contribution of TYR gene to coat color variation in llamas, we sequenced TYR exons and their flanking regions and genotyped animals with diluted, non-diluted, and white coat, including three blue-eyed white individuals. Moreover, we analyzed mRNA expression levels in skin biopsies by qPCR. TYR coding region presented nine SNPs, of which three were non-synonymous, c.428A > G, c.859G > T, and c.1490G > T. We also identified seven polymorphisms in non-coding regions, including two microsatellites, an homopolymeric repeat, and five SNPs: one in the promoter region (c.1-26C > T), two in the 3'-UTR, and two flanking the exons. Although no complete association was found between coat color and SNPs, c.1-26C > T was partially associated to diluted phenotypes. Additionally, the frequency of the G allele from c.428A > G was significantly higher in white compared to non-diluted. Results from qPCR showed that expression levels of TYR in white llamas were significantly lower (p < 0.05) than those in diluted and non-diluted phenotypes. Screening for variation in regulatory regions of TYR did not reveal polymorphisms that explain such differences. However, data from this study showed that TYR expression levels play a role in llama pigmentation.

Comparative FISH-Mapping of MC1R, ASIP, and TYRP1 in New and Old World Camelids and Association Analysis With Coat Color Phenotypes in the Dromedary (Camelus dromedarius)

Melanocortin 1 receptor (MC1R), the agouti signaling protein (ASIP), and tyrosinase related protein 1 (TYRP1) are among the major regulators of pigmentation in mammals. Recently, MC1R and ASIP sequence variants were associated with white and black/dark brown coat colors, respectively, in the dromedary. Here we confirmed this association by independent sequencing and mutation discovery of MC1R and ASIP coding regions and by TaqMan genotyping in 188 dromedaries from Saudi Arabia and United States, including 38 black, 53 white, and 97 beige/brown/red animals. We showed that heterozygosity for a missense mutation c.901C > T in MC1R is sufficient for the white coat color suggesting a possible dominant negative effect. Likewise, we confirmed that the majority of black dromedaries were homozygous for a frameshift mutation in ASIP exon 2, except for 4 animals, which were heterozygous. In search for additional mutations underlying the black color, we identified another frameshift mutation in ASIP exon 4 and 6 new variants in MC1R including a significantly associated SNP in 3'UTR. In pursuit of sequence variants that may modify dromedary wild-type color from dark-reddish brown to light beige, we identified 4 SNPs and one insertion in TYRP1 non-coding regions. However, none of these were associated with variations in wild-type colors. Finally, the three genes were cytogenetically mapped in New World (alpaca) and Old World (dromedary and Bactrian camel) camelids. The MC1R was assigned to chr21, ASIP to chr19 and TYRP1 to chr4 in all 3 species confirming extensive conservation of camelid karyotypes. Notably, while the locations of ASIP and TYRP1 were in agreement with human-camelid comparative map, mapping MC1R identified a new evolutionary conserved synteny segment between camelid chromosome 21 and HSA16. The findings contribute to coat color genomics and the development of molecular tests in camelids and toward the chromosome level reference assemblies of camelid genomes.

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.

Identification of a novel microRNA important for melanogenesis in alpaca (Vicugna pacos)

The molecular mechanisms underlying the formation of coat colors in animals are poorly understood. Recent studies have demonstrated that microRNA play important roles in the control of melanogenesis and coat color in mammals. In a previous study, we characterized the miRNA expression profiles in alpaca skin with brown and white coat color and identified a novel miRNA (named lpa-miR-nov-66) that is expressed significantly higher in white skin compared to brown skin. The present study was conducted to determine the functional roles of this novel miRNA in the regulation of melanogenesis in alpaca melanocytes. lpa-miR-nov-66 is predicted to target the soluble guanylate cyclase (sGC) gene based on presence of a binding site in the sGC coding sequence (CDS). Overexpression of lpa-miR-nov-66 in alpaca melanocyes upregulated the expression of sGC both at the mRNA and protein level. Overexpression of lpa-miR-nov-66 in melanocyes also resulted in decreased expression of key melanogenic genes including tyrosinase (TYR), tyrosinase related protein 1 (TYRP1), and microphthalmia transcription factor (MITF). Our ELISA assays showed increased cyclic guanosine monophosphate (cGMP) but decreased cyclic adenosine monophosphate (cAMP) production in melanocytes overexpressing lpa-miR-nov-66. In addition, overexpression of lpa-miR-nov-66 also reduced melanin production in cultured melanocytes. Results support a role of lpa-miR-nov-66 in melanocytes by directly or indirectly targeting , which regulates melanogenesis via the cAMP pathway.

Expression and tissue distribution of hepatocyte growth factor (HGF) and its receptor (c-Met) in alpacas (Vicugna pacos) skins associated with white and brown coat colors

Hepatocyte growth factor (HGF)/c-Met signaling has been considered as a key pathway in both melanocyte development and melanogenesis. To understand better the expression patterns and tissue distribution characterization of HGF and its receptor c-Met in skin of white versus brown alpaca (Vicugna pacos), we detected the tissue distribution of HGF and c-Met using immunohistochemistry and analyzed the expression patterns by using Western blot and quantitative real time PCR (qPCR). Immunohistochemistry analysis demonstrated that HGF staining robustly increased in the dermal papilla and mesenchymal cells of white alpaca skin compared with that of brown. However, c-Met staining showed strongly positive result, particularly inhair matrix and root sheath in brown alpaca skin. Western blot and qPCR results suggested that HGF and c-Met were expressed at significantly high levels in white and brown alpaca skins, respectively, and protein and transcripts possessed the same expression pattern in white and brown alpaca skins. The results suggested that HGF/c-Met signaling functions in alpaca coat color formation offer essential theoretical basis for further exploration of the role of HGF/c-Met signaling in pigment formation.

The alpaca melanocortin 1 receptor: gene mutations, transcripts, and relative levels of expression in ventral skin biopsies

The objectives of the present study were to characterize the MC1R gene, its transcripts and the single nucleotide polymorphisms (SNPs) associated with coat color in alpaca. Full length cDNA amplification revealed the presence of two transcripts, named as F1 and F2, differing only in the length of their 5′-terminal untranslated region (UTR) sequences and presenting a color specific expression. Whereas the F1 transcript was common to white and colored (black and brown) alpaca phenotypes, the shorter F2 transcript was specific to white alpaca. Further sequencing of the MC1R gene in white and colored alpaca identified a total of twelve SNPs; among those nine (four silent mutations (c.126C>A, c.354T>C, c.618G>A, and c.933G>A); five missense mutations (c.82A>G, c.92C>T, c.259A>G, c.376A>G, and c.901C>T)) were observed in coding region and three in the 3′UTR. A 4 bp deletion (c.224 227del) was also identified in the coding region. Molecular segregation analysis uncovered that the combinatory mutations in the MC1R locus could cause eumelanin and pheomelanin synthesis in alpaca. Overall, our data refine what is known about the MC1R gene and provides additional information on its role in alpaca pigmentation.

The genetic inheritance of the blue-eyed white phenotype in alpacas (Vicugna pacos)

White-spotting patterns in mammals can be caused by mutations in the gene KIT, whose protein is necessary for the normal migration and survival of melanocytes from the neural crest. The alpaca (Vicugna pacos) blue-eyed white (BEW) phenotype is characterized by 2 blue eyes and a solid white coat over the whole body. Breeders hypothesize that the BEW phenotype in alpacas is caused by the combination of the gene causing gray fleece and a white-spotting gene. We performed an association study using KIT flanking and intragenic markers with 40 unrelated alpacas, of which 17 were BEW. Two microsatellite alleles at KIT-related markers were significantly associated (P < 0.0001) with the BEW phenotype (bew1 and bew2). In a larger cohort of 171 related individuals, we identify an abundance of an allele (bew1) in gray animals and the occurrence of bew2 homozygotes that are solid white with pigmented eyes. Association tests accounting for population structure and familial relatedness are consistent with a proposed model where these alleles are in linkage disequilibrium with a mutation or mutations that contribute to the BEW phenotype and to individual differences in fleece color.

Identification and characterization of microRNAs in white and brown alpaca skin

AB BACKGROUND: MicroRNAs (miRNAs) are small, non-coding 21-25 nt RNA molecules that play an important role in regulating gene expression. Little is known about the expression profiles and functions of miRNAs in skin and their role in pigmentation. Alpacas have more than 22 natural coat colors, more than any other fiber producing species. To better understand the role of miRNAs in control of coat color we performed a comprehensive analysis of miRNA expression profiles in skin of white versus brown alpacas.

RESULTS

Two small RNA libraries from white alpaca (WA) and brown alpaca (BA) skin were sequenced with the aid of Illumina sequencing technology. 272 and 267 conserved miRNAs were obtained from the WA and BA skin libraries, respectively. Of these conserved miRNAs, 35 and 13 were more abundant in WA and BA skin, respectively. The targets of these miRNAs were predicted and grouped based on Gene Ontology and KEGG pathway analysis. Many predicted target genes for these miRNAs are involved in the melanogenesis pathway controlling pigmentation. In addition to the conserved miRNAs, we also obtained 22 potentially novel miRNAs from the WA and BA skin libraries.

CONCLUSION

This study represents the first comprehensive survey of miRNAs expressed in skin of animals of different coat colors by deep sequencing analysis. We discovered a collection of miRNAs that are differentially expressed in WA and BA skin. The results suggest important potential functions of miRNAs in coat color regulation.

Single nucleotide polymorphisms in the Melanocortin 1 Receptor gene are linked with lightness of fibre colour in Peruvian Alpaca (Vicugna pacos)

Melanocortin 1 receptor (alpha melanocyte stimulating hormone receptor) (MC1R) is a gene-controlling melanogenesis in mammals. However, it is not well characterized in alpacas and its association with colour is not known. The aim of this study was to look for polymorphisms in the MC1R gene in Peruvian Huacaya alpacas and to analyse the relationship between MC1R single nucleotide polymorphisms (SNPs) and the variations in the instrumental measurement of colour of alpaca fibre. Sixty alpaca fibre samples from black, brown, cream and white animals (15 for each colour) were used to extract DNA from hair bulbs. Colour was measured with a spectrophotometer to obtain quantitative values (CieL*a*b*). Sixteen samples, four of each colour group, were sequenced. Eighteen SNP mutations, 10 not previously described, were found in these 16 sequences. Three of them were chosen (c.82A>G, c.865C>T, c.901C>T) to analyse genotypes by PCR-RFLP in the other 44 fibre samples and to determine the association of mutations with instrumental colour. These three polymorphisms showed association with fibre lightness (P < 0.05), although there was no correlation with colour groups.