Novel variants in the KIT and PAX3 genes in horses with white‐spotted coat colour phenotypes

Variants in the EDNRB, KIT, MITF, PAX3 and TRPM1 genes are known to cause white spotting phenotypes in horses, which can range from the common white markings up to completely white horses. In this study, we investigated these candidate genes in 169 horses with white spotting phenotypes not explained by the previously described variants. We identified a novel missense variant, PAX3:p.Pro32Arg, in Appaloosa horses with a splashed white phenotype in addition to their leopard complex spotting patterns. We also found three novel variants in the KIT gene. The splice site variant c.1346+1G>A occurred in a Swiss Warmblood horse with a pronounced depigmentation phenotype. The missense variant p.Tyr441Cys was present in several part‐bred Arabians with sabino‐like depigmentation phenotypes. Finally, we provide evidence suggesting that the common and widely distributed KIT:p.Arg682His variant has a very subtle white‐increasing effect, which is much less pronounced than the effect of the other described KIT variants. We termed the new KIT variants W18–W20 to provide a simple and unambiguous nomenclature for future genetic testing applications.     

A complex structural variant at the KIT locus in cattle with the Pinzgauer spotting pattern

A specific white spotting phenotype, termed finching or line‐backed spotting, is known for all Pinzgauer cattle and occurs occasionally in Tux‐Zillertaler cattle, two Austrian breeds. The so‐called Pinzgauer spotting is inherited as an autosomal incompletely dominant trait. A genome‐wide association study using 27 white spotted and 16 solid‐coloured Tux‐Zillertaler cattle, based on 777k SNP data, revealed a strong signal on chromosome 6 at the KIT locus. Haplotype analyses defined a critical interval of 122 kb downstream of the KIT coding region. Whole‐genome sequencing of a Pinzgauer cattle and comparison to 338 control genomes revealed a complex structural variant consisting of a 9.4‐kb deletion and an inversely inserted duplication of 1.5 kb fused to a 310‐kb duplicated segment from chromosome 4. A diagnostic PCR was developed for straightforward genotyping of carriers for this structural variant (KIT^PINZ) and confirmed that the variant allele was present in all Pinzgauer and most of the white spotted Tux‐Zillertaler cattle. In addition, we detected the variant in all Slovenian Cika, British Gloucester and Spanish Berrenda en negro cattle with similar spotting patterns. Interestingly, the KIT^PINZ variant occurs in some white spotted animals of the Swiss breeds Evolèner and Eringer. The introgression of the KIT^PINZ variant confirms admixture and the reported historical relationship of these short‐headed breeds with Austrian Tux‐Zillertaler and suggests a mutation event, occurring before breed formation.     

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

White spotting phenotypes in horses can range in severity from the common white markings up to completely white horses. EDNRB, KIT, MITF, PAX3 and TRPM1 represent known candidate genes for such phenotypes in horses. For the present study, we re‐investigated a large horse family segregating a variable white spotting phenotype, for which conventional Sanger sequencing of the candidate genes’ individual exons had failed to reveal the causative variant. 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 ~1.9‐kb deletion spanning exons 10–13 of the KIT gene (chr3:77,740,239_77,742,136del1898insTATAT). In continuity with previously named equine KIT variants we propose to designate the newly identified deletion variant W22. We had access to 21 horses carrying the W22 allele. Four of them were compound heterozygous W20/W22 and had a completely white phenotype. Our data suggest that W22 represents a true null allele of the KIT gene, whereas the previously identified W20 leads to a partial loss of function. These findings will enable more precise genetic testing for depigmentation phenotypes in horses.     

Two variants in the KIT gene as candidate causative mutations for a dominant white and a white spotting phenotype in the donkey

White spotting phenotypes have been intensively studied in horses, and although similar phenotypes occur in the donkey, little is known about the molecular genetics underlying these patterns in donkeys. White spotting in donkeys can range from only a few white areas to almost complete depigmentation and is characterised by a loss of pigmentation usually progressing from a white spot in the hip area. Completely white‐born donkeys are rare, and the phenotype is characterised by the complete absence of pigment resulting in pink skin and a white coat. A dominant mode of inheritance has been demonstrated for spotting in donkeys. Although the mode of inheritance for the completely white phenotype in donkeys is not clear, the phenotype shows similarities to dominant white in horses. As variants in the KIT gene are known to cause a range of white phenotypes in the horse, we investigated the KIT gene as a potential candidate gene for two phenotypes in the donkey, white spotting and white. A mutation analysis of all 21 KIT exons identified a missense variant in exon 4 (c.662A>C; p.Tyr221Ser) present only in a white‐born donkey. A second variant affecting a splice donor site (c.1978+2T>A) was found exclusively in donkeys with white spotting. Both variants were absent in 24 solid‐coloured controls. To the authors’ knowledge, this is the first study investigating genetic mechanisms underlying white phenotypes in donkeys. Our results suggest that two independent KIT alleles are probably responsible for white spotting and white in donkeys.     

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.     

Genome‐wide variability and selection signatures in Italian island cattle breeds

In the present study, a sample of 88 animals belonging to four local (Modicana, Sarda, Sardo‐Bruna and Sardo‐Modicana) and one cosmopolitan (Italian Brown Swiss) cattle breeds were genotyped with a medium density SNP beadchip and compared to investigate their genetic diversity and the existence of selection signatures. A total of 43 012 SNPs distributed across all 29 autosomal chromosomes were retained after data quality control. Basic population statistics, Wright fixation index and runs of homozygosity (ROH) analyses confirmed that the Italian Brown Swiss genome was shaped mainly by selection, as underlined by the low values of heterozygosity and minor allele frequency. As expected, local cattle exhibited a large within‐breed genetic heterogeneity. The F_ST comparison revealing the largest number of significant SNPs was Sardo‐Bruna vs. Sardo‐Modicana, whereas the smallest was observed for Italian Brown Swiss vs. Sardo‐Modicana. Modicana exhibited the largest number of detected ROHs, whereas the smallest was observed for Sardo‐Modicana. Signatures of selection were detected in genomic regions that harbor genes involved in milk production traits for Italian Brown Swiss and fitness traits for local breeds. According to the results of multi‐dimensional scaling and the admixture analysis the Sardo‐Bruna is more similar to the Sarda than to the Italian Brown Swiss breed. Moreover, the Sardo‐Modicana is genetically closer to the Modicana than to the Sarda breed. Results of the present work confirm the usefulness of single nucleotide polymorphisms in deciphering the genetic architecture of livestock breeds.     

The genetics of brown coat color and white spotting in domestic yaks (Bos grunniens)

Domestic yaks (Bos grunniens) exhibit two major coat color variations: a brown vs. wild‐type black pigmentation and a white spotting vs. wild‐type solid color pattern. The genetic basis for these variations in color and distribution remains largely unknown and may be complicated by a breeding history involving hybridization between yaks and cattle. Here, we investigated 92 domestic yaks from China using a candidate gene approach. Sequence variations in MC1R, PMEL and TYRP1 were surveyed in brown yaks; TYRP1 was unassociated with the coloration and excluded. Recessive mutations from MC1R, or p.Gln34*, p.Met73Leu and possibly p.Arg142Pro, are reported in bovids for the first time and accounted for approximately 40% of the brown yaks in this study. The remaining 60% of brown individuals correlated with a cattle‐derived deletion mutation from PMEL (p.Leu18del) in a dominant manner. Degrees of white spotting found in yaks vary from color sidedness and white face, to completely white. After examining the candidate gene KIT, we suggest that color‐sided and all‐white yaks are caused by the serial translations of KIT (Cs₆ or Cs₂₉) as reported for cattle. The white‐faced phenotype in yaks is associated with the KIT haplotype S^wf. All KIT mutations underlying the serial phenotypes of white spotting in yaks are identical to those in cattle, indicating that cattle are the likely source of white spotting in yaks. Our results reveal the complex genetic origins of domestic yak coat color as either native in yaks through evolution and domestication or as introduced from cattle through interspecific hybridization.     

Allelic heterogeneity at the equine KIT locus in dominant white (W) horses

White coat color has been a highly valued trait in horses for at least 2,000 years. Dominant white (W) is one of several known depigmentation phenotypes in horses. It shows considerable phenotypic variation, ranging from ~50% depigmented areas up to a completely white coat. In the horse, the four depigmentation phenotypes roan, sabino, tobiano, and dominant white were independently mapped to a chromosomal region on ECA 3 harboring the KIT gene. KIT plays an important role in melanoblast survival during embryonic development. We determined the sequence and genomic organization of the ~82 kb equine KIT gene. A mutation analysis of all 21 KIT exons in white Franches-Montagnes Horses revealed a nonsense mutation in exon 15 (c.2151C>G, p.Y717X). We analyzed the KIT exons in horses characterized as dominant white from other populations and found three additional candidate causative mutations. Three almost completely white Arabians carried a different nonsense mutation in exon 4 (c.706A>T, p.K236X). Six Camarillo White Horses had a missense mutation in exon 12 (c.1805C>T, p.A602V), and five white Thoroughbreds had yet another missense mutation in exon 13 (c.1960G>A, p.G654R). Our results indicate that the dominant white color in Franches-Montagnes Horses is caused by a nonsense mutation in the KIT gene and that multiple independent mutations within this gene appear to be responsible for dominant white in several other modern horse populations.     

Genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle of Canada based on 15 bovine microsatellite markers.

The genetic diversity among Canadienne, Brown Swiss, Holstein, and Jersey cattle was estimated from relationships determined by genotyping 20 distantly related animals in each breed for 15 microsatellites located on separate chromosomes. The Canadienne, Holstein, and Jersey cattle had an average of six alleles per loci compared with five alleles for Brown Swiss. Furthermore, a number of potentially breed-specific alleles were identified. The allele size variance among breeds was similar, but varied considerably among loci. All of the loci studied were equally heterozygous, as were Brown Swiss, Canadienne, and Holstein cattle (0.68-0.69) whereas Jersey cattle showed lower heterozygosity (0.59). The within-breed estimates of genetic distance were greater than zero and significant. The genetic distance between Canadienne and Holstein (0.156), Brown Swiss (0.243), and Jersey (0.235) was negligible, suggesting close relationship. Concurrently, Brown Swiss and Holstein (0.211) cattle also demonstrated close relationship. In contrast, the Jersey breed was genetically distant from the Brown Swiss and Holstein cattle (0.427 and 0.320, respectively). The characterization of Canadienne cattle, as part of the genetic resource conservation effort currently underway in Canada, underscores the difficulty in scientifically establishing unique breeds. Therefore, the need to consider all relevant morphological characteristics and production performance in combination with available cultural, historical, pedigree, and molecular information becomes relevant when identifying breeds for conservation.     

Colour-sidedness in Gloucester cattle is associated with a complex structural variant impacting regulatory elements downstream of KIT

Colour-sidedness is a striking coat colour pattern found in a number of cattle breeds, typically characterised by a white stripe that extends along the back, head and belly of the animal. This dominant phenotype is caused by two related translocations (Cs₆ and Cs₂₉) that alter a region downstream of the KIT gene. Gloucester cattle are native to the UK and are known for an unusual colour-sided pattern that does not extend to the head. We carried out whole-genome sequencing of two Gloucester bulls as well as colour-sided Irish Moiled, British White and Pustertaler Sprinzen for comparison. We found that the Gloucester cattle also have a complex structural variant downstream of KIT, which overlaps the regions involved in Cs₆ and Cs_29. All three alleles potentially disrupt a number of putative regulatory elements downstream of KIT. These results complement and expand on the recently published work focused on the Pinzgauer breed from Austria, a carrier of the same colour-sided pattern as seen in Gloucester cattle.     

A Function Accounting for Training Set Size and Marker Density to Model the Average Accuracy of Genomic Prediction

Prediction of genomic breeding values is of major practical relevance in dairy cattle breeding. Deterministic equations have been suggested to predict the accuracy of genomic breeding values in a given design which are based on training set size, reliability of phenotypes, and the number of independent chromosome segments (). The aim of our study was to find a general deterministic equation for the average accuracy of genomic breeding values that also accounts for marker density and can be fitted empirically. Two data sets of 5′698 Holstein Friesian bulls genotyped with 50 K SNPs and 1′332 Brown Swiss bulls genotyped with 50 K SNPs and imputed to ∼600 K SNPs were available. Different k-fold (k = 2–10, 15, 20) cross-validation scenarios (50 replicates, random assignment) were performed using a genomic BLUP approach. A maximum likelihood approach was used to estimate the parameters of different prediction equations. The highest likelihood was obtained when using a modified form of the deterministic equation of Daetwyler et al. (2010), augmented by a weighting factor (w) based on the assumption that the maximum achievable accuracy is . The proportion of genetic variance captured by the complete SNP sets () was 0.76 to 0.82 for Holstein Friesian and 0.72 to 0.75 for Brown Swiss. When modifying the number of SNPs, w was found to be proportional to the log of the marker density up to a limit which is population and trait specific and was found to be reached with ∼20′000 SNPs in the Brown Swiss population studied.     

Molecular genetics of coat colour variations in White Galloway and White Park cattle

White Galloway cattle exhibit three different white coat colour phenotypes, that is, well marked, strongly marked and mismarked. However, mating of individuals with the preferred well or strongly marked phenotype also results in offspring with the undesired mismarked and/or even fully black coat colour. To elucidate the genetic background of the coat colour variations in White Galloway cattle, we analysed four coat colour relevant genes: mast/stem cell growth factor receptor (KIT), KIT ligand (KITLG), melanocortin 1 receptor (MC1R) and tyrosinase (TYR). Here, we show that the coat colour variations in White Galloway cattle and White Park cattle are caused by a KIT gene (chromosome 6) duplication and aberrant insertion on chromosome 29 (Cs₂₉) as recently described for colour‐sided Belgian Blue. Homozygous (Cs₂₉/Cs₂₉) White Galloway cattle and White Park cattle exhibit the mismarked phenotype, whereas heterozygous (Cs₂₉/wt₂₉) individuals are either well or strongly marked. In contrast, fully black individuals are characterised by the wild‐type chromosome 29. As known for other cattle breeds, mutations in the MC1R gene determine the red colouring. Our data suggest that the white coat colour variations in White Galloway cattle and White Park cattle are caused by a dose‐dependent effect based on the ploidy of aberrant insertions and inheritance of the KIT gene on chromosome 29.     

De novo stop-lost germline mutation in FGFR3 causes severe chondrodysplasia in the progeny of a Holstein bull

Fifteen cases of chondrodysplasia characterized by disproportionate dwarfism occurred in the progeny of a single Holstein bull. A de novo mutation event in the germline of the sire was suspected as cause. Whole-genome sequencing revealed a single protein-changing variant in the stop codon of FGFR3 gene on chromosome 6. Sanger sequencing of EDTA blood proved that this variant occurred de novo and segregates perfectly with the observed phenotype in the affected cattle family. FGFR3 is an important regulator gene in bone formation owing to its key role in the bone elongation induced by FGFR3-dimers. The detected paternally inherited stop-lost variant in FGFR3 is predicted to add 93 additional amino acids to the protein’s C-terminus. This study provides a second example of a dominant FGFR3 stop-lost variant as a pathogenic mutation of a severe form of chondrodysplasia. Even though FGFR3 is known to be associated with dwarfism and growth disorders in human and sheep, this study is the first to describe FGFR3-associated chondrodysplasia in cattle.     

Genome‐wide assessment of diversity and differentiation between original and modern Brown cattle populations

Identifying genomic regions involved in the differences between breeds can provide information on genes that are under the influence of both artificial and natural selection. The aim of this study was to assess the genetic diversity and differentiation among four different Brown cattle populations (two original vs. two modern populations) and to characterize the distribution of runs of homozygosity (ROH) islands using the Illumina Bovine SNP50 BeadChip genotyping data. After quality control, 34 735 SNPs and 106 animals were retained for the analyses. Larger heterogeneity was highlighted for the original populations. Patterns of genetic differentiation, multidimensional scaling, and the neighboring joining tree distinguished the modern from the original populations. The F_ST‐outlier identified several genes putatively involved in the genetic differentiation between the two groups, such as stature and growth, behavior, and adaptability to local environments. The ROH islands within both the original and the modern populations overlapped with QTL associated with relevant traits. In modern Brown (Brown Swiss and Italian Brown), ROH islands harbored candidate genes associated with milk production traits, in evident agreement with the artificial selection conducted to improve this trait in these populations. In original Brown (Original Braunvieh and Braunvieh), we identified candidate genes related with fat deposition, confirming that breeding strategies for the original Brown populations aimed to produce dual‐purpose animals. Our study highlighted the presence of several genomic regions that vary between Brown populations, in line with their different breeding histories.     

A novel KIT deletion variant in a German Riding Pony with white‐spotting coat colour phenotype

White spotting phenotypes in horses may be caused by developmental alterations impairing melanoblast differentiation, survival, migration and/or proliferation. Candidate genes for white‐spotting phenotypes in horses include EDNRB, KIT, MITF, PAX3 and TRPM1. We investigated a German Riding Pony with a sabino‐like phenotype involving extensive white spots on the body together with large white markings on the head and almost completely white legs. We obtained whole genome sequence data from this horse. The analysis revealed a heterozygous 1273‐bp deletion spanning parts of intron 2 and exon 3 of the equine KIT gene (Chr3: 79 579 925–79 581 197). We confirmed the breakpoints of the deletion by PCR and Sanger sequencing. Knowledge of the functional impact of similar KIT variants in horses and other species suggests that this deletion represents a plausible candidate causative variant for the white‐spotting phenotype. We propose the designation W28 for the mutant allele.     

Identification of small and large genomic candidate variants in bovine pulmonary hypoplasia and anasarca syndrome

The pulmonary hypoplasia and anasarca syndrome (PHA) is a congenital lethal disorder, which until now has been reported in cattle and sheep. PHA is characterized by extensive subcutaneous fetal edema combined with hypoplasia or aplasia of the lungs and dysplasia of the lymphatic system. PHA is assumed to be of genetic etiology. This study presents the occurrence of PHA in two different cattle breeds and their genetic causation. Two PHA cases from one sire were observed in Slovenian Cika cattle. Under the assumption of monogenic inheritance, genome-wide homozygosity mapping scaled down the critical regions to 3% of the bovine genome including a 43.6 Mb-sized segment on chromosome 6. Whole-genome sequencing of one case, variant filtering against controls and genotyping of a larger cohort of Cika cattle led to the detection of a likely pathogenic protein-changing variant perfectly associated with the disease: a missense variant on chromosome 6 in ADAMTS3 (NM_001192797.1: c.1222C>T), which affects an evolutionary conserved residue (NP_001179726.1: p.(His408Tyr)). A single PHA case was found in Danish Holstein cattle and was whole-genome sequenced along with its parents. However, as there was no plausible private protein-changing variant, mining for structural variation revealed a likely pathogenic trisomy of the entire chromosome 20. The identified ADAMTS3 associated missense variant and the trisomy 20 are two different genetic causes, which shows a compelling genetic heterogeneity for bovine PHA.     

Arachnomelia in Brown Swiss cattle maps to chromosome 5

Arachnomelia in Brown Swiss cattle is a monogenic autosomal recessive inherited congenital disorder of the skeletal system giving affected calves a spidery look (OMIA ID 000059). Over a period of 20 years 15 cases were sampled in the Swiss and Italian Brown cattle population. Pedigree data revealed that all affected individuals trace back to a single acknowledged carrier founder sire. A genome scan using 240 microsatellites spanning the 29 bovine autosomes showed homozygosity at three adjacent microsatellite markers on bovine Chr 5 in all cases. Linkage analysis confirmed the localization of the arachnomelia mutation in the region of the marker ETH10. Fine-mapping and haplotype analysis using a total of 34 markers in this region refined the critical region of the arachnomelia locus to a 7.19-Mb interval on bovine Chr 5. The disease-associated IBD haplotype was shared by 36 proven carrier animals and allows marker-assisted selection. As the corresponding human and mouse chromosome segments do not contain any clear functional candidate genes for this disorder, the mutation causing arachnomelia in the Brown Swiss cattle might help to identify an unknown gene in bone development. Electronic supplementary material  The online version of this article (doi:) contains supplementary material, which is available to authorized users.     

Association of the expression levels in the skeletal muscle and a SNP in the CDC10 gene with growth‐related traits in Japanese Black beef cattle

Growth performance, as well as marbling, is the main breeding objective in Japanese Black (JB) cattle, the major beef breed in Japan. The septin 7 (CDC10) gene, involved in cellular proliferation, is located within a genomic region of a quantitative trait locus for growth‐related traits. In this study, we first showed that the expression levels of the CDC10 gene in the skeletal muscle were higher in JB steers with extremely high growth performance than in JB steers with extremely low growth, using real‐time PCR. Further, a single nucleotide polymorphism (SNP), NC_007302.5:g.63264949G>C, was detected in the promoter region of the CDC10 gene and genotyped in three Japanese cattle breeds (known as ‘Wagyu’ in Japan) and the Brown Swiss dairy cattle breed. All four cattle populations showed a moderate genetic diversity at the SNP of the CDC10 gene. An association analysis indicated that the SNP was associated with growth‐related traits in JB cattle. These findings suggest possible effects of the expression levels in the skeletal muscle and the SNP of the CDC10 gene on growth‐related traits in JB cattle. The CDC10 SNP may be useful for effective marker‐assisted selection to increase beef productivity in JB beef cattle.     

KIT and melanoma predisposition in pigs: sequence variants and association analysis

KIT mutations have been detected in different cancer subtypes, including melanoma. The gene also has been extensively studied in farm animals for its prominent role in coat color. The present work aimed at detecting KIT variants in a porcine model of cutaneous melanoma, the melanoblastoma‐bearing Libechov Minipig (MeLiM). By sequencing exons and intron borders, 36 SNPs and one indel were identified. Of 10 coding SNPs, three were non‐synonymous mutations, likely to affect the protein conformation. A promising variant, located in exon 19 (p.Val870Ala), was genotyped in a MeLiM × Duroc cross, and an association analysis was conducted on several melanoma‐related traits. This variant showed a significant association with melanoma development, tumor ulceration and cutaneous invasion. In conclusion, although the KIT gene would not be a major causal gene for melanoma development in pig, its genetic variation could be influencing this trait.