Genetic evaluation for semen characteristics in a crossbreeding project involving Saudi and Spanish V-line rabbits

Abstract

Two parallel schemes of crosses were practised involving the Spanish maternal line called V-line (V) and Saudi Gabali (S) rabbits. The first scheme began by crossing Saudi Gabali bucks with V-line does to get the F₁ cross (½S½V), then does of this F₁ were backcrossed with bucks of V-line to get ¾V¼S, then the progenies of the backcross were inter se mated to get (¾V¼S)², followed by two more generations of inter se mating of (¾V¼S)² to get a new synthetic maternal line named Saudi 2. The second scheme began by crossing V-line bucks with Saudi does to get the F₁ cross (½V½S), then does of this F₁ cross were backcrossed with Saudi bucks to get ½S½V, then progeny of this backcross were inter se mated to get (¾S¼V)², followed by two generations of inter se mating of (¾S¼V)² to form a new synthetic paternal line named Saudi 3. A generalised least-square procedure was used to estimate differences between lines in direct and maternal additive effects, direct and maternal heterosis, and direct recombination losses. A total number of 2497 ejaculates of 642 bucks were evaluated for ejaculate volume (VE), pH of semen, sperm concentration (SC), percentage of motile spermatozoa (MS), percentage of spermatozoa with abnormal form (AS), percentage of dead spermatozoa (DS) and libido score (LS). The estimates of direct additive effect for DS (14.2%), SC (12.4%) and AS (12.0%) were moderate and favourable to Saudi rabbits. Maternal additive effects for VE (11.7%) were in favour of genes of V-line carried by the dams of bucks, while these estimates were in favour of Saudi genes for DS (18.6%), AS (13.2%), MS (10.8%) and SC (10.0%). Positive estimates of direct heterosis recorded for SC (13.5%), VE (10.6%) and MS (10.5%) and the negative estimates recorded for AS (−21.5%) and DS (−20.3%) were moderate and well disposed to semen parameters of crossbred bucks. Positive estimates of maternal heterosis for VE (24.0%), MS (21.8%), SC (10.3%) and LS (10.2%) of bucks, and the negative estimates for percentage of DS (−14.7%) and AS (−9.6%) were also significantly favourable. Estimates of direct recombination losses for the majority of semen traits were always not significant.

Introduction

Semen characteristics are very important if we consider that one single male affects the fertility and prolificacy of many females, especially when artificial insemination (AI) is performed as a routine in rabbit farms (Alvariño, 2000). On the other hand, semen characteristics are important not only due to their relationship with male fertility and prolificacy (which are, in general, good) but also because they determine the number of doses that can be obtained from one ejaculate, and therefore, the efficacy of using bucks for AI. The objective of the AI centres is to obtain a greater number of doses at lower cost, and at the same time it should allow semen production to be increased while maintaining a high level of male fertility and prolificacy. For rabbits to be raised in large-scale commercial production in hot climate areas, libido of bucks could affect negatively both conception and kindling rates. Therefore, bucks must be tested for libido behaviour in order to be used efficiently in any breeding programme.

Results on semen quantity and quality of crossbred males from sire lines could lead to a change in rabbit meat production schemes by using crossbred males in the terminal cross instead of purebred males (Egea De Prado and Roy Perez, 1992; Brun et al., 2002; Al-Sobayil and Khalil, 2002). Alvariño et al. (1996) reported that differences in semen parameters among paternal breeds could affect conception rates. Recently, García-Tomás et al. (2006a and b) indicated that the superiority of crossbred bucks has not been proved for Spanish sire lines and semen traits.

Vicente et al. (2000) in Spain observed considerable differences in seminal characteristics among males of rabbit lines. Brun et al. (2002) in France have reported significant heterotic effects on sperm concentration (SC), mass motility and percentage of motile sperm per ejaculate. Also, García-Tomás et al. (2006a) in Spain found a significant but negligible percentage of heterosis for sperm normalcy, a high percentage of heterosis for the two traits related to the presence of proximal and distal cytoplasmatic droplets in the spermatozoa and a non-significant heterosis effect for pH. To date, publications concerning genetic analysis of semen parameters for crossbred rabbits raised in hot climate countries are very scarce. The objective of the present study was mainly to estimate direct and maternal heterosis, direct and maternal genetic additive effects, and recombination losses for semen characteristics of bucks in a crossbreeding project carried out in a hot climate country involving Spanish V-line (V) and Saudi Gabali (S) rabbits, with the intention to synthesise new maternal and paternal lines suitable for hot climates.

Material and methods

Animals and breeding plan

The rabbits used in this study represent one Saudi breed (Saudi Gabali, S) and one exotic breed (Spanish V-line, V). Saudi Gabali is a desert breed raised under the Saudi desert conditions especially in the Najd area and rabbits of this breed are characterised by a large litter size of six to eight young, a heavy mature body weight of 3.2 to 3.8 kg and the ability to survive and adapt to produce and reproduce in a hot environment. The V-line is a maternal rabbit line selected for number of young weaned per litter (Estany et al., 1989) for 21 generations. The size of the line is around 120 does and 25 bucks.

Eighty pedigree does and 16 pedigree bucks of V-line rabbits were imported from Universidad Politécnica de Valencia, Valencia, Spain, in September 2000. A 4-year crossbreeding project involving these two breeds was started in September 2000 in the experimental rabbitry, College of Agriculture and Veterinary Medicine, El-Qassim region to develop a maternal line and a paternal line of rabbits in Saudi Arabia.

The breeding plan in the project permitted connected production of 11 genetic groups as shown in Table 1. Two parallel schemes of crosses were undertaken. The first scheme began by crossing Saudi Gabali bucks with V-line does to get the F₁ cross (½S½V), then does of this F₁ were backcrossed with bucks of V-line to get ¾V¼S, then progeny of the backcross were inter se mated to get (¾V¼S)² and, finally, followed by two generations of inter se mating of the previous progeny (¾V¼S)² to get a new synthetic maternal line named Saudi 2. The second scheme began by crossing V-line bucks with Saudi does to get the F₁ cross (½V½S), then does of this F₁ were backcrossed with Saudi bucks to get ¾S¼V, then progeny of this backcross were inter se mated to get (¾S¼V)² and, finally, followed by two inter se mating generations of (¾S¼V)² of the previous progeny to form a new synthetic paternal line named Saudi 3. The bucks were randomly assigned to mate the does naturally with the restriction to avoid matings of animals with common grandparents.

Table 1 Genetic group of the bucks and their parents and coefficients of the matrix which relate the EFs with the genetic crossbreeding parameters

EF = estimable function, V = V-line, S = Saudi Gabali.

^†G_i solution for the ith genetic group of buck.

^‡Defined as the difference between direct (maternal) additive effects for the V-line and Saudi Gabali breed.

Housing and feeding

Rabbits were raised in a semi-closed rabbitry. Breeding does and bucks were housed separately in individual wire cages. All cages are equipped with feeding hoppers and drinking nipples. In the rabbitry, the environmental conditions were monitored; temperature ranged from 20°C to about 32°C, the relative humidity ranged from 20% to 50% and the photoperiod was 16L : 8D. Young rabbits were weaned at 4 weeks of age. Rabbits were fed a commercial pelleted diet during the whole period. On a dry-matter basis, the diet contained 18.5% crude protein, 8.0% crude fibre, 3.0% ether extract and 6.5% ash. Feed and water were available ad libitum.

Semen collection and libido score

Data collected for semen characteristics were available for bucks of all generations. The distribution of bucks used and their ejaculates produced in each genetic group are presented in Table 2. A total number of 2497 ejaculates collected from 642 bucks were evaluated for semen characteristics. Bucks were housed in individual cages and they were accustomed to ejaculate in an artificial vagina. The semen samples were collected from the bucks at an early age of 5 to 6 months. Three to five ejaculates per buck were collected in the morning with a collection interval of 2 weeks. Once the ejaculate had been collected, the volume (ml) was measured using a graduated tube. The ejaculates with abnormal colours (in case of the presence of blood or urine, etc.) or with gel were not evaluated. The samples were placed in an incubator (37°C) for a period of 10 to 15 min to prevent cold shock. pH was measured by a pH meter and each ejaculate was evaluated manually and examined by microscopy. For measuring motility of spermatozoa, semen was diluted with 2.9% sodium citrate dihydrate solution (37°C) and was recorded on a subjective scale of 0% to 100% after viewing several microscopic fields. Percentages of live, dead and abnormal sperm were recorded as described by Chemineau et al. (1991). Duplicate smears from each ejaculate were stained with eosin–igrosin stain. A total of 200 spermatozoa were examined randomly (100 in each of the two smears). Sperm cell concentration (×10⁶/ml) was quantified by direct cell count using the improved Neubauer haemocytometer (Boussit, 1989).

Table 2 Distribution of bucks used and semen ejaculates obtained in different genetic groups

V = V-line, S = Saudi Gabali.

Each ejaculate was examined for volume of ejaculate (VE), pH of semen, SC, percentage of motile spermatozoa (MS), percentage of spermatozoa with abnormal form (AS) and percentage of dead spermatozoa (DS). At the time of semen collection for each ejaculate per buck, libido score (LS) of bucks was measured subjectively by giving the buck a score from 0 to 5. The score of 1 was given to the bucks that showed low libido, and the score of 5 to those that exhibited strong libido (i.e. a sexually very active buck).

Model of analysis and estimation of crossbreeding genetic parameters

The animal model used in analysing the semen parameters was

where y was the vector of observed semen parameter for bucks; b was the vector of fixed effects of the genetic group of buck (11 levels), year-season of semen collection (1-year season every 3 months); u_a was the vector of random additive effect of the bucks and sires and dams of bucks; u_p was the vector of random effects of the permanent non-additive effect of the bucks; X, Z_a and Z_p were the incidence matrices relating records to the fixed effects, additive genetic effects and permanent environment, respectively; and e was the vector of random residual effects.

The variance components of the random effects were estimated by REML using analytical gradients methodology implemented in the VCE software (Neumaier and Groeneveld, 1998). These parameters were used to solve the corresponding mixed model equations. The solutions obtained for estimable functions of the genetic groups (represented as the ĝ vector) and the matrix of their estimated (co) variances (Cⁱⁱ) were used to get the estimates of the crossbreeding genetic parameters of the lines, applying the procedure of generalised least-squares (GLS) using the CBE package of Wolf (1996). These estimable parameters represent the differences between direct genetic effects of the lines (), differences between maternal genetic effects of the lines (), individual heterosis (H ^I), maternal heterosis (H ^M) and recombination losses (R ^I) as stated by Dickerson (1992) and Baselga et al. (2003). The interpretation of ĝ in terms of the genetic parameters of the crosses are given in Table 1, where the combination of parameters that explain the effects of different genetic groups of buck are , , H ^I, H ^M and R ^I (the solution for the buck genetic group number 2 was equalled to 0). Thus, we have five parameters to estimate, included in the vector that we call z, thus

Since z explains ĝ through the following model:

with X being the matrix of incidence relating the solutions of buck genetic group effects with the parameters of crossbreeding (Table 1), the GLS solution for z, which we call , will be unique:

where = the variance–covariance matrix of needed to test the significance of the components of z.

Results and discussion

Overall actual means for semen parameters

Results in Table 3 describe the semen performance of bucks raised in the project. Most traits showed remarkable variability particularly in SC and percentages of MS. The ranges for semen parameters obtained (Table 3) were similar to those obtained by Kuzminsky et al. (1996), Minelli et al. (1999), Alvariño (2000), Arroita et al. (2000), Al-Sobayil and Khalil (2002) and Castellini et al. (2003). García-Tomás et al. (2006b) and Vicente et al. (2000) found marked differences in the overall means only for VE but not for the other semen characteristics in line R raised in Spain.

Table 3 Summary statistics for semen characteristics and libido score

VE = volume of ejaculate, pH = pH of semen, SC = concentration of sperm, MS = motility of sperm, AS = abnormal sperm, DS = dead sperm, LS = libido score.

Proportion of the phenotypic variance due to genetic additive effects (h²) and male permanent environmental and genetic non-additive effects (p²)

We define h ² and p ² as the ratios of the additive variance and environmental permanent variance to the phenotypic variance, all of them estimated by REML analysis (Table 4). These parameters, it must be understood, are extended to the two lines and crosses involved in this experiment. Thus, the h ², for example, has not the common meaning of heritability referring to a unique population in genetic equilibrium. The estimates of h ² and p ² were moderate, ranging from 0.12 to 0.18 and from 0.15 to 0.28, respectively. They are lower than those estimates presented by Brun et al. (2002) in France and by García-Tomás et al. (2006b) in Spain.

Table 4 Estimates of the proportion of the phenotypic variance due to genetic additive effects (h ²) and to permanent environment and genetic non-additive effects (p ²) with their standard errors (±s.e.) for semen characteristics and libido score

VE = volume of ejaculate, pH = pH of semen, SC = concentration of sperm, MS = motility of sperm, AS = abnormal sperm, DS = dead sperm, LS = libido score.

Direct and maternal additive effects

The estimates of direct and maternal additive effects for semen parameters were mostly moderate and in favour of Saudi Gabali breed as shown in Table 5.

Table 5 Estimates of differences between V-line and Saudi Gabali breed in direct and maternal additive effects and their standard errors (±s.e.) for semen characteristics and libido score

VE = volume of ejaculate, pH = pH of semen, SC = concentration of sperm, MS = motility of sperm, AS = abnormal sperm, DS = dead sperm, LS = libido score.

^†Percentage of the difference referred to the average of the values for V-line and Saudi Gabali breed; NS = non-significant; *significant at the 5% level.

Moderate and significant estimates of differences in direct additive effects for SC (12.4%), percentage of DS (14.2%) and percentage of AS (12.0%) indicate that in the bucks, the Saudi genes had higher value than the V-line genes (Table 5). Also, García-Tomás et al. (2006a) stated that differences in direct genetic effects between two sire lines (C and R) were significant and relevant for some semen-production traits (e.g. concentration and total number of spermatozoa per ejaculate) and some semen-quality traits (e.g. percentages of sperm viability, percentage of spermatozoa with normal apical ridge, percentage of sperm morphological abnormalities of neck–mid-piece and percentage of sperm with proximal cytoplasmatic droplet) and those differences were of high magnitude (about 50% of the actual mean) and in favour of line C for SC and the total number of spermatozoa per ejaculate.

As shown in Table 5, differences in the maternal additive effect for VE (11.7%) and LS (6.9%) were in favour of genes of the V-line carried by the dams of bucks, while these estimates were in favour of Saudi genes for SC (10.0%), percentage of AS (13.2%), percentage of DS (18.6%) and percentage of MS (10.8%). Also, García-Tomás et al. (2006a) with two sire lines (C and R) observed differences between the two lines in maternal genetic effects and these differences were important for VE (14%) and in favour of line C but in favour of line R for SC (40%), total number of spermatozoa per ejaculate (23%), mass and individual motility and percentage of spermatozoa with presence of morphological abnormalities of neck-piece. Brun et al. (2002) in France stated that differences between reciprocal crosses for individual motility (trait related to energy metabolism) were of considerable importance; they explained this concept on the basis of sex-linked effects as a result of maternal transmission of mitochondria (cell organites involved in energy metabolism).

Direct and maternal effects for sperms concentration and sperms motility (Table 5) were negative and in favour of Saudi genes. Brun et al. (2002) explained the maternal effects on this trait and they suggested that the differences in maternal genetics effects for the traits related to energy metabolism (sperm motility) could be explained as sex-linked effects on account of maternal transmission of the mitochondria. On the other hand, there is an antagonist relationship between motility and pH explained by the metabolic activity of the spermatozoa, which release lactic acid and consequently reduce pH. In our case, this relationship is not shown since the sign of the differences is in concordance with this result but the differences are not significant. Concentration and motility have been found to be positively correlated by Battaglini et al. (1992), Bencheikh (1995), Brun et al. (2002) and recently by García-Tomás et al. (2006b). This could explain the concordance in the sign of the differences between direct and maternal effects for both traits. It should also be indicated that maternal genetic effects were always of the same sign than direct genetic effects (Table 5), which agree with Brun et al. (2002) for volume and percent of MS and do not agree with García-Tomás et al. (2006b) for concentration and presence of spermatozoa with morphological abnormalities.

Direct and maternal heterosis

Positive estimates of direct heterosis recorded for VE (10.6%), SC (13.5%), motility of sperm (10.5%) and LS (6.2%) were favourable (Table 6). One of the explanations for positive heterotic effects in the percentage of sperm motility could be that sexual maturation in crossbred males was faster than in purebred males. In other species, such as pigs, it seems that crossbred males could be more advanced in sexual maturity. However, at the adult stage, differences among purebred and crossbred bucks in semen parameters tend to disappear (Wilson et al., 1977; Buchanan and Johnson, 1984). In general, crossbreeding among breeds raised in hot climates was associated with heterotic effects in semen characteristics as stated by El-Ezz et al. (1985) in Egypt and Al-Sobayil and Khalil (2002) in Saudi Arabia.

Table 6 Estimates of direct and maternal heterosis and their standard errors (±s.e.) for semen characteristics and libido score

VE = volume of ejaculate, pH = pH of semen, SC = concentration of sperm, MS = motility of sperm, AS = abnormal sperm, DS = dead sperm, LS = libido score.

^†Percentage of the heterosis referred to the average of the values for V-line and Saudi Gabali breed; NS = non-significant; *significant at the 5% level.

The estimates of direct heterosis for percentages of abnormal sperm (−21.5%; P < 0.05) and dead sperm (−20.3%; P < 0.05) were negative and moderate. Brun et al. (2002) reported high variability in the estimates of heterosis in function of the seminal trait since they observed a 6.8% heterosis in mass motility and a 4.1% in percentage of MS along with high values of heterosis in SC (37.5%), total number of spermatozoa per ejaculate (37.6%) and number of MS per ejaculate (42.3%). García-Tomás et al. (2006a) found high variabilities in the estimates of direct heterosis for several semen characteristics, being practically negligible for sperm normalcy (about 2%) but very high for the percentage of spermatozoa with of cytoplasmatic droplet (57%).

Estimates of maternal heterosis for different traits were significant and moderate for six traits out of seven (Table 6). The significant and positive estimates of maternal heterosis for VE (24.0%), SC (10.3%), percentages of motile sperm (21.8%) and libido of bucks (10.2%) and the negative and significant estimates for percentages of abnormal sperm (−9.6%) and dead sperm (−14.7%) show the interest of using crossbred dams to produce the crossbred bucks. It is also interesting to note that the values of the individual and maternal heterosis are very positive to determine the semen characters of the newly founded synthetic lines (Saudi 2 and 3 lines, Table 1).

Recombination losses

The recombination losses were always non-significant for all traits studied (Table 7) except that the VE that was significant but with low magnitude. These results indicate that recombination losses for these traits were of limited importance in crossbred bucks when compared with direct heterosis. In this concept, we found that the estimates of direct heterosis for the majority of the studied traits (Table 6) were generally the largest.

Table 7 Estimates of recombination losses (R ^I) and their standard errors (±s.e.) for semen characteristics and libido score

VE = volume of ejaculate, pH = pH of semen, SC = concentration of sperm, MS = motility of sperm, AS = abnormal sperm, DS = dead sperm, LS = libido score.

^†NS = non-significant; *significant at the 5% level.

Conclusions

Significant estimates of direct and maternal additive effects for most semen characters showed the superiority of the Saudi Gabali breed compared with the V-line. Moreover, the majority of semen parameters were improved in crossbred bucks relative to purebred bucks of V-line or Saudi Gabali.

The signs and values of individual and maternal heterosis were favourable to the expected semen characters of crossbred bucks, and of bucks of synthetic lines derived from crosses between V-line and Saudi Gabali.