Whole genome sequencing reveals a novel deletion variant in the KIT gene in horses with white spotted coat colour phenotypes

Predicted genetic burden and frequency of phenotype-associated variants in the horse

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Population Analysis Identifies 15 Multi-Variant Dominant White Haplotypes in Horses

The influence of a horse's appearance on health, sentimental and monetary value has driven the desire to understand the etiology of coat color. White markings on the coat define inclusion for multiple horse breeds, but they may disqualify a horse from registration in other breeds. In domesticated horses (Equus caballus), 35 KIT alleles are associated with or cause depigmentation and white spotting. It is a common misconception among the general public that a horse can possess only two KIT variants. To correct this misconception, we used BEAGLE 5.4-phased NGS data to identify 15 haplotypes possessing two or more KIT variants previously associated with depigmentation phenotypes. We sourced photos for 161 horses comprising 12 compound genotypes with three or more KIT variants and employed a standardized method to grade depigmentation, yielding average white scores for each unique compound genotype. We found that 7 of the 12 multi-variant haplotypes resulted in significantly more depigmentation relative to the single-variant haplotypes (ANOVA). It is clear horses can possess more than two KIT variants, and future work aims to document phenotypic variations for each compound genotype.

Spotting the Pattern: A Review on White Coat Color in the Domestic Horse

Traits such as shape, size, and color often influence the economic and sentimental value of a horse. Around the world, horses are bred and prized for the colors and markings that make their unique coat patterns stand out from the crowd. The underlying genetic mechanisms determining the color of a horse's coat can vary greatly in their complexity. For example, only two genetic markers are used to determine a horse's base coat color, whereas over 50 genetic variations have been discovered to cause white patterning in horses. Some of these white-causing mutations are benign and beautiful, while others have a notable impact on horse health. Negative effects range from slightly more innocuous defects, like deafness, to more pernicious defects, such as the lethal developmental defect incurred when a horse inherits two copies of the Lethal White Overo allele. In this review, we explore, in detail, the etiology of white spotting and its overall effect on the domestic horse to Spot the Pattern of these beautiful (and sometimes dangerous) white mutations.

Geometric morphometrics of face profile across horse breeds and within Arabian horses

Horse traits under selection are largely quantitative and affected by multiple genes. Horse face shape is an example of a continuous trait, which due to the reliance on observational assessments, is classified into; "dished", "straight", and "roman-nosed". This categorization is often inadequate to convey the full spectrum of the face shape variation especially for genetic studies. The first objective of the current study was to use geometric morphometric methods to quantitatively phenotype face shapes and examine its variation across horse breeds. The second objective was to analyze the face shape variation within Arabian horses since face shape is (1) favored, valued, and genetically selected in certain lineages (e.g. Egyptian), (2) is evaluated by registries and scored in shows, and (3) in its extreme forms pose health concerns. We digitized landmarks on lateral profile photos, particularly on the dorsal curvature of the rostrum, and subjected these landmarks to Generalized Procrustes Analysis to generate independent shape and size variables which were statistically compared across breeds and within Arabians. Horse breeds varied in nasal curvature, ranging from extremely concave to extremely convex, with over 70 % of horse breeds exhibiting intermediate concavity (i.e., straight profile). Interestingly, Arabian horses possessed the highest diversity in face profile and individuals clustered into three distinct shape sub-groups (one dished and two straight profile clusters). Our quantitative phenotyping method can be the basis of future genetic studies of facial profile within Arabian lineages as a favored traits and potentially manage its extreme forms as a likely genetic disease.

A de novo 2.3 kb structural variant in MITF explains a novel splashed white phenotype in a Thoroughbred family

Splashed white in horses is characterized by extensive white patterning on the legs, face and abdomen and may be accompanied by deafness. To date, seven variants in microphthalmia-associated transcription factor (MITF) and two variants in Paired Box 3 (PAX3) have been identified to explain this phenotype. A splashed white Thoroughbred stallion, whose sire and dam were not patterned, was hypothesized to have a de novo variant leading to his white coat pattern. A whole-genome sequencing candidate gene approach identified two single nucleotide variants (SNVs) in SOX10, four SNVs in MITF and a 2.3 kb deletion in MITF with the alternative allele present in this stallion but absent in the other 18 horses analyzed. All six SNVs were annotated as modifiers and were not further considered. The deletion in MITF (NC_009159.3:g.21555811_21558139delinsAAAT) encompasses exon 9 encoding a part of the helix-loop-helix domain required for DNA binding. Sanger sequencing and parentage testing confirmed that this deletion was a de novo mutation of maternal origin. Consistent with the published nomenclature, we denote this likely causal variant as SW8. Genotyping three of this stallion's offspring identified SW8 only in the nearly all-white foal that was confirmed deaf by brainstem auditory evoked response testing. This foal was also a compound heterozygote for dominant white variants (W20/W22), but to date, W variants alone have not been connected to deafness. SW8 marks the fourth de novo MITF variant in horses reported to cause white patterning. The link between deafness and all MITF variants with and without other variants impacting melanocyte development and function needs to be further explored.

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.

5'UTR Variant in KIT Associated with White Spotting in Horses

Mutations in KIT, a gene that influences melanoblast migration and pigmentation, often result in mammalian white spotting. As of February 2023, over 30 KIT variants associated with white spotting were documented in Equus caballus (horse). Here we report an association of increased white spotting on the skin and coat with a variant in the 5'UTR of KIT (rs1149701677: g.79,618,649A>C). Horses possessing at least one alternate allele demonstrate phenotypic characteristics similar to other KIT mutations: clear borders around unpigmented regions on the body, face, and limbs. Using a quantitative measure of depigmentation, we observed an average white score of 10.70 among individuals with rs1149701677, while the average score of the control, homozygous reference sample was 2.23 (p=1.892e-11, n=109, t-test). The rs1149701677 site has a cross-species conservation score of 3.4, one of the highest scores across the KIT 5'UTR, implying regulatory importance for this site. Ensembl also predicted a "moderately impactful" functional effect for the rs1149701677 variant. We propose that this single nucleotide variant likely alters the regulation of KIT, which in turn may disrupt melanoblast migration causing an increase in white spotting on the coat. Alternatively, the rs1149701677 variant may be in linkage with another nearby variant with an as-yet-undiscovered functional impact. We propose to term this new allele "Holiday White" or W35 based on conventional nomenclature.

Breed Distribution and Allele Frequencies of Base Coat Color, Dilution, and White Patterning Variants across 28 Horse Breeds

Since domestication, horses have been selectively bred for various coat colors and white spotting patterns. To investigate breed distribution, allele frequencies, and potential lethal variants for recommendations on genetic testing, 29 variants within 14 genes were investigated in 11,281 horses from 28 breeds. The recessive chestnut ea allele in melanocortin 1 receptor (MC1R) (p.D84N) was identified in four breeds: Knabstrupper, Paint Horse, Percheron, and Quarter Horse. After filtering for relatedness, ea allele frequency in Knabstruppers was estimated at 0.035, thus illustrating the importance of testing for mate selection for base coat color. The Rocky Mountain Horse breed had the highest allele frequency for two of the dilution variants under investigation (Za.f. = 0.32 and Cha.f. = 0.026); marker-assisted selection in this breed could aid in the production of horses with desirable dilute coats with less severe ocular anomalies caused by the silver (Z) allele. With regard to white patterning, nine horses homozygous for the paired box 3 (PAX3) splashed white 2 (SW2) allele (p.C70Y) and six horses homozygous for the KIT proto-oncogene, receptor tyrosine kinase (KIT) sabino 1 (SB1) allele (ECA3g.79544206A>T) were identified, thus determining they are rare and confirming that homozygosity for SW2 is not embryonic lethal. The KIT dominant white 20 (W20) allele (p.R682H) was identified in all but three breeds: Arabian (n = 151), Icelandic Horse (n = 66), and Norwegian Fjord Horse (n = 90). The role of W20 in pigmentation across breeds is not well understood; given the different selection regimes of the breeds investigated, these data provide justification for further evaluating the functional role of this allele in pigmentation. Here, we present the largest dataset reported for coat color variants in horses to date, and these data highlight the importance of breed-specific studies to inform on the proper use of marker-assisted selection and to develop hypotheses related to pigmentation for further testing in horses.

A KIT Variant Associated with Increased White Spotting Epistatic to MC1R Genotype in Horses (Equus caballus)

Over 40 identified genetic variants contribute to white spotting in the horse. White markings and spotting are under selection for their impact on the economic value of an equine, yet many phenotypes have an unknown genetic basis. Previous studies also demonstrate an interaction between MC1R and ASIP pigmentation loci and white spotting associated with KIT and MITF. We investigated two stallions presenting with a white spotting phenotype of unknown cause. Exon sequencing of the KIT and MITF candidate genes identified a missense variant in KIT (rs1140732842, NC_009146.3:g.79566881T>C, p.T391A) predicted by SIFT and PROVEAN as not tolerated/deleterious. Three independent observers generated an Average Grade of White (AGW) phenotype score for 147 individuals based on photographs. The KIT variant demonstrates a significant QTL association to AGW (p = 3.3 × 10−12). Association with the MC1R Extension locus demonstrated that, although not in LD, MC1R e/e (chestnut) individuals had higher AGW scores than MC1R E/- individuals (p = 3.09 × 10−17). We also report complete linkage of the previously reported KIT W19 allele to this missense variant. We propose to term this variant W34, following the standardized nomenclature for white spotting variants within the equine KIT gene, and report its epistatic interaction with MC1R.

A genome-wide epistatic network underlies the molecular architecture of continuous color variation of body extremities

Deciphering the molecular architecture of coat coloration for a better understanding of the biological mechanisms underlying pigmentation still remains a challenge. We took advantage of a rabbit French experimental population in which both a pattern and a gradient of coloration from white to brown segregated within the himalayan phenotype. The whole experimental design was genotyped using the high density Affymetrix® AxiomOrcun™ SNP Array and phenotyped into 6 different groups ordered from the lighter to the darker. Genome-wide association analyses pinpointed an oligogenic determinism, under recessive and additive inheritance, involving genes already known in melanogenesis (ASIP, KIT, MC1R, TYR), and likely processed pseudogenes linked to ribosomal function, RPS20 and RPS14. We also identified (i) gene-gene interactions through ASIP:MC1R affecting light cream/beige phenotypes while KIT:RPS responsible of dark chocolate/brown colors and (ii) a genome-wide epistatic network involving several others coloration genes such as POT1 or HPS5. Finally, we determined the recessive inheritance of the English spotting phenotype likely involving a copy number variation affecting at least the end of the coding sequence of the KIT gene. Our analyses of coloration as a continuous trait allowed us to go beyond much of the established knowledge through the detection of additional genes and gene-gene interactions that may contribute to the molecular architecture of the coloration phenotype.

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.

Identification of W13 in the American Miniature Horse and Shetland Pony Populations

Coat color is a trait of economic significance in horses. Variants in seven genes have been documented to cause white patterning in horses. Of the 34 variants that have been identified in KIT proto-oncogene, receptor tyrosine kinase (KIT), 27 have only been reported in a single individual or family and thus not all are routinely offered for genetic testing. Therefore, to enable proper use of marker-assisted selection, determining breed specificity for these alleles is warranted. Screening 19 unregistered all-white Shetland ponies for 16 white patterning markers identified 14 individuals whose phenotype could not be explained by testing results. In evaluating other known dominant white variants, 14 horses were heterozygous for W13. W13 was previously only reported in two quarter horses and a family of Australian miniature horses. Genotyping known white spotting variants in 30 owner-reported white animals (25 Miniature Horses and five Shetland ponies) identified two additional W13/N American Miniature Horses. The estimated allele frequency of W13 in the American Miniature Horse was 0.0063 (79 N/N, 1 W13/N) and the allele was not detected in a random sample (n = 59) of Shetland ponies. No homozygous W13 individuals were identified and W13/N ponies had a similar all-white coat with pink skin phenotype, regardless of the other white spotting variants present, demonstrating that W13 results in a Mendelian inherited dominant white phenotype and homozygosity is likely lethal. These findings document the presence of W13 in the American Miniature Horse and Shetland pony populations at a low frequency and illustrate the importance of testing for this variant in additional breeds.

Two Variants of KIT Causing White Patterning in Stock-Type Horses

Over 30 polymorphisms in the KIT Proto-Oncogene Receptor Tyrosine Kinase (KIT) gene have been implicated in white spotting patterns ranging from small areas to full dermal depigmentation in the horse. We performed a candidate-gene exon sequencing approach on KIT and MITF, 2 known causatives of white spotting patterns, within 2 families of horses of unknown white spotting. Family 1 (Fam1, N = 5) consisted of a Quarter Horse stallion and 4 offspring with white spotting pattern ranging from legs, lower ventral, and head regions with jagged borders, to almost complete white. The second family (Fam2, N = 7) consisted of 6 half-sibling American Paint Horse/Quarter Horse and their dam, demonstrating unpigmented limbs with belly spots and an extensive white patterning on the face. This approach resulted in 2 variants significantly associated with familial phenotypes, where Fam1 variant is an indel leading to a frameshift mutation, and Fam2 a non-synonymous SNP. We validated the variants within an unrelated population of horses (Fam2 variant, P = 0.00271944) as well as for protein functional impact with ExPASy, Protter, Phyre2, SMART, PROVEAN, SIFT, and I-TASSER, confirming the reported associations. Fam1 associated variant, deemed W31, alters the protein sequence, leading to an early stop codon truncating the normal amino acid sequence from 972 to just 115 amino acids. Fam2 associated variant, deemed W32, may have a subtle impact on receptor function or could be in linkage with a non-coding or regulatory change creating the mild spotting pattern observed in this family.

Genetic Testing in the Horse

Genetic testing in horses began in the 1960s, when parentage testing using blood group markers became the standard. In the 1990s, parentage testing shifted from evaluating blood groups to DNA testing. The development of genetics and genomics in both human and veterinarian medicine, along with continued technological advances in the last 2 decades, has helped unravel the causal variants for many horse traits. Genetic testing is also now possible for a variety of phenotypic and disease traits and is used to assist in breeding and clinical management decisions. This article describes the genetic tests that are currently available for horses.

MC1R and KIT Haplotypes Associate With Pigmentation Phenotypes of North American Yak (Bos grunniens)

Small numbers of domestic yak (Bos grunniens) were imported to North America in the late 19th century indirectly from the Qinghai-Tibetan Plateau. Coat color of yak is of interest for fiber production, aesthetics, and as a potential indicator of recent hybridization with cattle. North American yak are classified into 3 major coat color patterns depending upon the presence and extent of white markings. They are further classified by nose pigmentation (black or gray). The aim of this study was to identify loci involved in white patterning and nose pigmentation of North American yak. Genotyping by mass spectrometry of markers identified through Sanger and whole-genome sequencing revealed a 388 kb haplotype of KIT associated in a semi-dominant manner with white coloration in this population of yak. This KIT haplotype is similar to both a haplotype found in white-faced Chinese yak and to haplotypes found in cattle but is divergent from other Bos species such as bison, gaur, and banteng. Melanocortin 1 receptor (MC1R) was implicated as a dominant determinant of black nose color with a single haplotype containing 2 missense mutations perfectly associated with the phenotype. The MC1R haplotype associated with black nose pigment is also similar to cattle haplotypes. No cattle studied, however, shared either of the 2 haplotypes associated with color in yak, suggesting these alleles were introgressed into yak before they were imported to North America. These results provide molecular insight into the history of North American yak and information from which breeders can determine possible color outcomes of matings.

A quantitative trait locus on chromosome 2 was identified that accounts for a substantial proportion of phenotypic variance of the yellow plumage color in chicken

Chicken plumage color is an important economical trait in poultry breeding, as triple-yellow indigenous broilers are preferred over western commercial broilers in the Chinese market. However, the studies on the pigmentation of plumage coloration are relatively rare at present. Here, we performed a genome-wide mapping study on an F₂ intercross, whose 2 founders were one hybrid commercial line “High Quality chicken Line A” that originated from the Anak red chicken and one indigenous line “Huiyang Beard” chicken that is a classical “triple-yellow” Chinese indigenous breed. Moreover, we used an automatic colorimeter that can quantitatively assess the colorations in L∗, a∗, and b∗ values. One major quantitative trait locus (QTL) on chromosome 2 was thus identified by both genome-wide association and linkage analyses, which could explain 10 to 20% of the total phenotypic variance of the b∗ measurements of the back plumage color. Using linkage analysis, 2 additional QTL on chromosome 1 and 20 were also found to be significantly associated with the plumage coloration in this cross. With additional samples from Anak red and Huiyang Beard chickens as well as pooled resequencing data from the 2 founders of this cross, we then further narrowed down the QTL regions and identified several candidate genes, such as CABLES1, CHST11, BCL2L1, and CHD22. As the effects of QTL found in this study were substantial, quantitatively measuring the coloration rather than the descriptive measurements provides stronger statistical power for the analyses. In addition, this major QTL on chromosome 2 that was associated with feather pigmentation at the genome-wide level will facilitate the future chicken breeding for yellow plumage color. In conclusions, we mapped 3 associated QTL on chromosome 1, 2, and 20. The candidate genes identified in this study shed light in the genetic basis of yellow plumage color in chicken.

A Genome-Wide Association Analysis in Noriker Horses Identifies a SNP Associated With Roan Coat Color

The roan coat color in horses is characterized by dispersed white hair and dark points. This phenotype segregates in a broad range of horse breeds, while the underlying genetic background is still unknown. Previous studies mapped the roan locus to the KIT gene on equine chromosome 3 (ECA3). However, this association could not be validated across different horse breeds. Performing a genome-wide association analysis (GWAS) in Noriker horses, we identified a single nucleotide polymorphism (SNP) (ECA3:g.79,543.439 A > G) in the intron 17 of the KIT gene. The G -allele of the top associated SNP was present in other roan horses, namely Quarter Horse, Murgese, Slovenian, and Belgian draught horse, while it was absent in a panel of 15 breeds, including 657 non-roan horses. In further 379 gray Lipizzan horses, eight animals exhibited a heterozygous genotype (A/G). Comparative whole-genome sequence analysis of the KIT region revealed two deletions in the downstream region (ECA3:79,533,217_79,533,224delTCGTCTTC; ECA3:79,533,282_79,533,285delTTCT) and a 3 bp deletion combined with 17 bp insertion in intron 20 of KIT (ECA3:79,588,128_79,588,130delinsTTATCTCTATAGTAGTT). Within the Noriker sample, these loci were in complete linkage disequilibrium (LD) with the identified top SNP. Based upon pedigree information and historical records, we were able to trace back the genetic origin of roan coat color to a baroque gene pool. Furthermore, our data suggest allelic heterogeneity and the existence of additional roan alleles in ponies and breeds related to the English Thoroughbred. In order to study the roan phenotype segregating in those breeds, further association and verification studies are required.

Genomic Analysis Suggests KITLG is Responsible for a Roan Pattern in two Pakistani Goat Breeds

The roan coat color pattern is described as the presence of white hairs intermixed with pigmented hairs. This kind of pigmentation pattern has been observed in many domestic species, including the goat. The molecular mechanisms and inheritance that underlie this pattern are known for some species and the KITLG gene has been shown associated with this phenotype. To date, no research effort has been carried out to find the gene(s) that control(s) roan coat color pattern in goats. In the present study, after genotyping with the GoatSNP50 BeadChip, 35 goats that showed a roan pattern and that belonged to two Pakistan breeds (Group A) were analyzed and then compared to 740 goats of 39 Italian and Pakistan goat breeds that did not have the same coat color pattern (Group B). Runs of homozygosity-based and XP-EHH analyses were used to identify unique genomic regions potentially associated with the roan pattern. A total of 3 regions on chromosomes 5, 6, and 12 were considered unique among the group A versus group B comparisons. The A region > 1.7 Mb on chromosome 5 was the most divergent between the two groups. This region contains six genes, including the KITLG gene. Our findings support the hypothesis that the KITLG gene may be associated with the roan phenotype in goats.

A non-synonymous SNP in exon 3 of the KIT gene is responsible for the classic grey phenotype in alpacas (Vicugna pacos)

The alpaca classic grey phenotype is of particular interest to the industry. Until now, there were only indirect data suggesting that the KIT gene was involved in the classic grey phenotype. All exons of KIT in three black and three classic silvergrey alpacas were sequenced. Five non-synonymous SNPs were observed. There was only one SNP found that was present only in the silvergrey alpacas, and this was also the only SNP predicted to be damaging. This variant results in a change of a glycine (Gly) to an arginine (Arg) at amino acid position 126 (c.376G>A), occurring in the second Ig-like domain of the extracellular domain of KIT. Basic protein modelling predicted that this variant is likely destabilising. Therefore, an additional 488 alpacas were genotyped for this SNP using the tetra-primer amplification refractory mutation system PCR (Tetra-primer ARMS-PCR). All classic grey alpacas were observed to be heterozygous, and 99.3% of non-grey dark base colour alpacas were found to be homozygous for the wildtype allele in this position. These results confirm that the classic grey phenotype in alpacas is the result of a c.376G>A (p.Gly126Arg) SNP in exon 3 of KIT. These data also support the hypothesis that the grey phenotype is autosomal dominant and that the mutation is most likely homozygous lethal.

Characterization and expression analysis of KIT and MITF-M genes in llamas and their relation to white coat color

The llama (Lama glama) is a fiber-producing species that presents a wide range of coat colors, among which white is one of the most important for the textile industry. However, there is little information about the molecular mechanisms that control the white phenotype in this species. In domestic mammals, a white coat is usually produced by mutations in the KIT proto-oncogene receptor tyrosine kinase (KIT) and microphthalmia-associated transcription factor (MITF) genes. In this work we have sequenced and described the coding regions of KIT and MITF-M, the melanocyte-specific isoform, and the two transcriptional variants MITF-M(-) and MITF-M(+). Moreover, we studied the expression of these genes in the skin of white and colored llamas. Although no variants were revealed to be associated with white coat color, significant differences between phenotypes were observed in the expression levels of KIT and MITF-M. Interestingly, white llamas expressed less MITF-M(+) than did colored ones, which is consistent with a consequent reduction in the synthesis of melanin. Even though our results indicate that downregulation of KIT and MITF-M expression is involved in white phenotype production in llamas, the causative gene of white coat color remains unknown.

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.

A Frameshift Mutation in KIT is Associated with White Spotting in the Arabian Camel

While the typical Arabian camel is characterized by a single colored coat, there are rare populations with white spotting patterns. White spotting coat patterns are found in virtually all domesticated species, but are rare in wild species. Theories suggest that white spotting is linked to the domestication process, and is occasionally associated with health disorders. Though mutations have been found in a diverse array of species, fewer than 30 genes have been associated with spotting patterns, thus providing a key set of candidate genes for the Arabian camel. We obtained 26 spotted camels and 24 solid controls for candidate gene analysis. One spotted and eight solid camels were whole genome sequenced as part of a separate project. The spotted camel was heterozygous for a frameshift deletion in KIT (c.1842delG, named KIT^W1 for White spotting 1), whereas all other camels were wild-type (KIT⁺/KIT⁺). No additional mutations unique to the spotted camel were detected in the EDNRB, EDN3, SOX10, KITLG, PDGFRA, MITF, and PAX3 candidate white spotting genes. Sanger sequencing of the study population identified an additional five KIT^W1/KIT⁺ spotted camels. The frameshift results in a premature stop codon five amino acids downstream, thus terminating KIT at the tyrosine kinase domain. An additional 13 spotted camels tested KIT⁺/KIT⁺, but due to phenotypic differences when compared to the KIT^W1/KIT⁺ camels, they likely represent an independent mutation. Our study suggests that there are at least two causes of white spotting in the Arabian camel, the newly described KIT^W1 allele and an uncharacterized mutation.