Swedish Warmblood horses (SWB) are bred for show jumping and/or dressage with young horse test scores as indicator traits. This study aimed to investigate possible candidate genes and regions of importance for evaluated and linearly scored young horse test traits. A single-step genome-wide association study (ssGWAS) was done using the BLUPF90 suite of programs for factors scores from factor analysis of traits assessed at young horse tests together with height at withers. The ssGWAS included 20,814 SWB with factors scores for four factors for evaluated traits. A total of 6436 of these horses also had factor scores for 13 factors for linearly scored traits. Genotypes from a 670K SNP array were available for 380 of the horses in this study. All genotyped horses had factor scores for evaluated traits, and 379 also had factors scores for linearly scored traits. Significant SNPs associated with three factors related to size were located on ECA3 within or nearby a well-known region, including the genes ligand dependent nuclear receptor corepressor like (LCORL), non-SMC condensin I complex subunit G (NCAPG), DDB1 and CUL4 Associated Factor 16 (DCAF16), and the Family with Sequence Similarity 184 Member B (FAM184B). Significant SNPs were also detected for two factors for evaluated traits representing conformation and jumping, and four factors for linearly scored traits related to body length, neck conformation, walk and trot (hindleg position and activity), respectively. Among nearby genes, calcium/calmodulin-dependent protein kinase type 1D (CAMK1D) for the factor for linearly scored traits related to neck conformation and GLI Family Zinc Finger 2 (GLI2) for the factor for evaluated jumping traits, were most promising. For these, top associated SNPs were detected within the genes, and the known gene functions seems to be related to the phenotypes. In conclusion, ssGWAS is beneficial to detect plausible candidate genes/regions for desired traits in warmblood horses.

Assessment protocols to describe the various aspects of conformation, gait and jumping traits on a linear scale were introduced at young horse tests for Swedish Warmblood horses in 2013. The traits scored on a linear scale are assumed to be less subjective and more easily compared across populations than the traditional evaluated traits that are scored relative to the breeding goal. However, the resulting number of traits is considerable, and several of the traits are correlated. The aim of this study was to investigate the interrelationship between the different evaluated and linearly scored traits in Swedish Warmbloods using factor analysis. In total, 20,935 horses born 1996-2017 had information on evaluated traits, and 5450 of these also had linearly scored trait records assessed since 2014 when the protocol was updated. A factor analysis with varimax rotation was performed separately for evaluated and linearly scored traits using the Psych package in R. Height at withers was included in both analyses. A total of four factors for evaluated traits and 14 factors for linearly scored traits were kept for further analysis. Missing values for individual traits in horses with linearly scored trait records were imputed based on correlated traits before factor scores were calculated using factor loadings. Genetic parameters for, and correlations between, the resulting underlying factors were estimated using multiple-trait animal models in the BLUPF90 package. Heritability estimates were on a similar level as for the traits currently used in the genetic evaluation, ranging from 0.05 for the factor for linearly scored traits named L.behaviour (dominated by traits related to behaviour) to 0.59 for the factor for evaluated traits named E.size (dominated by height at withers and conformation). For both types of traits, separate factors were formed for jumping and gait traits, as well as for body size. High genetic correlations were estimated between such corresponding factors for evaluated traits and factors for linearly scored traits. In conclusion, factor analysis could be used to reduce the number of traits to be included in multiple-trait genetic evaluation or in genomic analysis for warmblood horses. It can also contribute to a better understanding of the interrelationships among the assessed traits and be useful to decide on subgroups of traits to be used in several multiple-trait evaluations on groups of original traits.

Background Warmblood fragile foal syndrome (WFFS) is a monogenetic defect caused by a recessive lethal missense point mutation in the procollagen-lysine, 2-oxoglutarate 5-dioxygenase 1 gene (PLOD1, c.2032G>A). The majority of homozygous WFFS horses are aborted during gestation. Clinical signs of affected horses include fragile skin, skin and mucosa lacerations, hyperextension of the articulations, and hematomas. In spite of its harmful effect, a relatively high frequency of WFFS carriers has been found in Warmblood horses, suggesting a heterozygote advantage. Thus, in this study our aims were to: (1) estimate the frequency of WFFS carriers in the Swedish Warmblood breed (SWB), (2) estimate the effect of WFFS carrier genotype on performance traits in two SWB subpopulations bred for different disciplines, and (3) simulate the potential effects of balancing selection and different selection strategies on the frequency of carriers. Methods In total, 2288 SWB sport horses born between 1971 and 2020 were tested for the WFFS mutation and had estimated breeding values (EBV) for ten traditional evaluating and 50 linear descriptive traits. Results The frequency of WFFS carriers calculated from a pool of 511 randomly selected SWB horses born in 2017 was equal to 7.4% and ranged from 0.0 to 12.0% among the whole set of tested SWB horses, starting from 1971 till 2020. The effect of the WFFS carrier genotype was significant for several EBV mainly related to movements and dressage traits and especially for horses not bred for the show jumping discipline. Using simulation, we showed that balancing selection can maintain a recessive lethal allele in populations such as the SWB breed over generations and that the frequency is expected to slowly decrease in absence of balancing selection. Finally, we showed that selection against carrier sires can result in a more rapid decrease of the frequency of the mutant allele over time. Conclusion Further research is needed to confirm the apparent association between equine performance and the WFFS carrier genotype. Identification of such associations or new causative mutations for horse performance traits can serve as new tools in horse breeding to select for healthy, sustainable, and better performing horses.

In many European warmblood studbooks, clear specialization toward either jumping or dressage horses is evident. The Swedish Warmblood (SWB) is also undergoing such specialization, creating a possible need for separate breeding programs and a discipline-specific Young Horse Test (YHT). This study investigated how far specialization of the SWB breed has proceeded and the potential consequences. Individuals in a population of 122,054 SWB horses born between 1980 and 2020 were categorized according to pedigree as jumping (J), dressage (D), allround (AR), or thoroughbred (Th). Data on 8,713 J horses and 6,477 D horses assessed for eight traits in YHT 1999-2020 were used to estimate genetic parameters within and between J and D horses and between different periods. Future scenarios in which young horses are assessed for either jumping or dressage traits at YHT were also analyzed. More than 80% of horses born in 1980-1985 were found to be AR horses, while 92% of horses born in 2016-2020 belonged to a specialized category. The average relationship within J or D category was found to increase during the past decade, whereas the relationship between these categories decreased. Heritability estimates for gait traits were 0.42-0.56 for D horses and 0.25-0.38 for J horses. For jumping traits, heritability estimates were 0.17-0.26 for J horses and 0.10-0.18 for D horses. Genetic correlations between corresponding traits assessed in J and D horses were within the range 0.48-0.81, with a tendency to be lower in the late study period. In the future scenarios, heritability and genetic variance both decreased for traits that were not assessed in all horses, indicating that estimation of breeding value and genetic progress for these traits could be affected by a specialized YHT. However, ranking of sires based on estimated breeding values (EBVs) and accuracy of EBVs was only slightly altered for discipline-specific traits. With continued specialization in SWB, specialization of the YHT should thus be considered.

Additional file 4: Table S3. Descriptive statistics of the traditional evaluating traits in the SJ horses. The following values are reported for each traditional evaluating trait in the SJ horses: number of observations, mean, standard deviation, minimum, 25th percentile, 75th percentile and maximum.

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.

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.

AbstractRecessive lethal variants often segregate at low frequencies in animal populations, such that two randomly selected individuals are unlikely to carry the same mutation. However, the likelihood of an individual inheriting two copies of a recessive lethal mutation is dramatically increased by inbreeding events. Such occurrences are particularly common in domestic animal populations, which are often characterised by high rates of inbreeding and low effective population sizes. To date there have been no published investigations into the presence of specific variants at high frequencies in domestic horse populations. 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 scanned genotype data from Thoroughbred horses(n= 526) for adjacent single nucleotide polymorphisms (SNPs) at high heterozygote frequencies, but with a complete absence of homozygotes. Two SNPs that matched these criteria were mapped to an intronic region in theLY49Bgene, indicating that a closely linked mutation may cause lethality in homozygous state. Despite a complete absence of homozygotes, almost 35% of Thoroughbreds included in these analyses were heterozygous for both SNPs. A similar loss or absence of homozygotes was observed in genotype data from other domestic horse breeds (n= 2030). Variant analysis of whole-genome sequence data (n= 90) identified two SNPs in the 3’UTR region of theLY49Bgene that may result in loss of function. Analysis of transcriptomic data from equine embryonic tissue revealed thatLY49Bis expressed in the trophoblast during placentation stage of development.In this study, a region in theLY49Bgene was identified as a strong candidate for harbouring a variant causing lethality in homozygous state. These findings suggest thatLY49Bmay 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.Author SummaryRecessive lethal mutations may reach high frequencies in livestock populations due to selective breeding practices, resulting in reduced fertility rates. In this study, we characterise recessive lethal mutations at high frequencies in the Thoroughbred horse population, a breed with high rates of inbreeding and low genetic diversity. We identified a haplotype in theLY49Bgene that shows strong evidence of being homozygous lethal, despite having high frequencies of heterozygotes in Thoroughbreds and other domestic horse breeds. Two 3’UTR variants were identified as most likely to cause loss of function in theLY49Bgene, resulting in lethality. This finding provides novel insights into the potential importance ofLY49Bin equine development. Additionally, this study may assist with breeding strategies to improve fertility rates in the Thoroughbred and other domestic horse breeds.

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 (F ST 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.

Breeding against canine hip dysplasia (HD) may benefit from the importation of foreign sires. When foreign sires are evaluated on a different HD scale, this may diminish the efficacy. Using stochastic simulations, we evaluated genetic change and inbreeding levels for different scenarios of importing sires with high genetic merit for HD. Population size and genetic parameters (e.g. heritability, accuracy of selection, genetic correlation) were based on actual data for HD in Golden retrievers and Labrador retrievers in the UK and Sweden. For countries with different HD scales and an estimated breeding value (EBV) evaluation in place, the importation was useful if imported sires had EBV rankings in the top 50% and if genetic correlations between EBV systems were above 0.85. When importing sires with EBV rankings in the top 10%, moderate accuracies of EBVs (>0.40) and moderately strong genetic correlations (>0.70) were needed. Selection against HD without the importation of sires may increase inbreeding levels, while the importation of sires can decrease inbreeding levels. For national genetic evaluation and selection programmes, importing sires with high genetic merit can be an effective breeding strategy, but care is needed to estimate reliable EBVs.

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.

The equestrian sport horse Swedish Warmblood (SWB) originates from versatile cavalry horses. Most modern SWB breeders have specialized their breeding either towards show jumping or dressage disciplines. The aim of this study was to explore the genomic structure of SWB horses to evaluate the presence of genomic subpopulations, and to search for signatures of selection in subgroups of SWB with high or low breeding values (EBVs) for show jumping. We analyzed high density genotype information from 380 SWB horses born in the period 2010–2011, and used Principal Coordinates Analysis and Discriminant Analysis of Principal Components to detect population stratification. Fixation index and Cross Population Extended Haplotype Homozygosity scores were used to scan the genome for potential signatures of selection. In accordance with current breeding practice, this study highlights the development of two separate breed subpopulations with putative signatures of selection in eleven chromosomes. These regions involve genes with known function in, e.g., mentality, endogenous reward system, development of connective tissues and muscles, motor control, body growth and development. This study shows genetic divergence, due to specialization towards different disciplines in SWB horses. This latter evidence can be of interest for SWB and other horse studbooks encountering specialized breeding.

Additional file 2. List of genes overlapping with the identified CNVRs.

A linear profiling protocol was introduced in 2013 at tests for 3-year-old Swedish Warmblood horses. In this protocol, traits are subjectively described on a nine-point linear scale from one biological extreme to the other. This complements the traditional scoring where horses are evaluated in relation to the breeding objective. This study aimed to investigate the suitability of the linear information for genetic evaluation. Data on 22 conformation traits, 17 movement traits, 14 jumping traits and one temperament trait from 3,410 horses tested between 2013 and 2016 were analysed using an animal model. For conformation traits, the heritabilities ranged from 0.10 for description of hock joint from behind to 0.52 for shape of the neck. For movement traits, the highest heritability (0.54) was estimated for elasticity in trot and the lowest (0.08) for energy in walk. The heritabilities for jumping traits ranged from 0.05 for the ability to focus on the assignment to 0.57 for scope. Genetic correlations between linear traits and corresponding traditionally scored traits were strong (-0.37 to in many cases <-0.9). The results show that the linear information is suitable for genetic evaluation and can be a useful tool for breeders.

In order to reduce the prevalence of inherited diseases in pedigree dogs, the feasibility of implementation of an international breeding program was investigated. One prerequisite is a strong genetic correlation between countries and our objective was to estimate this correlation for canine hip dysplasia (HD) across three countries to evaluate the feasibility of an international genetic evaluation. Data were provided by the Société Centrale Canine (SCC, France), Svenska Kennelklubben (SKK, Sweden) and The Kennel Club (KC, UK) on Golden retriever and Labrador retriever dogs. Trivariate analysis on the three different modes of scoring HD in France, Sweden and the UK was performed using a mixed linear animal model. Heritability, genetic correlation, number of common sires, genetic similarity, selection differentials and accuracy of selection were calculated. The estimated heritabilities of Golden retrievers (Labrador retrievers) for HD scores were 0.28 (0.15), 0.28 (0.29) and 0.41 (0.34) in France, Sweden and the UK, respectively. The feasibility of performing a genetic evaluation of HD across countries was indicated by the favourable genetic correlations estimated between score modes (ranged from 0.48 to 0.99). The accuracy of selection for the most recent birth year cohorts of male dogs was not improved by international evaluation compared to national evaluation. Improvement in genetic progress can however be achieved by selection across populations in different countries, particularly for small populations, which were indicated by the large difference between selection differentials based on the national and international evaluations.

There is substantial international influence on the Swedish Warmblood (SWB) sport horse population. The SWB Association suspects that imported horses are pre-selected based on their anticipated performance level, which could lead to biased estimated breeding values (EBVs) for stallions. This study examined different strategies to handle records for imported horses in the genetic evaluation. The stallions were evaluated for 10 traits using 3 different EBV estimation methods that were compared based on ranking comparisons of stallions, accuracy, correlations and absolute differences between EBVs. The results showed that the stallions' EBVs were affected by imported offspring which, as a group, had higher average scores than horses born in Sweden. Sire EBVs differed by up to 0.6 genetic standard deviation of the trait between methods for stallions with >50% imported offspring. Excluding imported offspring lowered the accuracies and caused larger re-ranking compared with including a fixed effect of origin in the model.

Breeding programmes for warmblood sport horses are similar in the Nordic countries Sweden, Denmark, Finland and Norway, and stallions of same origin are used. The aim was to investigate whether a joint Nordic genetic evaluation based on lifetime competition performance is feasible and beneficial for breeding competitive sport horses in the Nordic countries. Results for almost 45,000 horses in show jumping and 30,000 horses in dressage were available. The larger populations in Sweden and Denmark contributed with 85% of the results. Heritabilities and genetic correlations between performances in the different countries were estimated, and comparisons of accuracies of estimated breeding values (EBVs) and number of stallions with EBVs based on national or joint data were studied. The heritabilities ranged between 0.25 and 0.42 for show jumping and between 0.14 and 0.55 for dressage. The genetic correlations between competition performances in the Nordic countries were estimated to 0.63-1.00. EBVs based on joint data increased accuracies for EBVs for stallions by 38-81% and increased the number of available stallions with EBVs by 40-288%, compared to EBVs based on national data only. A joint Nordic genetic evaluation for sport horses is recommended.

Young horse performance test data from two warmblood riding horse populations, Norwegian warmblood (NWB) and Swedish warmblood (SWB), were analysed to examine whether including information from a related studbook would increase the accuracy of the genetic evaluations within a population. Ten conformation and performance traits from 31,588 horses, 774 NWB and 30,814 SWB were analysed separately and jointly using single trait animal models. Heritabilities were moderate to high, and varied from 0.15 (conformation, joint data) to 0.74 (jumping technique, NWB data). The genetic similarity (GS) between populations was 31%, with the SWB, as expected given the size of the populations, contributing most to the GS (98%). Genetic correlations between the same traits in the two populations were 0.43-0.90 but with large standard errors (0.2-0.3). Including information from the other population increased the average accuracy of estimated breeding values for common stallions, on average 4% for SWB and 110% for NWB.