Study of using marker assisted selection on a beef cattle breeding program by model comparison

A data set of a commercial Nellore beef cattle selection program was used to compare breeding models that assumed or not markers effects to estimate the breeding values, when a reduced number of animals have phenotypic, genotypic and pedigree information available. This herd complete data set was composed of 83,404 animals measured for weaning weight (WW), post-weaning gain (PWG), scrotal circumference (SC) and muscle score (MS), corresponding to 116,652 animals in the relationship matrix. Single trait analyses were performed by MTDFREML software to estimate fixed and random effects solutions using this complete data. The additive effects estimated were assumed as the reference breeding values for those animals. The individual observed phenotype of each trait was adjusted for fixed and random effects solutions, except for direct additive effects. The adjusted phenotype composed of the additive and residual parts of observed phenotype was used as dependent variable for models' comparison. Among all measured animals of this herd, only 3160 animals were genotyped for 106 SNP markers. Three models were compared in terms of changes on animals' rank, global fit and predictive ability. Model 1 included only polygenic effects, model 2 included only markers effects and model 3 included both polygenic and markers effects. Bayesian inference via Markov chain Monte Carlo methods performed by TM software was used to analyze the data for model comparison. Two different priors were adopted for markers effects in models 2 and 3, the first prior assumed was a uniform distribution (U) and, as a second prior, was assumed that markers effects were distributed as normal (N). Higher rank correlation coefficients were observed for models 3_U and 3_N, indicating a greater similarity of these models animals' rank and the rank based on the reference breeding values. Model 3_N presented a better global fit, as demonstrated by its low DIC. The best models in terms of predictive ability were models 1 and 3_N. Differences due prior assumed to markers effects in models 2 and 3 could be attributed to the better ability of normal prior in handle with collinear effects. The models 2_U and 2_N presented the worst performance, indicating that this small set of markers should not be used to genetically evaluate animals with no data, since its predictive ability is restricted. In conclusion, model 3_N presented a slight superiority when a reduce number of animals have phenotypic, genotypic and pedigree information. It could be attributed to the variation retained by markers and polygenic effects assumed together and the normal prior assumed to markers effects, that deals better with the collinearity between markers. (C) 2012 Elsevier B.V. All rights reserved.

Estimation methods of non-additive effects for characteristics of weight and scrotal circumference in crossbred beef cattle

The objective of this study was to investigate, in a population of crossbred cattle, the obtainment of the non-additive genetic effects for the characteristics weight at 205 and 390 days and scrotal circumference, and to evaluate the consideration of these effects in the prediction of breeding values of sires using different estimation methodologies. In method 1, the data were pre-adjusted for the non-additive effects obtained by least squares means method in a model that considered the direct additive, maternal and non-additive fixed genetic effects, the direct and total maternal heterozygosities, and epistasis. In method 2, the non-additive effects were considered covariates in genetic model. Genetic values for adjusted and non-adjusted data were predicted considering additive direct and maternal effects, and for weight at 205 days, also the permanent environmental effect, as random effects in the model. The breeding values of the categories of sires considered for the weight characteristic at 205 days were organized in files, in order to verify alterations in the magnitude of the predictions and ranking of animals in the two methods of correction data for the non-additives effects. The non-additive effects were not similar in magnitude and direction in the two estimation methods used, nor for the characteristics evaluated. Pearson and Spearman correlations between breeding values were higher than 0.94, and the use of different methods does not imply changes in the selection of animals.

Effect of supplementation of two sources and two levels of copper on lipid metabolism in Nellore beef cattle

This study was conducted with 35 Nellore beef cattle to determine the effect of supplementation of two levels and two copper sources (organic and inorganic) on metabolism of lipids and cholesterol of meat. The five treatments used were: Control: without copper supplementation, 110 or 140: 10 or 40 mg/kg DM (as Cu sulfate), O10 or O40: 10 or 40 mg/kg DM (as Cu proteinate). In general, the copper supplementation changed the fatty acid profile of meat (p < 0.05), with a higher proportion of unsaturated fatty acids and reduction of saturated fatty acids. There was no effect of supplementation on blood cholesterol and triglycerides, however; in general, there was a reduction in cholesterol concentration in the L dorsi (p < 0.05) compared to the control treatment through the reduction (p < 0.05) in the concentrations of GSH and GSH/GSSG ratio. The Cu supplementation did have an influence on metabolism of lipids. The production of healthier meat is beneficial to public health by reducing the risk of cardiovascular disease. (C) 2012 Elsevier Ltd. All rights reserved.

Trends of the genetic connectedness measures among Nelore beef cattle herds

Validity of comparisons between expected breeding values obtained from best linear unbiased prediction procedures in genetic evaluations is dependent on genetic connectedness among herds. Different cattle breeding programmes have their own particular features that distinguish their database structure and can affect connectedness. Thus, the evolution of these programmes can also alter the connectedness measures. This study analysed the evolution of the genetic connectedness measures among Brazilian Nelore cattle herds from 1999 to 2008, using the French Criterion of Admission to the group of Connected Herds (CACO) method, based on coefficients of determination (CD) of contrasts. Genetic connectedness levels were analysed by using simple and multiple regression analyses on herd descriptors to understand their relationship and their temporal trends from the 19992003 to the 20042008 period. The results showed a high level of genetic connectedness, with CACO estimates higher than 0.4 for the majority of them. Evaluation of the last 5-year period showed only a small increase in average CACO measures compared with the first 5 years, from 0.77 to 0.80. The percentage of herds with CACO estimates lower than 0.7 decreased from 27.5% in the first period to 16.2% in the last one. The connectedness measures were correlated with percentage of progeny from connecting sires, and the artificial insemination spread among Brazilian herds in recent years. But changes in connectedness levels were shown to be more complex, and their complete explanation cannot consider only herd descriptors. They involve more comprehensive changes in the relationship matrix, which can be only fully expressed by the CD of contrasts.

Genotype by environment interaction for birth and weaning weights of composite beef cattle in different regions of Brazil

The objectives of the present study were to characterize and define homogenous production environments of composite beef cattle in Brazil in terms of climatic and geographic variables using multivariate exploratory techniques and to use them to assess the presence of G x E for birth weight (BW) and weaning weight (WW). Data from animals born between 1995 and 2008 on 36 farms located in 27 municipalities of the Brazilian states were used. Fifteen years of climate observations (mean minimum and maximum annual temperature and mean annual rainfall) and geographic (latitude, longitude and altitude) data were obtained for each municipality where the farms were located for characterization of the production environments. 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 using deviance information criterion (DIC). Genetic parameters were estimated by multi-trait analysis considering the same trait to be different in each cluster. According to the values of DIC, the inclusion of sire-cluster effect did not improve the fit of the genetic evaluation model for BW and WW. Estimates of genetic correlations among clusters ranged from -0.02 to 0.92. The low genetic correlation among the most studied regions permits us to suggest that a separate genetic evaluation for some regions should be undertaken. (C) 2012 Elsevier B.V. All rights reserved.

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.

Interaction of dietary vitamin D-3 and sunlight exposure on B. indicus cattle: Animal performance, carcass traits, and meat quality

Attempts to improve beef tenderness through supplementation with dietary vitamin D-3 have been challenged by null results and negative impacts on animal performance and carcass traits. Because vitamin D-3 is also synthesised by the animal via ultraviolet radiation from sunlight, the effectiveness of supplementation with dietary vitamin D-3 may be modulated by the degree of exposure of the animal to sunlight. Hence, this work aimed to verify whether dietary vitamin D-3 modifies meat quality without negatively affecting animal performance and carcass traits in B. indicus beef cattle that were either exposed to or protected from natural sunlight. Forty-two (411 +/- 38 kg) Nellore-type castrated males were fed a high-concentrate diet for 45 days after assignment to a treatment group. The treatments comprised combinations of three levels of vitamin D3 [ViTD - none (V0) or 2 x 10(6) IU of vitamin D-3 administered for either 2 (V2) or 8 (V8) consecutive days pre-slaughter] and two shading conditions (SHADE - unshaded or shaded). The post-mortem (pm) measurements were taken in the Longissimus thoracis et lumborum muscle. The animal performance and carcass traits were unaffected by ViTD or SHADE The V2 treatment increased the Myofibrillar Fragmentation Index in shaded animals compared to unshaded ones. Animals under shade had higher muscle calcium concentration. There was no effect of either ViTD or SHADE on the shear force. The L* values were higher at 24 h pm than at 0 and 1 h pm, with no differences among the animals in the ViTD or SHADE groups. Higher a* values were observed among animals in the V8 group than in the V0 group, and higher b* values were observed among animals in the V8 group than in the V2 or V0 groups, which were not different. In conclusion, ViTD and SHADE did not affect animal performance, carcass traits or shear force, whereas animals receiving a lower ViTD dosage and SHADE exhibited altered myofibrillar fragmentation. ViTD affected the colour parameters, and changes in the lightness of the beef related to the time pm were found in meat from animals under SHADE. (c) 2012 Elsevier B.V. All rights reserved.

Genetic associations between flight speed and growth traits in Nellore cattle

The aim of the present study was to estimate genetic parameters for flight speed and its association with growth traits in Nellore beef cattle. The flight speed (FS) of 7,402 yearling animals was measured, using a device composed of a pair of photoelectric cells. Time interval data (s) were converted to speed (m/s) and faster animals were regarded as more reactive. The growth traits analyzed were weaning weight (WW), ADG from weaning to yearling age, and yearling scrotal circumference (SC). The (co)variance components were estimated using REML in a multitrait analysis applying an animal model. The model included random direct additive genetic and residual effects, fixed effects of contemporary groups, age of dam (classes), and age of animal as covariable. For WW, the model also included maternal genetic and permanent environmental random effects. The direct heritability estimate for FS was 0.26 +/- 0.05 and direct heritability estimates for WW, SC, and ADG were 0.30 +/- 0.01, 0.48 +/- 0.02, and 0.19 +/- 0.01, respectively. Estimates of the genetic correlation between FS and the growth traits were -0.12 +/- 0.07 (WW), -0.13 +/- 0.08 (ADG), and -0.11 +/- 0.07 (SC). Although the values were low, these correlations showed that animals with better temperaments (slower FS) tended to present better performance. It is possible to infer that long-term selection for weight and scrotal circumference can promote a positive genetic response in the temperament of animals. Nevertheless, to obtain faster genetic progress in temperament, it would be necessary to perform direct selection for such trait. Flight speed is an easily measured indicator of temperament and can be included as a selection criterion in breeding programs for Nellore cattle.

Genetic relationship between growth and reproductive traits in Nellore cattle

The objective of this study was to evaluate the genetic relationship between postweaning weight gain (PWG), heifer pregnancy (HP), scrotal circumference (SC) at 18 months of age, stayability at 6 years of age (STAY) and finishing visual score at 18 months of age (PREC), and to determine the potential of these traits as selection criteria for the genetic improvement of growth and reproduction in Nellore cattle. The HP was defined as the observation that a heifer conceived and remained pregnant, which was assessed by rectal palpation at 60 days. The STAY was defined as whether or not a cow calved every year up to the age of 6 years, given that she was provided the opportunity to breed. The Bayesian linear-threshold analysis via the Gibbs sampler was used to estimate the variance and covariance components applying a multitrait model. Posterior mean estimates of direct heritability were 0.15 +/- 0.00, 0.42 +/- 0.02, 0.49 +/- 0.01, 0.11 +/- 0.01 and 0.19 +/- 0.00 for PWG, HP, SC, STAY and PREC, respectively. The genetic correlations between traits ranged from 0.17 to 0.62. The traits studied generally have potential for use as selection criteria in genetic breeding programs. The genetic correlations between all traits show that selection for one of these traits does not imply the loss of the others.

Random regression analyses using B-spline functions to model growth of Nellore cattle

The objective of this study was to estimate (co)variance components using random regression on B-spline functions to weight records obtained from birth to adulthood. A total of 82 064 weight records of 8145 females obtained from the data bank of the Nellore Breeding Program (PMGRN/Nellore Brazil) which started in 1987, were used. The models included direct additive and maternal genetic effects and animal and maternal permanent environmental effects as random. Contemporary group and dam age at calving (linear and quadratic effect) were included as fixed effects, and orthogonal Legendre polynomials of age (cubic regression) were considered as random covariate. The random effects were modeled using B-spline functions considering linear, quadratic and cubic polynomials for each individual segment. Residual variances were grouped in five age classes. Direct additive genetic and animal permanent environmental effects were modeled using up to seven knots (six segments). A single segment with two knots at the end points of the curve was used for the estimation of maternal genetic and maternal permanent environmental effects. A total of 15 models were studied, with the number of parameters ranging from 17 to 81. The models that used B-splines were compared with multi-trait analyses with nine weight traits and to a random regression model that used orthogonal Legendre polynomials. A model fitting quadratic B-splines, with four knots or three segments for direct additive genetic effect and animal permanent environmental effect and two knots for maternal additive genetic effect and maternal permanent environmental effect, was the most appropriate and parsimonious model to describe the covariance structure of the data. Selection for higher weight, such as at young ages, should be performed taking into account an increase in mature cow weight. Particularly, this is important in most of Nellore beef cattle production systems, where the cow herd is maintained on range conditions. There is limited modification of the growth curve of Nellore cattle with respect to the aim of selecting them for rapid growth at young ages while maintaining constant adult weight.

Effects of two estradiol esters (benzoate and cypionate) on the induction of synchronized ovulations in Bos indicus cows submitted to a timed artificial insemination protocol

The effects of estradiol benzoate (EB) and estradiol cypionate (EC) on induction of ovulation after a synchronized LH surge and on fertility of Bos indicus females submitted to timed AI (TAI) were evaluated. In Experiment 1, ovariectomized Nelore heifers were used to evaluate the effect of EB (n = 5) and EC (n = 5) on the circulating LH profile. The LH surge timing (19.6 and 50.5 h; P = 0.001), magnitude (20.5 and 9.4 ng/mL; P = 0.005), duration (8.6 and 16.5 h; P = 0.001), and area under the LH curve (158.6 and 339.4 ng/mL; P = 0.01) differed between the EB and EC treatments, respectively. In Experiment 2 (follicular responses; n = 60) and 3 (pregnancy per AI; P/AI; n = 953) suckled Bos indicus beef cows submitted to an estradiol/progesterone-based synchronization protocol were assigned to receive one of two treatments to induce synchronized ovulation: 1 mg of EB im 24 h after progesterone (P4) device removal or 1 mg of EC im at P4 device removal. There was no difference (P > 0.05) between EB and EC treatments on follicular responses (maximum diameter of the ovulatory follicle, 13.1 vs. 13.9 mm; interval from progesterone device removal to ovulation, 70.2 vs. 68.5 h; and ovulation rate, 77.8 vs. 82.8%, respectively). In addition, P/AI was similar (P < 0.22) between the cows treated with EB (57.5%; 277/482) and EC (61.8%; 291/471). In conclusion, despite pharmacologic differences, both esters of estradiol administered either at P4 device removal (EC) or 24 h later (EB) were effective in inducing an LH surge which resulted in synchronized ovulations and similar P/AI in suckled Bos indicus beef cows submitted to TAI. (C) 2012 Elsevier Inc. All rights reserved.

MUC1 gene polymorphism in three Nelore lines selected for growth and its association with growth and carcass traits

The objective of this study was to describe the VNTR polymorphism of the mucin 1 gene (MUC1) in three Nelore lines selected for yearling weight to determine whether allele and genotype frequencies of this polymorphism were affected by selection for growth. In addition, the effects of the polymorphism on growth and carcass traits were evaluated. Birth, weaning and yearling weights, rump height, Longissimus muscle area, backfat thickness, and rump fat thickness, were analyzed. A total of 295 Nelore heifers from the Beef Cattle Research Center, Instituto de Zootecnia de Sertozinho, were used, including 41 of the control line, 102 of the selection line and 152 of the traditional. The selection and traditional lines comprise animals selected for higher yearling weight, whereas control line animals are selected for yearling weight close to the average. Five alleles were identified, with allele 1 being the most frequent in the three lines, especially in the lines selected for higher means for yearling weight. Heterozygosity was significantly higher in the control line. Association analyses showed significant effects of allele 1 on birth weight and weaning weight while the allele 3 exert significant effects on yearling weight and back fat thickness. Despite these findings, application of this marker to marker-assisted selection requires more consistent results based on the genotyping of a larger number of animals in order to increase the accuracy of the statistical analyses.

Pedigree analysis and inbreeding depression on growth traits in Brazilian Marchigiana and Bonsmara breeds

The study of population structure by pedigree analysis is useful to identify important circumstances that affect the genetic history of populations. The intensive use of a small number of superior individuals may reduce the genetic diversity of populations. This situation is very common for the beef cattle breeds. Therefore, the objectives of the present study were to analyze the pedigree and possible inbreeding depression on traits of economic interest in the Marchigiana and Bonsmara breeds and to test the inclusion of the individual inbreeding coefficient (F-i) or individual increases in inbreeding coefficient (Delta F-i) in the genetic evaluation model for the quantification of inbreeding depression. The complete pedigree file of the Marchigiana breed included 29,411 animals born between 1950 and 2003. For the Bonsmara breed, the pedigree file included 18,695 animals born between 1988 and 2006. Only animals with at least 2 equivalent generations of known pedigree were kept in the analyses of inbreeding effect on birth weight, weaning weight measured at about 205 d, and BW at 14 mo in the Marchigiana breed, and on birth weight, weaning weight, and scro-tal circumference measured at 12 mo in the Bonsmara breed. The degree of pedigree knowledge was greater for Marchigiana than for Bonsmara animals. The average generation interval was 7.02 and 3.19 for the Marchigiana and Bonsmara breed, respectively. The average inbreeding coefficient was 1.33% for Marchigiana and 0.26% for Bonsmara. The number of ancestors explaining 50% of the gene pool and effective population size computed via individual increase in coancestry were 13 and 97.79 for Marchigiana and 41 and 54.57 for Bonsmara, respectively. These estimates indicate reduction in genetic variability in both breeds. Inbreeding depression was observed for most of the growth traits. The model including Delta F-i can be considered more adequate to quantify inbreeding depression. The inclusion of F-i or Delta F-i in the genetic evaluation model may not result in better fit to the data. A genetic evaluation with simultaneous estimation of inbreeding depression can be performed in Marchigiana and Bonsmara breeds, providing additional information to producers and breeders.

A longitudinal study of Neospora caninum infection on three dairy farms in Brazil

Neospora caninum is an apicomplexan protozoan parasite that is one of the most important infectious causes of abortion in both dairy and beef cattle in many countries. The objectives of this longitudinal study were to determine the prevalence, rates of vertical and horizontal transmission of N. caninum and hazard for culling of N. caninum-seropositive animals in three Brazilian dairy herds. Blood samples from all animals were collected nine times at each of the three farms over a two-year period. Serum was tested for antibodies against N. caninum using the indirect fluorescent antibody test with a cutoff value of 1:100. The percentage of N. caninum-positive samples at each sampling time ranged at Farm 1 from 3.32% to 11.71%, at Farm 11 from 3.90% to 22.06% and at Farm III from 3.90% to 22.06%. The number of positive serological reactions varied in relation to the number of repeated samples taken from individual animals at each farm. In all herds, there was a high degree (P < 0.05) of association between the N. caninum serological status of dams and daughters. The seropositive conversion rate was estimated as 0.37%, 3.00% and 6.94% per 100 cow-years at Farms I, II and III, respectively. The seronegative conversion rate was estimated as 31.58% and 11.11% per 100 cow-years at Farms 1 and III, respectively. In all herds, there was no difference (P > 0.05) in the culling rate between the cattle that were seropositive cattle and seronegative for N. caninum infection. The results from this study confirm the importance of vertical transmission in the epidemiology of N. caninum. Although a few positive seroconversions indicated horizontal transmission, it does not appear to be the major route of infection for N. caninum. (C) 2012 Elsevier B.V. All rights reserved.

Degree of multicollinearity and variables involved in linear dependence in additive-dominant models

The objective of this work was to assess the degree of multicollinearity and to identify the variables involved in linear dependence relations in additive-dominant models. Data of birth weight (n=141,567), yearling weight (n=58,124), and scrotal circumference (n=20,371) of Montana Tropical composite cattle were used. Diagnosis of multicollinearity was based on the variance inflation factor (VIF) and on the evaluation of the condition indexes and eigenvalues from the correlation matrix among explanatory variables. The first model studied (RM) included the fixed effect of dam age class at calving and the covariates associated to the direct and maternal additive and non-additive effects. The second model (R) included all the effects of the RM model except the maternal additive effects. Multicollinearity was detected in both models for all traits considered, with VIF values of 1.03 - 70.20 for RM and 1.03 - 60.70 for R. Collinearity increased with the increase of variables in the model and the decrease in the number of observations, and it was classified as weak, with condition index values between 10.00 and 26.77. In general, the variables associated with additive and non-additive effects were involved in multicollinearity, partially due to the natural connection between these covariables as fractions of the biological types in breed composition.