Genetic diversity and signatures of selection in Icelandic horses and Exmoor ponies.

BACKGROUND: The Icelandic horse and Exmoor pony are ancient, native breeds, adapted to harsh environmental conditions and they have both undergone severe historic bottlenecks. However, in modern days, the selection pressures on these breeds differ substantially. The aim of this study was to assess genetic diversity in both breeds through expected (HE) and observed heterozygosity (HO) and effective population size (Ne). Furthermore, we aimed to identify runs of homozygosity (ROH) to estimate and compare genomic inbreeding and signatures of selection in the breeds. RESULTS: HO was estimated at 0.34 and 0.33 in the Icelandic horse and Exmoor pony, respectively, aligning closely with HE of 0.34 for both breeds. Based on genomic data, the Ne for the last generation was calculated to be 125 individuals for Icelandic horses and 42 for Exmoor ponies. Genomic inbreeding coefficient (FROH) ranged from 0.08 to 0.20 for the Icelandic horse and 0.12 to 0.27 for the Exmoor pony, with the majority of inbreeding attributed to short ROHs in both breeds. Several ROH islands associated with performance were identified in the Icelandic horse, featuring target genes such as DMRT3, DOCK8, EDNRB, SLAIN1, and NEURL1. Shared ROH islands between both breeds were linked to metabolic processes (FOXO1), body size, and the immune system (CYRIB), while private ROH islands in Exmoor ponies were associated with coat colours (ASIP, TBX3, OCA2), immune system (LYG1, LYG2), and fertility (TEX14, SPO11, ADAM20). CONCLUSIONS: Evaluations of genetic diversity and inbreeding reveal insights into the evolutionary trajectories of both breeds, highlighting the consequences of population bottlenecks. While the genetic diversity in the Icelandic horse is acceptable, a critically low genetic diversity was estimated for the Exmoor pony, which requires further validation. Identified signatures of selection highlight the differences in the use of the two breeds as well as their adaptive trait similarities. The results provide insight into genomic regions under selection pressure in a gaited performance horse breed and various adaptive traits in small-sized native horse breeds. This understanding contributes to preserving genetic diversity and population health in these equine populations.

An intronic copy number variation in Syntaxin 17 determines speed of greying and melanoma incidence in Grey horses.

The Greying with age phenotype in horses involves loss of hair pigmentation whereas skin pigmentation is not reduced, and a predisposition to melanoma. The causal mutation was initially reported as a duplication of a 4.6 kb intronic sequence in Syntaxin 17. The speed of greying varies considerably among Grey horses. Here we demonstrate the presence of two different Grey alleles, G2 carrying two tandem copies of the duplicated sequence and G3 carrying three. The latter is by far the most common allele, probably due to strong selection for the striking white phenotype. Our results reveal a remarkable dosage effect where the G3 allele is associated with fast greying and high incidence of melanoma whereas G2 is associated with slow greying and low incidence of melanoma. The copy number expansion transforms a weak enhancer to a strong melanocyte-specific enhancer that underlies hair greying (G2 and G3) and a drastically elevated risk of melanoma (G3 only). Our direct pedigree-based observation of the origin of a G2 allele from a G3 allele by copy number contraction demonstrates the dynamic evolution of this locus and provides the ultimate evidence for causality of the copy number variation of the 4.6 kb intronic sequence.

Effect of an endothelial regulatory module on plasma proteomics in exercising horses.

Elite performing exercise requires an intricate modulation of the blood pressure to support the working muscles with oxygen. We have previously identified a genomic regulatory module that associates with differences in blood pressures of importance for elite performance in racehorses. This study aimed to determine the effect of the regulatory module on the protein repertoire. We sampled plasma from 12 Coldblooded trotters divided into two endothelial regulatory module haplotype groups, a sub-elite performing haplotype (SPH) and an elite performing haplotype (EPH), each at rest and exercise. The haplotype groups and their interaction were interrogated in two analyses, i) individual paired ratio analysis for identifying differentially abundant proteins of exercise (DAPE) and interaction (DAPI) between haplotype and exercise, and ii) unpaired ratio analysis for identifying differentially abundant protein of haplotype (DAPH). The proteomics analyses revealed a widespread change in plasma protein content during exercise, with a decreased tendency in protein abundance that is mainly related to lung function, tissue fluids, metabolism, calcium ion pathway and cellular energy metabolism. Furthermore, we provide the first investigation of the proteome variation due to the interaction between exercise and related blood pressure haplotypes, which this difference was related to a faster switch to the lipoprotein and lipid metabolism during exercise for EPH. The molecular signatures identified in the present study contribute to an improved understanding of exercise-related blood pressure regulation.

An endothelial regulatory module links blood pressure regulation with elite athletic performance.

The control of transcription is crucial for homeostasis in mammals. A previous selective sweep analysis of horse racing performance revealed a 19.6 kb candidate regulatory region 50 kb downstream of the Endothelin3 (EDN3) gene. Here, the region was narrowed to a 5.5 kb span of 14 SNVs, with elite and sub-elite haplotypes analyzed for association to racing performance, blood pressure and plasma levels of EDN3 in Coldblooded trotters and Standardbreds. Comparative analysis of human HiCap data identified the span as an enhancer cluster active in endothelial cells, interacting with genes relevant to blood pressure regulation. Coldblooded trotters with the sub-elite haplotype had significantly higher blood pressure compared to horses with the elite performing haplotype during exercise. Alleles within the elite haplotype were part of the standing variation in pre-domestication horses, and have risen in frequency during the era of breed development and selection. These results advance our understanding of the molecular genetics of athletic performance and vascular traits in both horses and humans.

The genetics of gaits in Icelandic horses goes beyond DMRT3, with RELN and STAU2 identified as two new candidate genes.

BACKGROUND: In domesticated animals, many important traits are complex and regulated by a large number of genes, genetic interactions, and environmental influences. The ability of Icelandic horses to perform the gait 'pace' is largely influenced by a single mutation in the DMRT3 gene, but genetic modifiers likely exist. The aim of this study was to identify novel genetic factors that influence pacing ability and quality of the gait through a genome-wide association study (GWAS) and correlate new findings to previously identified quantitative trait loci (QTL) and mutations. RESULTS: Three hundred and seventy-two Icelandic horses were genotyped with the 670 K+ Axiom Equine Genotyping Array, of which 362 had gait scores from breeding field tests. A GWAS revealed several SNPs on Equus caballus chromosomes (ECA) 4, 9, and 20 that were associated (p < 1.0 × 10-5) with the breeding field test score for pace. The two novel QTL on ECA4 and 9 were located within the RELN and STAU2 genes, respectively, which have previously been associated with locomotor behavior in mice. Haplotypes were identified and the most frequent one for each of these two QTL had a large favorable effect on pace score. The second most frequent haplotype for the RELN gene was positively correlated with scores for tölt, trot, gallop, and canter. Similarly, the second most frequent haplotype for the STAU2 gene had favorable effects on scores for trot and gallop. Different genotype ratios of the haplotypes in the RELN and STAU2 genes were also observed in groups of horses with different levels of pacing ability. Furthermore, interactions (p < 0.05) were detected for the QTL in the RELN and STAU2 genes with the DMRT3 gene. The novel QTL on ECA4, 9, and 20, along with the effects of the DMRT3 variant, were estimated to account jointly for 27.4% of the phenotypic variance of the gait pace. CONCLUSIONS: Our findings provide valuable information about the genetic architecture of pace beyond the contribution of the DMRT3 gene and indicate genetic interactions that contribute to the complexity of this trait. Further investigation is needed to fully understand the underlying genetic factors and interactions.

Genetic diversity and signatures of selection in four indigenous horse breeds of Iran.

Indigenous Iranian horse breeds were evolutionarily affected by natural and artificial selection in distinct phylogeographic clades, which shaped their genomes in several unique ways. The aims of this study were to evaluate the genetic diversity and genomewide selection signatures in four indigenous Iranian horse breeds. We evaluated 169 horses from Caspian (n = 21), Turkmen (n = 29), Kurdish (n = 67), and Persian Arabian (n = 52) populations, using genomewide genotyping data. The contemporary effective population sizes were 59, 98, 102, and 113 for Turkmen, Caspian, Persian Arabian, and Kurdish breeds, respectively. By analysis of the population genetic structure, we classified the north breeds (Caspian and Turkmen) and west/southwest breeds (Persian Arabian and Kurdish) into two phylogeographic clades reflecting their geographic origin. Using the de-correlated composite of multiple selection signal statistics based on pairwise comparisons, we detected a different number of significant SNPs under putative selection from 13 to 28 for the six pairwise comparisons (FDR < 0.05). The identified SNPs under putative selection coincided with genes previously associated with known QTLs for morphological, adaptation, and fitness traits. Our results showed HMGA2 and LLPH as strong candidate genes for height variation between Caspian horses with a small size and the other studied breeds with a medium size. Using the results of studies on human height retrieved from the GWAS catalog, we suggested 38 new putative candidate genes under selection. These results provide a genomewide map of selection signatures in the studied breeds, which represent valuable information for formulating genetic conservation and improved breeding strategies for the breeds.

Refining the evolutionary tree of the horse Y chromosome.

The Y chromosome carries information about the demography of paternal lineages, and thus, can prove invaluable for retracing both the evolutionary trajectory of wild animals and the breeding history of domesticates. In horses, the Y chromosome shows a limited, but highly informative, sequence diversity, supporting the increasing breeding influence of Oriental lineages during the last 1500 years. Here, we augment the primary horse Y-phylogeny, which is currently mainly based on modern horse breeds of economic interest, with haplotypes (HT) segregating in remote horse populations around the world. We analyze target enriched sequencing data of 5 Mb of the Y chromosome from 76 domestic males, together with 89 whole genome sequenced domestic males and five Przewalski's horses from previous studies. The resulting phylogeny comprises 153 HTs defined by 2966 variants and offers unprecedented resolution into the history of horse paternal lineages. It reveals the presence of a remarkable number of previously unknown haplogroups in Mongolian horses and insular populations. Phylogenetic placement of HTs retrieved from 163 archaeological specimens further indicates that most of the present-day Y-chromosomal variation evolved after the domestication process that started around 4200 years ago in the Western Eurasian steppes. Our comprehensive phylogeny significantly reduces ascertainment bias and constitutes a robust evolutionary framework for analyzing horse population dynamics and diversity.

Genome-Wide Association Analyses of Osteochondrosis in Belgian Warmbloods Reveal Candidate Genes Associated With Chondrocyte Development.

Osteochondrosis (OC) is an important skeletal disease causing profound welfare concerns in horses. Although numerous studies have explored the genetics underlying OC in various breeds, the Belgian Warmblood (BW) remains unstudied despite having a concerning prevalence of 32.0%. As a result, this study aimed to conduct genome-wide association (GWA) analyses to identify candidate variants associated with OC in BWs. To achieve this, blood samples and radiographs were collected from 407 Belgian Warmbloods registered to one of two BW studbooks (Belgisch Warmbloedpaard and Zangersheide), and genotyping was performed using the 670K Axiom Equine Genotyping Array. GWA analyses using a principle component approach were then performed on OC status (OCS; presence or absence of OC at any joint), hock OC status (HOC) and stifle OC status (SOC). These analyses yielded significantly associated (P < .01) SNPs on Equus caballus chromosome (ECA) 3, ECA 12, and ECA 18 for OCS; however, no single nucleotide polymorphisms (SNPs) reached significance for HOC or SOC. Subsequent analysis of candidate genes within 500 kilobases of the significant SNPs revealed functions broadly related to cell differentiation and chondrocyte development. While this study represents another step forward in uncovering variants and biological pathways associated with OC, additional studies are needed to validate the newly identified candidate SNPs for OC in BWs. Further studies of OC in BWs, as well as other breeds, are critical in our efforts to fully understand the disease's etiopathogenesis and ultimately provide breeding programs better equipped to improve horse health and well-being.

Genetic diversity and selection in Puerto Rican horses.

Since the first Spanish settlers brought horses to America centuries ago, several local varieties and breeds have been established in the New World. These were generally a consequence of the admixture of the different breeds arriving from Europe. In some instances, local horses have been selectively bred for specific traits, such as appearance, endurance, strength, and gait. We looked at the genetics of two breeds, the Puerto Rican Non-Purebred (PRNPB) (also known as the "Criollo") horses and the Puerto Rican Paso Fino (PRPF), from the Caribbean Island of Puerto Rico. While it is reasonable to assume that there was a historic connection between the two, the genetic link between them has never been established. In our study, we started by looking at the genetic ancestry and diversity of current Puerto Rican horse populations using a 668 bp fragment of the mitochondrial DNA D-loop (HVR1) in 200 horses from 27 locations on the island. We then genotyped all 200 horses in our sample for the "gait-keeper" DMRT3 mutant allele previously associated with the paso gait especially cherished in this island breed. We also genotyped a subset of 24 samples with the Illumina Neogen Equine Community genome-wide array (65,000 SNPs). This data was further combined with the publicly available PRPF genomes from other studies. Our analysis show an undeniable genetic connection between the two varieties in Puerto Rico, consistent with the hypothesis that PRNPB horses represent the descendants of the original genetic pool, a mix of horses imported from the Iberian Peninsula and elsewhere in Europe. Some of the original founders of PRNRB population must have carried the "gait-keeper" DMRT3 allele upon arrival to the island. From this admixture, the desired traits were selected by the local people over the span of centuries. We propose that the frequency of the mutant "gait-keeper" allele originally increased in the local horses due to the selection for the smooth ride and other characters, long before the PRPF breed was established. To support this hypothesis, we demonstrate that PRNPB horses, and not the purebred PRPF, carry a signature of selection in the genomic region containing the DMRT3 locus to this day. The lack of the detectable signature of selection associated with the DMRT3 in the PRPF would be expected if this native breed was originally derived from the genetic pool of PRNPB horses established earlier and most of the founders already had the mutant allele. Consequently, selection specific to PRPF later focused on allels in other genes (including CHRM5, CYP2E1, MYH7, SRSF1, PAM, PRN and others) that have not been previously associated with the prized paso gait phenotype in Puerto Rico or anywhere else.

The Enrichment of Specific Hair Follicle-Associated Cell Populations in Plucked Hairs Offers an Opportunity to Study Gene Expression Underlying Hair Traits.

Gene expression differences can assist in characterizing important underlying genetic mechanisms between different phenotypic traits. However, when population-dense tissues are studied, the signals from scarce populations are diluted. Therefore, appropriately choosing a sample collection method that enriches a particular type of effector cells might yield more specific results. To address this issue, we performed a polyA-selected RNA-seq experiment of domestic horse (Equus ferus caballus) plucked-hair samples and skin biopsies. Then, we layered the horse gene abundance results against cell type-specific marker genes generated from a scRNA-seq supported with spatial mapping of laboratory mouse (Mus musculus) skin to identify the captured populations. The hair-plucking and skin-biopsy sample-collection methods yielded comparable quality and quantity of RNA-seq results. Keratin-related genes, such as KRT84 and KRT75, were among the genes that showed higher abundance in plucked hairs, while genes involved in cellular processes and enzymatic activities, such as MGST1, had higher abundance in skin biopsies. We found an enrichment of hair-follicle keratinocytes in plucked hairs, but detected an enrichment of other populations, including epidermis keratinocytes, in skin biopsies. In mammalian models, biopsies are often the method of choice for a plethora of gene expression studies and to our knowledge, this is a novel study that compares the cell-type enrichment between the non-invasive hair-plucking and the invasive skin-biopsy sample-collection methods. Here, we show that the non-invasive and ethically uncontroversial plucked-hair method is recommended depending on the research question. In conclusion, our study will allow downstream -omics approaches to better understand integumentary conditions in both health and disease in horses as well as other mammals.

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.

A QTL for conformation of back and croup influences lateral gait quality in Icelandic horses.

BACKGROUND: The back plays a vital role in horse locomotion, where the spine functions as a spring during the stride cycle. A complex interaction between the spine and the muscles of the back contribute to locomotion soundness, gait ability, and performance of riding and racehorses. Conformation is commonly used to select horses for breeding and performance in multiple horse breeds, where the back and croup conformation plays a significant role. The conformation of back and croup plays an important role on riding ability in Icelandic horses. However, the genes behind this trait are still unknown. Therefore, the aim of this study was to identify genomic regions associated with conformation of back and croup in Icelandic horses and to investigate their effects on riding ability. One hundred seventy-seven assessed Icelandic horses were included in the study. A genome-wide association analysis was performed using the 670 K+ Axiom Equine Genotyping Array, and the effects of different haplotypes in the top associated region were estimated for riding ability and additional conformation traits assessed during breeding field tests. RESULTS: A suggestive quantitative trait loci (QTL) for the score of back and croup was detected on Equus caballus (ECA) 22 (p-value = 2.67 × 10- 7). Haplotype analysis revealed two opposite haplotypes, which resulted in higher and lower scores of the back and croup, respectively (p-value < 0.001). Horses with the favorable haplotype were more inclined to have a well-balanced backline with an uphill conformation and had, on average, higher scores for the lateral gaits tölt (p-value = 0.02) and pace (p-value = 0.004). This genomic region harbors three genes: C20orf85, ANKRD60 and LOC100056167. ANKRD60 is associated with body height in humans. C20orf85 and ANKRD60 are potentially linked to adolescent idiopathic scoliosis in humans. CONCLUSIONS: Our results show that the detected QTL for conformation of back and croup is of importance for quality of lateral gaits in Icelandic horses. These findings could result in a genetic test to aid in the selection of breeding horses, thus they are of major interest for horse breeders. The results may also offer a gateway to comparative functional genomics by potentially linking both motor laterality and back inclination in horses with scoliosis in humans.

An investigation into factors influencing basal eye temperature in the domestic horse (Equus caballus) when measured using infrared thermography in field conditions.

Infrared thermography (IRT) is a popular technology used for the detection of thermal changes given its non-invasive nature and lack of direct contact with the individual. Accordingly, the maximal eye temperature (MaxET) measured with IRT has been extensively applied in equine research. However, there is little information available about the potential limitations of the MaxET in field studies. Thus, the aims of this study were to 1) quantify the individual variation of MaxET in field conditions and the effects of individual, breed, body size (height at withers), eye side, sex and age, 2) determine the effects of environment and operator, and 3) explore the relationship between MaxET and rectal temperature (RT) at rest. To accomplish these aims, 791 MaxET measures from 32 horses were collected in Sweden in five different months and five farms over a period of 12 months. There was an effect of individual on IRT (P < .05) and individual MaxET varied from 29.4 to 37.6 °C. IRT was also affected (P < .05) by breed and sex (maximal difference 1.1 °C and 0.3 °C, respectively) but not by eye side, age and height at withers. There were significant effects of month and farm (maximal differences; 2.4 and 2.3 °C, respectively), between outdoor and indoor measurements (0.8 °C) and also between operators (0.2 °C). There were no correlations between MaxET and RT. These results demonstrate that in horses observed at rest in their home environment, MaxET is affected by endogenous (sex and breed) and environmental factors (farm, location and month of the year) and shows no relationship to RT. We strongly suggest that IRT technology should be used with great caution in field studies and only under conditions where these factors can be accurately accounted for.

A genome-wide scan for candidate lethal variants in Thoroughbred horses.

Domestic animal populations are often characterised by high rates of inbreeding and low effective population sizes due to selective breeding practices. These practices can result in otherwise rare recessive deleterious alleles drifting to high frequencies, resulting in reduced fertility rates. This study aimed to identify potential recessive lethal haplotypes in the Thoroughbred horse breed, a closed population that has been selectively bred for racing performance. In this study, we identified a haplotype in the LY49B gene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Variant analysis of whole-genome sequence data identified two SNPs in the 3'UTR of the LY49B gene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed that LY49B is expressed in the trophoblast during placentation stage of development. These findings suggest that LY49B may have an essential, but as yet unknown function in the implantation stage of equine development. Further investigation of this region may allow for the development of a genetic test to improve fertility rates in horse populations. Identification of other lethal variants could assist in improving natural levels of fertility in horse populations.

Genetics of Skin Disease in Horses.

Equine skin diseases are common, causing increased costs and reduced welfare of affected horses.Genetic testing, if available, can complement early detection, disease diagnosis, and clinical treatment and offers horse breeders the possibility to rule out carrier status. The mechanisms of complex disease can be investigated by using the latest state-of-the-art genomic technologies. Genome-based strategies may also serve as an efficient and cost-effective strategy for the management of the disease severity levels, with particular interest in complex traits such as insect bite hypersensitivity, chronic progressive lymphedema, and melanoma.

Benefits and risks of barefoot harness racing in Standardbred trotters.

There is a lack of research on the benefits and risks of shoeing conditions in harness racing. Thus, our objectives were to: (a) investigate whether velocity times (VT; s/km) are affected by racing unshod (N = 76,932 records on 5,247 horses); (b) determine the potential risks of galloping, being penalized, and disqualification when competing unshod (N = 111,755 records on 6,423 horses); and (c) identify additional environmental factors that affect VT and risks. VT was found to be significantly influenced by shoeing condition (e.g., unshod, shod front, shod hind, or fully shod), but also by sex, age, season, track, track condition, start method, start position, distance, and driver-horse performance level (p < 2e-16). The risks of galloping and disqualification were significantly influenced by shoeing condition, sex, age, season, track, start method, start position, or driver-horse performance level (p ≤ .05). Horses racing unshod had 0.7 s/km lower VT than fully shod horses and showed better performance when racing on neutral tracks during the late summer than horses with other shoeing conditions during the same period. However, racing unshod increased the relative risks of galloping and disqualification by 15%-35% in all seasons. Horses shod only on the hind hooves showed better performance than fully shod horses, without higher risks associated with competing unshod.

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.

Characterization of A Homozygous Deletion of Steroid Hormone Biosynthesis Genes in Horse Chromosome 29 as A Risk Factor for Disorders of Sex Development and Reproduction.

Disorders of sex development (DSD) and reproduction are not uncommon among horses, though knowledge about their molecular causes is sparse. Here we characterized a ~200 kb homozygous deletion in chromosome 29 at 29.7-29.9 Mb. The region contains AKR1C genes which function as ketosteroid reductases in steroid hormone biosynthesis, including androgens and estrogens. Mutations in AKR1C genes are associated with human DSDs. Deletion boundaries, sequence properties and gene content were studied by PCR and whole genome sequencing of select deletion homozygotes and control animals. Deletion analysis by PCR in 940 horses, including 622 with DSDs and reproductive problems and 318 phenotypically normal controls, detected 67 deletion homozygotes of which 79% were developmentally or reproductively abnormal. Altogether, 8-9% of all abnormal horses were homozygous for the deletion, with the highest incidence (9.4%) among cryptorchids. The deletion was found in ~4% of our phenotypically normal cohort, ~1% of global warmblood horses and ponies, and ~7% of draught breeds of general horse population as retrieved from published data. Based on the abnormal phenotype of the carriers, the functionally relevant gene content, and the low incidence in general population, we consider the deletion in chromosome 29 as a risk factor for equine DSDs and reproductive disorders.

Genome-wide association study for insect bite hypersensitivity susceptibility in horses revealed novel associated loci on chromosome 1.

Equine insect bite hypersensitivity (IBH) is a pruritic skin allergy caused primarily by biting midges, Culicoides spp. IBH susceptibility has polygenic inheritance and occurs at high frequencies in several horse breeds worldwide, causing increased costs and reduced welfare of affected horses. The aim of this study was to identify and validate single nucleotide polymorphisms (SNPs) associated with equine IBH susceptibility. After quality control, 33,523 SNPs were included in a Bayesian genome-wide association study on 177 affected and 178 unaffected Icelandic horses. We report associated regions in E. caballus (ECA) 1, 3, 15 and 18, overlapping with known IBH QTLs in horses, and novel regions containing several genes, together explaining 11.46% of the total genetic variance. For validation, three SNPs on ECA 1 and ECA X (explaining the largest percentage of genetic variance) within 1-mb genomic windows for IBH were genotyped in an independent population of 280 Exmoor ponies. The associated genomic region (152-153 mb) on ECA 1 was confirmed in Exmoor ponies and contains the AQR gene involved in splicing processes and a long non-coding RNA. This study confirms the polygenic nature of IBH susceptibility and suggests a role of transcriptional regulatory mechanisms (e.g., alternative splicing) for IBH predisposition in these horse breeds.