Population study of the Pura Raza Español Horse regarding its coat colour

Methodological approaches to assessing the genetic predisposition of osteochondrosis of the tarsocrural joint in the Pura Raza Española horses

Osteochondrosis (OC) is a developmental orthopaedic disease of significant concern in numerous sport horse breeds, with significant international relevance. Using digital radiographs, we assessed the occurrence of hock (tarsocrural joint) OC in 3 048 Pura Raza Española (PRE) horses which took part in a morpho-functional test, in three specific locations in the tarsus limbs: the Distal Intermediate Ridge of the Tibia (DIRT), the lateral trochlear ridges of the talus (LTT), and the medial trochlear ridges of the talus (MTT). An incidence rate of 13.3% was found for hock OC in the analysed sample, with the highest incidence rate observed in DIRT (10.0%) and the lowest in MTT (0.2%). Estimates of genetic predisposition to hock OC were carried out using three genetic approaches: 1a) a binomial threshold model based on the presence or absence of OC, 1b) a multinomial threshold model, on a scale from 0 (absence) to 3 (maximum), and 2) a linear model. The effects considered in the models included sex, genetic origin and stud class. All the analyses were based on the Bayesian inference methodology, using the THRGIBBS3F90 software. The binomial threshold model yielded the most suitable results, with an estimated heritability for Overall hock OC of 0.71 ± 0.055 on the underlying scale (0.53 on the observed scale), ranging in different locations from 0.48 ± 0.087 (LTT) to 0.66 ± 0.063 (DIRT) on the underlying scale (0.10 and 0.38 on the observed scale, respectively). The highest significative genetic correlation was observed between Overall and DIRT (0.97) for approach 1a, and the lowest significant genetic correlation was between Overall and LTT (0.49), for approach 2. This study contributes valuable insights into the genetic predisposition towards, as well as for the potential for selective breeding against, hock OC in PRE horses, and provides a basis for future research and breeding programmes aimed at minimising the occurrence of hock OC and promoting the overall health of this breed.

Genomic tools for early selection among Thoroughbreds and Polo Argentino horses for practicing polo

The Polo Argentino (PA) horse is a recognized breed, developed originally by mixing crossbred and Thoroughbred (TB) horses to play polo. Early PA selection is difficult due to unreliable performance estimations. This study investigated the usefulness of genomic markers previously linked to morphological and functional traits as a tool for the early selection of PA. To this, we genotyped 520 PA and 30 TB horses using the Equine GGPArray (Illumina, n = 71,778 SNPs). Analyses included a genetic characterization of six genetic markers associated with behavioral (DRD4), muscular development (MSTN), and body size (LCORL, HMGA6, ZFAT, and LASP1) genes. Genetic differences in the DRD4, MSTN, and LCORL SNP were found between the two breeds, in the last two FST index between breeds was 0.13 and 0.6, respectively (p < 0.01). In DRD4, G allele was the more prevalent in PA (0.56 vs 0.45 in TB, p < 0.05), but no differences were observed between the genotypes associated with phenotypes. In MSTN, heterozygous genotypes were the most common in PA (48 %), with a significant decrease in AA (Hardy-Weinberg p < 0.05), suggesting a negative selection against it in polo horses. In body size, HMGA2 was monomorphic in all horses, while ZFAT and LASP1 SNP showed higher variability. Interestingly, 99 % of PA showed a TT genotype in LCORL (only 66 % in TB), demonstrating selection for smaller horses. Our results suggest that empirical selection in PA has generated an incipient genomic differentiation in discrete traits which could be used as a marker-assisted selection tool for early selection of polo horses.

Quantitative analysis of parent-of-origin effect in reproductive and morphological selection criteria in the Pura Raza Española horse

It is generally assumed that parents make a genetically equal contribution to their offspring, but this assumption might not always hold. This is because the expression of a gene can be blocked by methylation during gametogenesis, and the degree of methylation can depend on the origin of the parental gene (imprinting) or by preferential management associated with genetic merit. The first consequences of this for quantitative genetics is that the mean phenotypes of reciprocal heterozygotes need no longer be the same, as would be expected according to Mendelian heritage. We analysed three mare reproductive traits (reproductive efficiency, age at first foaling and foaling number) and three morphological traits (height at withers, thoracic circumference, and scapula-ischial length) in the Pura Raza Española (PRE) horse population, which possesses a deep and reliable pedigree, making it a perfect breed for analysing the quantitative effect of parent-of-origin. The number of animals analysed ranged from 44,038 to 144,191, all of them with both parents known. The model comparison between a model without parent-of-origin effects and three different models with parent-of-origin effects revealed that both maternal and paternal gametic effects influence all the analysed traits. The maternal gametic effect had a higher influence on most traits, accounting for between 3% and 11% of the total phenotypic variance, while the paternal gametic effect accounted for a higher proportion of variance in one trait, age at first foaling (4%). As expected, the Pearson's correlations between additive breeding values of models that consider parent-of-origin and that do not consider parent-of-origin were very high; however, the percentage of coincident animals slightly decreases when comparing animals with the highest estimated breeding values. Ultimately, this work demonstrates that parent-of-origin effects exist in horse gene transmission from a quantitative point of view. Additionally, including an estimate of the parent-of-origin effect within the PRE horse breeding program could be a great tool for a better parent's selection and that could be of interest for breeders, as this value will determine whether the animals acquire genetic categories and are much more highly valued.

Phenotypic and Genetic Study of the Presence of Hair Whorls in Pura Raza Español Horses

Hair whorls are a hereditary feature in horses that may be associated with temperament and coat color. Hair whorls are described as changes in the hair pattern and may take various forms, such as circular and linear whorls. We first carried out a frequency analysis of hair whorls (circular and linear). Next, a Generalized Non-Linear Model was computed to assess the significance of some potential influencing factors, and a genetic parameter estimation was performed. ENDOG software v4.8 was used to estimate the inbreeding coefficient of all the animals analyzed. It was more common to find horses with circular hair whorls than with linear whorls. The heritability ranges obtained were, in general, medium-high for both circular whorls (0.20 to 0.90) and linear whorls (0.44 to 0.84). High positive correlations were found on the between left and right positions, indicating a tendency to symmetry in certain locations. The laterality of hair whorls was also evidenced, with the biggest concentration on the left-hand side, particularly in gray horses, showing circular whorls below the central line of eyes, which has been associated in a previous paper with a calmer and more docile temperament.

Fine-tuning genomic and pedigree inbreeding rates in equine population with a deep and reliable stud book: the case of the Pura Raza Española horse

BACKGROUND

Estimating inbreeding, which is omnipresent and inevitable in livestock populations, is a primary goal for management and animal breeding especially for those interested in mitigating the negative consequences of inbreeding. Inbreeding coefficients have been historically estimated by using pedigree information; however, over the last decade, genome-base inbreeding coefficients have come to the forefront in this field. The Pura Raza Española (PRE) horse is an autochthonous Spanish horse breed which has been recognised since 1912. The total PRE population (344,718 horses) was used to estimate Classical (F), Ballou's ancestral, Kalinowski's ancestral, Kalinowski's new and the ancestral history coefficient values. In addition, genotypic data from a selected population of 805 PRE individuals was used to determine the individual inbreeding coefficient using SNP-by-SNP-based techniques (methods of moments -F_HOM-, the diagonal elements of the genomic -F_G-, and hybrid matrixes -F_H-) and ROH measures (F_RZ). The analyse of both pedigree and genomic based inbreeding coefficients in a large and robust population such as the PRE horse, with proven parenteral information for the last 40 years and a high degree of completeness (over 90% for the last 70 years) will allow us to understand PRE genetic variability better and the correlations between the estimations will give the data greater reliability.

RESULTS

The mean values of the pedigree-based inbreeding coefficients ranged from 0.01 (F for the last 3 generations -F3-) to 0.44 (ancestral history coefficient) and the mean values of genomic-based inbreeding coefficients varied from 0.05 (F_RZ for three generations, F_H and F_HOM) to 0.11 (F_RZ for nine generations). Significant correlations were also found between pedigree and genomic inbreeding values, which ranged between 0.58 (F3 with F_HOM) and 0.79 (F with F_RZ). In addition, the correlations between F_RZ estimated for the last 20 generations and the pedigree-based inbreeding highlight the fact that fewer generations of genomic data are required when comparing total inbreeding values, and the opposite when ancient values are calculated.

CONCLUSIONS

Ultimately, our results show that it is still useful to work with a deep and reliable pedigree in pedigree-based genetic studies with very large effective population sizes. Obtaining a satisfactory parameter will always be desirable, but the approximation obtained with a robust pedigree will allow us to work more efficiently and economically than with massive genotyping.

Genetic inbreeding depression load for fertility traits in Pura Raza Española mares

Fertility is a key factor in the economic success of horse farms. However, it has received little attention due to the difficulty of measuring fertility objectively. Since its studbook creation (1912), the Pura Raza Española (PRE) breed has been a closed population and become high in-bred resulting in inbreeding depression (poor phenotypic values). Nevertheless, heterogeneous effects of inbreeding depression have been detected among founders and nonfounders. The aims of this study were (1) to analyze the genetic parameters for reproductive traits in mares of the PRE horse breed and (2) to estimate, for the first time, the inbreeding depression load associated with common ancestors of the breed. A total of 22,799 mares were analyzed. Heritability estimates ranged from 0.05 (interval between first and second foaling) to 0.16 (age at first foaling), whereas inbreeding depression load ratios ranged from 0.06 (parturition efficiency at 6th foaling) to 0.17 (age at first foaling), for a partial inbreeding coefficient of 10%. Although heritability is related to the variability expressed in the population, inbreeding depression load ratios measure the potential variability, whether expressed in the population or not. Most correlations between additive and inbreeding depression load genetic values were significant (P < 0.001) and of low to moderate magnitude. Our results confirm that individual inbreeding depression loads allow us to select horses that have a genetic value resistant to the deleterious effects of inbreeding.

500 years of breeding in the Carthusian Strain of Pura Raza Español horse: An evolutional analysis using genealogical and genomic data

The Carthusian horse is a Pura Raza Español (PRE) strain (CS), bred as a closed population since its creation more than 500 years ago. The aim of this study was to analyse for the first time its population structure and situation of variability combining both genealogical (GEL) and genomic (GEN) data. The GEL data comprised 348,429 pedigree records (56,105 CS horses), while the GEN analysis included the high-density genotypes (670,804 SNPs) of 287 horses. Pedigree completeness demonstrated its accuracy, showing a good correlation of GEL (F) and GEN (F_ROH ) inbreeding coefficient in the case of PRE subpopulations partially related and non-related to Carthusian strain (0.68) but a lower value in the 100% Carthusian horses (0.42), due to the high weight of founders not detected by GEL analysis. GEN (PCA, AMOVA, and Admixture) and GEL analysis showed a good differentiation of subpopulations, but also a high level of introgression of the CS in the breed during past decades. A recent change in this trend was noteworthy, with a considerable reduction in CS variability and a genetic bottleneck (effective population sizes of 31.57 and 30.20 in GEL and GEN analysis, respectively, in last generation). The PRE has maintained its variability, and a considerable difference in estimated N_e by GEL (60.77) and GEN (188.0) data was observed. Using two sources of complementary information, it was found the existence of an ancient PRE strain with a unique genetic landmark, practically free from the influence of other equine populations.

Genetic inbreeding depression load for morphological traits and defects in the Pura Raza Española horse

BACKGROUND

Inbreeding is caused by mating between related individuals and is associated with reduced fitness and performance (inbreeding depression). Several studies have detected heterogeneity in inbreeding depression among founder individuals. Recently, a procedure was developed to predict hidden inbreeding depression load that is associated with founders using the Mendelian sampling of non-founders. The objectives of this study were to: (1) analyse the population structure and general inbreeding, and (2) test this recent approach for predicting hidden inbreeding depression load for four morphological traits and two morphology defects in the Pura Raza Española (PRE) horse breed.

RESULTS

The regression coefficients that were calculated between trait performances and inbreeding coefficients demonstrated the existence of inbreeding depression. In total, 58,772,533 partial inbreeding coefficients (Fij) were estimated for the whole PRE population (328,706 horses). We selected the descendants of horses with a Fij ≥ 6.25% that contributed to at least four offspring and for which morphological traits were measured for the subsequent analysis of inbreeding depression load (639 horses). A pedigree was generated with the last five generations (5026 animals) used as the reference population (average inbreeding coefficient of 8.39% and average relatedness coefficient of 10.76%). Heritability estimates ranged from 0.08 (cresty neck) to 0.80 (height at withers), whereas inbreeding depression load ratios ranged from 0.01 (knock knee) to 0.40 (length of shoulder), for an inbreeding coefficient of 10%. Most of the correlations between additive and inbreeding depression load genetic values and correlations between inbreeding depression load genetic values for the different traits were positive or near 0.

CONCLUSIONS

Although the average inbreeding depression loads presented negative values, a certain percentage of the animals showed neutral or even positive values. Thus, high levels of inbreeding do not always lead to a decrease in mean phenotypic value or an increase in morphological defects. Hence, individual inbreeding depression loads could be used as a tool to select the most appropriate breeding animals. The possibility of selecting horses that have a high genetic value and are more resistant to the deleterious effects of inbreeding should help improve selection outcomes.

Survey of Risk Factors and Genetic Characterization of Ewe Neck in a World Population of Pura Raza Español Horses

Simple Summary

Ewe Neck is a common morphological defect of the Pura Raza Español (PRE) population, which seriously affects the horse’s development. In this PRE population (35,267 PRE), a total of 9693 animals (27.12% of total) was Ewe Neck-affected. It has been demonstrated that genetic and risk factors (sex, age, geographical area, coat color, and stud size) are involved, being more prevalent in the males, 4–7 years old, chestnut coat, from small studs (less than 5 mares), and raised in North America. The morphological traits height at chest, length of back, head-neck junction, and bottom neck-body junction and the body indices, head index, and thoracic index were those most closely related with the appearance of this morphological defect. The additional genetic base of Ewe Neck in PRE, which presents low-moderate heritability (h2: 0.23–0.34), shows that the prevalence of this defect could be effectively reduced by genetic selection.

Abstract

Ewe Neck is a relatively common morphological defect in Pura Raza Español (PRE) horses and other Baroque type horse breeds, which adversely affects the breeding industry; (1) objectives: to establish the within-breed prevalence, possible associated factors, and heritability of Ewe Neck in PRE horses; (2) methods: the database included evaluations of 35,267 PRE horses. The Ewe Neck defect, 16 morphological traits, and 4 body indices were recorded. A Bayesian genetic animal model included the following systematic effects: sex, age, coat color, geographical area of the stud, and birth stud size were used; (3) results: in this PRE population, a total of 27.12% was affected. All the risk factors studied were significantly associated with the Ewe Neck score. The heritability coefficient for Ewe Neck score ranged from 0.23 to 0.34. Morphological traits (height at chest, length of back, head-neck junction, and bottom neck-body junction) and the indices (head and thoracic index) were those most closely related with the appearance of Ewe Neck; (4) conclusions: Ewe Neck is a relatively frequent defect in PRE horses, associated with risk factors and other morphological traits, with a moderate level of heritability. Breeding to select against this condition may therefore be beneficial in this breed.

Genetic Structure Analysis of the Pura Raza Español Horse Population through Partial Inbreeding Coefficient Estimation

The aim of this work was to analyze genetic parameters such as the inbreeding coefficient (F), relatedness coefficient (AR) and partial inbreeding coefficient (Fij) of the whole PRE population, and the ancestors which account for 50% of the total genetic variability of the current population, from genealogical information. The average F of the whole PRE population (328,706 animals) has decreased from 8.45% to 7.51% in the least 20 years. The Fij was estimated for the whole PRE population, resulting in a database of 58,772,533 records containing one record for each Fij that each animal receives from a certain common ancestor (CA). A total of 10,244 CAs contributed to the Fij with an average of 5370 descendants, with each descendant having an average of 170 CAs. Over the generations, the number of CAs has increased, while the proportion of Fij by each one has decreased. In addition, the contributions of the more influential ancestors have changed. The increased census, the limited use of artificial insemination and our increased knowledge about inbreeding depression and the animals' breeding values allow breeders to select horses more for their functionality and conformation than for their pedigree reputation, which is the cause of all these changes.