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

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.

Transmission ratio distortion detection by neutral genetic markers in the Pura Raza Española horse breed

Transmission Ratio Distortion (TRD) is a genetic phenomenon widely demonstrated in several livestock species, but barely in equine species. The TRD occurs when certain genotypes are over- or under-represented in the offspring of a particular mating and can be caused by a variety of factors during gamete formation or during embryonic development. For this study, 126 394 trios consisting of a stallion, mare, and offspring were genotyped using a panel of 17 neutral microsatellite markers recommended by the International Society for Animal Genetics for paternity tests and individual identification. The number of alleles available for each marker ranges from 13 to 18, been 268 the total number of alleles investigated. The TRDscan v.2.0 software was used with the biallelic procedure to identify regions with distorted segregation ratios. After completing the analysis, a total of 12 alleles (out of 11 microsatellites) were identified with decisive evidence for genotypic TRD; 3 and 9 with additive and heterosis patterns, respectively. In addition, 19 alleles (out of 10 microsatellites) were identified displaying allelic TRD. Among them, 14 and 5 were parent-unspecific and stallion-mare-specific TRD. Out of the TRD regions, 24 genes were identified and annotated, predominantly associated with cholesterol metabolism and homeostasis. These genes are often linked to non-specific symptoms like impaired fertility, stunted growth, and compromised overall health. The results suggest a significant impact on the inheritance of certain genetic traits in horses. Further analysis and validation are needed to better understand the TRD impact before the potential implementation in the horse breeding programme strategies.

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.

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.

Estimation of the Genetic Components of (Co)variance and Preliminary Genome-Wide Association Study for Reproductive Efficiency in Retinta Beef Cattle

In this study, we analyzed the variation of reproductive efficiency, estimated as the deviation between the optimal and real parity number of females at each stage of the cow's life, in 12,554 cows belonging to the Retinta Spanish cattle breed, using classical repeatability and random regression models. The results of the analyses using repeatability model and the random regression model suggest that reproductive efficiency is not homogeneous throughout the cow's life. The h² estimate for this model was 0.30, while for the random regression model it increased across the parities, from 0.24 at the first calving to 0.51 at calving number 9. Additionally, we performed a preliminary genome-wide association study for this trait in a population of 252 Retinta cows genotyped using the Axiom Bovine Genotyping v3 Array. The results showed 5 SNPs significantly associated with reproductive efficiency, located in two genomic regions (BTA4 and BTA28). The functional analysis revealed the presence of 5 candidate genes located within these regions, which were previously involved in different aspects related to fertility in cattle and mice models. This new information could give us a better understanding of the genetic architecture of reproductive traits in this species, as well as allow us to accurately select more fertile cows.

Unlocking Horse Y Chromosome Diversity

The present equine genetic variation mirrors the deep influence of intensive breeding programs during the last 200 years. Here, we provide a comprehensive current state of knowledge on the trends and prospects on the variation in the equine male-specific region of the Y chromosome (MSY), which was assembled for the first time in 2018. In comparison with the other 12 mammalian species, horses are now the most represented, with 56 documented MSY genes. However, in contrast to the high variability in mitochondrial DNA observed in many horse breeds from different geographic areas, modern horse populations demonstrate extremely low genetic Y-chromosome diversity. The selective pressures employed by breeders using pedigree data (which are not always error-free) as a predictive tool represent the main cause of this lack of variation in the Y-chromosome. Nevertheless, the detailed phylogenies obtained by recent fine-scaled Y-chromosomal genotyping in many horse breeds worldwide have contributed to addressing the genealogical, forensic, and population questions leading to the reappraisal of the Y-chromosome as a powerful genetic marker to avoid the loss of biodiversity as a result of selective breeding practices, and to better understand the historical development of horse breeds.

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.

Basal Reactivity Evaluated by Infrared Thermography in the "Caballo de Deporte Español" Horse Breed According to Its Coat Color

Horses have been valued for their diversity of coat color since prehistoric times. In particular, the pleiotropic effect that coat color genes have on behavior determines the way the horse perceives and reacts to its environment. The primary aim of this study was to evaluate the influence of coat color on basal reactivity assessed with infrared thermography as eye temperature at rest (ETR), determine their relation with the results obtained by these horses in Show Jumping competitions and to estimate the genetic parameters for this variable to test its suitability for genetic selection. A General Linear Model (GLM) and Duncan post-hoc analysis indicated differences in ETR due to coat color, sex, age, location, and breed-group factors. A Spearman’s rank correlation of 0.11 (p < 0.05) was found with ranking, indicating that less reactive horses were more likely to achieve better rankings. Heritability values ranged from 0.17 to 0.22 and were computed with a model with genetic groups and a model with residual variance heterogeneity. Breeding values were higher with the last genetic model, thus demonstrating the pleiotropic effect of coat color. These results indicate that ETR has a suitable genetic basis to be used in the breeding program to select for basal reactivity due to coat color.

Effects of Selection on Breed Contribution in the Caballo de Deporte Español

The equine breeding industry for sport’s performance has evolved into a fairly profitable economic activity. In particular, the Caballo de Deporte Español (CDE) is bred for different disciplines with a special focus on Show Jumping. The main aim of this study was to determine the effects of 15 years of selection and to find out whether it has been effective and adequate regarding the CDE main breeding objectives. The whole pedigree of 19,045 horses registered as CDE was used, comprising 47,884 animals (18,799 males and 29,085 females). An analysis performed to check for the pedigree completeness level yielded a number of equivalent complete generations (t) equal to 1.95, an average generation interval (GI) of 10.87 years, mean inbreeding coefficient (F) of 0.32%, an average relatedness coefficient (AR) of 0.09% and an effective population size (Ne) of 204. For the analyses, animals were divided into fourteen breed groups. Additionally, in order to study the evolution of these breeds over time and their influence on CDE pedigree, five different periods were considered according to the year of birth of the animals. Performance data used in the analyses were the estimated breeding values (EBV) of the Show Jumping sport discipline of 12,197 horses in the CDE pedigree, available from the 2020 routine genetic evaluations of the CDE breeding program (starting in 2004). Dressage and Eventing EBV values were also assessed. Results showed values of F higher than expected under random mating; this pointed to some degree of inbred matings. With regard to the evolution of breeding values, we found that, in general, EBVs of offspring were higher than the EBVs of parents. Notwithstanding, there is still a need for improvement in population management and the coordination of the breeders to get higher responses but controlling the loss of genetic diversity in the CDE breed.

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.