1. Introduction
Of the countries in the European Union, Spain has the second-highest number of dairy sheep, with a population of up to 16 million in 2018 [
1]. Dairy sheep farms are an essential factor in the maintenance of rural development through environmentally friendly and sustainable production systems, which provide a broad diversity of high-quality products, including artisanal mature cheese and suckling lamb meat [
2]. For many years, sheep milk production in Spain has been focused on the extensive exploitation of local sheep breeds, such as Churra, Latxa (which consists of two varieties, Cara Negra and Cara Rubia) and Manchega, which are highly adapted to the different local conditions, maximizing the use of local resources. In the last 20 years, however, some foreign specialized dairy sheep breeds, such as Assaf, Awassi, and Lacaune, are, to some extent, replacing the local breeds in semi-intensive and intensive farms [
3] because of their higher milk yield potential. In particular, Assaf has the most significant impact on the Spanish and Portuguese dairy sheep industry [
4].
For all the above-mentioned breeds, the genetic improvement for milk production is performed by the corresponding breeder association through official breeding programs, in which artificial insemination (AI) becomes a fundamental tool of progeny testing. AI plays a significant role because it is directly related to fertility and is the primary reproductive tool used for genetic improvement by introducing genes from superior sires into flocks [
5]. Hence, AI is one of the most powerful biotechnological reproductive methods used in most breeding programs of livestock species.
In relation to AI, it is important to study the characteristics of semen production and semen quality. The number of doses produced per ram ejaculate depends on the volume, sperm concentration, and motility of the sperm, which are mainly affected by environmental and genetic factors [
6].
To increase the efficiency of the dairy sheep AI centers, studies focused on the assessment of the environmental factors affecting semen production and the estimation of genetic parameters for seminal production and quality traits are required. In cattle, different authors have reported moderate heritabilities for some semen traits [
5,
7], suggesting the potential for genetic improvement not only for male reproduction traits but also for female reproductive traits, which are generally difficult to improve through direct selection [
8].
In sheep, the limited studies addressing the estimation of genetic parameters for ram semen traits in different breeds show a wide variation [
6,
9,
10], highlighting the importance of studying these traits for individual breeds. The estimation of genetic parameters is relevant not only to assess the possibility of the genetic improvement of semen traits but also to understand the interactions among the different traits and, therefore, help define the most efficient selection strategy and avoid undesired consequences of selection [
5]. Finally, the inclusion of seminal traits in the male selection indices could be an aspect to consider and investigate in sheep breeds.
Given these considerations, the present study aimed to estimate genetic parameters for traits related to semen production and quality in Spanish dairy sheep breeds. For this purpose, four local breeds (Churra, Manchega, Latxa Cara Negra, and Latxa Cara Rubia) and one foreign highly productive breed (Assaf) were analyzed. Although the analyses here presented have been performed individually for each breed, this study provides an opportunity to assess the genetic parameters for semen traits across different breeds.
3. Results
The descriptive statistics of the semen traits for the five breeds are shown in
Table 1. The ASS breed had the highest mean value for two of the three traits studied, with values of 1.36 mL and 4.87 for the VOL and MOT traits, respectively. The ASS breed had the most extensive standard deviations for the VOL (0.54) and SC (1320.50) traits among the five breeds. The MAN breed showed the highest mean value for the SC (4513.80) trait and the largest standard deviation for the MOT trait (0.44). The LCR and LCN breeds presented very similar results for each of the studied traits. The phenotypic distributions observed for each trait and breed are provided in
Supplementary Material (Figures S1–S3). As can be observed, the distributions for VOL and SC were very similar across breeds and close to normality, whereas the MOT distributions were not normal.
The estimates of the genetic, permanent environmental, and residual variances and repeatability are presented in
Table 2. In general, the additive genetic variances estimated for the seminal traits were of similar magnitude among breeds, except for those estimated for the LCR breed, for all the traits, and the LCN breed, for SC and MOT traits, which were remarkably lower. The ASS breed showed the highest additive genetic variance for the VOL (0.030) and SC (281,870) traits, whereas the highest additive genetic variance for the MOT trait was found in CHU (0.017). The LCN and LCR breeds obtained the lowest additive genetic variances for all the traits. The repeatability estimates were similar among the breeds for VOL and SC, ranging between 0.211–0.291 for VOL and 0.270–0.303 for SC, whereas for MOT, the estimates obtained for the LCN and LCR (0.096 and 0.077, respectively) and MAN (0.079) breeds were notably lower than those estimated for ASS and CHU (0.197 and 0.151, respectively).
Regarding the genetic parameters presented in
Table 3 and considering the three traits under study, we observed that the ASS, LCR, and MAN breeds showed similar heritability estimates (
Table 3), with SC being the most heritable trait (0.20, 0.11, and 0.15, respectively). The VOL trait showed intermediate heritability estimates (0.12, 0.08, and 0.11, respectively), and the MOT trait showed estimates close to zero (0.03, 0.01, and 0.03, respectively). In the CHU breed, the SC trait showed similar heritability estimate (0.18) to those found in ASS and MAN, whereas the estimates obtained for VOL and MOT were substantially higher than those found in the other breeds (0.20 and 0.11, respectively). The heritability estimate for the VOL trait in CHU was the highest heritability reported in this work. Finally, the LCN breed showed similar trends of heritability estimates to those of CHU, although with lower values, with the VOL trait showing the highest heritability value (0.16) within this breed.
Among the genetic correlations estimated in this work (
Table 3), only those estimated for the MAN breed showed acceptable accuracy, as very large standard errors (SE) were observed for the rest of breeds estimates. Focusing on MAN, moderate and negative correlations between the VOL trait and the two other traits were observed (−0.49 and −0.29 with SC and MOT, respectively), whereas a moderate positive correlation was estimated between the SC and MOT traits (0.32).
In the ASS and CHU breeds, moderate and positive phenotypic correlations were found for the SC-VOL and SC-MOT trait pairs (with a range from 0.16 to 0.23), with close to zero negative estimates between MOT and VOL; in the LCN and LCR breeds, the MOT-SC phenotypic correlation was positive (0.16 and 0.18, respectively), whereas all the phenotypic correlations in MAN were low (
Table 3).
4. Discussion
Dairy sheep are traditionally reared in Mediterranean areas on a wide range of production systems. In the last 30 years, dairy sheep systems have kept their dual purpose and have maintained an income characterized by a 65–75% contribution from milk and 25–35% from meat. The efficiency of selection for milk traits is now whetting interest in new breeding goals, which are relevant to satisfying consumer demand with other products, such as quality-labeled cheeses, and “healthy” food products [
3].
Milk yield is one of the most important economic traits in dairy sheep, but in practice (for genetic improvement), fertility is also a critical trait for the economic efficiency of dairy sheep flocks [
14]. According to different studies in other species, semen quantity and quality are the main factors affecting the fertility of rams used for AI [
15,
16]. AI is becoming increasingly essential within sheep breeding programs, allowing the dissemination of the genetic gain to the whole population. Reproductive efficiency has become a key factor to consider because improved reproductive success could result in greater genetic progress for selected traits [
17]. Therefore, the selection of rams should be focused not only on milk production traits but also on semen traits in which the estimation of genetic parameters is a required step.
The dataset compiled through this study has provided the opportunity to analyze different seminal traits rarely studied in the ovine species and to compare phenotypic and genetic parameters among five of the most important Spanish dairy sheep breeds. Traits such as volume, sperm concentration, and sperm motility are relevant traits to analyze in the ejaculate of rams.
In general, the mean values reported here for these traits were very similar to those obtained in other sheep breeds [
6,
18]. The mean value observed for the MOT trait in the MAN breed (3.61) was low compared with that of the rest of breeds analyzed here and also stood out compared with French breeds [
6]. We think that the deviation from the normality and the subjectivity of the MOT has affected the accuracy of the genetic parameters estimated here for this trait, as will be discussed below.
Seminal traits are affected by environmental, management, and genetic effects and physiological status [
6]. In our case, three AI centers located in different regions of Spain provided semen production data, and the main environmental factors affecting sperm production were collected at each center. In this study, because the analyses were performed individually for each breed and the different centers analyzed different breeds, the reported differences might be due not only to the breed factor but also to the different criteria to consider an ejaculate suitable for AI at each AI center. For example, for the MAN breed, a sequential evaluation of ejaculates is carried out in such a way that a minimum volume is required for an ejaculate; those ejaculates that pass this criterion are evaluated for motility, and once they pass this criterion, the concentration is measured. Therefore, all the ejaculates with concentration data will have volume and motility data above the respective thresholds. On the other hand, the AI center that analyzed the samples from LCN and LCR breeds implements more restrictive criteria to declare an ejaculate suitable for AI (a mass motility value of 4 and 0.3 for volume are required as minimum values). Hence, these breeds show a relatively low variability compared with the rest of breeds, especially for the SC and MOT traits. In the case of the breeds analyzed in the same AI centers, CHU-ASS and LCN-LCR, the differences reported between them are more likely to reflect genuine differences due to the breed effect. In any case, we acknowledge that the use of these quality thresholds can introduce certain biases on the genetic parameter estimates reported here (e.g., a lower heritability due to the artificially reduced phenotypic variance). Hence, the ideal situation for this kind of study would be that all semen quality measures were recorded without any censorship in all the AI centers; but this state is unrealistic in practical conditions. The impact of the thresholds will depend on the proportion of samples that are below the corresponding threshold in each breed and also the change in the variable distributions when including or discarding some of the data.
The effects of flock and the interaction of collection season–year were significant for all the traits and breeds. Our results are in agreement with those obtained by Anel et al. [
13] in the Churra breed, where the farm effect was found to have a major influence on ewes’ fertility. The season and year of collection were also found to be significant environmental factors, probably related to the seasonal reproductive behavior of the ovine species. As expected, all breeds analyzed in this work showed the highest values for the three studied traits at the favorable reproductive season (i.e., the months from September to February). This agrees with observations reported in other sheep breeds [
19]. These two factors, whose relevance to semen traits has been highlighted by other authors, involve many environmental factors, such as animal care, climatic changes, temperature, feed quality, and length of day [
20,
21]. Another critical effect influencing semen production and quality detected here was the age of the male at the collection time. This factor was significant for all the traits and breeds, except for VOL in ASS and CHU. In this case, the age factor showed a variable effect on the studied traits depending on the breed. For example, the highest values for the three studied traits were observed, for the ASS, CHU, and MAN breeds, in rams with age between 37 and 60 months; and for the LCN and LCR breeds, for rams included in two different age groups (10–24 months and 49–60 months). The common effect of age across all the breeds was an age-related decrease in semen trait values (rams with more than 60 months), which is in agreement with observations reported in Lacaune and Manech Tête Rousse breeds [
6]. Additionally, in Austrian Simmental bulls, a decrease in semen quality traits, without effect on the VOL trait, has also been reported to be associated with increased age [
22]. Finally, both the mounting effect (in ASS, CHU, LCN, and LCR breeds) and mounting regime effect (only in MAN breed) were significant for all the traits. In our analysis, a maximum of two ejaculates per animal was collected every 2–5 min (except in the MAN breed, in which only one ejaculate per animal was routinely collected). Comparing the results obtained here in the different mountings, the mean values for VOL and SC traits in mount 2 were, in general, lower than in mount 1, whereas slightly higher MOT mean values were observed in mount 2 compared to mount 1. David et al. [
6] also found a significant influence of the mounting effect in the Lacaune and Manech Tête Rousse breeds, reporting higher motility and a lower number of spermatozoa with a shorter interval between two collections. They also suggested that there could be an optimum interval, which ensures a higher quantity of highly motile spermatozoa, but this question needs to be further investigated. For the MAN breed, with only one mount per collection, the mounting regime effect, measured as the interval in days between two consecutive semen collections, was important. Thus, an interval of 3 days was found to be optimum for the three semen traits assessed in this study. Lower intervals (1–2 days) resulted in a decrease in semen production and quality, whereas intervals above 3 days did not increase the trait values.
The highest additive genetic variances were found for the VOL and SC traits in the ASS breed and the MOT trait in the CHU breed. The large additive genetic variances observed in ASS can be explained by the fact that this breed has been recently created by the absorption of several Spanish breeds, the large current census of the Assaf population and its recent breeding program (initiated in 2005), compared with that of the other studied breeds.
The coefficients of repeatability describing the degree of similarity between consecutive measurements affected by random environmental factors were also estimated. The SC trait showed the highest repeatability value in all the breeds studied. The repeatability values found in this work were lower than those estimated in other sheep breeds [
6] and bovine breeds [
20,
21]. The low to moderate repeatability values obtained for the analyzed semen traits suggest that first collection records are not good indicators of the future reproductive yield of AI rams.
The heritability estimates obtained in this study are similar to those described in the literature. For example, David et al. [
6] reported similar results in two French sheep breeds ranging from 0.07 (for the MOT trait) to 0.27 (for the VOL and SC traits) in the Lacaune breed and from 0.02 (for the MOT trait) to 0.33 (for the VOL trait) in the Manech Tête Rousse breed. We think the deviation from the normality observed for the MOT trait, and the subjectivity linked to the measurements of this trait, can explain why this is the trait for which the lowest heritability estimates were observed (four out of the five breeds had heritability estimates lower than 0.04, only CHU reached an estimate higher than 0.1). Our data are also comparable to those published for other species. For example, Li et al. [
16] reported similar heritabilities in pigs (Landrace boars), ranging from 0.11 to 0.23 (for the MOT and SC traits, respectively), and Karoui et al. [
7], showed a similar pattern, obtaining heritabilities from 0.16 to 0.22 (for the MOT and VOL traits, respectively) in Holstein bulls. In contrast, other researchers found higher heritabilities for the MOT trait, surpassing those found in the VOL or SC traits. For example, in Ethiopian sheep, Rege et al. [
9] found values of 0.32 and 0.27 for the MOT trait in 9 and 12-month-old rams, respectively. For Holstein bulls, Druet et al. [
5] obtained a heritability of 0.43, and Li et al. [
16], found values of 0.26 and 0.42 in Yorkshire and Duroc pig breeds, respectively. Heritability estimates found in the present study and the variability found in different breeds and species point towards the opportunity for genetic improvement of some semen traits through selective breeding. In any case, from a practical point of view, before considering the improvement of these traits through selection, it would be necessary to estimate the genetic correlations between the semen traits and the current selection objectives of the breeding programs of these breeds, which are mainly focused on dairy-related traits.
Understanding the relationship between the different seminal traits is very informative for the efficient development of selection strategies aimed at the genetic improvement of the considered target traits. Generally, similar values of phenotypic correlations to those estimated here for VOL, SC, and MOT were found in the literature [
6]. The low phenotypic correlations obtained for some of the breeds in this and other studies [
9] may indicate that some environmental factors, particularly when considering a wide range of time, influencing the studied traits are not controllable by the AI centers (e.g., seasonality temperature, dryness). It should also be considered that for some of these environmental effects, for example, related to the mounting regimen, ovine AI centers may implement a more variable routine than other species AI centers (e.g., cattle, pig).
This influence of non-controllable environmental factors might explain the high standard error of the genetic correlations reported for most of the breeds in this study. Only for the MAN breed did genetic correlation estimates among the semen traits show adequate accuracy (low SE). This may be due to the higher number of animals analyzed in the MAN breed (1200 vs. 500 in the other breeds) and to a more suitable family structure for the determination of genetic relationships. These are key factors influencing the precision of the genetic correlation estimates [
23]. Hence, focusing on the MAN estimates, the moderate negative genetic correlation (−0.49) observed between VOL and SC agrees with the average estimates reported, for both young and adult animals, for this pair of traits in Lacaune and Manech Tête Rousse breeds [
6] and Ethiopian highland sheep [
9]. This agrees with the expected relationship between traits of this nature (volume/concentration), which would determine the need for considering both traits when trying to improve semen production and quality in rams. For the other two trait pairs involving the MOT trait, which showed a non-normal distribution, the genetic correlations showed remarkably high SE for MAN, and also for the other breeds. When comparing the published genetic correlations for SC-MOT and VOL-MOT, we found important discrepancies among studies, even different for breeds and age groups analyzed in the same study [
6,
9]. Whereas for the VOL-MOT pair, the genetic correlation estimate reported here for MAN (−0.29) shows the same direction but higher magnitude than those reported in French breeds (range: −0.03 to −0.09) [
6], this estimate clearly differs from that reported in Ethiopian sheep (0.51) [
9]. The moderate positive genetic correlation that we found for the SC-MOT pair (0.32) agrees with the estimates reported for young animals of the Manech Tête Rousse (0.44) [
6] and Ethiopian highland breeds (0.58) [
9], but not in the Lacune breed (0.04–0.07) [
6]. All these discrepancies are likely to be explained by the high level of subjectivity when measuring the MOT trait in AI centers (as shown in
Figures S1–S3), which is also directly related to the low heritability estimates reported here for this trait, not only in MAN but in all the studied breeds. Because of the relevance of the MOT trait in sheep fertility [
24], we consider that if the moderate positive correlation reported here and by Davis et al. [
6] for the SC-MOT pair (0.32 and 0.44, respectively) is confirmed in MAN by future studies, the MOT trait could be improved by indirect selection using the SC traits as a proxy.