Analysis of longitudinal data of beef cattle raised on pasture from northern Brazil using nonlinear models

This study was undertaken to aim of estimating the genetic parameters and trends for asymptotic weight (A) and maturity rate (k) of Nellore cattle from northern Brazil. The data set was made available by the Brazilian Association of Zebu Breeders and collected between the years of 1997 and 2007. The Von Bertalanffy, Brody, Gompertz, and logistic nonlinear models were fitted by the Gauss?CNewton method to weight?Cage data of 45,895 animals collected quarterly of the birth to 750?days old. The curve parameters were analyzed using the procedures GLM and CORR. The estimation of (co)variance components and genetic parameters was obtained using the MTDFREML software. The estimated heritability coefficients were 0.21?±?0.013 and 0.25?±?0.014 for asymptotic weight and maturity rate, respectively. This indicates that selection for any trait shall results in genetic progress in the herd. The genetic correlation between A and k was negative (?0.57?±?0.03) and indicated that animals selected for high maturity rate shall result in low asymptotic weight. The Von Bertalanffy function is adequate to establish the mean growth patterns and to predict the adult weight of Nellore cattle. This model is more accurate in predicting the birth weight of these animals and has better overall fit. The prediction of adult weight using nonlinear functions can be accurate when growth curve parameters and their (co)variance components are estimated jointly. The model used in this study can be applied to the prediction of mature weight in herds where a portion of the animals are culled before they reach the adult age.     

Optimal beef cattle diets formulated by nonlinear programming

A beef diet model based on National Research Council recommendations is significantly nonlinear for feed ingredients, daily gain and weight of cattle. Solving a diet model has been difficult, but advances in nonlinear programming now allow solutions that are quick and easy. This study developed a nonlinear programming method for optimally planning a feeding program by choosing feeds, daily gains and selling weight. Two types of diets are important for this purpose:optimal-return diets and least-cost-gain diets. For a given weight of cattle, an optimal-return diet chooses feeds and daily gain to maximize returns above feed costs. A least-cost-gain diet chooses feeds and daily gain to minimize feed plus yardage costs per kilogram of gain. In an optimal feeding program, a sequence of optimal-return diets is fed to increasing weights of cattle. Feed costs plus yardage per kilogram of gain rise to equal the actual selling price at the optimal selling weight, and the cattle are sold. Cattle feeders and researchers with access to a microcomputer can maximize net returns from a feeding program.     

从肉牛到牛肉的经济指标分析

【目的】为了解现阶段肉牛到牛肉之间的联系现状和产出能力,掌握主要经济指标,确定计算肉牛产值的依据。 【方法】以张掖市肉牛与牛肉之间的联系现状为例,从中选取１个能够代表现阶段生产水平的随机样本（n=11777）,利用统计分析方法进行分类和分组,列出次数分布表与分布图,研究宰前重、胴体重、净肉重、屠宰率和净肉率5个变量的频率分布,并根据平均数() 和标准差（S）及其分布规律推算出总体平均数u的概率置信区间代表群体生产水平,作为计算肉牛产值的依据。【结果】现阶段进行肉牛生产的主要品种为西门塔尔、安格斯和荷斯坦,公牛育肥产肉占84.54%。其中西门塔尔占74.61%居主导地位,安格斯占16.08%,荷斯坦占9.31%。西门塔尔牛以宰前重在500~700kg的小架子牛育肥产肉占56.12%居主导地位,依次为大架子育肥占39.28%、小牛肉占2.57%和成年牛占2.03%。肉牛宰前重、胴体重、净肉重和屠宰率、净肉率5个指标样本分组频率与正态分布理论概率相关度高、绝对差值小,概率相近。 【结论】现阶段生产牛肉的主要品种有西门塔尔、安格斯和荷斯坦,以公牛育肥为主。主导品种西门塔尔牛生产周期有缩短趋势,小牛肉生产起步。 西门塔尔牛群的宰前重、胴体重、净肉重和屠宰率、净肉率均值为：684.94kg、385.16kg、316.85kg和56.47%、45.71%时,区间为：581.33≦x≦788.55、321.9≦x≦448.40、266.09≦x≦367.61、53.85≦x≦59.09、42.57≦x≦48.85。 大群体宰前重、胴体重和屠宰率的均值为：678.13kg、379.54kg和55.81%%,区间为572.63≦x≦783.63 、311.42≦x≦447.66、52.34≦x≦59.28 。 反映肉牛生产水平的5个经济指标的概率分布为正态分布,平均值（）一倍标准差左右取值概率为0.6827,可用概率密度函数公式进行相应的计算。根据正态分布特点,平均值左右一倍标准差范围内取值的概率最大为,平均值即为群体盈亏平衡点,且主导品种的宰前重、胴体重和屠宰率较大群体分别高6.81kg、5.62kg和0.66%。 祁连牧歌实业有限公司为张掖市规模最大肉牛屠宰加工企业,屠宰牛源为周边规模养殖场育肥牛,该结果能够代表张掖市域内肉牛大群体生产水平,指标平均值（）可以作为计算肉牛产值的统计依据。 生产中可以此为标准对经营情况进行分析,当个体牛只或经营单位的宰前重、胴体重、净肉重和屠宰率、净肉率高于平均值时说明该牛只育肥状况较好或该经营单位经营状态良好,大概率保持盈利。     

HACCP模式在肉牛安全生产中的应用

HACCP是迄今为止食品安全生产方面最有效、最科学的质量管理方法之一。本文根据HACCP的技术规程,对肉牛生产的各流程进行分析,确定了肉牛生产过程中的关键控制点,并提出了控制措施。     

一例肉牛违法饲喂瘦肉精的案例分析

非法生产、销售、使用盐酸克伦特罗（瘦肉精）,严重危害广大人民群众的身体健康。严厉打击在饲料和饮用水中添加或以其他方式饲喂和使用瘦肉精,一直以来是各级政府和畜牧兽医部门的工作重点。今年4月份,岫岩满族自治县动物卫生监督管理局开展瘦肉精专项整治活动,在对一牛场的尿样检测中,有2头牛的尿样经检测为阳性,经过大量的工作,此案件依法予以处罚。     

Sire X environment interactions in beef cattle weaning weight field data

Weaning weight field records, supplied by the American Polled Hereford Association, were used to examine sire X environment interactions. Sire X herd/region and sire X contemporary group/herd interactions were evaluated from a data set containing 19,503 records. Sire X region interaction was evaluated from a data set containing 8,659 records. The genetic correlations of sire progeny performance across contemporary groups/herd were .59 and .37 across herds and contemporary groups/region. The average genetic correlation of sire progeny performance across regions was .64. Heritability of weaning weight was .11 across regions, .17 within region and .28 within herd. Mixed-model sire analyses of Polled Hereford weaning weight field records should include sire X herd/region and sire X contemporary group/herd random effects to reduce the sire X environment effects particular to any herd or contemporary group, and to account for the distribution of sire progeny across herds and contemporary groups in the estimation of prediction error variance. It may be necessary to perform separate sire analyses for some regions to evaluate the breeding values of sires in regions where rank changes are likely to occur.     

全株甘蔗饲养肉牛技术模式

甘蔗是一种以生产蔗糖为目的的经济作物 ,在我省热区的种植业中发挥着极其重要的作用 ,但在糖价波动较大的情况下 ,由于甘蔗的利用途径单一 ,蔗农的经济收入因制糖业的影响而起伏。通过研究 ,发现糖用甘蔗饲养肉牛具有一定的应用前景。直接利用全株甘蔗饲喂黄牛 10 0d而不补饲精料的试验结果显示 ,牛群没有出现病理反应 ;在饲喂全株甘蔗的基础上 ,在市场牛价不低于 6 0 0元 /kg活重时每天补饲0 5kg浓缩饲料 ,云南黄牛日增重达 0 5 5kg,可给蔗农带来较好的经济收益。这些研究结果为全株糖用甘蔗的饲用技术开发提供了依据。1　甘蔗的饲用特…     

Economics of deworming beef cattle

From the practitioner level and for the producer, there is far too much confusion concerning internal parasites in cattle and the justification of control measures. By coordinating research efforts between parasitologists and other scientists, awareness of internal parasites will be raised and the total impact of parasitism can better be evaluated.     

Development of a model to predict dietary metabolizable energy from digestible energy in beef cattle

Understanding the utilization of feed energy is essential for precision feeding in beef cattle production. We aimed to assess whether predicting the metabolizable energy (ME) to digestible energy (DE) ratio (MDR), rather than a prediction of ME with DE, is feasible and to develop a model equation to predict MDR in beef cattle. We constructed a literature database based on published data. A meta-analysis was conducted with 306 means from 69 studies containing both dietary DE and ME concentrations measured by calorimetry to test whether exclusion of the y-intercept is adequate in the linear relationship between DE and ME. A random coefficient model with study as the random variable was used to develop equations to predict MDR in growing and finishing beef cattle. Routinely measured or calculated variables in the field (body weight, age, daily gain, intake, and dietary nutrient components) were chosen as explanatory variables. The developed equations were evaluated with other published equations. The no-intercept linear equation was found to represent the relationship between DE and ME more appropriately than the equation with a y-intercept. The y-intercept (−0.025 ± 0.0525) was not different from 0 (P = 0.638), and Akaike and Bayesian information criteria of the no-intercept model were smaller than those with the y-intercept. Within our growing and finishing cattle data, the animal’s physiological stage was not a significant variable affecting MDR after accounting for the study effect (P = 0.213). The mean (±SE) of MDR was 0.849 (±0.0063). The best equation for predicting MDR (n = 106 from 28 studies) was 0.9410 ( ± 0.02160) +0.0042 ( ± 0.00186) × DMI (kg) – 0.0017 ( ± 0.00024) × NDF(% DM) – 0.0022 ( ± 0.00084) × CP(% DM). We also presented a model with a positive coefficient for the ether extract (n = 80 from 22 studies). When using these equations, the observed ME was predicted with high precision (R² = 0.92). The model accuracy was also high, as shown by the high concordance correlation coefficient (>0.95) and small root mean square error of prediction (RMSEP), <5% of the observed mean. Moreover, a significant portion of the RMSEP was due to random bias (> 93%), without mean or slope bias (P > 0.05). We concluded that dietary ME in beef cattle could be accurately estimated from dietary DE and its conversion factor, MDR, predicted by the dry matter intake and concentration of several dietary nutrients, using the 2 equations developed in this study.     

A practical longitudinal model for evaluating growth in Gelbvieh cattle

Genetic evaluation of growth in Gelbvieh beef cattle was examined by multiple-trait (MTM) and random regression (RRM) analysis. The data set comprised 541,108 animals with 1,120,086 records. Approximately 15% of the animals in the data set had at least one record measured outside of the accepted MTM age ranges for weaning weight (Wwt) and yearling weight (Ywt). Fourteen percent of Wwt records and 19% of Ywt records were measured outside the accepted ranges for MTM analysis, and thus were excluded from MTM evaluations. Two RRM evaluations were performed using cubic Legendre polynomials (RRML) and linear splines (RRMS) with three knots at 1, 205, and 365 d of age. Data Set 1 (d1) utilized all available records, whereas Data Set 2 (d2) included only records measured within MTM ranges (1 d, 160 to 250 d, and 320 to 410 d). The RRML models did not reach convergence until diagonalization was imposed. After diagonalization, it was found that all longitudinal models required fewer iterations to converge than the MTM. Correlations between the MTM, RRML-d2, and RRMS-d2 evaluations were >or=0.99 for all three traits, indicating that these models were equivalent when predicting breeding values from data within the MTM age ranges. Correlations between MTM, RRML-d1, and RRMS-d1 were >0.99 for Bwt and >0.95 for Wwt and Ywt. The lower correlations for Wwt and Ywt indicate that the added information does affect breeding value prediction. The RRM has the capability to incorporate records measured at all ages into genetic evaluations at a computing cost similar to the MTM.     

Evaluation of cytoplasmic genetic effects in beef cattle using an animal model

Mixed model techniques were used to evaluate the importance of cytoplasmic genetic effects on beef cattle performance. Birth weight (BWT), preweaning average daily gain (ADG), weaning weight (WT205), postweaning gain (PG), ultrasonic backfat thickness (FAT) and predicted milk yield (MILK) data were collected in two herds of Hereford cattle located at Plymouth and Raleigh, North Carolina. Cytoplasmic lines were determined based on the foundation female in the maternal lineage of each animal. An animal model was used to account for all nuclear genetic variation among animals within herds. Direct breeding values were estimated for all animals with records and their parents for all traits. For MILK, permanent environmental effects were estimated for animals with multiple records. For preweaning traits, maternal breeding values and permanent maternal environmental effects also were estimated. In all analyses, F-tests for cytoplasmic effects were not significant. Probability values approached significance (P = .15 to P = .10) only for PG and FAT at Plymouth. Assumptions regarding the ratios of genetic and environmental variances and covariances had no effect on F-tests. Results contrast with earlier analyses of the same data in which nuclear genetic effects were accounted for by including sire and maternal grandsire in the statistical model. This study failed to show that cytoplasmic genetic effects were important sources of variation in performance; residual additive genetic effects were confounded with cytoplasmic lines for these herds. Because cytoplasmic sources may be regarded as founder effects, further research is needed in other populations.     

Utilization of beef from different cattle phenotypes to produce a guaranteed tender beef product

Cattle (n = 303) were visually selected from four feed yards to represent six phenotypes (English [EN; n = 50], 3/4 English-1/4 Brahman [ENB; n = 52], 1/2 English-1/2 Exotic [ENEX; n = 56], 1/2 English-1/4 Exotic-1/4 Brahman [ENEXB; n = 47], 3/4 Exotic-1/4 Brahman [EXB; n = 49], and 1/2 Exotic-1/4 English-1/4 Brahman [EXENB; n = 49]). Carcasses were processed at a commercial beef packing facility, and strip loins were collected after 48-h chilling. Strip loins were aged for 14 d at 2 degrees C and frozen at -20 degrees C for 3 to 5 d before three 2.5-cm-thick steaks were cut for Warner-Bratzler shear force (WBSF) determinations and sensory evaluations. Phenotype EN had the highest (P < 0.05) adjusted fat thickness, and EXB had adjusted fat thickness that was lower (P < 0.05) than all other phenotypes except EXENB. Carcasses of EN and ENB had smaller (P < 0.05) longissimus muscle areas than phenotypes ENEX, EXB, and EXENB. Phenotype EN produced carcasses with the highest (P < 0.05) numerical yield grade, whereas carcasses originating from phenotype EXB had lower (P < 0.05) numerical yield grades than all other phenotypes except ENEX. No differences (P > 0.05) were found among phenotypes for mean WBSF values or sensory panel ratings for initial and sustained tenderness, initial and sustained juiciness, beef flavor characteristics, and overall mouthfeel. More than 90% of steaks from carcasses of all phenotypes had WBSF values less than 3.6 kg when cooked to an internal cooked temperature of 70 degrees C. Results from this study indicated that all phenotypes represented in this study could be managed to produce tender beef.     

Joint longitudinal modeling of age of dam and age of animal for growth traits in beef cattle

Two methods to jointly model age of dam (AOD) and age of animal in random regression analyses of growth in Gelbvieh cattle were examined. The first method (M1) was analogous to the multiple-trait analysis and consisted of AOD as a nested class variable and a cubic polynomial regression on age nested within birth, weaning, and yearly weights. The second method (M2) used two-dimensional splines, with age knots at 150, 205, 270, 340, and 390 d. The AOD knots were placed at 725, 1,464, and 2,189 d. These selected knots were used to form a two-dimensional grid containing 15 knots, each representing a specific age and AOD combination. A data set containing Gelbvieh growth records was split along contemporary groups into two data sets. Data set 1 contained 316,078 records and was used for prediction by mixed-model equations. Data set 2 contained 164,167 records and was used for cross validation. In the complete data set, only 90 and 30% of animals with birth weight had records on weaning and yearling weights, respectively. Models were evaluated based on R2, average squared error (ASE), percent bias, and plots of solutions. The ASE for weights associated with birth weight, weaning weight, and yearling weight for M1 were 15, 505, and 703 kg2. With M2, large jumps in fixed-effect estimates were observed outside the two-dimensional grid. To eliminate this problem, weighted one-dimensional splines were used for extrapolation beyond the two-dimensional grid. For M2 with weighted spline extrapolation, the ASE were 15, 542, and 777 kg2 for birth weight, weaning weight, and yearling weight, respectively. Creation of optimal two-dimensional splines is difficult when data are clustered. Despite such difficulties, the two-dimensional spline was capable of jointly and continuously modeling AOD and age of animal.     

Two methods for parameter estimation using multiple-trait models and beef cattle field data

Two methods are presented for estimating variances and covariances from beef cattle field data using multiple-trait sire models. Both methods require that the first trait have no missing records and that the contemporary groups for the second trait be subsets of the contemporary groups for the first trait; however, the second trait may have missing records. One method uses pseudo expectations involving quadratics composed of the solutions and the right-hand sides of the mixed model equations. The other method is an extension of Henderson's Simple Method to the multiple trait case. Neither of these methods requires any inversions of large matrices in the computation of the parameters; therefore, both methods can handle very large sets of data. Four simulated data sets were generated to evaluate the methods. In general, both methods estimated genetic correlations and heritabilities that were close to the Restricted Maximum Likelihood estimates and the true data set values, even when selection within contemporary groups was practiced. The estimates of residual correlations by both methods, however, were biased by selection. These two methods can be useful in estimating variances and covariances from multiple-trait models in large populations that have undergone a minimal amount of selection within contemporary groups.     

Illustration of the maternal animal model used for genetic evaluation of beef cattle

National cattle evaluation programs for weaning weight in most beef breed associations involve implementation of the maternal animal model to predict direct and maternal EPD. With this model, direct breeding values are predicted for all animals with records or pedigree ties to animals with records, or both. Even though maternal genetic value is expressed only in animals that become dams, these effects are transmitted by all parents and inherited from parents by all animals, leading to maternal breeding values being predicted for all animals as well. A small example data set was simulated involving 12 parents, 8 nonparents, and 13 animals with weaning weight records. The pedigree was developed to include paternal and maternal half-sib families, full-sibs, and some inbreeding, similar to field populations of beef cattle. Assembly of the mixed model equations and solutions for the maternal animal model are illustrated explicitly to assist animal breeding students in their understanding of the properties of the maternal animal model and to explicitly implement the model. Model parameters and moments, fixed contemporary group solutions, adjustment of breeding values for merit of mates, interpretation of maternal permanent environmental effect solutions, and alternatives for the assembly of the equations are shown. This example should lead to increased student and producer understanding of genetic improvement programs for weaning weight in beef cattle.     

澳大利亚2019年肉牛产业发展分析与借鉴

澳大利亚是世界上牛肉出口大国之一,肉牛产业是其重要的支柱产业。我国是澳大利亚重要的牛肉出口国之一,随着我国与澳大利亚自由贸易协议的签订及实施,近几年我国对澳大利亚的牛肉进口量持续快速攀升,2019年跃居澳大利亚牛肉出口国第一位。本文收集汇总了澳大利亚官方发布的相关数据,对2019年澳大利亚肉牛产业发展现状、贸易情况进行追踪,通过分析以期为我国肉牛产业持续健康发展提出借鉴。     

Leptin as a predictor of carcass composition in beef cattle

Our objective was to determine if serum concentrations of leptin could be used to predict carcass composition and merit in feedlot finished cattle. Two different groups of crossbred Bos taurus steers and heifers were managed under feedlot conditions near Miles City, MT. The first group consisted of 88 1/2 Red Angus, 1/4 Charolais, and 1/4 Tarentaise composite gene combination steers (CGC) harvested at the ConAgra processing facility in Greeley, CO. The second group (Lean Beef Project; LB) consisted of 91 F2 steers and heifers born to Limousin, Hereford, or Piedmontese by CGC F1 cows crossed to F1 bulls of similar breed composition and harvested at a local processing facility in Miles City, MT. Blood samples were collected approximately 24 h before harvest (CGC) or approximately 3 d before and at harvest (LB). No differences in serum concentrations of leptin were detected (P > 0.10) between Hereford, Limousin, or Piedmontese F2 calves nor between LB steers and heifers. Positive correlations (P < 0.01) existed between serum leptin and marbling score (r = 0.35 and 0.50), fat depth measured between the 12th and 13th rib (r = 0.34 and 0.46), kidney, pelvic, and heart fat (KPH) (r = 0.42 and 0.46), and quality grade (r = 0.36 and 0.49) in CGC and LB cattle, respectively. Serum leptin was also positively correlated with calculated yield grade for CGC steers (r = 0. 19; P = 0. 10) and LB cattle (r = 0.52; P < 0.01). Longissimus area was not correlated with serum leptin in CGC steers (r = 0.12; P > 0.10). However, a negative correlation existed between longissimus area and serum leptin in the LB cattle (r = -0.45; P < 0.01). Serum concentrations of leptin were significantly associated with carcass composition (marbling, back fat depth, and KPH fat) and quality grade in both groups of cattle studied and may provide an additional indicator of fat content in feedlot cattle.     

Sulfur-induced polioencephalomalacia in a herd of rotationally grazed beef cattle

Central nervous system disease occurred in a herd of rotationally grazed beef cattle consuming water containing 3400 ppm sulfate. Clinical signs, pathologic findings, and high water sulfate levels confirmed the diagnosis of sulfur-induced polioencephalomalacia. The incidence of disease reduced when the herd was switched to a low sulfate water source.     

Estimates of covariance functions for growth of Nelore cattle applying a parametric correlation structure to model within-animal correlations

A total of 20,065 weights recorded on 3016 Nelore animals were used to estimate covariance functions for growth from birth to 630 days of age, assuming a parametric correlation structure to model within-animal correlations. The model of analysis included fixed effects of contemporary groups and age of dam as quadratic covariable. Mean trends were taken into account by a cubic regression on orthogonal polynomials of animal age. Genetic effects of the animal and its dam and maternal permanent environmental effects were modelled by random regressions on Legendre polynomials of age at recording. Changes in direct permanent environmental effect variances were modelled by a polynomial variance function, together with a parametric correlation function to account for correlations between ages. Stationary and nonstationary models were used to model within-animal correlations between different ages. Residual variances were considered homogeneous or heterogeneous, with changes modelled by a step or polynomial function of age at recording. Based on Bayesian information criterion, a model with a cubic variance function combined with a nonstationary correlation function for permanent environmental effects, with 49 parameters to be estimated, fitted best. Modelling within-animal correlations through a parametric correlation structure can describe the variation pattern adequately. Moreover, the number of parameters to be estimated can be decreased substantially compared to a model fitting random regression on Legendre polynomial of age.