Genome-wide association study for feed efficiency in collective cage-raised rabbits under full and restricted feeding.

Feed efficiency (FE) is one of the most economically and environmentally relevant traits in the animal production sector. The objective of this study was to gain knowledge about the genetic control of FE in rabbits. To this end, GWASs were conducted for individual growth under two feeding regimes (full feeding and restricted) and FE traits collected from cage groups, using 114 604 autosome SNPs segregating in 438 rabbits. Two different models were implemented: (1) an animal model with a linear regression on each SNP allele for growth trait; and (2) a two-trait animal model, jointly fitting the performance trait and each SNP allele content, for FE traits. This last modeling strategy is a new tool applied to GWAS and allows information to be considered from non-genotyped individuals whose contribution is relevant in the group average traits. A total of 189 SNPs in 17 chromosomal regions were declared to be significantly associated with any of the five analyzed traits at a chromosome-wide level. In 12 of these regions, 20 candidate genes were proposed to explain the variation of the analyzed traits, including genes such as FTO, NDUFAF6 and CEBPA previously associated with growth and FE traits in monogastric species. Candidate genes associated with behavioral patterns were also identified. Overall, our results can be considered as the foundation for future functional research to unravel the actual causal mutations regulating growth and FE in rabbits.

Variability, repeatability and phenotypic relationships of several characteristics of production and semen quality in rabbit.

A total of 2140 ejaculates from 156 adult males pertaining to four groups of bucks were analysed and the following traits were recorded: pH, ejaculate volume (V), mass motility (Mm), individual motility (Mi), concentration (Cn), total number of spermatozoa per ejaculate (TSE), percentage of sperm viability (Vi), percentage of sperm with acrosome integrity (NAR), percentage of sperm normalcy (Nr), percentage of sperm morphological abnormalities of head (H), neck-midpiece (Nm) and tail (T) and presence of proximal and distal cytoplasmic droplet (Dp, Dd). Principal component (PC) analysis and phenotypic correlations were performed in order to examine the relationships between qualitative and quantitative traits of rabbit semen. The repeatability of the traits measured was also estimated. Phenotypic correlations between sperm traits were estimated as the residual correlation from an analysis of variance, including the effects of: genetic type of the male, order of the ejaculate, day of collection and the permanent non-additive random effect of the male to which the observation corresponds. The repeatability of these traits was analysed separately, in a set of univariate analyses, using VCE software. The previously defined mixed model was used for this analysis. The principal component analysis was performed using the Princomp procedure of the SAS v.8 package. The first four PC explained 62% of total variation: 23%, 18%, 12% and 8%, respectively. Percentage of sperm viability, NAR, Nr, T and Nm were the predominant variables in the first PC. Mass motility, Mi, pH, Cn and TSE were located in the second. Percentage of sperm with the presence of proximal and distal cytoplasmic droplet, NAR, Vi and Nr were the predominant traits in the third and V defined the fourth. Repeatability of semen quality traits was moderate for most of the traits. The highest values (about 0.45) corresponded to V, Cn, Nr and Dp and the lowest (about 0.10) to H and Nm. Repeatability of the first ejaculate was greater than that of the second for pH, Mm, Mi, H, Dp and Dd.

Lactating performance, water and feed consumption of rabbit does reared under a Mediterranean summer circadian cycle of temperature v. comfort temperature conditions.

The general aim of this research was to study the effect of high ambient temperature on the performance of does during lactation, specifically the following factors: average daily feed (ADFI) and water (ADWI) intakes, daily milk yield (DMY); milk composition: dry matter (DM), CP and gross energy (GE); doe BW (DW); individual kit weaning weight (IWW) and litter survival rate during lactation (SR). The study was undertaken comparing the performance of two groups of contemporary does reared under the same management, feeding regime and environmental conditions, except the environmental temperature and humidity. A total of 80 females were randomly allocated, at 60 days of age, into two identical and continuous rooms. In one room, the temperature was maintained permanently within the thermo-neutral zone (between 18°C to 22°C); thus, environmental conditions in this room were considered as comfort conditions. In the second room, the environmental temperature pattern simulated the daily temperature cycles that were characteristic of the summer in Mediterranean countries (24°C at 0800 h, increasing up to 29°C until 1100 h; maintenance at 29°C to 31°C for 4 h and decreasing to about 24°C to 26°C around 1700 h until 0800 h of the following day), which were considered as thermal stress conditions. Females followed a semi-intensive reproductive rhythm, first artificial insemination at 4.5 months of age, with subsequent 42-day reproductive cycles. Traits were recorded from a total of 138 lactations. Does were controlled up to the 5th lactation. Data were analyzed using linear and linear mixed models. High ambient temperature led to a lower ADFI (-9.4%), DW (-6.2%) and IWW (-8%), but it did not affect ADWI. No significant difference was found either for DMY, milk composition (DM, CP and GE) and SR during the lactation period. Heat stress was moderate, and does were able to adapt to it behaviorally by decreasing feed intake (to reduce heat production), but also live weight, allowing them to preserve milk yield and composition for assuring litter survival. On the other hand, water consumption could not be the main animal mechanism to overcome heat stress.

Daily exposure to summer temperatures affects the motile subpopulation structure of epididymal sperm cells but not male fertility in an in vivo rabbit model.

High temperatures have negative effects on sperm quality leading to temporary or permanent sterility. The aim of the study was to assess the effect of long exposure to summer circadian heat stress cycles on sperm parameters and the motile subpopulation structure of epididymal sperm cells from rabbit bucks. Twelve White New Zealand rabbit bucks were exposed to a daily constant temperature of the thermoneutral zone (from 18 °C to 22 °C; control group) or exposed to a summer circadian heat stress cycles (30 °C, 3 h/day; heat stress group). Spermatozoa were flushed from the epididymis and assessed for sperm quality parameters at recovery. Sperm total motility and progressivity were negatively affected by high temperatures (P < 0.05), as were also specific motility parameters (curvilinear velocity, linear velocity, mean velocity, straightness coefficient, linearity coefficient, wobble coefficient, and frequency of head displacement; P < 0.05, but not the mean amplitude of lateral head displacement). Heat stress significantly increased the percentage of less-motile sperm subpopulations, although the percentage of the high-motile subpopulation was maintained, which is consistent with the fact that no effect was detected on fertility rates. However, prolificacy was reduced in females submitted to heat stress when inseminated by control bucks. In conclusion, our results suggest that environmental high temperatures are linked to changes in the proportion of motile sperm subpopulations of the epididymis, although fertility is still preserved despite the detrimental effects of heat stress. On the other hand, prolificacy seems to be affected by the negative effects of high temperatures, especially by altering female reproduction.

Comparison of behaviour, performance and mortality in restricted and ad libitum-fed growing rabbits.

The objective of this study was to determine whether rabbits fed in a restricted regimen (75%) showed increased competition for feeding, drinking and use of specific areas of the cages as compared with those provided feed ad libitum. This evaluation was carried out by measuring their space utilisation in the cage, the incidence of agonistic behaviour and rates of mortality. In total, 504 rabbits between 31 and 66 days of age were used in this study. A total of 200 heavy-weight rabbits and 56 light-weight rabbits were randomly housed in 32 cages, each cage containing eight rabbits: 25 cages housing heavy rabbits and seven cages housing the light-weight ones. They were all fed ad libitum (AD). In addition, a total of 208 heavy-weight rabbits and 40 light-weight rabbits were randomly housed in 31 cages, each of them containing eight rabbits: 26 cages housing heavy weight rabbits and five cages housing light-weight ones. They were all fed a restricted diet (R) regimen. The restriction was calculated to be 75% of the feed consumed by the AD group. The total space available in the cage was 3252 cm(2), with a stocking density of 24.6 animals/m(2). Animals between 32 and 60 days of age from 20 different cages were observed nine times per week (morning or afternoon) by means of scan and focal sampling by one observer. During each period, cages were assessed for 5 min, registering every minute the position of all the animals in relation to Area A (feeder), Area B (central part) or Area C (back and drinker area). The incidence of agonistic behaviour such as displacement, biting and jumping on each other was also assessed. Performance variables such as daily gain and feed conversion ratio, in addition to general health status and mortality rates, were recorded for all rabbits. When the rabbits were under restricted feeding, the competition for feed and drink increased with clear signs of agonistic behaviour such as biting, displacement and animals jumping on top of each other. Although this competition was maintained during the entire growing period, the BW homogeneity between animals in the same cage was similar in both cases, suggesting that all animals could consume similar quantities of feed. The possible advantages of a restricted diet, such as better feed conversion ratio, were observed in this study only in the last few weeks of the growing period.

Elliptical selection experiment for the estimation of genetic parameters of the growth rate and feed conversion ratio in rabbits.

Two elliptical selection experiments were performed in two contemporary sire lines of rabbits (C and R) in order to optimize the experimental design for estimating the genetic parameters of the growth rate (GR) and feed conversion ratio (FCR). Twelve males and 19 females from line C, and 13 males and 23 females from line R, were selected from an ellipse defined by a quadratic index based on these traits. Data from 160 rabbits of each of the parental generations of lines C and R and their offspring (275 and 266 animals, respectively) were used for the analysis. A Bayesian framework was adopted for inference. Marginal posterior distributions of the genetic parameters were obtained by Gibbs sampling. An animal model including batch, parity order, litter size, and common environmental litter effects was assumed. Posterior means (posterior standard deviations) for heritabilities of GR and FCR were estimated to be 0.31 (0.10) and 0.31 (0.10), respectively, in line C and 0.21 (0.08) and 0.25 (0.12) in line R. Posterior means of the proportion of the variance due to common litter environmental effects were 0.14 (0.06) and 0.21 (0.06) for GR and FCR, respectively, in line C and 0.17 (0.06) and 0.22 (0.06) in line R. Posterior means of genetic correlation between both traits were -0.49 (0.25) in line C and -0.47 (0.32) in line R, indicating that selection for GR was expected to result in a similar correlated response in FCR in both lines.

Effect of heat intensity and persistency on prolificacy and preweaning kit growth at different stages of the rabbit production cycle.

The aim of this research was to assess the effect of temperature intensity and variation both throughout the day and between days within different periods of the reproductive cycle and the lactation of the rabbit. This information would help in establishing optimal patterns of environmental temperature control in rabbitries. The traits analyzed were total number of kits born (TB), number of kits alive at weaning (NW), and average individual weight at weaning (AvgWW). For each trait, several mixed models were fitted to the data, differing only in the number and type of temperature descriptors included in the vector of fixed effects. Those descriptors were the average daily mean, maximum, and range of temperatures (AvgTmean, AvgTmax, and AvgRg, respectively) and the CV of daily mean temperature (CVTmean). All were calculated for periods in which important physiological processes related to the studied traits occur. High environmental temperature was found to have a detrimental effect on prolificacy and preweaning growth of the kits. When the average daily mean reached 20°C, it produced a linear decay of TB of around 0.1 kit/°C. The most sensitive period for TB could cover from spermatogenesis to embryo implantation. However, the high correlation between descriptors calculated for different periods makes it difficult to assign an effect to each specific period and therefore to the specific physiological process occurring in that period. The effect on NW was smaller and quadratic, with an optimum value between 18°C and 21°C. Weaning weight was the most strongly affected trait similar to NW. It also showed a quadratic response to AvgTmean, with an optimum value in the same temperature interval as NW and a strong decline in weaning weight with temperatures higher than 21°C (-14 g/°C). There were no differences on the effect of heat at different stages of lactation on NW and AvgWW. The impact of high environmental temperatures on prolificacy is alleviated if a drop in temperature is produced during the day. Thus, the effect of AvgRg was relevant and positive for TB and NW. However, it had a quadratic negative effect for AvgWW at late lactation. The temperature variation between days within a period has a positive effect on TB and AvgWW but a negative effect on NW when it is produced at middle lactation, whereas there is no effect at other stages of lactation.

Use of infrared thermography to assess the influence of high environmental temperature on rabbits.

The aim of this work was to ascertain if infrared thermography (IRT) can be used on rabbits to assess differences in surface body temperature when they are subjected to two different environmental temperatures outside the comfort zone. Rabbits housed in room A were maintained at a temperature of below 30°C and rabbits in room B at a temperature of above 32°C for a year. Faeces were collected six times during the year to assess stress by means of faecal cortisol metabolites (FCM). The assessment of IRT was carried out to assess maximum and minimum temperatures on the eyes, nose and ears. FCM concentration was higher in room B than A, to confirm that stress conditions were higher in room B. Significant differences in IRT were found between the animals housed in both rooms. It was observed that it was more difficult for animals from room B to maintain a regular heat loss. Although all the body zones used to assess temperature with IRT gave statistical differences, the correlations found between the eyes, nose and ears were moderate, suggesting that they were giving different information. In addition, differences up to 3.36°C were found in the eye temperature of rabbits housed in the same room, with a clear effect of their position in relation to extractors and heating equipments. Therefore, IRT could be a good tool to assess heat stress in animals housed on typical rabbit farm buildings, giving a measure of how the animal is perceiving a combination of humidity, temperature and ventilation. Some face areas were better for analysing images. Minimum temperature on eyes and temperatures on nose are suggested to assess heat losses and critical areas of the farm for heat stress in rabbits.

Genetic basis of semen traits and their relationship with growth rate in rabbits.

This work aims to estimate the genetic parameters of seminal and production traits in a paternal line of rabbits selected for ADG during the fattening period. The considered traits were male libido (Lib) defined as successful mounting of an artificial vagina; presence of urine (Ur) and calcium carbonate deposits (Ca) in the ejaculate; semen pH; individual sperm motility (IM); the suitability for AI of the ejaculate (Sui), which involves the subjective combination of several quality traits; the average ejaculate volume (Vol); sperm concentration (Conc); and the average sperm production per ejaculate (Prod = Vol × Conc). The genetic relationship between all of these traits with ADG is also provided. Male libido and seminal data came either from routine evaluations of the ejaculates in an AI center or from 2 experiments in which bucks from the same population were used. Two consecutive ejaculates per male and per week were collected, leaving 7 d within weekly collections. A linear tri-trait model was used to analyze Conc, Vol, and ADG, whereas linear and threshold-linear 2-trait models were used to analyze male libido and the remaining seminal traits with ADG. A Bayesian approach was adopted for inference. Approximately 38% of ejaculates were rejected for AI primarily due to low IM scores. Variables related to the quality of the ejaculate (Ur, Ca, pH, IM, Sui) and Lib were found to be lowly heritable (h(2) ranged from 0.04 to 0.11), but repeatable. This indicates performance of bucks for seminal quality traits and libido in AI centers would be more strongly affected by management practices rather than genetic selection. Semen production traits exhibited moderate values of h(2) (0.22, 0.27, and 0.23 for Conc, Vol, and Prod, respectively), suggesting the possibility of effective selection for these traits. A moderate to high negative genetic correlation (r(g); posterior mean; highest posterior density at 95%, HPD(95%)) was estimated between Conc and Vol (-0.53, HPD(95%) = -0.76, -0.27). The ADG was estimated to have an h(2) of 0.16, to have a low, positive r(g) with Conc (0.21, HPD(95%) = -0.03, 0.48), to have a low, negative r(g) with Vol (-0.19, HPD(95%) = -0.47, 0.08), and to be genetically uncorrelated with all remaining traits analyzed. Therefore, selection for increasing ADG in paternal lines is expected to have no detrimental effects on Ur, Ca, pH, IM, Sui, and Lib and little to no effect on Conc, Vol, and Prod.

Different ways to model biological relationships between fertility and pH of the semen in rabbits.

This work aimed to study the relationship between pH of the semen and fertility (Fert, defined as the success or failure of conception), which is of special interest because pH of the semen can be considered a global marker of the expression of some seminal quality traits. Different methods used to model the relationship between Fert and pH are presented here: 1) ignoring genetic and environmental correlations and including pH either as a covariate or as a cross-classified effect on fertility, 2) a bivariate mixed model, and 3) recursive bivariate mixed models. A total of 653 pH records and 6,365 Fert records after AI were used. Crossbreed does from 2 maternal lines were artificially inseminated with buck semen from a paternal line in a commercial environment. A negative, and almost linear, effect of pH on Fert was detected. The posterior median of pH and Fert heritabilities, and the highest posterior density interval at 95% (in parentheses) were approximately 0.18 (0.05, 0.29) and approximately 0.10 (0.02, 0.20) across all the models, respectively. Genetic correlations between traits were negative, but the highest posterior density interval at 95% included zero [i.e., -0.31 (-0.91, 0.33) in the bivariate mixed model and -0.17 (-0.99, 0.48) and -0.44 (-0.99, 0.10) in the recursive bivariate mixed models including pH as a covariate or as a cross-classified effect, respectively]. All models predicted Fert data reasonably well (i.e., 76 and 62% correct predictions for success and failure, respectively). No differences in the prediction of the EBV for male fertility were encountered between models, showing a good concordance in the animals ranked by their EBV (the correlation between EBV in all models was close to 1). Thus, no differences in results were obtained considering, or not considering, genetic and environmental correlations between pH and Fert and assuming, or not assuming, recursiveness between each trait. This is because the magnitude of the effect of pH on Fert was not large enough; therefore, the same results were obtained even though the models were of different complexity.

Using the product threshold model for estimating separately the effect of temperature on male and female fertility.

Animals under environmental thermal stress conditions have reduced fertility due to impairment of some mechanisms involved in their reproductive performance that are different in males and females. As a consequence, the most sensitive periods of time and the magnitude of effect of temperature on fertility can differ between sexes. The objective of this study was to estimate separately the effect of temperature in different periods around the insemination time on male and on female fertility by using the product threshold model. This model assumes that an observed reproduction outcome is the result of the product of 2 unobserved variables corresponding to the unobserved fertilities of the 2 individuals involved in the mating. A total of 7,625 AI records from rabbits belonging to a line selected for growth rate and indoor daily temperature records were used. The average maximum daily temperature and the proportion of days in which the maximum temperature was greater than 25°C were used as temperature descriptors. These descriptors were calculated for several periods around the day of AI. In the case of males, 4 periods of time covered different stages of the spermatogenesis, the transit through the epididymus of the sperm, and the day of AI. For females, 5 periods of time covered the phases of preovulatory follicular maturation including day of AI and ovulation, fertilization and peri-implantational stage of the embryos, embryonic and early fetal periods of gestation, and finally, late gestation until birth. The effect of the different temperature descriptors was estimated in the corresponding male and female liabilities in a set of threshold product models. The temperature of the day of AI seems to be the most relevant temperature descriptor affecting male fertility because greater temperature records on the day of AI caused a decrease in male fertility (-6% in male fertility rate with respect to thermoneutrality). Departures from the thermal zone in temperature descriptors covering several periods before AI until early gestation had a negative effect on female fertility, with the pre- and peri-implantational period of the embryos being especially sensitive (from -5 to -6% in female fertility rate with respect to thermoneutrality). The latest period of gestation was unaffected by the temperature. Overall, magnitude and persistency of the temperatures reached in the conditions of this study do not seem to be great enough to have a large effect on male and female rabbit fertility.

Interaction of genotype × artificial insemination conditions for male effect on fertility and prolificacy.

Failures in fertilization or embryogenesis have been shown to be partly the result of poor semen quality. When AI is practiced, fertilization rate depends on the number and quality of spermatozoa in the insemination dose around the time of application. Individual variation in the male effect on fertility (success or failure to conceive; Fert) and prolificacy (total number of kids born per litter; TB) could also depend on these factors, and it could be better observed under limited conditions of AI, such as decreased sperm concentration, small or null preselection of ejaculates for any semen quality trait, or a long storage period of the AI doses. The aim of this research was to determine if an interaction existed between male genotype and the AI conditions for male effects on Fert and TB after AI was performed under different conditions. Fertility and TB were assumed to be different traits and were analyzed in 2 sets of independent analyses. In the first step, the different conditions were determined uniquely by the sperm dosage. Artificial insemination was performed at 10 and 40 × 10(6) spermatozoa/mL. In the second step, the different conditions were determined by all the factors involved in the AI process as a whole (conditions and duration of the storage period of the dose, genetic type of the female, and environmental conditions on the farm). Data from AI from the former experiment were analyzed with data from AI performed under different conditions. Threshold and linear 2-trait models were assumed for Fert and TB, respectively. The sperm dosage had a clear effect on Fert and TB, which favored the greater dosage (+0.13% and +1.25 kids born, respectively). Prolificacy was more sensitive to sperm reduction than was fertility. Male heritabilities for Fert were 0.09 for both sperm dosages, and were 0.08 and 0.06 for male TB with a smaller and larger sperm dosage, respectively. No genotype × sperm dosage interaction was found. Therefore, the same response to selection to improve male Fert and TB could be achieved at any sperm concentration. However, an interaction between male genotype and the AI conditions as a whole seemed to exist, indicating that the AI conditions for selection for Fert and TB could be modified to maximize genetic progress. Consequently, the optimization of a breeding program for male Fert and TB under a given set of semen utilization conditions is achievable.

Crossbreeding parameters for growth and feed consumption traits from a five diallel mating scheme in rabbits.

Twenty-five genetic groups of young rabbits originated from a diallel-crossbreeding scheme among 5 selected lines: 3 maternal and 2 terminal sire lines belonging to 2 Spanish research institutions. A sample of 2,773 young rabbits from 525 litters was controlled during the fattening period lasting from 5 to 9 wk. Growth and feed consumption traits were evaluated throughout different biweekly batches. A Bayesian approach was used for inference from an animal model with common litter effects. On average, genetic groups coming from lines selected for growth rate were heavier (+58.9 g at 32 d and +315.5 g at 60 d), had greater growth rate (9.24 and 8.15 g/d from individual or cage analysis, respectively) and feed intake (+13.24 g/d), and showed better feed conversion ratio (-0.21 g of intake/g of gain), than the genetic groups originated from crosses among lines selected for litter size. Crossbreeding parameters were estimated from the samples of the marginal posterior distribution of the genetic group effect according to the Dickerson model. Maternal genetic and individual heterosis effects were null or very low. Direct genetic effects mainly regulated the expression of growth traits.

Individual efficiency for the use of feed resources in rabbits.

A Bayesian procedure, which allows consideration of the individual variation in the feed resource allocation pattern, is described and implemented in 2 sire lines of rabbit (Caldes and R). The procedure is based on a hierarchical Bayesian scheme, where the first stage of the model consists of a multiple regression model of feed intake on metabolic BW and BW gain. In a second stage, an animal model was assumed including batch, parity order, litter size, and common environmental litter effects. Animals were reared during the fattening period (from weaning at 32 d of age to 60 d of age) in individual cages on an experimental farm, and were fed ad libitum with a commercial diet. Body weight (g) and cumulative feed intake (g) were recorded weekly. Individual BW gain (g) and average BW (ABW, g) were calculated from these data for each 7-d period. Metabolic BW (g(0.75)) was estimated as ABW(0.75). The number of animals actually measured was 444 and 445 in the Caldes and R lines, respectively. Marginal posterior distributions of the genetic parameters were obtained by Gibbs sampling. Posterior means (posterior SD) for heritabilities for partial coefficients of regression of feed intake on metabolic BW and feed intake on BW gain were estimated to be 0.35 (0.17) and 0.40 (0.17), respectively, in the Caldes line and 0.26 (0.19) and 0.27 (0.14), respectively, in line R. The estimated posterior means (posterior SD) for the proportion of the phenotypic variance due to common litter environmental effects of the same coefficients of regression were respectively, 0.39 (0.14) and 0.28 (0.13) in the Caldes line and 0.44 (0.22) and 0.49 (0.14) in line R. These results suggest that efficiency of use of feed resources could be improved by including these coefficients in an index of selection.

Crossbreeding parameter estimation for functional longevity in rabbits using survival analysis methodology.

A complete diallel cross involving 3 maternal lines of rabbit was performed to estimate cross-breeding parameters for functional longevity. This trait was defined as the ability to delay involuntary culling. The lines considered, A, V, and Prat, had all been selected by litter size at weaning for a long period. Data were related to a total of 653 does belonging to the 9 genetic types from the diallel cross; does were reared and bred on the same commercial farm. Survival analysis was performed using a Cox proportional hazard model. The model incorporated time-dependent factors, such as year-season, litter size, and the interaction between cycle and physiological status of the female; time-independent factors, such as the genetic type of the doe; and sire and dam random factors. Crossbreeding parameters were estimated from the solutions obtained for the type of doe and its estimated variance-covariance matrix, using a generalized least squares procedure. The estimated parameters were the differences between lines in direct genetic effects and maternal genetic effects and individual heterosis. Relevant differences were observed in direct genetic effects between lines A and Prat but not in any maternal effects. Heterosis was found to be significant and favorable between lines A and Prat, and between the lines V and Prat. The magnitude of this effect was variable but important, especially in the first cross. Results stress the importance of using crosses between specialized lines to produce does for intensive meat rabbit production.

Genetics of litter size in three maternal lines of rabbits: repeatability versus multiple-trait models.

Variance components were estimated in 3 lines of rabbits selected for litter size at weaning (A, Prat, and V) to test one of the assumptions of the models used for selection: that litter size data at different parities are repeated measurements of the same trait. Multiple-trait analyses were performed for the total number of kits born (TB), the number of kits born alive (BA), and the number of kits weaned (NW) per litter. Estimates were obtained by REML in multivariate analyses, including all of the information of the selection criteria, under a repeatability model or a multiple-trait model, considering litter size at the first, second, and subsequent parities as different traits. Models included the fixed effects of the physiological status of the female and the year-season of mating day, buck and doe random permanent environmental effects, and doe additive genetic effects. Results indicated that prolificacy was determined mainly by doe components and that the service sire had a very small effect. Heritabilities for the first and second parities were greater than the estimates obtained under the repeatability model (0.04 to 0.14 for the repeatability model). In the A and V lines, similar values of heritability were found at the first and second parities, but in the Prat line heritability at the second parity was always greater than at the first and greater parities (values of 0.21, 0.17, and 0.15 for TB, BA, and NW, respectively, in second parities of the Prat line). Genetic correlations between the same traits at different parities were approximately 0.8 for all traits in line A, but much lower in the other 2 lines. On average, the values were 0.64 for TB, 0.48 for BA, and 0.39 for NW between the first and second parities, and 0.65 for TB, 0.56 for BA, and 0.45 for NW between the first and third and greater parities. Genetic correlations between the second and greater parities showed the greatest values (approximately 0.8) in lines A and Prat for all traits, but they were lower in line V (0.63 for BA and 0.37 for NW). The heterogeneity of heritabilities and genetic correlations between parities lower than 0.9 suggests that litter size at different parities could be considered as different traits when genetic evaluations are performed. However, when the accuracies of predicted breeding values under a multiple-trait and a repeatability model were calculated, assuming the first to be the true model, the values obtained were nearly the same for all traits in all lines.

Survival analysis in two lines of rabbits selected for reproductive traits.

The objective of the study was to analyze the reproductive longevity of 2 selected lines of rabbits. The first one was the Prat line, a line selected in Spain on litter size at weaning, and raised in overlapping generations. The second one was a French line, the A1077 line, selected on litter size at birth and individual weight at 63 d of age, managed in batches, and artificially inseminated with discrete generations. Reproductive longevity was measured beginning at the first successful mating, assessed by a pregnancy diagnosis in the Prat line, and at the first kindling in the A1077 line. In the A1077 line, culling for infertility occurred after 3 unsuccessful artificial inseminations. The trait analyzed, defined as the doe length of productive life (LPL), was the time in days between date of the first positive pregnancy diagnosis and date of culling or death in the Prat line. In the A1077 line, the trait was the number of AI after the first kindling. Effects included in the model were year-season, litter size at birth, reproductive cycle or physiological status x cycle interaction, age at first mating, batch (only for the A1077 line), and additive genetic value of the animal as a random effect. Survival analyses were carried out with a Cox model for the Prat line and a discrete model for the A1077 line. The estimated heritability values for LPL were around 0.16 in the Prat and A1077 lines with a model including physiological status x cycle interaction effect. Removing this effect from the model led to an increase in estimated genetic variance with h2 = 0.24 and 0.19 in the Prat and A1077 lines, respectively. Including the traits LPL and number of AI from first fertile mating or AI in selection programs could increase reproductive longevity and decrease the replacement rate.

Genetic parameters of fertility in two lines of rabbits with different reproductive potential.

A Bayesian analysis with a threshold model was performed for fertility defined as a binary trait (1 = successful mating, 0 = unsuccessful mating) in two populations of rabbits of different reproductive potential and different genetic origin: Line P selected for litter size and Line C selected for growth rate. There were 20,793 records of natural mating (86.2% successful) in Line C between 1983 and 2003, and 17,548 records (80.5% successful) in Line P, between 1992 and 2003. Data related to 5,388 and 3,848 females and 1,021 and 685 males in Lines C and P, respectively. The pedigree included 6,409 and 4,533 individuals in Lines C and P, respectively. The binary response was modeled under a probit approach. The model for the latent variable included male and female additive genetic effects, male and female permanent environmental effects, and the year-season and physiological status of the female (nulliparous, multiparous lactating, or multiparous nonlactating) as systematic effects. Means (standard deviation in parentheses) of the estimated marginal posterior distribution (EMPD) of male heritability were 0.013 (0.006) and 0.010 (0.008) in Lines C and P, respectively, and those of EMPD of female heritability were 0.056 (0.013) and 0.062 (0.018) in Lines C and P, respectively. Means of the EMPD of the proportion of the phenotypic variance due to environmental male and female effects were, respectively, 0.031 (0.007) and 0.128 (0.018) in Line C and 0.053 (0.010) and 0.231 (0.024) in Line P. Means (standard deviations in parentheses) of the EMPD of genetic correlation between male and female fertility were 0.733 (0.197) in Line C and 0.434 (0.381) in Line P. The posterior distribution of genetic correlations presents a huge dispersion, and the estimates should be taken with caution because of the almost negligible estimate of the male genetic component. Results indicate that little genetic variation exists for female fertility, and practically none for male fertility. It would, therefore, be possible to improve reproductive performance by including female fertility in a breeding program, but response to selection would be very small.