Our long-term objective is to develop breeding strategies for improving feed efficiency in dairy cattle. In this study, phenotypic data were pooled across multiple research stations to facilitate investigation of the genetic and nongenetic components of feed efficiency in Holstein cattle. Specifically, the heritability of residual feed intake (RFI) was estimated and heterogeneous relationships between RFI and traits relating to energy utilization were characterized across research stations. Milk, fat, protein, and lactose production converted to megacalories (milk energy; MilkE), dry matter intakes (DMI), and body weights (BW) were collected on 6,824 lactations from 4,893 Holstein cows from research stations in Scotland, the Netherlands, and the United States. Weekly DMI, recorded between 50 to 200 d in milk, was fitted as a linear function of MilkE, BW0.75, and change in BW (ΔBW), along with parity, a fifth-order polynomial on days in milk (DIM), and the interaction between this polynomial and parity in a first-stage model. The residuals from this analysis were considered to be a phenotypic measure of RFI. Estimated partial regression coefficients of DMI on MilkE and on BW0.75 ranged from 0.29 to 0.47 kg/Mcal for MilkE across research stations, whereas estimated partial regression coefficients on BW0.75 ranged from 0.06 to 0.16 kg/kg0.75. Estimated partial regression coefficients on ΔBW ranged from 0.06 to 0.39 across stations. Heritabilities for country-specific RFI were based on fitting second-stage random regression models and ranged from 0.06 to 0.24 depending on DIM. The overall heritability estimate across all research stations and all DIM was 0.15±0.02, whereas an alternative analysis based on combining the first- and second-stage model as 1 model led to an overall heritability estimate of 0.18±0.02. Hence future genomic selection programs on feed efficiency appear to be promising; nevertheless, care should be taken to allow for potentially heterogeneous variance components and partial relationships between DMI and other energy sink traits across environments when determining RFI.

Feed efficiency, as defined by the fraction of feed energy or dry matter captured in products, has more than doubled for the US dairy industry in the past 100 yr. This increased feed efficiency was the result of increased milk production per cow achieved through genetic selection, nutrition, and management with the desired goal being greater profitability. With increased milk production per cow, more feed is consumed per cow, but a greater portion of the feed is partitioned toward milk instead of maintenance and body growth. This dilution of maintenance has been the overwhelming driver of enhanced feed efficiency in the past, but its effect diminishes with each successive increment in production relative to body size and therefore will be less important in the future. Instead, we must also focus on new ways to enhance digestive and metabolic efficiency. One way to examine variation in efficiency among animals is residual feed intake (RFI), a measure of efficiency that is independent of the dilution of maintenance. Cows that convert feed gross energy to net energy more efficiently or have lower maintenance requirements than expected based on body weight use less feed than expected and thus have negative RFI. Cows with low RFI likely digest and metabolize nutrients more efficiently and should have overall greater efficiency and profitability if they are also healthy, fertile, and produce at a high multiple of maintenance. Genomic technologies will help to identify these animals for selection programs. Nutrition and management also will continue to play a major role in farm-level feed efficiency. Management practices such as grouping and total mixed ration feeding have improved rumen function and therefore efficiency, but they have also decreased our attention on individual cow needs. Nutritional grouping is key to helping each cow reach its genetic potential. Perhaps new computer-driven technologies, combined with genomics, will enable us to optimize management for each individual cow within a herd, or to optimize animal selection to match management environments. In the future, availability of feed resources may shift as competition for land increases. New approaches combining genetic, nutrition, and other management practices will help optimize feed efficiency, profitability, and environmental sustainability.

Hierarchical mixed effects models have been demonstrated to be powerful for predicting genomic merit of livestock and plants, on the basis of high-density single-nucleotide polymorphism (SNP) marker panels, and their use is being increasingly advocated for genomic predictions in human health. Two particularly popular approaches, labeled BayesA and BayesB, are based on specifying all SNP-associated effects to be independent of each other. BayesB extends BayesA by allowing a large proportion of SNP markers to be associated with null effects. We further extend these two models to specify SNP effects as being spatially correlated due to the chromosomally proximal effects of causal variants. These two models, that we respectively dub as ante-BayesA and ante-BayesB, are based on a first-order nonstationary antedependence specification between SNP effects. In a simulation study involving 20 replicate data sets, each analyzed at six different SNP marker densities with average LD levels ranging from r 2 ¼ 0.15 to 0.31, the antedependence methods had significantly (P , 0.01) higher accuracies than their corresponding classical counterparts at higher LD levels (r 2 . 0. 24) with differences exceeding 3%. A cross-validation study was also conducted on the heterogeneous stock mice data resource (http:// mus.well.ox.ac.uk/mouse/HS/) using 6-week body weights as the phenotype. The antedependence methods increased cross-validation prediction accuracies by up to 3.6% compared to their classical counterparts (P , 0.001). Finally, we applied our method to other benchmark data sets and demonstrated that the antedependence methods were more accurate than their classical counterparts for genomic predictions, even for individuals several generations beyond the training data.

The objectives of this work were to assess alternative linear reaction norm (RN) models for genetic evaluation of Angus cattle in Brazil. That is, we investigated the interaction between genotypes and continuous descriptors of the environmental variation to examine evidence of genotype by environment interaction (G×E) in post-weaning BW gain (PWG) and to compare the environmental sensitivity of national and imported Angus sires. Data were collected by the Brazilian Angus Improvement Program from 1974 to 2005 and consisted of 63,098 records and a pedigree file with 95,896 animals. Six models were implemented using Bayesian inference and compared using the Deviance Information Criterion (DIC). The simplest model was M(1), a traditional animal model, which showed the largest DIC and hence the poorest fit when compared with the 4 alternative RN specifications accounting for G×E. In M(2), a 2-step procedure was implemented using the contemporary group posterior means of M(1) as the environmental gradient, ranging from -92.6 to +265.5 kg. Moreover, the benefits of jointly estimating all parameters in a 1-step approach were demonstrated by M(3). Additionally, we extended M(3) to allow for residual heteroskedasticity using an exponential function (M(4)) and the best fitting (smallest DIC) environmental classification model (M(5)) specification. Finally, M(6) added just heteroskedastic residual variance to M(1). Heritabilities were less at harsh environments and increased with the improvement of production conditions for all RN models. Rank correlations among genetic merit predictions obtained by M(1) and by the best fitting RN models M(3) (homoskedastic) and M(5) (heteroskedastic) at different environmental levels ranged from 0.79 and 0.81, suggesting biological importance of G×E in Brazilian Angus PWG. These results suggest that selection progress could be optimized by adopting environment-specific genetic merit predictions. The PWG environmental sensitivity of imported North American origin bulls (0.046 ± 0.009) was significantly larger (P < 0.05) than that of local sires (0.012 ± 0.013). Moreover, PWG of progeny of imported sires exceeded that of native sires in medium and superior production levels. On the other hand, Angus cattle locally selected in Brazil tended to be more robust to environmental changes and hence be more suitable when production environments for potential progeny is uncertain.

Reliable biomarkers predictive of productive herd life (time in herd after birth of first calf) have heretofore not been discovered in dairy cattle. However, circulating concentrations of anti-Müllerian hormone (AMH) are positively associated with number of follicles or antral follicle count (AFC), ovarian function, and fertility, and approximately 25% of cows have a relatively low AFC and low AMH concentrations. The present study tested the hypothesis that heifers with the lowest AMH concentrations have suboptimal fertility and are removed from a herd for poor reproductive performance at a greater rate, and therefore have a shorter productive herd life compared with age-matched herdmates with higher AMH. To test this hypothesis, 11- to 15-mo-old Holstein heifers (n=281) were subjected to a single measurement of AMH. All heifers not removed from the herd had the opportunity to complete 2 lactations and start their third lactation after calving. During this time, performance and health parameters for each individual were recorded daily by herd managers. Results showed that the quartile of heifers with the lowest AMH concentration also had, on average, a shorter productive herd life (by 196 d), a reduced survival rate after birth of the first calf, the lowest level of milk production (first lactation), the lowest total percentage of cows pregnant (across all lactations), the highest culling rates (first and second lactations and overall), and the highest culling rate for poor reproduction (first lactation) compared with age-matched herdmates with higher AMH. We concluded that a single determination of AMH concentration in young adult dairy heifers may be a simple diagnostic method to predict herd longevity, and AMH may be a useful phenotypic marker to improve longevity of dairy cows.