Equine vitiligo-like depigmentation in grey horses is related to genes involved in immune response and tumor metastasis.

BACKGROUND: In horses, the autoimmune disease vitiligo is characterized by the loss of melanocytes and results in patchy depigmentation of the skin around the eyes, muzzle and the perianal region. Vitiligo-like depigmentation occurs predominantly in horses displaying the grey coat colour and is observed at a prevalence level of 26.0-67.0% in grey horses compared with only 0.8-3.5% in non-grey horses. While the polygenetic background of this complex disease is well documented in humans, the underlying candidate genes for this skin disorder in horses remain unknown. In this study we aim to perform a genome-wide association study (GWAS) for identifying putative candidate loci for vitiligo-like depigmentation in horses. METHODS: In the current study, we performed a GWAS analysis using high-density 670 k single nucleotide polymorphism (SNP) data from 152 Lipizzan and 104 Noriker horses, which were phenotyped for vitiligo-like depigmentation by visual inspection. After quality control 376,219 SNPs remained for analyses, the genome-wide Bonferroni corrected significance level was p < 1.33e-7. RESULTS: We identified seven candidate genes on four chromosomes (ECA1, ECA13, ECA17, ECA20) putatively involved in vitiligo pathogenesis in grey horses. The highlighted genes PHF11, SETDB2, CARHSP1 and LITAFD, are associated with the innate immune system, while the genes RCBTB1, LITAFD, NUBPL, PTP4A1, play a role in tumor suppression and metastasis. The antagonistic pathogenesis of vitiligo in relation to cancer specific enhanced cell motility and/or metastasis on typical melanoma predilection sites underlines a plausible involvement of RCBTB1, LITAFD, NUBPL, and PTP4A1. CONCLUSIONS: The proposed candidate genes for equine vitiligo-like depigmentation, indicate an antagonistic relation between vitiligo and tumor metastasis in a horse population with higher incidence of melanoma. Further replication and expression studies should lead to a better understanding of this skin disorder in horses.

High-resolution population structure and runs of homozygosity reveal the genetic architecture of complex traits in the Lipizzan horse.

BACKGROUND: The sample ascertainment bias due to complex population structures remains a major challenge in genome-wide investigations of complex traits. In this study we derived the high-resolution population structure and levels of autozygosity of 377 Lipizzan horses originating from five different European stud farms utilizing the SNP genotype information of the high density 700 k Affymetrix Axiom™ Equine genotyping array. Scanning the genome for overlapping runs of homozygosity (ROH) shared by more than 50% of horses, we identified homozygous regions (ROH islands) in order to investigate the gene content of those candidate regions by gene ontology and enrichment analyses. RESULTS: The high-resolution population network approach revealed well-defined substructures according to the origin of the horses (Austria, Slovakia, Croatia and Hungary). The highest mean genome coverage of ROH (SROH) was identified in the Austrian (SROH = 342.9), followed by Croatian (SROH = 214.7), Slovakian (SROH = 205.1) and Hungarian (SROH = 171.5) subpopulations. ROH island analysis revealed five common islands on ECA11 and ECA14, hereby confirming a closer genetic relationship between the Hungarian and Croatian as well as between the Austrian and Slovakian samples. Private islands were detected for the Hungarian and the Austrian Lipizzan subpopulations. All subpopulations shared a homozygous region on ECA11, nearly identical in position and length containing among other genes the homeobox-B cluster, which was also significantly (p < 0.001) highlighted by enrichment analysis. Gene ontology terms were mostly related to biological processes involved in embryonic morphogenesis and anterior/posterior specification. Around the STX17 gene (causative for greying), we identified a ROH island harbouring the genes NR4A3, STX17, ERP44 and INVS. Within further islands on ECA14, ECA16 and ECA20 we detected the genes SPRY4, NDFIP1, IMPDH2, HSP90AB1, whereas SPRY4 and HSP90AB1 are involved in melanoma metastasis and survival rate of melanoma patients in humans. CONCLUSIONS: We demonstrated that the assessment of high-resolution population structures within one single breed supports the downstream genetic analyses (e.g. the identification of ROH islands). By means of ROH island analyses, we identified the genes SPRY4, NDFIP1, IMPDH2, HSP90AB1, which might play an important role for further studies on equine melanoma. Furthermore, our results highlighted the impact of the homeobox-A and B cluster involved in morphogenesis of Lipizzan horses.

Population Networks Associated with Runs of Homozygosity Reveal New Insights into the Breeding History of the Haflinger Horse.

Within the scope of current genetic diversity analyses, population structure and homozygosity measures are independently analyzed and interpreted. To enhance analytical power, we combined the visualization of recently described high-resolution population networks with runs of homozygosity (ROH). In this study, we demonstrate that this approach enabled us to reveal important aspects of the breeding history of the Haflinger horse. We collected high-density genotype information of 531 horses originating from 7 populations which were involved in the formation of the Haflinger, namely 32 Italian Haflingers, 78 Austrian Haflingers, 190 Noriker, 23 Bosnian Mountain Horses, 20 Gidran, 33 Shagya Arabians, and 155 Purebred Arabians. Model-based cluster analysis identified substructures within Purebred Arabian, Haflinger, and Noriker that reflected distinct genealogy (Purebred Arabian), geographic origin (Haflinger), and coat color patterns (Noriker). Analysis of ROH revealed that the 2 Arabian populations (Purebred and Shagya Arabians), Gidran and the Bosnian Mountain Horse had the highest genome proportion covered by ROH segments (306-397 Mb). The Noriker and the Austrian Haflinger showed the lowest ROH coverage (228, 282 Mb). Our combined visualization approach made it feasible to clearly identify outbred (admixture) and inbred (ROH segments) horses. Genomic inbreeding coefficients (FROH) ranged from 10.1% (Noriker) to 17.7% (Purebred Arabian). Finally it could be demonstrated, that the Austrian Haflinger sample has a lack of longer ROH segments and a deviating ROH spectrum, which is associated with past bottleneck events and the recent mating strategy favoring out-crosses within the breed.

Phenotypic and Genetic Analysis of the Leopard Complex Spotting in Noriker Horses.

Genetic analyses of coat colors are frequently restricted to subjectively categorized phenotype information. The aim of this study was to develop a method to numerically quantify the variability of leopard complex (LP) spotting phenotypes introducing tools from image analysis. Generalized Procrustes analysis eliminates systematic errors due to imaging process. The binarization of normalized images and the application of principal component analysis (PCA) on the derived pixel matrices, transform pixel information into numerical data space. We applied these methods on 90 images to ascertain the specific leopard patterns within the Noriker breed. Furthermore, we genotyped a representative sample of 191 Noriker horses for the known LP spotting associated loci. Ninety-seven percentage of the genotyped leopard spotted horses were heterozygous for LP and had at least one copy of the PATN1 allele. However, the remaining pattern variation was great, indicating other genetic factors influencing the expression of LP spotting. Based upon this data, we estimated effect sizes of the modifier PATN1, and additional factors including sex, age, base color, and spotting phenotype of parents. The PCA of the pixel matrix resulted in 2 significant components accounting for 51% of the variation. Applying a linear model, we identified significant effects for age groups and base color on the first and second components, while for sex and parents' LP phenotype significant effects were found on 4 additional components.

Complex inheritance of melanoma and pigmentation of coat and skin in Grey horses.

The dominant phenotype of greying with age in horses, caused by a 4.6-kb duplication in intron 6 of STX17, is associated with a high incidence of melanoma and vitiligo-like skin depigmentation. However, the progressive greying and the incidence of melanoma, vitiligo-like depigmentation, and amount of speckling in these horses do not follow a simple inheritance pattern. To understand their inheritance, we analysed the melanoma grade, grey level, vitiligo grade, and speckling grade of 1,119 Grey horses (7,146 measurements) measured in six countries over a 9-year period. We estimated narrow sense heritability (h(2)), and we decomposed this parameter into polygenic heritability (h(2) (POLY)), heritability due to the Grey (STX17) mutation (h(2) (STX17)), and heritability due to agouti (ASIP) locus (h(2) (ASIP)). A high heritability was found for greying (h(2) = 0.79), vitiligo (h(2) = 0.63), and speckling (h(2) = 0.66), while a moderate heritability was estimated for melanoma (h(2) = 0.37). The additive component of ASIP was significantly different from zero only for melanoma (h(2) (ASIP) = 0.02). STX17 controlled large proportions of phenotypic variance (h(2) (STX17) = 0.18-0.55) and overall heritability (h(2) (STX17)/h(2) = 0.28-0.83) for all traits. Genetic correlations among traits were estimated as moderate to high, primarily due to the effects of the STX17 locus. Nevertheless, the correlation between progressive greying and vitiligo-like depigmentation remained large even after taking into account the effects of STX17. We presented a model where four traits with complex inheritance patterns are strongly influenced by a single mutation. This is in line with evidence of recent studies in domestic animals indicating that some complex traits are, in addition to the large number of genes with small additive effects, influenced by genes of moderate-to-large effect. Furthermore, we demonstrated that the STX17 mutation explains to a large extent the moderate to high genetic correlations among traits, providing an example of strong pleiotropic effects caused by a single gene.

Genetic diversity, population structure and subdivision of local Balkan pig breeds in Austria, Croatia, Serbia and Bosnia-Herzegovina and its practical value in conservation programs.

BACKGROUND: At present the Croatian Turopolje pig population comprises about 157 breeding animals. In Austria, 324 Turopolje pigs originating from six Croatian founder animals are registered. Multiple bottlenecks have occurred in this population, one major one rather recently and several more older and moderate ones. In addition, it has been subdivided into three subpopulations, one in Austria and two in Croatia, with restricted gene flow. These specificities explain the delicate situation of this endangered Croatian lard-type pig breed. METHODS: In order to identify candidate breeding animals or gene pools for future conservation breeding programs, we studied the genetic diversity and population structure of this breed using microsatellite data from 197 individuals belonging to five different breeds. RESULTS: The genetic diversity of the Turopolje pig is dramatically low with observed heterozygosities values ranging from 0.38 to 0.57. Split into three populations since 1994, two genetic clusters could be identified: one highly conserved Croatian gene pool in Turopoljski Lug and the"Posavina" gene pool mainly present in the Austrian population. The second Croatian subpopulation in Lonjsko Polje in the Posavina region shows a constant gene flow from the Turopoljski Lug animals. CONCLUSIONS: One practical conclusion is that it is necessary to develop a "Posavina" boar line to preserve the "Posavina" gene pool and constitute a corresponding population in Croatia. Animals of the highly inbred herd in Turopoljski Lug should not be crossed with animals of other populations since they represent a specific phenotype-genotype combination. However to increase the genetic diversity of this herd, a program to optimize its sex ratio should be carried out, as was done in the Austrian population where the level of heterozygosity has remained moderate despite its heavy bottleneck in 1994.

DPF3, A Putative Candidate Gene For Melanoma Etiopathogenesis in Gray Horses.

Melanoma prevalence in gray horses reaches up to 50% and more. Several studies have documented a genetic melanoma predisposition which is referred to the 4.6 kb duplication in intron 6 of STX17 and its surrounding haplotype. However, the genetic background and mechanisms responsible for differences in etiopathogenesis of equine dermal melanomatosis still remain unknown. In the current study, we performed a genome wide association analysis in 141 Lipizzan horses and subsequently identified one candidate gene on chromosome 24 putatively involved in melanoma pathogenesis in gray horses. The associated SNP was located in the intronic region of DPF3, a gene which is involved in humans in cell growth, proliferation, apoptosis and motility of cancer cells. The replication study in 1210 horses from seven breeds demonstrated, that the G/G genotype of the DPF3 associated SNP exhibits putative melanoma suppression effects. As a conclusion DPF3 represents a candidate gene, which might play an essential role for gray horses coping with high genetic melanoma related tumor load.

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.

Genome-Wide Association Studies Based on Equine Joint Angle Measurements Reveal New QTL Affecting the Conformation of Horses.

The evaluation of conformation traits is an important part of selection for breeding stallions and mares. Some of these judged conformation traits involve joint angles that are associated with performance, health, and longevity. To improve our understanding of the genetic background of joint angles in horses, we have objectively measured the angles of the poll, elbow, carpal, fetlock (front and hind), hip, stifle, and hock joints based on one photograph of each of the 300 Franches-Montagnes (FM) and 224 Lipizzan (LIP) horses. After quality control, genome-wide association studies (GWASs) for these traits were performed on 495 horses, using 374,070 genome-wide single nucleotide polymorphisms (SNPs) in a mixed-effect model. We identified two significant quantitative trait loci (QTL) for the poll angle on ECA28 (p = 1.36 × 10-7), 50 kb downstream of the ALX1 gene, involved in cranial morphology, and for the elbow joint on ECA29 (p = 1.69 × 10-7), 49 kb downstream of the RSU1 gene, and 75 kb upstream of the PTER gene. Both genes are associated with bone mineral density in humans. Furthermore, we identified other suggestive QTL associated with the stifle joint on ECA8 (p = 3.10 × 10-7); the poll on ECA1 (p = 6.83 × 10-7); the fetlock joint of the hind limb on ECA27 (p = 5.42 × 10-7); and the carpal joint angle on ECA3 (p = 6.24 × 10-7), ECA4 (p = 6.07 × 10-7), and ECA7 (p = 8.83 × 10-7). The application of angular measurements in genetic studies may increase our understanding of the underlying genetic effects of important traits in equine breeding.

Genome-Wide Homozygosity Patterns and Evidence for Selection in a Set of European and Near Eastern Horse Breeds.

Intensive artificial and natural selection have shaped substantial variation among European horse breeds. Whereas most equine selection signature studies employ divergent genetic population structures in order to derive specific inter-breed targets of selection, we screened a total of 1476 horses originating from 12 breeds for the loss of genetic diversity by runs of homozygosity (ROH) utilizing a 670,000 single nucleotide polymorphism (SNP) genotyping array. Overlapping homozygous regions (ROH islands) indicating signatures of selection were identified by breed and similarities/dissimilarities between populations were evaluated. In the entire dataset, 180 ROH islands were identified, whilst 100 islands were breed specific, all other overlapped in 36 genomic regions with at least one ROH island of another breed. Furthermore, two ROH hot spots were determined at horse chromosome 3 (ECA3) and ECA11. Besides the confirmation of previously documented target genes involved in selection for coat color (MC1R, STX17, ASIP), body size (LCORL/NCAPG, ZFAT, LASP1, HMGA2), racing ability (PPARGC1A), behavioral traits (GRIN2B, NTM/OPCML) and gait patterns (DMRT3), several putative target genes related to embryonic morphogenesis (HOXB), energy metabolism (IGFBP-1, IGFBP-3), hair follicle morphogenesis (KRT25, KRT27, INTU) and autophagy (RALB) were highlighted. Furthermore, genes were pinpointed which might be involved in environmental adaptation of specific habitats (UVSSA, STXBP4, COX11, HLF, MMD).

Repeatability, reproducibility and consistency of horse shape data and its association with linearly described conformation traits in Franches-Montagnes stallions.

Linear description (LD) of conformation traits was introduced in horse breeding to minimise subjectivity in scoring. However, recent studies have shown that LD traits show essentially the same problems as traditionally scored traits, such as data converging around the mean value with very small standard deviations. To improve the assessment of conformation traits of horses, we investigated the application of the recently described horse shape space model based upon 403 digitised photographs of 243 Franches-Montagnes (FM) stallions and extracted joint angles based on specific landmark triplets. Repeatability, reproducibility and consistency of the resulting shape data and joint angles were assessed with Procrustes ANOVA (Rep) and intra-class correlation coefficients (ICC). Furthermore, we developed a subjective score to classify the posture of the horses on each photograph. We derived relative warp scores (PCs) based upon the digitised photos conducting a principal component analysis (PCA). The PCs of the shapes and joint angles were compared to the posture scores and to the linear description data using linear mixed effect models including significant posture scores as random factors. The digitisation process was highly repeatable and reproducible for the shape (Rep = 0.72-0.99, ICC = 0.99). The consistency of the shape was limited by the age and posture (p < 0.05). The angle measurements were highly repeatable within one digitiser. Between digitisers, we found a higher variability of ICC values (ICC = 0.054-0.92), indicating digitising error in specific landmarks (e.g. shoulder point). The posture scores were highly repeatable (Fleiss' kappa = 0.713-0.857). We identified significant associations (p(X2) < 0.05) with traits describing the withers height, shoulder length and incline, overall leg conformation, walk and trot step length. The horse shape data and angles provide additional information to explore the morphology of horses and therefore can be applied to improve the knowledge of the genetic architecture of LD traits.

Y Chromosome Uncovers the Recent Oriental Origin of Modern Stallions.

The Y chromosome directly reflects male genealogies, but the extremely low Y chromosome sequence diversity in horses has prevented the reconstruction of stallion genealogies [1, 2]. Here, we resolve the first Y chromosome genealogy of modern horses by screening 1.46 Mb of the male-specific region of the Y chromosome (MSY) in 52 horses from 21 breeds. Based on highly accurate pedigree data, we estimated the de novo mutation rate of the horse MSY and showed that various modern horse Y chromosome lineages split much later than the domestication of the species. Apart from few private northern European haplotypes, all modern horse breeds clustered together in a roughly 700-year-old haplogroup that was transmitted to Europe by the import of Oriental stallions. The Oriental horse group consisted of two major subclades: the Original Arabian lineage and the Turkoman horse lineage. We show that the English Thoroughbred MSY was derived from the Turkoman lineage and that English Thoroughbred sires are largely responsible for the predominance of this haplotype in modern horses.

Identification of genetic variation on the horse y chromosome and the tracing of male founder lineages in modern breeds.

The paternally inherited Y chromosome displays the population genetic history of males. While modern domestic horses (Equus caballus) exhibit abundant diversity within maternally inherited mitochondrial DNA, no significant Y-chromosomal sequence diversity has been detected. We used high throughput sequencing technology to identify the first polymorphic Y-chromosomal markers useful for tracing paternal lines. The nucleotide variability of the modern horse Y chromosome is extremely low, resulting in six haplotypes (HT), all clearly distinct from the Przewalski horse (E. przewalskii). The most widespread HT1 is ancestral and the other five haplotypes apparently arose on the background of HT1 by mutation or gene conversion after domestication. Two haplotypes (HT2 and HT3) are widely distributed at high frequencies among modern European horse breeds. Using pedigree information, we trace the distribution of Y-haplotype diversity to particular founders. The mutation leading to HT3 occurred in the germline of the famous English Thoroughbred stallion "Eclipse" or his son or grandson and its prevalence demonstrates the influence of this popular paternal line on modern sport horse breeds. The pervasive introgression of Thoroughbred stallions during the last 200 years to refine autochthonous breeds has strongly affected the distribution of Y-chromosomal variation in modern horse breeds and has led to the replacement of autochthonous Y chromosomes. Only a few northern European breeds bear unique variants at high frequencies or fixed within but not shared among breeds. Our Y-chromosomal data complement the well established mtDNA lineages and document the male side of the genetic history of modern horse breeds and breeding practices.

A cis-acting regulatory mutation causes premature hair graying and susceptibility to melanoma in the horse.

In horses, graying with age is an autosomal dominant trait associated with a high incidence of melanoma and vitiligo-like depigmentation. Here we show that the Gray phenotype is caused by a 4.6-kb duplication in intron 6 of STX17 (syntaxin-17) that constitutes a cis-acting regulatory mutation. Both STX17 and the neighboring NR4A3 gene are overexpressed in melanomas from Gray horses. Gray horses carrying a loss-of-function mutation in ASIP (agouti signaling protein) had a higher incidence of melanoma, implying that increased melanocortin-1 receptor signaling promotes melanoma development in Gray horses. The Gray horse provides a notable example of how humans have cherry-picked mutations with favorable phenotypic effects in domestic animals.

Comparative histopathology of grey-horse-melanoma and human malignant melanoma.

Equine melanoma shows striking features particularly with regard to clinical development in grey horses: in contrast to malignant melanoma in humans and in solid coloured horses that are characterized by early onset of metastasis, pigment cell tumours display almost benign clinical features in ageing grey horses. Through evolution, grey horses appear to be in a favourable position in regard to the biological behaviour of melanomas. Yet unknown factors inhibiting or retarding early melanoma metastasis may be responsible for this phenomenon. In this study, immunostaining profiles and histopathologic patterns of equine vs. human melanotic tumours were compared. In addition, the expression of melanoma markers currently used in human melanoma detection and characterization were evaluated for their applicability in equine melanoma diagnosis. Immunohistopathologic investigations revealed that benign grey horse melanomas share common features with human blue nevi and with human malignant desmoplastic melanomas, whereas their resemblance to other types of human cutaneous malignant melanomas is less pronounced. Our data equally underline that S-100, proliferating cell nuclear antigen (PCNA), HMB-45, Ki-67, T-311 and CD44 can serve as reliable markers for horse melanomas. Further investigations aiming at identifying factors retarding metastasis in affected grey horses are needed, as they may contribute to the development of novel treatment strategies for human malignant melanoma.