Mutations in MITF and PAX3 cause "splashed white" and other white spotting phenotypes in horses

Three Novel KIT Polymorphisms Found in Horses with White Coat Color Phenotypes

This paper reports three novel KIT variants likely responsible for previously unexplained white patterning phenotypes observed in three groups of horses. White spots and markings may have substantial consequences on the value and health of domesticated horses. This study aims to elucidate the genetic mechanisms underlying depigmented coat colors to aid in producing prosperous herds. Aligned whole genome sequences were manually screened to identify three polymorphisms in a family of Anglo-Arabian horses (N = 7), a family of Warmblood horses (N = 5), and a single stock-type mare with unexplained white markings. Sanger sequencing confirmed the presence of the variants, and in silico predictive programs were used to predict the functional impacts of each. We propose to term the novel variants W37, W38, and W39, respectively, per convention. The W37 polymorphism was always observed in the presence of one W35 allele, suggesting complete linkage. All three variants were predicted to alter or remove the KIT protein active domain, repressing typical protein folding and impacting pathways that upregulate pigmentation. The severe predicted impact on biological function suggests that these variants may cause increased white spotting, providing a possible explanation for the depigmentation phenotypes observed in affected individuals.

A Ball Python Colour Morph Implicates MC1R in Melanophore-Xanthophore Distribution and Pattern Formation

Reptiles showcase an extensive array of skin colours and patterns, yet little is known about the genetics of reptile colouration. Here, we investigate the genetic basis of the Clown colour morph found in captive-bred ball pythons (Python regius) to study skin pigmentation and patterning in snakes. We obtained samples by crowdsourcing shed skin from commercial breeders and hobbyists. We applied a case-control design, whole-genome pool sequencing, variant annotation, histological analyses, and electron microscopy imaging. We identified a missense mutation in a transmembrane region of the melanocortin-1 receptor (MC1R) associated with the Clown phenotype. In classic avian and mammalian model species, MC1R is known for controlling the type and amount of melanin produced. In contrast, our results suggest that MC1R signalling might play a key role in pattern formation in ball pythons, affecting xanthophore-melanophore distribution. This work highlights the varied functions of MC1R across different vertebrate lineages and promotes a novel model system to study reptile colouration.

PAX3 haploinsufficiency in Maine Coon cats with dominant blue eyes and hearing loss resembling the human Waardenburg syndrome

This study investigated the dominant blue eyes (DBE) trait linked to hearing impairment and variable white spotting in Maine Coon cats. Fifty-eight animals descending from 2 different DBE lineages, the Dutch and the Topaz lines, were sampled. They comprised 48 cats from the Dutch bloodline, including 9 green-eyed and 31 blue-eyed cats, with some individuals exhibiting signs of deafness, and 8 stillborn kittens. Samples from the Topaz lineage included 10 blue-eyed animals. A brainstem auditory evoked response test revealed a reduced to absent response to auditory stimuli and absent physiological waveforms in all of the 8 examined DBE animals. We sequenced the genome of 2 affected cats from the Dutch line and searched for variants in 19 candidate genes for the human Waardenburg syndrome and pigmentary disorders. This search yielded 9 private protein-changing candidate variants in the genes PAX3, EDN3, KIT, OCA2, SLC24A5, HERC2, and TYRP1. The genotype-phenotype cosegregation was observed for the PAX3 variant within all animals from the Dutch lineage. The mutant allele was absent from 461 control genomes and 241 additionally genotyped green-eyed Maine Coons. We considered the PAX3 variant as the most plausible candidate-a heterozygous nonsense single base pair substitution in exon 6 of PAX3 (NC_051841.1:g.205,787,310G>A, XM_019838731.3:c.937C>T, XP_019694290.1:p.Gln313*), predicted to result in a premature stop codon. PAX3 variants cause auditory-pigmentary syndrome in humans, horses, and mice. Together with the comparative data from other species, our findings strongly suggest PAX3:c.937C>T (OMIA:001688-9685) as the most likely candidate variant for the DBE, deafness, and minimal white spotting in the Maine Coon Dutch line. Finally, we propose the designation of DBERE (Rociri Elvis Dominant Blue Eyes) allele in the domestic cat.

A PAX3 insertion in the Celestial breed and certain feline breeding lines with dominant blue eyes

During the last 60 years many inherited traits in domestic outbred cats were selected and retained giving birth to new breeds characterised by singular coat or morphological phenotypes. Among them, minimal white spotting associated with blue eyes was selected by feline breeders to create the Altai, Topaz, and Celestial breeds. Various established breeds also introduced this trait in their lineages. The trait, that was confirmed as autosomal dominant by breeding data, was first described in domestic cats from Kazakhstan and Russia, in British shorthair and British longhair from Russia, and in Maine Coon cats from the Netherlands, suggesting different founding effects. Using a genome-wide association study we identified a single region on chromosome C1 that was associated with the minimal white spotting and blue eyes phenotype (also called DBE by breeders for dominant blue eyes) in the French Celestial breed. Within that region we identified Paired Box 3 (PAX3) as the strongest candidate gene, since PAX3 is a key regulator of MITF (Melanocyte-Inducing Transcription Factor) and PAX3 variants have been previously identified in various species showing white spotting with or without blue eyes including the mouse and the horse. Whole genome sequencing of a Celestial cat revealed an endogenous retrovirus LTR (long terminal repeat) insertion within PAX3 intron 4 known to contain regulatory sequences (conserved non-coding element [CNE]) involved in PAX3 expression. The insertion is in the vicinity of CNE2 and CNE3. All 52 Celestial and Celestial-mixed cats with a DBE phenotype presented the insertion, that was absent in their 22 non-DBE littermates and in 87 non-DBE cats from various breeds. The outbred Celestial founder was also heterozygous for the insertion. Additionally, the variant was found in nine DBE Maine Coon cats related to the Celestial founder and four DBE Siberian cats with an uncertain origin. Segregation of the variant in the Celestial breed is consistent with dominant inheritance and does not appear to be associated with deafness. We propose that this NC_018730.3:g.206974029_206974030insN[395] variant represents the DBECEL (Celestial Dominant Blue Eyes) allele in the domestic cat.

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.

Island demographics and trait associations in white-tailed deer

When a population is isolated and composed of few individuals, genetic drift is the paramount evolutionary force and results in the loss of genetic diversity. Inbreeding might also occur, resulting in genomic regions that are identical by descent, manifesting as runs of homozygosity (ROHs) and the expression of recessive traits. Likewise, the genes underlying traits of interest can be revealed by comparing fixed SNPs and divergent haplotypes between affected and unaffected individuals. Populations of white-tailed deer (Odocoileus virginianus) on islands of Saint Pierre and Miquelon (SPM, France) have high incidences of leucism and malocclusions, both considered genetic defects; on the Florida Keys islands (USA) deer exhibit smaller body sizes, a polygenic trait. Here we aimed to reconstruct island demography and identify the genes associated with these traits in a pseudo case-control design. The two island populations showed reduced levels of genomic diversity and a build-up of deleterious mutations compared to mainland deer; there was also significant genome-wide divergence in Key deer. Key deer showed higher inbreeding levels, but not longer ROHs, consistent with long-term isolation. We identified multiple trait-related genes in ROHs including LAMTOR2 which has links to pigmentation changes, and NPVF which is linked to craniofacial abnormalities. Our mixed approach of linking ROHs, fixed SNPs and haplotypes matched a high number (~50) of a-priori body size candidate genes in Key deer. This suite of biomarkers and candidate genes should prove useful for population monitoring, noting all three phenotypes show patterns consistent with a complex trait and non-Mendelian inheritance.

Different Founding Effects Underlie Dominant Blue Eyes (DBE) in the Domestic Cat

During the last twenty years, minimal white spotting associated with blue eyes was selected by feline breeders to create the Altai, Topaz, and Celestial breeds. Additionally, certain breeders introduced this trait in their lineages of purebred cats. The trait has been called "dominant blue eyes (DBE)" and was confirmed to be autosomal dominant in all lineages. DBE was initially described in outbred cats from Kazakhstan and Russia and in two purebred lineages of British cats from Russia, as well as in Dutch Maine Coon cats, suggesting different founding effects. We have previously identified two variants in the Paired Box 3 (PAX3) gene associated with DBE in Maine Coon and Celestial cats; however, the presence of an underlying variant remains undetermined in other DBE breeding lines. Using a genome-wide association study, we identified a single region on chromosome C1 that was associated with DBE in British cats. Within that region, we identified PAX3 as the strongest candidate gene. Whole-genome sequencing of a DBE cat revealed an RD-114 retrovirus LTR (long terminal repeat) insertion within PAX3 intron 4 (namely NC_018730.3:g.206975776_206975777insN[433]) known to contain regulatory sequences. Using a panel of 117 DBE cats, we showed that this variant was fully associated with DBE in two British lineages, in Altai cats, and in some other DBE lineages. We propose that this NC_018730.3:g.206975776_206975777insN[433] variant represents the DBEALT (Altai Dominant Blue Eye) allele in the domestic cat. Finally, we genotyped DBE cats from 14 lineages for the three PAX3 variants and showed that they were not present in four lineages, confirming genetic heterogeneity of the DBE trait in the domestic cat.

Coloration in Equine: Overview of Candidate Genes Associated with Coat Color Phenotypes

Variation in coat color among equids has attracted significant interest in genetics and breeding research. The range of colors is primarily determined by the type, concentration, and distribution of melanin pigments, with the balance between eumelanin and pheomelanin influenced by numerous genetic factors. Advances in genomic and sequencing technologies have enabled the identification of several candidate genes that influence coat color, thereby clarifying the genetic basis of these diverse phenotypes. In this review, we concisely categorize coat coloration in horses and donkeys, focusing on the biosynthesis and types of melanin involved in pigmentation. Moreover, we highlight the regulatory roles of some key candidate genes, such as MC1R, TYR, MITF, ASIP, and KIT, in coat color variation. Moreover, the review explores how coat color relates to selective breeding and specific equine diseases, offering valuable insights for developing breeding strategies that enhance both the esthetic and health aspects of equine species.

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.

The Hoof Color of Australian White Sheep Is Associated with Genetic Variation of the MITF Gene

Studying the characteristics of mammalian hoof colors is important for genetic improvements in animals. A deeper black hoof color is the standard for breeding purebred Australian White (AUW) sheep and this phenotype could be used as a phenotypic marker of purebred animals. We conducted a genome-wide association study (GWAS) analysis using restriction site associated DNA sequencing (RAD-seq) data from 577 Australian White sheep (black hoof color = 283, grey hoof color = 106, amber hoof color = 186) and performed association analysis utilizing the mixed linear model in EMMAX. The results of GWAS demonstrated that a specific single-nucleotide polymorphism (SNP; g. 33097911G>A) in intron 14 of the microphthalmia-associated transcription factor (MITF) gene was significantly associated with the hoof color in AUW sheep (p = 9.40 × 10-36). The MITF gene plays a key role in the development, differentiation, and functional regulation of melanocytes. Furthermore, the association between this locus and hoof color was validated in a cohort of 212 individuals (black hoof color = 122, grey hoof color = 38, amber hoof color = 52). The results indicated that the hoof color of AUW sheep with GG, AG, and AA genotypes tended to be black, grey, and amber, respectively. This study provided novel insights into hoof color genetics in AUW sheep, enhancing our comprehension of the genetic mechanisms underlying the diverse range of hoof colors. Our results agree with previous studies and provide molecular markers for marker-assisted selection for hoof color in sheep.

A gene regulatory network combining Pax3/7, Sox10 and Mitf generates diverse pigment cell types in medaka and zebrafish

Neural crest cells generate numerous derivatives, including pigment cells, and are a model for studying how fate specification from multipotent progenitors is controlled. In mammals, the core gene regulatory network for melanocytes (their only pigment cell type) contains three transcription factors, Sox10, Pax3 and Mitf, with the latter considered a master regulator of melanocyte development. In teleosts, which have three to four pigment cell types (melanophores, iridophores and xanthophores, plus leucophores e.g. in medaka), gene regulatory networks governing fate specification are poorly understood, although Mitf function is considered conserved. Here, we show that the regulatory relationships between Sox10, Pax3 and Mitf are conserved in zebrafish, but the role for Mitf is more complex than previously emphasized, affecting xanthophore development too. Similarly, medaka Mitf is necessary for melanophore, xanthophore and leucophore formation. Furthermore, expression patterns and mutant phenotypes of pax3 and pax7 suggest that Pax3 and Pax7 act sequentially, activating mitf expression. Pax7 modulates Mitf function, driving co-expressing cells to differentiate as xanthophores and leucophores rather than melanophores. We propose that pigment cell fate specification should be considered to result from the combinatorial activity of Mitf with other transcription factors.

Two Novel Variants in MITF and PAX3 Associated With Splashed White Phenotypes in Horses

Mutations causing depigmentation are relatively common in Equus caballus (horse). Over 40 alleles in multiple genes are associated with increased white spotting (as of February 2023). The splashed white phenotype, a coat spotting pattern described as appearing like the horse has been splashed with white paint, was previously associated with variants in the PAX3 and MITF genes. Both genes encode transcription factors known to control melanocyte migration and pigmentation. We report two novel mutations, a stop-gain mutation in PAX3 (XM_005610643.3:c.927C>T, ECA6:11,196,181, EquCab3.0) and a missense mutation in a binding domain of MITF (NM_001163874.1:c.993A>T, ECA16:21,559,940, EquCab3.0), each with a strong association with increased depigmentation in Pura Raza Española horses (P = 1.144E-11, N = 30, P = 4.441E-16, N = 39 respectively). Using a quantitative method to score depigmentation, the PAX3 and MITF mutations were found to have average white scores of 25.50 and 24.45, respectively, compared to the average white coat spotting score of 1.89 in the control set. The functional impact for each mutation was predicted to be moderate to extreme (I-TASSER, SMART, Variant Effect Predictor, SIFT). We propose to designate the MITF mutant allele as Splashed White 9 and the PAX3 mutant allele as Splashed White 10 per convention.

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.

MITF Contributes to the Body Color Differentiation of Sea Cucumbers Apostichopus japonicus through Expression Differences and Regulation of Downstream Genes

Melanin, which is a pigment produced in melanocytes, is an important contributor to sea cucumber body color. MITF is one of the most critical genes in melanocyte development and melanin synthesis pathways. However, how MITF regulates body color and differentiation in sea cucumbers is poorly understood. In this study, we analyzed the expression level and location of MITF in white, purple, and green sea cucumbers and identified the genes regulated by MITF using chromatin immunoprecipitation followed by sequencing. The mRNA and protein expression levels of MITF were all highest in purple morphs and lowest in white morphs. In situ hybridization indicated that MITF mRNA were mainly expressed in the epidermis. We also identified 984, 732, and 1191 peaks of MITF binding in green, purple, and white sea cucumbers, which were associated with 727, 557, and 887 genes, respectively. Our findings suggested that MITF contributed to the body color differentiation of green, purple, and white sea cucumbers through expression differences and regulation of downstream genes. These results provided a basis for future studies to determine the mechanisms underlying body color formation and provided insights into gene regulation in sea cucumbers.

PMEL p.Leu18del dilutes coat color of Kumamoto sub-breed of Japanese Brown cattle

Background

Coat color is important for registration and maintenance of livestock. Standard coat color of Kumamoto sub-breed of Japanese Brown cattle is solid brown, but individuals with diluted coat color have been observed recently. In this study, we attempted to identify polymorphism(s) responsible for coat color dilution by whole genome analysis.

Results

One of the diluted cattle possessed 7302 exonic polymorphisms which could affect genes’ function. Among them, 14 polymorphisms in 10 coat color-related genes were assumed to be specific for the diluted cattle. Subsequent genotyping with three diluted cattle and 74 standard cattle elucidated that PMEL p.Leu18del was the causative polymorphism for coat color dilution in this sub-breed. Individuals with del/del type of this polymorphism showed diluted coat color, but coat color of heterozygotes were intermediate with various dilution rates.

Conclusions

Coat color dilution of Kumamoto sub-breed was caused by PMEL p.Leu18del. The causative del allele has been detected in several genetically distant cattle breeds, suggesting that PMEL p.Leu18del can be used as a DNA marker to control cattle coat color.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-022-08916-8.

De-novo and genome-wide meta-analyses identify a risk haplotype for congenital sensorineural deafness in Dalmatian dogs

Congenital sensorineural deafness (CSD) has been reported to affect up to 30% of Dalmatian dogs world-wide and while unilaterally deaf dogs can live a close to normal life, dogs suffering bilateral deafness are frequently euthanized. Extreme-white coat patterning as encoded by the gene Melanocyte Inducing Transcription Factor (MITF) has long been postulated as the major risk factor for CSD in the Dalmatian breed. While attempts to identify causative risk variants associated with CSD have been numerous, no genome-wide association study has positively identified MITF as a risk locus for either bilateral or unilateral deafness in the Dalmatian breed to date. In this study, we identified an association with CSD on CFA20 in the vicinity of MITF within Australian Dalmatian dogs. Although not genome-wide significant, the association signal was validated by reanalysing publicly available data and merging the wider data resource with the local data to improve statistical power. The merged data, representing three major global populations of Dalmatian dogs, enabled us to identify a single, well-defined genome-wide significant risk haplotype for CSD. The haplotype was formed by three genome-wide significant associated markers (BICF2G630233852T>C, BICF2G630233861T>C, BICF2G630233888G>A) on CFA20 with 62% of bilaterally deaf dogs homozygous for the risk haplotype (CCA), while 30% of bilaterally deaf and 45% of hearing dogs carried one copy of the risk haplotype. Animals homozygous or heterozygous for the low-risk haplotype were less likely to be unilaterally deaf. While the association between the risk haplotype and deafness is incomplete, animals homozygous for the risk haplotype were 10-times more likely to be bilaterally deaf. Although the underlying causative variants are yet to be discovered, results from this study can now assist with reducing deafness in Dalmatian dogs.

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.

Genome-Wide Analysis Identifies Candidate Genes Encoding Beak Color of Duck

Beak color diversity is a broadly occurring phenomenon in birds. Here, we used ducks to identify candidate genes for yellow, black, and spotted beaks. For this, an F2 population consisting of 275 ducks was genotyped using whole genome resequencing containing 12.6 M single-nucleotide polymorphisms (SNPs) and three beak colors. Genome-wide association studies (GWAS) was used to identify the candidate and potential SNPs for three beak colors in ducks (yellow, spotted, and black). The results showed that 2753 significant SNPs were associated with black beaks, 7462 with yellow, and 17 potential SNPs with spotted beaks. Based on SNP annotation, MITF, EDNRB2, members of the POU family, and the SLC superfamily were the candidate genes regulating pigmentation. Meanwhile, isoforms MITF-M and EDNRB2 were significantly different between black and yellow beaks. MITF and EDNRB2 likely play a synergistic role in the regulation of melanin synthesis, and their mutations contribute to phenotypic differences in beak melanin deposition among individuals. This study provides new insights into genetic factors that may influence the diversity of beak color.

Selection signatures and formation of the Samosir goat breed through the cultures of the Batak Toba Tribe in Samosir Island, Indonesia

Background and Aim:

The Samosir goat has a high cultural value and is a source of germplasm in Indonesia. This study aimed to reveal the history and selection signatures of the Samosir goat.

Materials and Methods:

A total of 25 goats were divided into seven subpopulations of Indonesian goat breeds. Deoxyribonucleic acid (DNA) from blood samples was isolated with the use of the gSYNC™ DNA Mini Kit (Geneaid, Taipei, Taiwan). Cytb gene amplification was performed by the polymerase chain reaction (PCR) method, and the PCR products were sequenced. A phylogenetic tree was constructed by the neighbor-joining method using MEGA 11 software. A questionnaire was used to collect information related to the history and breeding practices of the Samosir goat on Samosir Island.

Results:

Samosir goats are divided into four groups based on their coat color: Completely white, white with brown spots, white with black spots, and white with brown and black spots. The body form of the Samosir goat is similar to that of the Kacang goat. The space below a traditional Toba Batak house is used as a goat pen. The genetic difference between the Samosir goat and the Kacang goat based on the Cytb gene was approximately 0.1%.

Conclusion:

Phylogenetic analysis between Samosir goats and other indigenous Indonesian goats revealed that Samosir goats form a single clade, with a very close genetic distance from other local goats, such as the Kacang goat. The Toba Batak culture on Samosir Island has significantly influenced the selection and formation of the Samosir goat breed.

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.

Canine coat pigmentation genetics: a review

Our understanding of canine coat colour genetics and the associated health implications is developing rapidly. To date, there are 15 genes with known roles in canine coat colour phenotypes. Many coat phenotypes result from complex and/or epistatic genetic interactions among variants within and between loci, some of which remain unidentified. Some genes involved in canine pigmentation have been linked to aural, visual and neurological impairments. Consequently, coat pigmentation in the domestic dog retains considerable ethical and economic interest. In this paper we discuss coat colour phenotypes in the domestic dog, the genes and variants responsible for these phenotypes and any proven coat colour-associated health effects.

Prevalence of Varied Coat Coloration in a Yellow-Throated Marten (Martes flavigula) Population

Simple Summary

Abnormal coloration is very rare in any given population of wildlife; however, our research identified a yellow-throated marten population with a high ratio of this phenomenon for the first time. Across the main distribution of the species with relevant observational data, we observed abnormally-colored martens in only Northeast Tiger and Leopard National Park. Abnormal coloration had a variety of forms and individuals with white paws that accounted for a larger proportion of the overall population than normal individuals. This shows heritable variation in the region, which is worthy of further research.

Abstract

Mammalian coat color is determined by heritable variations such as disease, nutrition, and hormone levels. Variation in animal coat color is also considered an environmental indicator and provides clues for the study of population genetics and biogeography. Records of abnormal coloration in the wild are rare, not only because it is often selected against, but also because of the difficulties in detection of the phenomenon. We used long-term camera-trapping data to first report abnormal coat coloration in yellow-throated marten (Martes flavigula) in China. Six types of abnormal coloration were found only in the Northeast Tiger and Leopard National Park, Northeast China, which were not reported in other regions in China. A total of 268 videos of Martes flavigula contained normal coloration, 455 videos of individuals of the species contained abnormal coloration, 437 contained the ‘gloving’ type (martens with de-pigmented front toes, paws or lower forelimbs), while the remaining other 18 videos contained five types (different degrees of white-spotting and dilution). The higher relative abundance index (0.428, ‘gloving’ to 0.329, normal) and wide distribution area of the ‘gloving’ type indicated that this abnormal coat coloration type is usual in Northeast China, which may reflect genetic variability in the local population. These records will contribute to further research on animal coat color and its corresponding adaptive strategy.

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 Basis of Striking Fin Shapes and Colors in the Fighting Fish

Resolving the genomic basis underlying phenotypic variations is a question of great importance in evolutionary biology. However, understanding how genotypes determine the phenotypes is still challenging. Centuries of artificial selective breeding for beauty and aggression resulted in a plethora of colors, long-fin varieties, and hyper-aggressive behavior in the air-breathing Siamese fighting fish (Betta splendens), supplying an excellent system for studying the genomic basis of phenotypic variations. Combining whole-genome sequencing, quantitative trait loci mapping, genome-wide association studies, and genome editing, we investigated the genomic basis of huge morphological variation in fins and striking differences in coloration in the fighting fish. Results revealed that the double tail, elephant ear, albino, and fin spot mutants each were determined by single major-effect loci. The elephant ear phenotype was likely related to differential expression of a potassium ion channel gene, kcnh8. The albinotic phenotype was likely linked to a cis-regulatory element acting on the mitfa gene and the double-tail mutant was suggested to be caused by a deletion in a zic1/zic4 coenhancer. Our data highlight that major loci and cis-regulatory elements play important roles in bringing about phenotypic innovations and establish Bettas as new powerful model to study the genomic basis of evolved changes.

Detecting the footprint of selection on the genomes of Murciano-Granadina goats

Artificial selection is one of the major forces modifying the genetic composition of livestock populations. Identifying genes under selection could be useful to elucidate their impact on phenotypic variation. We aimed to identify genomic regions targeted by selection for dairy and pigmentation traits in Murciano-Granadina goats. Performance of a selection scan based on the integrated haplotype score test in a population of 1183 Murciano-Granadina goats resulted in the identification of 77 candidate genomic regions/SNPs. The most significant selective sweeps mapped to chromosomes 1 (69.86 Mb), 4 (41.80-49.95 Mb), 11 (65.74 Mb), 12 (31.24 and 52.51 Mb), 17 (34.76-37.67 Mb), 22 (31.75 Mb), and 26 (26.69-31.05 Mb). By using previously generated RNA-Seq data, we built a catalogue of 6414 genes that are differentially expressed across goat lactation (i.e. 78 days post-partum, early lactation; 216 days post-partum, late lactation; 285 days post-partum, dry period). Interestingly, 183 of these genes mapped to selective sweeps and several of them display functions related with lipid, protein, and carbohydrate metabolism, insulin signaling, cell proliferation, as well as mammary development and involution. Of particular interest are the CSN3 and CSN1S2 genes, which encode two major milk proteins. Additionally, we found three pigmentation genes (GLI3, MC1R, and MITF) co-localizing with selective sweeps. Performance of a genome-wide association study and Sanger sequencing and TaqMan genotyping experiments revealed that the c.801C>G (p.Cys267Trp) polymorphism in the melanocortin 1 receptor (MC1R) gene is the main determinant of the black (GG or GC genotypes) and brown (CC genotypes) colorations of Murciano-Granadina goats.

Deubiquitination of MITF-M Regulates Melanocytes Proliferation and Apoptosis

Microphthalmia-associated transcription factor-M (MITF-M) is the key gene in the proliferation and differentiation of melanocytes, which undergoes an array of post-translation modifications. As shown in our previous study, deubiquitinase USP13 is directly involved in melanogenesis. However, it is still ambiguous that the effect of USP13-mediated MITF-M expression on melanocytes proliferation and apoptosis. Herein, we found that MITF-M overexpressing melanocytes showed high cell proliferation, reduced apoptosis, and increased melanin levels. Besides, melanin-related genes, TYR, DCT, GPNMB, and PMEL, were significantly up-regulated in MITF-M overexpressing melanocytes. Furthermore, Exogenous USP13 significantly upregulated the endogenous MITF-M protein level, downregulated USP13 significantly inhibited MITF-M protein levels, without altering MITF-M mRNA expression. In addition, USP13 upregulation mitigated the MITF-M degradation and significantly increased the half-life of MITF-M. Also, USP13 stabilized the exogenous MITF protein levels. In conclusion, the MITF-M level was regulated by USP13 deubiquitinase in melanocytes, affecting melanocytes proliferation and apoptosis. This study provides the theoretical basis for coat color transformation that could be useful in the development of the new breed in fur animals.

A De Novo MITF Deletion Explains a Novel Splashed White Phenotype in an American Paint Horse

Splashed white is a coat color pattern in horses characterized by extensive white patterning on the legs, belly, and face often accompanied by blue eyes and deafness. Three mutations in microphthalmia-associated transcription factor (MITF) and two mutations in Paired Box 3 (PAX3) have been identified that explain splashed white patterns (SW1-SW5). An American Paint Horse stallion with a splashed white phenotype and blue eyes, whose parents were not white patterned, was negative for the 5 known splashed white variants and other known white spotting alleles. This novel splashed white phenotype (SW6) was hypothesized to be caused by a de novo mutation in MITF or PAX3. Analysis of whole-genome sequencing using the EquCab3.0 reference genome for comparison identified an 8.7 kb deletion in MITF on ECA16 (NC_009159.3:g.21551060-21559770del). The deletion encompassed part of intron 7 through the 3' UTR of exon 9 of MITF, including the helix-loop-helix DNA-binding domain (ENSECAT00000006375.3). This variant is predicted to truncate protein and impair binding to DNA. Sanger sequencing confirmed the stallion was heterozygous for the MITF deletion. No single nucleotide polymorphisms (SNPs) or structural variants were identified in PAX3 or any of the other candidate genes that were unique to the stallion or predicted to affect protein function. Genotyping five of the stallion's splashed white offspring, including one all white foal, found that they were also heterozygous for the deletion. Given the role of MITF in producing white pattern phenotypes, and the predicted deleterious effect of this mutation, this 8.7 kb deletion is the likely causal variant for SW6.

An atlas of CNV maps in cattle, goat and sheep

Copy number variation (CNV) is the most prevalent type of genetic structural variation that has been recognized as an important source of phenotypic variation in humans, animals and plants. However, the mechanisms underlying the evolution of CNVs and their function in natural or artificial selection remain unknown. Here, we generated CNV region (CNVR) datasets which were diverged or shared among cattle, goat, and sheep, including 886 individuals from 171 diverse populations. Using 9 environmental factors for genome-wide association study (GWAS), we identified a series of candidate CNVRs, including genes relating to immunity, tick resistance, multi-drug resistance, and muscle development. The number of CNVRs shared between species is significantly higher than expected (P<0.00001), and these CNVRs may be more persist than the single nucleotide polymorphisms (SNPs) shared between species. We also identified genomic regions under long-term balancing selection and uncovered the potential diversity of the selected CNVRs close to the important functional genes. This study provides the evidence that balancing selection might be more common in mammals than previously considered, and might play an important role in the daily activities of these ruminant species.

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.

Transcript Profiles of Stria Vascularis in Models of Waardenburg Syndrome

Background

Waardenburg syndrome is an uncommon genetic condition characterized by at least some degree of congenital hearing loss and pigmentation deficiencies. However, the genetic pathway affecting the development of stria vascularis is not fully illustrated.

Methods

The transcript profile of stria vascularis of Waardenburg syndrome was studied using Mitf-M mutant pig and mice models. Therefore, GO analysis was performed to identify the differential gene expression caused by Mitf-M mutation.

Results

There were 113 genes in tyrosine metabolism, melanin formation, and ion transportations showed significant changes in pig models and 191 genes in mice models. In addition, there were some spice's specific gene changes in the stria vascularis in the mouse and porcine models. The expression of tight junction-associated genes, including Cadm1, Cldn11, Pcdh1, Pcdh19, and Cdh24 genes, were significantly higher in porcine models compared to mouse models. Vascular-related and ion channel-related genes in the stria vascularis were also shown significantly difference between the two species. The expression of Col2a1, Col3a1, Col11a1, and Col11a2 genes were higher, and the expression of Col8a2, Cd34, and Ncam genes were lower in the porcine models compared to mouse models.

Conclusions

Our data suggests that there is a significant difference on the gene expression and function between these two models.

Identification of key microRNAs affecting melanogenesis of breast muscle in Muchuan black-boned chickens by RNA sequencing

1. Melanin content is considered an important indicator of meat quality in black-boned chickens, which have a high market value. To understand the complex physiological processes underlying muscle melanogenesis in this chicken, differentially expressed miRNAs (DEMs) were detected between black muscle (BM) and white muscle (WM) of chickens using high-throughput sequencing technology. Six small RNA libraries were constructed, and more than 16.75 million clean reads were obtained for each library. 2. A total of 582 known miRNAs and 65 novel miRNAs were identified from the six chicken sequence libraries. A total of 19 DEMs were identified between the two groups, of which nine were upregulated and 10 were downregulated. Furthermore, the DEMs were predicted to target 572 genes. 3. Certain DEMs (such as miR-204, miR-133b, and miR-12 229-3p) and their target genes may play an important role in muscle melanogenesis of chickens. These findings provide a foundation for clarifying the miRNA regulatory mechanisms involved in muscle pigmentation in avian species.

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.

Equine Genotyping Arrays

High-quality genomic tools have been integral in understanding genomic architecture and function in the modern-day horse. The equine genetics community has a long tradition of pooling resources to develop genomic tools. Since the equine genome was sequenced in 2006, several iterations of high throughput genotyping arrays have been developed and released, enabling rapid and cost-effective genotyping. This review highlights the design considerations of each iteration, focusing on data available during development and outlining considerations in selecting the genetic variants included on each array. Additionally, we outline recent applications of equine genotyping arrays as well as future prospects and applications.

Genetic Diversity and Signatures of Selection in a Native Italian Horse Breed Based on SNP Data

Simple Summary

The Bardigiano horse is a native Italian breed bred for living in rural areas, traditionally used in agriculture. The breed counts about 3000 horses, and it is nowadays mainly used for recreational purposes. The relatively small size and the closed status of the breed raise the issue of monitoring genetic diversity. We therefore characterized the breed’s genetic diversity based on molecular data. We showed a critical reduction of genetic variability mainly driven by past bottlenecks. We also highlighted homozygous genomic regions that might be the outcome of directional selection in recent years, in line with the conversion of Bardigiano horses from agricultural to riding purposes.

Abstract

Horses are nowadays mainly used for sport and leisure activities, and several local breeds, traditionally used in agriculture, have been exposed to a dramatic loss in population size and genetic diversity. The loss of genetic diversity negatively impacts individual fitness and reduces the potential long-term survivability of a breed. Recent advances in molecular biology and bioinformatics have allowed researchers to explore biodiversity one step further. This study aimed to evaluate the loss of genetic variability and identify genomic regions under selection pressure in the Bardigiano breed based on GGP Equine70k SNP data. The effective population size based on Linkage Disequilibrium (N_e) was equal to 39 horses, and it showed a decline over time. The average inbreeding based on runs of homozygosity (ROH) was equal to 0.17 (SD = 0.03). The majority of the ROH were relatively short (91% were ≤2 Mbp long), highlighting the occurrence of older inbreeding, rather than a more recent occurrence. A total of eight ROH islands, shared among more than 70% of the Bardigiano horses, were found. Four of them mapped to known quantitative trait loci related to morphological traits (e.g., body size and coat color) and disease susceptibility. This study provided the first genome-wide scan of genetic diversity and selection signatures in an Italian native horse breed.

The expression of miR-129-5p and its target genes in the skin of goats

Coat color is one of the major quality traits of animals, and miR-129-5p acts as an important regulator for melanin biosynthesis in mammals. In this study, real-time PCR and western blotting were used to examine the expression of miR-129-5p and its targets genes in the skin of different coat color goats. The results showed that the expression of miR-129-5p in the skin samples of Inner Mongolia cashmere goats (IMCG) was higher than that of Dazu black goat (DBG). Also, the target genes (tyrosinase (TYR), frizzled 6 (FZD6) and glycogen synthase kinase 3β (GSK3β)) of miR-129-5p was highly expressed in the skin samples of DBG. The expression of miR-129-5p firstly increased and then decreased with age in F1 hybrid generation of DBG and IMCG. In addition, the expression of TYR decreased with age, while the expression of MITF increased with age but then decreased. The expression of FZD6 and GSK3β in the skin samples of F1 of different ages were irregular. Our results indicated that miR-129-5p mainly affects the formation of coat color of goats by decreasing the expression of TYR. This study suggests that miR-129-5p can act as a suppressor in the formation of coat color to lay the foundation for studying the effect of miR-129-5p on melanin synthesis.

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.

Integrated Analysis of mRNA Expression, CpG Island Methylation, and Polymorphisms in the MITF Gene in Ducks (Anas platyrhynchos)

Microphthalmia-associated transcription factor (MITF) is a key regulator for the development and function of melanocytes in skin, eye, and plumage pigmentations. Thus, the MITF was selected as a candidate gene associated with plumage coloration in ducks. This study analyzed the mRNA expression, promoter methylation, and polymorphisms in the MITF gene in ducks with different plumage colors (Putian Black, Putian White, Liancheng White, and Longsheng Jade-green). No expression of the MITF melanin-specific isoform (MITF-M) was detected in white feather bulbs. By contrast, the mRNA expression levels of MITF-M were high in black feather bulbs. Bioinformatics analysis showed that two CpG islands were present in the promoter region of the MITF gene. The methylation level of the second CpG island was significantly lower in black feather bulbs than in white feather bulbs. However, the methylation level of the first CpG island was not different among the feather bulbs with various colors except Liancheng White feather bulbs. The methylation status of the whole CpG island significantly and negatively correlated with the mRNA expression of MITF-M (P < 0.05). Furthermore, four novel SNPs (single nucleotide polymorphisms) were identified in the 5′UTR, exon 4, intron 7, and intron 8 of the MITF gene. Allele T in g.39807T>G and allele G in g.40862G>A were the predominant alleles only found in Putian White, whereas the variant A allele in g.32813G>A exhibited a high allele frequency in Liancheng White. Collectively, these results contributed to the understanding of the function of the MITF gene in duck plumage coloration.

Frameshift Variant in MFSD12 Explains the Mushroom Coat Color Dilution in Shetland Ponies

Mushroom is a unique coat color phenotype in Shetland Ponies characterized by the dilution of the chestnut coat color to a sepia tone and is hypothesized to be a recessive trait. A genome wide association study (GWAS), utilizing the Affymetrix 670K array (MNEc670k) and a single locus mixed linear model analysis (EMMAX), identified a locus on ECA7 for further investigation (P_corrected = 2.08 × 10⁻¹⁰). This locus contained a 3 Mb run of homozygosity in the 12 mushroom ponies tested. Analysis of high throughput Illumina sequencing data from one mushroom Shetland pony compared to 87 genomes from horses of various breeds, uncovered a frameshift variant, p.Asp201fs, in the MFSD12 gene encoding the major facilitator superfamily domain containing 12 protein. This variant was perfectly concordant with phenotype in 96 Shetland Ponies (P = 1.15 × 10⁻²²), was identified in the closely related Miniature Horse for which the mushroom phenotype is suspected to occur (f_mu = 0.02), and was absent in 252 individuals from seven additional breeds not reported to have the mushroom phenotype. MFSD12 is highly expressed in melanocytes and variants in this gene in humans, mice, and dogs impact pigmentation. Given the role of MFSD12 in melanogenesis, we propose that p.Asp201fs is causal for the dilution observed in mushroom ponies.