Seven novel KIT mutations in horses with white coat colour phenotypes

Genomic analysis reveals the association of KIT and MITF variants with the white spotting in swamp buffaloes

BACKGROUND

Swamp-type buffaloes with varying degrees of white spotting are found exclusively in Tana Toraja, South Sulawesi, Indonesia, where spotted buffalo bulls are highly valued in accordance with the Torajan customs. The white spotting depigmentation is caused by the absence of melanocytes. However, the genetic variants that cause this phenotype have not been fully characterized. The objective of this study was to identify the genomic regions and variants responsible for this unique coat-color pattern.

RESULTS

Genome-wide association study (GWAS) and selection signature analysis identified MITF as a key gene based on the whole-genome sequencing data of 28 solid and 39 spotted buffaloes, while KIT was also found to be involved in the development of this phenotype by a candidate gene approach. Alternative candidate mutations included, in addition to the previously reported nonsense mutation c.649 C > T (p.Arg217*) and splice donor mutation c.1179 + 2T > A in MITF, a nonsense mutation c.2028T > A (p.Tyr676*) in KIT. All these three mutations were located in the genomic regions that were highly conserved exclusively in Indonesian swamp buffaloes and they accounted largely (95%) for the manifestation of white spotting. Last but not the least, ADAMTS20 and TWIST2 may also contribute to the diversification of this coat-color pattern.

CONCLUSIONS

The alternative mutations identified in this study affect, at least partially and independently, the development of melanocytes. The presence and persistence of such mutations may be explained by significant financial and social value of spotted buffaloes used in historical Rambu Solo ceremony in Tana Toraja, Indonesia. Several de novo spontaneous mutations have therefore been favored by traditional breeding for the spotted buffaloes.

Single nucleotide polymorphisms for parentage testing of horse breeds in Korea

OBJECTIVE

In this study, we aimed to evaluate the usability single nucleotide polymorphisms (SNPs) for parentage testing of horse breeds in Korea.

METHODS

The genotypes of 93 horse samples (38 Thoroughbred horses, 17 Jeju horses, 20 Quarter horses, and 18 American miniature horses) were determined using 15 microsatellite (Ms) markers (AHT4, AHT5, ASB2, ASB17, ASB23, CA425, HMS1, HMS2, HMS3, HMS6, HMS7, HTG4, HTG10, LEX3, and VHL20) and 101 SNP markers.

RESULTS

Paternity tests were performed using 15 Ms markers and 101 SNP markers in Thoroughbred horses and Quarter horses. AHT5, ASB2, ASB17, ASB23, CA425, HMS7, HTG10, and LEX3 did not follow Mendelian inheritance in Thoroughbred horses, whereas in Quarter horses, only AHT4, ASB2, and HMS2 showed Mendelian inheritance, consequently, paternity was not established. Meanwhile, 31 markers, including MNEc_2_2_ 2_98568918_BIEC2_502451, in Thoroughbred horses, and 30 markers, including MNEc_ 2_30_7430735_BIEC2_816793, in Quarter horses did not conform with Mendelian inheritance and therefore, could not be used for establishing parentage.

CONCLUSION

The possibility of replacing Ms markers with SNP markers for paternity testing in horses was confirmed. However, further research using more samples is necessary.

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.

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.

Digital Phenotyping Reveals Phenotype Diversity and Epistasis among White Spotting Alleles in the American Paint Horse

White spotting is an iconic feature of the American Paint Horse. The American Paint Horse Association (APHA) is dedicated to recording pedigree and performance of this stock-type breed, while preserving its distinctive coat color and conformation. Here, the depigmented proportion of the coat (% white coat) was measured using digital photograph analysis of 1195 registered American Paint Horses. Genotypes for nine white-spotting polymorphisms commonly found in Paint Horses, and two pigment-producing loci MCIR and ASIP genes, were also provided by the APHA. White-coat percent significantly increased in horses with more white-spotting alleles present, regardless of the number of loci bearing those alleles, likely due to a strong additive genetic effect at each white-spotting locus, as well as an additive epistatic effect among white spotting loci. Paint Horses with a chestnut base coat color (genotype e/e at MC1R) possessed a significantly higher white coat percentage, suggesting confirming an epistatic interaction between pigmentation signaling genes and loci for white spotting. The APHA registry categories of Regular versus Solid Paint-Bred also differed in their median white coat percentage (p < 0.0001), but not in the overall ranges of this phenotype, reenforcing the importance of the regional patterns of the depigmentation in the definition of the desired APHA phenotype. Multi-locus phenotype prediction models for white-coat percentage performed only moderately well, and improvements in the sample size and the number of loci genotyped will likely be needed before such an approach could be used practically by APHA breeders. In the future, models that enable phenotype prediction based on genotypes, and automated phenotype assessment could increase the production of valuable visual traits in the American Paint Horse population and improve the APHA member experience during the registration process.

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.

Genetic characterization of Macaca arctoides: A highlight of key genes and pathways

When compared to the approximately 22 other macaque species, Macaca arctoides has many unique phenotypes. These traits fall into various phenotypic categories, including genitalia, coloration, mating, and olfactory traits. Here we used a previously identified whole genome set of 690 outlier genes to look for possible genetic explanations of these unique traits. Of these, 279 genes were annotated miRNAs, which are non-coding. Patterns within the remaining outliers in coding genes were investigated using GO (n = 370) and String (n = 383) analysis, which showed many interconnected immune-related genes. Further, we compared the outliers to candidate pathways associated with M. arcotides' unique phenotypes, revealing 10/690 outlier genes that overlapped these four pathways: hedgehog signaling, WNT signaling, olfactory, and melanogenesis. Of these, genes in all pathways except olfactory had higher F_ST values than the rest of the genes in the genome based on permutation tests. Overall, our results point to many genes each having a small impact on phenotype, working in tandem to cause large systemic changes. Additionally, these results may indicate pleiotropy. This seems to be especially true with the development and coloration of M. arctoides. Our results highlight that development, melanogenesis, immune function, and miRNAs may be heavily involved in M. arctoides' evolutionary history.

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.

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.

Genomic scan revealed KIT gene underlying white/gray plumage color in Chinese domestic geese

Goose is an important type of domesticated poultry. The wild geese that are regarded as the ancestors of modern domestic geese present gray plumage. Domesticated, geese have both white and gray feathers. To elucidate the genetic mechanisms underlying the formation of white and gray plumage in geese, we resequenced the whole genome of 18 geese from six populations including white and gray goose breeds. The average sequencing depth per individual was 9.81× and the average genome coverage was 96.8%. A total of 346 genes were detected in the top 1% of F_ST scores of gray- and white-feathered geese, and a significant F_ST site was located in the intron region within the KIT gene, the 18 bp deletion in KIT having the strongest potential association with white feathers. It has been reported that a number of genes are associated with plumage colors in birds. However, no studies have identified the relationship between KIT and plumage color in birds at present, although the white coat can be attributed to mutations in KIT in some mammals. Our study showed that that KIT is a plausible candidate gene for white/gray plumage color in Chinese domestic geese.

Selection in the Finnhorse, a native all-around horse breed

Selection by breeders modifies the morphology, behaviour and performance of domesticated species. Here, we examined signs of selection in Finnhorse, the only native horse breed in Finland. We first searched divergent genomic regions between Finnhorses and other breeds, as well as between different breeding sections of the Finnhorse with data from Illumina Equine SNP70 BeadChip, and then studied several of the detected regions in more detail. We found altogether 35 common outlier SNPs between Finnhorses and other breeds using two different selection tests. Many of the SNPs were located close to genes affecting coat colour, performance, size, sugar metabolism, immune response and olfaction. We selected genes affecting coat colour (KIT, MITF, PMEL), performance (MSTN) and locomotion (DMRT3) for a more detailed examination. In addition, we looked for, and found, associations with height at withers and SNPs located close to gene LCORL. Among the four breeding sections of Finnhorses (harness trotters, riding horses, draught horses and pony‐sized horses), a single SNP located close to the DMRT3 gene was significantly differentiated and only between harness trotters and pony‐sized horses.

Genome-wide analysis of homozygosity regions in european simmental bulls

The study of Runs of Homozygosity (ROH) is a useful approach for the characterization of the genome of livestock populations. Due to their high relationship with autozygosity, ROH allow to make inference about population genetic history, to estimate the level of inbreeding, to assess within breed heterogeneity and to detect the footprints of selection on livestock genomes. Aim of this study was to investigate the distribution of runs of homozygosity in bulls belonging to five European Simmental populations and to assess the relationship between three production traits (milk yield, fat and protein contents) and autozygosity. ROH count, distribution and ROH-based coefficient of inbreeding (F_ROH ) were calculated for 3,845 Simmental bulls of five different European countries: Austria (AT), Switzerland (CH), Czech Republic (CZ), Germany (DE) and Italy (IT). Average values of ROH number per animal, and total genome length covered by ROH were 77.8 ± 20.7 and 205 ± 74.4 Mb, respectively. Bulls from AT, DE and IT exhibited similar ROH characteristics. Swiss animals showed the highest (12.6%), while CZ the lowest (4.6%) F_ROH coefficient. The relationship between ROH occurrence and milk production traits was investigated through a genome-wide ROH-traits association analysis (GWRA). A total of 34 regions previously associated with milk traits (yield and/or composition) were identified by GWRA. Results of the present research highlight a mixed genetic background in the 5 European Simmental populations, with the possible presence of three subgroups. Moreover, a strong relationship between autozygosity and production traits has been detected.

Coat Color Roan Shows Association with KIT Variants and No Evidence of Lethality in Icelandic Horses

Roan (Rn) horses show a typical seasonal change of color. Their body is covered with colored and white hair. We performed a descriptive statistical analysis of breeding records of Icelandic horses to challenge the hypothesis of roan being lethal in utero under homozygous condition. The roan to non-roan ratio of foals from roan × roan matings revealed homozygous roan Icelandic horses to be viable. Even though roan is known to be inherited in a dominant mode and epistatic to other coat colors, the causative mutation is still unknown. Nevertheless, an association between roan phenotype and the KIT gene was shown for different horse breeds. In the present study, we identified KIT variants by Sanger sequencing, and show that KIT is also associated with roan in the Icelandic horse breed.

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.

Discovery of selection-driven genetic differences of Duroc, Landrace, and Yorkshire pig breeds by EigenGWAS and Fst analyses

Pigs are one of the earliest domesticated animals and multiple breeds have been developed to meet the various demands of consumers. EigenGWAS is a novel strategy to identify candidate genes that underlying population genetic differences and to infer candidate regions under selection as well. In this study, EigenGWAS and F_st analyses were performed using the public re-sequencing data of three typical commercial pig breeds, Duroc, Landrace and Yorkshire. The intersection of genome-wide significant SNPs detected by EigenGWAS and top-ranked 1% SNPs of F_st results were treated as signals under selection. Using the data of all three breeds, 3062 signals under selection were detected and the nearby genomic regions within 300 kb upstream and downstream covered 6.54% of whole genome. Pairs of breeds were analysed along with the pathway analysis. The gene function enrichment results indicated that many candidate genes located in the genomic regions of the signals under selection were associated with biological processes related to growth, metabolism, reproduction, sensory perception, etc. Among the candidate genes, the FSHB, AHR, PTHLH, KDR and FST genes were reported to be associated with reproductive performance; the KIT, KITLG, MITF, MC1R and EDNRB genes were previously identified to affect coat colour; the RETREG1, TXNIP, BMP5, PPARD and RBP4 genes were reported to be associated with lipid metabolism and growth traits. The identified genetic differences across the three commercial breeds will advance understanding of the artificial selection history of pigs and the signals under selection will suggest potential uses in pig genomic breeding programmes.

Biological adaptations in the Arctic cervid, the reindeer (Rangifer tarandus)

The reindeer is an Arctic species that exhibits distinctive biological characteristics, for which the underlying genetic basis remains largely unknown. We compared the genomes of reindeer against those of other ruminants and nonruminant mammals to reveal the genetic basis of light arrhythmicity, high vitamin D metabolic efficiency, the antler growth trait of females, and docility. We validate that two reindeer vitamin D metabolic genes (CYP27B1 and POR) show signs of positive selection and exhibit higher catalytic activity than those of other ruminants. A mutation upstream of the reindeer CCND1 gene endows an extra functional binding motif of the androgen receptor and thereby may result in female antlers. Furthermore, a mutation (proline-1172→threonine) in reindeer PER2 results in loss of binding ability with CRY1, which may explain circadian arrhythmicity in reindeer.

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 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.

Selection signatures in goats reveal copy number variants underlying breed-defining coat color phenotypes

Domestication and human selection have formed diverse goat breeds with characteristic phenotypes. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. Using whole genome sequencing of DNA pools (pool-seq) from 20 genetically diverse modern goat breeds and bezoars, we identified 2,239 putative selection signatures. In two Pakistani goat breeds, Pak Angora and Barbari, we found selection signatures in a region harboring KIT, a gene involved in melanoblast development, migration, and survival. The search for candidate causative variants responsible for these selective sweeps revealed two different copy number variants (CNVs) downstream of KIT that were exclusively present in white Pak Angora and white-spotted Barbari goats. Several Swiss goat breeds selected for specific coat colors showed selection signatures at the ASIP locus encoding the agouti signaling protein. Analysis of these selective sweeps revealed four different CNVs associated with the white or tan (A^Wt), Swiss markings (A^sm), badgerface (A^b), and the newly proposed peacock (A^pc) allele. RNA-seq analyses on skin samples from goats with the different CNV alleles suggest that the identified structural variants lead to an altered expression of ASIP between eumelanistic and pheomelanistic body areas. Our study yields novel insights into the genetic control of pigmentation by identifying six functionally relevant CNVs. It illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats.

Domestic animals have been selected for hundreds or sometimes even thousands of years for traits that were appreciated by their human owners. This process correlated with the fixation of causative genetic variants controlling breed-specific traits within regions of reduced genetic diversity, so called selection signatures or selective sweeps. We conducted a comprehensive screen for selection signatures in 20 phenotypically and genetically diverse modern goat breeds and identified a total of 2,239 putative selection signatures in our dataset. Follow-up experiments on selection signatures harboring known candidate genes for coat color revealed six different copy number variants (CNVs). Two of these CNVs were located in the 3’-flanking region of KIT and associated with a completely white coat color phenotype in Pak Angora goats and a white-spotted coat color phenotype in Barbari goats, respectively. The other four CNVs were located at the ASIP locus. They were associated with four different types of coat color patterning in seven Swiss goat breeds. Their functional effect is mediated by region-specific quantitative changes in ASIP mRNA expression. Our study illustrates how structural changes of the genome have contributed to phenotypic evolution in domestic goats.

Genome-wide association analysis reveals QTL and candidate mutations involved in white spotting in cattle

BACKGROUND

White spotting of the coat is a characteristic trait of various domestic species including cattle and other mammals. It is a hallmark of Holstein-Friesian cattle, and several previous studies have detected genetic loci with major effects for white spotting in animals with Holstein-Friesian ancestry. Here, our aim was to better understand the underlying genetic and molecular mechanisms of white spotting, by conducting the largest mapping study for this trait in cattle, to date.

RESULTS

Using imputed whole-genome sequence data, we conducted a genome-wide association analysis in 2973 mixed-breed cows and bulls. Highly significant quantitative trait loci (QTL) were found on chromosomes 6 and 22, highlighting the well-established coat color genes KIT and MITF as likely responsible for these effects. These results are in broad agreement with previous studies, although we also report a third significant QTL on chromosome 2 that appears to be novel. This signal maps immediately adjacent to the PAX3 gene, which encodes a known transcription factor that controls MITF expression and is the causal locus for white spotting in horses. More detailed examination of these loci revealed a candidate causal mutation in PAX3 (p.Thr424Met), and another candidate mutation (rs209784468) within a conserved element in intron 2 of MITF transcripts expressed in the skin. These analyses also revealed a mechanistic ambiguity at the chromosome 6 locus, where highly dispersed association signals suggested multiple or multiallelic QTL involving KIT and/or other genes in this region.

CONCLUSIONS

Our findings extend those of previous studies that reported KIT as a likely causal gene for white spotting, and report novel associations between candidate causal mutations in both the MITF and PAX3 genes. The sizes of the effects of these QTL are substantial, and could be used to select animals with darker, or conversely whiter, coats depending on the desired characteristics.

The Genetic Basis of Piebald Coat Colour in Hucul Horses in the Context of White Markings and Crypto-Tobiano as a Breeding Problem in Poland

The aim of the study was to analyse the genetic basis of piebald coat colour in Hucul horses and to verify their coat colour in breeding records. Tests were performed with DNA purified from the whole blood samples of 242 Hucul horses with different coat colour patterns. DNA was analysed to identify an inversion in ECA3 (PCR). The results confirmed that the inversion on ECA3 is a direct factor determining piebald (tobiano) colour in the analysed Hucul horses. No inversion was observed in any of the solid coloured horses, but it was present in all the piebald ones. It was also identified in 18% (11 of 61) of the horses from the group of horses qualified in the passport as solid coloured with white markings. In fact, these horses had the tobiano gene that is phenotypically identifiable as crypto-tobiano, which may give the false impression of having white markings and lead to error when describing a horse. This is an important issue, in particular with regard to the breed standard, which eliminates Hucul horses with white markings from breeding.

High-density genotyping reveals signatures of selection related to acclimation and economically important traits in 15 local sheep breeds from Russia

Background

Domestication and centuries of selective breeding have changed genomes of sheep breeds to respond to environmental challenges and human needs. The genomes of local breeds, therefore, are valuable sources of genomic variants to be used to understand mechanisms of response to adaptation and artificial selection. As a step toward this we performed a high-density genotyping and comprehensive scans for signatures of selection in the genomes from 15 local sheep breeds reared across Russia.

Results

Results demonstrated that the genomes of Russian sheep breeds contain multiple regions under putative selection. More than 50% of these regions matched with intervals identified in previous scans for selective sweeps in sheep genomes. These regions contain well-known candidate genes related to morphology, adaptation, and domestication (e.g., KITLG, KIT, MITF, and MC1R), wool quality and quantity (e.g., DSG@, DSC@, and KRT@), growth and feed intake (e.g., HOXA@, HOXC@, LCORL, NCAPG, LAP3, and CCSER1), reproduction (e.g., CMTM6, HTRA1, GNAQ, UBQLN1, and IFT88), and milk-related traits (e.g., ABCG2, SPP1, ACSS1, and ACSS2). In addition, multiple genes that are putatively related to environmental adaptations were top-ranked in selected intervals (e.g., EGFR, HSPH1, NMUR1, EDNRB, PRL, TSHR, and ADAMTS5). Moreover, we observed that multiple key genes involved in human hereditary sensory and autonomic neuropathies, and genetic disorders accompanied with an inability to feel pain and environmental temperatures, were top-ranked in multiple or individual sheep breeds from Russia pointing to a possible mechanism of adaptation to harsh climatic conditions.

Conclusions

Our work represents the first comprehensive scan for signatures of selection in genomes of local sheep breeds from the Russian Federation of both European and Asian origins. We confirmed that the genomes of Russian sheep contain previously identified signatures of selection, demonstrating the robustness of our integrative approach. Multiple novel signatures of selection were found near genes which could be related to adaptation to the harsh environments of Russia. Our study forms a basis for future work on using Russian sheep genomes to spot specific genetic variants or haplotypes to be used in efforts on developing next-generation highly productive breeds, better suited to diverse Eurasian environments.

Electronic supplementary material

The online version of this article (10.1186/s12864-019-5537-0) contains supplementary material, which is available to authorized users.

Selection signatures in four German warmblood horse breeds: Tracing breeding history in the modern sport horse

The study of selection signatures helps to find genomic regions that have been under selective pressure and might host genes or variants that modulate important phenotypes. Such knowledge improves our understanding of how breeding programmes have shaped the genomes of livestock. In this study, 942 stallions were included from four, exemplarily chosen, German warmblood breeds with divergent historical and recent selection focus and different crossbreeding policies: Trakehner (N = 44), Holsteiner (N = 358), Hanoverian (N = 319) and Oldenburger (N = 221). Those breeds are nowadays bred for athletic performance and aptitude for show-jumping, dressage or eventing, with a particular focus of Holsteiner on the first discipline. Blood samples were collected during the health exams of the stallion preselections before licensing and were genotyped with the Illumina EquineSNP50 BeadChip. Autosomal markers were used for a multi-method search for signals of positive selection. Analyses within and across breeds were conducted by using the integrated Haplotype Score (iHS), cross-population Extended Haplotype Homozygosity (xpEHH) and Runs of Homozygosity (ROH). Oldenburger and Hanoverian showed very similar iHS signatures, but breed specificities were detected on multiple chromosomes with the xpEHH. The Trakehner clustered as a distinct group in a principal component analysis and also showed the highest number of ROHs, which reflects their historical bottleneck. Beside breed specific differences, we found shared selection signals in an across breed iHS analysis on chromosomes 1, 4 and 7. After investigation of these iHS signals and shared ROH for potential functional candidate genes and affected pathways including enrichment analyses, we suggest that genes affecting muscle functionality (TPM1, TMOD2-3, MYO5A, MYO5C), energy metabolism and growth (AEBP1, RALGAPA2, IGFBP1, IGFBP3-4), embryonic development (HOXB-complex) and fertility (THEGL, ZPBP1-2, TEX14, ZP1, SUN3 and CFAP61) have been targeted by selection in all breeds. Our findings also indicate selection pressure on KITLG, which is well-documented for influencing pigmentation.

Comparative transcriptome analysis reveals the genetic basis of coat color variation in Pashmina goat

The genetics of coat color variation remains a classic area. Earlier studies have focused on a limited number of genes involved in color determination; however, the complete set of trait determinants are still not well known. In this study, we used high-throughput sequencing technology to identify and characterize intricate interactions between genes that cause complex coat color variation in Changthangi Pashmina goats, producer of finest and costly commercial animal fiber. We systematically identified differentially expressed mRNAs and lncRNAs from black, brown and white Pashmina goat skin samples by using RNA-sequencing technique. A pairwise comparison of black, white and brown skin samples yielded 2479 significantly dysregulated genes (2422 mRNA and 57 lncRNAs). Differentially expressed genes were enriched in melanin biosynthesis, melanocyte differentiation, developmental pigmentation, melanosome transport activities GO terms. Our analysis suggested the potential role of lncRNAs on color coding mRNAs in cis and trans configuration. We have also developed online data repository as a component of the study to provide a central location for data access, visualization and interpretation accessible through http://pcd.skuastk.org/.

Whole-genome sequencing reveals a large deletion in the MITF gene in horses with white spotted coat colour and increased risk of deafness

White spotting phenotypes in horses are highly valued in some breeds. They are quite variable and may range from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for white spotting phenotypes in horses. For the present study, we investigated an American Paint Horse family segregating a phenotype involving white spotting and blue eyes. Six of eight horses with the white-spotting phenotype were deaf. We obtained whole-genome sequence data from an affected horse and specifically searched for structural variants in the known candidate genes. This analysis revealed a heterozygous ~63-kb deletion spanning exons 6-9 of the MITF gene (chr16:21 503 211-21 566 617). We confirmed the breakpoints of the deletion by PCR and Sanger sequencing. PCR-based genotyping revealed that all eight available affected horses from the family carried the deletion. The finding of an MITF variant fits well with the syndromic phenotype involving both depigmentation and an increased risk for deafness and corresponds to human Waardenburg syndrome type 2A. Our findings will enable more precise genetic testing for depigmentation phenotypes in horses.

Scans for signatures of selection in Russian cattle breed genomes reveal new candidate genes for environmental adaptation and acclimation

Domestication and selective breeding has resulted in over 1000 extant cattle breeds. Many of these breeds do not excel in important traits but are adapted to local environments. These adaptations are a valuable source of genetic material for efforts to improve commercial breeds. As a step toward this goal we identified candidate regions to be under selection in genomes of nine Russian native cattle breeds adapted to survive in harsh climates. After comparing our data to other breeds of European and Asian origins we found known and novel candidate genes that could potentially be related to domestication, economically important traits and environmental adaptations in cattle. The Russian cattle breed genomes contained regions under putative selection with genes that may be related to adaptations to harsh environments (e.g., AQP5, RAD50, and RETREG1). We found genomic signatures of selective sweeps near key genes related to economically important traits, such as the milk production (e.g., DGAT1, ABCG2), growth (e.g., XKR4), and reproduction (e.g., CSF2). Our data point to candidate genes which should be included in future studies attempting to identify genes to improve the extant breeds and facilitate generation of commercial breeds that fit better into the environments of Russia and other countries with similar climates.

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.

Spotted phenotypes in horses lost attractiveness in the Middle Ages

Horses have been valued for their diversity of coat colour since prehistoric times; this is especially the case since their domestication in the Caspian steppe in ~3,500 BC. Although we can assume that human preferences were not constant, we have only anecdotal information about how domestic horses were influenced by humans. Our results from genotype analyses show a significant increase in spotted coats in early domestic horses (Copper Age to Iron Age). In contrast, medieval horses carried significantly fewer alleles for these phenotypes, whereas solid phenotypes (i.e., chestnut) became dominant. This shift may have been supported because of (i) pleiotropic disadvantages, (ii) a reduced need to separate domestic horses from their wild counterparts, (iii) a lower religious prestige, or (iv) novel developments in weaponry. These scenarios may have acted alone or in combination. However, the dominance of chestnut is a remarkable feature of the medieval horse population.

Distribution of coat-color-associated alleles in the domestic horse population and Przewalski's horse

Considering the hidden mode of inheritance of some coat-color-associated alleles, we investigated the presence/absence of coat-color-associated alleles in 1093 domestic horses of 55 breeds and 20 specimens of Przewalski's horse. For coat-color genotyping, allele specific PCR, pyrosequencing and Li-Cor analyses were conducted on 12 coat-color-associated alleles of five genes. Our data provide deep insight into the distribution of coat-color-associated alleles within breeds. We found that the alleles for the basic colorations (bay, black, and chestnut) are widely distributed and occur in nearly all breeds. Alleles leading to dilutions or patterns are rare in domestic breeds and were not found in Przewalski's horse. Higher frequencies of these alleles are only found in breeds that are selected for their expressed phenotypes (e.g., Kinsky horse, Lewitzer, Tinker). Nevertheless, our study produced strong evidence that molecular testing of the coat color is necessary for well-defined phenotyping to avoid unexpected colorations of offspring that can result in legal action.

The Genetic Architecture of Domestication in Animals

Domestication has been essential to the progress of human civilization, and the process itself has fascinated biologists for hundreds of years. Domestication has led to a series of remarkable changes in a variety of plants and animals, in what is termed the “domestication phenotype.” In domesticated animals, this general phenotype typically consists of similar changes in tameness, behavior, size/morphology, color, brain composition, and adrenal gland size. This domestication phenotype is seen in a range of different animals. However, the genetic basis of these associated changes is still puzzling. The genes for these different traits tend to be grouped together in clusters in the genome, though it is still not clear whether these clusters represent pleiotropic effects, or are in fact linked clusters. This review focuses on what is currently known about the genetic architecture of domesticated animal species, if genes of large effect (often referred to as major genes) are prevalent in driving the domestication phenotype, and whether pleiotropy can explain the loci underpinning these diverse traits being colocated.

The Genetics of Deafness in Domestic Animals

Although deafness can be acquired throughout an animal's life from a variety of causes, hereditary deafness, especially congenital hereditary deafness, is a significant problem in several species. Extensive reviews exist of the genetics of deafness in humans and mice, but not for deafness in domestic animals. Hereditary deafness in many species and breeds is associated with loci for white pigmentation, where the cochlear pathology is cochleo-saccular. In other cases, there is no pigmentation association and the cochlear pathology is neuroepithelial. Late onset hereditary deafness has recently been identified in dogs and may be present but not yet recognized in other species. Few genes responsible for deafness have been identified in animals, but progress has been made for identifying genes responsible for the associated pigmentation phenotypes. Across species, the genes identified with deafness or white pigmentation patterns include MITF, PMEL, KIT, EDNRB, CDH23, TYR, and TRPM1 in dog, cat, horse, cow, pig, sheep, ferret, mink, camelid, and rabbit. Multiple causative genes are present in some species. Significant work remains in many cases to identify specific chromosomal deafness genes so that DNA testing can be used to identify carriers of the mutated genes and thereby reduce deafness prevalence.