Genomic prediction of residual feed intake in US Holstein dairy cattle.

Genomic selection is an important tool to introduce feed efficiency into dairy cattle breeding. The goals of the current research are to estimate genomic breeding values of residual feed intake (RFI) and to assess the prediction reliability for RFI in the US Holstein population. The RFI data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States, and were pre-adjusted to remove phenotypic correlations with milk energy, metabolic body weight, body weight change, and for several environmental effects. In the current analyses, genomic predicted transmitting abilities of milk energy and of body weight composite were included into the RFI model to further remove the genetic correlations that remained between RFI and these energy sinks. In the first part of the analyses, a national genomic evaluation for RFI was conducted for all the Holsteins in the national database using a standard multi-step genomic evaluation method and 60,671 SNP list. In the second part of the study, a single-step genomic prediction method was applied to estimate genomic breeding values of RFI for all cows with phenotypes, 5,252 elite young bulls, 4,029 young heifers, as well as their ancestors in the pedigree, using a high-density genotype chip. Theoretical prediction reliabilities were calculated for all the studied animals in the single-step genomic prediction by direct inversion of the mixed model equations. In the results, breeding values were estimated for 1.6 million genotyped Holsteins and 60 million ungenotyped Holsteins, The genomic predicted transmitting ability correlations between RFI and other traits in the index (e.g., fertility) are generally low, indicating minor correlated responses on other index traits when selecting for RFI. Genomic prediction reliabilities for RFI averaged 34% for all phenotyped animals and 13% for all 1.6 million genotyped animals. Including genomic information increased the prediction reliabilities for RFI compared with using only pedigree information. All bulls had low reliabilities, and averaged to only 16% for the top 100 net merit progeny-tested bulls. Analyses using single-step genomic prediction and high-density genotypes gave similar results to those obtained from the national evaluation. The average theoretical reliability for RFI was 18% among the elite young bulls under 5 yr old, being lower in the younger generations of elite bulls compared with older bulls. To conclude, the size of the reference population and its relationship to the predicted population remain as the limiting factors in the genomic prediction for RFI. Continued collection of feed intake data is necessary so that reliabilities can be maintained due to close relationships of phenotyped animals with breeding stock. Considering the currently low prediction reliability and high cost of data collection, focusing RFI data collection on relatives of elite bulls that will have the greatest genetic contribution to the next generation will give more gains and profit.

Defining the optimal period length and stage of growth or lactation to estimate residual feed intake in dairy cows.

Residual feed intake (RFI) is an estimate of animal feed efficiency, calculated as the difference between observed and expected feed intake. Expected intake typically is derived from a multiple regression model of dry matter intake on energy sinks, including maintenance and growth in growing animals, or maintenance, gain in body reserves, and milk production in lactating animals. The best period during the production cycle of a dairy cow to estimate RFI is not clear. Here, we characterized RFI in growing Holstein heifers (RFIGrowth; ∼10 to 14 mo of age; n = 226) and cows throughout a 305-d lactation (RFILac-Full; n = 118). The goals were to characterize relationships between RFI estimated at different production stages of the dairy cow; determine effects of selection for efficiency during growth on subsequent lactation and feed efficiency; and identify the most desirable testing scheme for RFILac-Full. For RFIGrowth, intake was predicted from multiple linear regression of metabolizable energy (ME) intake on mid-test body weight (BW)0.75 and average daily gain (ADG). For RFILac-Full, predicted intake was based on regression of BW0.75, ADG, and energy-corrected milk yield. Mean energy intake of the least and most efficient growing heifers (±0.5 standard deviations from mean RFIGrowth of 0) differed by 3.01 Mcal of ME/d, but the groups showed no difference in mid-test BW or ADG. Phenotypic correlation between RFIGrowth and RFI of heifers estimated in the first 100 d in milk (RFILac100DIM; n = 130) was 0.37. Ranking of these heifers as least (mean + 0.5 standard deviations), middle, or most efficient (mean - 0.5 standard deviations) based on RFIGrowth resulted in 43% maintaining the same ranking by RFILac100DIM. On average, the most efficient heifers ate 3.27 Mcal of ME/d less during the first 100 DIM than the least efficient heifers, but exhibited no differences in average energy-corrected milk yield, ADG, or BW. The correlation between RFILac100DIM and RFILac-Full was 0.72. Thus, RFIGrowth may serve as an indicator trait for RFI during lactation, and selection for heifers exhibiting low RFIGrowth should improve overall herd feed efficiency during lactation. Correlation analysis between RFILac-Full (10 to 305 DIM) and subperiod estimates of RFI during lactation indicated a test period of 64 to 70 d in duration occurring between 150 to 220 DIM provided a reliable approximation (r ≥ 0.90) of RFILac-Full among the test periods evaluated.

High-density genome-wide association study for residual feed intake in Holstein dairy cattle.

Improving feed efficiency (FE) of dairy cattle may boost farm profitability and reduce the environmental footprint of the dairy industry. Residual feed intake (RFI), a candidate FE trait in dairy cattle, can be defined to be genetically uncorrelated with major energy sink traits (e.g., milk production, body weight) by including genomic predicted transmitting ability of such traits in genetic analyses for RFI. We examined the genetic basis of RFI through genome-wide association (GWA) analyses and post-GWA enrichment analyses and identified candidate genes and biological pathways associated with RFI in dairy cattle. Data were collected from 4,823 lactations of 3,947 Holstein cows in 9 research herds in the United States. Of these cows, 3,555 were genotyped and were imputed to a high-density list of 312,614 SNP. We used a single-step GWA method to combine information from genotyped and nongenotyped animals with phenotypes as well as their ancestors' information. The estimated genomic breeding values from a single-step genomic BLUP were back-solved to obtain the individual SNP effects for RFI. The proportion of genetic variance explained by each 5-SNP sliding window was also calculated for RFI. Our GWA analyses suggested that RFI is a highly polygenic trait regulated by many genes with small effects. The closest genes to the top SNP and sliding windows were associated with dry matter intake (DMI), RFI, energy homeostasis and energy balance regulation, digestion and metabolism of carbohydrates and proteins, immune regulation, leptin signaling, mitochondrial ATP activities, rumen development, skeletal muscle development, and spermatogenesis. The region of 40.7 to 41.5 Mb on BTA25 (UMD3.1 reference genome) was the top associated region for RFI. The closest genes to this region, CARD11 and EIF3B, were previously shown to be related to RFI of dairy cattle and FE of broilers, respectively. Another candidate region, 57.7 to 58.2 Mb on BTA18, which is associated with DMI and leptin signaling, was also associated with RFI in this study. Post-GWA enrichment analyses used a sum-based marker-set test based on 4 public annotation databases: Gene Ontology, Kyoto Encyclopedia of Genes and Genomes (KEGG) pathways, Reactome pathways, and medical subject heading (MeSH) terms. Results of these analyses were consistent with those from the top GWA signals. Across the 4 databases, GWA signals for RFI were highly enriched in the biosynthesis and metabolism of amino acids and proteins, digestion and metabolism of carbohydrates, skeletal development, mitochondrial electron transport, immunity, rumen bacteria activities, and sperm motility. Our findings offer novel insight into the genetic basis of RFI and identify candidate regions and biological pathways associated with RFI in dairy cattle.

The effects of feeding mixed tocopherol oil on whole-blood respiratory burst and neutrophil immunometabolic-related gene expression in lactating dairy cows.

The 4 major tocopherol isoforms differ in their biochemical reactivity and cellular effects due to basic chemical structural differences. Alpha-tocopherol has been well studied regarding effects on bovine polymorphonuclear leukocyte (PMN) function and its involvement in respiratory burst. However, no studies to date have identified the effects of supplementing a mixed tocopherol oil (Tmix) particularly enriched in non-α tocopherol isoforms (i.e., γ- and δ-isoforms) on fundamental immunometabolic changes in dairy cows. Therefore, the objectives of this study were to determine whether short-term feeding of vegetable oil-derived Tmix alters specific biomarkers of metabolism, whole-blood leukocyte populations, respiratory burst, immunometabolic-related gene expression of PMN, or gene expression of isolated PMN when challenged with lipopolysaccharides (LPS). Clinically healthy multiparous lactating Holstein cows (n = 12; 179 ± 17 d in milk, 40.65 ± 3.68 kg of milk yield) were fed Tmix (620 g/d) for 7 consecutive days. Jugular blood (EDTA anticoagulant) was collected from all cows on d 0 before treatment initiation and again on d 7 after Tmix feeding. Total stimulated respiratory burst activity (RBA) and leukocyte populations were assessed in whole blood, and tocopherol isoform concentrations, metabolites, and hormones were measured in plasma. For gene expression analysis, isolated PMN from cows before and after Tmix feeding were incubated with LPS at a final concentration of either 0.0 or 1.5 µg/mL. Feeding of Tmix for 7 d increased the concentrations of α- and γ-tocopherol. The Tmix did not alter plasma insulin but decreased cholesterol. The Tmix did not alter whole-blood RBA or the leukocyte populations. The LPS challenge increased the expression of proinflammatory genes TNFA and IL6. However, Tmix treatment did not alter the patterns of LPS-affected expression of genes (e.g., TNFA, ITGB2, PPARA, and RXRA) associated with the immune or metabolic response. In conclusion, short-term feeding of Tmix may have no negative effect on animal health as Tmix increased α- and γ-tocopherol concentrations in blood and did not impair whole-blood RBA or alter leukocyte populations. The data provide further support that the α- and γ-tocopherol isoforms do not interfere with normal immune or metabolic function.

Genome-wide association analyses based on a multiple-trait approach for modeling feed efficiency.

Genome-wide association (GWA) of feed efficiency (FE) could help target important genomic regions influencing FE. Data provided by an international dairy FE research consortium consisted of phenotypic records on dry matter intakes (DMI), milk energy (MILKE), and metabolic body weight (MBW) on 6,937 cows from 16 stations in 4 counties. Of these cows, 4,916 had genotypes on 57,347 single nucleotide polymorphism (SNP) markers. We compared a GWA analysis based on the more classical residual feed intake (RFI) model with one based on a previously proposed multiple trait (MT) approach for modeling FE using an alternative measure (DMI|MILKE,MBW). Both models were based on a single-step genomic BLUP procedure that allowed the use of phenotypes from both genotyped and nongenotyped cows. Estimated effects for single SNP markers were small and not statistically important but virtually identical for either FE measure (RFI vs. DMI|MILKE,MBW). However, upon further refining this analysis to develop joint tests within nonoverlapping 1-Mb windows, significant associations were detected between either measure of FE with a window on each of Bos taurus autosomes BTA12 and BTA26. There was, as expected, no overlap between detected genomic regions for DMI|MILKE,MBW and genomic regions influencing the energy sink traits (i.e., MILKE and MBW) because of orthogonal relationships clearly defined between the various traits. Conversely, GWA inferences on DMI can be demonstrated to be partly driven by genetic associations between DMI with these same energy sink traits, thereby having clear implications when comparing GWA studies on DMI to GWA studies on FE-like measures such as RFI.

Use of genotype × environment interaction model to accommodate genetic heterogeneity for residual feed intake, dry matter intake, net energy in milk, and metabolic body weight in dairy cattle.

Feed efficiency in dairy cattle has gained much attention recently. Due to the cost-prohibitive measurement of individual feed intakes, combining data from multiple countries is often necessary to ensure an adequate reference population. It may then be essential to model genetic heterogeneity when making inferences about feed efficiency or selecting efficient cattle using genomic information. In this study, we constructed a marker × environment interaction model that decomposed marker effects into main effects and interaction components that were specific to each environment. We compared environment-specific variance component estimates and prediction accuracies from the interaction model analyses, an across-environment analyses ignoring population stratification, and a within-environment analyses using an international feed efficiency data set. Phenotypes included residual feed intake, dry matter intake, net energy in milk, and metabolic body weight from 3,656 cows measured in 3 broadly defined environments: North America (NAM), the Netherlands (NLD), and Scotland (SAC). Genotypic data included 57,574 single nucleotide polymorphisms per animal. The interaction model gave the highest prediction accuracy for metabolic body weight, which had the largest estimated heritabilities ranging from 0.37 to 0.55. The within-environment model performed the best when predicting residual feed intake, which had the lowest estimated heritabilities ranging from 0.13 to 0.41. For traits (dry matter intake and net energy in milk) with intermediate estimated heritabilities (0.21 to 0.50 and 0.17 to 0.53, respectively), performance of the 3 models was comparable. Genomic correlations between environments also were computed using variance component estimates from the interaction model. Averaged across all traits, genomic correlations were highest between NAM and NLD, and lowest between NAM and SAC. In conclusion, the interaction model provided a novel way to evaluate traits measured in multiple environments in which genetic heterogeneity may exist. This model allowed estimation of environment-specific parameters and provided genomic predictions that approached or exceeded the accuracy of competing within- or across-environment models.

Reducing gut effects from Cryptosporidium parvum infection in dairy calves through prophylactic glucagon-like peptide 2 therapy or feeding of an artificial sweetener.

Glucagon-like peptide 2 (GLP-2) therapy was shown previously to reduce inflammation-related gut damage from coccidiosis in dairy calves, and feeding of artificial sweetener stimulates GLP-2 secretion from intestinal L cells. The purpose of this study was to determine whether GLP-2 treatment or artificial sweetener feeding beginning 1 wk before an experimental inoculation with the coccidian parasite Cryptosporidium parvum can reduce infection-related intestinal damage in Holstein bull calves. Newborn calves were assigned to 1 of 4 treatment groups of 6 calves each, including noninfected control calves injected s.c. every 12 h with control buffer (CON), infected control calves injected s.c. every 12 h with control buffer (INF), infected calves injected s.c. every 12 h with 50 µg/kg of body weight of GLP-2 (GLP2), and infected calves injected s.c. every 12 h with control buffer and supplemented in the diet with Sucram (Pancosma, Geneva, Switzerland) at 400 mg/kg of dry matter of milk replacer (SUC). Treatments were initiated on d 1, and calves in INF, GLP2, and SUC were orally dosed on d 8 with 12,500 C. parvum oocysts. Fecal scores were recorded daily, plasma was collected on d 1, 8, 12, 15, and 18 to evaluate markers of inflammation, and fecal samples were collected on d 1, 8, and every other day thereafter to determine the presence of oocysts. Calves were euthanized on d 18 for collection of intestinal tissues and histological and gene expression analyses. Relative to CON, calves in INF exhibited an increase in diarrhea severity, increased plasma serum amyloid A concentration on d 15 and 18, reduced intestinal villus height, increased villus apoptosis and crypt cell proliferation, and increased intestinal mRNA expression of MARVELD2 and GPX2. However, calves in SUC and GLP2 had reduced diarrhea severity and fecal C. parvum oocyst shedding, reduced plasma serum amyloid A concentration on d 15 and 18, and, depending on the intestinal segment, increased villus height, reduced crypt cell proliferation, and reduced mRNA expression of MARVELD2, GPX2, and other tight junction proteins relative to INF. Lastly, GLP2 and SUC exhibited increased intestinal mass-to-length ratio and decreased length-to-empty body weight ratio relative to INF. Our findings suggest that GLP-2 and Sucram treatments administered before a low-level C. parvum exposure may contribute to fewer effects on intestinal integrity, morphology, and inflammation in response to infection, and shorter, denser intestines.

Modeling genetic and nongenetic variation of feed efficiency and its partial relationships between component traits as a function of management and environmental factors.

Feed efficiency (FE), characterized as the fraction of feed nutrients converted into salable milk or meat, is of increasing economic importance in the dairy industry. We conjecture that FE is a complex trait whose variation and relationships or partial efficiencies (PE) involving the conversion of dry matter intake to milk energy and metabolic body weight may be highly heterogeneous across environments or management scenarios. In this study, a hierarchical Bayesian multivariate mixed model was proposed to jointly infer upon such heterogeneity at both genetic and nongenetic levels on PE and variance components (VC). The heterogeneity was modeled by embedding mixed effects specifications on PE and VC in addition to those directly specified on the component traits. We validated the model by simulation and applied it to a joint analysis of a dairy FE consortium data set with 5,088 Holstein cows from 13 research stations in Canada, the Netherlands, the United Kingdom, and the United States. Although no differences were detected among research stations for PE at the genetic level, some evidence was found of heterogeneity in residual PE. Furthermore, substantial heterogeneity in VC across stations, parities, and ration was observed with heritability estimates of FE ranging from 0.16 to 0.46 across stations.

The genetic and biological basis of feed efficiency in mid-lactation Holstein dairy cows.

The objective of this study was to identify genomic regions and candidate genes associated with feed efficiency in lactating Holstein cows. In total, 4,916 cows with actual or imputed genotypes for 60,671 single nucleotide polymorphisms having individual feed intake, milk yield, milk composition, and body weight records were used in this study. Cows were from research herds located in the United States, Canada, the Netherlands, and the United Kingdom. Feed efficiency, defined as residual feed intake (RFI), was calculated within location as the residual of the regression of dry matter intake (DMI) on milk energy (MilkE), metabolic body weight (MBW), change in body weight, and systematic effects. For RFI, DMI, MilkE, and MBW, bivariate analyses were performed considering each trait as a separate trait within parity group to estimate variance components and genetic correlations between them. Animal relationships were established using a genomic relationship matrix. Genome-wide association studies were performed separately by parity group for RFI, DMI, MilkE, and MBW using the Bayes B method with a prior assumption that 1% of single nucleotide polymorphisms have a nonzero effect. One-megabase windows with greatest percentage of the total genetic variation explained by the markers (TGVM) were identified, and adjacent windows with large proportion of the TGVM were combined and reanalyzed. Heritability estimates for RFI were 0.14 (±0.03; ±SE) in primiparous cows and 0.13 (±0.03) in multiparous cows. Genetic correlations between primiparous and multiparous cows were 0.76 for RFI, 0.78 for DMI, 0.92 for MBW, and 0.61 for MilkE. No single 1-Mb window explained a significant proportion of the TGVM for RFI; however, after combining windows, significance was met on Bos taurus autosome 27 in primiparous cows, and nearly reached on Bos taurus autosome 4 in multiparous cows. Among other genes, these regions contain ß-3 adrenergic receptor and the physiological candidate gene, leptin, respectively. Between the 2 parity groups, 3 of the 10 windows with the largest effects on DMI neighbored windows affecting RFI, but were not in the top 10 regions for MilkE or MBW. This result suggests a genetic basis for feed intake that is unrelated to energy consumption required for milk production or expected maintenance as determined by MBW. In conclusion, feed efficiency measured as RFI is a polygenic trait exhibiting a dynamic genetic basis and genetic variation distinct from that underlying expected maintenance requirements and milk energy output.

Short communication: Glucagon-like peptide-2 and coccidiosis alter tight junction gene expression in the gastrointestinal tract of dairy calves.

Tight junction (TJ) proteins are integral factors involved in gut barrier function, and therapy with glucagon-like peptide-2 (GLP-2) enhances gut integrity. Our aim was to assess effects of GLP-2 treatment on mRNA expression of 8 TJ complex proteins in the intestine of dairy calves not infected or infected with Eimeria bovis at 11±3d of age. Mucosal epithelium from jejunum, ileum, and cecum was collected at slaughter from Holstein bull calves assigned to 4 groups: noninfected, buffer-treated (n=5); noninfected, GLP-2 treated (n=4); E. bovis-infected, buffer-treated (n=5); and E. bovis-infected, GLP-2-treated (n=4). Infected calves were orally dosed with 100,000 to 200,000 sporulated E. bovis oocysts on d 0; GLP-2-treated calves received 50 µg of GLP-2/kg of body weight subcutaneously twice daily for 10d beginning on d 18; and buffer-treated calves received an equal injection volume of 0.01 M Na bicarbonate buffer. All calves were killed on d 28. The mRNA expression of coxsackie and adenovirus receptor (CXADR), claudins 1, 2, and 4 (CLDN1, CLDN2, and CLDN4), F11 receptor (F11R), junction adhesion molecule 2 (JAM2), occludin (OCLN), and tight junction protein ZO-1 (TJP1) was determined by real-time quantitative PCR. In jejunum and ileum, an interaction of E. bovis infection and GLP-2 treatment on gene expression was noted. In jejunum of noninfected calves, GLP-2 increased CXADR, CLDN2, OCLN, and TJP1 mRNA expression but had no effect on mRNA expression in infected calves. Treatment with GLP-2 also increased tight junction protein ZO-1 protein expression in jejunum of noninfected calves as determined by immunohistochemistry. In ileum, E. bovis decreased expression of JAM2, OCLN, and TJP1 in buffer-treated calves, and GLP-2 increased TJP1 expression in infected calves. In cecum, E. bovis infection reduced expression of CXADR, CLDN4, F11R, and OCLN, and GLP-2 therapy increased expression of CLDN4, F11R, OCLN, and TJP1. Results are consistent with studies in nonruminants showing decreased expression of TJ complex proteins in the intestinal tract during pathogen-induced diarrhea and increased TJ protein expression in intestinal tissues in response to GLP-2 treatment. In conclusion, E. bovis reduces gene expression of TJ proteins primarily in cecum of calves 28d postinfection, and GLP-2 increases expression of selected TJ genes in intestinal tissues. Use of GLP-2 to improve gut barrier function in ruminants during pathogen-induced diarrhea warrants additional study.

Heterogeneity in genetic and nongenetic variation and energy sink relationships for residual feed intake across research stations and countries.

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.

An alternative approach to modeling genetic merit of feed efficiency in dairy cattle.

Genetic improvement of feed efficiency (FE) in dairy cattle requires greater attention given increasingly important resource constraint issues. A widely accepted yet occasionally contested measure of FE in dairy cattle is residual feed intake (RFI). The use of RFI is limiting for several reasons, including interpretation, differences in recording frequencies between the various component traits that define RFI, and potential differences in genetic versus nongenetic relationships between dry matter intake (DMI) and FE component traits. Hence, analyses focusing on DMI as the response are often preferred. We propose an alternative multiple-trait (MT) modeling strategy that exploits the Cholesky decomposition to provide a potentially more robust measure of FE. We demonstrate that our proposed FE measure is identical to RFI provided that genetic and nongenetic relationships between DMI and component traits of FE are identical. We assessed both approaches (MT and RFI) by simulation as well as by application to 26,383 weekly records from 50 to 200 d in milk on 2,470 cows from a dairy FE consortium study involving 7 institutions. Although the proposed MT model fared better than the RFI model when simulated genetic and nongenetic associations between DMI and FE component traits were substantially different from each other, no meaningful differences were found in predictive performance between the 2 models when applied to the consortium data.

Hot topic: Brown marmorated stink bug odor compounds do not transfer into milk by feeding bug-contaminated corn silage to lactating dairy cattle.

Brown marmorated stink bug (BMSB; Halyomorpha halys) is an emerging invasive species of grave concern to agriculture as a polyphagous plant pest with potential negative effects on the dairy industry. The purpose of this study was to determine the risk of including BMSB-contaminated silage in lactating dairy cow rations. First, 6 dairies, either highly infested (n=3; 30 to 100 bugs per stalk) or not infested (n=3), were sampled to assess the prevalence of bug secretion compounds tridecane (major component) and E-2-decenal (stink odor component) in silage and milk. Second, using wild BMSB, a mini-silo dose-response experiment (adding 100, 50, 25, 10, and 1 freshly crushed bugs/0.5kg of chopped corn) was conducted to assess the effect of ensiling on BMSB stink odor compounds. Finally, synthetic BMSB stink odor compounds (10g of tridecane and 5g of E-2-decenal) were ruminally infused twice daily over 3 d, and samples of milk, urine, and rumen fluid were collected to evaluate disposition. Bug stink odor compounds were sampled by solid-phase microextraction (SPME) and analyzed by gas chromatography-mass spectrometry (GC-MS). Milk production and feed composition were unaffected when BMSB-contaminated silage was fed. Moreover, no E-2-decenal was detected in silage or milk (detection threshold = 0.00125µg/mL). The dose-response of tridecane in mini-silo samples exhibited a linear relationship (R(2)=0.78) with the amount of BMSB added; however, E-2-decenal was completely decomposed and undetectable in spiked mini-silos after ensiling. Both synthetic secretion compounds infused into rumen were undetectable in all milk and urine samples. E-2-Decenal was not detectable in rumen fluid, whereas tridecane was detected only at 15 min postinfusion but not present thereafter. Feed intake was unaffected by infusion treatment and BMSB secretion compounds (E-2-decenal and tridecane) were not observed in milk. E-2-Decenal and tridecane from the metathoracic gland of BMSB are not able to contaminate milk either due to the ensiling process or because of metabolism within the rumen. Concern over BMSB stink odor compounds contaminating the fluid milk supply, even on highly infested farms, is not warranted.

Transcriptional regulators transforming growth factor-ß1 and estrogen-related receptor-α identified as putative mediators of calf rumen epithelial tissue development and function during weaning.

Molecular mechanisms regulating rumen epithelial development remain largely unknown. To identify gene networks and regulatory factors controlling rumen development, Holstein bull calves (n=18) were fed milk replacer only (MRO) until 42 d of age. Three calves each were euthanized at 14 and 42 d of age for tissue collection to represent preweaning, and the remaining calves were provided diets of either milk replacer + orchard grass hay (MH; n=6) to initiate weaning without development of rumen papillae, or milk replacer + calf starter (MG; n=6) to initiate weaning and development of rumen papillae. At 56 and 70 d of age, 3 calves from the MH and MG groups were euthanized for collection of rumen epithelium. Total RNA and protein were extracted for microarray analysis and to validate detected changes in selected protein expression, respectively. As expected, calves fed MRO had no rumen papillae and development of papillae was greater in MG versus MH calves. Differentially expressed genes between the MRO diet at d 42 (preweaning) versus the MG or MH diets at d 56 (during weaning) were identified using permutation analysis of differential expression. Expression of 345 and 519 transcripts was uniquely responsive to MG and MH feeding, respectively. Ingenuity Pathway Analysis (Qiagen, Redwood City, CA) indicated that the top-ranked biological function affected by the MG diet was the cell cycle, and TFGB1, FBOX01, and PPARA were identified as key transcriptional regulators of genes responsive to the MG diet and associated with development of rumen papillae. Increased expressions of TGFB1 mRNA and protein in response to the MG diet were confirmed by subsequent analyses. The top-ranking biological function affected by the MH diet was energy production. Receptors for IGF-1 and insulin, ESRRA, and PPARD were identified by ingenuity pathway analysis as transcriptional regulators of genes responsive to the MH diet. Further analysis of TGFB1 and ESRRA mRNA expression in rumen epithelium obtained from a separate ontogenic study of Holstein calves (n=26) euthanized every 7d from birth to 42 d of age showed increases in transcript expression with advancing age, supporting their roles in mediating rumen epithelial development and function during weaning. Additional evaluation of gene expression in the rumen epithelium of adult cows ruminally infused with butyrate also suggested that observed changes in ESRRA mRNA expression in developing calf rumen may be mediated by increased butyrate concentration. Our results identify TGFB1 and ESRRA as likely transcriptional regulators of rumen epithelial development and energy metabolism, respectively, and provide targets for modulation of rumen development and function in the growing calf.

Glucagon-like peptide 2 therapy reduces negative effects of diarrhea on calf gut.

Damage to the intestinal epithelium reduces nutrient absorption and animal growth, and can have negative long-term health effects on livestock. Because the intestinotropic hormone glucagon-like peptide 2 (GLP-2) has been shown to contribute to gut integrity, reduce inflammation, and improve nutrient absorption, the present study was designed to determine whether administration of GLP-2 to calves with coccidiosis in the first month of life affects intestinal growth and mediates negative effects of the proinflammatory response. Holstein bull calves (n=19) were assigned to 4 treatment groups of 4 to 5 calves each: (1) infected with Eimeria bovis, GLP-2 treated; (2) noninfected, GLP-2 treated; (3) infected with E. bovis, buffer treated; and (4) noninfected, buffer treated. Infected calves received 100,000 to 200,000 sporulated E. bovis oocysts suspended in milk replacer on d 0 of the study. On d 18, calves in the GLP-2 groups received a subcutaneous injection of 50 µg of bovine GLP-2/kg of body weight twice daily for 10 d, and calves in the buffer-treated groups received an equivalent volume of sodium bicarbonate buffer only. On d 28, calves were slaughtered 2h after injection of 5-bromo-2'-deoxyuridine (BrdU). Intestinal tissues were measured and villus height, crypt depth, and BrdU immunostaining were evaluated in segments of the small intestine. Nitrotyrosine immunostaining, a measure of nitro-oxidative damage, was evaluated in the ileum and cecum. No GLP-2 treatment by E. bovis infection interaction was observed for any parameter measured, with the exception of nitrotyrosine immunostaining in the cecum. Large intestinal weight was greater in infected than noninfected calves and with GLP-2 treatment relative to buffer treatment. Calves that received GLP-2 also had greater small intestinal weight but no difference in cell proliferation, as assessed by BrdU labeling, relative to buffer-treated calves. No treatment effects were detected for villus height, crypt depth, or villus height:crypt depth ratio in segments of the small intestine. Protein tyrosine nitration was over 3-fold greater in the ileum and cecum of infected calves relative to noninfected calves, and GLP-2 therapy reduced tyrosine nitration in infected calves by 47% in the ileum and 69% in the cecum relative to buffer-treated calves. Treatment with GLP-2 promotes intestinal growth in neonatal calves and reduces the detrimental effects of nitro-oxidative stress in the ileocecum of calves with coccidiosis.

Short communication: Effect of heat stress during the dry period on gene expression in mammary tissue and peripheral blood mononuclear cells.

Heat stress (HT) during the dry period compromises mammary gland development, decreases future milk production, and impairs the immune status of dairy cows. Our objective was to evaluate the effect of cooling HT cows during the dry period on gene expression of the mammary gland and peripheral blood mononuclear cells (PBMC). Cows were dried off 46 d before their expected calving and assigned to 2 treatments, HT or cooling (CL). Cows in the CL group were cooled with sprinklers and fans whereas HT cows were not. After parturition, all cows were housed in a freestall barn with cooling. The PBMC were isolated at dry-off and at -20, 2, and 20 d relative to calving from a subset of cows (HT, n=9; CL, n=10), and mammary biopsies were taken at the same intervals (HT, n=7; CL, n=6) for RNA extraction. Gene expression was assessed using a custom multiplex gene expression assay based on traditional reverse transcription-PCR. Genes involved in prolactin (PRL) signaling [PRL receptor long form, PRL receptor short form, suppressor of cytokine signaling (SOCS)2, SOCS3, IGF2, IGF binding protein 5, and cyclin D1], fatty acid metabolism (acetyl-CoA carboxylase α (ACACA) and lipoprotein lipase (LPL)], and IGF1 were evaluated in mammary tissue, and genes related to fatty acid metabolism [ACACA, fatty acid synthase (FASN), and LPL], cytokine production [IL6, IL8, and tumor necrosis factor (TNF)], and IGF1 were evaluated in PBMC. No differences were observed in PRL signaling or fatty acid metabolism gene expression in the mammary gland. In PBMC, HT cows had greater mRNA expression of IGF1 and TNF during the transition period relative to CL and upregulated IL8 and downregulated FASN mRNA expression at 2 d relative to calving. We conclude that cooling HT cows during the dry period alters expression of genes involved in cytokine production and lipid metabolism in PBMC.

Genomic imputation and evaluation using high-density Holstein genotypes.

Genomic evaluations for 161,341 Holsteins were computed by using 311,725 of 777,962 markers on the Illumina BovineHD Genotyping BeadChip (HD). Initial edits with 1,741 HD genotypes from 5 breeds revealed that 636,967 markers were usable but that half were redundant. Holstein genotypes were from 1,510 animals with HD markers, 82,358 animals with 45,187 (50K) markers, 1,797 animals with 8,031 (8K) markers, 20,177 animals with 6,836 (6K) markers, 52,270 animals with 2,683 (3K) markers, and 3,229 nongenotyped dams (0K) with >90% of haplotypes imputable because they had 4 or more genotyped progeny. The Holstein HD genotypes were from 1,142 US, Canadian, British, and Italian sires, 196 other sires, 138 cows in a US Department of Agriculture research herd (Beltsville, MD), and 34 other females. Percentages of correctly imputed genotypes were tested by applying the programs findhap and FImpute to a simulated chromosome for an earlier population that had only 1,112 animals with HD genotypes and none with 8K genotypes. For each chip, 1% of the genotypes were missing and 0.02% were incorrect initially. After imputation of missing markers with findhap, percentages of genotypes correct were 99.9% from HD, 99.0% from 50K, 94.6% from 6K, 90.5% from 3K, and 93.5% from 0K. With FImpute, 99.96% were correct from HD, 99.3% from 50K, 94.7% from 6K, 91.1% from 3K, and 95.1% from 0K genotypes. Accuracy for the 3K and 6K genotypes further improved by approximately 2 percentage points if imputed first to 50K and then to HD instead of imputing all genotypes directly to HD. Evaluations were tested by using imputed actual genotypes and August 2008 phenotypes to predict deregressed evaluations of US bulls proven after August 2008. For 28 traits tested, the estimated genomic reliability averaged 61.1% when using 311,725 markers vs. 60.7% when using 45,187 markers vs. 29.6% from the traditional parent average. Squared correlations with future data were slightly greater for 16 traits and slightly less for 12 with HD than with 50K evaluations. The observed 0.4 percentage point average increase in reliability was less favorable than the 0.9 expected from simulation but was similar to actual gains from other HD studies. The largest HD and 50K marker effects were often located at very similar positions. The single-breed evaluation tested here and previous single-breed or multibreed evaluations have not produced large gains. Increasing the number of HD genotypes used for imputation above 1,074 did not improve the reliability of Holstein genomic evaluations.

Heat stress abatement during the dry period influences metabolic gene expression and improves immune status in the transition period of dairy cows.

Heat stress (HT) and photoperiod affect milk production and immune status of dairy cows. The objective was to evaluate the effects of HT abatement prepartum under controlled photoperiod on hepatic metabolic gene expression and cellular immune function of periparturient Holstein cows (n=21). Cows were dried off 46 d before expected calving date and assigned to treatments by mature equivalent milk production. The treatments were 1) HT and 2) cooling (CL), both imposed during a photoperiod of 14L:10D. Rectal temperature was measured twice daily, whereas respiration rate was measured 3 times/wk at 1500 h during the entire dry period. After calving, cows were housed in a freestall barn with cooling, and milk yield was recorded daily up to 140 d in milk. Liver samples were taken at dry off, -20, 2, and 20 d relative to calving by biopsy. Under a similar schedule, neutrophil function was determined in blood of cows on HT (n=12) and CL (n=9). Blood samples were taken on -46, -32, -18, 0, 14, 28, and 42 d relative to calving for measurement of metabolites and were collected twice daily from -7 to 2 d relative to calving for prolactin (PRL) analysis. The HT cows had greater concentrations of PRL at 0 d relative to calving (150 vs. 93; SEM=11 ng/mL) and had higher afternoon rectal temperatures (39.4 vs. 39.0; SEM=0.04°C) and elevated respiration rates (78 vs. 56; SEM=2 breaths/min) during the prepartum period compared with CL cows. Relative to HT cows, CL cows had greater hepatic expression of PRL-R, SOCS-3, and CAV-1 mRNA. Neutrophil oxidative burst was greater in CL cows relative to HT cows at 2 d (61 vs. 42; SEM=6%) and at 20 d (62 vs. 49; SEM=5%) relative to calving, and phagocytosis was greater in CL cows at 20 d (47 vs. 33; SEM=4%) relative to calving compared with HT cows. Humoral response, as measured by IgG secretion against ovalbumin challenge, was greater for CL cows at -32 d (0.44 vs. 0.33; SEM=0.05 OD) and -21 d (0.60 vs. 0.50±0.04 OD) relative to calving compared with HT cows. These results suggest that HT abatement during the dry period improved innate and acquired immune status as measured by neutrophil function and immunoglobulin secretion against ovalbumin challenge, and altered hepatic gene expression related to PRL signaling in the periparturient period or subsequent lactation.

Characterization of glucagon-like peptide 2 pathway member expression in bovine gastrointestinal tract.

Glucagon-like peptide 2 (GLP-2), secreted by enteroendocrine cells, has several physiological effects on the intestine of monogastric species, including promotion of growth of intestinal epithelium, reduction of epithelial cell apoptosis, and enhancement of intestinal blood flow, nutrient absorption, and epithelial barrier function. The regulatory functions of GLP-2 in the ruminant gastrointestinal tract (GIT) have not been well studied. The objectives of this investigation were to characterize the mRNA expression of 4 members of the GLP-2 pathway throughout the bovine GIT, including (1) proglucagon (GCG), the parent peptide from which GLP-2 is derived through cleavage by prohormone convertase; (2) prohormone convertase (PCSK1); (3) GLP-2 receptor (GLP2R); and (4) dipeptidyl peptidase IV (DPP4), the enzyme that inactivates GLP-2. Gene expression was evaluated in rumen, reticulum, omasum, abomasum, duodenum, jejunum, ileum, cecum, and rectum collected at slaughter from prepubertal heifers, mature cows in early, mid, and late lactation, and nonlactating cows (n=3 per stage) by a gene expression profiling assay. In addition, mRNA expression of 14 genes involved in nutrient transport, enzyme activity, blood flow, apoptosis, and proliferation were evaluated in the 9 GIT tissues for their association with GCG and GLP2R mRNA expression. Immunohistochemistry was used to localize GLP2R protein in tissues of the lower GIT. Results indicated that mRNA expression of GCG, PCSK1, GLP2R, and DPP4 varies across the 9 GIT tissues, with greatest expression in small and large intestines, and generally nondetectable levels in forestomachs. Expression of DPP4 and GLP2R mRNA varied by developmental stage or lactational state in intestinal tissues. Expression of GCG or GLP2R mRNA was correlated with molecular markers of proliferation, apoptosis, blood flow, enzyme activity, and urea transport, depending on the tissue examined, which suggests a potential for involvement of GLP-2 in these physiological processes in the ruminant GIT. The GLP2R protein was expressed in intestinal crypts of the bovine GIT, which is consistent with the distribution in monogastric species. Our findings support a functional role of the GLP-2 pathway in bovine GIT and the potential for use of GLP-2 as a therapy to improve intestinal function and nutrient absorption in ruminants.

An improved algorithm of white matter hyperintensity detection in elderly adults.

Automated segmentation of the aging brain raises significant challenges because of the prevalence, extent, and heterogeneity of white matter hyperintensities. White matter hyperintensities can be frequently identified in magnetic resonance imaging (MRI) scans of older individuals and among those who have Alzheimer's disease. We propose OASIS-AD, a method for automatic segmentation of white matter hyperintensities in older adults using structural brain MRIs. OASIS-AD is an approach evolved from OASIS, which was developed for automatic lesion segmentation in multiple sclerosis. OASIS-AD is a major refinement of OASIS that takes into account the specific challenges raised by white matter hyperintensities in Alzheimer's disease. In particular, OASIS-AD combines three processing steps: 1) using an eroding procedure on the skull stripped mask; 2) adding a nearest neighbor feature construction approach; and 3) applying a Gaussian filter to refine segmentation results, creating a novel process for WMH detection in aging population. We show that OASIS-AD performs better than existing automatic white matter hyperintensity segmentation approaches.