Genetic evaluation for latent variables derived from factor analysis in broilers.

1. The aim of this study was to investigate the associations between several carcass, performance and meat quality traits in broilers through factor analysis and use the latent variables (i.e. factors) as pseudo-phenotypes in genetic evaluations.2. Factors were extracted using the principal components method and varimax rotation algorithm. Genetic parameters were estimated via Bayesian inference under a multiple-trait animal model.3. All factors taken together explained 71% of the original variance of the data. The first factor, denominated as 'weight', was associated with carcass and body weight traits; and the second factor, defined as 'tenderness', represented traits related to water-holding capacity and shear force. The third factor, 'colour', was associated with traits related to meat colour, whereas the fourth, referenced as 'viscera', was related to heart, liver and abdominal fat.4. The four biological factors presented moderate to high heritability (ranging from 0.35 to 0.75), which may confer genetic gains in this population.5. In conclusion, it seems possible to reduce the number of traits in the genetic evaluation of broilers using latent variables derived from factor analysis.

Genotypic and allelic frequencies of gene polymorphisms associated with meat tenderness in Nellore beef cattle.

The objectives of this study were to characterize the allelic and genotypic frequencies of polymorphisms in the µ-calpain and calpastatin genes, and to assess their association with meat tenderness and animal growth in Nellore cattle. We evaluated 605 Nellore animals at 24 months of age, on average, at slaughter. The polymorphisms were determined for the molecular markers CAPN316, CAPN530, CAPN4751, CAPN4753, and UOGACAST1. Analyses of meat tenderness at 7, 14, and 21 days of maturation were performed in samples of longissimus thoracis obtained between the 12th and 13th rib and sheared using a Warner Bratzler Shear Force. Significant effects were observed for meat tenderness at days 7, 14, and 21 of maturation for the marker CAPN4751, at day 21 for the marker CAPN4753, and at days 14 and 21 for the marker UOGCAST1. For genotypic combinations of markers, the results were significant for the combination CAPN4751/UOGCAST1 in the three maturation periods and CAPN4753/UOGCAST1 at days 14 and 21 of maturation.

Multibreed genetic evaluation in bovines using simulated data employing a composite population.

Genetic evaluations in Brazil are performed within each animal breed; however, with the wide range of extant genetic groups in the country and the increased use of genetic crossing as a form of rapid meat production, more elaborate programs that can jointly evaluate animals of different genetic groups are needed. Genetic evaluation of a composite breed is difficult because of the variation in the genetic composition of a given herd, as well as the inclusion of non-additive genetic effects among breeds that can be important for selecting traits in certain breed combinations. Newer models include additive and non-additive effects; however, few studies have investigated these aspects in tropical breeds. The aim of this study was to simulate genetic values to compare different models. Non-inclusion of maternal effects in models leads to overestimation of variance and direct heritability. Estimates of the biological additive effects are influenced by the maternal effects; however, estimates of the non-additive effects are minimally influenced by the maternal effects and are well estimated in any situation. The studied models effectively predict the direct genetic values.

Estimation of genetic parameters and genetic changes for growth characteristics of Santa Ines sheep.

Studying genetic parameters and genetic changes in Santa Ines sheep is important, because it is the commonest breed in Brazil. This study obtained genetic data from 37,735 pedigree records of lambs over 12 years (2003-2014) from 33 flocks in 10 Brazilian States; 11,851 records of performance were available. (Co)variance components, genetic parameters and breeding values estimates were obtained by derivative-free restricted maximum likelihood in a univariate analysis that included maternal additive genetic and maternal permanent environmental effects. Birth weight, weaning weight, weight at 180 days of age, weight at 270 days of age, average daily weight gain in the following states: from birth to weaning, from weaning to 6 months, from 6 months to 9 months, and from weaning to 9 months; presence of hair in fur and leg muscularity were assessed. (Co)variance component values increased in the weight traits with age. A significant maternal effect was found in the pre-weaned stage that decreased in the post-weaned stage. High values were estimated for the maternal permanent environmental effect, possibly because of the extensive grassland that was available. High total heritability values were estimated for all of the traits evaluated. Significant, positive correlations were found between direct and maternal additive genetic traits with a gradual decrease as the lambs gained independence from their mothers. The genetic trends observed were irregular and incremental. Significant genetic variance suggests that direct selection for pre-weaning traits results in indirect selection of maternal abilities, and individual selection of any post-weaning trait results in rapid genetic improvement.

Inclusion of cytoplasmic lineage effect and direct-maternal genetic covariance for genetic evaluation of growth traits in Nellore cattle.

We evaluated the impact of cytoplasmic lineage effects (Lc) for growth traits on genetic evaluation, including the genetic covariance between direct and maternal effects (σam). Pedigree data from 496,190 Nellore animals and observations on birth weight (BW, N = 243,391), weaning weight (WW, N = 431,681), and post-weaning weight gain adjusted to 345 days (PWG, N = 172,131) were analyzed. Four univariate models were used to obtain estimates of (co)variance components using the restricted maximum likelihood method in the BLUPF90 program. Model 1 included Lc and σam. Model 2 included Lc and σam was set to zero. Model 3 did not include Lc. Model 4 did not include Lc and σam was set to zero. These models considered the effects of the Lc as random. Phenotypic variance obtained through cytoplasmic lineage effects was determined for all traits, ranging from 0.07 to 0.15, 0.15 to 0.03, and 0.05 to 0.03% for BW, WW, and PWG, respectively, for models 1 and 2. Correlations between direct and maternal genetic components were positive for WW and negative for BW and PWG. No differences were observed for genetic parameter estimates or animal ranking with the inclusion of σam. For BW, the likelihood ratio suggested that model 1 best fits the data, while model 4 was the most appropriate for WW and PWG. Thus, these models are recommended for genetic evaluations. Despite the low magnitude of cytoplasmic lineages, this effect could predict breeding value and improve the selection of animals for BW in this Nellore population.

Genetic variation of the weaning weight of beef cattle as a function of accumulated heat stress.

The objective of this study was to identify the genetic variation in the weaning weight (WW) of beef cattle as a function of heat stress. The WWs were recorded at approximately 205 days of age in three Brazilian beef cattle populations: Nelore (93,616), Brangus (18,906) and Tropical Composite (62,679). In view of the cumulative nature of WW, the effect of heat stress was considered as the accumulation of temperature and humidity index units (ACTHI) from the animal's birth to weaning. A reaction norm model was used to estimate the (co)variance components of WW across the ACTHI scale. The accumulation of THI units from birth to weaning negatively affected the WW. The definition of accumulated THI units as an environmental descriptor permitted to identify important genetic variation in the WW as a function of heat stress. As evidence of genotype by environment interaction, substantial heterogeneity was observed in the (co)variance components for WW across the environmental gradient. In this respect, the best animals in less stressful environments are not necessarily the best animals in more stressful environments. Furthermore, the response to selection for WW is expected to be lower in more stressful environments.

Comparison of bivariate and multivariate joint analyses on the selection loss of beef cattle.

For genetic evaluation of beef cattle, univariate or bivariate analyses are often performed as an alternative to decrease the complexity of matrices and mathematical models compared to multivariate analysis, which considers a larger number of joint traits. The use of bivariate methods to calculate genetic predictors may cause bias in the estimation of breeding values and, as a consequence, reclassification of the rank of top-selected sires, resulting in a loss of genetic gain in future generations. The objective of this study was to compare the bivariate and multivariate joint methods of genetic evaluation, verifying the selection loss, and reclassification of the ranking of the best animals with different selection intensities. Records of 431,224 Nellore breed animals were evaluated for birth weight, weaning weight, post-weaning gain, muscle score, scrotal circumference, and selection index. The pedigree file consisted of 505,848 animals, including 218,727 males and 287,121 females. The predicted breeding values were obtained using the program PEST 2, and the complete pedigree analysis was performed by the PopReport software. The results showed that, for the four different selection intensities considered (TOP 10 and 1, 10, and 30%), selection loss and reclassification of animals in ranking, were detected for all traits evaluated when the two methods of analysis were compared.

Use of weaning management group as a random effect for a more robust estimation of genetic parameters for post-weaning traits in Nellore cattle.

Data from 69,525 animals were used to compare two types of analyses, one of them having the weaning management group (WEMANG) included as an effect in the contemporary group (F_WEMANG) and the other considering the weaning management group as a random effect, not related to the mathematical model (R_WEMANG) for post-weaning traits. The components of (co)variance were estimated for pre-weaning traits (birth weight and weaning weight) and for post-weaning traits [scrotal circumference (SC), weight gain from weaning to 18 months of age (WG) and muscle score (MUSC)] in Nellore cattle, based on a complete animal model. Heritability of SC, WG and MUSC for the F_WEMANG model was equal to 0.46 ± 0.02, 0.38 ± 0.03 and 0.26 ± 0.01, and for the R_WEMANG model it was 0.45 ± 0.02, 0.31 ± 0.03 and 0.25 ± 0.01, respectively. Genetic correlations between all the studied traits varied between 0.07 ± 0.01 and 0.77 ± 0.03 in F_WEMANG and between 0.02 ± 0.01 and 0.76 ± 0.04 in R_WEMANG. The R_ WEMANG model allowed a decrease in the number of contemporary groups as well as an increase in the number of observations per group without significant alterations in heritability coefficients, for the post-weaning traits. Consequently, the analysis became more robust and avoided having contemporary groups with low variability.

Genotype by environment interaction for post-weaning weight gain, scrotal circumference, and muscling score of composite beef cattle in different regions of Brazil.

The objectives of this study were to characterize and define homogenous production environments of composite beef cattle in Brazil in terms of climatic and geographic variables by using multivariate exploratory techniques; to evaluate the presence of genotype by environment interaction (GxE) for post-weaning weight gain (PWG), yearling scrotal circumference (SC), and yearling muscling (MUS). Hierarchical and nonhierarchical cluster analysis was used to group farms located in regions with similar environmental variables into clusters. Six clusters of farms were formed. The effect of sire-cluster interaction was tested by single-trait analysis. Genetic parameters were estimated by multi-trait analysis considering the same trait to be different in each cluster. The effect of sire-cluster interaction was significant (P <0.01) for PWG and MUS. Estimates of genetic correlations among clusters ranged from 0.31 to 0.93 for PWG, 0.64 to 0.89 for SC, and 0.18 to 0.80 for MUS. These results indicate the need for a genetic analysis on a regional basis or inclusion of the GxE effect in the statistical model to permit appropriate evaluation of the animals.

Sire effects on carcass and meat quality traits of young Nellore bulls.

Meat quality is being increasingly demanded by consumers in recent years. Several factors can affect meat quality, ranging from animal traits such as breed and genetic heritage to pre- and post-slaughter processes. This study investigated the influence of Nellore bulls on carcass and meat quality traits. We used 475 young uncastrated males, the progeny of 54 bulls, to evaluate characteristics of the following carcass traits: hot carcass weight, rib-eye area, and fat thickness. We also evaluated the following beef quality traits: marbling, color, drip loss, cooking loss, and shear force at 0, 7, and 14 days of aging. Bulls had a significant influence (P≤0.05) on rib-eye area, fat thickness, marbling, drip loss at 14 days of aging and color at all aging periods. Based on these results, the use of bulls with high breeding values for these traits can provide important advances in carcass traits and meat quality in breeding programs of Nellore cattle that are raised in tropical conditions.

Quantitative genetic study of age at subsequent rebreeding in Nellore cattle by using survival analysis.

The continuous trait age at subsequent rebreeding (ASR) was evaluated using survival analysis in Nellore breed cows that conceived for the first time at approximately 14 months of age. This methodology was chosen because the restricted breeding season produces censored data. The dataset contained 2885 records of ASR (in days). Records of females that did not produce calves in the following year after being exposed to a sire were considered censored (48.3% of the total). The statistical model used was a Weibull mixed survival model, which included fixed effects of contemporary groups (CG) and period and a random effect of individual animal. The effect of contemporary groups on ASR was significant (P < 0.01). Heritabilities obtained for ASR were 0.03 and 0.04 in logarithmic and original scales, respectively. These results indicate that the genetic selection response for subsequent reproduction of 2-year-old Nellore breed females is not expected to be effective based on survival analysis. Furthermore, these results suggest that environmental improvement is fundamental to this important trait. It should be highlighted that an increase in the average date of birth can produce an adverse effect in the future, since this cannot be compensated by genetic improvement.

Genome-wide association study for feedlot average daily gain in Nellore cattle (Bos indicus).

The genome-wide association study (GWAS) results are presented for average daily gain (ADG) in Nellore cattle. Phenotype of 720 male Bos indicus animals with information of ADG in feedlots and 354,147 single-nucleotide polymorphisms (SNPs) obtained from a database added by information from Illumina Bovine HD (777,962 SNPs) and Illumina BovineSNP50 (54,609) by imputation were used. After quality control and imputation, 290,620 SNPs remained in the association analysis, using R package Genome-wide Rapid Association using Mixed Model and Regression method GRAMMAR-Gamma. A genomic region with six significant SNPs, at Bonferroni-corrected significance, was found on chromosome 3. The most significant SNP (rs42518459, BTA3: 85849977, p = 9.49 × 10(-8)) explained 5.62% of the phenotypic variance and had the allele substitution effect of -0.269 kg/day. Important genes such as PDE4B, LEPR, CYP2J2 and FGGY are located near this region, which is overlapped by 12 quantitative trait locus (QTLs) described for several production traits. Other regions with markers with suggestive effects were identified in BTA6 and BTA10. This study showed regions with major effects on ADG in Bos indicus in feedlots. This information may be useful to increase the efficiency of selecting this trait and to understand the physiological processes involved in its regulation.

Genotype by environment interaction and model comparison for growth traits of Santa Ines sheep.

The objectives of the present study were to compare alternative models for the genetic evaluation and assess the importance of genotype by environment interaction (G×E) in the estimation of genetic parameters and genetic evaluation of birth weight (BW), weight at 60 days of age (W60) and weight at 180 days of age (W180) of Santa Ines sheep. Data comprise 7622 BW, 4673 W60 and 2830 W180 records from animals born in 44 Brazilian herds. Four models were used for the analyses: animal model (AM) with homogeneous residual variance (1), or heterogeneous residual variance (2), hierarchical reaction norms model (HRNM) with homogeneous (1) or heterogeneous residual variance (2). The models that best fit the BW, W60 and W180 data were AM2, HRNM1 and HRNM2 respectively. Thus, models for genetic evaluation that consider heterogeneity of variances are recommended to evaluate growth traits of sheep. The correlation between intercept and slope of the HRNM was higher than 0.70 for all traits studied, indicating that animals with higher average breeding values responded better to improvement in environmental conditions, a fact characterizing the scale effect of G×E. Therefore, G×E is an important factor to be considered in the estimation of genetic parameters and genetic evaluation of growth traits of sheep.

Heritability estimate and genetic correlations of reproductive features in Nellore bulls, offspring of super precocious, precocious and normal cows under extensive farming conditions.

The present work aimed to estimate heritability and genetic correlations of reproductive features of Nellore bulls, offspring of mothers classified as superprecocious (M1), precocious (M2) and normal (M3). Twenty one thousand hundred and eighty-six animals with average age of 21.29 months were used, evaluated through the breeding soundness evaluation from 1999 to 2008. The breeding soundness features included physical semen evaluation (progressive sperm motility and sperm vigour), semen morphology (major, minor and total sperm defects), scrotal circumference (SC), testicular volume (TV) and SC at 18 months of age (SC18). The components of variance, heritability and genetic correlations for and between the features were estimated simultaneously by restricted maximum likelihood, with the use of the vce software system vs 6. The heritability estimates were high for SC18, SC and TV (0.43, 0.63 and 0.54; 0.45, 0.45 and 0.44; 0.42, 0.45 and 0.41, respectively for the categories of mothers M1, M2 and M3) and low for physical and morphological semen aspects. The genetic correlations between SC18 and SC were high, as well as between these variables with TV. High and positive genetic correlations were recorded among SC18, SC and TV with the physical aspects of the semen, although no favourable association was verified with the morphological aspects, for the three categories of mothers. It can be concluded that the mother's sexual precocity did not affect the heritability of their offspring reproduction features.

Association of single nucleotide polymorphisms in the bovine leptin and leptin receptor genes with growth and ultrasound carcass traits in Nellore cattle.

Given the important role of leptin in metabolism, we looked for a possible association of leptin and leptin receptor polymorphisms with carcass and growth traits in Nellore cattle. We examined associations of leptin and leptin receptor SNPs with ultrasound carcass (longissimus dorsi muscle area (ribeye area), backfat thickness and rump fat thickness and growth traits (weaning weight adjusted to 210 days of age, yearling weight adjusted to 550 days of age, weight gain of weaning to yearling and scrotal circumference adjusted to 550 days of age) of 2162 Bos primigenius indicus (Nellore) animals. Allele and genotypic frequencies were calculated for each marker. Allele substitution, additive and dominance effects of the polymorphisms were also evaluated. Some alleles of the molecular markers had low frequencies, lower than 1%, in the sample analyzed, although the same polymorphisms described for B. p. taurus cattle were found. Due to very low allelic frequencies, the E2JW, A59V and UASMS2 markers were not included in the analysis, because they were almost fixed. E2FB was found to be significantly associated with weight gain, ribeye area and backfat thickness. The promoter region markers, C963T and UASMS1, were also found to be significantly associated with ribeye area. T945M was significantly associated with weight gain. We conclude that the leptin and receptor gene markers would be useful for marker-assisted selection.

Single nucleotide polymorphisms in CAPN and leptin genes associated with meat color and tenderness in Nellore cattle.

We analyzed single nucleotide polymorphisms in calpain, leptin, leptin receptor, and growth hormone receptor genes and their association with color, drip and cooking losses of longissimus muscle at 7, 14 and 21 days postmortem in 638 purebred Nellore bulls slaughtered between 22 and 26 months of age. Meat samples were vacuum-packed and aged at 4°C. The single nucleotide polymorphisms T945M, GHR2, E2FB, and CAPN4751 were evaluated. All genotypic classes were observed; however, the T/T genotype of T945M and E2FB was found at a low frequency. A significant association of E2FB with drip loss (a measure of water-holding capacity) was detected at seven days of meat aging. CAPN4751 had an additive effect on red and yellow color intensities. The T allele of CAPN4751 was found to be positively associated with improved meat color, but not with meat tenderness, differing from a previous report indicating that it is associated with meat tenderness. We conclude that the potential for use of CAPN4751 as a marker for these meat quality traits requires further research.

Genetic parameters for postweaning traits in Braunvieh cattle.

Genetic parameters for traits related to postweaning growth in Braunvieh cattle, reared under tropical and sub-tropical conditions in Brazil, were studied. Weight traits were weight at 365 days of age (W365, N = 4055), at 450 days (W450, N = 3453), and at 550 days (W550, N = 1946), while weight gains were gain from weaning to 365 days of age (WGW365, N = 3060), from weaning to 450 days (WGW450, N = 2764), from weaning to 550 days (WGW550, N = 1531), from 365 to 550 days of age (WG365550, N = 1528), from 365 to 450 days (WG365450, N = 2401), and from 450 to 550 days (WG450550, N = 1563). A full animal model was used for estimating the variance components, using the MTDFREML software. The dataset contained 18,688 animals with phenotypic measures and 35,188 animals in the relationship matrix. Heritability estimates for postweaning weights decreased with age. For W365, W450 and W550, respectively, the direct heritability estimates were 0.29 +/- 0.061, 0.25 +/- 0.057, 0.16 +/- 0.060, maternal heritability was 0.20 +/- 0.035, 0.18 +/- 0.035, 0.13 +/- 0.052, and total heritability was 0.30, 0.35, 0.26. In this breed, maternal influence was found to be important up to 550 days of age. The greater genetic correlations between weights were observed for weights measured at shorter intervals. A large environmental effect was observed for weight gain between weaning and 550 days; this effect was greater for the gains between 365 and 550 days.

Additive genetic relationship of longevity with fertility and production traits in Nellore cattle based on bivariate models.

Survival or longevity is an economically important trait in beef cattle. The main inconvenience for its inclusion in selection criteria is delayed recording of phenotypic data and the high computational demand for including survival in proportional hazard models. Thus, identification of a longevity-correlated trait that could be recorded early in life would be very useful for selection purposes. We estimated the genetic relationship of survival with productive and reproductive traits in Nellore cattle, including weaning weight (WW), post-weaning growth (PWG), muscularity (MUSC), scrotal circumference at 18 months (SC18), and heifer pregnancy (HP). Survival was measured in discrete time intervals and modeled through a sequential threshold model. Five independent bivariate Bayesian analyses were performed, accounting for cow survival and the five productive and reproductive traits. Posterior mean estimates for heritability (standard deviation in parentheses) were 0.55 (0.01) for WW, 0.25 (0.01) for PWG, 0.23 (0.01) for MUSC, and 0.48 (0.01) for SC18. The posterior mean estimates (95% confidence interval in parentheses) for the genetic correlation with survival were 0.16 (0.13-0.19), 0.30 (0.25-0.34), 0.31 (0.25- 0.36), 0.07 (0.02-0.12), and 0.82 (0.78-0.86) for WW, PWG, MUSC, SC18, and HP, respectively. Based on the high genetic correlation and heritability (0.54) posterior mean estimates for HP, the expected progeny difference for HP can be used to select bulls for longevity, as well as for post-weaning gain and muscle score.

Association of SNPs on CAPN1 and CAST genes with tenderness in Nellore cattle.

We examined whether single-nucleotide polymorphisms (SNPs) in the calpain (CAPN) and calpastatin (CAST) genes, described from Bos primigenius taurus, are polymorphic in Nellore cattle. We also looked for a possible association of linkage disequilibrium of this polymorphism with tenderness of the longissimus dorsi muscle after 7, 14 and 21 days of postmortem aging in 638 purebred Nellore bulls. Meat tenderness was measured as Warner-Bratzler shear force. Additive and dominance effects were tested for SNPs of the three genotypic classes; the substitution effect was tested for SNPs with missing genotypic classes. Genotypic and gene frequencies were also calculated for the different SNPs. An increase in tenderness was observed from 7 to 21 days; the average values for shear force at 7, 14 and 21 days of aging were 5.92 +/- 0.06, 4.92 +/- 0.05, and 4.38 +/- 0.04 kg, respectively. All markers showed polymorphism, but there was no CC genotype for CAPN316, and few animals showed the AA genotype for CAPN530. The alleles CAPN4751, UOGCAST1, and WSUCAST were found to have additive and dominance effects for shear force at 7, 14 and 21 days, while CAPN316 showed a substitution effect for shear force at 7 and 21 days. An additive-by-additive epistatic interaction was observed between CAPN4751 and markers on the CAST gene. In conclusion, these markers should be considered for use in breeding programs.

Comparison of different models to estimate genetic parameters for carcass traits in a commercial broiler line.

The effect of genetic and non-genetic factors for carcass, breast meat and leg weights, and yields of a commercial broiler line were investigated using the restricted maximum likelihood method, considering four different animal models, including or excluding maternal genetic effect with covariance between direct and maternal genetic effects, and maternal permanent environmental effect. The likelihood ratio test was used to determine the most adequate model for each trait. For carcass, breast, and leg weight, and for carcass and breast yield, maternal genetic and permanent environmental effects as well as the covariance between direct and maternal genetic effects were significant. The estimates of direct and maternal heritability were 0.17 and 0.04 for carcass weight, 0.26 and 0.06 for breast weight, 0.22 and 0.02 for leg weight, 0.32 and 0.02 for carcass yield, and 0.52 and 0.04 for breast yield, respectively. For leg yield, maternal permanent environmental effect was important, in addition to direct genetic effects. For that trait, direct heritability and maternal permanent environmental variance as a proportion of the phenotypic variance were 0.43 and 0.02, respectively. The results indicate that ignoring maternal effects in the models, even though they were of small magnitude (0.02 to 0.06), tended to overestimate direct genetic variance and heritability for all traits.