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.

Adipose triglyceride lipase protein abundance and translocation to the lipid droplet increase during leptin-induced lipolysis in bovine adipocytes

Proper regulation of lipid metabolism is critical for preventing the development of metabolic diseases. It is clear that leptin plays a critical role in the regulation of energy homeostasis by regulating energy intake. However, leptin can also regulate energy homeostasis by inducing lipolysis in adipocytes, but it is unclear how the major lipases are involved in leptin-stimulated lipolysis. Therefore, the objectives of this study were to determine if (1) leptin acts directly to induce lipolysis in bovine adipocytes, (2) the potential lipases involved in leptin-induced lipolysis in bovine adipocytes, and (3) increases translocation of adipose triglyceride lipase (ATGL) and hormone sensitive lipase (HSL) during leptin-stimulated lipolysis in bovine stromal vascular cell-derived adipocytes. As hypothesized, leptin induced a lipolytic response (P = 0.02) in isolated adipocytes which was accompanied by an increase in phosphorylation of signal transducer and activator of transcription (STAT)3 (P = 0.03), a well-documented secondary messenger of leptin, and ATGL protein abundance (P < 0.01). Protein abundance of STAT3, perilipin, HSL, and phosphorylation of HSL by PKA and AMPK were not altered during leptin-stimulated lipolysis (P > 0.05). Immunostaining techniques were employed to determine the location of HSL and ATGL. Both lipases translocated to the lipid droplet after 2 h of exposure to isoproterenol (P < 0.02). However, only ATGL was translocated to the lipid droplet during leptin-stimulated lipolysis (P = 0.04), indicating ATGL may be the active lipase in leptin-stimulated lipolysis. In summary, leptin stimulates lipolysis in bovine adipocytes. The lack of phosphorylated HSL and translocation of HSL to the lipid droplet during leptin-stimulated lipolysis suggest minimal activity by PKA. Interestingly, leptin-stimulated lipolysis is accompanied by an increase in ATGL protein abundance and translocation to the lipid droplet, indicating its involvement in leptin-stimulated lipolysis either due to an increase in protein abundance or through a novel lipolytic cascade.

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.

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.

Effects of acute heat stress on lipid metabolism of bovine primary adipocytes

Heat stress (HS) affects numerous physiological processes including nutrient partitioning and lipid metabolism. Objectives of this study were to evaluate how acute HS affects lipid metabolism in subcutaneous adipose tissue of dairy cattle. Adipose tissue biopsies were performed on Holstein cows for bovine primary adipocyte isolation and cultured at either 42°C (HS) or 37°C (thermal neutral, TN). Adipocytes were incubated with increasing isoproterenol (ISO), and with increasing concentrations of insulin in the presence of ISO to evaluate changes in lipolysis. Incorporation of radioactive acetate into lipids was measured as an indicator of lipogenesis. Abundance and phosphorylation of several lipolytic and lipogenic proteins were also measured. Adipocytes exposed to HS had an elevated maximal response to ISO and were more sensitive to lipolytic stimulation by ISO compared with cells cultured at TN. Thermal treatment did not affect the antilipolytic effects of insulin in the presence of ISO. Lipogenesis measured as acetate incorporation was not altered by HS, but a temperature by insulin interaction was observed for the regulation of acetyl CoA carboxylase, such that the presence of insulin resulted in a reduction in phosphorylation of acetyl CoA carboxylase in adipocytes cultured at TN but not HS conditions. Results of this study demonstrate that acute HS has a direct effect on the regulation of lipolysis and the rate-limiting enzyme of lipogenesis in isolated bovine adipocytes.

The impact of 3 strategies for incorporating polled genetics into a dairy cattle breeding program on the overall herd genetic merit

Dehorning in cattle has been associated with behavioral, physiological, and neuroendocrine responses indicative of pain. Unaddressed, the pain associated with a routine production procedure could contribute to a negative public perception of livestock production practices. Alternative considerations of dehorning include the selection of polled cattle within herds, thereby avoiding pain and production loss. As polledness results from an autosomal dominant pattern of inheritance, genetic selection for polled cattle could reduce the prevalence of the horned trait. Herein we discuss 3 strategies to incorporate polled genetics into a cow herd and the estimated impact on the overall genetic merit of the herd. Furthermore, the availability and genetic merit of polled artificial insemination bulls in the United States is summarized. Both Holstein and Jersey dairy bulls registered with the National Association of Animal Breeders from December 2010 through April 2013 were queried. Polled bulls were identified as either being homozygous (PP) or heterozygous (Pp) and the average net merit (NM) predicted transmitting ability (PTA) of each sire group was calculated. The percentage of polled calves born each year over a 10-yr period was calculated for the following 3 scenarios: (A) various percentages of horned cows were randomly mated to Pp bulls, (B) various percentages of horned cows were preferentially mated to Pp bulls, and (C) horned cows were selectively mated to PP bulls, heterozygous cows to Pp bulls, and homozygous polled cows to horned bulls. Additionally, the change in NM PTA of the cow herd was calculated over the same period. The highest percentage of polled animals (87%) was achieved in scenario C. An evaluation of the herd NM PTA highlights the trade-offs associated with increasing polled genetics. Given the current genetic merit of horned and polled bulls, increasing the percentage of polled calves will decrease the NM PTA in Holstein, but may have minimal impact in Jersey herds. Decisions regarding selective breeding to increase polled genetics will need to be evaluated in the context of production objectives, cost of dehorning, and impact on overall genetic merit.

Regulation of feed efficiency in dairy cattle

Abstract

The efficiency with which dairy cows convert feed to milk impacts farm profitability and environmental stewardship. Thus, feed efficiency is a particularly critical trait to the dairy industry as feed costs increase. Feed efficiency is a complex trait influenced by many components, including milk output, nutrient intake, body weight and body weight change. Recent studies show that feed efficiency is low to moderately heritable, indicating it should respond to genetic selection. Progress is being made towards understanding the genetic regulation of feed efficiency, estimation of genomic breeding values and incorporation of feed efficiency into multi-trait selection programmes. However, underlying biological processes that contribute to variation in feed efficiency are poorly understood, particularly in lactating dairy cattle. Selection experiments in beef cattle, pigs, chickens and mice demonstrate that feed efficiency can be altered by selection. Heat production, activity or feeding behaviour, digestibility, muscle turnover and mitochondrial function are associated with variation in feed efficiency in many of these selection lines, but data from lactating dairy cows are currently lacking. A more detailed understanding of underlying biological processes and their relative importance to variation in feed efficiency is crucial for the development of novel strategies to further improve feed efficiency in dairy cattle. The objective of this review is to summarize our current understanding of the genetic and biological regulation of feed efficiency in lactating dairy cows. Insights from the evaluation of feed efficiency in beef cattle and other species are also considered.

Genetic parameters for energy balance, feed efficiency, and related traits in Holstein cattle

Objectives of the current study were to estimate genetic parameters in Holstein cows for energy balance (EB) and related traits including dry matter intake (DMI), body weight (BW), body condition score (BCS), energy-corrected milk (ECM) production, and gross feed efficiency (GFE), defined as the ratio of total ECM yield to total DMI over the first 150 d of lactation. Data were recorded for the first half of lactation on 227 and 175 cows in their first or later lactation, respectively. Random regression models were fitted to longitudinal data. Also, each trait was averaged over monthly intervals and analyzed by single and multivariate animal models. Heritability estimates ranged from 0.27 to 0.63, 0.12 to 0.62, 0.12 to 0.49, 0.63 to 0.72, and 0.49 to 0.53 for DMI, ECM yield, EB, BW, and BCS, respectively, averaged over monthly intervals. Daily heritability estimates ranged from 0.18 to 0.30, 0.10 to 0.26, 0.07 to 0.22, 0.43 to 0.67, and 0.25 to 0.38 for DMI, ECM yield, EB, BW, and BCS, respectively. Estimated heritability for GFE was 0.32. The genetic correlation of EB at 10d in milk (DIM) with EB at 150 DIM was -0.19, suggesting the genetic regulation of this trait differs by stage of lactation. Positive genetic correlations were found among DMI, ECM yield, and BW averaged over monthly intervals, whereas correlations of these traits with BCS depended upon stage of lactation. Total ECM yield for the lactation was positively correlated with DMI, but a negative genetic correlation between total ECM yield and EB was found. However, the genetic correlation between total ECM yield and EB in the first month of lactation was -0.02, indicating that total production is not genetically correlated with EB during the first month of lactation, when negative EB is most closely associated with diminished fitness. The genetic correlation between GFE and EB ranged from -0.73 to -0.99, indicating that selection for more efficient cows would favor a lower energy status. However, the genetic correlation between EB in the first month of lactation and GFE calculated from 75 to 150 DIM was not significant, indicating that the unfavorable correlation between GFE and EB in early lactation may be minimized with alternative definitions of efficiency. Thus, EB, GFE and related traits will likely respond to genetic selection in Holstein cows. However, the impact of selection for improved feed efficiency on EB must be carefully considered to avoid potential negative consequences of further reductions in EB at the onset of lactation.

Adipose tissue angiopoietin-like protein 4 messenger RNA changes with altered energy balance in lactating Holstein cows

Negative energy balance at the onset of lactation is unfavorably associated with fitness traits in high-producing dairy cows. Angiopoietin-like protein 4 (ANGPTL4) is an adipokine that has been associated with the regulation of lipid metabolism through the inhibition of lipoprotein lipase activity and regulation of lipolysis. Expression of ANGPTL4 messenger RNA (mRNA) increases during early lactation, but its regulation with changing energy status is currently unknown. Accordingly, the objective of this study was to determine whether ANGPTL4 mRNA abundance is responsive to declining energy balance induced by the transition from pregnancy to lactation, feed restriction, and GH administration in lactating dairy cows. The mRNA abundance of leptin, adiponectin, and adiponectin receptor 2 were also measured to compare adipokine mRNA profiles during changes in energy metabolism. Repeated adipose tissue biopsies were taken from different cows during transition from late pregnancy to lactation (n = 26), feed restriction (n = 19), and GH administration (n = 20). As expected, milk yield increased with the onset of lactation and GH administration (P < 0.01) but declined during feed restriction. Energy balance declined in each experiment, resulting in negative energy balance at the onset of lactation and after feed restriction. Abundance of ANGPTL4 mRNA expression increased 2- to 6-fold with declining energy balance in each experiment. Leptin mRNA declined with feed restriction, and adiponectin mRNA decreased with the onset of lactation. The consistency and magnitude of the increase in ANGPTL4 mRNA across multiple models of altered energy balance identifies it as an adipokine that is uniquely responsive to changes in energy balance in the lactating dairy cow.

Phosphorylation of perilipin is associated with indicators of lipolysis in Holstein cows

Perilipin is a regulatory protein that coats the lipid droplet in adipocytes. In the basal state, perilipin inhibits lipolysis by restricting access of hormone sensitive lipase (HSL) to the lipid droplet. In contrast, during stimulated lipolysis, phosphorylated perilipin interacts with HSL such that the catalytic activity of HSL on its acylglycerol substrate is enhanced. However, the regulation and function of perilipin in vivo has not been defined clearly across comparative animal models. Consequently, this study was undertaken to determine if changes in perilipin mRNA, protein, or phosphorylation state are associated with in vivo indicators of lipolysis in the dairy cow as a model of lipolysis induced by the marked metabolic demands of lactation. Semiquantitative western blotting and quantitative PCR were used to quantify total and phosphorylated HSL and perilipin in adipose tissue obtained from cows in early [5-14 days in milk (DIM), n=11] and mid (176-206 DIM, n=9) lactation. As expected, circulating NEFA and glycerol concentrations, and phosphorylated HSL were greater in early versus mid lactation, indicative of greater lipolytic activity in early lactation. Furthermore, phosphorylated, but not total perilipin abundance, was greater in early lactation when the metabolic demand for energy is greater than in mid lactation. Finally, the abundance of phosphorylated perilipin was positively correlated with circulating glycerol and NEFA concentrations during both early and mid lactation. Collectively, these data support the hypothesis that phosphorylated perilipin is a critical determinant of lipolytic activity stemming from the metabolic demands of lactation.

Relationships of a transforming growth factor-beta2 single nucleotide polymorphism and messenger ribonucleic acid abundance with bone and production traits in chickens

Osteoporosis is a serious problem for the laying hen industry with economic, production, and welfare consequences. Transforming growth factor-beta2 (TGFbeta2) has been implicated as an important factor in coupling bone resorption and formation in bone remodeling. The current study was designed to determine if TGFbeta2 was associated with variation in bone mineralization in chickens, using 2 complementary experimental approaches. First, an intronic single nucleotide polymorphism (SNP) present in TGFbeta2 was investigated in an F(2) population to determine its association with bone, growth, and egg traits of importance to the layer and broiler industries. The TGFbeta2 SNP was significantly associated (P < 0.05) with bone mineral density and content. However, these associations became nonsignificant when BW was included as a covariate in analyses. The TGFbeta2 SNP was also significantly associated (P < 0.05) with BW from 1 to 6 wk of age and egg production from 46 to 55 wk of age. To further explore the relationship between TGFbeta2 and bone strength, bone marrow TGFbeta2 mRNA abundance was compared between broiler and layer chickens at 15, 35, and 60 wk of age. Bone and egg traits were measured along with mRNA abundance at each age and found to differ significantly between lines. The TGFbeta2 mRNA abundance was approximately 4-fold greater in broiler compared with layer hens at 15 wk of age but was similar between lines at later ages. Thus, even though the TGFbeta2 SNP will likely not be an effective marker for improving bone strength independently of changes in BW, further research is warranted to investigate the relationship of TGFbeta2 mRNA abundance to bone strength in laying hens.

Postnatal changes in the expression of genes located in the callipyge region in sheep skeletal muscle

The expression of five genes surrounding the callipyge (CLPG) mutation was analysed in skeletal muscles from lambs at one prenatal and two postnatal ages that coincide with the onset and establishment of muscle hypertrophy. Genotype-specific changes in transcript abundance were detected for paternal allele-specific DLK1 and PEG11 (the official symbol of the latter is RTL1) and the maternal allele-specific MEG3, PEG11AS and MEG8 when the mutation was inherited in cis. There were differences in the temporal and muscle-specific effects on expression between the maternal allele-specific genes and paternal allele-specific genes. Maternal inheritance of the CLPG allele had a significant effect on the expression of MEG3 and MEG8 at prenatal and postnatal ages, whereas paternal inheritance of DLK1 and PEG11 only affected postnatal expression. Genotype-specific changes in PEG11AS expression were detected only in prenatal muscle. Maternal inheritance of the mutation caused similar changes in MEG3 and MEG8 expression in the semimembranosus, which undergoes hypertrophy, and the supraspinatus, which does not hypertrophy. Paternal inheritance of the mutation caused changes in PEG11 expression in both muscles, although the magnitude of expression in semimembranosus was more than 100-fold greater than in supraspinatus. DLK1 expression was upregulated in callipyge animals at both postnatal ages in the semimembranosus, but there was no effect of genotype on DLK1 expression in the supraspinatus at any age. Increased DLK1 expression was likely the primary cause of muscle hypertrophy, but a contribution of PEG11 to the phenotype cannot be ruled out based on gene expression.