Search

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

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

Abstract 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

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

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

Abstract 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

Abstract 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

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.

Abstract Background Copy Number Variation (CNV) is a common form of genetic variation underlying animal evolution and phenotypic diversity across a wide range of species. In the mammalian genome, high frequency of CNV differentiation between breeds may be candidates for population-specific selection. However, CNV differentiation, selection and its population genetics have been poorly explored in horses. Results We investigated the patterns, population variation and gene annotation of CNV using the Axiom® Equine Genotyping Array (670,796 SNPs) from a large cohort of individuals (N = 1755) belonging to eight European horse breeds, varying from draught horses to several warmblood populations. After quality control, 152,640 SNP CNVs (individual markers), 18,800 segment CNVs (consecutive SNP CNVs of same gain/loss state or both) and 939 CNV regions (CNVRs; overlapping segment CNVs by at least 1 bp) compared to the average signal of the reference (Belgian draught horse) were identified. Our analyses showed that Equus caballus chromosome 12 (ECA12) was the most enriched in segment CNV gains and losses (~ 3% average proportion of the genome covered), but the highest number of segment CNVs were detected on ECA1 and ECA20 (regardless of size). The Friesian horses showed private SNP CNV gains (> 20% of the samples) on ECA1 and Exmoor ponies displayed private SNP CNV losses on ECA25 (> 20% of the samples). The Warmblood cluster showed private SNP CNV gains located in ECA9 and Draught cluster showed private SNP CNV losses located in ECA7. The length of the CNVRs ranged from 1 kb to 21.3 Mb. A total of 10,612 genes were annotated within the CNVRs. The PANTHER annotation of these genes showed significantly under- and overrepresented gene ontology biological terms related to cellular processes and immunity (Bonferroni P-value

Abstract Background A growing demand for improved physical skills and mental attitude in modern sport horses has led to strong selection for performance in many warmblood studbooks. The aim of this study was to detect genomic regions with low diversity, and therefore potentially under selection, in Swedish Warmblood horses (SWB) by analysing high-density SNP data. To investigate if such signatures could be the result of selection for equestrian sport performance, we compared our SWB SNP data with those from Exmoor ponies, a horse breed not selected for sport performance traits. Results The genomic scan for homozygous regions identified long runs of homozygosity (ROH) shared by more than 85% of the genotyped SWB individuals. Such ROH were located on ECA4, ECA6, ECA7, ECA10 and ECA17. Long ROH were instead distributed evenly across the genome of Exmoor ponies in 77% of the chromosomes. Two population differentiation tests (FST and XP-EHH) revealed signatures of selection on ECA1, ECA4, and ECA6 in SWB horses. Conclusions Genes related to behaviour, physical abilities and fertility, appear to be targets of selection in the SWB breed. This study provides a genome-wide map of selection signatures in SWB horses, and ground for further functional studies to unravel the biological mechanisms behind complex traits in horses.

Abstract Background Since the 1950s, the Norwegian–Swedish Coldblooded trotter (NSCT) has been intensively selected for harness racing performance. As a result, the racing performance of the NSCT has improved remarkably; however, this improved racing performance has also been accompanied by a gradual increase in inbreeding level. Inbreeding in NSCT has historically been monitored by using traditional methods that are based on pedigree analysis, but with recent advancements in genomics, the NSCT industry has shown interest in adopting molecular approaches for the selection and maintenance of this breed. Consequently, the aims of the current study were to estimate genomic-based inbreeding coefficients, i.e. the proportion of runs of homozygosity (ROH), for a sample of NSCT individuals using high-density genotyping array data, and subsequently to compare the resulting rate of genomic-based F (FROH) to that of pedigree-based F (FPED) coefficients within the breed. Results A total of 566 raced NSCT were available for analyses. Average FROH ranged from 1.78 to 13.95%. Correlations between FROH and FPED were significant (P

Abstract Background Horses have been strongly selected for speed, strength, and endurance-exercise traits since the onset of domestication. As a result, highly specialized horse breeds have developed with many modern horse breeds often representing closed populations with high phenotypic and genetic uniformity. However, a great deal of variation still exists between breeds, making the horse particularly well suited for genetic studies of athleticism. To identify genomic regions associated with athleticism as it pertains to trotting racing ability in the horse, the current study applies a pooled sequence analysis approach using a unique Nordic horse model. Results Pooled sequence data from three Nordic horse populations were used for FST analysis. After strict filtering, FST analysis yielded 580 differentiated regions for trotting racing ability. Candidate regions on equine chromosomes 7 and 11 contained the largest number of SNPs (n = 214 and 147, respectively). GO analyses identified multiple genes related to intelligence, energy metabolism, and skeletal development as potential candidate genes. However, only one candidate region for trotting racing ability overlapped a known racing ability QTL. Conclusions Not unexpected for genomic investigations of complex traits, the current study identified hundreds of candidate regions contributing to trotting racing ability in the horse. Likely resulting from the cumulative effects of many variants across the genome, racing ability continues to demonstrate its polygenic nature with candidate regions implicating genes influencing both musculature and neurological development.

Abstract Background Although harness racing is of high economic importance to the global equine industry, significant genomic resources have yet to be applied to mapping harness racing success. To identify genomic regions associated with harness racing success, the current study performs genome-wide association analyses with three racing performance traits in the Norwegian-Swedish Coldblooded Trotter using the 670 K Axiom Equine Genotyping Array. Results Following quality control, 613 horses and 359,635 SNPs were retained for further analysis. After strict Bonferroni correction, nine genome-wide significant SNPs were identified for career earnings. No genome-wide significant SNPs were identified for number of gallops or best km time. However, four suggestive genome-wide significant SNPs were identified for number of gallops, while 19 were identified for best km time. Multiple genes related to intelligence, energy metabolism, and immune function were identified as potential candidate genes for harness racing success. Conclusions Apart from the physiological requirements needed for a harness racing horse to be successful, the results of the current study also advocate learning ability and memory as important elements for harness racing success. Further exploration into the mental capacity required for a horse to achieve racing success is likely warranted.

Abstract Background Many common and relevant diseases affecting equine welfare have yet to be tested regarding structural variants such as copy number variations (CNVs). CNVs make up a substantial proportion of total genetic variability in populations of many species, resulting in more sequence differences between individuals than SNPs. Associations between CNVs and disease phenotypes have been established in several species, but equine CNV studies have been limited. Aim of this study was to identify CNVs and to perform a genome-wide association (GWA) study in Friesian horses to identify genomic loci associated with insect bite hypersensitivity (IBH), a common seasonal allergic dermatitis observed in many horse breeds worldwide. Results Genotypes were obtained using the Axiom® Equine Genotyping Array containing 670,796 SNPs. After quality control of genotypes, 15,041 CNVs and 5350 CNV regions (CNVRs) were identified in 222 Friesian horses. Coverage of the total genome by CNVRs was 11.2% with 49.2% of CNVRs containing genes. 58.0% of CNVRs were novel (i.e. so far only identified in Friesian horses). A SNP- and CNV-based GWA analysis was performed, where about half of the horses were affected by IBH. The SNP-based analysis showed a highly significant association between the MHC region on ECA20 and IBH in Friesian horses. Associations between the MHC region on ECA20 and IBH were also detected based on the CNV-based analysis. However, CNVs associated with IBH in Friesian horses were not often in close proximity to SNPs identified to be associated with IBH. Conclusions CNVs were identified in a large sample of the Friesian horse population, thereby contributing to our knowledge on CNVs in horses and facilitating our understanding of the equine genome and its phenotypic expression. A clear association was identified between the MHC region on ECA20 and IBH in Friesian horses based on both SNP- and CNV-based GWA studies. These results imply that MHC contributes to IBH sensitivity in Friesian horses. Although subsequent analyses are needed for verification, nucleotide differences, as well as more complex structural variations like CNVs, seem to contribute to IBH sensitivity. IBH should be considered as a common disease with a complex genomic architecture.

Abstract Background Miniature size in horses represents an extreme reduction of withers height that originated after domestication. In some breeds, it is a highly desired trait representing a breed- or subtype-specific feature. The genomic changes that emerged due to strong-targeted selection towards this distinct type remain unclear. Results Comparisons of whole-genome sequencing data from two Miniature Shetland ponies and one standard-sized Shetland pony, performed to elucidate genetic determinants for miniature size, revealed four synergistic variants, limiting withers height to 34.25 in. (87 cm). Runs of homozygosity regions were detected spanning these four variants in both the Miniature Shetland ponies and the standard-sized Shetland pony. They were shown to be characteristic of the Shetland pony breed, resulting in a miniature type under specific genotypic combinations. These four genetic variants explained 72% of the size variation among Shetland ponies and related breeds. The length of the homozygous regions indicate that they arose over 1000 years ago. In addition, a copy number variant was identified in DIAPH3 harboring a loss exclusively in ponies and donkeys and thus representing a potential height-associated variant. Conclusion This study reveals main drivers for miniature size in horses identified in whole genome data and thus provides relevant candidate genes for extremely short stature in mammals.

Abstract Background To date, genome-scale analyses in the domestic horse have been limited by suboptimal single nucleotide polymorphism (SNP) density and uneven genomic coverage of the current SNP genotyping arrays. The recent availability of whole genome sequences has created the opportunity to develop a next generation, high-density equine SNP array. Results Using whole genome sequence from 153 individuals representing 24 distinct breeds collated by the equine genomics community, we cataloged over 23 million de novo discovered genetic variants. Leveraging genotype data from individuals with both whole genome sequence, and genotypes from lower-density, legacy SNP arrays, a subset of ~5 million high-quality, high-density array candidate SNPs were selected based on breed representation and uniform spacing across the genome. Considering probe design recommendations from a commercial vendor (Affymetrix, now Thermo Fisher Scientific) a set of ~2 million SNPs were selected for a next-generation high-density SNP chip (MNEc2M). Genotype data were generated using the MNEc2M array from a cohort of 332 horses from 20 breeds and a lower-density array, consisting of ~670 thousand SNPs (MNEc670k), was designed for genotype imputation. Conclusions Here, we document the steps taken to design both the MNEc2M and MNEc670k arrays, report genomic and technical properties of these genotyping platforms, and demonstrate the imputation capabilities of these tools for the domestic horse.

Abstract Background The syndrome Multiple Congenital Ocular Anomalies (MCOA) is a congenital eye disorder in horses. Both the MCOA syndrome and the Silver coat colour in horses are caused by the same missense mutation in the premelanosome protein (PMEL) gene. Horses homozygous for the Silver mutation (TT) are affected by multiple ocular defects causing visual impairment or blindness. Horses heterozygous for the Silver mutation (CT) have less severe clinical signs, usually cysts arising from the ciliary body iris or retina temporally. It is still unknown if the vision is impaired in horses heterozygous for the Silver mutation. A recent study reported that Comtois horses carrying the Silver mutation had significantly deeper anterior chambers of the eye compared to wild-type horses. This could potentially cause refractive errors. The purpose of the present study was to investigate if Icelandic horses with the Silver mutation have refractive errors compared to wild-type horses. One hundred and fifty-two Icelandic horses were included in the study, 71 CT horses and five TT horses. All horses were genotyped for the missense mutation in PMEL. Each CT and TT horse was matched by a wild-type (CC) horse of the same age ± 1 year. Skiascopy and a brief ophthalmic examination were performed in all horses. Association between refraction and age, eye, genotype and sex was tested by linear mixed-effect model analysis. TT horses with controls were not included in the statistical analyses as they were too few. Results The interaction between age and genotype had a significant impact on the refractive state (P = 0.0001). CT horses older than 16 years were on average more myopic than wild-type horses of the same age. No difference in the refractive state could be observed between genotypes (CT and CC) in horses younger than 16 years. TT horses were myopic (−2 D or more) in one or both eyes regardless of age. Conclusion Our results indicate that an elderly Icelandic horse (older than 16 years) carrying the Silver mutation is more likely to be myopic than a wild-type horse of the same age.

[This corrects the article DOI: 10.1371/journal.pone.0202584.].

Skeletal atavism in Shetland ponies is a heritable disorder characterized by abnormal growth of the ulna and fibula that extend the carpal and tarsal joints, respectively. This causes abnormal skeletal structure and impaired movements, and affected foals are usually killed. In order to identify the causal mutation we subjected six confirmed Swedish cases and a DNA pool consisting of 21 control individuals to whole genome resequencing. We screened for polymorphisms where the cases and the control pool were fixed for opposite alleles and observed this signature for only 25 SNPs, most of which were scattered on genome assembly unassigned scaffolds. Read depth analysis at these loci revealed homozygosity or compound heterozygosity for two partially overlapping large deletions in the pseudoautosomal region (PAR) of chromosome X/Y in cases but not in the control pool. One of these deletions removes the entire coding region of the SHOX gene and both deletions remove parts of the CRLF2 gene located downstream of SHOX. The horse reference assembly of the PAR is highly fragmented, and in order to characterize this region we sequenced bacterial artificial chromosome (BAC) clones by single-molecule real-time (SMRT) sequencing technology. This considerably improved the assembly and enabled size estimations of the two deletions to 160−180 kb and 60−80 kb, respectively. Complete association between the presence of these deletions and disease status was verified in eight other affected horses. The result of the present study is consistent with previous studies in humans showing crucial importance of SHOX for normal skeletal development.

The Colombian paso horse, the most important horse breed in Colombia, performs specific and particular gaits (paso fino, trocha, and Colombian trot), which display different footfall patterns and stride frequencies. The breed has been selected for gait and conformation for more than 50 years and we hypothesize that this selection has led to kinematic differences of the gaits that can be explained by different genetic variants. Hence, the aims of the study were: 1. To identify if there are any differences in the kinematic and genetic variants between the Colombian paso horse's gaits. 2. To evaluate if and how much the gait differences were explained by the nonsense mutation in the DMRT3 gene and 3. To evaluate these results for selecting and controlling the horses gait performance. To test our hypotheses, kinematic data, microsatellites and DMRT3 genotypes for 187 Colombian paso horses were analyzed. The results indicated that there are significant kinematic and DMRT3 differences between the Colombian paso horse's gaits, and those parameters can be used partially to select and control the horses gait performance. However, the DMRT3 gene does not play a major role in controlling the trocha and the Colombian trot gaits. Therefore, modifying genes likely influence these gaits. This study may serve as a foundation for implementing a genetic selection program in the Colombian paso horse and future gene discovery studies for locomotion pattern in horses.

The Swedish-Norwegian Coldblooded trotter (CBT) is a local breed in Sweden and Norway mainly used for harness racing. Previous studies have shown that a mutation from cytosine (C) to adenine (A) in the doublesex and mab-3 related transcription factor 3 (DMRT3) gene has a major impact on harness racing performance of different breeds. An association of the DMRT3 mutation with early career performance has also been suggested. The aim of the current study was to investigate this proposed association in a randomly selected group of CBTs. 769 CBTs (485 raced, 284 unraced) were genotyped for the DMRT3 mutation. The association with racing performance was investigated for 13 performance traits and three different age intervals: 3 years, 3 to 6 years, and 7 to 10 years of age, using the statistical software R. Each performance trait was analyzed for association with DMRT3 using linear models. The results suggest no association of the DMRT3 mutation with precocity (i.e. performance at 3 years of age). Only two traits (race time and number of disqualifications) were significantly different between the genotypes, with AA horses having the fastest times and CC horses having the highest number of disqualifications at 3 years of age. The frequency of the AA genotype was significantly lower in the raced CBT sample compared with the unraced sample and less than 50% of the AA horses participated in a race. For the age intervals 3 to 6 and 7 to 10 years the AA horses also failed to demonstrate significantly better performance than the other genotypes. Although suggested as the most favorable genotype for racing performance in Standardbreds and Finnhorses across all ages, the AA genotype does not appear to be associated with superior performance, early or late, in the racing career of CBTs.

While susceptibility to hypersensitive reactions is a common problem amongst humans and animals alike, the population structure of certain animal species and breeds provides a more advantageous route to better understanding the biology underpinning these conditions. The current study uses Exmoor ponies, a highly inbred breed of horse known to frequently suffer from insect bite hypersensitivity, to identify genomic regions associated with a type I and type IV hypersensitive reaction. A total of 110 cases and 170 controls were genotyped on the 670K Axiom Equine Genotyping Array. Quality control resulted in 452,457 SNPs and 268 individuals being tested for association. Genome-wide association analyses were performed using the GenABEL package in R and resulted in the identification of two regions of interest on Chromosome 8. The first region contained the most significant SNP identified, which was located in an intron of the DCC netrin 1 receptor gene. The second region identified contained multiple top SNPs and encompassed the PIGN, KIAA1468, TNFRSF11A, ZCCHC2, and PHLPP1 genes. Although additional studies will be needed to validate the importance of these regions in horses and the relevance of these regions in other species, the knowledge gained from the current study has the potential to be a step forward in unraveling the complex nature of hypersensitive reactions.

Horses were domesticated from the Eurasian steppes 5,000-6,000 years ago. Since then, the use of horses for transportation, warfare, and agriculture, as well as selection for desired traits and fitness, has resulted in diverse populations distributed across the world, many of which have become or are in the process of becoming formally organized into closed, breeding populations (breeds). This report describes the use of a genome-wide set of autosomal SNPs and 814 horses from 36 breeds to provide the first detailed description of equine breed diversity. F(ST) calculations, parsimony, and distance analysis demonstrated relationships among the breeds that largely reflect geographic origins and known breed histories. Low levels of population divergence were observed between breeds that are relatively early on in the process of breed development, and between those with high levels of within-breed diversity, whether due to large population size, ongoing outcrossing, or large within-breed phenotypic diversity. Populations with low within-breed diversity included those which have experienced population bottlenecks, have been under intense selective pressure, or are closed populations with long breed histories. These results provide new insights into the relationships among and the diversity within breeds of horses. In addition these results will facilitate future genome-wide association studies and investigations into genomic targets of selection.

Equine Multiple Congenital Ocular Anomalies (MCOA) syndrome is a heritable eye disorder mainly affecting silver colored horses. Clinically, the disease manifests in two distinct classes depending on the horse genotype. Horses homozygous for the mutant allele present with a wide range of ocular defects, such as iris stromal hypoplasia, abnormal pectinate ligaments, megaloglobus, iridociliary cysts and cataracts. The phenotype of heterozygous horses is less severe and predominantly includes iridociliary cysts, which occasionally extend into the temporal retina. In order to determine the genetic cause of MCOA syndrome we sequenced the entire previously characterized 208 kilobase region on chromosome 6 in ten individuals; five MCOA affected horses from three different breeds, one horse with the intermediate Cyst phenotype and four unaffected controls from two different breeds. This was performed using Illumina TruSeq technology with paired-end reads. Through the systematic exclusion of all polymorphisms barring two SNPs in PMEL, a missense mutation previously reported to be associated with the silver coat colour and a non-conserved intronic SNP, we establish that this gene is responsible for MCOA syndrome. Our finding, together with recent advances that show aberrant protein function due to the coding mutation, suggests that the missense mutation is causative and has pleiotrophic effect, causing both the horse silver coat color and MCOA syndrome.

Several studies based on a variety of genetic markers have attempted to establish the origins of horse domestication. Thus far a discrepancy between the results of mitochondrial DNA analysis, which show high levels of diversity, and results from the Y-chromosome, with almost no genetic variability, has been identified. Most previous work on the horse Y-chromosome has focused on widespread, popular breeds or local Asian breeds. It is possible that these breeds represent a reduced set of the genetic variation present in the species. Additional genetic variation may be present in local breeds and ancient feral populations, such as the Retuertas horse in Spain. In this study we analyzed the Y-chromosome of the Retuertas horse, a feral horse population on the Iberian Peninsula that is at least several hundred years old, and whose genetic diversity and morphology suggests that it has been reproductively isolated for a long time. Data from the Retuertas horse was compared to another 11 breeds from the region (Portugal, Spain and France) or likely of Iberian origin, and then to data from 15 more breeds from around the globe. We sequenced 31 introns, Zinc finger Y-chromosomal protein (ZFY) and anonymous Y-linked fragments and genotyped 6 microsatellite loci found on the Y-chromosome. We found no sequence variation among all individuals and all breeds studied. However, fifteen differences were discovered between our data set and reference sequences in GenBank. We show that these likely represent errors within the deposited sequences, and suggest that they should not be used as comparative data for future projects.

Intense selective pressures applied over short evolutionary time have resulted in homogeneity within, but substantial variation among, horse breeds. Utilizing this population structure, 744 individuals from 33 breeds, and a 54,000 SNP genotyping array, breed-specific targets of selection were identified using an F(ST)-based statistic calculated in 500-kb windows across the genome. A 5.5-Mb region of ECA18, in which the myostatin (MSTN) gene was centered, contained the highest signature of selection in both the Paint and Quarter Horse. Gene sequencing and histological analysis of gluteal muscle biopsies showed a promoter variant and intronic SNP of MSTN were each significantly associated with higher Type 2B and lower Type 1 muscle fiber proportions in the Quarter Horse, demonstrating a functional consequence of selection at this locus. Signatures of selection on ECA23 in all gaited breeds in the sample led to the identification of a shared, 186-kb haplotype including two doublesex related mab transcription factor genes (DMRT2 and 3). The recent identification of a DMRT3 mutation within this haplotype, which appears necessary for the ability to perform alternative gaits, provides further evidence for selection at this locus. Finally, putative loci for the determination of size were identified in the draft breeds and the Miniature horse on ECA11, as well as when signatures of selection surrounding candidate genes at other loci were examined. This work provides further evidence of the importance of MSTN in racing breeds, provides strong evidence for selection upon gait and size, and illustrates the potential for population-based techniques to find genomic regions driving important phenotypes in the modern horse.

An equine SNP genotyping array was developed and evaluated on a panel of samples representing 14 domestic horse breeds and 18 evolutionarily related species. More than 54,000 polymorphic SNPs provided an average inter-SNP spacing of ∼43 kb. The mean minor allele frequency across domestic horse breeds was 0.23, and the number of polymorphic SNPs within breeds ranged from 43,287 to 52,085. Genome-wide linkage disequilibrium (LD) in most breeds declined rapidly over the first 50-100 kb and reached background levels within 1-2 Mb. The extent of LD and the level of inbreeding were highest in the Thoroughbred and lowest in the Mongolian and Quarter Horse. Multidimensional scaling (MDS) analyses demonstrated the tight grouping of individuals within most breeds, close proximity of related breeds, and less tight grouping in admixed breeds. The close relationship between the Przewalski's Horse and the domestic horse was demonstrated by pair-wise genetic distance and MDS. Genotyping of other Perissodactyla (zebras, asses, tapirs, and rhinoceros) was variably successful, with call rates and the number of polymorphic loci varying across taxa. Parsimony analysis placed the modern horse as sister taxa to Equus przewalski. The utility of the SNP array in genome-wide association was confirmed by mapping the known recessive chestnut coat color locus (MC1R) and defining a conserved haplotype of ∼750 kb across all breeds. These results demonstrate the high quality of this SNP genotyping resource, its usefulness in diverse genome analyses of the horse, and potential use in related species.