Genetics of swayback in American Saddlebred horses

Predicted genetic burden and frequency of phenotype-associated variants in the horse

Disease-causing variants have been identified for less than 20% of suspected equine genetic diseases. Whole genome sequencing (WGS) allows rapid identification of rare disease causal variants. However, interpreting the clinical variant consequence is confounded by the number of predicted deleterious variants that healthy individuals carry (predicted genetic burden). Estimation of the predicted genetic burden and baseline frequencies of known deleterious or phenotype associated variants within and across the major horse breeds have not been performed. We used WGS of 605 horses across 48 breeds to identify 32,818,945 variants, demonstrate a high predicted genetic burden (median 730 variants/horse, interquartile range: 613-829), show breed differences in predicted genetic burden across 12 target breeds, and estimate the high frequencies of some previously reported disease variants. This large-scale variant catalog for a major and highly athletic domestic animal species will enhance its ability to serve as a model for human phenotypes and improves our ability to discover the bases for important equine phenotypes.

Assessing genetic diversity and defining signatures of positive selection on the genome of dromedary camels from the southeast of the Arabian Peninsula

Dromedary camels (Camelus dromedarius) are members of the Camelini tribe within the Camelidae family. They are distributed throughout North Africa, the Arabian Peninsula and Southeast Asia. This domestic species is characterized by its superior adaptability to the harsh desert environment. In this study, whole autosomal data of 29 dromedary samples from the Southeast Arabian Peninsula in Oman; 10 from Muscat, 14 from Al-Batinah, and 5 from Al-Sharqiya, were investigated to assess their genetic relationship and to define candidate signatures of positive selection. A minimal genetic distinction that separates Muscat dromedaries from the other two populations was observed, with a degree of genetic admixture between them. Using the de-correlated composite of multiple signals (DCMS) approach, a total of 47 candidate regions within the autosomes of these dromedary populations were defined with signatures of positive selection. These candidate regions harbor a total of 154 genes that are mainly associated with functional categories related to immune response, lipid metabolism and energy expenditure, optical and auditory functions, and long-term memory. Different functional genomic variants were called on the candidate regions and respective genes that warrant further investigation to find possible association with the different favorable phenotypes in dromedaries. The output of this study paves the way for further research efforts aimed at defining markers for use in genomic breeding programs, with the goal of conserving the genetic diversity of the species and enhancing its productivity.

New genomic insights into the conformation of Lipizzan horses

Conformation traits are important selection criteria in equine breeding, as they describe the exterior aspects of the horse (height, joint angles, shape). However, the genetic architecture of conformation is not well understood, as data of these traits mainly consist of subjective evaluation scores. Here, we performed genome-wide association studies on two-dimensional shape data of Lipizzan horses. Based on this data, we identified significant quantitative trait loci (QTL) associated with cresty neck on equine chromosome (ECA)16 within the MAGI1 gene, and with type, hereby differentiating heavy from light horses on ECA5 within the POU2F1 gene. Both genes were previously described to affect growth, muscling and fatty deposits in sheep, cattle and pigs. Furthermore, we pin-pointed another suggestive QTL on ECA21, near the PTGER4 gene, associated with human ankylosing spondylitis, for shape differences in the back and pelvis (roach back vs sway back). Further differences in the shape of the back and abdomen were suggestively associated with the RYR1 gene, involved in core muscle weakness in humans. Therefore, we demonstrated that horse shape space data enhance the genomic investigations of horse conformation.

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.

Investigating the population structure and genetic diversity of Arabian horses in Oman using SNP markers

Arabian horses were selected for metabolic efficiency, beauty, efficiency and endurance. Therefore, Bedouins have for centuries traced their prized horses' ancestries. With the establishment of the World Arabian Horse Organization (WAHO), registration of Arabian horses became centralized and countries worldwide registered them in its database. Most existing Arabian horses in Oman today were imported after the 1970s and are predominantly flat-racing Arabians. This work aimed at revealing the genetic background and diversity of Omani Arabian horses by comparing them with Arabian horses from a diverse genetic background. To that end, we genotyped 63 randomly sampled Arabian horses from Oman using the Illumina Equine SNP70. For comparison, SNP genotypes of 12 Saudi Arabian horses, 27 French, 77 Egyptian, 11 Polish and 36 US Arabians were included in the study. We additionally included 17 Thoroughbred horses and 21 horses representing large and small breeds as an outgroup. Our MDS analysis and phylogenetic analysis showed that the Arabian horses in Oman cluster primarily with French Arabian horses, with a few horses clustering within the Polish/US Arabians. The French Arabian horse cluster was the closest to the Thoroughbred horses. Amongst the Arabian horses, plink average genomic inbreeding levels were highest in the Egyptian Arabian (0.169) followed by the Saudi Arabian horses (0.137) and lowest in the Omani and French Arabian horses, -0.041 and -0.079 respectively. To our knowledge, this is the first report on the genetic background and diversity of Arabian horses in Oman. Our results demonstrated a definite subpopulation structure among Arabian horses and this information should advise future decision-making on Arabian horse breeding.

Genetics of Equine Orthopedic Disease

Orthopedic diseases are a common cause for limited exercise capacity in the horse. They often underlie genetic risk factors, which can affect bone, articular cartilage, tendons, ligaments, and adnexal structures among others. The genetic effects can directly interfere with tissue development and skeletal growth or can trigger degenerative or inflammatory processes. Many of these diseases of the locomotor system like osteochondrosis are complex and can be affected by multifactorial influences. For this reason, it is important for those performing diagnostic procedures to have a comprehensive knowledge of orthopedic diseases, their prevalence within breeds, and genetic background.

Inter- and intra-breed genome-wide copy number diversity in a large cohort of European equine breeds

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 < 0.05). We identified 80 CNVRs overlapping with known QTL for fertility, coat colour, conformation and temperament. We also report 67 novel CNVRs.

CONCLUSIONS

This work revealed that CNV patterns, in the genome of some European horse breeds, occurred in specific genomic regions. The results provide support to the hypothesis that high frequency private CNVs residing in genes may potentially be responsible for the diverse phenotypes seen between horse breeds.

Ten years of the horse reference genome: insights into equine biology, domestication and population dynamics in the post-genome era

The horse reference genome from the Thoroughbred mare Twilight has been available for a decade and, together with advances in genomics technologies, has led to unparalleled developments in equine genomics. At the core of this progress is the continuing improvement of the quality, contiguity and completeness of the reference genome, and its functional annotation. Recent achievements include the release of the next version of the reference genome (EquCab3.0) and generation of a reference sequence for the Y chromosome. Horse satellite‐free centromeres provide unique models for mammalian centromere research. Despite extremely low genetic diversity of the Y chromosome, it has been possible to trace patrilines of breeds and pedigrees and show that Y variation was lost in the past approximately 2300 years owing to selective breeding. The high‐quality reference genome has led to the development of three different SNP arrays and WGSs of almost 2000 modern individual horses. The collection of WGS of hundreds of ancient horses is unique and not available for any other domestic species. These tools and resources have led to global population studies dissecting the natural history of the species and genetic makeup and ancestry of modern breeds. Most importantly, the available tools and resources, together with the discovery of functional elements, are dissecting molecular causes of a growing number of Mendelian and complex traits. The improved understanding of molecular underpinnings of various traits continues to benefit the health and performance of the horse whereas also serving as a model for complex disease across species.

A Frameshift Mutation in KIT is Associated with White Spotting in the Arabian Camel

While the typical Arabian camel is characterized by a single colored coat, there are rare populations with white spotting patterns. White spotting coat patterns are found in virtually all domesticated species, but are rare in wild species. Theories suggest that white spotting is linked to the domestication process, and is occasionally associated with health disorders. Though mutations have been found in a diverse array of species, fewer than 30 genes have been associated with spotting patterns, thus providing a key set of candidate genes for the Arabian camel. We obtained 26 spotted camels and 24 solid controls for candidate gene analysis. One spotted and eight solid camels were whole genome sequenced as part of a separate project. The spotted camel was heterozygous for a frameshift deletion in KIT (c.1842delG, named KIT^W1 for White spotting 1), whereas all other camels were wild-type (KIT⁺/KIT⁺). No additional mutations unique to the spotted camel were detected in the EDNRB, EDN3, SOX10, KITLG, PDGFRA, MITF, and PAX3 candidate white spotting genes. Sanger sequencing of the study population identified an additional five KIT^W1/KIT⁺ spotted camels. The frameshift results in a premature stop codon five amino acids downstream, thus terminating KIT at the tyrosine kinase domain. An additional 13 spotted camels tested KIT⁺/KIT⁺, but due to phenotypic differences when compared to the KIT^W1/KIT⁺ camels, they likely represent an independent mutation. Our study suggests that there are at least two causes of white spotting in the Arabian camel, the newly described KIT^W1 allele and an uncharacterized mutation.

Host genetic influence on papillomavirus-induced tumors in the horse

The common equine skin tumors known as sarcoids have been causally associated with infection by bovine papillomavirus (BPV). Additionally, there is evidence for host genetic susceptibility to sarcoids. We investigated the genetic basis of susceptibility to sarcoid tumors on a cohort of 82 affected horses and 270 controls genotyped on a genome-wide platform and two custom panels. A Genome Wide Association Study (GWAS) identified candidate regions on six chromosomes. Bayesian probability analysis of the same dataset verified only the regions on equine chromosomes (ECA) 20 and 22. Fine mapping using custom-produced SNP arrays for ECA20 and ECA22 regions identified two marker loci with high levels of significance: SNP BIEC2-530826 (map position 32,787,619) on ECA20 in an intron of the DQA1 gene in the Major Histocompatibility Complex (MHC) class II region (p = 4.6e-06), and SNP BIEC2-589604 (map position 25,951,536) on ECA22 in a 200 kb region containing four candidate genes: PROCR, EDEM2, EIF6 and MMP24 (p = 2.14e-06). The marker loci yielded odds ratios of 5.05 and 4.02 for ECA20 and ECA22, respectively. Associations between genetic MHC class II variants and papillomavirus-induced tumors have been reported for human papillomavirus and cottontail rabbit papillomavirus infections. This suggests a common mechanism for susceptibility to tumor progression that may involve subversion of the host immune response. This study also identified a genomic region other than MHC that influenced papillomavirus-induced tumor development in the studied population.

Skeletal variation in Tennessee Walking Horses maps to the LCORL/NCAPG gene region

Conformation has long been a driving force in horse selection and breed creation as a predictor for performance. The Tennessee Walking Horse (TWH) ranges in size from 1.5 to 1.7 m and is often used as a trail, show, and pleasure horse. To investigate the contribution of genetics to body conformation in the TWH, we collected DNA samples, body measurements, and gait/training information from 282 individuals. We analyzed the 32 body measures with a principal component analysis. Principal component (PC)1 captured 28.5% of the trait variance, while PC2 comprised just 9.5% and PC3 6.4% of trait variance. All 32 measures correlated positively with PC1, indicating that PC1 describes overall body size. We genotyped 109 horses using the EquineSNP70 bead chip and marker association assessed the data using PC1 scores as a phenotype. Mixed-model linear analysis (EMMAX) revealed a well-documented candidate locus on ECA3 (raw P = 3.86 × 10(-9)) near the LCORL gene. A custom genotyping panel enabled fine-mapping of the PC1 body-size trait to the 3'-end of the LCORL gene (P = 7.09 × 10(-10)). This position differs from other reports suggesting single nucleotide polymorphisms (SNPs) upstream of the LCORL coding sequence regulate expression of the gene and, therefore, body size in horses. Fluorescent in situ hybridization analysis defined the position of a highly homologous 5 kb retrogene copy of LCORL (assigned to unplaced contigs of the EquCab 2.0 assembly) at ECA9 q12-q13. This is the first study to identify putative causative SNPs within the LCORL transcript itself, which are associated with skeletal size variation in horses.

The MC1R and ASIP Coat Color Loci May Impact Behavior in the Horse

Shared signaling pathways utilized by melanocytes and neurons result in pleiotropic traits of coat color and behavior in many mammalian species. For example, in humans polymorphisms at MC1R cause red hair, increased heat sensitivity, and lower pain tolerance. In deer mice, rats, and foxes, ASIP polymorphisms causing black coat color lead to more docile demeanors and reduced activity. Horse (Equus caballus) base coat color is primarily determined by polymorphisms at the Melanocortin-1 Receptor (MC1R) and Agouti Signaling Protein (ASIP) loci, creating a black, bay, or chestnut coat. Our goal was to investigate correlations between genetic loci for coat color and temperament traits in the horse. We genotyped a total of 215 North American Tennessee Walking Horses for the 2 most common alleles at the MC1R (E/e) and ASIP (A/a) loci using previously published PCR and RFLP methods. The horses had a mean age of 10.5 years and comprised 83 geldings, 25 stallions, and 107 mares. To assess behavior, we adapted a previously published survey for handlers to score horses from 1 to 9 on 20 questions related to specific aspects of temperament. We utilized principle component analysis to combine the individual survey scores into 4 factors of variation in temperament phenotype. A factor component detailing self-reliance correlated with genotypes at the ASIP locus; black mares (aa) were more independent than bay mares (A_) (P = 0.0063). These findings illuminate a promising and novel animal model for future study of neuroendocrine mechanisms in complex behavioral phenotypes.

A high density SNP array for the domestic horse and extant Perissodactyla: utility for association mapping, genetic diversity, and phylogeny studies

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.