Introduction
Late pregnancy and lactogenesis are linked to comprehensive endocrine and metabolic changes. Homeorhetic adaptations are necessary to accomplish the nutrient drain towards the mammary gland. They are associated with decreasing peripheral insulin sensitivity (IS) in peripheral tissues other than the mammary gland in dairy cattle [1] just as in other mammalian species like humans [2]. Nevertheless, genetic selection for milk synthesis in dairy cattle resulted in a particular distinctive and prolonged negative energy balance (NEB) which is often associated with metabolic diseases like fatty liver [3].
Messenger molecules derived from adipose tissue (AT), i.e. adipokines, are important metabolic regulators and modulate IS [4]. This is well known for monogastric species in contrast to ruminants. The adipokines are related to energy metabolism (e.g. adiponectin (ADIPOQ) and leptin (LEP) and to inflammation (e.g. tumor-necrosis factor-a (TNF-a) and interleukin-6 (IL-6)) [5]. LEP and ADIPOQ increase IS whereby LEP is positively and ADIPOQ is negatively linked with body fat content [6,7]. LEP mainly exerts its effects through the leptin receptor (LEPR) isoform b (LEPRB) [8] and ADIPOQ via ADIPOQ receptor 1 (ADIPOR1) and ADIPOQ receptor 2 (ADIPOR2), as observed in humans [9]. Adiponectin reduces gluconeogenesis and decreases hepatic glucose release. Beta-oxidation of fatty acids is increased by ADIPOQ in liver and skeletal muscle. In the latter glucose uptake is also improved [7,9]. LEP is negatively associated with energy intake. It stimulates beta-oxidation and reduces blood glucose and blood lipids [6]. Activation of inflammatory pathways in AT and release of inflammatory cytokines from AT can affect IS directly at the level of insulin sensing, e.g. reduced phosphorylation of insulin receptor and insulin receptor substrate-1 lead to decreased glucose uptake as shown for TNF-a [10] and IL-6 [11].
The transcription factor peroxisome proliferator-activated receptor gamma (PPARc) is a central regulator of genes related to insulin sensitivity, including ADIPOQ, IL-6 and TNF-a [12,13]. In dairy cattle it was shown that PPARc gene expression is linked to lipogenic gene expression in the mammary gland [14] and that the PPARc agonist 2,4-thiazolidinedione increases plasma glucose [15]. In cattle, two isoforms of PPARc (PPARc1 and PPARc2) sharing 90% identity are described [16]. In general, PPARc2 is the predominant isoform in AT; PPARc plays an important role in adipocyte differentiation and improves IS e.g. by increasing glucose uptake and by the induction of genes related to fatty acid metabolism and fatty acid uptake. In addition, the isoform PPARc2 is needed for adiponectin gene transcription [17].
Fat from different areas of the body displays distinct structural and functional properties including adipokine secretion. In humans and rodents, visceral (v.c.) AT feature higher lipogenic and lipolytic activities and produce more pro-inflammatory cytokines, while subcutaneous (s.c.) AT are the main source of LEP and ADIPOQ, and thus have disparate roles, amongst others in IS as reviewed in [18]. So far, studies in dairy cows were mostly limited to s.c. fat which can readily be biopsied.
Linoleic acid and their derivatives like conjugated linoleic acids (CLA) are ligands for PPARs [19,20]. In 3T3L1 adipocytes trans- 10, cis-12 CLA decreases ADIPOQ in both PPARc dependent and independent mechanisms in vitro [21]. CLA is also a ligand for free fatty acid receptor 1 and activation of this receptor increases glucose stimulated insulin secretion in mice [22]. From rodent studies using high doses (up to 1% of the diet) of trans-10, cis-12 CLA, dramatic loss of AT, inflammation, insulin resistance and hepatic steatosis are reported [23]; however, when using lower doses, the effects were largely limited to reduced AT lipogenesis. In dairy cattle, CLA may act energy-sparing as the trans-10, cis-12 CLA isomer is related to milk fat depression whereas whole body fat seems largely unaffected during both short and long-term treatments with this CLA isomer [24]. Milk fat depression may improve NEB in dairy cows but the effects are not conclusive yet as discussed by Ho¨tger et al. [25]. This group has shown that the main isomers (trans-10, cis-12 CLA, cis-9, trans-11 CLA) decrease endogenous glucose production and milk fat synthesis but increase plasma glucose, whereby insulin secretion in cows during their second lactation was not affected.
The understanding of the expression and regulation of IS related genes in relevant tissues involved in energy metabolism in dairy cattle may help to develop strategies to reduce the incidence for metabolic diseases in future. Therefore we aimed to characterise the time course of the expression of IS related genes during lactation in different AT of v.c. and s.c. localisation in addition to muscle, liver and mammary gland. Based on the potential of CLA to spare energy in addition to affect glucose metabolism and possibly IS, we supplemented primiparous and pluriparous Holstein cows’ diet with a 1:1 mixture of trans-10, cis- 12 CLA and cis-9, trans-11 CLA which is registered as feed supplement in the EU. The main focus aims were (1) to characterise the temporal expression pattern in s.c. AT and liver biopsies from pluriparous cows from late pregnancy to late lactation, (2) to compare the mRNA expression between different s.c. and v.c. AT depots and in other metabolically relevant tissues available through timed slaughter of primiparous cows and (3) to test for potential CLA effects in cows of both parity groups. The target genes we selected herein are related to IS, as substantiated by a vast body of literature, and we therefore aimed to test not only tissue mRNA expression but also to estimate systemic IS using the ‘‘Revised Quantitative Insulin Sensitivity Check Index’’ (RQUICKI) that was suggested for use in dairy cows [26]. We aimed to identify the AT depots that are mostly affected throughout lactation and to define the predominant genes out of the analysed panel of genes involved in the regulation IS throughout lactation.
Materials and Methods
All animals were housed in free stall barns at the experimental station of the Friedrich Loeffler Institute, Federal Research Institute for Animal Health, Braunschweig, Germany
Ethics statement The experiments were approved by the competent authority, the lower saxony state office for consumer protection and food safety (LAVES, file no. 33.11.42502-04-071/07, Oldenburg, Germany). The regulations of the German Animal Welfare Act (TierSchG) in its respective edition were met.
Animals, diets and treatments
Trial 1. German Holstein pluriparous cows (n = 21) were fed a partial mixed ration consisting of 38% grass silage, 25% corn silage and 37% concentrate on a dry matter basis for ad libitum intake over the entire experimental period. Diets were formulated according to the recommendations of the German Society of Nutrition Physiology [27]. More details about the composition of the diet, the experiment, and the outcomes in terms of performance (including the calculated energy balance), metabolite concentrations and milk FFA profile from this experiment are described by Pappritz et al. [28]. On day 1 in milk (DIM) 1 the animals were allocated to 2 groups: a CLA group (n = 11) and a control group (n = 10). During the supplementation period (from 1 to 182 DIM), the animals received 4 kg of additional concentrate containing the fat supplements. The animals of the CLA group received 100 g/d encapsulated CLA (Lutrell Pure, BASF, Ludwigshafen, Germany). The control group received 100 g/d of a control fat supplement (Silafat, BASF) in which CLA were substituted by stearic acid to form an isoenergetic control diet using a fatty acid with the same number of carbon atoms as in CLA. In the CLA group, the animals consumed 7.6 g/day each of the trans-10, cis-12 and cis-9, trans-11 CLA isomers. Liver and s.c. AT from tail head were biopsied on d 21 before calving and on DIM 1, 21, 70, 105 182 196, 224, and 252. The biopsy procedures were previously described in detail by Saremi et al. [29]. For gene expression analyses, samples were restricted to the pluriparous cows of the control group (all sampling dates) and of the CLA group (d 221, 21, 105, 196, and 252) as the focus of the study was related to changes throughout lactation and the capacity of handling the samples was limited. Jugular blood samples were collected on d –21, –14, –7, 1, 7, 14, 21, 35, 49, 70, 105, 140, 182, 189, 196, 210, 224, 238, and 252 relative to calving.
Trial 2. German Holstein primiparous cows (n = 25, average age at calving was 23 months) were fed a partial mixed ration (ad libitum) according to the recommendations of the German Society of Nutrition Physiology [27]. The ration consisted of 63% silage (60% corn silage, 40% grass silage) and 37% concentrate on a dry matter basis. Five cows were slaughtered on 1 DIM and the remaining animals (n = 20) were randomly allocated to either receive 100 g/d of the control fat supplement (Silafat, BASF) or of CLA (Lutrell Pure, BASF) starting from 1 DIM. All experimental details and main outcomes in terms of performance, adipose depot and organ weights, body composition, body fat mobilization, protein accretion, and energy utilization were reported by von Soosten et al. [30,31]. Five cows per group were slaughtered at 42 and 105 DIM. Samples from 3 v.c. AT (omental, mesenteric, and retroperitoneal), 3 s.c. AT (tail head, withers, and sternum), liver, M. semitendinosus, pancreas, and mammary gland parenchyma were taken about 25 min after stunning. Blood samples of about 40 mL were collected from V. jugularis after the morning milking on d –21, –14, –7, –3, 1 (n = 25), 7, 14, 28, 42 (n = 20), and 105 (n = 10) relative to parturition.
The tissue samples from both trials were snap-frozen in liquid nitrogen and stored at –80uC for RNA extraction and analysis. The blood samples were collected in different tubes to obtain sodium fluoride-EDTA plasma, heparin plasma or serum. Plasma samples were stored at –20uC and sera at –80uC until analysed.
Relative quantification by real-time PCR
Preparation of samples including RNA extraction and cDNA synthesis together with real-time PCR and reference genes characteristics, selection, and measurements are described in detail by Saremi et al. [29,32]. Relative quantification of the target genes was performed as summarized in Table 1. Amplicon standard curve was used except for PPARc2 in Trial 1 and 2 and LEPRB and TNF-a in Trial 1 where cDNA standards curve was used. The standard curve was used to correct data based on PCR efficiency per run. Two mL cDNA (diluted 1 to 4) as template and 5 mL SYBR Green JumpStart Taq Readymix (Sigma-Aldrich, Steinheim, Germany) in a total volume of 10 mL were run in an Mx3000P real-time PCR cycler (Stratagene, Amsterdam, Netherlands and Agilent, Santa Clara, CA). Based on biological (abundance) and/or technical reasons (total RNA in use for cDNA synthesis), not all genes were quantified in each tissue in Trial 1 (ADIPOQ and PPARc2 in liver, IL-6 and LEPRB in s.c. AT were not evaluated) and in Trial 2 (ADIPOQ in liver and muscle, PPARc in muscle and mammary gland, LEP in liver were not evaluated).
Reference gene stability and data analysis
In each individual case, comparing treatment effects within a single tissue or comparing different tissues, an adequate panel of reference genes was selected according to Saremi et al. [29,32].
Trial 1. The reference genes used were low density lipoprotein receptor-related protein 10 (LRP10), glyceraldehyde-phosphatedehydrogenase (GAPDH), and RNA Polymerase II (POLR2A) for s.c. AT and LRP10, POLR2A, and eukariotic translation initiation factor 3, subunit K (EIF3K) for liver tissue according to Saremi et al. [29]. The geometric mean of the reference gene abundances was used for normalisation. Data are presented as ratio of the copy numbers of genes of interest and the geometric mean of the corresponding reference genes.
Trial 2. The quantification cycle values were imported in qBASEplus version 2.0 (Biogazelle, Ghent, Belgium) and all subsequent calculations and data quality controls were done based on this software [33]. Reference genes were tested per tissue and the lowest V value obtained for the number of reference genes needed for normalisation as described by Saremi et al. [32]. Briefly, for normalisation of the efficiency corrected data the reference genes EIF3K, LRP10, POLR2A, Emerin, Marvel domain containing 1 and Hippocalcin-like 1 were selected and used, depending on the tissues being compared [32] were used. The data used for statistics was generated by geNormplus as a part of qBASEplus version 2.0 (Biogazelle, Ghent, Belgium).
Analysis of non esterified fatty acids, ß-hydroxy butyric acid, LEP, and insulin in blood samples
Glucose was measured in sodium fluoride-EDTA plasma, for quantification of insulin, non esterified fatty acids (NEFA), and ß- hydroxy butyric acid (BHBA) heparinised plasma was used. The plasma concentrations of glucose, NEFA, and BHBA were determined by an automatic analyser (Cobas Mira Plus System from Roche Diagnostica Ltd, Basel, Switzerland) using commercial test kits (Glucose: Glucose Hexokinase Fluid 5+1, MTI Diagnostics GmbH, Idstein, Germany; NEFA: HR(2) R1+R2 Set, WAKO Chemicals GmbH, Neuss, Germany; BHBA: RANBUT, RB 1008, Randox Laboratories GmbH, Wu¨lfrath, Germany). LEP concentrations were determined in serum by ELISA [34]. The intra and inter assay coefficients of variations were 6.3% and 13.9%, respectively, the limit of detection was 0.3 ng/mL. The concentrations of insulin in plasma were measured using a commercially available double antibody radioimmunoassay (DSL-1600, Diagnostic Systems Laboratories, Inc., TX). The intra and inter assay coefficients of variation were 6.3% and 8.8%, respectively
Calculations and statistical analyses
The RQUICKI was calculated based on the concentrations of plasma glucose, NEFA, and insulin [35] as previously used to estimate IS in dairy cows [26] by the following equation:
All statistical analyses were performed using SPSS (version 20, SPSS Inc., Chicago, IL). For Trial 1, the mixed model procedure was used. Treatment (control and CLA) was considered as fixed factor, sampling dates (time) as repeated effect, and the respective interaction was included in the model. For serum LEP, but not for mRNA analyses which were limited to the cows from Trial 1, parity (cows or heifers) was considered as fixed effect together with its interaction with treatment. For each time series the prepartum value was considered as a covariate in the model to check for the effect of possible differences between treatment groups prepartum. The covariate was never significant. The covariance structures compound symmetry, unstructured and heterogeneous first-order autoregressive were tested; latter was used for serum data while compound symmetry was used for mRNA data. Bonferroni correction was used for correcting multiple comparisons.
For Trial 2 in which the repeated design for tissue samples was not applicable, all data were tested for homogeneity of variances (P#0.1). Accordingly the general linear model with the fixed effects of treatment, dates, and the respective interaction, or non parametric tests (Kruskal-Wallis or Mann-Whitney) were applied (P#0.05). For correlation analyses Pearsons product-moment coefficients were calculated, presentation of the coefficients of correlations is limited to significant correlations (P#0.05). All data are presented as mean 6 SEM. For both trials P#0.05 was considered significant.
Results
Trial 1 During the transition from late pregnancy to early lactation (transition period) in pluriparous cows, the mRNA abundance of the target genes was decreased in the majority of the cases in s.c. AT (Figures 1, 2, and 3). In detail, the mRNA abundance for LEP, LEPR, ADIPOQ, ADIPOR1, ADIPOR2, PPARc, and PPARc2 was decreased at that time in s.c. AT. The mRNAs for PPARc, PPARc2 and LEPR in s.c. AT remained lower throughout the analysed time period up to 252 DIM as compared to ante partum (a.p.). Similarly, the amount of PPARc and IL-6 mRNA in liver was reduced. In contrast, the quantity of LEPR and LEPRB as well as ADIPOR2 mRNA in liver increased during the transition period. The abundance of LEPR and LEPRB mRNA in liver was higher up to day 196 DIM than a.p., and decreased thereafter to values comparable to the situation a.p., but remained higher at day 252 in case of LEPR. The mRNA quantity for ADIPOR2 in liver was higher at 21 DIM vs. 221 DIM and decreased to the end of the analysed period. No peripartal and lactation related differences were observed for TNF-a in liver and s.c. AT. The pattern of serum LEP followed roughly the mRNA abundance throughout the transition period but these changes were not significant (Figure 2).
The coefficients of correlation calculated between the mRNA abundance of the target genes and blood glucose, insulin, NEFA and BHBA, RQICKI as well as with back fat thickness, body condition score and live weight (Table 2) ranged between r = 0.440 (LEP mRNA in s.c. AT and insulin) and r = 0.180 (LEPRB mRNA in liver and NEFA). Most correlations were related to ADIPOR1 mRNA in s.c. AT (7 out of 8 variables) which were negative for NEFA and BHB but there was no relation with back fat thickness. In contrast, significant correlations of ADIPOR2 mRNA in s.c. AT were limited to body condition score (BCS) and liver weight, whereas ADIPOR2 mRNA in liver had more correlations (6 out of 8 variables). These correlations were negative for glucose and insulin and liver weight but positive for NEFA and BHB. Serum LEP correlated positively with back fat thickness, comparable with LEP mRNA in s.c. AT (r = 0.225 vs. r = 0.245). LEPR in s.c. AT was positively associated with back fat thickness, but LEPR and LEPRB mRNA in liver were negatively correlated with back fat thickness. PPARc mRNA in liver was negatively correlated with NEFA, BHB, insulin and back fat thickness but positively associated with the RQUICKI index. Supplementation with CLA only sporadically increased or decreased the expression of few genes at discrete time points in liver and s.c. AT; however, no consistent changes were observed. Treatment with CLA increased the mRNA abundance of LEPRB and TNF-a mRNA abundance in liver, whereas PPARc mRNA was lower in CLA cows than in the controls at defined time points throughout lactation (Figure 3). Lower values than in the control group were also observed at 196 DIM i.e. beyond the supplementation period in case of PPARc2 in s.c. AT of CLA cows (Figure 3). When testing CLA effects on insulin, glucose, and NEFA, and RQUICKI values, CLA effects were evident: the insulin concentrations were greater and RQUICKI data were lower (8% or 17% during supplementation and depletion period, respectively) than in the control group (Figure 4).
Trial 2
Comparing the temporal changes of the different target genes across the first 105 days of lactation in different adipose depots, liver, muscle, and mammary gland of primiparous cows (Table 3 and 4), the most changes with time were observed for ADIPOR1. In mammary gland and all fat depots except omental AT, the receptor mRNA abundance increased with DIM. In contrast, the amount of this mRNA was reduced in muscle at d 42 and d 105, and at d 42 in liver tissue when compared against d 1. The mRNA abundance of ADIPOR2 increased with DIM in retroperitoneal AT; in mammary gland the values were higher at d 105 compared to d 42. The ADIPOQ mRNA abundance increased in two v.c. AT (omental and retroperitoneal) with DIM. LEPRB mRNA abundance was decreased with DIM in all v.c. AT and in the mammary gland. No changes were found for LEP mRNA in primiparous cows during lactation. The mRNA abundance for PPARc and PPARc2 changed in some tissues but not consistently. For PPARc an increase over time was detected in omental AT, and for PPARc2 in retroperitoneal AT. In addition, an increase was observed in liver for both transcripts. The amount of TNF-a mRNA increased with DIM in all s.c. AT and in liver; in mammary gland; in retroperitoneal AT, a decrease with DIM was observed. The mRNA of IL-6 was decreased in s.c. AT from tail head and withers but increased in liver over the time.
Summarising the associations between blood metabolites, body fat and liver mass with ADIPOQ mRNA abundance in different AT depots, negative correlations were observed between ADIPPOQ mRNA in the mesenteric AT depot and omental and mesenteric fat mass (r = 20.433 and r = 20.484, respectively) as well as with total fat mass (i.e. mass of s.c. and all v.c. depots taken together; r = 20.432). In addition, ADIPOQ mRNA in retroperitoneal AT was negatively correlated with NEFA (r = 20.428). Correlation analysis between the mRNAs of ADIPOR1 and ADIPOR2 in different AT depots and adipose tissue masses as well as liver mass revealed negative correlations between ADIPOR1 in mesenteric AT and AT masses (omental: r = 2 0.462, mesenteric: r = 20.461, retroperitoneal: r = 20.467, mass of s.c. fat and all v.c. depots: r = 20.480), but a positive correlation with liver mass was observed for ADIPOR1 in two v.c. depots (mesenteric: r = 0.407 and retroperitoneal: r = 0.456) and in the s.c. depot from sternum (r = 0.633). For ADIPOR2 correlations were limited to the s.c. depots, mainly the one from the withers area: positive associations were established with omental (r = 0.479), retroperitoneal (r = 0.491), total s.c. AT mass and total AT mass of all depots considered (r = 0.426 and r = 0.494).
Supplementation with CLA only had few and inconsistent effects on gene expression in primiparous cows (Table 3 and 4). The mRNA abundance of TNF-a mRNA was reduced by CLA at 42 DIM in mesenteric AT, and the ADIPOQ and ADIPOR2 mRNA abundance were decreased at 105 DIM in omental and retroperitoneal AT, respectively. The mRNA abundance of PPARc2 was reduced by CLA only in the mammary gland at both 42 and 105 DIM.
Differences in gene expression of target genes between individual adipose tissue depots and in other tissues
The mRNA abundance of the target genes was compared between different AT from Trial 2 irrespective of treatment and DIM (Figure 5). Differences between the different AT depots were seen for all target mRNAs except ADIPOR2. In addition, the mRNA abundance of most target genes was different when comparing only within v.c. depots or within s.c. depots.
By the real-time PCR protocol applied, ADIPOQ mRNA was detectable in the mammary gland only in one third of the samples. The abundance was about 900 times lower than in AT. For ADIPOR1, a 2 fold lower abundance was observed in the liver than in the average of all AT whereas in muscle the values were 3 fold higher than in AT. The abundance of ADIPOR2 mRNA in muscle was higher than in liver. In both tissues the abundance was mostly in the same range as in the different AT depots although liver and muscle had lower and higher ADIPOR2 mRNA abundance in comparison to retroperitoneal and mesenteric AT, respectively. Compared to AT, there was about 100 fold less PPARc2 mRNA in liver, 79 fold less in muscle and 17 times less in mammary gland. PPARc mRNA abundance in liver was 109 fold lower in comparison to the mean of AT. LEP mRNA abundance was low in muscle and mammary gland in comparison to the mean of fat depots (171 and 353 fold lower, respectively). In liver, 34 times more LEPRB mRNA abundance than in AT was observed. About 5 fold more TNF-a mRNA was observed in liver than in AT, but IL-6 mRNA abundance was 8.3 times lower in liver than in AT.
Discussion
Gene expression of the ADIPOQ and LEP system during lactation
In general, the adipokines LEP and ADIPOQ are considered as positive regulators of IS as demonstrated in monogastrics [4]. Skeletal muscle lipid stores are linked to insulin secretion and peripheral insulin sensitivity. LEP improves insulin sensitivity by increasing the activity of AMP-activated protein kinase (AMPK) in muscle which reduces the activity of acetyl-CoA carboxylase, therewith reducing lipid stores by increasing fatty acid oxidation in skeletal muscle [36]. On the other hand it should be taken into account that LEP also increases pro-inflammatory protein secretion like TNF-a by macrophages [37] but reduces proinflammatory protein secretion indirectly by activation of AMPK e.g. in adipose tissue [38]. In contrast to LEP, ADIPOQ acts antiinflammatory e.g. by reducing TNF-a secretion. A main target of ADIPOQ is AMPK. Its activation in liver is associated with reduced gluconeogenesis by the reduction of phosphoenolpyruvate carboxykinase and glucose-6-phosphatase expression and increased fatty acid oxidation by inactivation of acetyl-CoA carboxylase [39]. LEP is positively and ADIPOQ is negatively linked with body fat content [6,7]. Effects of LEP are mainly mediated by LEPRB [8] and ADIPOQ via ADIPOR1 and ADIPOR2, as observed in humans [9]. As confirmed herein, the blood concentrations of LEP in dairy cows are low in early lactation [40] and thus may affect the secretion of proinflammatory proteins as discussed above. Correlation analyses revealed deviant regulation of LEPR and LEPRB mRNA between back fat thickness and BCS in s.c. AT (positively correlated) vs. liver (negatively correlated). In liver LEPR mRNA was also negatively correlated with glucose and insulin but positively with NEFA and BHB. This might be an indication of the importance of LEP in regulating liver lipid content and gluconeogenesis [41]. In type 2 diabetic mice LEP improves hepatic insulin sensitivity by reducing gluconeogenesis which should be associated with the negative correlations between LEPR mRNA vs. glucose and insulin. On the other hand, possibly more relevant for dairy cows, LEP increases fatty acid oxidation in liver and reduces hepatic lipotoxicity [41], and decreased LEP at parturition might contribute to the increased vulnerability to hepatic steatosis. The observed upregulation of the receptor mRNA for LEPR and LEPRB in liver as well as the positive correlation with NEFA and in case of LEPR mRNA also with BHBA indicates a counterbalancing of LEP effects at the receptor level in this situation.
For ADIPOQ less information is available in cattle but decreased concentrations were documented using semiquantitative Western blots [42,43] and a bovine specific ELISA [43].
Western blots [42,43] and a bovine specific ELISA [43]. Koltes and Spurlock [44] observed a decrease of the ADIPOQ mRNA abundance in s.c. AT from tail head throughout the transition period. This result is comparable to our observation in Trial 1 in pluriparous cows and may affect signalling by ADIPOQ receptors which were down regulated in s.c. fat after parturition. For the LEPR mRNA in s.c. AT, our present results are not in line with data published by our group and others [45,46]. In these two studies, an increase in LEPRB and LEPR mRNA, respectively was observed in s.c. AT and explained by Thorn et al. [46] at least partly by the low insulin concentrations during early lactation. As expected, we also observed low insulin concentrations after calving.
In humans, ADIPOQ and its receptors are inversely related to body fat mass; both receptors and serum ADIPOQ are downregulated by insulin in lean men [9,47]. The reduction in mRNA abundance of the ADIPOQ system we observed during the transition period in s.c. AT of pluriparous cows in Trial 1 therefore indicates an uncoupling between energy balance and the ADIPOQ system as they were not regulated contrary to the energy status and insulin values as shown in humans and mice [9,47]. ADIPOQ mRNA was positively correlated with back fat thickness and BCS, as well as ADIPOR1/2 mRNA in s.c. AT thus emphasising the differential regulation in cows as compared to humans or mice at least at the mRNA level. Based on the number of correlations between the analysed parameters and ADIPOR1 mRNA in s.c. AT of pluriparous cows this receptor might be more relevant for paracrine/autocrine ADIPOQ effects in this AT depot than ADIPOR2; the latter seems to be closer linked to liver of pluriparous dairy cows since ADIPOR1 was positively but ADIPOR2 mRNA was negatively correlated with the mass of this organ. ADIPOR2 mRNA is more abundant than ADIPOR1 mRNA in liver of cattle [48], which may underline its importance in fatty acid oxidation by PPARa signalling [49]. ADIPOQ is associated with lipogenesis and differentiation during adipogenesis at least in murine 3T3-L1 adipocytes [50] and we observed a positive correlation with back fat thickness and BCS in pluriparous cows. The decreased ADIPOQ mRNA abundance in AT during early lactation might therefore be related to reduced lipogenesis and adipogenesis during NEB in discrete AT depots. The regulation of both receptors in liver differs from s.c. AT as the ADIPOR1 mRNA abundance was stable and the one for ADIPOR2 increased in liver during the transition period. Differential regulation of both receptors in these two tissues is underlined by different correlation patterns of both receptor mRNAs in s.c. AT and liver. More correlations with the different variables (e.g. negative correlations with NEFA and BHBA) were found for ADIPOR1 in s.c. AT and for ADIPOR2 (e.g. negative correlation with glucose and insulin, positive correlation with NEFA) in liver. Signalling by ADIPOR2 in liver during negative energy balance may underline a positive link to fatty acid oxidation and energy metabolism in cattle as suggested for mice [49]. Based on ADIPOQ signal transduction in liver as discovered in a mouse model with targeted disruption of the ADIPORs [49], reduced ADIPOQ serum concentrations during the transition period may enhance gluconeogenesis. The mRNA of ADIPOR2 increased after parturition in liver of pluriparous cows in this study. Similar observations were reported by Loor et al. [51]. In their study the expression of the receptor mRNA was further increased by feed restriction before calving. In addition, growth hormone (GH) might be related to this observation: GH knowingly peaks around parturition in dairy cattle [52] and ADIPOR2 but not ADIPOR1 was shown to be upregulated by GH in liver of mice [53]. Signalling by ADIPOR2 is associated with less oxidative stress in liver [49], and oxidative stress is increased by NEFA as shown in human hepatocytes [54]. During NEB in early lactation, NEFA concentrations increase [3]. Based on the fact that targeted ADIPOR2 disruption in liver of mice increased oxidative stress markers [49], upregulation of ADIPOR2 may help to mitigate oxidative stress in liver of dairy cows, that is increased due to lipid mobilization in the periparturient period [55].
In the primiparous cows from Trial 2, most of the temporal changes occurred for ADIPOR1 mRNA in several adipose depots, in particular in the s.c. AT depots in which the abundance increased with DIM. As mentioned earlier, this might also be related to the effects of ADIPOQ on adipogenesis and lipogenesis [50] which are down-regulated during NEB with increased lipolysis [56]. As discussed by Dagon et al. [57] and Rossmeisl et al. [58], AMPK is an important mediator of hormonal and nutritional effects on AT and is thus related to the control of body fat mass. Its activation involves the down-regulation of lipogenic and adipogenic genes like PPARc. Both ADIPORs are linked to AMPK activation but ADIPOR1 to a higher extent than ADIPOR2, being closer associated with PPARa signalling, at least in liver [49]. Therefore, an increase in abundance of ADIPOR1 until peak lactation might be related to an increase in AMPK activity, controlling lipolysis tightly in retroperitoneal AT and in different s.c. AT of primiparous cows. Based on a recent study, the energy sensor AMPK may inhibit lipolysis and hormone sensitive lipase transport to the lipid droplet [59], which is in line with our observation of increasing ADIPOR1 mRNA abundance in AT. Therefore, AMPK may protect the adipocyte against energy depletion as discussed by Gauthier et al. [60]. Antilipolytic effects of AMPK were also addressed for dairy cows based on increasing phosphorylation of AMPKa1 after parturition [61]. Summarising, increasing ADIPOR1 mRNA throughout lactation might be associated with reduced lipolysis in dairy cows by increasing AMPK activity which would be in line with positive effects of ADIPOQ on lipid accumulation in AT [50] and lowest ADIPOQ concentrations after parturition in dairy cows [43]. Interestingly, in the primiparous cows we observed negative correlations between mesenteric ADIPOQ mRNA and the omental and mesenteric AT mass comparable to ADIPOR1 mRNA in this depot which was also negatively correlated with retroperitoneal AT mass as well as total AT mass. This observation could be related to the stage of lactation (sampling dates: DIM 1, 42, 105) and negative ADIPOQ effects on gluconeogenesis [39] as well as the cows’ performance level. AT from cows at DIM 100 of two families with differences in their capability to store body fat and to secrete milk revealed positive correlations between mesenteric ADIPOQ mRNA and omental as well as mesenteric AT mass [48]. In humans, mesenteric AT has an important role in regulating IS [62], and possibly also in cows. High AT mass of this depot might be unfavourable by its secretion products on gluconeogenesis and hepatic IS in early lactation of high yielding dairy cows.
Gene expression of PPARc and PPARc2 during lactation
The transcription factor PPARc is a key regulator of IS and adipogenesis [63]. Targeted deletion of PPARc in murine liver decreased hepatic triglyceride content and impaired systemic IS [64]. Therefore, regulation of PPARc around calving may affect the insulin responsiveness of tissues. Like in other species, two PPARc isoforms have been identified in cattle i.e. PPARc1 and PPARc2 [16]. In vitro, the isoform PPARc2 has stronger adipogenic activity than the isoform PPARc1 [65] and might at least partly explain the positive associations we found between PPARc and PPARc2 vs. back fat thickness and BCS.
In the pluriparous cows of Trial 1, PPARc was reduced in s.c. AT and in the liver from late pregnancy to early lactation. We suggest that the expression patterns of PPARc and PPARc2 mRNA, with highest values prepartum compared to postpartum in AT, are possibly associated to the genes we analysed throughout lactation. Lower values in PPARc mRNA abundance in liver were also reported 3 wk after parturition as compared to the first wk of lactation [66]. The nadir observed for PPARc expression in AT but also in liver at wk 3 postpartum is in line with this report. In contrast, PPARa is upregulated after parturition [67]. As discussed by Lee et al. [68], the uppermost expression for PPARc is found in AT, where it is related to adipogenesis and IS. In liver, it is linked to fatty acid uptake, lipid storage, and reduced gluconeogenesis. We observed negative correlations between PPARc and NEFA, BHBA but, albeit weaker, also to insulin concentrations. This constellation in particular with regard to NEFA and BHBA may mitigate hepatic fatty acid uptake at least in pluriparous dairy cows. High abundance of PPARc mRNA in s.c. AT during late pregnancy might be related to energy accretion by AT. In liver, the decreased abundance of PPARc during early lactation may also result in increased gluconeogenic capacity.
In Trial 2, no changes in PPARc and PPARc2 were observed in s.c. AT after parturition in primiparous cows, which is in line with the results obtained in s.c. AT from Trial 1 with pluriparous cows. In addition to the potential effects of PPARc on gluconeogenesis, the incidence for fatty liver might be reduced by decreasing expression of PPARc2 after parturition: inactivation of PPARc2 may improve fatty liver induced by high fat diet in mice [69]. This suggestion is in line with our observation about the negative correlations between PPARc vs. NEFA and BHBA in pluriparous cows as discussed above.
For the mammary gland, Bionaz and Loor [70] reported that PPARc is upregulated at the onset of lactation in comparison to 15 DIM. Comparing 1, 42 and 105 DIM in our study, PPARc2 mRNA remained constant in mammary gland, similar to another study [71].
Gene expression of the pro-inflammatory proteins IL-6 and TNF-a
The adipokines IL-6 and TNF-a are related to insulin resistance [72]. In general, TNF-a reduces IS by phosphorylation of insulin receptor substrate which in turn blocks insulin signalling [73]. In consequence, plasma NEFA concentrations rise thus further enhancing insulin receptor substrate phosphorylation indirectly [74]. TNF-a increases factors related to fatty acid uptake in liver but does not affect genes associated with fatty acid oxidation as demonstrated in TNF-a knockout mice [75]. Accordingly, the intracellular rise of NEFA in non-alcoholic fatty liver increases the hepatic diacylglycerol content, which further reduces IS as discussed above [76].
In dairy cows, plasma TNF-a peaks transiently at parturition and decreases thereafter below prepartum values [77]; high TNFa serum concentrations are associated with insulin resistance and fatty liver 1 to 2 wk after calving [78]. For serum IL-6, no differences were observed at different stages of lactation [79]. Limited data is available about the mRNA expression of both adipokines in bovine AT, but liver TNF-a mRNA abundance is higher than the amount of IL-6 mRNA in liver [80].
In Trial 1, the TNF-a mRNA abundance remained fairly constant throughout lactation in liver and in s.c. AT. Besides confirming earlier reports focusing on the transition period [66,81], we also provided data about TNF-a mRNA expression in later stages of lactation. In contrast to TNFa, IL-6 mRNA in liver decreased during the transition period. IL-6 not only has adverse effects on liver metabolism but also stimulates hepatocyte proliferation and prevents cell damage [82]. The observed decline towards early lactation might be related to a higher susceptibility for metabolic health disorders during early lactation [83]. The abundance of IL-6 mRNA in s.c. AT of pluriparous cows was too low for valid quantification by our protocol and was not detectable in each sample.
In primiparous cows (Trial 2) we observed an increase in TNF-a mRNA abundance in all s.c. AT, at least at one of the sampling dates after parturition, which is different from the studies of Sadri et al. [81] and van Dorland et al. [66] reporting no changes during early lactation in pluriparous cows. In retroperitoneal AT TNF-a mRNA was less abundant at 105 than on 1 DIM. In contrast to pluriparous cows we observed an increase of TNF-a mRNA abundance in liver. Therefore, we suggest that TNF-a mRNA expression differs between primiparous and pluriparous cows and between different AT of primiparous cows. Differences between the different AT depots might be related to local regulation of IS in either a paracrine or autocrine manner.
The IL-6 mRNA in s.c. AT from tail head and withers was reduced from 1 to 105 DIM. Therefore, the regulation of the proinflammatory proteins IL-6 and TNF-a seems different at least at the level of the mRNA in these AT depots. In contrast, the IL-6 mRNA in liver increased during this time. This observation might be related to an improved protective and regenerative capacity related to IL-6 at this stage of lactation [82] in primiparous, and, to a lower extent, in pluriparous cows, where only a numerical increase was observed. On the other hand, IL-6 is associated with reduced IS [72]. Network analysis in cows revealed that IL-6 is related to many liver specific pathways e.g. lipoprotein metabolism and fatty acid oxidation and increases during ketosis [84]. Therefore, lower abundance in liver around parturition and early lactation might be advantageous in relation to IS but the protection of the liver against oxidative stress during this time period might be reduced as mentioned previously.
As discussed by Alluwaimi [85], TNF-a is suggested to maintain and regulate immunological functions within the mammary gland and the decrease from 1 to 42 or 105 DIM, respectively, may reflect the adaptation of the gland to galactopoesis as compared to colostrogenesis and lactogenesis.
Effects of CLA supplementation on gene expression
In general the effects of CLA supplementation were scarce and not consistent between genes, tissues and time points throughout lactation. However, in case of differences between control and CLA cows, the abundance of the respective target mRNA was always decreased, with only one exception (LEPRB in liver of pluriparous cows). CLA reduced ADIPOQ mRNA in omental and ADIPOR2 mRNA in retroperitoneal AT at 105 DIM. In addition, lower TNF-a mRNA abundance was observed in mesenteric AT of the CLA group, although only at 42 DIM. Comparable to a study of Sigl et al. [86] using primiparous Brown Swiss cows, we detected no CLA effect on hepatic PPARc mRNA abundance in primiparous German Holstein cows, whereas CLA treatment in pluriparous cows resulted in sporadically decreased hepatic PPARc mRNA abundance. CLA are ligands for PPARs [19,20]. Nevertheless, the effects of CLA in primiparous dairy cattle may differ from the effects of the PPARc agonist 2,4-thiazolidinedione observed in bulls [87]. However, the effects of 2,4-thiazolidinedione were not consistent in bulls comparable to our CLA effects. Therefore, inconsistency for PPARc agonists like 2,4-thiazolidinedione and CLA between different tissues or AT depots in cattle may exists. This result is in contrast to the situation in mice, were broad consistency was observed regarding to the effects of both compounds on glucose and lipid metabolism [88].
We suppose that the impact of CLA supplementation is more important for the mammary gland than for all other analysed tissues. Treatment with CLA decreased PPARc2 mRNA at both time points investigated in primiparous cows, accounting for 2 out of 5 CLA effects we observed in total. PPARc agonists have been demonstrated to stimulate lipid synthesis in a mammary epithelial cell line [14]. We observed lower PPARc2 mRNA abundance at 42 and at 105 DIM in the animals treated with CLA than in the cows receiving the control fat supplement. The mammary gland is an important target for CLA; as discussed by Bauman et al. [24], milk fat depression by CLA treatment is well established and was observed in the animals from the two trials herein, too [26,29]. PPARc was suggested by Kadegowda et al. [14] to be an important regulator of milk fat synthesis. It was shown that in parallel to PPARc, lipogenic genes like stearoyl-CoA desaturase, diacylglycerol O-acyltransferase 2 and fatty acid synthase were down regulated during the peripartal period [89]. We speculate that reduced receptor capacity of the PPARc2 isoform might be responsible for milk fat depression by CLA.
CLA supplementation and systemic insulin sensitivity
For estimating the systemic IS, hyperinsulinaemic euglycaemic glucose clamps (HEC) are considered as ‘‘gold standard’’, however, HEC can not realistically be accomplished at all days sampling was targeted. We used a surrogate marker, i.e. RQUICKI that can be accessed from individual more frequently obtained blood samples. The validity of the RQUICKI index was confirmed in dairy cows using data from previous clamp studies [90]. No CLA effects were observed on NEFA and BHBA concentrations [28], nevertheless, the RQUICKI index was reduced by 8% during the supplementation period and by 17% during the depletion period after 182 DIM in CLA treated cows vs. control cows indicating a reduction of IS by CLA. This was mainly attributable to greater insulin concentrations in CLA treated animals. Comparable effects on insulin concentrations using a similar CLA mixture were also observed in a mouse model, in which the results were explained by an increase in ß-cell mass and number [91]. In contrast to our results, Ho¨tger et al. [25] observed no increase of insulin during supplementation with CLA for 9 weeks in dairy cows during their second lactation. The role of CLA as a ligand for the free fatty acid receptor 1 increasing glucose stimulated insulin secretion [22] might at least partly account for the insulin-increasing effects of CLA we observed in pluriparous cows. In addition, signalling by this receptor could also increase incretins, leading to insulin secretion [92]. Differences in the expression of free fatty acid receptor 1 on pancreatic ß-cells and on enteroendocrine cells between primiparous vs. pluriparous cows could not be excluded. Parity seems important in this context since we did not see the CLA effect on insulin or RQUICKI when studying an additional set of 10 primiparous cows in Trial 1 (5 cows each supplemented with either CLA or control fat) from which only blood samples were assessed (data not shown).
Conclusion
Studying primiparous and pluriparous cows, we observed changes in the mRNA abundance of different modulators of IS like the ADIPOQ and LEP system in different tissues focusing on AT and liver. Besides TNF-a, the abundance of all mRNAs of interest decreased after parturition in s.c. AT of pluriparous cows but not in liver where most of the mRNAs increased except IL-6. By our study we have shown in primiparous cows that individual AT depots of cows may differentially influence the regulation of IS. Based on the number of changes in different tissues, the impact of different genes on metabolism during lactation may also vary. Most changes were observed in case of ADIPOQR1 and TNF-a. Most of the differences occurred in retroperitoneal AT. Differences between different AT depots, even of s.c. origin should be considered for further studies. Finally, we showed that PPARc2, as a regulator of milk fat synthesis in cattle, is reduced in mammary gland by CLA supplementation which might be a direct link to CLA induced milk fat depression.
Acknowledgments
The authors would like to thank the technical assistance of Inga Hofs, Isabella Israel, Birgit Mielenz and Karin Strack (Institute of Animal Science, Physiology & Hygiene Unit, University of Bonn). We also appreciate Iris Gockel-Bo¨hner (Institute of Animal Science, Physiology & Hygiene Unit, University of Bonn) for her help during sampling of the cows.
Author Contributions
Conceived and designed the experiments: SD GB JR HS MM. Performed the experiments: BS SW PF AK JR MM SH. Analyzed the data: BS PF MM. Wrote the paper: BS MM HS.
This article was originally published in PLoS ONE 9(1): e86211. doi:10.1371/journal.pone.0086211. This is an Open Access article distributed under the terms of the Creative Commons Attribution License.
References
1. Bell AW, Bauman DE (1997) Adaptations of glucose metabolism during pregnancy and lactation. J Mammary Gland Biol Neoplasia
2: 265–278. 2. Lain KY, Catalano PM (2007) Metabolic changes in pregnancy. Clin Obstet Gynecol 50: 938–948.
3. LeBlanc S (2010) Monitoring metabolic health of dairy cattle in the transition period. J Reprod Dev 56 Suppl: S29–35.
4. Antuna-Puente B, Feve B, Fellahi S, Bastard JP (2008) Adipokines: The missing link between insulin resistance and obesity. Diabetes & Metabolism 34: 2–11.
5. Karastergiou K, Mohamed-Ali V (2010) The autocrine and paracrine roles of adipokines. Mol Cell Endocrinol 318: 69–78.
6. Ahima RS, Osei SY (2008) Adipokines in obesity. Front Horm Res 36: 182–197.
7. Kadowaki T, Yamauchi T, Kubota N, Hara K, Ueki K, et al. (2006) Adiponectin and adiponectin receptors in insulin resistance, diabetes, and the metabolic syndrome. The Journal of Clinical Investigation 116: 1784–1792.
8. Bartha T, Sayed-Ahmed A, Rudas P (2005) Expression of leptin and its receptors in various tissues of ruminants. Domest Anim Endocrinol 29: 193–202.
9. Kadowaki T, Yamauchi T (2005) Adiponectin and adiponectin receptors. Endocr Rev 26: 439–451.
10. Hotamisligil GS (2000) Molecular mechanisms of insulin resistance and the role of the adipocyte. Int J Obes Relat Metab Disord 24 Suppl 4: S23–27.
11. Nieto-Vazquez I, Fernandez-Veledo S, de Alvaro C, Rondinone CM, Valverde AM, et al. (2007) Protein-tyrosine phosphatase 1B-deficient myocytes show increased insulin sensitivity and protection against tumor necrosis factor-alphainduced insulin resistance. Diabetes 56: 404–413.
12. van Doorn M, Kemme M, Ouwens M, van Hoogdalem EJ, Jones R, et al. (2006) Evaluation of proinflammatory cytokines and inflammation markers as biomarkers for the action of thiazolidinediones in Type 2 diabetes mellitus patients and healthy volunteers. Br J Clin Pharmacol 62: 391–402.
13. Tontonoz P, Spiegelman BM (2008) Fat and beyond: the diverse biology of PPARgamma. Annu Rev Biochem 77: 289–312.
14. Kadegowda AKG, Bionaz M, Piperova LS, Erdman RA, Loor JJ (2009) Peroxisome proliferator-activated receptor-c activation and long-chain fatty acids alter lipogenic gene networks in bovine mammary epithelial cells to various extents. J Dairy Sci 92: 4276–4289.
15. Schoenberg KM, Overton TR (2011) Effects of plane of nutrition and 2,4- thiazolidinedione on insulin responses and adipose tissue gene expression in dairy cattle during late gestation. J Dairy Sci 94: 6021–6035.
16. Sundvold H, Brzozowska A, Lien S (1997) Characterisation of bovine peroxisome proliferator-activated receptors gamma 1 and gamma 2: genetic mapping and differential expression of the two isoforms. Biochem Biophys Res Commun 239: 857–861.
17. Hammarstedt A, Andersson CX, Rotter Sopasakis V, Smith U (2005) The effect of PPARc ligands on the adipose tissue in insulin resistance. Prostaglandins Leukot Essent Fatty Acids 73: 65–75.
18. Wronska A, Kmiec Z (2012) Structural and biochemical characteristics of various white adipose tissue depots. Acta Physiol (Oxf) 205: 194–208.
19. Bensinger SJ, Tontonoz P (2008) Integration of metabolism and inflammation by lipid-activated nuclear receptors. Nature 454: 470–477.
20. Reynolds CM, Roche HM (2010) Conjugated linoleic acid and inflammatory cell signalling. Prostaglandins Leukot Essent Fatty Acids 82: 199–204.
21. Miller JR, Siripurkpong P, Hawes J, Majdalawieh A, Ro H-S, et al. (2008) The trans-10, cis-12 isomer of conjugated linoleic acid decreases adiponectin assembly by PPARc-dependent and PPARc-independent mechanisms. J Lipid Res 49: 550–562.
22. Schmidt J, Liebscher K, Merten N, Grundmann M, Mielenz M, et al. (2011) Conjugated linoleic acids mediate insulin release through islet G protein-coupled receptor FFA1/GPR40. J Biol Chem 286: 11890–11894.
23. Vyas D, Kadegowda AK, Erdman RA (2012) Dietary conjugated linoleic Acid and hepatic steatosis: species-specific effects on liver and adipose lipid metabolism and gene expression. J Nutr Metab 2012: 932928.
24. Bauman DE, Harvatine KJ, Lock AL (2011) Nutrigenomics, rumen-derived bioactive fatty acids, and the regulation of milk fat synthesis. Annu Rev Nutr 31: 299–319.
25. Hotger K, Hammon HM, Weber C, Gors S, Troscher A, et al. (2013) Supplementation of conjugated linoleic acid in dairy cows reduces endogenous glucose production during early lactation. J Dairy Sci 96: 2258–2270.
26. Holtenius P, Holtenius K (2007) A model to estimate insulin sensitivity in dairy cows. Acta Vet Scand 49: 29.
27. GfE (2001) (German Society of Nutrition Physiology). Ausschuss fu¨ r Bedarfsnormen der Gesellschaft fu¨ r Erna¨hrungsphysiologie. Nr. 8. Empfehlungen zur Energie- und Na¨hrstoffversorgung der Milchku¨he und Aufzuchtrinder (Recommendations of energy and nutrient supply for dairy cows and breeding cattle). DLG-Verlag, Frankfurt am Main, Germany.
28. Pappritz J, Meyer U, Kramer R, Weber E-M, Jahreis G, et al. (2011) Effects of long-term supplementation of dairy cow diets with rumen-protected conjugated linoleic acids (CLA) on performance, metabolic parameters and fatty acid profile in milk fat. Arch Anim Nutr 65: 89–107.
29. Saremi B, Al-Dawood A, Winand S, Muller U, Pappritz J, et al. (2012) Bovine haptoglobin as an adipokine: serum concentrations and tissue expression in dairy cows receiving a conjugated linoleic acids supplement throughout lactation. Vet Immunol Immunopathol 146: 201–211.
30. von Soosten D, Meyer U, Piechotta M, Flachowsky G, Danicke S (2012) Effect of conjugated linoleic acid supplementation on body composition, body fat mobilization, protein accretion, and energy utilization in early lactation dairy cows. J Dairy Sci 95: 1222–1239.
31. von Soosten D, Meyer U, Weber EM, Rehage J, Flachowsky G, et al. (2011) Effect of trans-10, cis-12 conjugated linoleic acid on performance, adipose depot weights, and liver weight in early-lactation dairy cows. J Dairy Sci 94: 2859– 2870.
32. Saremi B, Sauerwein H, Da¨nicke S, Mielenz M (2012) Technical note: Identification of reference genes for gene expression studies in different bovine tissues focusing on different fat depots. J Dairy Sci 95: 3131–3138.
33. Hellemans J, Mortier G, De Paepe A, Speleman F, Vandesompele J (2007) qBase relative quantification framework and software for management and automated analysis of real-time quantitative PCR data. Genome Biology 8: R19.
34. Sauerwein H, Heintges U, Hennies M, Selhorst T, Daxenberger A (2004) Growth hormone induced alterations of leptin serum concentrations in dairy cows as measured by a novel enzyme immunoassay. Livestock Prod Sci 87: 189– 195.
35. Perseghin G, Caumo A, Caloni M, Testolin G, Luzi L (2001) Incorporation of the fasting plasma FFA concentration into QUICKI improves its association with insulin sensitivity in nonobese individuals. J Clin Endocrinol Metab 86: 4776–4781.
36. Minokoshi Y, Kim YB, Peroni OD, Fryer LG, Muller C, et al. (2002) Leptin stimulates fatty-acid oxidation by activating AMP-activated protein kinase. Nature 415: 339–343.
37. Loffreda S, Yang SQ, Lin HZ, Karp CL, Brengman ML, et al. (1998) Leptin regulates proinflammatory immune responses. FASEB J 12: 57–65.
38. Gauthier MS, O’Brien EL, Bigornia S, Mott M, Cacicedo JM, et al. (2011) Decreased AMP-activated protein kinase activity is associated with increased inflammation in visceral adipose tissue and with whole-body insulin resistance in morbidly obese humans. Biochem Biophys Res Commun 404: 382–387.
39. Yamauchi T, Kamon J, Minokoshi Y, Ito Y, Waki H, et al. (2002) Adiponectin stimulates glucose utilization and fatty-acid oxidation by activating AMPactivated protein kinase. Nat Med 8: 1288–1295.
40. Leury BJ, Baumgard LH, Block SS, Segoale N, Ehrhardt RA, et al. (2003) Effect of insulin and growth hormone on plasma leptin in periparturient dairy cows. Am J Physiol Regul Integr Comp Physiol 285: R1107–1115.
41. Toyoshima Y, Gavrilova O, Yakar S, Jou W, Pack S, et al. (2005) Leptin improves insulin resistance and hyperglycemia in a mouse model of type 2 diabetes. Endocrinology 146: 4024–4035.
42. Giesy SL, Yoon B, Currie WB, Kim JW, Boisclair YR (2012) Adiponectin deficit during the precarious glucose economy of early lactation in dairy cows. Endocrinology 153: 5834–5844.
43. Mielenz M, Mielenz B, Singh SP, Kopp C, Heinz J, et al. (2013) Development, validation, and pilot application of a semiquantitative Western blot analysis and an ELISA for bovine adiponectin. Domest Anim Endocrinol 44: 121–130.
44. Koltes DA, Spurlock DM (2012) Adipose tissue angiopoietin-like protein 4 messenger RNA changes with altered energy balance in lactating Holstein cows. Domest Anim Endocrinol 43: 307–316.
45. Lemor A, Hosseini A, Sauerwein H, Mielenz M (2009) Transition period-related changes in the abundance of the mRNAs of adiponectin and its receptors, of visfatin, and of fatty acid binding receptors in adipose tissue of high-yielding dairy cows. Domest Anim Endocrinol 37: 37–44.
46. Thorn SR, Ehrhardt RA, Butler WR, Quirk SM, Boisclair YR (2008) Insulin regulates hepatic leptin receptor expression in early lactating dairy cows. Am J Physiol Regul Integr Comp Physiol 295: R1455–R1462.
47. Faraj M, Beauregard G, Tardif A, Loizon E, Godbout A, et al. (2008) Regulation of leptin, adiponectin and acylation-stimulating protein by hyperinsulinaemia and hyperglycaemia in vivo in healthy lean young men. Diabetes Metab 34: 334–342.
48. Mielenz M, Kuhla B, Hammon HM (2013) Abundance of adiponectin system and G-protein coupled receptor GPR109A mRNA in adipose tissue and liver of F2 offspring cows of Charolais x German Holstein crosses that differ in body fat accumulation. J Dairy Sci 96: 278–289.
49. Yamauchi T, Nio Y, Maki T, Kobayashi M, Takazawa T, et al. (2007) Targeted disruption of AdipoR1 and AdipoR2 causes abrogation of adiponectin binding and metabolic actions. Nat Med 13: 332–339.
50. Fu Y, Luo N, Klein RL, Garvey WT (2005) Adiponectin promotes adipocyte differentiation, insulin sensitivity, and lipid accumulation. J Lipid Res 46: 1369– 1379.
51. Loor JJ, Dann HM, Guretzky NA, Everts RE, Oliveira R, et al. (2006) Plane of nutrition prepartum alters hepatic gene expression and function in dairy cows as assessed by longitudinal transcript and metabolic profiling. Physiol Genomics 27: 29–41.
52. Schams DF, Graf B, Graule M, Abele M, Prokopp S (1991) Hormonal changes during lactation in cows of three different breeds. Livest Prod Sci 27: 285–296.
53. Qin Y, Tian YP (2010) Hepatic adiponectin receptor R2 expression is upregulated in normal adult male mice by chronic exogenous growth hormone levels. Mol Med Rep 3: 525–530.
54. Soardo G, Donnini D, Domenis L, Catena C, De Silvestri D, et al. (2011) Oxidative stress is activated by free fatty acids in cultured human hepatocytes. Metab Syndr Relat Disord 9: 397–401.
55. Bernabucci U, Ronchi B, Lacetera N, Nardone A (2005) Influence of body condition score on relationships between metabolic status and oxidative stress in periparturient dairy cows. J Dairy Sci 88: 2017–2026.
56. Khan MJ, Hosseini A, Burrell S, Rocco SM, McNamara JP, et al. (2013) Change in subcutaneous adipose tissue metabolism and gene network expression during the transition period in dairy cows, including differences due to sire genetic merit. J Dairy Sci 96: 2171–2182.
57. Dagon Y, Avraham Y, Berry EM (2006) AMPK activation regulates apoptosis, adipogenesis, and lipolysis by eIF2alpha in adipocytes. Biochem Biophys Res Commun 340: 43–47.
58. Rossmeisl M, Flachs P, Brauner P, Sponarova J, Matejkova O, et al. (2004) Role of energy charge and AMP-activated protein kinase in adipocytes in the control of body fat stores. Int J Obes Relat Metab Disord 28 Suppl 4: S38–44.
59. Daval M, Diot-Dupuy F, Bazin R, Hainault I, Viollet B, et al. (2005) Antilipolytic action of AMP-activated protein kinase in rodent adipocytes. J Biol Chem 280: 25250–25257.
60. Gauthier MS, Miyoshi H, Souza SC, Cacicedo JM, Saha AK, et al. (2008) AMP-activated protein kinase is activated as a consequence of lipolysis in the adipocyte: potential mechanism and physiological relevance. J Biol Chem 283: 16514–16524.
61. Locher LF, Rehage J, Khraim N, Meyer U, Danicke S, et al. (2012) Lipolysis in early lactation is associated with an increase in phosphorylation of adenosine monophosphate-activated protein kinase (AMPK)a1 in adipose tissue of dairy cows. J Dairy Sci 95: 2497–2504.
62. Yang YK, Chen M, Clements RH, Abrams GA, Aprahamian CJ, et al. (2008) Human mesenteric adipose tissue plays unique role versus subcutaneous and omental fat in obesity related diabetes. Cell Physiol Biochem 22: 531–538.
63. Rosen ED, MacDougald OA (2006) Adipocyte differentiation from the inside out. Nat Rev Mol Cell Biol 7: 885–896.
64. Matsusue K, Haluzik M, Lambert G, Yim S-H, Gavrilova O, et al. (2003) Liverspecific disruption of PPARc in leptin-deficient mice improves fatty liver but aggravates diabetic phenotypes. J Clin Invest 111: 737–747.
65. Mueller E, Drori S, Aiyer A, Yie J, Sarraf P, et al. (2002) Genetic analysis of adipogenesis through peroxisome proliferator-activated receptor gamma isoforms. J Biol Chem 277: 41925–41930.
66. van Dorland HA, Sadri H, Morel I, Bruckmaier RM (2011) Coordinated gene expression in adipose tissue and liver differs between cows with high or low NEFA concentrations in early lactation. J Anim Physiol Anim Nutr (Berl) 96: 137–147.
67. Loor JJ, Dann HM, Everts RE, Oliveira R, Green CA, et al. (2005) Temporal gene expression profiling of liver from periparturient dairy cows reveals complex adaptive mechanisms in hepatic function. Physiol Genomics 23: 217–226.
68. Lee CH, Olson P, Evans RM (2003) Minireview: lipid metabolism, metabolic diseases, and peroxisome proliferator-activated receptors. Endocrinology 144: 2201–2207.
69. Yamazaki T, Shiraishi S, Kishimoto K, Miura S, Ezaki O (2011) An increase in liver PPARc2 is an initial event to induce fatty liver in response to a diet high in butter: PPARc2 knockdown improves fatty liver induced by high-saturated fat. J Nutr Biochem 22: 543–553.
70. Bionaz M, Loor J (2008) Gene networks driving bovine milk fat synthesis during the lactation cycle. BMC Genomics 9: 366.
71. Komatsu M, Ohfusa H, Aizawa T, Hashizume K (2007) Adiponectin inversely correlates with high sensitive C-reactive protein and triglycerides, but not with insulin sensitivity, in apparently healthy Japanese men. Endocr J 54: 553–558.
72. Bruun JM, Verdich C, Toubro S, Astrup A, Richelsen B (2003) Association between measures of insulin sensitivity and circulating levels of interleukin-8, interleukin-6 and tumor necrosis factor-alpha. Effect of weight loss in obese men. Eur J Endocrinol 148: 535–542.
73. Boura-Halfon S, Zick Y (2009) Phosphorylation of IRS proteins, insulin action, and insulin resistance. Am J Physiol Endocrinol Metab 296: E581–591.
74. Boden G (2005) Free fatty acids and insulin secretion in humans. Curr Diab Rep 5: 167–170.
75. Salles J, Tardif N, Landrier JF, Mothe-Satney I, Guillet C, et al. (2012) TNFa gene knockout differentially affects lipid deposition in liver and skeletal muscle of high-fat-diet mice. J Nutr Biochem 23: 1685–1693.
76. Jornayvaz FR, Shulman GI (2012) Diacylglycerol activation of protein kinase Ce and hepatic insulin resistance. Cell Metab 15: 574–584.
77. Schoenberg KM, Perfield KL, Farney JK, Bradford BJ, Boisclair YR, et al. (2011) Effects of prepartum 2,4-thiazolidinedione on insulin sensitivity, plasma concentrations of tumor necrosis factor-a and leptin, and adipose tissue gene expression. J Dairy Sci 94: 5523–5532.
78. Ohtsuka H, Koiwa M, Hatsugaya A, Kudo K, Hoshi F, et al. (2001) Relationship between serum TNF activity and insulin resistance in dairy cows affected with naturally occurring fatty liver. J Vet Med Sci 63: 1021–1025.
79. Hagiwara K, Yamanaka H, Hisaeda K, Taharaguchi S, Kirisawa R, et al. (2001) Concentrations of IL-6 in serum and whey from healthy and mastitic cows. Vet Res Commun 25: 99–108.
80. Vels L, Røntved CM, Bjerring M, Ingvartsen KL (2009) Cytokine and acute phase protein gene expression in repeated liver biopsies of dairy cows with a lipopolysaccharide-induced mastitis. J Dairy Sci 92: 922–934.
81. Sadri H, Bruckmaier RM, Rahmani HR, Ghorbani GR, Morel I, et al. (2010) Gene expression of tumour necrosis factor and insulin signalling-related factors in subcutaneous adipose tissue during the dry period and in early lactation in dairy cows. J Anim Physiol Anim Nutr (Berl) 94: e194–e202.
82. Cressman DE, Greenbaum LE, DeAngelis RA, Ciliberto G, Furth EE, et al. (1996) Liver failure and defective hepatocyte regeneration in interleukin-6- deficient mice. Science 274: 1379–1383.
83. Drackley JK (1999) ADSA Foundation Scholar Award. Biology of dairy cows during the transition period: the final frontier? J Dairy Sci 82: 2259–2273.
84. Loor JJ, Everts RE, Bionaz M, Dann HM, Morin DE, et al. (2007) Nutritioninduced ketosis alters metabolic and signaling gene networks in liver of periparturient dairy cows. Physiol Genomics 32: 105–116.
85. Alluwaimi AM (2004) The cytokines of bovine mammary gland: prospects for diagnosis and therapy. Res Vet Sci 77: 211–222.
86. Sigl T, Schlamberger G, Kienberger H, Wiedemann S, Meyer H, et al. (2010) Rumen-protected conjugated linoleic acid supplementation to dairy cows in late pregnancy and early lactation: effects on milk composition, milk yield, blood metabolites and gene expression in liver. Acta Vet Scand 52: 16.
87. Arevalo-Turrubiarte M, Gonzalez-Davalos L, Yabuta A, Garza JD, Davalos JL, et al. (2012) Effect of 2,4-thiazolidinedione on limousin cattle growth and on muscle and adipose tissue metabolism. PPAR Res 2012: 891841.
88. Wendel AA, Belury MA (2006) Effects of conjugated linoleic acid and troglitazone on lipid accumulation and composition in lean and Zucker diabetic fatty (fa/fa) rats. Lipids 41: 241–247.
89. Ji P, Osorio JS, Drackley JK, Loor JJ (2012) Overfeeding a moderate energy diet prepartum does not impair bovine subcutaneous adipose tissue insulin signal transduction and induces marked changes in peripartal gene network expression. J Dairy Sci 95: 4333–4351.
90. Kusenda M, Kaske M, Piechotta M, Locher L, Starke A, et al. (2013) Effects of dexamethasone-21-isonicotinate on peripheral insulin action in dairy cows 5 days after surgical correction of abomasal displacement. J Vet Intern Med 27: 200–206.
91. Poirier H, Rouault C, Cle´ment L, Niot I, Monnot MC, et al. (2005) Hyperinsulinaemia triggered by dietary conjugated linoleic acid is associated with a decrease in leptin and adiponectin plasma levels and pancreatic beta cell hyperplasia in the mouse. Diabetologia 48: 1059–1065.
92. Luo J, Swaminath G, Brown SP, Zhang J, Guo Q, et al. (2012) A potent class of GPR40 full agonists engages the enteroinsular axis to promote glucose control in rodents. PLoS One 7: e46300.
93. Yuen BS, McMillen IC, Symonds ME, Owens PC (1999) Abundance of leptin mRNA in fetal adipose tissue is related to fetal body weight. J Endocrinol 163: R11–14.
94. Hosseini A, Behrendt C, Regenhard P, Sauerwein H, Mielenz M (2011) Differential effects of propionate or beta-hydroxybutyrate on genes related to energy balance and insulin sensitivity in bovine white adipose tissue explants from a subcutaneous and a visceral depot. J Anim Physiol Anim Nutr (Berl) 96: 570–580.