Search

This technical note presents the program PaGELL v.1.5 (Parametric Genetic Evaluation of Lifespan in Livestock), a flexible software program to analyze (right-censored) longevity data in livestock populations, with a special emphasis on the genetic evaluation of the breeding stock. This software relies on a parametric generalization of the proportional hazard model; more specifically, the baseline hazard function follows a Weibull process and flexibility is gained by including an additional time-dependent effect with the number of change points defined by the user. The program can accommodate 3 different sources of variation (i.e., systematic, permanent environmental, and additive genetic effects) and both fixed and time-dependent patterns (only for systematic and permanent environmental effects). Analyses are performed within a Bayesian context by sampling from the joint posterior distribution of the model, and model fit can be easily determined by the calculation of the deviance information criterion. Although this software has already been used on field data sets, its performance has been double-checked on simulated data set, and results are presented in this technical note. PaGELL v.1.5 was written in Fortran 95 language and, after compiling with the GNU Fortran Compiler v.4.7 and later, it has been tested in Windows, Linux, and MacOS operating systems (both 32- and 64-bit platforms). [ABSTRACT FROM AUTHOR]

Although heteroskedasticity has been a main topic of interest in beef cattle during recent decades, the current availability of canalization models provided new insights for animal breeding programs. Within this context, birth BW (BWT) was analyzed in the Bruna dels Pirineus beef cattle breed by implementing canalization models that accounted for heterogeneous residual variances due to systematic, permanent environmental effects and additive genetic effects. Analyses were performed on BWT data from 8,130 calves born in 12 commercial breeding herds contributing to the yield recording scheme of the Bruna dels Pirineus breed. Analytical models accounted for direct additive genetic, permanent environmental, and 4 systematic effects (i.e., age of the dam, sex of the calf, birth type, and herd-year-season), and the same effects were evaluated as potential sources of variation in the residual term. Their relevance was checked by the deviance information criterion (DIC), and only residual additive genetic, permanent environmental, birth type, and herd-year-season remained in the operational model, all of them originating relevant reductions in the DIC parameter. Bruna dels Pirineus calves showed a moderate heritability of 0.30 (95% high posterior density, 0.19 to 0.40) for BWT; additional additive genetic variability was revealed in the residual term, this being positively correlated with the direct additive genetic component (0.44; 95% high posterior density, 0.37 to 0.54). Genetic trends were evaluated on both sources of additive genetic variance, and relevant patterns were identified in several herds. Although this breed did not evidence a homogeneous genetic trend for the whole population, herd-specific positive and negative trends were revealed, suggesting the plausibility of genetic selection for canalization on BWT in beef cattle breeds. These results must be viewed as a contribution to the canalization research field, providing relevant information for the breeding scheme of the Bruna dels Pirineus breed, as well as important insights about the genetic background of BWT for the beef industry worldwide. [ABSTRACT FROM AUTHOR]

Lambing interval is a relevant reproductive indicator for sheep populations under continuous mating systems, although there is a shortage of selection programs accounting for this trait in the sheep industry. Both the historical assumption of small genetic background and its unorthodox distribution pattern have limited its implementation as a breeding objective. In this manuscript, statistical performances of 3 alternative parametrizations [i.e., symmetric Gaussian mixed linear (GML) model, skew-Gaussian mixed linear (SGML) model, and piecewise Weibull proportional hazard (PWPH) model] have been compared to elucidate the preferred methodology to handle lambing interval data. More specifically, flock-by-flock analyses were performed on 31,986 lambing interval records (257.3 ± 0.2 d) from 6 purebred RipoUesa flocks. Model performances were compared in terms of deviance information criterion (DIC) and Bayes factor (BF). For all flocks, PWPH models were clearly preferred; they generated a reduction of 1,900 or more DIC units and provided BF estimates larger than 100 (i.e., PWPH models against linear models). These differences were reduced when comparing PWPH models with different number of change points for the baseline hazard function. In 4 flocks, only 2 change points were required to minimize the DIC, whereas 4 and 6 change points were needed for the 2 remaining flocks. These differences demonstrated a remarkable degree of heterogeneity across sheep flocks that must be properly accounted for in genetic evaluation models to avoid statistical biases and suboptimal genetic trends. Within this context, all 6 RipoUesa flocks revealed substantial genetic background for lambing interval with heritabilities ranging between 0.13 and 0.19. This study provides the first evidence of the suitability of PWPH models for lambing interval analysis, clearly discarding previous parametrizations focused on mixed linear models. [ABSTRACT FROM AUTHOR]

Little is known about mutational variability in livestock, among which only a few mutations with relatively large effects have been reported. In this manuscript, mutational variability was analyzed in 1,765 litter size records from 404 Ripollesa ewes to characterize the magnitude of this genetic source of variation and check the suitability of including mutational effects in genetic evaluations of this breed. Threshold animal models accounting for additive genetic mutations were preferred to models without mutational contributions, with an average difference in the deviance information criterion of more than 5 units. Moreover, the statistical relevance of the additive genetic mutation term was checked through a Bayes factor approach, which showed that the models with mutational variability were 8.5 to 22.7 times more probable than the others. The mutational heritability (percentage of the phenotypic variance accounted for by mutational variance) was 0.6 or 0.9%, depending on whether genetic dominance effects were accounted for by the analytical model. The inclusion of mutational effects in the genetic model for evaluating litter size in Ripollesa ewes called for some minor modifications in the genetic merit order of the individuals evaluated, which suggested that the continuous uploading of new additive mutations could be taken into account to optimize the selection scheme. This study is the first attempt to estimate mutational variances in a livestock species and thereby contribute to better characterization of the genetic background of productive traits of interest. [ABSTRACT FROM AUTHOR]

Although inbreeding (F) is a topic of major concern in animal breeding, estimates of inbreeding depression are usually obtained by modeling the overall F coefficient of each individual, without considering that the recessive (deleterious) genetic load of a given population may be unevenly distributed among the founder genomes. The founder-specific partial F coefficient is calculated as the identity-by-descent probability at any given autosomal locus related to a particular founder and allows a more detailed analysis of inbreeding depression on productive traits. Within this context, birth BW data from 2,459 Ripollesa lambs were analyzed under a hierarchical animal model without F-related covariates (model 0), with inbreeding depression modeled by the overall F coefficient (model F1), or by the partial F coefficient of 9 founders that made a relevant contribution to the population inbreeding (model F2). A straightforward empirical Bayes factor (BF) was developed for testing statistical relevance of each F-related covariate, in which greater-than-1 values favored the model including the covariate. The deviance information criterion (DIC) clearly supported model F1 (5,767.8) rather than model 0 (5,771.2), suggesting that inbreeding depression had a relevant influence on birth BW data. The linear effect of inbreeding depression was statistically relevant in model F1 (BF = 2.52 × 1035), with lamb birth BW declining by -13.6 g with each 1% F increase. The quadratic effect of inbreeding depression was almost null in model F1 (BF = 0.02), as suggested by the reduction in DIC (5,766.9) when this effect was removed from model F1. On the other hand, model F2 provided a similar DIC (5,767.9) value, with this parameter decreasing to 5,764.7 when nonrelevant founder-specific inbreeding depression effects were removed. Substantial heterogeneity in founder-specific inbreeding depression was reported by model F2, in which estimates for 4 of the 9 founders did not differ from zero (BF between 0.05 and 0.42), whereas 5 founders originated moderate (-8.2 g for each 1% F increase; BF = 1.42) to large inbreeding depression (-96.2 g for each 1% F increase; BF = 8.80 × 1019). The substantial variability between founder estimates suggested that inbreeding depression effects may mainly be due to a few alleles with major deleterious effects. These results contribute valuable information that should help to achieve more accurate management of inbreeding in the Ripollesa breed. [ABSTRACT FROM AUTHOR]

A performance data set from 376 Ripollesa purebred ewes of the experimental flock of the Universitat Autônoma of Barcelona was analyzed using a bivariate Bayesian threshold-linear model. The data set contained 1,598 litter size records and 1,699 days-to-lambing records. The model included the additive genetic effect of each animal and 3 nongenetic sources of variation: ewe age, year of lambing, and the permanent environmental effect characterized by the ewe. The flock was phenotypically selected for litter size since 1986, and replacement ewes and rams were selected from the progeny of the more prolific ewes, which had at least 3 deliveries recorded. The phenotypic trend for litter size was positive, whereas days to lambing followed an unclear pattern. Both traits had low heritabilities; 0.13 for litter size and 0.11 for days to lambing. Response to selection was evaluated through (a) the average breeding value of the ewe lambs chosen annually, and (b) the average breeding value of the overall flock. The first measurement suggested a positive trend for litter size, although it showed important oscillations. On the other hand, the average breeding value for the overall flock showed a stable positive tendency after yr 4 of selection, with estimates clearly different from zero after yr 11 of selection. A significant increase in the incidence of multiple births was observed, with a mode of approximately 10%. The correlated response in days to lambing did not show a significant trend. The effect of year of lambing also positively influenced both litter size and days to lambing, although important oscillations were observed between years. Results indicated that litter size in sheep can be effectively improved through phenotypic selection, even in small flocks; moreover, days to lambing could also be genetically improved, given the estimate obtained for its heritability. [ABSTRACT FROM AUTHOR]

The analysis focused on model fitting of 2 ewe reproductive traits, litter size, and days to lambing (interval between the introduction of the ram into the flock and the subsequent parturition of the ewes). The experimental data set of the Universitat Autònoma of Barcelona flock was used, including 1,598 records of litter size and 1,699 records of days to lambing from 376 Ripollesa ewes between 1986 and 2005. Univariate and bivariate models were considered as beginning points with linear or threshold approximation for litter size. Model fitting was evaluated in terms of goodness-of-fit and predictive ability, using the mean square error and the correlation between phenotypic and predicted records (ρyŷ) as reference parameters. The bivariate model was preferable for both variables, minimizing mean square error and maximizing ρyŷ. A threshold approximation for litter size was preferable over a linear approximation. Models were also compared with a simulation study, comparing the correlation coefficient between simulated and predicted breeding values (ρa,â). The bivariate threshold model was favored, with a ρyŷ of 0.677 and 0.834 for litter size and days to lambing, respectively. Correlation coefficients between simulated and predicted breeding values in the bivariate linear model were reduced slightly to 0.651 and 0.831, respectively, and they were lowest with linear univariate models (0.642 and 0.802). Although the bivariate models for ewe litter size and days to lambing were more accurate than the univariate models, the threshold approaches showed a greater advantage under the bivariate model. For the purpose of genetic evaluation of litter size in sheep, use of the threshold-linear model seems justified. In the Ripollesa breed, the evaluation of litter size can benefit from recording birth weight. [ABSTRACT FROM AUTHOR]

The Alberes population is a native bovine breed of Catalonia with an unclear origin, which historically some authors have assumed as being composed of two different colour varieties (black and fawn). Sixteen microsatellite loci were analysed, all of them included in the AIRE2066 European Concerted Action list. Overall expected and observed heterozygosities reached values of 0.649 and 0.662, respectively. Genetic differences among black and fawn varieties were not significant (FST = 0.007), and therefore the population is a single variety with a great colour gradation. On the contrary, we detected significant genetic differences among herds (FST = 0.026; p < 0.001), showing a genetic heterogeneity over short geographical distances. The number of migrants per generation among pairs of herds oscillates between 1.46 ( Roig and Freixe herds) and 5.62 ( Castanyers and Roig herds). Moreover, inbreeding and bottleneck situations can be rejected. The Alberes breed has been grouped within the Cantabrian trunk, closely related to the Asturiana de la Montaña and Alistana breeds, although some other breeds may also have influenced the population along its history. [ABSTRACT FROM AUTHOR]