The domestic chicken provides a widespread and relatively inexpensive source of dietary protein for humans. In addition to its role as a food animal, the chicken has a long history as a valuable model research organism 1 . These dual considerations led to the selection of chicken as the first agricultural animal model to be sequenced at the genome level 2 . While chickens have been used heavily for studies of developmental biology and immunology, a number of traits make them a viable model for studies of adipose biology, obesity and insulin resistance. Commercial broiler chickens, in particular, rapidly accumulate excess adipose tissue as a result of genetic selection for growth and are considered “obese” relative to leaner egg-laying or wild strains of chickens(rev. in 3 ). Chickens mimic the early stage of type 2 diabetes in humans, exhibiting both hyperglycemia (up to 200 mg/dL in the fasting state) and resistance to exogenous insulin 4 5 . Like humans, but unlike rodents or pigs, chickens rely on liver rather than adipose tissue for the majority of de novo lipid synthesis 6 7 8 . Most metabolic genes are conserved with humans, and a number of the quantitative trait loci (QTLs) that have been linked to fatness in chickens contain genes implicated in human susceptibility to obesity or diabetes 9 . Chickens also represent a model for studying mechanisms of adipocyte hyperplasia during development, a process that may exacerbate adult obesity. During at least the first several weeks after hatch, chicken adipose tissue expands more through adipocyte hyperplasia than hypertrophy, and an early increase in adipocyte number is a common feature of some lines genetically selected for excess adiposity 10 11 . Finally, the egg presents opportunities to directly manipulate the developmental milieu and study the consequences on adipose metabolism via in ovo injection.

Relatively little is known about regulation of adipose tissue deposition and metabolism in chicken. Because of its relative importance in lipogenesis, most studies have focused on the role of liver in adipose expansion. Several genetic lines of fat and lean chickens have been developed through phenotypic selection, most of which have both elevated plasma levels of very low density lipoprotein (VLDL) and lower levels of plasma glucose, reflecting the importance of hepatic lipogenesis and glucose consumption in fat accretion. Reciprocally, phenotypic selection for low plasma glucose simultaneously selects for fatness 12 . Both chicken and mammalian adipocytes develop through a sequence of molecular triggers including activation of CCAAT-enhancer-binding protein alpha (CEBPα) and peroxisome proliferator-activated receptor gamma (PPARγ) 13 . A clear point of divergence, however, is their responsiveness to insulin. Unlike in mammals, insulin has minimal effect on glucose uptake in chicken adipose tissue 14 . In fact, an avian homolog of the insulin-sensitive glucose transporter GLUT4 has not been identified in the current chicken genome database. Insulin does, however, stimulate uptake of acetate, which is the preferred substrate for de novo lipogenesis in chicken adipocytes, although the magnitude of the effect is relatively modest 15 . Insulin signaling appears to proceed through tissue specific cascades in chicken metabolic tissues. In liver, insulin elicits a signaling cascade that parallels the response in mammals, including tyrosine phosphorylation of insulin receptor β-subunit (IRβ), insulin receptor substrate-1 (IRS-1) and Src homology 2 domain-containing substrate (Shc) and activation of phosphatidylinositol 3-kinase (PI3K) 16 17 . The situation in skeletal muscle is more complex. Tyrosine phosphorylation of IRβ and IRS-1 and PI3K activity are not regulated by insulin, whereas events downstream of PI3K (e.g. Akt and P70S6K activation) are accordingly sensitive 18 . We recently reported that insulin also does not elicit a classical IRβ initiated cascade in chicken adipose tissue, including the downstream steps of Akt and P70S6K activation 19 . Insulin also does not inhibit lipolysis in chicken adipose tissue; glucagon, is the primary lipolytic hormone (rev. in 20 ).

In the present study we simultaneously characterized the effects of a short term (5 hours) fast or neutralization of insulin action (5 hours) on adipose tissue of young (16–17 day-old), fed commercial broiler chickens. The goals of this study were two-fold. First, we sought to identify pathways activated by feed restriction, reasoning that they may highlight potential strategies for control of fatness through either genetic selection or improved management practices. Simultaneously, we sought to understand the contribution of insulin, if any, into chicken adipose physiology. No experimental model of diabetes exist in chicken: total pancreatectomies are not achievable, and alloxan and streptozotocin are inefficient at destroying pancreatic chicken beta-cells (rev. in 5 ). The two treatments were compared to distinguish potential insulin-specific changes from those that could be mimicked by fasting through changes in nutrient availability. Both treatments were shown previously to elicit significant alterations in several plasma metabolic and endocrine parameters 18 ; in the studies reported herein, samples of abdominal adipose tissue were issued from the same experiment. Tissue metabolomics was combined with microarrays to bridge the gap between gene expression, metabolic and physiological responses, and to identify the composite effects of both fasting and insulin deprivation on chicken adipose tissue.

Expression levels of a total of 2016 genes were significantly altered by fasting and/or insulin neutralization when compared to fed controls based on an FDR adjusted p-value < 0.05 ( Additional file 1; Figure 1A). Sixty-nine percent of these genes showed a fold-change ≥ |1.5| (Figure 1B). The majority of changes in expression were attributable to fasting, with 917 up-regulated and 863 down-regulated genes in fasted vs. fed adipose tissue. Insulin neutralization altered expression of 92 genes, 72 of which were also differentially expressed with fasting (Figure 1A). All genes that were affected by both treatments changed in the same direction (i.e., up- or down-regulated in both groups). Real-time RT-PCR was employed to validate differential expression based on the microarray data. Eleven genes were selected based on fold-change or biological functions of interest (Table 1). Differential expression under fasting versus fed conditions was validated for all genes except pre-B-cell leukemia homeobox 3 (PBX3). Ten of the eleven genes were also differentially expressed in insulin neutralized compared to fed birds based on QPCR.

FDR p-value: * p < 0.05, ** p < 0.005, *** p < 0.0005.

Genes that were differentially expressed in at least one pairwise comparison were clustered to visualize the similarities between groups and to determine if insulin-neutralized expression profiles were more similar to fasted or to fed status. As shown in Figure 2A, samples within each of the three experimental groups clustered together. The dendrogram also showed that the fasting group was distant from fed and insulin-neutralized groups, which were closer to each other. To further visualize relationships between treatments with regard to gene expression, distinct clusters of genes were extracted and submitted to gene set enrichment analysis to identify GO terms and pathways that were significantly overrepresented among genes contained in these clusters. Seven clusters represented four general patterns of similarities between treatments (Figure 2B). Clusters 1, 3 and 4 consisted of genes with higher expression in fasting compared to both insulin-neutralized and fed conditions, with insulin-neutralized intermediate between fasted and fed. This set of genes was significantly enriched in GO terms related to protein and lipid catabolism and to cell signaling, including regulation of the stress-sensitive NFκB cascade (Table 2). These three clusters were also enriched in members of the KEGG pathways ubiquitin mediated proteolysis, sphingolipid metabolism, PPAR signaling, fatty acid metabolism and the peroxisome. The rate-limiting genes for fatty acid oxidation (ACOX1 and CPT1A), along with fatty acid binding proteins 5 and 6, are contained in these three clusters. Clusters 5 and 7 also contained genes with higher levels in fasted vs. the other two groups, but with comparable expression levels between insulin-neutralized and fed, and thus no clear effect of insulin loss. These two clusters were significantly enriched in pathways related to signaling and metabolism, including enzyme linked receptor protein signaling pathway (p = 0.0097) and in the KEGG pathways for glycerolipid metabolism and PPAR. Genes responsible for the latter enrichment include PPARΔ, which was recently shown to increase total oxidative metabolism in white adipose tissue 21 . Clusters 2 and 6 contained genes expressed at lowest levels in fasted chickens. Genes in cluster 2 were expressed at intermediate levels in the insulin-neutralized group relative to fed and fasted. This set of genes was significantly enriched in GO annotations related to monosaccharide catabolic process and glucose metabolism, and in genes comprising the KEGG pathways for carbohydrate metabolism, TCA cycle and glycolysis (Table 2). Finally, cluster 6 consisted of genes that were also lowest in fasting but showed no clear effect of insulin loss, with similar expression in fed and insulin-neutralized groups. This set of genes was significantly enriched for the KEGG pathways steroid biosynthesis (p = 0.0043), glyoxylate and dicarboxylate metabolism (P = 0.012) and pyruvate metabolism (P = 0.033), along with a number of genes involved in lipid biosynthesis, which was the highest scoring GO category (p = 0.6). Cluster 8 was a distinct, small cluster with variable expression within group and no significant GO or KEGG annotations.

Global biological responses to fasting and to insulin neutralization were further characterized using KEGG pathway matching, based on genes with statistically significant differential expression (FDR < 0.05) and absolute fold-change ≥1.5. Genes altered exclusively by fasting represented a wide range of cellular pathways, indicating significant effects of even a five hour fast on adipose function and metabolism in chicken. Fasting exerted significant effects on pathways related to carbohydrate, amino acid and lipid metabolism and synthesis. Within the categories related to lipid metabolism, fasting up-regulated expression of genes involved in fatty acid oxidation (e.g., acyl-CoA oxidase 1 (ACOX1), acetyl-CoA carboxylase beta (ACACB), carnitine palmitoyltransferase 1A (CPT1A)) and down-regulated expression of genes that control fatty acid, cholesterol and triacylglycerol synthesis (e.g., 1-acylglycerol-3-phosphate O-acyltransferase 9 (AGPAT9), ATP citrate lyase (ACLY), farnesyl diphosphate synthase (FDPS), acetyl-Coenzyme A carboxylase alpha (ACACA) and acetoacetyl-CoA synthetase (AACS)). Fasting also up-regulated expression of many genes involved in proteolysis and amino acid degradation. In addition to pathways highlighted by KEGG analysis, fasting down-regulated a number of genes (e.g., TGFβ, BMP) that mediate mesenchymal stem cell (MSC) commitment, an early step in the formation of new adipocytes ( Additional file 1). Finally, a number of phosphodiesterases were up-regulated with fasting, presumably in response to the increased plasma glucagon 18 and subsequent elevations in cyclic adenosine monophosphate (cAMP; Additional file 1). Collectively, these categories indicate that chicken adipose tissue responds to a relatively short duration (five hour) fast with sweeping changes in gene expression that suppress synthesis and storage of lipids and other macromolecules and up-regulate mobilization and metabolism of fatty acids and proteins.

Loss of insulin action also resulted in significant effects on adipose gene expression, the majority of which overlapped with the response to fasting (Table 3; Additional file 2). Several genes central to energy metabolism were affected. Diacylglycerol O-acyltransferase homolog 2 (DGAT2), which catalyzes the final and only committed step in triacylglycerol synthesis, was down-regulated (10.5-and 6.1-fold, respectively, fasted and insulin neutralized) in both treatment groups relative to the fed group. Conversely, acyl-Coenzyme A binding domain containing 5 (ACBD5) and pyruvate dehydrogenase kinase 4 (PDK4) were significantly up-regulated in both treatments relative to fed controls. ACBD5 is one of a family of long chain fatty acyl CoA trafficking proteins that play roles in both triglyceride synthesis and beta-oxidation 22 . PDK4, which was up-regulated vs. fed by ~ 17-fold with fasting and 6-fold with insulin neutralization, acts as a fuel switch by phosphorylating and inactivating pyruvate dehydrogenase, shifting metabolism from glycolysis to fatty acid oxidation 23 . Fasting and insulin neutralization also up-regulated expression of the type I angiotensin II receptor (AGTR1). Angiotensin II alters adipocyte lipid metabolism and insulin signaling 24 25 26 , and increased AGTR1 expression in adipose tissue is associated with enhanced insulin sensitivity 27 . Finally, a number of genes regulated by both fasting and insulin neutralization function in general processes related to protein synthesis.

FDR p-value <0.05.

A total of thirteen genes (four up-regulated and nine down-regulated) were differentially expressed only with insulin neutralization (Table 4). The most interesting of these responses were upregulation of GCG, which encodes preproglucagon (fold change = 2.91), in parallel with downregulation of the glucagon receptor (LOC425670, fold change = −2.77). Other genes uniquely affected by insulin have less clear relevance to adipose biology according to current knowledge.

FDR p-value <0.05.

Tissue metabolomic analysis was used to identify the metabolic intermediates that were altered by fasting and insulin neutralization. A total of 92 metabolites were detected based on signal-to-noise ratios ( Additional file 3). It is worth noting that glucose-6-phosphate content was similar in fasted or “diabetic” vs. fed status, despite a large range of plasma glucose levels (232–747 mg/100 ml). A total of 12 metabolites were significantly different between treatment groups based on p < 0.05 and an additional five were suggestive of significance (p < 0.10; Table 5). Tissue levels of amino acids were consistently lower in fasted vs. fed tissue, with statistically significant reductions in asparagine and glutamine (p < 0.01). Presumably, these effects were due to a change in the balance of protein synthesis/proteolysis and to the catabolism of carbon skeletons for energy in response to energy restriction, which is consistent with up-regulated expression of genes involved in amino acid catabolism (Figure 3). They may also reflect a decrease in plasma amino acid supply as suggested by the decrease in total plasma amino acid levels (evaluated by the content in total α-NH₂-non-protein nitrogen (αNH₂NPN), i.e., mostly total amino acids), as compared to fed controls 18 . In contrast to fasting, tissue amino acid levels tended to be increased in insulin-neutralized vs. fed, although only glutamine showed a statistically significant response. Comparison of insulin-neutralized vs. fasted chickens highlights the divergent effects of treatments on amino acids (Figure 4). Alanine, arginine, asparagine, glutamine, histidine, proline, serine, threonine and tyrosine levels were all significantly higher in insulin-neutralized vs. fasted, with differences ranging from 1.7- to 3.4-fold. Two metabolites related to glucose metabolism, D-glucono-1,5-lactone-6-phosphate and glycerol-3-phosphate, were lower in both fasted and insulin-neutralized treatments vs. fed, with the latter comparison nearing statistical significance (p < 0.1). D-glucono-1,5-lactone-6-phosphate is a product of glucose-6-phosphate dehydrogenase (G6PDH), an enzyme that, in mammals is insulin-sensitive and rate-limiting for pentose phosphate pathway activity and production of cellular NADPH, an important cofactor for lipid metabolism 28 . However, pentose phosphate pathway activity is intrinsically low in chicken and is not stimulated when lipogenesis is high; the production of cellular NADPH is more closely related to malic enzyme activity 29 . Glycerol-3-phosphate is a product of both glucose and pyruvate metabolism and is used in triacylglycerol synthesis. Lower levels with both treatments may reflect glycerol demand for fatty acid reesterification in light of the apparent increase in lipolysis in both treatment groups 18 .

* p < 0.10, ** p < 0.05, *** p < 0.01, **** p < 0.001.

Correlated patterns of gene expression and metabolite abundance were extracted using hierarchical clustering to interconnect treatment effects on transcripts and metabolites. Clusters 2 and 3 contained genes and metabolites with lower abundance in fasted vs. fed or insulin neutralized tissue. The two clusters differed with respect to the insulin neutralized group: cluster 3 contained analytes at intermediate levels between fasted and fed, while cluster 2 contained those at levels comparable to or greater than fed. Twelve of the 17 metabolites with statistically suggestive or significant effects of treatment, including all of the amino acids and amino acid derivatives ( Additional file 4), were present in cluster 2 along with a set of genes that included the p85α regulatory subunit of PI3 kinase (PIK3R1), as well as ME, malonyl CoA decarboxylase and ELOVL6. Cluster 3 contained several metabolites including both NAD + and NADPH and was significantly enriched in GO annotations related to carbohydrate metabolism and in the KEGG pathways TCA cycle, glycolysis/gluconeogenesis, pyruvate metabolism and steroid biosynthesis (p < 0.05; Additional file 4). Clusters 7 and 8 consisted of genes and metabolites with higher levels in fasted than in the other two treatment groups. These clusters were significantly enriched in GO categories PPAR signaling and negative regulation of cellular biosynthesis and also contained citrate and pyruvate (p < 0.05; Additional file 4).

Despite roles as both a domestic food animal of worldwide economic importance and a widely used model organism with relevance for human obesity and insulin resistance, few studies have examined regulation of gene expression in chicken adipose tissue. To our knowledge, no studies of nutritional regulation of chicken adipose tissue at the genomic level have been reported in the published literature. Likewise, although insulin is the most well-defined hormonal mediator of metabolism in mammalian adipose tissue, its role in chicken remains to be clarified. Therefore the current study addressed two objectives: 1) characterize the transcriptomic and metabolomic response to energy manipulation as a step toward enhanced understanding of adipose biology in chicken; and 2) identify the effects of insulin on chicken adipose tissue by including a group of birds in which insulin action was blocked by immunoneutralization with an anti-insulin antibody. We sought to both identify potential new targets for genetic selection or management strategies to reduce fat accumulation in commercial broilers and to further develop chicken as a model organism for studies of human obesity.

Although intrinsic lipogenic activity is low in chicken adipose tissue, genes involved in fatty acid synthesis and storage were suppressed and those in fatty acid mobilization and oxidation were up-regulated by fasting. The 40 down-regulated genes with fold changes greater than three were significantly enriched for the GO annotation lipid biosynthetic process (FDR <0.05), including genes that control triglyceride synthesis (DGAT2 and AGPAT9) and fatty acid synthesis (ACACA, ACLY and ME), elongation (ELOVL6), and desaturation (FADS1). AGPAT9 and DGAT2 catalyze the initial and final steps, respectively, of de novo triglyceride synthesis. ACLY is the main enzyme for synthesis of cytosolic acetyl-CoA, which is carboxylated to malonyl-CoA by ACACA, the rate-limiting step in fatty acid synthesis. Reducing equivalents for the conversion of malonyl-CoA to palmitate are supplied by malic enzyme (ME). ELOVL6 catalyzes elongation of palmitate to stearate and appears to play a key role in insulin sensitivity 30 31 . Finally, FADS1 is rate-limiting for polyunsaturated fatty acids (PUFA) biosynthesis and was recently implicated in control of fasting glucose homeostasis in humans 32 . Genes altered by fasting in adipose tissue in this study overlapped with those shown to be differentially expressed in chicken liver after 16 or 48 hours of fasting, including ACLY, ACOX1, BCAT1 and PDK4 33 . These authors used a different array platform than ours, which precludes precise quantitative comparisons. However, among the genes changed in both studies, the fold changes observed in adipose tissue were consistently greater than those in liver, despite the longer duration of fasting in that study. For example, PDK4 expression was up-regulated ~ 18-fold by a five hour fast in adipose tissue, but only ~ 1.5-fold after a 16 hour fast in liver. While differences in sensitivity between the two array platforms must be kept in mind, these data suggest that adipose tissue metabolism in chicken is at least as sensitive to energy status as hepatic metabolism. Our results indicate that both fatty acid synthesis and storage are dynamically regulated by energy status in chicken adipose tissue, despite its modest (~ 15%) contribution to the amount of stored fatty acids.

Both fasted and insulin-neutralized birds exhibited significant increases in plasma glucagon. Parallel elevations in plasma NEFA suggested that this resulted in significant lipolysis of stored triacylglycerol in both treatment groups. During fasting, a considerable percentage of the liberated fatty acids are re-esterified in adipocytes, and only a small fraction traditionally have been thought to be oxidized in the mitochondria of adipocytes through beta oxidation 34 . However, recent studies in mice and in human adipose tissue demonstrate that in some conditions fatty acid oxidation in white adipose tissue is considerable and may be an important determinant of obesity 35 36 37 . Consistent with this concept, we found significant increases in a number of key enzymes that mediate mobilization of fatty acids and their oxidation, including the rate-limiting enzymes in both mitochondrial and peroxisomal fatty acid oxidation (CPT1A and ACOX1, respectively). We measured tissue levels of beta-hydroxybutyrate, a ketone product of beta oxidation, to confirm that changes in gene expression had functional consequences and found them to be significantly elevated in adipose tissue of fasted vs. fed chickens. Levels were numerically but not statistically higher in insulin-neutralized adipose tissue (data not shown). Qualitatively, fasting-induced changes in gene expression resemble those induced by the fibrate class of drugs, which activate PPARα and promote fatty acid oxidation in white adipose tissue and are used clinically to treat hyperlipidemia 35 37 38 39 40 . These data suggest that dietary activation of PPARα, for example through supplementation with fatty acids that preferentially bind and activate this member of the PPAR family 41 , may be a means to attenuate fat deposition in commercial broilers. Such action may underlie the reduced abdominal fat mass reported in broilers that were fed diets rich in n-3 PUFA 42 .

Both fasting and insulin neutralization elicited marked upregulation of PDK4. PDK4 is a nutrient sensing fuel switch that phosphorylates and inactivates pyruvate dehydrogenase, which shifts fuel use from glucose to fatty acids and spares glucose for the brain during periods of fasting 23 43 44 45 . PDK4 also enhances glycerol synthesis in white adipose tissue by shunting pyruvate into glyceroneogenesis, at least in the fed state 46 . Hepatic and skeletal muscle expression of PDK4 is increased by fatty acids, acetyl CoA, NADH and the diabetic state and decreased by insulin and pyruvate (rev. in 23 47 ). Little is known about PDK4 in chicken, but a recent study suggests it acts as a glycogen sensor in muscle and thus plays comparable roles to those in mammals 48 . In mouse white adipose tissue, PDK4 expression was shown to be induced by activation of p38MAPK 49 , which we found to be significantly up-regulated with fasting and, to a lesser extent, with insulin neutralization 19 . Although PDK4 was up-regulated in both treatment groups, and both groups showed evidence of increased lipolysis 18 , only fasted chickens presented a gene expression signature and tissue beta-hydroxybutyrate levels that were clearly indicative of fatty acid oxidation. Although we did not measure malonyl-CoA levels, we predict that they were reduced with fasting, but not insulin neutralization, based on reduced expression of ACACA. Malonyl-CoA allosterically binds and inhibits CPT1A, minimizing fatty acid transport and subsequent oxidation in mitochondria 50 . With insulin neutralization, increased PDK4 may thus be more aligned with the demand for glycerol needed to re-esterify fatty acids liberated by lipolysis 46 . Additional experiments are needed to confirm that manipulation of PDK4 alters fatty acid oxidation in chicken adipose tissue and to delineate its relative contributions to fatty acid oxidation and glyceroneogenesis under varying metabolic states. If manipulation of PDK4 does alter fatty acid oxidation, our results highlight this pathway as a potential target for reducing fatness, which has relevance for both poultry and humans.

Microarray data indicate that the effects of fasting in chicken adipose tissue extend beyond metabolism. GO analysis highlighted pathways such as cell cycle and cytokine-cytokine receptor interaction that are most likely related to changes in the stromal vascular fraction, which contains proliferating preadipocytes and cells of the immune system. In particular, a number of genes that regulate multiple steps in adipogenesis were significantly altered by fasting. Chickens rapidly accumulate abdominal fat after hatch, and until approximately 7 weeks of age this is due more to formation of new adipocytes than to adipocyte hypertrophy 11 . Adipocytes arise from mesenchymal stem cells in a two stage process of lineage commitment to an adipocyte fate, followed by differentiation of fibroblast-like preadipocytes into mature fat-storing cells 51 . Members of both the Wnt (MSC lineage commitment) and TGFβ/BMP (MSC lineage commitment and preadipocyte competence) signaling pathways were significantly regulated by fasting. Fasting down-regulated expression of CEBPα and PPARγ, two transcription factors that orchestrate the cascade of gene expression changes that lead to terminal adipocyte differentiation 51 . Expression of other adipogenic mediators including fibroblast growth factor 2 (FGF2), fibroblast growth factor receptor 1 (FGFR1) 52 , and nuclear receptor corepressor 1 (NCOR1) 53 were also significantly regulated by fasting. Collectively, these changes suggest that adipocyte number in chickens is dynamically tied to energy status, at least in young chicks (such as those used herein) that are rapidly forming new adipocytes. An elegant study by Arner et al. concluded that adipocyte number in humans is a major determinant of adult fat mass and is determined during early childhood 54 . Less is known about this process in humans due to the limitations of sampling adipose tissue, particularly during development and from different abdominal depots. In light of what appears to be sensitive regulation of adipogenesis by nutritional state, chickens may thus be particularly valuable models in which to elucidate mechanisms of adipocyte hyperplasia during development that would inform the study of human obesity.

It is worth noting that, despite the uncertainty about insulin signaling in chicken adipose tissue, fasting altered the expression of several messengers encoding elements of the insulin signaling cascade. Expression of PIK3CB, which encodes the catalytic p110 subunit of PI3K, was up-regulated with fasting, while PIK3R1, which encodes the regulatory p85 subunit, was down-regulated. Such regulation could maintain some insulin signals despite a fall in plasma insulin level. CBLB and CRK, which mediate insulin signals that are associated with lipid rafts 55 , were also up-regulated with fasting. In mammals, this pathway stimulates glucose uptake independently of PI3K activation, which may shed light on the apparent refractoriness of PI3K activity to insulin that was described in chicken skeletal muscle 18 . Therefore, the potential effects of insulin on lipid storage and energy utilization appear to be defended in the fasting state, when insulin levels fall, by enhanced insulin sensitivity at the post-receptor level. Additional studies are needed to confirm this effect and to further explore the potential of PI3K-independent effects of insulin on glucose utilization in chicken adipose tissue.

Insulin is not considered to be a key regulator of glucose metabolism in chicken adipose tissue, although it does induce glucose disposal in chicken liver and muscle 14 . It is therefore not surprising that the majority of genes significantly altered by both insulin neutralization and fasting are not related to glucose metabolism and lipid synthesis. The main exception is DGAT2, which catalyzes the final step in esterification of fatty acids into triglycerides. In fact, DGAT2 showed the most extreme down-regulation (6.1- and 10.5-fold, insulin-neutralized and fasted, respectively) in each treatment group, which is surprising considering that other genes related to lipogenesis were only regulated by fasting. Suppression of DGAT2 expression may be due to feedback by lipolysis, which appeared to be increased in both treatment groups based on plasma NEFA levels. In general, our data indicate that insulin deprivation altered fatty acid and glucose metabolism in a manner comparable to fasting but to a lesser extent, such that most genes involved in these pathways did not exhibit statistically significant changes in expression. For example, cluster analysis (Figure 2) revealed that some genes upregulated by fasting were also increased by insulin neutralization (clusters 1, 3 and 4); these three clusters were enriched with genes in the KEGG pathways for fatty acid metabolism and PPAR signaling, including both ACOX1 and CPT1A, among others. Similarly, among genes that were downregulated by fasting, clustering discriminated a set of genes (Figure 2, cluster 2) with a trend to also be decreased (albeit to a lesser extent than in fasting) by insulin deprivation. Interestingly, this cluster was significantly enriched in functions related to carbohydrate metabolism, suggesting that insulin does play some role in chicken adipose glucose metabolism. Comparable trends appeared in the metabolomic data. For example, stearate and palmitate (the only fatty acids measured by our MS platform) were lower (although not significantly) in both fasted and insulin neutralized compared to fed birds ( Additional file 3). While the purpose of our study design was to determine the specific effects of insulin on chicken adipose tissue, we cannot exclude the possibility that some of the overlapping changes in gene expression were secondary to systemic factors, such as hyperglucagonemia present in both treatment groups 18 . In vitro experiments using primary adipocytes or adipose explants will be useful to confirm specific effects of insulin on genes identified herein.

Of the 13 changes in expression that were unique to insulin neutralization, the most interesting responses were up-regulation of GCG, which encodes preproglucagon (fold change = 2.91), and down-regulation of the glucagon receptor (LOC425670, fold change = −2.77). The proglucagon system in avians is more complex than in mammals. The avian preproglucagon locus encodes two distinct precursor proteins that yield different peptides through alternative posttranslational processing: the class A transcript (PGA) yields glucagon and glucagon-like peptide-1 (GLP-1), while the class B transcript (PGB) additionally produces glucagon-like peptide-2 (GLP-2) and is more like the mammalian transcript (rev. in 56 ). Adipose tissue expresses both transcripts, with PGA being slightly more abundant, and is the third highest preproglucagon expressing tissue in chicken, behind pancreas and the proventriculus 57 . We used transcript-specific QPCR to determine that only the PGB transcript was up-regulated by insulin neutralization (data not shown). Additional experiments are necessary to delineate which of the encoded peptides are up-regulated in parallel, but the coincident down-regulation of the glucagon receptor suggests a paracrine glucagon axis in chicken adipose tissue, and one that is regulated by insulin. In support of this concept, plasma glucagon (presumably derived largely from pancreas) was elevated comparably in both treatment groups 18 , while GCG expression in adipose tissue was only up-regulated by insulin neutralization.

Tissue metabolomic analysis highlighted effects of insulin neutralization that were divergent from fasting and not readily apparent from microarray data. Most of the tissue amino acids that were measured were higher with insulin-neutralization but lower with fasting when each group was compared to ad libitum fed controls. This pattern parallels the levels of αNH₂NPN levels in blood 18 . Low levels in fasted adipose tissue were most likely due to oxidation of the carbon skeletons for cellular energy through the tricarboxylic acid cycle (TCA) cycle and/or for glyceroneogenesis, in the absence of dietary glucose. Increased amino acid catabolism was reflected in the differential expression profiles of the fasted vs. fed comparison (Figure 3; Additional file 1). In the insulin neutralized group, however, glucose supply from food was maintained and preferentially oxidized for energy. Elevated amino acids in the insulin neutralized group may also reflect reduced utilization due to the lack of insulin’s anabolic effects, particularly on the proliferating cell population within adipose tissue. The metabolomics approach used here measured only metabolite pool sizes at the time that tissues were harvested, rather than the effect of fasting or insulin neutralization on the rates of metabolism through glycolysis and the TCA cycle. The latter would be much more informative with respect to the dynamic impact of treatment, but requires the use of isotopic labeling (e.g., by feeding ¹³C-labelled glucose) which was not performed in this study. Nonetheless, we were able to demonstrate significant effects on some metabolites that inform the parallel changes in gene expression, particularly in relation to amino acid metabolism. Combined clustering of metabolomic and gene expression together identified a set of genes correlated with many amino acid levels, including PIK3R1, ME and MCD.

In summary, we determined that adipose tissue metabolism in the chicken is regulated by energy status and, to a lesser extent by insulin. Although adipose tissue is not a primary site of lipogenesis in chicken, the rate-limiting genes for fatty acid synthesis were suppressed by fasting. Likewise, fasting appeared to increase aspects of insulin sensitivity based on expression profiles, despite the view that chicken adipose tissue is relatively insensitive to insulin. Consistent with this paradigm, insulin neutralization significantly altered the expression of only a few genes related to glucose and lipid metabolism. Nonetheless, a considerable number of genes were altered by insulin neutralization, many of which thus far have unclear roles in adipose biology. Expression profiles suggest that even short term fasting alters fat storage in broilers by enhancing the oxidation of fatty acids. The initiating events that trigger upregulation of the corresponding genes are unclear, but there is considerable evidence for activation of PPARa, LXRa, and potentially other transcription factors that are activated by fatty acid ligands. Further studies are warranted to identify these triggers because of their potential impact on fat storage. Our data also suggest that broiler chicks may be an informative model organism in which to investigate dietary effects on adipose development in light of what appears to be a relationship between energy intake and adipogenesis. The results of this study thus have dual benefit for both the poultry industry and for studies of obesity in humans.

VLDL: Very low density lipoprotein; CEBPα :CCAAT: Enhancer-binding protein alpha; PPAR: Peroxisome proliferator-activated receptor; IRβ: Insulin receptor β-subunit; IRS-1: Insulin receptor substrate-1; Shc: Src homology 2 domain-containing substrate; PI3K: Phosphatidylinositol 3-kinase; QPCR: Quantitative PCR; FDR: False discovery rate; GO: Gene ontology; PDGF: Platelet derived growth factor; MSC: Mesenchymal stem cell; Camp: Cyclic adenosine monophosphate; DGAT2: Diacylglycerol O-acyltransferase homolog 2; PDK4: Dehydrogenase kinase 4; ACBD5: Acyl-Coenzyme A binding domain containing 5; AGTR1: Type I angiotensin II receptor; PUFA: Polyunsaturated fatty acids; NEFA: Non-esterified fatty acid; ATGL: Adipose triglyceride lipase; PD: Pyruvate dehydrogenase; LXR: Liver X receptor; GCG: Preproglucagon; PGB: Preproglucagon class B transcript; BCAAs: Branched chain amino acids.

The authors declare that they have no competing interests.

BJ performed the experimental assays, analyzed genomic data and drafted the manuscript. BE, JG and SC generated and analyzed metabolomics data. SD provided technical support and AM provided statistical support. JS, SM-C and JD performed the original in vivo study and provided tissue samples. BV designed and supervised the study and personnel and edited the manuscript. All authors read and approved the final manuscript.