D.M. Spurlock

Changes in skeletal muscle gene expression following clenbuterol administration

BackgroundBeta-adrenergic receptor agonists (BA) induce skeletal muscle hypertrophy, yet specific mechanisms that lead to this effect are not well understood. The objective of this research was to identify novel genes and physiological pathways that potentially facilitate BA induced skeletal muscle growth. The Affymetrix platform was utilized to identify gene expression changes in mouse skeletal muscle 24 hours and 10 days after administration of the BA clenbuterol.ResultsAdministration of clenbuterol stimulated anabolic activity, as indicated by decreased blood urea nitrogen (BUN; P < 0.01) and increased body weight gain (P < 0.05) 24 hours or 10 days, respectively, after initiation of clenbuterol treatment. A total of 22,605 probesets were evaluated with 52 probesets defined as differentially expressed based on a false discovery rate of 10%. Differential mRNA abundance of four of these genes was validated in an independent experiment by quantitative PCR. Functional characterization of differentially expressed genes revealed several categories that participate in biological processes important to skeletal muscle growth, including regulators of transcription and translation, mediators of cell-signalling pathways, and genes involved in polyamine metabolism.ConclusionGlobal evaluation of gene expression after administration of clenbuterol identified changes in gene expression and overrepresented functional categories of genes that may regulate BA-induced muscle hypertrophy. Changes in mRNA abundance of multiple genes associated with myogenic differentiation may indicate an important effect of BA on proliferation, differentiation, and/or recruitment of satellite cells into muscle fibers to promote muscle hypertrophy. Increased mRNA abundance of genes involved in the initiation of translation suggests that increased levels of protein synthesis often associated with BA administration may result from a general up-regulation of translational initiators. Additionally, numerous other genes and physiological pathways were identified that will be important targets for further investigations of the hypertrophic effect of BA on skeletal muscle.

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.

Coordination of lipid droplet-associated proteins during the transition period of Holstein dairy cows.

Dairy cows often experience negative energy balance with the onset of lactation, and severe or prolonged negative energy balance can contribute to declines in overall fitness. Energy stores, in the form of adipose tissue triacylglycerides, are mobilized during times of energy deficit, and recent research has implicated several proteins associated with the lipid droplet as lipolytic regulators. The objective of this study was to determine if these novel proteins associated with lipolytic regulation are altered with the changing metabolic demands of lactation. Weekly blood samples were collected from 26 Holstein cows from 21 d before expected parturition through 28 d postpartum, and again at 150 d postpartum. Serum nonesterified fatty acids, glycerol, and β-hydroxybutyrate were measured. Energy balance was calculated from daily feed intake and milk yield, weekly body weight, and monthly milk component measurements. Adipose tissue biopsies were taken 21 d before expected parturition (-21 d) and at 5, 21, and 150 d postpartum. Semiquantitative Western blotting was used to measure abundance of hormone-sensitive lipase (HSL), phosphorylated HSL, perilipin, phosphorylated perilipin (PPLIN), adipose triglyceride lipase (ATGL), and comparative gene identity-58 (CGI-58). Abundance of ATGL was less at 5 and 21 d in milk (DIM) compared with -21 and 150 DIM, even though cows were in negative energy balance and experiencing increased rates of lipolysis in early lactation. In contrast, phosphorylated HSL and PPLIN increased with increasing lipolysis immediately after parturition. Additionally, PPLIN was negatively correlated with milk yield at 5, 21, and 150 d postpartum, and negatively correlated with feed intake and energy balance at 21 d postpartum. This result is consistent with the hypothesis that phosphorylation of perilipin is responsive to signals for increased triaclyglyceride mobilization. Finally, a consistent negative correlation between abundance of perilipin and CGI-58 proteins was observed throughout the transition period. These results confirm that novel lipolytic proteins in adipose tissue are regulated at the level of protein abundance and phosphorylation during the periparturient period and into mid lactation.

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.

Harnessing the genetics of the modern dairy cow to continue improvements in feed efficiency

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.

Random Forests approach for identifying additive and epistatic single nucleotide polymorphisms associated with residual feed intake in dairy cattle

Feed efficiency is an economically important trait in the beef and dairy cattle industries. Residual feed intake (RFI) is a measure of partial efficiency that is independent of production level per unit of body weight. The objective of this study was to identify significant associations between single nucleotide polymorphism (SNP) markers and RFI in dairy cattle using the Random Forests (RF) algorithm. Genomic data included 42,275 SNP genotypes for 395 Holstein cows, whereas phenotypic measurements were daily RFI from 50 to 150 d postpartum. Residual feed intake was defined as the difference between an animals feed intake and the average intake of its cohort, after adjustment for year and season of calving, year and season of measurement, age at calving nested within parity, days in milk, milk yield, body weight, and body weight change. Random Forests is a widely used machine-learning algorithm that has been applied to classification and regression problems. By analyzing the tree structures produced within RF, the 25 most frequent pairwise SNP interactions were reported as possible epistatic interactions. The importance scores that are generated by RF take into account both main effects of variables and interactions between variables, and the most negative value of all importance scores can be used as the cutoff level for declaring SNP effects as significant. Ranking by importance scores, 188 SNP surpassed the threshold, among which 38 SNP were mapped to RFI quantitative trait loci (QTL) regions reported in a previous study in beef cattle, and 2 SNP were also detected by a genome-wide association study in beef cattle. The ratio of number of SNP located in RFI QTL to the total number of SNP in the top 188 SNP chosen by RF was significantly higher than in all 42,275 whole-genome markers. Pathway analysis indicated that many of the top 188 SNP are in genomic regions that contain annotated genes with biological functions that may influence RFI. Frequently occurring ancestor-descendant SNP pairs can be explored as possible epistatic effects for further study. The importance scores generated by RF can be used effectively to identify large additive or epistatic SNP and informative QTL. The consistency in results of our study and previous studies in beef cattle indicates that the genetic architecture of RFI in dairy cattle might be similar to that of beef cattle.

Genetic parameters between feed-intake-related traits and conformation in 2 separate dairy populations—the Netherlands and United States

To include feed-intake-related traits in the breeding goal, accurate estimates of genetic parameters of feed intake, and its correlations with other related traits (i.e., production, conformation) are required to compare different options. However, the correlations between feed intake and conformation traits can vary depending on the population. Therefore, the objective was to estimate genetic correlations between 6 feed-intake-related traits and 7 conformation traits within dairy cattle from 2 countries, the Netherlands (NL) and the United States (US). The feed-intake-related traits were dry matter intake (DMI), residual feed intake (RFI), milk energy output (MilkE), milk yield (MY), body weight (BW), and metabolic body weight (MBW). The conformation traits were stature (ST), chest width (CW), body depth (BD), angularity (ANG), rump angle (RA), rump width (RW), and body condition score (BCS). Feed intake data were available for 1,665 cows in NL and for 1,920 cows in US, from 83 nutritional experiments (48 in NL and 35 in US) conducted between 1991 and 2011 in NL and between 2007 and 2013 in US. Additional conformation records from relatives of the animals with DMI records were added to the database, giving a total of 37,241 cows in NL and 28,809 in US with conformation trait information. Genetic parameters were estimated using bivariate animal model analyses. The model included the following fixed effects for feed-intake-related traits: location by experiment-ration, age of cow at calving modeled with a second order polynomial by parity class, location by year-season, and days in milk, and these fixed effects for the conformation traits: herd by classification date, age of cow at classification, and lactation stage at classification. Both models included additive genetic and residual random effects. The highest estimated genetic correlations involving DMI were with CW in both countries (NL=0.45 and US=0.61), followed by ST (NL=0.33 and US=0.57), BD (NL=0.26 and US=0.49), and BCS (NL=0.24 and US=0.46). The MilkE and MY were moderately correlated with ANG in both countries (0.33 and 0.47 in NL, and 0.36 and 0.48 in US). Finally, BW was highly correlated with CW (0.77 in NL and 0.84 in US) and with BCS (0.83 in NL and 0.85 in US). Feed-intake-related traits were moderately to highly genetically correlated with conformation traits (ST, CW, BD, and BCS) in both countries, making them potentially useful as predictors of DMI.

Regulation of lipid droplet-associated proteins following growth hormone administration and feed restriction in lactating Holstein cows

Lipid metabolism plays a crucial role in the adaptation of dairy cows to periods of energy insufficiency. The objective of the current study was to determine if lipolytic proteins are consistently regulated when energy mobilization is stimulated by different factors. We evaluated 2 models of altered energy balance in mid-lactation Holstein cows, including feed restriction (FR) and administration of bovine growth hormone (GH), by quantifying the abundance and (or) phosphorylation of hormone-sensitive lipase (HSL), perilipin (PLIN), and adipose triglyceride lipase (ATGL). For GH administration, adipose tissue and blood samples were collected 4d before and 3 and 7d after administration of GH (n=20 cows). Similarly, adipose and blood samples were obtained 6d before and 1 and 4d after initiation of FR (n=18 cows). Estimated net energy balance decreased and nonesterified fatty acid concentration increased in both experimental models. Decreased ATGL and PLIN protein abundance was observed with GH administration and FR. Additionally, the abundance of phosphorylated HSLSer565 decreased in both models. Decreased abundance of phosphorylated PLIN was observed with GH administration, but not FR. Decreased ATGL protein abundance appears to be a consistent response to energy insufficiency in lactating cows, as this response was also described with negative energy balance at the onset of lactation. In contrast, the abundance of PLIN protein and phosphorylation of HSL using antibodies targeting serine residue 565 of HSL (HSLSer565) were altered in the current research, but not at the onset of lactation. Our findings demonstrate that lipolysis is altered through the regulation of multiple proteins, and that this regulation differs according to physiological state in lactating cows.

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.

Considerations when combining data from multiple nutrition experiments to estimate genetic parameters for feed efficiency

Prior to genomic selection on a trait, a reference population needs to be established to link marker genotypes with phenotypes. For costly and difficult-to-measure traits, international collaboration and sharing of data between disciplines may be necessary. Our aim was to characterize the combining of data from nutrition studies carried out under similar climate and management conditions to estimate genetic parameters for feed efficiency. Furthermore, we postulated that data from the experimental cohorts within these studies can be used to estimate the net energy of lactation (NE(L)) densities of diets, which can provide estimates of energy intakes for use in the calculation of the feed efficiency metric, residual feed intake (RFI), and potentially reduce the effect of variation in energy density of diets. Individual feed intakes and corresponding production and body measurements were obtained from 13 Midwestern nutrition experiments. Two measures of RFI were considered, RFI(Mcal) and RFI(kg), which involved the regression of NE(L )intake (Mcal/d) or dry matter intake (DMI; kg/d) on 3 expenditures: milk energy, energy gained or lost in body weight change, and energy for maintenance. In total, 677 records from 600 lactating cows between 50 and 275 d in milk were used. Cows were divided into 46 cohorts based on dietary or nondietary treatments as dictated by the nutrition experiments. The realized NE(L) densities of the diets (Mcal/kg of DMI) were estimated for each cohort by totaling the average daily energy used in the 3 expenditures for cohort members and dividing by the cohorts total average daily DMI. The NE(L) intake for each cow was then calculated by multiplying her DMI by her cohorts realized energy density. Mean energy density was 1.58 Mcal/kg. Heritability estimates for RFI(kg), and RFI(Mcal) in a single-trait animal model did not differ at 0.04 for both measures. Information about realized energy density could be useful in standardizing intake data from different climate conditions or management systems, as well as investigating potential genotype by diet interactions.