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

Reproductive traits in Pura Raza Española mares manifest inbreeding depression from low levels of homozygosity

Inbreeding depression is a genetic phenomenon associated with the loss of fitness and mean phenotypic performance due to mating between relatives. Historically, inbreeding coefficients have been estimated from pedigree information. However, the onset of genomic selection programs provides large datasets of individuals genotyped using SNP arrays, enabling more precise assessment of an individual's genomic-level inbreeding using genomic data. One of the traits most sensitive to issues stemming from increased inbreeding is reproduction. This is particularly important in equine, in which fertility is only moderate compared to other livestock species. To explore this further, we evaluated the effect of inbreeding on five reproductive traits (age at first foaling (AFF), average interval between foalings (AIF), total number of foalings (NF), productive life (PL) and reproductive efficiency (RE)) in Pura Raza Español mares using genomic data. Residual predicted phenotypes were obtained by purging these traits through the REML (wgResidual) and ssGREML (gResidual) approaches in reproductive data of 29,847 PRE mares using the BLUPF90+ program. Next, we used pedigree-based (Fped) and ROH-based genomic (FROH) inbreeding coefficients derived from 1018 animals genotyped with 61,271 SNPs to estimate the inbreeding depression (linear regression). Our results indicated significant levels of inbreeding depression for all reproductive traits, with the exception of the AIF trait when Fped was used. However, all traits were negatively affected by the increase in genomic inbreeding, and FROH was found to capture more inbreeding depression than Fped. Likewise, REML models (ssGREML) using genomic data for estimated predicted residual phenotypes resulted in higher variance explained by the model compared with the models not using genomics (REML). Finally, a segmented regression analysis was conducted to evaluate the effect of inbreeding depression, revealing that the levels of genealogical and genomic homozygosity do not manifest uniformly in reproductive traits. In contrast, the levels of inbreeding depression ranged from low to high as homozygosity increased. This analysis also showed that reproductive traits are very sensitive to inbreeding depression, even with relatively low levels of homozygosity.

Designing an early selection morphological traits index for reproductive efficiency in Pura Raza Española mares

The low rate of reproductive efficiency in horses may be linked to the equine industry's practice of maintaining breeding stock that stands out for its athletic or morphological performance but exhibits poor reproductive efficiency. In this study, the age at first foaling, age at last foaling, interval between first and second foaling, average interval between foaling, foaling number and reproductive efficiency, and their relationship with the morphological records in 19,758 Pura Raza Española (PRE) mares were analyzed. After a partial least squares analysis height at withers, lateral hock angle, dorsal-sternal diameter, perimeter of anterior cannon bone, angle of shoulder, thoracic perimeter, hip-stifle distance, and angle of croup were the most important traits related with reproductive traits in PRE mares. A multivariate animal model with both morphological and reproductive traits was applied, including age, geographical area, coat color, and average stud size in the decade of the mare's first foaling. The results indicate that reproductive selection is feasible, and is enhanced by the use of morphological traits, given the moderate to high magnitudes of heritability values in both types of traits, which oscillated between 0.1 (interval between first and second foaling) and 0.95 (height at withers). The resulting genetic parameters were used to develop a series of selection indexes based on morphological or morphological and reproductive combination traits to improve reproductive efficiency traits in PRE mares and thus compute the expected genetic response (EGR) for the different strategies. The increase in EGR, when comparing EGR using only reproductive traits as selection criteria vs. using a combined index of both morphological and reproductive traits, oscillated between 4.0% for the age at first foaling to 46.8% for the interval between first and second foaling. In addition, two precocity indexes showed positive EGR when the age at last foaling and the interval between first and second foaling, both with and without morphological traits, were used as selection criteria. Therefore, this analysis reveals that the preselection of reproductive traits based on morphological traits is possible in PRE mares. Ultimately, this knowledge will help breeders achieve genetic progress in reproductive traits, leading to healthier and more successful breeding outcomes in horses.

Orthopedic Diseases in the Pura Raza Española Horse: The Prevalence and Genetic Parameters of Angular Hoof Deviations

Abnormalities in hoof shape are usually connected with limb conformation defects. The role of angular hoof deviations is important for longevity in sports competitions and is increasingly recognized as a factor associated with lameness in performance horses. In this paper, we measured the prevalence of four defects related to the angulation of the hoof in the Pura Raza Española horse (PRE): splay-footed forelimb (SFF), pigeon-toed forelimb (PTF), splay-footed rear limb (SFR), and pigeon-toed rear limb (PTR). A total of 51,134 animals were studied, of which only 15.75% did not have any of the four angular hoof defects investigated, while 26.61%, 23.76%, 79.53%, and 3.86% presented SFF, PTF, SFR, and PTR, respectively. Angular defects were evaluated using two different models; model A was a linear scale composed of three categories, where 0 corresponded to the absence of defects, 1 to a minor presence of the defect and 2 to the highest degree of the defect. Model B was composed of two categories, where 0 corresponded to the absence of defects and 1 to the presence of defects, joining classes 1 and 2. We measured the factors influencing the appearance of these defects: age, inbreeding coefficient, sex, and birth stud size. The heritability of each defect was also estimated using a multivariate animal model, using the Gibbsf90+ software from the BLUPF90 family, resulting in heritability estimates of 0.18 (s.d. = 0.009), 0.20 (s.d. = 0.010), 0.11 (s.d. = 0.009), and 0.31 (s.d. = 0.010) for SFF, PTF, SFR, and PTR defects, respectively, for model A, and 0.17 (s.d. = 0.008), 0.19 (s.d. = 0.009), 0.11 (s.d. = 0.009), and 0.29 (s.d. = 0.009) for SFF, PTF, SFR, and PTR defects, respectively, for model B. Finally, the genetic correlation between the diameter of the superficial digital flexor tendon (SDFT) and the proportionality index (PI) in relation to the higher or lower prevalence of the defects was analyzed. We concluded that diameter of SDFT development is strongly correlated with inward toe conditions (PTF, PTR; P≠0 ≥ 0.95), while PI is associated with outward toe defects (SFF, SFR; P≠0 ≥ 0.95).

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.

The Social and Reproductive Challenges Faced by Free-Roaming Horse (Equus caballus) Stallions

In captivity, intact male horses, due to their sexual drive, are usually socially isolated from other horses. This lifestyle strongly contrasts with that experienced by horses living in free-roaming, feral, or semi-feral conditions, where adult stallions have several roles in their social group, with successful reproduction being their primary drive. Reproductive skew in wild populations is high; many stallions will fail to reproduce at all, while others achieve high levels of reproductive success, siring a large number of foals. Successful stallions are those with particular characteristics and abilities that facilitate harem formation and tenure, allowing them to successfully take over a harem or establish a new one, protect mares from rival stallions, employ appropriate social behaviour to maintain group cohesion, and avoid kin-mating, for example through kin recognition mechanisms. Whilst the life of free-living stallions is far from stress-free, they retain ancestral adaptations to selection pressures (such as predation and competition) exhibited by their natural environment over thousands of years. Here, we discuss the challenges faced by free-living horse stallions, the roles they play in social groups, and their resulting social needs. By understanding these pressures and how stallions react to them, we highlighted the importance of the social environment for the stallion. It is hoped that a better understanding of wild stallions’ lives will lead to their needs being more clearly met in captivity, reducing stereotypical behaviour and improving welfare.

Analysis of Inbreeding Effects on Survival at Birth of Pannon White Rabbits Using the Inbreeding-Purging Model

Mating between related animals is an inevitable consequence of a closed population structure especially when it coincides with a small population size. As a result, inbreeding depression may be encountered especially when considering fitness traits. However, under certain circumstances, the joint effects of inbreeding and selection may at least partly purge the detrimental genes from the population. In the course of this study, our objective was to determine the status of purging and to quantify the magnitude of the eliminated genetic load for the survival at birth of Pannon White rabbit kits maintained in a closed nucleus population. The evolution of the survival at birth was evaluated by applying the PurgeR R package based on the inbreeding-purging model. In the period from 1992 to 2017, 22.718 kindling records were analyzed. According to the heuristic approach, the purging coefficient reached the maximum possible value of 0.5 when estimating between 1992 and 1997. Based on the expected fitness over generations and on the expressed opportunity of purging, the beneficial effects of purging could be expected after 10 generations. The proportion of the purged genetic load could be between 20% and 60%. While the results obtained are not entirely conclusive, they do raise the possibility that some of the inbreeding load was caused, at least in part, by genes that could be successfully removed from the population by purging.

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.

Champing at the Bit for Improvements: A Review of Equine Welfare in Equestrian Sports in the United Kingdom

Equestrian sports, including racing (e.g., flat, steeple-chasing, harness or donkey derby); show-jumping; cross-country; dressage; polo; polocrosse; endurance; carriage driving; vaulting and hunting; are hugely popular in the UK, and they involve a significant number of people, both as participants and spectators, and tens of thousands of equids. In this paper, we discuss animal welfare as a complex and disputed issue, clarifying what the term means and how it can be measured. We review many aspects of welfare risk to equids used for sport, addressing issues encountered throughout their lives, including housing, feeding, veterinary intervention, shoeing, handling, training, breeding and equipment. This is followed by a unique exploration of the institutions and social processes influencing equine welfare. The institutional components comprise the rules of competition, the equids, attributes of the stakeholders and the space where participants strive to achieve a common purpose. We endeavour to untangle the most significant elements that create barriers or provide opportunities for equine welfare improvement. We expose the challenges faced by a broad range of stakeholders with differing ethics, attitudes and values. Evidently, there are many welfare risks to which equids used in sports continue to be exposed. It is also evident that significant improvements have occurred in recent times, but there remains a barrier to reducing the risks to an acceptable level. We conclude with recommendations regarding a process for change, involvement of stakeholders and management of knowledge to improve equine welfare that involves identifying and prioritising the risk factors and ultimately leading to interventions, further research and/or education.

Fine-Scale Analysis of Runs of Homozygosity Islands Affecting Fertility in Mares

The loss of genetic variability in livestock populations bred under strict selection processes is a growing concern, as it may lead to increased inbreeding values and lower fertility, as a consequence of the “inbreeding depression” effect. This is particularly important in horses, where inbreeding levels tend to rise as individuals become more and more closely related. In this study, we evaluated the effect of increased inbreeding levels on mare fertility by combining an SNP-based genomic approach using runs of homozygosity and the estimation of genetic breeding values for reproductive traits in a large population of Pura Raza Española mares. Our results showed a negative correlation between whole-genome homozygosity and fertility estimated breeding values (EBVs) at the genome level (ρ = −0.144). However, the analysis at chromosome level revealed a wide variability, with some chromosomes showing higher correlations than others. Interestingly, the correlation was stronger (−0.241) when we repeated the analysis in a reduced dataset including the 10% most and least fertile individuals, where the latter showed an increase in average inbreeding values (F_ROH) of around 30%. We also found 41 genomic regions (ROHi, runs of homozygosity islands) where homozygosity increased 100-fold, 13 of which were significantly associated with fertility after cross-validation. These regions encompassed 17 candidate genes previously related to oocyte and embryo development in several species. Overall, we demonstrated the relationship between increased homozygosis at the genomic level and fertility in mares. Our findings may help to deal with the occurrence of inbreeding depression, as well as further our understanding of the mechanisms underlying fertility in mares.

The effect of the Robertsonian translocation 1/29 on the fertility of beef cattle reared under extensive conditions: a 30 years retrospective study

The Robertsonian translocation 1/29 (rob(1;29)) is the most worldwide widespread chromosomal abnormality in domestic animals. Previous studies have demonstrated its negative effect on fertility in dairy herds, but not in beef cattle extensively bred. In this study, we analyzed the effect of rob(1;29) in a Retinta cattle breed dataset gathered during the last 30 years. The data presented herein includes rob(1;29) analysis of 11,505 cows from 251 herds, pedigree information of 24,790 animals, and 67,457 calving records. Fertility was evaluated using estimated breeding values for the reproductive efficiency (Re), calculated as the percentage ratio between the number of calvings of an individual with the number expected in an optimal situation. Our results showed that cows carrying the heterozygote genotype showed a significant decrease in their Re (-5.10%, p<0.001). No decrease was detected in free rob(1;29) animals and homozygous carriers. In addition, the incidence of rob(1;29) in the breed fertility was decreased to very low values after 30 years of avoiding selection bulls carrier as stallions. The effect of rob(1;29) in cattle fertility is only significant when the prevalence of carrier individuals is high. Selecting against the disease only by the paternal side reduced the incidence to negligible values.

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.

Bayesian inference of the inbreeding load variance for fertility traits in Brown Swiss cattle

Our study investigated the inbreeding load for fertility traits in the Italian Brown Swiss dairy cattle breed. Fertility traits included continuous traits (i.e., interval from calving to first service, days open, and calving interval) and categorical traits (i.e., calving rate at first insemination and nonreturn date at d 56). We included only records of the first 3 parities of cows that calved between 2010 and 2018. We traced up the pedigree of the cows with records as far as possible, ending up with a total of 73,246 animals. The final data set consisted of 59,864 records from 34,921 cows. We analyzed all models using a Bayesian approach that included a covariate with total inbreeding in addition to systematic, permanent environment, additive genetic, and inbreeding load effects. We then evaluated the trends in heritabilities and ratios of the inbreeding load using a continuum of partial inbreeding coefficients from 0.001 to 0.100 as reference. Posterior estimates of heritabilities tended to decrease across the continuum, whereas ratios of the inbreeding load tended to increase, more noticeably in categorical traits (calving rate at first insemination and nonreturn date at d 56). From the results obtained, we confirmed the presence of heterogeneity in inbreeding depression. We then predicted the inbreeding load effects, which had a low reliability of prediction, explained by having only 513 ancestors generating inbreeding. However, reliability of prediction was high enough for some of the individuals, obtaining a favorable prediction of inbreeding load for a relevant percentage, which improved the phenotypic performance of their inbred descendants. These results make it feasible to implement breeding and management strategies that select ancestors with a favorable inbreeding load prediction. In addition, it opens the possibility to define a global index for the expected consequences of the inbreeding generated by each individual.