Methionine supplementation during a hydrogen peroxide challenge alters components of insulin signaling and antioxidant proteins in subcutaneous adipose explants from dairy cows.

Enhanced postruminal supply of methionine (Met) during the peripartal period alters protein abundance of insulin, AA, and antioxidant signaling pathways in subcutaneous adipose tissue (SAT). Whether SAT is directly responsive to supply of Met and can induce molecular alterations is unknown. Our objective was to examine whether enhanced Met supply during an oxidative stress challenge in vitro alters insulin, AA, inflammation, and antioxidant signaling-related protein networks. Four late-lactation Holstein cows (average 27.0 kg of milk per day) were used for SAT collection. Tissue was incubated in duplicate for 4 h in a humidified incubator with 5% CO2 at 37°C according to the following experimental design: control medium with an "ideal" profile of essential AA (CTR; Lys:Met 2.9:1), CTR plus 100 µM H2O2 (HP), or CTR with greater Met supply plus 100 µM H2O2 (HPMET; Lys:Met 2.5:1). Molecular targets associated with insulin signaling, lipolysis, antioxidant nuclear factor, erythroid 2 like 2 (NFE2L2), inflammation, and AA metabolism were determined through reverse-transcription quantitative PCR and western blotting. Data were analyzed using the MIXED procedure of SAS 9.4 (SAS Institute Inc.). Among proteins associated with insulin signaling, compared with CTR, HP led to lower abundance of phosphorylated AKT serine/threonine kinase (p-AKT) and solute carrier family 2 member 4 (SLC2A4; insulin-induced glucose transporter). Although incubation with HPMET restored abundance of SLC2A4 to levels in the CTR and upregulated abundance of fatty acid synthase (FASN) and phosphorylated 5'-prime-AMP-activated protein kinase (p-AMPK), it did not alter p-AKT, which remained similar to HP. Among proteins associated with AA signaling, compared with CTR, challenge with HP led to lower abundance of phosphorylated mechanistic target of rapamycin (p-MTOR), and HPMET did not restore abundance to CTR levels. Among inflammation-related targets studied, incubation with HPMET led to greater protein abundance of nuclear factor kappa B subunit p65 (NFKB-RELA). The response in NFKB observed with HPMET was associated with a marked upregulation of the antioxidant transcription regulator NFE2L2 and the antioxidant enzyme glutathione peroxidase 1 (GPX1). No effects of treatment were detected for mRNA abundance of proinflammatory cytokines or antioxidant enzymes, underscoring the importance of post-transcriptional regulation. Overall, data indicated that short-term challenge with H2O2 was particularly effective in reducing insulin and AA signaling. Although a greater supply of Met had little effect on those pathways, it seemed to restore the protein abundance of the insulin-induced glucose transporter. Overall, the concomitant upregulation of key inflammation and antioxidant signaling proteins when a greater level of Met was supplemented to oxidant-challenged SAT highlighted the potential role of this AA in regulating the inflammatory response and oxidant status. Further studies should be conducted to assess the role of postruminal supply of Met and other AA in the regulation of immune, antioxidant, and metabolic systems in peripartal cow adipose tissue.

Genome-wide association study of resistance/susceptibility to infectious bovine keratoconjunctivitis in Brazilian Hereford cattle.

Genome-wide association studies were conducted to identify the more informative genomic regions and SNPs, as well as to identify candidate genes associated with infectious bovine keratoconjunctivitis (IBK) resistance/susceptibility in Hereford cattle. A Bayes B statistical approach was initially applied in genome-wide association studies by using deregressed estimated breeding values for IBK resistance/susceptibility. To estimate the combined effect of a genomic region that is potentially associated with QTL, 2504 non-overlapping 1-Mb windows that varied in SNP number were defined, with the most informative 24 windows including 427 SNPs and explaining more than 20% of the estimated genetic variance for IBK resistance/susceptibility. These regions were explored with respect to their biological functions through functional analysis to map potential candidate genes. The significant SNPs were mapped on chromosomes 1, 3, 5, 6, 7, 8, 10, 11, 12, 13, 14, 15, 18, 20, 23, and 28, and candidate genes were detected as related to the IBK. Most informative SNPs in term of genetic variance were located in proximity of genes related to phenotypic expression of lesions and biological processes associated to the IBK. Knowledge about phenotypic and genomic variation generated in the present study can be used to on design selection strategies to improve the resistance to IBK of Hereford cattle herds.

Short communication: Effect of glucose infusion dose and stage of lactation on glucose tolerance test kinetics in lactating dairy cows.

The objective for this study was to determine the effect of glucose dose and days following peak milk yield on plasma glucose, serum insulin, and plasma nonesterified fatty acids (NEFA) kinetics during an intravenous glucose tolerance test (IVGTT) in lactating dairy cattle. Six lactating Holstein dairy cows (3 primiparous and 3 multiparous) were assigned to 2 squares and received 0.092, 0.15, or 0.3 g of glucose/kg of body weight (BW) during an IVGTT at 74 and 221 d in milk (DIM), representing early (post-peak) lactation and mid lactation, respectively. Treatments were applied in a replicated Latin square design using contiguous 7-d periods within each stage of lactation. Milk production and dry matter intake were determined daily during the first 6 d of each period. The IVGTT was performed on d 7. For the IVGTT, cows were prepared with indwelling catheters in each jugular vein, and blood samples were collected at -15, -10, 5, 10, 15, 20, 30, 45, 60, 90, and 120 min relative to the glucose infusion. Samples were analyzed for plasma glucose, serum insulin, and plasma NEFA concentrations. Increasing the glucose dose during the IVGTT increased plasma glucose area under the curve (AUC), decreased glucose half-life, and increased maximal plasma glucose concentrations in plasma during the IVGTT. Greater glucose dose during the IVGTT elevated serum insulin AUC and increased nadir NEFA concentrations. Maximal plasma glucose concentration during the IVGTT was lower, whereas maximum NEFA concentration, NEFA AUC, and NEFA clearance rate were greater at 74 than at 221 DIM. Only glucose half-life was responsive to stage of lactation × glucose dose effects during the IVGTT, and the decrease in glucose half-life with increasing glucose dose was greater at 74 than at 221 DIM. Glucose AUC was greater and NEFA AUC lower for cows at 74 than at 221 DIM. For the doses tested, a glucose dose greater than 0.092 g/kg of BW resulted in peak blood glucose concentration that exceeded the previously reported renal glucose excretion threshold of 8.3 mM. There is a need for accompanying data to determine if this is the case for the glucose doses evaluated in this experiment. Based on maximal peak glucose concentrations and effects on glucose half-life, we identify 0.092 g of glucose/kg of BW (0.46 g/kg of metabolic body weight) as the preferred dose for the IVGTT for cows at 74 and 221 DIM in this study.

Multiple trait and random regression models using linear splines for genetic evaluation of multiple breed populations.

This study aimed to verify if random regression models using linear splines (RRMLS) are suitable for identifying genetic parameters in multiple-breed populations and also to investigate whether an interaction exists between the breeding value (BV) of sires and their progeny breed group. Ten populations were simulated by crossing 2 breeds with distinct genetic variance and nonzero segregation variance. To obtain the genetic parameters, 2 models were used: a multiple-trait model (MULT), in which the trait was considered distinct when evaluated in each group (1/2P1 + 1/2P2, 5/8P1 + 3/8P2, and 3/4P1 + 1/4P2), and a RRMLS with the spline polynomial knots adjusted to these same groups. The genetic parameters estimated through MULT and RRMLS did not differ from the simulated values. The correlations between BV (simulated and estimated) of animals were high and varied from 0.74 to 0.76, which indicates the efficiency of using MULT and RRMLS for predicting BV. Using field data, the traits age at first calving (AFC), first lactation length (LL), and 305-d milk yield (MY-305) from a multiple-breed population of Holstein-Gyr cattle were analyzed. The BV of animals were modeled through RRMLS with 3, 5, and 7 knots, distributed in accordance with the fraction of Holstein breed in each progeny breed group. It was verified that RRMLS with 7 knots for adjusting mean trajectories and genetic effects, with homogeneous residual variance, best fit AFC and LL. For MY-305, the best fit for mean trajectory and genetic effects was the RRMLS with 5 knots and with homogeneous residual variance. The posterior means of heritability varied from 0.21 to 0.48, 0.21 to 0.38, and 0.10 to 0.33 for AFC, LL, and MY-305, respectively. Estimates from genetic parameters obtained by using RRMLS with field data showed that this model is a useful tool for genetic evaluations of populations formed by a great number of breed groups. An interaction occurred between the BV of sires and their progeny breed group, and the genetic parameters for AFC, LL, and MY-305 traits from a multiple-breed population depend on breed composition of the progeny from which the evaluations are based.

Method for the estimation of genetic merit of animals with uncertain paternity under Bayesian inference.

The use of controlled mating or artificial insemination is impracticable in the case of large herds, mainly because of labour costs and the need to delimit areas during the breeding period. However, the exclusion of information from animals with uncertain paternity reduces genetic progress. The objectives of this study were as follows: (i) propose an iterative empirical Bayesian procedure to implement the hierarchical animal model (ITER); (ii) calculate the posterior probabilities of paternity by the maximum likelihood method following the concepts; (iii) compare an average numerator relationship matrix (ANRM), Bayesian hierarchical (HIER) models and ITER. Records of Nellore animals born between 1984 and 2006 from the zootechnical archive of Agropecuária Jacarezinho Ltda were used. For data consistency, records of contemporary groups (CGs) with fewer than three animals and animals whose records were 3.5 standard deviations above or below the mean of their CG were eliminated. After editing the data, 62,212 animals in the file and 12,876 animals in pedigree file were maintained, respectively. Spearman and Pearson correlations between the posterior mean of the genetic effects of animals were calculated to compare the ranking of animals for selection. Simulated data were used to confirm the veracity of the model. The correlations between ITER and HIER and between ITER and ANRM were similar evaluating different files, which decreased at the same proportion when only high-ranked animals were evaluated. In conclusion, the model proposed herein is a suitable computational alternative to improve the prediction of breeding values of animals in genetic evaluations using large databases, including animals with uncertain paternity.

Analyses of reaction norms reveal new chromosome regions associated with tick resistance in cattle.

Despite single nucleotide polymorphism (SNP) availability and frequent cost reduction has allowed genome-wide association studies even in complex traits as tick resistance, the use of this information source in SNP by environment interaction context is unknown for many economically important traits in cattle. We aimed at identifying putative genomic regions explaining differences in tick resistance in Hereford and Braford cattle under SNP by environment point of view as well as to identify candidate genes derived from outliers/significant markers. The environment was defined as contemporary group means of tick counts, since they seemed to be the most appropriate entities to describe the environmental gradient in beef cattle. A total of 4363 animals having tick counts (n=10 673) originated from 197 sires and 3966 dams were used. Genotypes were acquired on 3591 of these cattle. From top 1% SNPs (410) having the greatest effects in each environment, 75 were consistently relevant in all environments, which indicated SNP by environment interaction. The outliers/significant SNPs were mapped on chromosomes 1, 2, 5, 6, 7, 9, 11, 13, 14, 15, 16, 18, 21, 23, 24, 26 and 28, and potential candidate genes were detected across environments. The presence of SNP by environment interaction for tick resistance indicates that genetic expression of resistance depends upon tick burden. Markers with major portion of genetic variance explained across environments appeared to be close to genes with different direct or indirect functions related to immune system, inflammatory process and mechanisms of tissue destruction/repair, such as energy metabolism and cell differentiation.

Expression of paraoxonase types 1, 2 and 3 in reproductive tissues and activity of paraoxonase type 1 in the serum and seminal plasma of bulls.

The paraoxonases types 1, 2 and 3 (PON1, PON2 and PON3, respectively) are enzymes that degrade lipid peroxides, preventing oxidative damages relevant for male reproductive function. This study determined the expression of those three paraoxonases in reproductive tissues of bulls and evaluated correlations among the activity of PON1 in the serum and seminal plasma with breeding soundness parameters in bulls. The expression of PON1, PON2 and PON3 was characterised by RT-PCR in samples of testicular parenchyma, vesicular glands and epididymis collected from three slaughtered bulls. All three paraoxonases were expressed in the testicular parenchyma, PON2 and PON3 were both expressed in the epididymis head and PON3 was also expressed in the epididymis tail. The PON1 activity was determined in samples of serum and seminal plasma from 110 bulls submitted to breeding soundness evaluation. There was a strong correlation (r = .90) between the activity of the PON1 in both serum and seminal plasma (p < .0001). The PON1 activity in the seminal plasma was positively correlated with ejaculate's colour, sperm mass activity (p = .04), motility, vigour and viability (all p < .01). Thus, PON1 may be a potential marker for sperm motility and viability in bulls.

Genotype by environment interaction effects in genetic evaluation of preweaning gain for Line 1 Hereford cattle from Miles City, Montana.

It has long been recognized that genotype × environment interaction potentially influences genetic evaluation of beef cattle. However, this recognition has largely been ignored in systems for national cattle evaluation. The objective of this investigation was to determine if direct and maternal genetic effects on preweaning gain would be reranked depending on an environmental gradient as determined by year effects. Data used were from the 76-yr selection experiment with the Line 1 Hereford cattle raised at Miles City, MT. The data comprised recorded phenotypes from 7,566 animals and an additional 1,862 ancestral records included in the pedigree. The presence of genotype × environment interaction was examined using reaction norms wherein year effects on preweaning gain were hypothesized to linearly influence the EBV. Estimates of heritability for direct and maternal effects, given the average environment, were 10 ± 2 and 26 ± 3%, respectively. In an environment that is characterized by the 5th (95th) percentile of the distribution of year effects, the corresponding estimates of heritability were 18 ± 3 (22 ± 3%) and 30 ± 3% (30 ± 3%), respectively. Rank correlations of direct and maternal EBV appropriate to the 5th and 95th percentiles of the year effects were 0.67 and 0.92, respectively. In the average environment, the genetic trends were 255 ± 1 g/yr for direct effects and 557 ± 3 g/yr for maternal effects. In the fifth percentile environment, the corresponding estimates of genetic trend were 271 ± 1 and 540 ± 3 g/yr, respectively, and in the 95th percentile environment, they were 236 ± 1 and 578 ± 3 g/yr, respectively. Linear genetic trends in environmental sensitivity were observed for both the direct (-8.06 × 10 ± 0.49 × 10) and maternal (8.72 × 10 ± 0.43 × 10) effects. Therefore, changing systems of national cattle evaluation to more fully account for potential genotype × environment interaction would improve the assessment of breeding stock, particularly for direct effects. Estimates of environmental sensitivity parameters could also facilitate identification of genetic limitations to production.

Genomewide association study for production and meat quality traits in Canchim beef cattle.

The commercial value of the bovine carcass is determined by a set of traits, such as weight, yield, back fat thickness, and marbling; therefore, the genetic improvement of growth, meat, and carcass quality traits is an important tool to add value to the supply chain. Genomewide association studies (GWAS) enable the identification of loci that control phenotypic expression of quantitative traits (QTL). Therefore, the objective of this work was to perform a GWAS to identify genomic regions and genes associated with growth, carcass traits, and meat quality in Canchim beef cattle. These traits were yearling weight (YW), rib eye area (REA), back fat thickness (BFT), and marbling (MARB). To increase sample size and marker density, genotype imputation was performed, and only markers imputed with greater than 95% accuracy were used. Genomewide association study was performed using a Bayesian approach, by the Bayes B statistical method, incorporating genotypes and phenotypes from 614 animals from both the Canchim breed and the MA genetic group (offspring of Charolais bulls and one-half Canchim + one-half Zebu cows). This investigation identified 1 and 4 genomic regions explaining 0.23 and 7.35% of the genetic variance for REA and YW, respectively. These regions harbor a total of 19 genes, 7 of which were classified for biological functions by functional analysis. Significant associations were not observed for BFT and MARB. The identification of QTL that had been previously described in the literature reinforces associations found in this study.

Economic selection indexes for Hereford and Braford cattle raised in southern Brazil.

Economic selection indexes (EI) are considered the best way to select the most profitable animals for specific production systems. Nevertheless, in Brazil, few genetic evaluation programs deliver such indexes to their breeders. The aims of this study were to determine the breeding goals (BG) and economic values (EV, in US$) for typical beef cattle production systems in southern Brazil, to propose EI aimed to maximize profitability, and to compare the proposed EI with the currently used empirical index. Bioeconomic models were developed to characterize 3 typical production systems, identifying traits of economic impact and their respective EV. The first was called the calf-crop system and included the birth rate (BR), direct weaning weight (WWd), and mature cow weight (MCW) as selection goals. The second system was called the full-cycle system, and its breeding goals were BR, WWd, MCW, and carcass weight (CW). Finally, the third was called the stocking and finishing system, which had WWd and CW as breeding goals. To generate the EI, we adopted the selection criteria currently measured and used in the empirical index of PampaPlus, which is the genetic evaluation program of the Brazilian Hereford and Braford Association. The comparison between the EI and the current PampaPlus index was made by the aggregated genetic-economic gain per generation (Δ). Therefore, for each production system an index was developed using the derived economic weights, and it was compared with the current empirical index. The relative importance (RI) for BR, WWd, and MCW for the calf-crop system was 68.03%, 19.35%, and 12.62%, respectively. For the full-cycle system, the RI for BR, WWd, MCW, and CW were 69.63%, 7.31%, 5.01%, and 18.06%, respectively. For the stocking and finishing production system, the RI for WWd and CW was 34.20% and 65.80%, respectively. The Δ for the calf-crop system were US$6.12 and US$4.36, using the proposed economic and empirical indexes, respectively. Respective values were US$19.87 and US$18.22 for the full-cycle system and US$20.52 and US$18.52 in the stocking and finishing system. The efficiency of the proposed EI had low sensitivity to changes in the values of the economic and genetic parameters. The 3 EI generated higher Δ when using the proposed economic weight compared to the Δ provided by a PampaPlus index, suggesting the use of proposed EI to obtain greater economic profitability in relation to the current empirical PampaPlus index.

Design of selection schemes to include tick resistance in the breeding goal for Hereford and Braford cattle.

Ticks are one of the main causes of losses in cattle, causing economic impact by reducing productivity and fertility and by transmission of diseases. The objective of this study was to analyze the genetic gains obtained through different strategies to include traditional (EBV) or genomic EBV (GEBV) for tick count (TC) in selection indexes for Hereford and Braford cattle. Besides TC, we also considered traits currently included in the Delta G Breeding Program Index (DGI): preweaning gain, weaning conformation, weaning precocity, weaning muscling, postweaning gain, yearling conformation, yearling precocity, yearling muscling, and scrotal circumference. Genetic gain per generation (ΔG) was evaluated using the current DGI and including TC in 8 alternative scenarios with TC relative weightings of 10, 50, or 100% and using phenotype or GEBV. Genomic EBV accuracy () ranged between 0.1 and 0.9. As expected, increasing increases the accuracy of the index () for all scenarios in which GEBV were considered. As the relative weight of TC was increased to 50%, greater ΔG differences in relation to the baseline DGI ($53.03) scenario were observed when the GEBV information was included with equal to or greater than 0.7 only for TC (ΔG between $61.06 and $74.26) or equal to or greater than 0.5 for all traits (ΔG between $56.03 and $83.36). To achieve these accuracies for traits with low heritability, a large calibration data set would be required. Focusing only on TC, the availability of genomic information would be desirable to avoid the need to count ticks and the exposure of animals to parasitism risks. However, for = 0.7, the respective numbers for Hereford and Braford would be 4,703 and 6,522 animals. As expected, when comparing the relative index weights of 10, 50, and 100% for TC, the highest response to selection per generation (RS) for TC was in the scenario was with 100% relative weight and GEBV for this trait (SR = -0.09 SD with = 0.9). This would be the recommended scenario to form tick-resistant lines in Hereford and Braford cattle. However, with 50% relative weight for TC, including GEBV information for TC only or for all traits in index ( = 0.9), it should yield 93 or 84% of RS, respectively, compared to that obtained with full emphasis on TC (100% relative weight) and GEBV information. This indicates that in the presence of highly accurate GEBV, despite slightly slower gain for TC, indexes with 50% relative weight for TC are interesting alternatives to jointly improve tick resistance and other relevant traits.

Use of molecular markers to improve relationship information in the genetic evaluation of beef cattle tick resistance under pedigree-based models.

The selection of genetically superior individuals is conditional upon accurate breeding value predictions which, in turn, are highly depend on how precisely relationship is represented by pedigree. For that purpose, the numerator relationship matrix is essential as a priori information in mixed model equations. The presence of pedigree errors and/or the lack of relationship information affect the genetic gain because it reduces the correlation between the true and estimated breeding values. Thus, this study aimed to evaluate the effects of correcting the pedigree relationships using single-nucleotide polymorphism (SNP) markers on genetic evaluation accuracies for resistance of beef cattle to ticks. Tick count data from Hereford and Braford cattle breeds were used as phenotype. Genotyping was carried out using a high-density panel (BovineHD - Illumina® bead chip with 777 962 SNPs) for sires and the Illumina BovineSNP50 panel (54 609 SNPs) for their progenies. The relationship between the parents and progenies of genotyped animals was evaluated, and mismatches were based on the Mendelian conflicts counts. Variance components and genetic parameters estimates were obtained using a Bayesian approach via Gibbs sampling, and the breeding values were predicted assuming a repeatability model. A total of 460 corrections in relationship definitions were made (Table 1) corresponding to 1018 (9.5%) tick count records. Among these changes, 97.17% (447) were related to the sire's information, and 2.8% (13) were related to the dam's information. We observed 27.2% (236/868) of Mendelian conflicts for sire-progeny genotyped pairs and 14.3% (13/91) for dam-progeny genotyped pairs. We performed 2174 new definitions of half-siblings according to the correlation coefficient between the coancestry and molecular coancestry matrices. It was observed that higher-quality genetic relationships did not result in significant differences of variance components estimates; however, they resulted in more accurate breeding values predictions. Using SNPs to assess conflicts between parents and progenies increases certainty in relationships and consequently the accuracy of breeding value predictions of candidate animals for selection. Thus, higher genetic gains are expected when compared to the traditional non-corrected relationship matrix.

Genome-enabled prediction for tick resistance in Hereford and Braford beef cattle via reaction norm models.

Very few studies have been conducted to infer genotype × environment interaction (G×E) based in genomic prediction models using SNP markers. Therefore, our main objective was to compare a conventional genomic-based single-step model (HBLUP) with its reaction norm model extension (genomic 1-step linear reaction norm model [HLRNM]) to provide EBV for tick resistance as well as to compare predictive performance of these models with counterpart models that ignore SNP marker information, that is, a linear animal model (ABLUP) and its reaction norm extension (1-step linear reaction norm model [ALRNM]). Phenotypes included 10,673 tick counts on 4,363 Hereford and Braford animals, of which 3,591 were genotyped. Using the deviance information criterion for model choice, ABLUP and HBLUP seemed to be poorer fitting in comparison with their respective genomic model extensions. The HLRNM estimated lower average and reaction norm genetic variability compared with the ALRNM, whereas ABLUP and HBLUP seemed to be poorer fitting in comparison with their respective genomic reaction norm model extensions. Heritability and repeatability estimates varied along the environmental gradient (EG) and the genetic correlations were remarkably low between high and low EG, indicating the presence of G×E for tick resistance in these populations. Based on 5-fold -means partitioning, mean cross-validation estimates with their respective SE of predictive accuracy were 0.66 (SE 0.02), 0.67 (SE 0.02), 0.67 (SE 0.02), and 0.66 (SE 0.02) for ABLUP, HBLUP, HLRNM, and ALRNM, respectively. For 5-fold random partitioning, HLRNM (0.71 ± 0.01) was statistically different from ABLUP (0.67 ± 0.01). However, no statistical significance was reported when considering HBLUP (0.70 ± 0.01) and ALRNM (0.70 ± 0.01). Our results suggest that SNP marker information does not lead to higher prediction accuracies in reaction norm models. Furthermore, these accuracies decreased as the tick infestation level increased and as the relationship between animals in training and validation data sets decreased.

Multiple-breed reaction norm animal model accounting for robustness and heteroskedastic in a Nelore-Angus crossed population.

Our objective was to genetically characterize post-weaning weight gain (PWG), over a 345-day period after weaning, of Brangus-Ibagé (Nelore×Angus) cattle. Records (n=4016) were from the foundation herd of the Embrapa South Livestock Center. A Bayesian approach was used to assess genotype by environment (G×E) interaction and to identify a suitable model for the estimation of genetic parameters and use in genetic evaluation. A robust and heteroscedastic reaction norm multiple-breed animal model was proposed. The model accounted for heterogeneity of residual variance associated with effects of breed, heterozygosity, sex and contemporary group; and was robust with respect to outliers. Additive genetic effects were modeled for the intercept and slope of a reaction norm to changes in the environmental gradient. Inference was based on Monte Carlo Markov Chain of 110 000 cycles, after 10 000 cycles of burn-in. Bayesian model choice criteria indicated the proposed model was superior to simpler sub-models that did not account for G×E interaction, multiple-breed structure, robustness and heteroscedasticity. We conclude that, for the Brangus-Ibagé population, these factors should be jointly accounted for in genetic evaluation of PWG. Heritability estimates increased proportionally with improvement in the environmental conditions gradient. Therefore, an increased proportion of differences in performance among animals were explained by genetic factors rather than environmental factors as rearing conditions improved. As a consequence response to selection may be increased in favorable environments.

Genomic prediction for tick resistance in Braford and Hereford cattle.

One of the main animal health problems in tropical and subtropical cattle production is the bovine tick, which causes decreased performance, hide devaluation, increased production costs with acaricide treatments, and transmission of infectious diseases. This study investigated the utility of genomic prediction as a tool to select Braford (BO) and Hereford (HH) cattle resistant to ticks. The accuracy and bias of different methods for direct and blended genomic prediction was assessed using 10,673 tick counts obtained from 3,435 BO and 928 HH cattle belonging to the Delta G Connection breeding program. A subset of 2,803 BO and 652 HH samples were genotyped and 41,045 markers remained after quality control. Log transformed records were adjusted by a pedigree repeatability model to estimate variance components, genetic parameters, and breeding values (EBV) and subsequently used to obtain deregressed EBV. Estimated heritability and repeatability for tick counts were 0.19 ± 0.03 and 0.29 ± 0.01, respectively. Data were split into 5 subsets using k-means and random clustering for cross-validation of genomic predictions. Depending on the method, direct genomic value (DGV) prediction accuracies ranged from 0.35 with Bayes least absolute shrinkage and selection operator (LASSO) to 0.39 with BayesB for k-means clustering and between 0.42 with BayesLASSO and 0.45 with BayesC for random clustering. All genomic methods were superior to pedigree BLUP (PBLUP) accuracies of 0.26 for k-means and 0.29 for random groups, with highest accuracy gains obtained with BayesB (39%) for k-means and BayesC (55%) for random groups. Blending of historical phenotypic and pedigree information by different methods further increased DGV accuracies by values between 0.03 and 0.05 for direct prediction methods. However, highest accuracy was observed with single-step genomic BLUP with values of 0.48 for -means and 0.56, which represent, respectively, 84 and 93% improvement over PBLUP. Observed random clustering cross-validation breed-specific accuracies ranged between 0.29 and 0.36 for HH and between 0.55 and 0.61 for BO, depending on the blending method. These moderately high values for BO demonstrate that genomic predictions could be used as a practical tool to improve genetic resistance to ticks and in the development of resistant lines of this breed. For HH, accuracies are still in the low to moderate side and this breed training population needs to be increased before genomic selection could be reliably applied to improve tick resistance.

Response of lactating cows to live yeast supplementation during summer.

Dairy cows experiencing heat stress have reduced intake and increased reliance on glucose, making feeding strategies capable of improving diet digestibility plausible for improving postrumen nutrient flow and performance. The effect of yeast on digestion and performance of lactating cows during the warm summer months of southeastern Brazil was evaluated. Cows were individually fed in tie stalls and temperature-humidity index was above 68 during 75.6% of the experiment. Twenty-eight Holstein cows (207±87 d in milk) received a standard diet for 14 d and then a treatment for 70 d, in a covariate-adjusted, randomized block design with repeated measures over time. Treatments were yeast (Saccharomyces cerevisiae) or control. Yeast was top dressed to the diet in the morning, equivalent to 25×10(10) cfu of live cells and 5×10(10) cfu of dead cells. The diet contained corn silage (37.7%), Tifton silage (7.1%), raw soybeans (4.1%), soybean meal (16.5%), finely ground corn (20.7%), and citrus pulp (11.9%). Yeast increased milk (26.7 vs. 25.4 kg/d) and solids yield (3.06 vs. 2.92 kg/d), especially lactose. Response in milk yield was consistent over time and started at d 5. The daily intake of digestible OM, total-tract digestibility of nutrients, urinary allantoin excretion, chewing pattern throughout the day, and dry matter intake did not respond to yeast. A trend was observed for increased plasma glucose with yeast (62.9 vs. 57.3mg/dL), lowered respiratory frequency (48 vs. 56 breaths/min), and increased plasma niacin content (1.31 vs. 1.22 µg/mL), though cows had similar rectal temperature. Ruminal lactate and butyrate as proportions of ruminal organic acids were reduced by yeast, but no effects on other organic acids, ruminal pH, or protozoa content were detected. Plasma urea N over 24h was increased by yeast. On d 72 to 74, citrus pulp was abruptly replaced with finely ground corn to induce acidosis. The increased load of starch increased dry matter intake between 0700 and 1300 h, jugular blood partial pressure of CO2, HCO3-, and base excess, and decreased blood pH for both treatments. The yeast treatment had a higher blood pH compared with the control, 7.34, and 7.31, respectively. Yeast supplementation improved lactation performance of dairy cows under heat stress. Improvement in lactation performance apparently involved the regulation of body homeothermia, rather than improved digestibility.

Bovine genetic resistance effects on biological traits of Rhipicephalus (Boophilus) microplus.

This study aimed to verify the influence of bovine genetic resistance on biological traits of the Rhipicephalus (Boophilus) microplus tick. Genetic resistance or susceptibility was determined according to breeding values for tick counts, predicted using a dataset of 9007 Hereford and Braford (Hereford×Zebu) bovines naturally infested and raised under extensive production systems in southern Brazil. From a total of 974 Braford heifers born in 2008, 20 were classified as genetically tick-resistant and 20 classified as genetically tick-susceptible, and used to obtain the ticks samples used in this study. The 40 heifers were exposed to four subsequent artificial infestations with approximately 20,000 larvae at 14-day intervals. From the 19th to 23rd day of each infestation tick counts were performed on the left body side of the heifers. Engorged ticks were manually collected on the day of highest observed burden after each infestation. Tick counts on susceptible heifers were 5.5, 10.5, 11.1 and 6.9 times larger than on resistant heifers, respectively, after the first, second, third and fourth artificial infestations. In the third infestation, ticks from resistant heifers showed lower egg production index (P<0.0001) than ticks from susceptible heifers. In the fourth infestation, ticks from susceptible group showed higher egg mass weight (P<0.05) and nutrient index (P<0.0001) than ticks from resistant heifers. Tick initial weights showed a positive association with egg production index in susceptible heifers (P<0.05) and a negative association in the resistant group (P<0.05), suggesting a host defense mechanism that reduces the conversion efficiency of ingested blood to eggs in engorged ticks from resistant cattle. This shows that bovine genetic tick resistance, in addition to affecting the number of ticks carried by the animals, also affected the egg mass weight, egg production and nutrient indexes of ticks. The results of the present study imply that the selection of resistant animals could be used as a strategic tool for tick control in production systems, reducing infestation levels on cattle and environment.

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.

Milk yield and composition from Angus and Angus-cross beef cows raised in southern Brazil.

This study assessed milk yield and composition of Angus and Angus-cross beef cows raised in southern Brazil. A total of 128 records were collected in 2 consecutive calving seasons from cows between 3 and 5 yr of age of 4 breed compositions: Angus (ANAN), Caracu × Angus (CRAN), Hereford × Angus (HHAN), and Nelore × Angus (NEAN). These cows were mated to Brangus (BN) or Braford (BO) bulls and managed under extensive grazing conditions in southern Brazil. Milk production of these cows was assessed by 2 procedures: indirectly by the calf weigh-suckle-weigh procedure (WD) and directly by machine milking (MM). Lactation curves were estimated using nonlinear regression and the following related traits were derived: peak yield (PY), peak week (PW), total yield at 210 d (TY210), and lactation persistence (PERS). Milk composition and calf weaning weight adjusted to 210 d (WW210) were also determined. The MM technique was considered more accurate because of lower standard errors of estimated means, greater statistical power, and greater correlation between TY210 and WW210 (0.50) compared to WD (0.36). Considering the more precise evaluation by MM, the CRAN and NEAN cows had greater TY210 (1070 and 1116 kg, respectively) and PY (8.1 and 7.8 kg, respectively) compared to ANAN and HHAN cows, which had 858 and 842 kg for TY210 and 6.6 and 6.3 kg for PY, respectively. The NEAN cows had the latest PW at 10.8 wk. Late-calving cows had 21% lower TY210 compared to cows that calved earlier. Milk composition was influenced by cow genotype, with CRAN and NEAN cows producing milk with greater fat (3.8 and 3.9%, respectively) and protein (3.2 and 3.1%, respectively) content compared to ANAN and HHAN cows. Regardless of the genotype, fat, protein, and total solids increased in concentration from beginning to end of lactation, while lactose content decreased. Crossbreeding of Angus with adapted breeds of taurine or indicine origin can be effective in increasing milk yield and nutrient content and, consequently, producing heavier calves at weaning under extensive grazing in southern Brazil and other similar subtropical climate regions.

Bioeconomic model and selection indices in Aberdeen Angus cattle.

A bioeconomic model was developed to calculate economic values for biological traits in full-cycle production systems and propose selection indices based on selection criteria used in the Brazilian Aberdeen Angus genetic breeding programme (PROMEBO). To assess the impact of changes in the performance of the traits on the profit of the production system, the initial values ​​of the traits were increased by 1%. The economic values for number of calves weaned (NCW) and slaughter weight (SW) were, respectively, R$ 6.65 and R$ 1.43/cow/year. The selection index at weaning showed a 44.77% emphasis on body weight, 14.24% for conformation, 30.36% for early maturing and 10.63% for muscle development. The eighteen-month index showed emphasis of 77.61% for body weight, 4.99% for conformation, 11.09% for early maturing, 6.10% for muscle development and 0.22% for scrotal circumference. NCW showed highest economic impact, and SW had important positive effect on the economics of the production system. The selection index proposed can be used by breeders and should contribute to greater profitability.