Phenotypic and genetic analysis of reproductive traits in horse populations with different breeding purposes

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.

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.

Genomic screening of allelic and genotypic transmission ratio distortion in horse

The phenomenon in which the expected Mendelian inheritance is altered is known as transmission ratio distortion (TRD). The TRD analysis relies on the study of the transmission of one of the two alleles from a heterozygous parent to the offspring. These distortions are due to biological mechanisms affecting gametogenesis, embryo development and/or postnatal viability, among others. In this study, TRD phenomenon was characterized in horses using SNP-by-SNP model by TRDscan v.2.0 software. A total of 1,041 Pura Raza Español breed horses were genotyped with 554,634 SNPs. Among them, 277 horses genotyped in trios (stallion-mare-offspring) were used to perform the TRD analysis. Our results revealed 140 and 42 SNPs with allelic and genotypic patterns, respectively. Among them, 63 displayed stallion-TRD and 41 exhibited mare-TRD, while 36 SNPs showed overall TRD. In addition, 42 SNPs exhibited heterosis pattern. Functional analyses revealed that the annotated genes located within the TRD regions identified were associated with biological processes and molecular functions related to spermatogenesis, oocyte division, embryonic development, and hormonal activity. A total of 10 functional candidate genes related to fertility were found. To our knowledge, this is the most extensive study performed to evaluate the presence of alleles and functional candidate genes with transmission ratio distortion affecting reproductive performance in the domestic horse.

Reproductive Performance of Tunisian Arabian Stallions: A Study on the Variance and Estimation of Heritability

A consistently high level of stallion fertility plays an economically important role in modern horse breeding. To better understand the factors affecting variation in stallion fertility, we have performed a statistical analysis study on some breeding factors: year of breeding, stud farm, age of the stallion, number of covered mares per stallion, reproduction methods, and age of the mare. This work was conducted on 94 purebred Arabian stallions in four different regions of Tunisia. The results showed an increase in the number of stallions during the study period, ranging from 11.33% in 2011 to 13.29% in 2018. Sidi Thabet’s stud farm contained the highest number of purebred Arabian stallions. The majority of stallions were between 15 and 21 years old and had covered 1 to 20 mares; 95.19% of stallions were used in natural mating (Nat); 50.36% had low fertility, 17.69% had medium fertility, and 32.3% had excellent fertility according to fertility standards. Depending on the year and stud, there was a variation in fertility per cycle (FERPCE) and end-of-season fertility (FERPSE) of the stallions. The highest average FERPCE and FERPSE values were obtained using artificial insemination with fresh semen (AIF). Analysis of FERPCE and FERPSE showed that the model used in our study explained 40.21% of total variability observations for FERPCE and 42.1% for FERPSE. The used statistical model showed that the breeding year, the stud, the age of the stallion, the number of covered mares by stallions and the method of reproduction significantly affected both FERPCE and FERPSE (with p = 0.001). Low to moderate heritability estimations for FERPCE (hs2 = 0.08) and FERPSE (hes2 = 0.36) were obtained.

A genome-wide association study of mare fertility in the Pura Raza Español horse

Despite the economic importance of fertility for the horse industry, few efforts have been made to achieve a better understanding of the genetic mechanisms underlying its control. This is probably due to the difficulty of obtaining reliable phenotypes and the complexity of modelling the environmental and management factors. This work is novel in that we propose to use reproductive efficiency (RE) as an indicator of mare fertility. To achieve this, we performed a genome-wide association study in the Pura Raza Español horse aimed at identifying genomic variants, regions, and candidate genes associated with fertility in mares. The dataset included 819 animals genotyped with the Affymetrix Axiom™ Equine 670 K single-nucleotide polymorphisms (SNPs) Genotyping Array and the deregressed breeding values for RE trait, obtained using a ssBLUP model, employed as pseudo-phenotypic data. Our results showed 28 SNPs potentially associated with RE, which explained 87.19% of the genetic variance and 6.61% of the phenotypic variance. Those results were further validated in BayesB, showing a correlation between observed and predicted RE of 0.57. In addition, 15 candidate genes (HTRA3, SPIRE1, APOE, ERCC1, FOXA3, NECTIN-2, KLC3, RSPH6A, PDPK1, MEIOB, PAQR4, NM3, PKD1, PRSS21, IFT140) previously related to fertility in mammals were associated with the markers and genomic regions significantly associated with RE. To our knowledge, this is the first genome-wide association study performed on mare fertility.

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.

Systematic Analysis of Breed, Methodological, and Geographical Impact on Equine Sperm Progressive Motility

Simple Summary

Semen quality is an important indicator of reproductive health and fertility. With adverse temporal trends in human semen quality over the past 50 years paralleled in male animals, there is increasing concern for the causes and implications of perturbed male reproductive health. The evaluation of equine progressive motility (PM), a parameter closely associated with fertility, provides information on the fertilising capacity of equine ejaculate and current reproductive health of the equine stallion. Using systematic analysis, recent trends in equine PM were determined from 696 estimates from 280 individual studies. Temporal trends indicate equine PM has not significantly changed between the years 1990 and 2018. Significant breed, methodological, and geographical variations observed in equine PM may considerably influence actual and reported stallion fertility potential. Information on stallion PM meaningfully contributes to the wider literature on semen quality and provides avenues for future stallion fertility research. This systematic analysis presents the wider challenges associated with semen quality assessment, particularly within the equine sector, and provides recommendations to promote consistency across industry and research.

Abstract

Over the past five decades, there has been increasing evidence to indicate global declines in human semen quality. Parallel adverse trends measured in male animals indicate a potential environmental aetiology. This study evaluated the progressive motility (PM) of stallion ejaculate through a systematic review and meta-analysis. A total of 696 estimates of equine PM from 280 studies, which collected semen samples between the years 1990 and 2018, were collated for meta-analysis. The method of motility analysis, breed, season of collection, and geographical location were extracted. Simple linear regression determined temporal trends in stallion PM. Studies using microscopy estimated PM to be significantly greater compared to computer-automated methods (p ≤ 0.001). For Arabian breeds, PM was consistently higher than other breeds. Over time, there was a significant decline in PM for studies from Europe (n = 267) but a significant increase for studies from North America (n = 259). Temporal trends indicate the fertilising capacity of equine ejaculate has remained consistently high in the last three decades. That being so, variations observed suggest methodological, geographical, and individual stallion differences may significantly influence actual and reported stallion fertility potential.

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.

Breeding Strategies to Optimize Effective Population Size in Low Census Captive Populations: The Case of Gazella cuvieri

Simple Summary

Small-sized populations frequently undergo a significant loss of genetic variability that can lead to their extinction. Therefore, research on animal breeding has focused on mating systems for minimizing the disappearance of genetic variability. Minimizing the average coancestry of offspring has been described as the best strategy for this purpose. Traditionally, the preservation of genetic variability has been approached via breeding strategies for increasing the effective population size (Ne). The main objective of this study was to compare, via computer simulations, the performance of different breeding schemes to limit the losses of genetic diversity in small populations. This objective was achieved by monitoring the evolution of the effective size obtained by different strategies across 20 generations with a starting point of two pedigree real populations of Gazella cuvieri. The results showed that minimizing average coancestry in a cohort did not maximize the effective size as compared with new strategies that were designed for this study. Furthermore, the best strategy may vary for each population and should be studied individually.

Abstract

Small-sized animal populations can undergo significant loss of genetic variability that can lead to their extinction. Therefore, studies on animal breeding have focused on mating systems for minimizing the disappearance of genetic variability. The main objective of this study was to compare, using computer simulations, the performance of different breeding schemes to limit the loss of genetic diversity in small-sized populations. This objective was achieved by monitoring the evolution of the effective population size obtained by 23 strategies throughout 20 generations in two populations of Gazella cuvieri. The scenarios were designed with different assumptions, in both reference subpopulations, regarding: the use of parents coancestry or offspring coancestry, the use of their increases or the coefficients themselves, and the number of males and females involved. Computations were performed using an experimental module of Endog v4.9 developed for this purpose. The results of the study showed that strategies for minimizing the coancestry of the parents were better in the short term; however, these strategies were worse in the long term. Minimizing the average coancestry of the offspring was a better approach in the long term. Nevertheless, in both populations, the best results were obtained when both the coancestry of the parents and the coancestry of the offspring were weighted at 5% each and neither males nor females were assumed to contribute to the next generation. In any case, not all strategies had the same evolutionary pattern throughout generations in both populations. The current results show that neither traditional nor new strategies have any general use. Therefore, it is important to carefully test these strategies before applying them to different populations with different breeding needs under different conditions, such as different generation intervals, and different natural breeding systems such as monogamy or polygyny.

Impaired Reproductive Function in Equines: From Genetics to Genomics

Fertility is one of the key factors in the economic and productive success of the equine industry. Despite this, studies on the genetic causes affecting reproductive performance are scarce, especially in mares, where the genetic architecture of the reproductive traits is extremely complex. Today, with the increasing availability of new genomic methodologies for this species, we are presented with an interesting opportunity to understand the genetic basis of equine reproductive disorders. These include, among others, novel techniques for detecting chromosomal abnormalities, whose association with infertility in horses was established over 50 years ago; new sequencing technologies permitting an accurate detection of point mutations influencing fertility, as well as the study of inbreeding and molecular homozygosity, which has been widely suggested as one of the main causes of low reproductive performance in horses. Finally, over the last few years, reproductive performance has also been associated with copy number variants and candidate genes detected by genome-wide association studies on fertility traits. However, such studies are still scarce, probably because they depend on the existence of large and accurate phenotypic datasets of reproductive and/or fertility traits, which are still difficult to obtain in equines.

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.