Register      Login
Animal Production Science Animal Production Science Society
Food, fibre and pharmaceuticals from animals
RESEARCH ARTICLE (Open Access)

Identification of six genomic regions associated with pigment intensity in alpacas

Naomi Gray A , Ishani Shah A , David Groth A and Kylie A. Munyard https://orcid.org/0000-0002-5113-8646 A *
+ Author Affiliations
- Author Affiliations

A Curtin Medical School and Curtin Health Innovation Research Institute, Faculty of Health Sciences, Curtin University, Perth, WA, Australia.

* Correspondence to: K.Munyard@curtin.edu.au

Handling Editor: Sue Hatcher

Animal Production Science 63(11) 1052-1062 https://doi.org/10.1071/AN22467
Submitted: 15 December 2022  Accepted: 10 April 2023   Published: 11 May 2023

© 2023 The Author(s) (or their employer(s)). Published by CSIRO Publishing. This is an open access article distributed under the Creative Commons Attribution-NonCommercial-NoDerivatives 4.0 International License (CC BY-NC-ND)

Abstract

Context: Currently, there is a growing consumer demand for more ecologically sustainable practices in the textile industry. Fabric dyeing is highly pollutive, and one way to avoid dyeing is to use naturally coloured fibres. Alpacas exhibit a wide range of fleece colours, thereby making them a good source of fibre for sustainable textile production. Our understanding of the colour genetics of alpacas is improving, but there is still no explanation for all the variation seen in alpaca coat colours.

Aims: To identify a region or regions in the alpaca genome that contribute to differences in pigment intensity.

Methods: Colorimetric analysis using L*a*b* colourspace of fibre from white and black alpacas, that had been genotyped using the Neogen Australasia alpaca coat colour test, was used to classify the samples into intense and dilute groups for each colour. Illumina short-read genotyping by sequencing of the DNA from these alpacas was used to identify single nucleotide polymorphisms that were subsequently used in a case–control genome-wide association study (GWAS) comparing the extreme dilute and intensely pigmented animals.

Key results: Intense eumelanic fibre is darker (P = 0.0003), less red (P = 0.004), and more blue (P = 0.001) than is dilute eumelanic fibre. Intense pheomelanic fibre is darker (P = 1 × 10−7), more red (P = 3 × 10−20), and more yellow (P = 2 × 10−6) than is dilute pheomelanic fibre. The GWAS showed six regions of genome-wide significance. After manual inspection of these six regions, the best candidate region was upstream of KITLG, a gene previously associated with pigment intensity in dogs.

Conclusions: In combination with ASIP genotype, a regulatory mutation in a region upstream of KITLG in alpacas potentially has a major effect on pigment intensity in the species.

Implications: Successful identification of a marker for pigment intensity will allow breeders to select more precisely for breeding and production animals that will assist them in supplying the desired fibre colours to the textile industry.

Keywords: alpaca, ASIP, CNV, colorimetry, colour, genotyping by sequencing, GWAS, KITLG, pigment intensity.


References

Aigner B, Besenfelder U, Müller M, Brem G (2000) Tyrosinase gene variants in different rabbit strains. Mammalian Genome 11, 700–702.
Tyrosinase gene variants in different rabbit strains.Crossref | GoogleScholarGoogle Scholar |

Altschul SF, Gish W, Miller W, Myers EW, Lipman DJ (1990) Basic local alignment search tool. Journal of Molecular Biology 215, 403–410.
Basic local alignment search tool.Crossref | GoogleScholarGoogle Scholar |

Andersson LS, Wilbe M, Viluma A, Cothran G, Ekesten B, Ewart S, Lindgren G (2013) Equine multiple congenital ocular anomalies and silver coat colour result from the pleiotropic effects of mutant PMEL. PLoS ONE 8, e75639
Equine multiple congenital ocular anomalies and silver coat colour result from the pleiotropic effects of mutant PMEL.Crossref | GoogleScholarGoogle Scholar |

Barsh GS (1996) The genetics of pigmentation: from fancy genes to complex traits. Trends in Genetics 12, 299–305.
The genetics of pigmentation: from fancy genes to complex traits.Crossref | GoogleScholarGoogle Scholar |

Bauer A, Kehl A, Jagannathan V, Leeb T (2018) A novel MLPH variant in dogs with coat colour dilution. Animal Genetics 49, 94–97.
A novel MLPH variant in dogs with coat colour dilution.Crossref | GoogleScholarGoogle Scholar |

Baxter LL, Watkins-Chow DE, Pavan WJ, Loftus SK (2019) A curated gene list for expanding the horizons of pigmentation biology. Pigment Cell & Melanoma Research 32, 348–358.
A curated gene list for expanding the horizons of pigmentation biology.Crossref | GoogleScholarGoogle Scholar |

Baye LM, Patrinostro X, Swaminathan S, Beck JS, Zhang Y, Stone EM, Sheffield VC, Slusarski DC (2011) The N-terminal region of centrosomal protein 290 (CEP290) restores vision in a zebrafish model of human blindness. Human Molecular Genetics 20, 1467–1477.
The N-terminal region of centrosomal protein 290 (CEP290) restores vision in a zebrafish model of human blindness.Crossref | GoogleScholarGoogle Scholar |

Beck BB, Phillips JB, Bartram MP, Wegner J, Thoenes M, Pannes A, Sampson J, Heller R, Göbel H, Koerber F, Neugebauer A, Hedergott A, Nürnberg G, Nürnberg P, Thiele H, Altmüller J, Toliat MR, Staubach S, Boycott KM, Valente EM, Janecke AR, Eisenberger T, Bergmann C, Tebbe L, Wang Y, Wu Y, Fry AM, Westerfield M, Wolfrum U, Bolz HJ (2014) Mutation of POC1B in a severe syndromic retinal ciliopathy. Human Mutation 35, 1153–1162.
Mutation of POC1B in a severe syndromic retinal ciliopathy.Crossref | GoogleScholarGoogle Scholar |

Bedell MA, Brannan CI, Evans EP, Copeland NG, Jenkins NA, Donovan PJ (1995) DNA rearrangements located over 100 kb 5’ of the Steel (Sl)-coding region in Steel-panda and Steel-contrasted mice deregulate Sl expression and cause female sterility by disrupting ovarian follicle development. Genes & Development 9, 455–470.
DNA rearrangements located over 100 kb 5’ of the Steel (Sl)-coding region in Steel-panda and Steel-contrasted mice deregulate Sl expression and cause female sterility by disrupting ovarian follicle development.Crossref | GoogleScholarGoogle Scholar |

Besmer P (1991) The kit ligand encoded at the murine Steel locus: a pleiotropic growth and differentiation factor. Current Opinion in Cell Biology 3, 939–946.
The kit ligand encoded at the murine Steel locus: a pleiotropic growth and differentiation factor.Crossref | GoogleScholarGoogle Scholar |

Blaszczyk WM, Arning L, Hoffmann K-P, Epplen JT (2005) A Tyrosinase missense mutation causes albinism in the Wistar rat. Pigment Cell Research 18, 144–145.
A Tyrosinase missense mutation causes albinism in the Wistar rat.Crossref | GoogleScholarGoogle Scholar |

Bonaventure J, Domingues MJ, Larue L (2013) Cellular and molecular mechanisms controlling the migration of melanocytes and melanoma cells. Pigment Cell & Melanoma Research 26, 316–325.
Cellular and molecular mechanisms controlling the migration of melanocytes and melanoma cells.Crossref | GoogleScholarGoogle Scholar |

Brancalion L, Haase B, Wade CM (2022) Canine coat pigmentation genetics: a review. Animal Genetics 53, 3–34.
Canine coat pigmentation genetics: a review.Crossref | GoogleScholarGoogle Scholar |

Broad Institute (2019) Picard toolkit. Available at https://broadinstitute.github.io/picard/

Brunberg E, Andersson L, Cothran G, Sandberg K, Mikko S, Lindgren G (2006) A missense mutation in PMEL17 is associated with the Silver coat color in the horse. BMC Genetics 7, 46
A missense mutation in PMEL17 is associated with the Silver coat color in the horse.Crossref | GoogleScholarGoogle Scholar |

Bultman SJ, Michaud EJ, Woychik RP (1992) Molecular characterization of the mouse agouti locus. Cell 71, 1195–1204.
Molecular characterization of the mouse agouti locus.Crossref | GoogleScholarGoogle Scholar |

Camacho-Hübner A, Richard C, Beermann F (2002) Genomic structure and evolutionary conservation of the tyrosinase gene family from Fugu. Gene 285, 59–68.
Genomic structure and evolutionary conservation of the tyrosinase gene family from Fugu.Crossref | GoogleScholarGoogle Scholar |

Chen S, Zhou Y, Chen Y, Gu J (2018) fastp: an ultra-fast all-in-one FASTQ preprocessor. Bioinformatics 34, i884–i890.
fastp: an ultra-fast all-in-one FASTQ preprocessor.Crossref | GoogleScholarGoogle Scholar |

Costin G-E, Hearing VJ (2007) Human skin pigmentation: melanocytes modulate skin color in response to stress. The FASEB Journal 21, 976–994.
Human skin pigmentation: melanocytes modulate skin color in response to stress.Crossref | GoogleScholarGoogle Scholar |

Cruz A, Yucra A, Gutiérrez GA, Burgos A, Morante R, Gutiérrez JP, Cervantes I, Wurzinger M (2021) Colorimetry analysis of coat color and its relationship with fiber traits in alpacas. Animal 15, 100219
Colorimetry analysis of coat color and its relationship with fiber traits in alpacas.Crossref | GoogleScholarGoogle Scholar |

Danecek P, Bonfield JK, Liddle J, Marshall J, Ohan V, Pollard MO, Whitwham A, Keane T, McCarthy SA, Davies RM, Li H (2021) Twelve years of SAMtools and BCFtools. GigaScience 10, giab008
Twelve years of SAMtools and BCFtools.Crossref | GoogleScholarGoogle Scholar |

Drögemüller C, Philipp U, Haase B, Günzel-Apel A-R, Leeb T (2007) A noncoding melanophilin gene (MLPH) SNP at the splice donor of exon 1 represents a candidate causal mutation for coat color dilution in dogs. Journal of Heredity 98, 468–473.
A noncoding melanophilin gene (MLPH) SNP at the splice donor of exon 1 represents a candidate causal mutation for coat color dilution in dogs.Crossref | GoogleScholarGoogle Scholar |

D’Mello SAN, Finlay GJ, Baguley BC, Askarian-Amiri ME (2016) Signaling pathways in melanogenesis. International Journal of Molecular Sciences 17, 1144
Signaling pathways in melanogenesis.Crossref | GoogleScholarGoogle Scholar |

Feeley NL, Munyard KA (2009) Characterisation of the melanocortin-1 receptor gene in alpaca and identification of possible markers associated with phenotypic variations in colour. Animal Production Science 49, 675–681.
Characterisation of the melanocortin-1 receptor gene in alpaca and identification of possible markers associated with phenotypic variations in colour.Crossref | GoogleScholarGoogle Scholar |

Feeley NL, Bottomley S, Munyard KA (2011) Three novel mutations in ASIP associated with black fibre in alpacas (Vicugna pacos). The Journal of Agricultural Science 149, 529–538.
Three novel mutations in ASIP associated with black fibre in alpacas (Vicugna pacos).Crossref | GoogleScholarGoogle Scholar |

Fleis RI, Scott DW (2010) The microanatomy of healthy skin from alpacas (Vicugna pacos). Journal of Veterinary Diagnostic Investigation 22, 716–719.
The microanatomy of healthy skin from alpacas (Vicugna pacos).Crossref | GoogleScholarGoogle Scholar |

Floriot S, Duchesne A, Grohs C, Hozé C, Deloche MC, Fayolle G, Vilotte JL, Boichard D, Fritz S, Boussaha M (2021) A missense mutation in the FZD7 gene is associated with dilution of the red areas of the coat in Montbéliarde cattle. Animal Genetics 52, 351–355.
A missense mutation in the FZD7 gene is associated with dilution of the red areas of the coat in Montbéliarde cattle.Crossref | GoogleScholarGoogle Scholar |

Frank EN, Hick MVH, Gauna CD, Lamas HE, Renieri C, Antonini M (2006) Phenotypic and genetic description of fibre traits in South American domestic camelids (llamas and alpacas). Small Ruminant Research 61, 113–129.
Phenotypic and genetic description of fibre traits in South American domestic camelids (llamas and alpacas).Crossref | GoogleScholarGoogle Scholar |

Furumura M, Sakai C, Abdel-Malek Z, Barsh GS, Hearing VJ (1996) The interaction of agouti signal protein and melanocyte stimulating hormone to regulate melanin formation in mammals. Pigment Cell Research 9, 191–203.
The interaction of agouti signal protein and melanocyte stimulating hormone to regulate melanin formation in mammals.Crossref | GoogleScholarGoogle Scholar |

GHD (2022) Alpaca: a compact business case. Available at https://agrifutures.com.au/product/alpaca-a-compact-business-case/

Graf J, Hodgson R, van Daal A (2005) Single nucleotide polymorphisms in the MATP gene are associated with normal human pigmentation variation. Human Mutation 25, 278–284.
Single nucleotide polymorphisms in the MATP gene are associated with normal human pigmentation variation.Crossref | GoogleScholarGoogle Scholar |

Grichnik JM, Burch JA, Burchette J, Shea CR (1998) The SCF/KIT pathway plays a critical role in the control of normal human melanocyte homeostasis. Journal of Investigative Dermatology 111, 233–238.
The SCF/KIT pathway plays a critical role in the control of normal human melanocyte homeostasis.Crossref | GoogleScholarGoogle Scholar |

Gutiérrez-Gil B, Wiener P, Williams JL (2007) Genetic effects on coat colour in cattle: dilution of eumelanin and phaeomelanin pigments in an F2-Backcross Charolais × Holstein population. BMC Genetics 8, 56
Genetic effects on coat colour in cattle: dilution of eumelanin and phaeomelanin pigments in an F2-Backcross Charolais × Holstein population.Crossref | GoogleScholarGoogle Scholar |

Hoekstra HE (2006) Genetics, development and evolution of adaptive pigmentation in vertebrates. Heredity 97, 222–234.
Genetics, development and evolution of adaptive pigmentation in vertebrates.Crossref | GoogleScholarGoogle Scholar |

Hédan B, Cadieu E, Botherel N, Dufaure de Citres C, Letko A, Rimbault M, Drögemüller C, Jagannathan V, Derrien T, Schmutz S, Leeb T, André C (2019) Identification of a missense variant in MFSD12 involved in dilution of phaeomelanin leading to white or cream coat color in dogs. Genes-Basel 10, 386
Identification of a missense variant in MFSD12 involved in dilution of phaeomelanin leading to white or cream coat color in dogs.Crossref | GoogleScholarGoogle Scholar |

Ito S, Wakamatsu K (2003) Quantitative analysis of eumelanin and pheomelanin in humans, mice, and other animals: a comparative review. Pigment Cell Research 16, 523–531.
Quantitative analysis of eumelanin and pheomelanin in humans, mice, and other animals: a comparative review.Crossref | GoogleScholarGoogle Scholar |

Kimura S, Hatakeyama T, Koutaka T, Kubo K, Morita S, Eguchi K, Saitoh K, Yamauchi K, Imai S, Kashimura A, Inenaga T, Matsumoto H (2022) PMEL p.Leu18del dilutes coat color of Kumamoto sub-breed of Japanese Brown cattle. BMC Genomics 23, 694
PMEL p.Leu18del dilutes coat color of Kumamoto sub-breed of Japanese Brown cattle.Crossref | GoogleScholarGoogle Scholar |

Lander ES, Botstein D (1989) Mapping mendelian factors underlying quantitative traits using RFLP linkage maps. Genetics 121, 185–199.
Mapping mendelian factors underlying quantitative traits using RFLP linkage maps.Crossref | GoogleScholarGoogle Scholar |

Langmead B, Salzberg SL (2012) Fast gapped-read alignment with Bowtie 2. Nature Methods 9, 357–359.
Fast gapped-read alignment with Bowtie 2.Crossref | GoogleScholarGoogle Scholar |

Ligges U, Maechler M (2003) scatterplot3d – an R package for visualizing multivariate data. Journal of Statistical Software 8, 1–20.
scatterplot3d – an R package for visualizing multivariate data.Crossref | GoogleScholarGoogle Scholar |

Lu D, Willard D, Patel IR, Kadwell S, Overton L, Kost T, Luther M, Chen W, Woychik RP, Wilkison WO, Cone RD (1994) Agouti protein is an antagonist of the melanocyte-stimulating-hormone receptor. Nature 371, 799–802.
Agouti protein is an antagonist of the melanocyte-stimulating-hormone receptor.Crossref | GoogleScholarGoogle Scholar |

Lupton CJ, McColl A, Stobart RH (2006) Fiber characteristics of the Huacaya alpaca. Small Ruminant Research 64, 211–224.
Fiber characteristics of the Huacaya alpaca.Crossref | GoogleScholarGoogle Scholar |

Ly BCK, Dyer EB, Feig JL, Chien AL, Del Bino S (2020) Research techniques made simple: cutaneous colorimetry: a reliable technique for objective skin color measurement. Journal of Investigative Dermatology 140, 3–12.e1.
Research techniques made simple: cutaneous colorimetry: a reliable technique for objective skin color measurement.Crossref | GoogleScholarGoogle Scholar |

Marees AT, de Kluiver H, Stringer S, Vorspan F, Curis E, Marie-Claire C, Derks EM (2018) A tutorial on conducting genome-wide association studies: quality control and statistical analysis. International Journal of Methods in Psychiatric Research 27, e1608
A tutorial on conducting genome-wide association studies: quality control and statistical analysis.Crossref | GoogleScholarGoogle Scholar |

Mariat D, Taourit S, Guérin G (2003) A mutation in the MATP gene causes the cream coat colour in the horse. Genetics Selection Evolution 35, 119
A mutation in the MATP gene causes the cream coat colour in the horse.Crossref | GoogleScholarGoogle Scholar |

Miller CT, Beleza S, Pollen AA, Schluter D, Kittles RA, Shriver MD, Kingsley DM (2007) cis-regulatory changes in kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans. Cell 131, 1179–1189.
cis-regulatory changes in kit ligand expression and parallel evolution of pigmentation in sticklebacks and humans.Crossref | GoogleScholarGoogle Scholar |

Munyard K (2011) ‘Inheritance of white colour in alpacas: identifying the genes involved.’ (Rural Industries Research and Development Corporation)

Neuffer SJ, Cooper CD (2022) Zebrafish syndromic albinism models as tools for understanding and treating pigment cell disease in humans. Cancers 14, 1752
Zebrafish syndromic albinism models as tools for understanding and treating pigment cell disease in humans.Crossref | GoogleScholarGoogle Scholar |

Newton JM, Cohen-Barak O, Hagiwara N, Gardner JM, Davisson MT, King RA, Brilliant MH (2001) Mutations in the human orthologue of the mouse underwhite gene (uw) underlie a new form of oculocutaneous albinism, OCA4. The American Journal of Human Genetics 69, 981–988.
Mutations in the human orthologue of the mouse underwhite gene (uw) underlie a new form of oculocutaneous albinism, OCA4.Crossref | GoogleScholarGoogle Scholar |

Oetting WS, Brilliant MH, King RA (1996) The clinical spectrum of albinism in humans. Molecular Medicine Today 2, 330–335.
The clinical spectrum of albinism in humans.Crossref | GoogleScholarGoogle Scholar |

ØYehaug L, Plahte E, VÅGe DI, Omholt SW (2002) The regulatory basis of melanogenic switching. Journal of Theoretical Biology 215, 449–468.
The regulatory basis of melanogenic switching.Crossref | GoogleScholarGoogle Scholar |

Pillaiyar T, Manickam M, Jung S-H (2017) Recent development of signaling pathways inhibitors of melanogenesis. Cell Signal 40, 99–115.
Recent development of signaling pathways inhibitors of melanogenesis.Crossref | GoogleScholarGoogle Scholar |

Poplin R, Ruano-Rubio V, DePristo MA, Fennell TJ, Carneiro MO, Van der Auwera GA, Kling DE, Gauthier LD, Levy-Moonshine A, Roazen D, Shakir K, Thibault J, Chandran S, Whelan C, Lek M, Gabriel S, Daly MJ, Neale B, MacArthur DG, Banks E (2018) Scaling accurate genetic variant discovery to tens of thousands of samples. bioRxiv 201178
Scaling accurate genetic variant discovery to tens of thousands of samples.Crossref | GoogleScholarGoogle Scholar |

Powell AJ, Moss MJ, Tree LT, Roeder BL, Carleton CL, Campbell E, Kooyman DL (2008) Characterization of the effect of Melanocortin 1 Receptor, a member of the hair color genetic locus, in alpaca (Lama pacos) fleece color differentiation. Small Ruminant Research 79, 183–187.
Characterization of the effect of Melanocortin 1 Receptor, a member of the hair color genetic locus, in alpaca (Lama pacos) fleece color differentiation.Crossref | GoogleScholarGoogle Scholar |

Purcell S, Neale B, Todd-Brown K, Thomas L, Ferreira MAR, Bender D, Maller J, Sklar P, de Bakker PIW, Daly MJ, Sham PC (2007) PLINK: a tool set for whole-genome association and population-based linkage analyses. The American Journal of Human Genetics 81, 559–575.
PLINK: a tool set for whole-genome association and population-based linkage analyses.Crossref | GoogleScholarGoogle Scholar |

R Core Team (2022) ‘R: a language and environment for statistical computing.’ (R Foundation for Statistical Computing: Vienna, Austria). Available at http://www.R-project.org/

Rangwala SH, Kuznetsov A, Ananiev V, Asztalos A, Borodin E, Evgeniev V, Joukov V, Lotov V, Pannu R, Rudnev D, Shkeda A, Weitz EM, Schneider VA (2021) Accessing NCBI data using the NCBI sequence viewer and genome data viewer (GDV). Genome Research 31, 159–169.
Accessing NCBI data using the NCBI sequence viewer and genome data viewer (GDV).Crossref | GoogleScholarGoogle Scholar |

Robbins LS, Nadeau JH, Johnson KR, Kelly MA, Roselli-Rehfuss L, Baack E, Mountjoy KG, Cone RD (1993) Pigmentation phenotypes of variant extension locus alleles result from point mutations that alter MSH receptor function. Cell 72, 827–834.
Pigmentation phenotypes of variant extension locus alleles result from point mutations that alter MSH receptor function.Crossref | GoogleScholarGoogle Scholar |

Robinson JT, Thorvaldsdóttir H, Winckler W, Guttman M, Lander ES, Getz G, Mesirov JP (2011) Integrative genomics viewer. Nature Biotechnology 29, 24–26.
Integrative genomics viewer.Crossref | GoogleScholarGoogle Scholar |

Rural Industries Research and Development Corporation (2016) Alpaca market assessment 2016. Rural Industries Research and Development Corporation. Available at https://www.agrifutures.com.au/product/project-summary-alpaca-market-assessment-2016/

Sayers EW, Bolton EE, Brister JR, Canese K, Chan J, Comeau DC, Connor R, Funk K, Kelly C, Kim S, Madej T, Marchler-Bauer A, Lanczycki C, Lathrop S, Lu Z, Thibaud-Nissen F, Murphy T, Phan L, Skripchenko Y, Tse T, Wang J, Williams R, Trawick BW, Pruitt KD, Sherry ST (2022) Database resources of the national center for biotechnology information. Nucleic Acids Research 50, D20–D26.
Database resources of the national center for biotechnology information.Crossref | GoogleScholarGoogle Scholar |

Schmutz SM, Berryere TG, Ciobanu DC, Mileham AJ, Schmidtz BH, Fredholm M (2004) A form of albinism in cattle is caused by a tyrosinase frameshift mutation. Mammalian Genome 15, 62–67.
A form of albinism in cattle is caused by a tyrosinase frameshift mutation.Crossref | GoogleScholarGoogle Scholar |

Shah I, Gray N, Groth D, Brooks S, Munyard K (2023) A copy number variant near KITLG is associated with the roan pattern in alpacas. Animal Production Science

Slavney AJ, Kawakami T, Jensen MK, Nelson TC, Sams AJ, Boyko AR (2021) Five genetic variants explain over 70% of hair coat pheomelanin intensity variation in purebred and mixed breed domestic dogs. PLoS ONE 16, e0250579
Five genetic variants explain over 70% of hair coat pheomelanin intensity variation in purebred and mixed breed domestic dogs.Crossref | GoogleScholarGoogle Scholar |

Song X, Xu C, Liu Z, Yue Z, Liu L, Yang T, Cong B, Yang F (2017) Comparative transcriptome analysis of mink (Neovison vison) skin reveals the key genes involved in the melanogenesis of black and white coat colour. Scientific Reports 7, 12461
Comparative transcriptome analysis of mink (Neovison vison) skin reveals the key genes involved in the melanogenesis of black and white coat colour.Crossref | GoogleScholarGoogle Scholar |

Tam V, Patel N, Turcotte M, Bossé Y, Paré G, Meyre D (2019) Benefits and limitations of genome-wide association studies. Nature Reviews Genetics 20, 467–484.
Benefits and limitations of genome-wide association studies.Crossref | GoogleScholarGoogle Scholar |

Vaez M, Follett SA, Bed’hom B, Gourichon D, Tixier-Boichard M, Burke T (2008) A single point-mutation within the melanophilin gene causes the lavender plumage colour dilution phenotype in the chicken. BMC Genetics 9, 7
A single point-mutation within the melanophilin gene causes the lavender plumage colour dilution phenotype in the chicken.Crossref | GoogleScholarGoogle Scholar |

Van Buren SL, Minor KM, Grahn RA, Mickelson JR, Grahn JC, Malvick J, Colangelo JR, Mueller E, Kuehnlein P, Kehl A (2020) A third MLPH variant causing coat color dilution in dogs. Genes 11, 639
A third MLPH variant causing coat color dilution in dogs.Crossref | GoogleScholarGoogle Scholar |

Van der Auwera GA, O’Connor BD (2020) ‘Genomics in the cloud: using docker, GATK, and WDL in terra.’ (O’Reilly Media, Inc)

Wakamatsu K, Munyard K, Oddie C, Ito S (2021) Photobleached oxidative degradation of melanins: chemical characterization of melanins present in alpaca fiber. Photochemistry and Photobiology 97, 1493–1497.
Photobleached oxidative degradation of melanins: chemical characterization of melanins present in alpaca fiber.Crossref | GoogleScholarGoogle Scholar |

Wang S, Dvorkin D, Da Y (2012) SNPEVG: a graphical tool for GWAS graphing with mouse clicks. BMC Bioinformatics 13, 319
SNPEVG: a graphical tool for GWAS graphing with mouse clicks.Crossref | GoogleScholarGoogle Scholar |

Wang Y, Viennet C, Robin S, Berthon J-Y, He L, Humbert P (2017) Precise role of dermal fibroblasts on melanocyte pigmentation. Journal of Dermatological Science 88, 159–166.
Precise role of dermal fibroblasts on melanocyte pigmentation.Crossref | GoogleScholarGoogle Scholar |

Wehrle-Haller B (2003) The role of kit-ligand in melanocyte development and epidermal homeostasis. Pigment Cell Research 16, 287–296.
The role of kit-ligand in melanocyte development and epidermal homeostasis.Crossref | GoogleScholarGoogle Scholar |

Weich K, Affolter V, York D, Rebhun R, Grahn R, Kallenberg A, Bannasch D (2020) Pigment intensity in dogs is associated with a copy number variant upstream of KITLG. Genes 11, 75
Pigment intensity in dogs is associated with a copy number variant upstream of KITLG.Crossref | GoogleScholarGoogle Scholar |

Wickham H (2016) ‘ggplot2: elegant graphics for data analysis.’ (Springer-Verlag: New York, NY, USA)

Wu S, Li J, Ma T, Li J, Li Y, Jiang H, Zhang Q (2021) MiR-27a regulates WNT3A and KITLG expression in Cashmere goats with different coat colors. Animal Biotechnology 32, 205–212.
MiR-27a regulates WNT3A and KITLG expression in Cashmere goats with different coat colors.Crossref | GoogleScholarGoogle Scholar |

Yoo JC, Lim Ty, Park JS, Hah Y-S, Park N, Hong S-G, Park J-Y, Yoon T-J (2013) SYT14L, especially its C2 domain, is involved in regulating melanocyte differentiation. Journal of Dermatological Science 72, 246–251.
SYT14L, especially its C2 domain, is involved in regulating melanocyte differentiation.Crossref | GoogleScholarGoogle Scholar |