A survey of polymorphisms detected from sequences of popular beef breeds.

The genome sequence was obtained from 270 sires used in the Germplasm Evaluation (GPE) project. These bulls included 154 purebred AI sires from GPE Cycle VII breeds (Hereford, Angus, Simmental, Limousin, Charolais, Gelbvieh, and Red Angus), 83 F crosses of those breeds, and 33 AI sires from 8 other breeds. The exome capture sequence targeting coding regions of the genome was obtained from 176 of these bulls. Sequence reads were mapped to the UMD 3.1 bovine genome assembly; a mean of 2.5-fold (x) coverage per bull was obtained from the genomic sequence, and the targeted exons were covered at a mean of 20.0x. Over 28.8 million biallelic sequence variants were detected where each allele was present in at least 3 different bulls. These included 22.0 million previously reported variants and 94.1% of the 774,660 autosomal and BTA X SNP on the BovineHD BeadChip assay (HD). More than 92% of the variants detected in targeted exons were also detected from the low-coverage genome sequence. Less than 1% of the variants detected from the combined genome and exome sequence occurred in annotated protein-coding sequences and 5' and 3' untranslated regions (UTR) surrounding the 19,994 annotated protein coding regions. Variation was detected in the coding sequence or UTR of 96.8% of the genes: loss-of-function variants were predicted for 3,298 genes, 14,973 contained nonsynonymous variants, 11,276 had variation in UTR, and 17,721 genes contained synonymous variants. Minor allele frequencies (MAF) were <0.05 for 47.8% of the coding sequence and UTR variants, and MAF distributions were skewed toward low MAF. In contrast, 11.1% of the HD SNP detected in these bulls had MAF < 0.05, and the distribution was skewed toward higher MAF. Genes involved in immune system processes and immune response were overrepresented among those genes containing high MAF loss-of-function and nonsynonymous polymorphisms. Detected variants were submitted to the National Center for Biotechnology Information genetic variation database (dbSNP) under the handle MARC, batch GPE_Bull_GenEx.

Genomewide association study of reproductive efficiency in female cattle.

Reproductive efficiency is of economic importance in commercial beef cattle production, as failure to achieve pregnancy reduces the number of calves marketed per cow exposed. Identification of genetic markers with predictive merit for reproductive success would facilitate early selection of sires with daughters having improved reproductive rate without increasing generation intervals. To identify regions of the genome harboring variation affecting reproductive success, we applied a genomewide association study (GWAS) approach based on the >700,000 SNP marker assay, using a procedure based on genotyping multianimal pools of DNA to increase the number of animals that could be genotyped with available resources. Cows from several populations were classified according to reproductive efficiency, and DNA was pooled within population and phenotype prior to genotyping. Populations evaluated included a research population at the U.S. Meat Animal Research Center, 2 large commercial ranch populations, and a number of smaller populations (<100 head) across the United States. We detected 2 SNP with significant genomewide association (P ≤ 1.49 × 10(-7)), on BTA21 and BTA29, 3 SNP with suggestive associations (P ≤ 2.91 × 10(-6)) on BTA5, and 1 SNP with suggestive association each on BTA1 and BTA25. In addition to our novel findings, we confirmed previously published associations for SNP on BTA-X and all autosomes except 3 (BTA21, BTA22, and BTA28) encompassing substantial breed diversity including Bos indicus and Bos taurus breeds. The study identified regions of the genome associated with reproductive efficiency, which are being targeted for further analysis to develop robust marker systems, and demonstrated that DNA pooling can be used to substantially reduce the cost of GWAS in cattle.

Use of robust multivariate linear mixed models for estimation of genetic parameters for carcass traits in beef cattle.

Assumptions of normality of residuals for carcass evaluation may make inferences vulnerable to the presence of outliers, but heavy-tail densities are viable alternatives to normal distributions and provide robustness against unusual or outlying observations when used to model the densities of residual effects. We compare estimates of genetic parameters by fitting multivariate Normal (MN) or heavy-tail distributions (multivariate Student's t and multivariate Slash, MSt and MS) for residuals in data of hot carcass weight (HCW), longissimus muscle area (REA) and 12th to 13th rib fat (FAT) traits in beef cattle using 2475 records from 2007 to 2008 from a large commercial operation in Nebraska. Model comparisons using deviance information criteria (DIC) favoured MSt over MS and MN models, respectively. The posterior means (and 95% posterior probability intervals, PPI) of v for the MSt and MS models were 5.89 ± 0.90 (4.35, 7.86) and 2.04 ± 0.18 (1.70, 2.41), respectively. Smaller values of posterior densities of v for MSt and MS models confirm that the assumption of normally distributed residuals is not adequate for the analysis of the data set. Posterior mean (PM) and posterior median (PD) estimates of direct genetic variances were variable with MSt having the highest mean value followed by MS and MN, respectively. Posterior inferences on genetic variance were, however, comparable among the models for FAT. Posterior inference on additive heritabilities for HCW, REA and FAT using MN, MSt and MS models indicated similar and moderate heritability comparable with the literature. Posterior means of genetic correlations for carcass traits were variable but positive except for between REA and FAT, which showed an antagonistic relationship. We have demonstrated that genetic evaluation and selection strategies will be sensitive to the assumed model for residuals.

Deletion on chromosome 5 associated with decreased reproductive efficiency in female cattle.

Reproductive efficiency is arguably the most economically important trait in commercial beef cattle production, as failure to achieve pregnancy reduces the number of calves marketed per cow exposed to breeding. Identification of variation in the genome with predictive merit for reproductive success would facilitate accurate prediction of daughter pregnancy rate in sires enabling effective selection of bulls producing daughters with improved fertility. Copy number variation (CNV) is one form of structural variation that results in abnormal copies of DNA in the genome. Our lab previously reported a region between 25 and 70 Mb on chromosome 5 associated with reproductive efficiency. To further evaluate this region for genomic variations such as CNV, we initially applied a genomewide association approach based on the >700,000 SNP marker assay and PennCNV analysis to 68 Bos indicus cross females from a large commercial ranch in central Florida. A genomic segment located on chromosome 5, spanning the region of approximately 58.37 to 58.44 Mb (Bovine UMD3.1 assembly) was identified as containing a deletion associated with decreased reproductive efficiency. To verify this deletion, a quantitative real-time PCR test was developed to evaluate additional females in the central Florida and U.S. Meat Animal Research Center (USMARC) populations. The frequency of the homozygous deletion was 7% (23/319) in the central Florida population (Bos indicus cross) for females with 2 consecutive failed breeding seasons and 0% in the USMARC population (predominantly Bos taurus) of low reproductive females. To date, we have not identified the deletion in Bos taurus cattle, suggesting that the deletion is specific to Bos indicus cattle. From these data, we have identified a deletion on chromosome 5 associated with reproductive efficiency in Bos indicus-influenced cattle.

A simulation-based approach for evaluating and comparing the environmental footprints of beef production systems.

A methodology was developed and used to determine environmental footprints of beef cattle produced at the U.S. Meat Animal Research Center (MARC) in Clay Center, NE, with the goal of quantifying improvements achieved over the past 40 yr. Information for MARC operations was gathered and used to establish parameters representing their production system with the Integrated Farm System Model. The MARC farm, cow-calf, and feedlot operations were each simulated over recent historical weather to evaluate performance, environmental impact, and economics. The current farm operation included 841 ha of alfalfa and 1,160 ha of corn to produce feed predominately for the beef herd of 5,500 cows, 1,180 replacement cattle, and 3,724 cattle finished per year. Spring and fall cow-calf herds were fed on 9,713 ha of pastureland supplemented through the winter with hay and silage produced by the farm operation. Feedlot cattle were backgrounded for 3 mo on hay and silage with some grain and finished over 7 mo on a diet high in corn and wet distillers grain. For weather year 2011, simulated feed production and use, energy use, and production costs were within 1% of actual records. A 25-yr simulation of their current production system gave an average annual carbon footprint of 10.9±0.6 kg of CO2 equivalent units per kg BW sold, and the energy required to produce that beef (energy footprint) was 26.5±4.5 MJ/kg BW. The annual water required (water footprint) was 21,300±5,600 L/kg BW sold, and the water footprint excluding precipitation was 2,790±910 L/kg BW. The simulated annual cost of producing their beef was US$2.11±0.05/kg BW. Simulation of the production practices of 2005 indicated that the inclusion of distillers grain in animal diets has had a relatively small effect on environmental footprints except that reactive nitrogen loss has increased 10%. Compared to 1970, the carbon footprint of the beef produced has decreased 6% with no change in the energy footprint, a 3% reduction in the reactive nitrogen footprint, and a 6% reduction in the real cost of production. The water footprint, excluding precipitation, has increased 42% due to greater use of irrigated corn production. This proven methodology provides a means for developing the production data needed to support regional and national full life cycle assessments of the sustainability of beef.

Application of multivariate heavy-tailed distributions to residuals in the estimation of genetic parameters of growth traits in beef cattle.

Assumptions of normality in most animal breeding applications may make inferences vulnerable to the presence of outliers. Heavy-tail densities are viable alternatives to normal distributions and provide robustness against unusual or outlying observations when used to model the densities of residual effects. Our objective is to compare estimates of genetic parameters by fitting multivariate normal (MN) or heavy-tail distributions [multivariate Student's t (MSt) and multivariate slash (MS)] for residuals in data of body birth weight (BBW), weaning (WW), and yearling (YW) weight traits in beef cattle. A total of 17,019 weight records for BBW, WW, and YW from 1998 through 2010 from a large commercial cow/calf operation in the sand hills of Nebraska were analyzed. Models included fixed effects of contemporary group and sire breed whereas animal and maternal effects were random and the degrees of freedom (v) was treated as unknown for MSt and MS. Model comparisons using deviance information criteria (DIC) favored MSt over MS and MN models, respectively. The posterior means [and 95% posterior probability intervals (PPI)] of v for the MSt and MS models were 5.28 (4.80, 5.85) and 1.88 (1.76, 2.00), respectively. Smaller values of posterior densities of v for MSt and MS models confirm that the assumption of normally distributed residuals is not adequate for the analysis of BBW, WW, and YW datasets. Posterior mean (PM) and posterior median (PD) estimates of direct and maternal genetic variances were the same and posterior densities of these parameters were found to be symmetric. The 95% PPI estimates from MN and MSt models for BBW did not overlap, which indicates significant difference between PM estimates from MN or MSt models. The observed antagonistic relationship between additive direct and additive maternal effects indicated that genetic evaluation and selection strategies will be sensitive to the assumed model for residuals.

The accuracies of DNA-based estimates of genetic merit derived from Angus or multibreed beef cattle training populations.

Several organizations have developed prediction models for molecular breeding values (MBV) for quantitative growth and carcass traits in beef cattle using Bovine SNP50 genotypes and phenotypic or EBV data. Molecular breeding values for Angus cattle have been developed by IGENITY, Pfizer Animal Genetics, and a collaboration between researchers from Iowa State University and the University of Missouri-Columbia (ISU/UMC). The U.S. Meat Animal Research Center (USMARC; Clay Center, NE) has also developed MBV for 16 cattle breeds using 2 multibreed populations, the Germplasm Evaluation (GPE) Program and the 2,000 Bull Project (2K(ALL)), and 2 single breed subpopulations of the 2,000 Bull Project, Angus (2K(AN)) and Hereford (2K(HH)). In this study, these MBV were assessed relative to commercial ranch EBV estimated from the progeny phenotypes of Angus bulls naturally mated in multisire breeding pastures to commercial cows: 121 for USMARC MBV, 99 for ISU/UMC MBV, and 29 for IGENITY and Pfizer MBV (selected based on number of progeny carcass records). Five traits were analyzed: weaning weight (WW), HCW, marbling score (MS), rib-eye muscle area (RE), and, for IGENITY and Pfizer only, feedlot ADG. The average accuracies of MBV across traits were 0.38 ± 0.05 for IGENITY, 0.61 ± 0.12 for Pfizer, 0.46 ± 0.12 for ISU/UMC, 0.16 ± 0.04 for GPE, 0.26 ± 0.05 for 2K(ALL), 0.24 ± 0.04 for 2K(AN), and 0.02 ± 0.12 for 2K(HH). Angus-based MBV (IGENITY, Pfizer, ISU/UMC, and 2K(AN)) explained larger proportions of genetic variance in this population than GPE, 2K(ALL), or 2K(HH) MBV for the same traits. In this data set, IGENITY, Pfizer, and ISU/UMC MBV were predictive of realized performance of progeny, and incorporation of that information into national genetic evaluations would be expected to improve EPD accuracy, particularly for young animals.

Y are you not pregnant: identification of Y chromosome segments in female cattle with decreased reproductive efficiency.

Reproductive efficiency is of economic importance in commercial beef cattle production, since failure to achieve pregnancy reduces the number of calves marketed. Identification of genetic markers with predictive merit for reproductive success would facilitate early selection of females and avoid inefficiencies associated with sub-fertile cows. To identify regions of the genome harboring variation affecting reproductive success, we applied a genome-wide association approach based on the >700,000 SNP marker assay. To include the largest number of individuals possible under the available budget, cows from several populations were assigned to extremes for reproductive efficiency, and DNA was pooled within population and phenotype before genotyping. Surprisingly, pools prepared from DNA of low reproductive cattle returned fluorescence intensity data intermediate between fertile females and males for SNP mapped to the Y chromosome (i.e., male sex chromosome). The presence of Y-associated material in low reproductive heifers or cows was confirmed by Y-directed PCR, which revealed that 21 to 29% of females in the low reproductive category were positive by a Y chromosome PCR test normally used to sex embryos. The presence of the Y chromosome anomaly was further confirmed with application of additional Y-specific PCR amplicons, indicating the likelihood of the presence of some portion of male sex chromosome in female cattle in various beef cattle herds across the U.S. Discovery of this Y anomaly in low reproductive females may make an important contribution to management of reproductive failures in beef cattle operations.

Genetic evaluation of retail product percentage in Simmental cattle.

The objective of this study was to estimate genetic parameters required for genetic evaluation of retail product percentage (RPP) in Simmental cattle. Carcass weight (HCW), subcutaneous fat thickness (FAT), longissimus muscle area (REA) and kidney, pelvic, and heart fat (KPH) records were available to compute RPP on steers (n = 5171) and heifers (n = 1400) from the American Simmental Association database; animals were sired by 561 Simmental bulls and out of 5886 crossbred dams. Genetic parameters were estimated using residual maximal likelihood and a four trait animal model for the components of RPP including fixed harvest contemporary group effects, random animal genetic effects, and a linear covariate for age at harvest. Heritability estimates were 0.51 +/- 0.05, 0.36 +/- 0.05, 0.46 +/- 0.05, and 0.18 +/- 0.05 for HCW, FAT, REA and KPH respectively. Non-zero genetic correlations were estimated between HCW and REA (rg = 0.51 +/- 0.06) and between REA and FAT (rg = -0.43 +/- 0.08), but other genetic correlation estimates among the component traits were low. As a linear function of its components, heritability and genetic correlations involving RPP were estimated using index methods. The heritability estimate for RPP was 0.41, and genetic correlations were -0.17, -0.83, 0.67, and 0.01 with HCW, FAT, REA and KPH respectively. Therefore, RPP was strongly associated with muscle and fat deposition, but essentially independent of carcass weight and internal body cavity fat. Genetic evaluation of RPP would be straightforward using multiple trait index methods and genetic regression, although the inclusion of KPH would be of marginal value.

DNA-based paternity analysis and genetic evaluation in a large, commercial cattle ranch setting.

Deoxyribonucleic acid-based tests were used to assign paternity to 625 calves from a multiple-sire breeding pasture. There was a large variability in calf output and a large proportion of young bulls that did not sire any offspring. Five of 27 herd sires produced over 50% of the calves, whereas 10 sires produced no progeny and 9 of these were yearling bulls. A comparison was made between the paternity results obtained when using a DNA marker panel with a high (0.999), cumulative parentage exclusion probability (P(E)) and those obtained when using a marker panel with a lower P(E) (0.956). A large percentage (67%) of the calves had multiple qualifying sires when using the lower resolution panel. Assignment of the most probable sire using a likelihood-based method based on genotypic information resolved this problem in approximately 80% of the cases, resulting in 75% agreement between the 2 marker panels. The correlation between weaning weight, on-farm EPD based on pedigrees inferred from the 2 marker panels was 0.94 for the 24 bulls that sired progeny. Partial progeny assignments inferred from the lower resolution panel resulted in the generation of EPD for bulls that actually sired no progeny according to the high-P(E) panel, although the Beef Improvement Federation accuracies of EPD for these bulls were never greater than 0.14. Simulations were performed to model the effect of loci number, minor allele frequency, and the number of offspring per bull on the accuracy of genetic evaluations based on parentage determinations derived from SNP marker panels. The SNP marker panels of 36 and 40 loci produced EPD with accuracies nearly identical to those EPD resulting from use of the true pedigree. However, in field situations where factors including variable calf output per sire, large sire cohorts, relatedness among sires, low minor allele frequencies, and missing data can occur concurrently, the use of marker panels with a larger number of SNP loci will be required to obtain accurate on-farm EPD.

Genetic evaluation of beef carcass data using different endpoint adjustments.

Carcass data from 6,795 Simmental-sired animals born from 1992 to 2001 were used to determine whether adjustment to a constant age, back-fat, HCW, or marbling score would result in differences in heritability of the carcass traits and, correspondingly, if EPD calculated using those variance components and adjustments would result in sire reranking. The endpoints were age (EPA), backfat (EPF), HCW (EPC), or marbling (EPM). The traits analyzed were 12th-rib backfat (FAT), HCW, marbling (MRB), LM area (LMA), and percentage retail cuts (PRC). The data were analyzed using an animal model, where contemporary group was included as a fixed effect and was composed of slaughter date, sex, and herd. Random effects included in the model were direct genetic and residual. Estimates of heritability ranged from 0.12 to 0.14, 0.32 to 0.34, and 0.26 to 0.27 for FAT, HCW, and LMA, respectively, for the corresponding endpoints. Heritability for MRB was estimated to be 0.27 at all endpoints. For PRC, estimates of heritability were more variable, with estimates of 0.23 +/- 0.05, 0.32 +/- 0.05, 0.21 +/- 0.05, and 0.20 +/- 0.04 for EPA, EPF, EPC, and EPM, respectively. However, because the EPF and EPC adjustments adjust for a component trait of PRC (FAT and HCW, respectively), they may be altering the trait to one different from PRC. Spearman rank correlations between EPD within a trait using EPA compared with the other endpoints were >0.90 (P < 0.01) for FAT, HCW, MRB, and LMA. For PRC, Spearman rank correlations with EPA EPD were 0.73 (P < 0.01), 0.93 (P < 0.01), and 0.95 (P < 0.01) for EPF, EPC, and EPM, respectively. For most traits and endpoints, there was little reranking among sires when alternative endpoints were used. However, adjusting PRC to EPF appears to result in a greater heritability and substantial re-ranking of sires, potentially due to the adjustment changing the trait to one other than PRC.

Phenotypic ranges and relationships among carcass and meat palatability traits for fourteen cattle breeds, and heritabilities and expected progeny differences for Warner-Bratzler shear force in three beef cattle breeds.

Carcass and Warner-Bratzler shear force (WBSF) data from strip loin steaks were obtained from 7,179 progeny of Angus, Brahman, Brangus, Charolais, Gelbvieh, Hereford, Limousin, Maine-Anjou, Red Angus, Salers, Shorthorn, Simbrah, Simmental, and South Devon sires. Trained sensory panel (TSP) evaluations were obtained on 2,320 steaks sampled from contemporary groups of progeny from one to five sires of each breed. Expected progeny differences for marbling and WBSF were developed for 103 Simmental sires from 1,295 progeny, 23 Shorthorn sires from 310 progeny, and 69 Hereford sires from 1,457 progeny. Pooled phenotypic residual correlations, including all progeny, showed that marbling was lowly correlated with WBSF (-0.21) and with TSP overall tenderness (0.18). The residual correlation between WBSF and TSP tenderness was -0.68, whereas residual correlations for progeny sired by the three Bos indicus breeds were only slightly different than for progeny sired by Bos taurus breeds. The phenotypic range of mean WBSF among sires across breeds was 6.27 kg, and the phenotypic range among breed means was 3.93 kg. Heritability estimates for fat thickness, marbling score, WBSF, and TSP tenderness, juiciness, and flavor were 0.19, 0.68, 0.40, 0.37, 0.46, and 0.07, respectively. Ranges in EPD for WBSF and marbling were -0.41 to +0.26 kg and +0.48 to -0.22, respectively, for Simmentals; -0.41 to +0.36 kg and 0.00 to -0.32, respectively, for Shorthorns; and -0.48 to +0.22 kg and +0.40 to -0.24, respectively, for Herefords. More than 20% of steaks were unacceptable in tenderness. Results of this study demonstrated that 1) selection for marbling would result in little improvement in meat tenderness; 2) heritability of marbling, tenderness, and juiciness are high; and 3) sufficient variation exists in WBSF EPD among widely used Simmental, Shorthorn, and Hereford sires to allow for genetic improvement in LM tenderness.

Evaluation of Simmental carcass EPD estimated using live and carcass data.

This study was conducted to compare carcass EPD predicted using yearling live animal data and/or progeny carcass data, and to quantify the association between the carcass phenotype of progeny and the sire EPD. The live data model (L) included scan weight, ultrasound fat thickness, longissimus muscle area, and percentage of intramuscular fat from yearling (369 d of age) Simmental bulls and heifers. The carcass data model (C) included hot carcass weight, fat thickness, longissimus muscle area, and marbling score from Simmental-sired steers and cull heifers (453 d of age). The combined data model (F) included live animal and carcass data as separate but correlated traits. All data and pedigree information on 39,566 animals were obtained from the American Simmental Association, and all EPD were predicted using animal model procedures. The genetic model included fixed effects of contemporary group and a linear covariate for age at measurement, and a random animal genetic effect. The EPD from L had smaller variance and range than those from either C or F. Further, EPD from F had highest average accuracy. Correlations indicated that evaluations from C and F were most similar, and L would significantly (P < 0.05) re-rank sires compared with models including carcass data. Progeny (n = 824) with carcass data collected subsequent to evaluation were used to quantify the association between progeny phenotype and sire EPD using a model including contemporary group, and linear regressions for age at slaughter and the appropriate sire EPD. The regression coefficient was generally improved for sire EPD from L when genetic regression was used to scale EPD to the appropriate carcass trait basis. The EPD from C and F had similar linear associations with progeny phenotype, although EPD from F may be considered optimal because of increased accuracy. These data suggest that carcass EPD based on a combination of live and carcass data predict differences in progeny phenotype at or near theoretical expectation.

Genetic parameters for carcass traits and their live animal indicators in Simmental cattle.

The objective of this study was to estimate parameters required for genetic evaluation of Simmental carcass merit using carcass and live animal data. Carcass weight, fat thickness, longissimus muscle area, and marbling score were available from 5,750 steers and 1,504 heifers sired by Simmental bulls. Additionally, yearling ultrasound measurements of fat thickness, longissimus muscle area, and estimated percentage of intramuscular fat were available on Simmental bulls (n = 3,409) and heifers (n = 1,503). An extended pedigree was used to construct the relationship matrix (n = 23,968) linking bulls and heifers with ultrasound data to steers and heifers with carcass data. All data were obtained from the American Simmental Association. No animal had both ultrasound and carcass data. Using an animal model and treating corresponding ultrasound and carcass traits separately, genetic parameters were estimated using restricted maximum likelihood. Heritability estimates for carcass traits were 0.48 +/- 0.06, 0.35 +/- 0.05, 0.46 +/- 0.05, and 0.54 +/- 0.05 for carcass weight, fat thickness, longissimus muscle area, and marbling score, respectively. Heritability estimates for bull (heifer) ultrasound traits were 0.53 +/- 0.07 (0.69 +/- 0.09), 0.37 +/- 0.06 (0.51 +/- 0.09), and 0.47 +/- 0.06 (0.52 +/- 0.09) for fat thickness, longissimus muscle area, and intramuscular fat percentage, respectively. Heritability of weight at scan was 0.47 +/- 0.05. Using a bivariate weight model including scan weight of bulls and heifers with carcass weight of slaughter animals, a genetic correlation of 0.77 +/- 0.10 was obtained. Models for fat thickness, longissimus muscle area, and marbling score were each trivariate, including ultrasound measurements on yearling bulls and heifers, and corresponding carcass traits of slaughter animals. Genetic correlations of carcass fat thickness with bull and heifer ultrasound fat were 0.79 +/- 0.13 and 0.83 +/- 0.12, respectively. Genetic correlations of carcass longissimus muscle area with bull and heifer ultrasound longissimus muscle area were 0.80 +/- 0.11 and 0.54 +/- 0.12, respectively. Genetic correlations of carcass marbling score with bull and heifer ultrasound intramuscular fat percentage were 0.74 +/- 0.11 and 0.69 +/- 0.13, respectively. These results provide the parameter estimates necessary for genetic evaluation of Simmental carcass merit using both data from steer and heifer carcasses, and their ultrasound indicators on yearling bulls and heifers.

An autoregressive repeatability animal model for test-day records in multiple lactations.

Test-day (TD) models are becoming a standard for genetic evaluation of production traits in dairy cattle. Various approaches to model covariances between TD records include random regression, autoregressive repeatability, orthogonal polynomials, and models based on character processing. The applicability of these models is mainly associated with the number of parameters to estimate, incorporation of multiple lactations, and the accuracy of correlations generated by the cow's repeated expression of milking performance (TD yields) within and across lactations. We define and evaluate a multiple-lactation, autoregressive-repeatability model that disentangles environmental effects due to cow within and between lactations. Simulated records either included or ignored a long-term environmental effect between lactations. Our autoregressive TD animal model correctly detected presence and the absence of this effect and accurately recovered the assumed variance components and correlations underlying the data (10 parameters for three lactations). Estimates of variance components and autocorrelation coefficients were obtained using DFREML-simplex methodology. Given the value of this approach to reduce the size of residual variance components, autoregressive animal models are a preferable alternative to classical methods based on cumulative lactation yield to improve milk production in dairy cattle.

Genetic antagonism between body weight and milk production in beef cattle.

Korean cattle have an unusually short suckling period (4 mo) due to poor milking ability, and this is a hindrance to growth of calves. Therefore, Korean cattle breeders have shown interest in genetic improvement of milking ability. In this study, body weight (birth weight, weaning weight, and yearling weight) and five daily milk yields by period in Korean cattle (Hanwoo) were analyzed using a two-trait sire and maternal grandsire mixed model. The milk yields used were actually measured at sequential intervals from 1 to 4 mo after calving. Posterior means of the parameters were estimated using Gibbs sampling. Heritability estimates (0.25 to 0.26) for daily milk yield at weaning were larger than those with other periods. Genetic impact on daily milk yield, especially at weaning, was emphasized in order to lengthen the suckling period of Korean cattle. Genetic correlation estimates between BW and daily milk yield were all negative (-0.08 to -0.16 for birth weight, -0.04 to -0.21 for weaning weight, and -0.12 to -0.19 for yearling weight), whereas environmental correlation estimates were all positive (0.20 to 0.39 for birth weight, 0.34 to 0.51 for weaning weight, and 0.30 to 0.45 for yearling weight). The negative estimates of genetic correlation between weight and milk yield implied genetic antagonism between direct and maternal effects for weaning weight of beef cattle.

Genetic analysis of Holstein cattle populations in Brazil and the United States.

Genetic relationships between Brazilian and US Holstein cattle populations were studied using first-lactation records of 305-d mature equivalent (ME) yields of milk and fat of daughters of 705 sires in Brazil and 701 sires in the United States, 358 of which had progeny in both countries. Components of(co)variance and genetic parameters were estimated from all data and from within herd-year standard deviation for milk (HYSD) data files using bivariate and multivariate sire models and DFREML procedures distinguishing the two countries. Sire (residual) variances from all data for milk yield were 51 to 59% (58 to 101%) as large in Brazil as those obtained from half-sisters in the average US herd. Corresponding proportions of the US variance in fat yield that were found in Brazil were 30 to 41% for the sire component of variance and 48 to 80% for the residual. Heritabilities for milk and fat yields from multivariate analysis of all the data were 0.25 and 0.22 in Brazil, and 0.34 and 0.35 in the United States. Genetic correlations between milk and fat were 0.79 in Brazil and 0.62 in the United States. Genetic correlations between countries were 0.85 for milk, 0.88 for fat, 0.55 for milk in Brazil and fat in the US, and 0.67 for fat in Brazil and milk in the United States. Correlated responses in Brazil from sire selection based on the US information increased with average HYSD in Brazil. Largest daughter yield response was predicted from information from half-sisters in low HYSD US herds (0.75 kg/kg for milk; 0.63 kg/kg for fat), which was 14% to 17% greater than estimates from all US herds because the scaling effects were less severe from heterogeneous variances. Unequal daughter response from unequal genetic (co)variances under restrictive Brazilian conditions is evidence for the interaction of genotype and environment. The smaller and variable yield expectations of daughters of US sires in Brazilian environments suggest the need for specific genetic improvement strategies in Brazilian Holstein herds. A US data file restricting daughter information to low HYSD US environments would be a wise choice for across-country evaluation. Procedures to incorporate such foreign evaluations should be explored to improve the accuracy of genetic evaluations for the Brazilian Holstein population.

Variation in the volume of zebra finch song control nuclei is heritable: developmental and evolutionary implications.

In many songbird species, females prefer males that sing a larger repertoire of syllables. Males with more elaborate songs have a larger high vocal centre (HVC) nucleus, the highest structure in the song production pathway. HVC size is thus a potential target of sexual selection. Here we provide evidence that the size of the HVC and other song production nuclei are heritable across individual males within a species. In contrast, we find that heritabilities of other nuclei in a song-learning pathway are lower, suggesting that variation in the sizes of these structures is more closely tied to developmental and environmental differences between individuals. We find that evolvability, a statistical measure that predicts response to selection, is higher for the HVC and its target for song production, the robustus archistriatalis (RA), than for all other brain volumes measured. This suggests that selection based on the functions of these two structures would result in rapid major shifts in their anatomy. We also show that the size of each song control nucleus is significantly correlated with the song related nuclei to which it is monosynaptically connected. Finally, we find that the volume of the telencephalon is larger in males than in females. These findings begin to join theoretical analyses of the role of female choice in the evolution of bird song to neurobiological mechanisms by which the evolutionary changes in behaviour are expressed.

Application of an autoregressive process to estimate genetic parameters and breeding values for daily milk yield in a tropical herd of Lucerna cattle and in United States Holstein herds.

The objectives of this study were to estimate from test day records the genetic and environmental (co)variance components, correlations, and breeding values to increase genetic gain in milk yield of Lucerna and US Holstein cattle. The effects of repeated observations (within cow) were explained by first-order autoregressive processes within and across lactations using an animal model. Estimates of variance components and correlation coefficients between test days were obtained using derivative-free REML methodology. The autoregressive structure significantly reduced the model error component by disentangling the short-term environmental effects. The additional information and the more heterogeneous environmental variances between lactations in the multiple-lactation test day model than in the first lactation model provided substantially larger estimates of additive genetic variance (0.62 kg2 for Lucerna; 14.73 kg2 for Holstein), heritability (0.13 for Lucerna; 0.42 for Holstein), and individual genetic merit. Rank correlations of breeding values from multiple lactations and from first lactations ranged from 0.18 to 0.37 for females and from 0.73 to 0.89 for males, respectively. Consequently, more selection errors and less genetic gain would be expected from selection decisions based on an analysis of first lactation only, and greater accuracy would be achieved from multiple lactations. Results indicated that substantial genetic gain was possible for milk yield in the Lucerna herd (34 kg/yr). Estimates of genetic variance for Holsteins were larger than previously reported, which portends more rapid genetic progress in US herds also; under our assumptions, increases would be from 173 to 197 kg/yr.