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.

Multiple-breed genetic models recently have been demonstrated to account for the heterogenous genetic variances that exist between different beef cattle breed groups. We extend these models to allow for residual heteroskedasticity (heterogeneous residual variances), specified as a function of fixed effects (e.g., sex, breed proportion, breed group heterozygosity) and random effects such as contemporary groups (CG). We additionally specify the residual distributions to be either Gaussian or based on heavier-tailed alternatives such as the Student's t or Slash densities. For each of these three residual densities using either homoskedastic (homogeneous variance) or heteroskedastic error specifications, we analyzed 22,717 postweaning gain records from a Nelore-Hereford population based on a Markov chain Monte Carlo animal model implementation. The heteroskedastic Student's t error model (with estimated df = 7.33 +/- 0.48) was clearly the best-fitting model based on a pseudo-Bayes factor criterion. Breed group heterozygosity and, to a lesser extent, calf sex seemed to be marginally important sources of residual heteroskedasticity. Specifically, the residual variance in F1 animals was estimated to be 0.70 +/- 0.16 times that for purebreds, whereas the male residual variance was estimated to be 1.13 +/- 0.09 times that for females. The CG effects were important random sources of residual heteroskedasticity (i.e., the coefficient of variation of CG-specific residual variances was estimated to be 0.72 +/- 0.06). Purebred Nelores were estimated to have a larger genetic variance (124.84 +/- 21.75 kg2) compared with Herefords (40.89 +/- 6.70 kg2) under the heteroskedastic Student's t error model; however, the converse was observed from results based on a homoskedastic Student's t error model (46.24 +/- 10.90 kg2 and 60.11 +/- 8.54 kg2, respectively). These results indicate that allowing for robustness to outliers and accounting for heteroskedasticity of residual variances has potentially important implications for variance component and genetic parameter estimates from data on multiple-breed populations.

Pigs from the F(2) generation of a Duroc x Pietrain resource population were evaluated to discover QTL affecting growth and composition traits. Body weight and ultrasound estimates of 10th-rib backfat, last-rib backfat, and LM area were serially measured throughout development. Estimates of fat-free total lean, total body fat, empty body protein, empty body lipid, and ADG from 10 to 22 wk of age were calculated, and random regression analyses were performed to estimate individual animal phenotypes representing intercept and linear rates of increase in these serial traits. A total of 510 F(2) animals were genotyped for 124 micro-satellite markers evenly spaced across the genome. Data were analyzed with line cross, least squares regression, interval mapping methods using sex and litter as fixed effects. Significance thresholds of the F-statistic for single QTL with additive, dominance, or imprinted effects were determined at the chromosome- and genome-wise levels by permutation tests. A total of 43 QTL for 22 of the 29 measured traits were found to be significant at the 5% chromosome-wise level. Of these 43 QTL, 20 were significant at the 1% chromosome-wise significance threshold, 14 of these 20 were also significant at the 5% genome-wise significance threshold, and 10 of these 14 were also significant at the 1% genome-wise significance threshold. A total of 22 QTL for the animal random regression terms were found to be significant at the 5% chromosome-wise level. Of these 22 QTL, 6 were significant at the 1% chromosome-wise significance threshold, 4 of these 6 were also significant at the 5% genome-wise significance threshold, and 3 of these 4 were also significant at the 1% genome-wise significance threshold. Putative QTL were discovered for 10th-rib and last-rib backfat on SSC 6, body composition traits on SSC 9, backfat and lipid composition traits on SSC 11, 10th-rib backfat and total body fat tissue on SSC 12, and linear regression of last-rib backfat and total body fat tissue on SSC 8. These results will facilitate fine-mapping efforts to identify genes controlling growth and body composition of pigs that can be incorporated into marker-assisted selection programs to accelerate genetic improvement in pig populations.

Background: Recent evidence suggests that some sex differences in brain and behavior might result from direct genetic effects, and not solely the result of the organizational effects of steroid hormones. The present study examined the potential role for sex-biased gene expression during development of sexually dimorphic singing behavior and associated song nuclei in juvenile zebra finches.