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

Parent of origin genetic effects on milk production traits in a population of Iranian Holstein cows

The aim was to estimate the relative contribution of imprinting effects from both paternal and maternal sides to phenotypic variation in milk production traits including 305 days milk yield (MY), average daily milk production (ADM), fat percentage (F%), protein percentage (P%), 305 days fat yield (FY), 305 days protein yield (PY), ratio of fat percentage to protein percentage (F:P) and somatic cell score (SCS) in Iranian Holstein cows. To do this, each trait was analysed with a series of four animal models, which were identical for fixed and additive genetic effects but differed for combinations of paternal and maternal imprinting effects. The log-likelihood ratio test (LRT) and Akaike's information criteria (AIC) were used to select the best model for each trait. Correlations between traits due to additive and imprinting effects were estimated by bivariate analyses. For all traits studied, fitting the imprinting effect led to a better data fit. Also, it resulted in a noticeable decrease in additive genetic variance from 8% (SCS) to 28% (F:P). A significant maternal imprinting effect was detected on all traits studied. Estimates of maternal imprinting heritability (h mi 2 $$ {h}_{\mathrm{mi}}^2 $$ ) were 0.07 ± 0.02, 0.04 ± 0.01, 0.06 ± 0.01, 0.05 ± 0.01, 0.5 ± 0.01, 0.09 ± 0.02, 0.07 ± 0.02 and 0.06 ± 0.01 for MY, ADM, F%, P%, FY, PY, F:P and SCS, respectively. For F:P, in addition to the maternal imprinting effect, a significant paternal imprinting component was also detected with a 7% contribution to phenotypic variance of F:P. Estimates of direct heritability (h a 2 $$ {h}_{\mathrm{a}}^2 $$ ) were 0.29 ± 0.02, 0.17 ± 0.01, 0.22 ± 0.02, 0.11 ± 0.01, 0.18 ± 0.02, 0.22 ± 0.02, 0.15 ± 0.04 and 0.06 ± 0.01 for MY, ADM, F%, P%, FY, PY, F:P and SCS, respectively. Maternal imprinting correlations (rmi) were in a wide range between -0.75 ± 0.15 (P%-SCS) and 0.95 ± 0.11 (MY-ADM). Additive genetic correlations (ra) ranged between -0.54 ± 0.05 (MY-P%) and 0.99 ± 0.01 (MY-ADM) and phenotypic correlations (rp) ranged from -0.30 ± 0.01 (MY-F%) to 0.93 ± 0.01 (MY-ADM). The Spearman's correlation between additive breeding values including and excluding imprinting effects deviated from unity especially for top-ranked animals implying re-ranking of top animals following the inclusion of imprinting effects in the model. Since including imprinting effects in the model resulted in better data fit and re-ranking of top animals, including these effects in the genetic evaluation models for milk production traits was recommended.

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.