- Research article
- Open Access
Transcriptional analysis of abdominal fat in genetically fat and lean chickens reveals adipokines, lipogenic genes and a link between hemostasis and leanness
BMC Genomics volume 14, Article number: 557 (2013)
This descriptive study of the abdominal fat transcriptome takes advantage of two experimental lines of meat-type chickens (Gallus domesticus), which were selected over seven generations for a large difference in abdominal (visceral) fatness. At the age of selection (9 wk), the fat line (FL) and lean line (LL) chickens exhibit a 2.5-fold difference in abdominal fat weight, while their feed intake and body weight are similar. These unique avian models were originally created to unravel genetic and endocrine regulation of adiposity and lipogenesis in meat-type chickens. The Del-Mar 14K Chicken Integrated Systems microarray was used for a time-course analysis of gene expression in abdominal fat of FL and LL chickens during juvenile development (1–11 weeks of age).
Microarray analysis of abdominal fat in FL and LL chickens revealed 131 differentially expressed (DE) genes (FDR≤0.05) as the main effect of genotype, 254 DE genes as an interaction of age and genotype and 3,195 DE genes (FDR≤0.01) as the main effect of age. The most notable discoveries in the abdominal fat transcriptome were higher expression of many genes involved in blood coagulation in the LL and up-regulation of numerous adipogenic and lipogenic genes in FL chickens. Many of these DE genes belong to pathways controlling the synthesis, metabolism and transport of lipids or endocrine signaling pathways activated by adipokines, retinoid and thyroid hormones.
The present study provides a dynamic view of differential gene transcription in abdominal fat of chickens genetically selected for fatness (FL) or leanness (LL). Remarkably, the LL chickens over-express a large number of hemostatic genes that could be involved in proteolytic processing of adipokines and endocrine factors, which contribute to their higher lipolysis and export of stored lipids. Some of these changes are already present at 1 week of age before the divergence in fatness. In contrast, the FL chickens have enhanced expression of numerous lipogenic genes mainly after onset of divergence, presumably directed by multiple transcription factors. This transcriptional analysis shows that abdominal fat of the chicken serves a dual function as both an endocrine organ and an active metabolic tissue, which could play a more significant role in lipogenesis than previously thought.
The chicken was the first avian species and domestic animal selected for complete genome sequencing and assembly . Subsequently, the chicken has emerged as a premier model in animal agriculture [2–4] and developmental biology . Although now recognized as a model organism for biomedical research , the chicken has not been extensively used for the study of human diseases, especially metabolic disorders (i.e., insulin resistance, diabetes, obesity and metabolic syndrome). Several unique features of avian metabolism make the chicken an interesting model for understanding the interactions between genetic and endocrine factors that contribute to development of obesity and related metabolic disorders. In particular, chickens normally exhibit “hyperglycemia” [7, 8], insulin resistance [8–11], hepatic de novo synthesis of lipids  and, like humans , abdominal (visceral) fatness is a polygenic trait [14–19]. Despite their relative insensitivity to insulin, acute immunoneutralization of insulin in the chicken provokes differential expression of more than a thousand genes in both liver and in skeletal muscle . In contrast, only 69 genes were differentially expressed (DE) in abdominal fat of chickens following insulin immunoneutralization, albeit short-term fasting produced a much larger change (1780 DE genes) in transcription of abdominal fat genes . This recent work also shows a rather large decrease in expression of lipogenic genes in abdominal fat of fasted chickens. A detailed examination of the insulin signaling cascade in adipose tissue of the chicken shows a distinct unresponsiveness to insulin . Collectively, these observations support the chicken as a unique model for the study of the genetic and biological mechanisms controlling fatness or leanness.
Most mammalian models of obesity exploit single gene mutations or use high-energy, high-fat diets to induce obesity . Our chicken models are two experimental lines of meat-type chickens that were divergently selected over seven generations for either high (FL) or low (LL) abdominal (visceral) fatness [24, 25]. These chickens exhibit a 2.5-fold difference in abdominal fat weight at 9 weeks (wk) of age, albeit their body weight and feed intake are similar . Furthermore, the FL chickens present hyperplasia and hypertrophy of adipocytes at an earlier age than do LL chickens [27, 28].
Differential abundance of lipogenic genes in liver of the FL and LL chickens was determined earlier by differential mRNA display , quantitative RT-PCR [30, 31] or targeted low-density array . Our preliminary analysis of the liver transcriptome in the FL and LL chickens during juvenile development revealed 1,805 differentially expressed (DE) genes . Quantitative trait loci (QTL) analyses of an FL x LL intercross identified a major QTL for abdominal fatness at the distal end of chromosome 5 (GGA 5) [16, 17, 33]. Further, the expression quantitative trait loci (eQTL) analysis of GGA 5, involving a three generation intercross of the FL x LL chickens, identified variations in expression of 660 hepatic genes that were correlated with abdominal fatness traits .
The present study has a dual purpose to explore the abdominal fat transcriptome of juvenile FL and LL chickens and to identify major gene networks controlling adiposity and lipogenesis in these divergently selected models. Using the Del-Mar 14K Chicken Integrated Systems cDNA microarray, we took transcriptional snapshots of gene expression in abdominal fat across two genotypes (four birds/genotype) and six ages during juvenile development (1–11 weeks of age). Interestingly, our time-course transcriptional analysis of abdominal fat revealed numerous DE genes (main effect of genotype or age × genotype interaction) that are involved in hemostasis (14 genes), adipokine signaling (8 genes), retinol metabolism (13 genes), and the synthesis (37 genes), oxidation (12 genes) and transport (12 genes) of lipids. The liver is widely considered as the major site of lipogenesis in chickens and other birds. However, the present transcriptional analysis of visceral adiposity has identified 37 lipogenic DE genes, including FASN, SCD, SREBF1, SREBF2 and THRSPA that are expressed higher in FL chickens. The greater abundance of thrombogenic enzymes and related protease inhibitors in abdominal fat of the LL chickens suggests enhanced proteolytic processing of adipokines and other endocrine factors, with local and/or humoral actions, that could contribute to their leaner phenotype. Although abdominal fat is generally considered as a passive depot for lipids, the present descriptive study in FL and LL chickens supports our idea that it does contribute to lipid synthesis and serves as an endocrine organ, which liberates a host of adipokines and endocrine factors with intrinsic and/or extrinsic activity.
Animals and tissue collection
The birds were bred and raised at the Institut National de la Recherche Agronomique (INRA) UE1295 Pôle d’Expérimentation Avicole de Tours, F-37380 Nouzilly, France. At hatching, FL and LL cockerels were wing-banded and vaccinated against Marek’s disease virus. Birds were reared together in floor pens (4.4 × 3.9 m) and provided ad libitum access to water and conventional starter feed for three weeks [3,050 kcal of metabolizable energy (ME)/kg and 22% protein] and thereafter with a grower ration (3,025 kcal ME/kg and 17.9% protein). Chicks were held under continuous light (24 h or LL) for the first two days after hatching, followed by a 14 h light/10 h dark cycle (14L:10D) for the remainder of the experiment. Infrared gas heaters provided supplemental heat and ambient temperature was decreased weekly from 32 C at hatching until 22 C was reached at 3 wk of age. Eight birds from each genotype were randomly selected at six ages (1, 3, 5, 7, 9, and 11 wk), weighed, bled into heparinized syringes, and killed by cervical dislocation. Abdominal fat was quickly dissected and weighed; a sample was immediately snap frozen in liquid nitrogen and stored at −75 C until further processing. All animal procedures were performed under the strict supervision of a French government veterinarian and in accordance with protocols approved by the French Agricultural Agency, the Scientific Research Agency, and the Institutional Animal Care and Use Committees at INRA, Nouzilly, France. These procedures were also in compliance with the United States Department of Agriculture (USDA) guidelines on the use of agricultural animals in research and approved by the University of Delaware Agricultural Animal Care and Use Committee.
Four birds per genotype and age were randomly selected from the total of eight birds sampled per genotype and age for microarray analysis of abdominal fat (Additional file 1). Total cellular RNA was extracted from abdominal fat using guanidine thiocyanate and CsCl gradient purification , followed by a separate step for DNase I treatment. The RNA concentration was determined with a NanoDrop ND-1000 spectrophotometer (NanoDrop Technologies; Wilmington, DE). RNA integrity was examined using an RNA 6000 Nano Assay kit and the Model 2100 Bioanalyzer (Agilent Technologies; Palo Alto, CA) to assess the quality of the RNA samples (RNA integrity number, RIN ≥ 9 was considered acceptable). Twenty μg of total RNA was indirectly labeled using SuperScript Plus Indirect cDNA Labeling System (Invitrogen, Carlsbad, CA). First strand cDNA synthesis was performed in a 30 μl final volume containing 1× first-strand buffer, 5 μg of anchored oligo(dT20), DTT, dNTP mix (including aminoallyl- and aminohexyl-modified nucleotides), 40 U of RNaseOUT and 800 U of SuperScript III reverse transcriptase with an incubation at 46 C for 3 h. The original RNA template was removed by NaOH hydrolysis, and followed by neutralization with HCl. The cDNA was purified using a low-elution volume spin cartridge (Invitrogen; Carlsbad, CA) and labeled with either Alexa Fluor® 555 or Alexa Fluor® 647 succinimidyl ester in the dark at room temperature for 2 h. After purification of labeled cDNA with a low-elution-volume spin cartridge, the efficiency of dye incorporation was determined using the Microarray Module on the NanoDrop ND-1000 spectrophotometer and the Base:Dye Ratio Calculator on the Invitrogen website .
Twenty-four Del-Mar 14K Chicken Integrated Systems microarrays (NCBI GEO Platform # GLP1731) were hybridized with 48 labeled samples using a balanced block design, where half of the birds from each genotype and age were labeled with Alexa Fluor® 647 (red dye) and the other half with Alexa Fluor® 555 (green dye; see Additional file 1 for details of the hybridization design). Hybridized slides were scanned using a GenePix 4000B scanner with GenePix Pro 4.1 software (Molecular Devices, Union City, CA) at wavelengths of 635 nm (Alexa Fluor® 647-labeling) and 532 nm (Alexa Fluor® 555-labeling) generating a combined TIFF image file for each slide. The laser power was set at 100% with the photomultiplier tube (PMT) setting being adjusted for each scan to produce a PMT count near unity. All slides were manually checked for quality and all spots with inadequacies in signal, background or morphology were eliminated from further analysis. The image analysis results were merged with Excel files in GenePix Report (GPR) format, which contains clone identification, spot location on slide, and most current gene name/function (based on BLASTX/BLASTN score).
The microarray GPR files were analyzed using the linear models for statistical analysis of microarray data (LIMMA, version 3.4.5) software  package in R (version 2.11.1) . Median intensities for each dye were Loess normalized (without background subtraction) within array and between array (“Aquantile” method) to correct for dye and slide biases. A two-way analysis of variation (ANOVA) was used on Loess normalized intensity values from this factorial design experiment to determine the main effect of genotype (G), main effect of age (A), and the interaction of age and genotype (A × G). The Benjamini-Hochberg procedure  was used to control the experiment-wise false discovery rate (FDR) from multiple testing procedures.
Quantitative RT-PCR analysis
Several DE and prior candidate genes were selected for verification of expression by quantitative RT-PCR (qRT-PCR) analysis. First-strand cDNA synthesis was performed by incubation of a 13 μl reaction volume (containing 1 μg of total DNase-treated RNA, 1 μl of 100 μM oligo dT20, 1 μl of 10 mM dNTP mix, and water to 13 μl total volume) for 5 min at 70 C and then placed on ice for 2 min. A master mix containing 5 μl of 5× first-strand synthesis buffer, 1 μl of 0.1 M dithiothreitol (DTT), 1 μl of RNaseOUT, and 200 U of SuperScript® III reverse transcriptase (Invitrogen, Carlsbad, CA) was added to the RNA in a final reaction volume of 20 μl. The cDNA was diluted to achieve a concentration of 50 ng/μl. Primers were designed for qRT-PCR using Primer Express® v2.0 software (Applied Biosystems, Foster City, CA). Detailed information for each primer pair including gene name, gene symbol, primer sequences (forward and reverse), GenBank accession number and amplicon size are provided in Additional file 2.
The qRT-PCR assay was performed in an ABI Prism Sequence Detection System 7900HT, using Power SYBR® green PCR master mix (Applied Biosystems, Foster City, CA) and 400 nM of each primer (forward and reverse; Sigma-Aldrich, St. Louis, MO) in duplicate wells. Disassociation curves of each sample were analyzed to validate specific amplification and verify absence of primer dimers. PCR products were analyzed using agarose gel electrophoresis to compare approximate product size to expected amplicon size. The Ct for each sample was normalized to the corresponding sample geometric mean of three housekeeping genes [protein kinase, AMP-activated, beta 2 non-catalytic subunit (PRKAB2), protein kinase, AMP-activated, gamma 1 non-catalytic subunit (PRKAG1), and serpin peptidase inhibitor, clade E (nexin, plasminogen activator inhibitor type 1), member 2 (PAI-2 or SERPINE2)]. These housekeeping genes were selected using the RefFinder website  as the most stably expressed genes (i.e., genes with the lowest M-value) in the experiment. The 2-(∆∆Ct) formula was used to calculate relative transcript abundance . The statistical analysis was performed using a general linear model procedure in SAS v9.3. The data (log2 transformed normalized expression values) was analyzed using a two-factor analysis of variance to determine significant effects of genotype (G), age (A), and the interaction of age x genotype (A × G). Pearson’s correlation coefficient (r) was used to compare log2 FL/LL expression ratios between the microarray and qRT-PCR analyses of select genes.
Body weight (BW, kg), abdominal fat weight (g), and relative abdominal fat content (percent of body weight, %BW) in juvenile FL and LL chickens are presented in Table 1. The BW of FL and LL cockerels was similar for all ages between 1 and 11 wk. The absolute and relative abdominal fat weights of the FL chickens were 2.5-fold higher (P≤0.01) than those of the LL at all ages between 3 and 11 wk of age.
Abdominal fat gene expression
Differentially expressed (DE) genes were defined as those having a significant false discovery rate (FDR)-adjusted P-value. The significance level was set at P≤0.05 for genotype (G) and the age (A) × G interaction; or P≤0.001 for A. The statistical analysis of this factorial design (2 genotypes × 6 ages) experiment provided DE gene sets for the main effects of G and A, or the interaction of A × G. The main effect of G was determined by comparing gene expression values of each genotype (FL vs. LL) averaged across the six juvenile ages (1–11 wk). Likewise, the main effect of age (A) was determined by comparing gene expression values of each age averaged across both genotypes. To distinguish differences between ages, five single-degree-of-freedom contrasts were made by comparing the average of each subsequent age (averaged across both genotypes) against the 1 wk average (1 vs. 3 wk, 1 vs. 5 wk, 1 vs. 7 wk, 1 vs. 9 wk and 1 vs. 11 wk). The Venn diagram (Figure 1) shows the overall number of DE genes for G (344 genes), A (3,223 genes) and the interaction of A × G (254 genes) and the intersection of these DE gene sets. The number of unique genes are indicated for G (119 DE genes), A (3,183 DE genes) and the A × G interaction (32 DE genes). There were 213 DE genes in common between G and the A × G interaction. Thirty-one DE genes were shared between A and G, whereas 28 DE genes were in common between A and the A × G interaction. Overall, 19 DE genes were found in common among all three effects.
All DE genes involved in the higher order A × G interaction were removed from the main effects of A and G lists. Therefore, the total number of DE genes presented in the A × G interaction list (254 DE genes) in Additional file 3 reflects the 32 unique genes for the A × G interaction plus those genes shared with the main effects of A (19 DE genes) and G (194 DE genes). The number of DE genes presented in the main effect of A (3,195 DE genes) and A (131 DE genes) lists represent the total number of DE genes for that effect minus those genes that belong to the A × G interaction (Additional file 3).
Ingenuity® Pathway Analysis (IPA®) of DE gene sets
Significant genes (cDNA clone IDs) from the microarray analysis were annotated using the GeneBase tool on our website , which provides protein IDs (from GenBank or Swiss-Prot databases) of microarray cDNA probes derived from BLASTX analysis. Lists of DE genes containing the protein ID and log2 ratio for each gene were then submitted to the Ingenuity® Knowledge Base  for functional annotation and mapping to canonical metabolic and regulatory pathways. “Analysis ready” genes were mapped by IPA for the genotype (100 DE genes), age (2,301 DE genes), and age x genotype interaction (212 DE genes) lists. The IPA® Upstream Regulator Analysis was used to identify transcription factor (TF) interaction networks, predicted activation or inhibition of TF, and their direct targets from DE gene sets.
A summary of the IPA “Diseases and Disorders” category under “Biological Functions” is presented in Table 2. The subcategories of major interest were “Developmental Disorder” (33 genes), “Hereditary Disorder” (71 genes), “Inflammatory Disease” (7 genes, out of which 6 were up regulated in LL chickens), “Metabolic Disease” (41 genes), and “Organismal Injury and Abnormalities” (31 genes). A group of 33 genes were classified as “inborn error of metabolism” in three of the above subcategories (Developmental Disorder, Hereditary Disorder, and Metabolic Disease) [see Additional file 4].
One gene interaction network identified by IPA was heavily populated with a large number of hemostatic genes, which were up regulated in abdominal fat of the LL chickens (Figure 2). These genes are involved in coagulation [F2, A2M, carboxypeptidase B2 (CPB2), fibrinogen alpha (FGA), PLG, protein C (PROC) and serine peptidase inhibitor, clade D, member 1 (SERPIND1)] and inflammation [CD163 and retinoic acid receptor responder 2 (RARRES2) or chemerin]. Another group of DE genes [taste receptor, type 1, member 1 (TAS1R1), motilin receptor (MLNR), vasoactive intestinal peptide receptor 1 (VIPR1), and omega-3 fatty acid receptor 1 (O3FAR1)] are G-coupled receptors linked through the chemokine ligand CXCL12. Three genes shown in this pathway are related to steroid metabolism [hydroxysteroid (17-β) dehydrogenase 2 (HSD17B2) and hydroxysteroid (17-β) dehydrogenase 7 (HSD17B7)] and action [nuclear receptor subfamily 5, group A, member 1 (NR5A1)]. The transcription factor HNF1A regulates several hemostatic genes (FGA, PLG, PROC and SERPINA1) in visceral fat of LL chickens. PPARG directly regulates three adipokines (RARRES2, RBP4 and SERPINA1) expressed at higher abundance of LL and several additional genes up regulated in the FL (AGT, CTSL1, MCM7, SIM1, TGFBR1 and TUBB), including three metabolic enzymes (FASN, PDE3B and PYGL).
The higher expression of select hemostatic genes found in abdominal fat of LL chickens was verified by qRT-PCR analysis (Figure 3): serine proteases [F2, coagulation factor IX (F9) and protein C (PROC)] and protease inhibitors [A2M, annexin A5 (ANXA5), and SERPIND1]. Thrombin (F2) was more abundant in abdominal fat of the LL at all ages, except 11 wk. The expression of PROC was 3-fold higher in the LL at 1 and 3 wk, and over 10-fold higher at 5 wk. The coagulation factor F9 (Christmas factor) was over expressed in visceral fat of the LL by 3-fold, 24-fold, and 29-fold at 1, 5 and 7 wk, respectively. The expression patterns of two serine proteases (F9 and PROC) were similar with the greatest differences at 1 and 5 wk. The qRT-PCR analysis shows similar expression patterns between some hemostatic factors and adipokines (Figure 4). For example, the expression of ANXA5, F2, adiponectin (ADIPOQ), adiponectin receptor 1 (ADIPOR1) and attractin (ATRN) were highest in abdominal fat of the LL at 9 wk. Similarly, expression profiles of A2M, retinol binding protein 4 (RBP4) and angiopoietin-like 4 (ANGPTL4) were greatest in the LL at 7 wk. The adipokine visfatin [or nicotinamide phosphoribosyltransferase (NAMPT)] was not differentially expressed in adipose tissue of juvenile FL and LL. Both ADIPOQ and ANGPTL4 were identified in the main effect of age (A) by microarray analysis, although the log2 expression ratios were only slightly higher in the FL. The qRT-PCR analysis shows that the expression of ADIPOQ was higher (P≤0.05) in the LL between 7–11 wk of age, while the abundance of ANGPTL4 was elevated at 1, 5, 7 and 11 wk of age, albeit only age (A) produced a significant main effect.
The top canonical pathways identified by IPA (Additional file 5) reflect the prevalence of hemostatic genes in adipose tissue of LL chickens. The IPA software provided functional assignments of DE genes to “Coagulation System” (7 genes), “Acute Phase Response Signaling” (15 genes) and “Intrinsic Prothrombin Activation” (6 genes) pathways. These adipose genes include serine proteases [F2, PLG, PROC, and complement factor B (CFB)], protease inhibitors [A2M, serine peptidase inhibitor clade A member 1 (SERPINA1), and SERPIND1] and transporters of retinol [retinol binding protein 4 (RBP4) and 7 (RBP7)]. The IPA functional category “Lipid Metabolism” (Additional file 6) shows high representation of numerous genes involved in “oxidation of lipid” (11/12 genes higher in LL chickens), “transport of lipid” (9/12 genes higher in LL chickens), “synthesis of lipid” (18/37 genes up regulated in FL chickens) and “metabolism of retinoid” (5/5 genes higher in LL chickens).
Higher expression of lipogenic genes in adipose tissue of FL chickens
The abdominal fat of FL chickens exhibits higher expression of lipogenic transcription factors [sterol regulatory element binding transcription factor 1 (SREBF1), thyroid hormone responsive Spot 14 protein (THRSP) and sirtuin 2 (SIRT2)] (Figure 5A). In contrast, additional regulators of transcription [THRSP-like (THRSPL) or MID1 interacting protein 1 (MID1IP1); the nuclear liver X receptor-β (LXR β or NR1H2); and the proto-oncogene jun (JUN)] were more abundant in abdominal fat of the LL. As shown in this IPA gene interaction network, SREBF1 directly up regulates several genes in the FL that are involved in lipid biosynthesis [FASN, stearoyl-CoA desaturase (SCD), fatty acid desaturase 2 (FADS2), sterol-C5-desaturase (SC5DL), mevalonate decarboxylase (MVD), 7-dehydrocholesterol reductase (DHCR7), 3-hydroxy-3-methylglutaryl-CoA reductase (HMGCR) and lanosterol synthase (LSS)] and ketogenesis [3-hydroxy-3-methylglutaryl-CoA synthase 2 (HMGCS2)]. Some of these genes are also targets of SIRT2 and THRSPA and differentially expressed in adipose tissue of the FL. In addition, SREBF1 directly affects numerous genes that are expressed higher in the LL [fatty acid desaturase 1 (FADS1), acetyl-CoA carboxylase alpha (ACACA), acetoacetyl-CoA synthetase (AACS), farnesyl-diphosphate farnesyltransferase 1 (FDFT1), solute carrier family 2 (SLC2A2; or facilitated glucose transporter 2, GLUT2), succinate-CoA ligase, alpha subunit (SUCLG1), and phosphomevalonate kinase (PMVK)]. Two JUN targets, prostaglandin D2 synthase (PTGDS) and MID1IP1 (which regulates transcription of ACACA), were over-expressed in adipose tissue of the LL. Insulin-like growth factor binding protein 4 (IGFBP4) is another target of JUN that was expressed at higher levels in FL adipose tissue. The IPA Upstream Regulator Analysis predicts that JUN and SREBF1 lead to activation (indicated by orange arrows) of numerous up-regulated target genes (red symbols) in abdominal fat of the FL chickens (Figure 5B and 5C). In contrast, the blue blunted lines show predicted inhibition and down-regulation of DE genes by JUN (TP53) and SREBF1 (PCK1 and SLC2A2) in the LL (green symbols). The yellow arrows indicate that the IPA analysis would expect these targets to be activated by JUN and SREBF1, rather than down regulated as shown by these green gene symbols. The majority of the up-regulated genes found in abdominal fat of the FL (red symbols) are enzymes involved in lipid metabolism.
The expression profiles of eight genes mainly associated with lipid metabolism were examined by qRT-PCR analysis (Figure 6). A main effect (P≤0.05) of genotype (G) was observed for FASN (4-fold increase in FL at wk 7), SCD (4-fold and 3-fold increase in FL at wk 3 and 7, respectively), and pyruvate dehydrogenase kinase, isozyme 4 (PDK4, over expressed in LL chickens from 7 to 11 wk). A significant age by genotype (A x G) interaction (P≤0.05) was observed for facilitated glucose transporter, member 1 (GLUT1), perilipin 2 (PLIN2) and lipoprotein lipase (LPL). A main effect of age (A; P≤0.05) was also observed for FASN, GLUT1, PLIN2, PDK4, LPL, facilitated glucose transporter, member 8 (GLUT8) and superoxide dismutase 3 (SOD3).
Another network populated by numerous DE genes (up regulated in the FL) that control lipid metabolism shows the interaction of four transcription regulators (SIRT2, PPARD, EGR1 and CUX1), also up regulated in visceral fat of the FL (Figure 7A). Peroxisome proliferator-activated receptor delta (PPARD) interacts directly with patatin-like phospholipase domain containing 2 (PNPLA2), long chain acyl-CoA dehydrogenase (ACADL), aminoacylase 1 (ACY1), aldehyde dehydrogenase 2 (ALDH2), peroxiredoxin 6 (PRDX6), fatty acid binding protein 7 (FABP7), sorbitol dehydrogenase (SORD) and chemokine (C-C motif) ligand 13 (CCL13). Early growth response 1 (EGR1) interacts with CCL13 and 3-hydroxy-3-methylglutaryl CoA reductase (HMGCR), the rate-limiting enzyme in biosynthesis of cholesterol, which is a target of the histone deacetylase sirtuin 2 (SIRT2). The ketogenic enzyme 3-hydroxy-3-methylglutary-Coenzyme A synthase 2 (HMGCS2) is a downstream target of both SIRT2 and PPARD. Three additional metabolic enzymes [lanosterol synthase (LSS) and 7-dehydrocholesterol reductase (DHCR7)] were also expressed at higher levels in the FL, while acetoacetyl-CoA synthetase (AACS), fatty acid desaturase 1 (FADS1), and phosphomevalonate kinase (PMVK) were more abundant in abdominal fat of the LL. The Ingenuity Upstream Regulator Analysis identified several additional targets of PPARD (Figure 7B), which were expressed at higher levels in either the FL (FASN, FLT1, SCD and TLR5) or the LL (ACAA, APOA4, BCL2, GPD2, SLC27A1, SLC27A2, UCP3 and VLDLR) chickens. The IPA software predicts that PPARD is inhibited (blue color) based on prior knowledge of PPARD action in mammals and the observed higher expression of PPARD-activated targets in the LL (green symbols). IPA predicts that SIRT2 should be activated (orange color) and has a direct positive action on five target genes (DHCR7, HMGCR, LSS, MVD and SC5DL), which were up regulated in adipose tissue of the FL chickens. However, the yellow-colored arrows indicate that IPA expected three target genes (AACS, FDFT1 and PMVK) to be up regulated in the FL, rather than the LL as we observed.
Ligand activated nuclear receptors and other transcription factors
Of special interest are genes involved in ligand-activated gene transcription (e.g., retinol and thyroid hormone signaling) which regulate lipid metabolism (Table 3). Functional annotation of DE genes by IPA analysis identified five genes [alcohol dehydrogenase 1C (ADH1C), alcohol dehydrogenase 5 (ADH5), cytochrome P450, family 2, subfamily E, polypeptide 1 (CYP2E1), RARRES2 and RBP4] related to “metabolism of retinoid” (Additional file 6). An additional four retinol-related genes [RBP7, nucleolar protein 7 (NOL7), transthyretin (TTR) and retinol dehydrogenase 1 (RDH1)] were found by microarray analysis (Additional file 3).
The qRT-PCR analysis of six candidate transcription factors is presented in Additional file 7. Four genes [peroxisome proliferator-activated receptor alpha (PPARA), peroxisome proliferator-activated receptor gamma (PPARG), PPARD and sterol regulatory element binding transcription factor 2 (SREBF2) showed only a main effect of age (A). A main effect of genotype (G) was observed for SREBF1 due to higher expression in the FL at 1 and 5 wk. Similarly, the abundance of retinoid X receptor, gamma (RXRG) was higher in adipose tissue of the FL at 1, 5 and 11 wk, which produced a main effect of genotype (G).
The abundance of eight additional genes controlling metabolism and signaling of thyroid hormone and retinol was examined by qRT-PCR analysis (Additional file 8). Type I iodothyronine deiodinase (DIO1), which converts the prohormone T4 to metabolically active T3, showed only a main effect of age, whereas type III iodothyronine deiodinase (DIO3) presented main effects of age and genotype due to a consistently higher abundance in abdominal fat of LL chickens (Additional file 8-A). In contrast, the transcriptional regulator THRSPA and thioredoxin interacting protein (TXNIP) showed a main effect of age (A) and genotype (G) with higher expression in visceral fat of FL chickens at five of the six ages examined. Four genes involved in retinol metabolism [beta-carotene 15, 15′-monooxygenase (BCMO1), beta-carotene oxygenase 2 (BCO2), retinol saturase (RETSAT)] and the retinoic acid-induced gene 3 (RAIG3) [or G protein-coupled receptor, family C, group 5, member C (GPRC5C)] were also examined by qRT-PCR analysis (Additional file 8-B). Although higher in the LL between 7 and 11 wk of age, BCMO1 produced only a main effect of age (A). The expression of BCO2 sharply increased with age (main effect) and was consistently higher in abdominal fat of the LL birds (main effect of genotype). Similarly, RAIG3 showed main effects of age and genotype, with higher expression in the LL at 7 wk of age. The abundance of RETSAT was higher in visceral fat of the FL at 3 and 9 wk. Furthermore, the retinoid ligand-activated transcription factor RXRG was up-regulated in the FL, especially at 11 wk of age (Additional file 7).
An array of DE and prior candidate genes was selected for verification of gene expression using qRT-PCR analysis (Table 3; Figures 3, 4 and 6; Additional files 7 and 8). Pearson’s correlation coefficient (r) of expression ratios (FL/LL) of 15 select genes subjected to both microarray and qRT-PCR analyses indicates a significant (P≤0.01) correlation between the two methods (r = 0.64). The exclusion of two genes with the lowest microarray FC estimate (ANGPTL4 and ADIPOQ) greatly increased the Pearson correlation coefficient (r = 0.79) and the significance level (P≤0.01).
Another gene interaction network (Figure 8) identified by IPA shows interactions of several ligand-activated nuclear receptors and transcription regulators [RXRG, CEPBZ, NR1H4 (farnesoid X receptor, FXR), THRA, THRSP, MID1IP1, nuclear receptor coactivator 1 (NCOA1), forkhead box J1 (FOXJ1), and CCCTC-binding factor (CTCF)]. The target genes of these upstream regulators were up regulated in abdominal fat of the FL [GH, DNER, CYP2C9, ALAS1, CRYAB, ICMT, GPC4, SERINC1, CAMK2B, HMOX2 and SNX7] or LL chickens (CYP4F2, FABP1 and ACACA).
A final non-redundant set of genes involved in lipid metabolism was identified by IPA from the G, A and A × G DE gene lists and then was used for Ingenuity® Upstream Regulator Analysis. This analysis illustrates the interaction of numerous ligand-activated nuclear receptors and other transcription factors (TF), providing predictions of either an activated (orange color) or inhibited (blue color) state (Figure 9). These predictions are based on prior knowledge of transcriptional responses (from human and murine studies in the literature) and the responses of downstream targets found in the DE gene data set. For example, this mechanistic network of transcription regulators indicates whether the TF (orange color) and target gene (red gene symbol) are both activated or if the activity of the TF is inhibited (blue color), which would be associated with increased expression in the LL (green gene symbol). This mechanistic network predicts inhibition (blue lines and symbols) of eight transcription factors (PPARA, RXRA, NR1H2, NR1H3, PPARD, PPARG, NROB2 and NR5A2) and the activation (orange lines and symbols) of an additional eight transcription factors (NR1H4, THRB, CEBPA, CEBPB, CREB1, PPARGC1B, SREBF1 and SREBF2) (Table 4). The gene targets are presented for two transcription factors (PPARA and CEBPA) predicted to be inhibited and four transcription factors (THRB, SREBF2, CEBPB and CREB1) that were predicted to be activated by the IPA Upstream Regulator Analysis. This mechanistic analysis shows that three transcription factors (PPARA, CEBPA and CEBPB) exert direct actions on target genes up-regulated in the LL, while three other transcription factors (THRB, SREBF2 and CREB1) mainly target up-regulated genes in the FL, which are involved in the synthesis, transport or metabolism of lipids.
The divergent FL and LL chickens were originally developed as experimental models to study genetic and endocrine mechanisms controlling excessive abdominal fatness . Indeed, juvenile FL and LL chickens exhibit a 2.5-fold difference in abdominal fatness between 3 and 11 wk of age while maintaining similar rates of growth (Table 1), feed intake, and energy metabolism . The divergence of adiposity between the FL and LL chickens occurs at 3 wk of age ; hyperplasia of adipocytes was found as early as 2 wk of age in the FL , followed by marked hypertrophy of adipocytes by 9 wk of age . The FL chickens appear to favor partitioning of energy and nutrients into abdominal fat, whereas the LL deposit more protein into skeletal muscle, especially breast muscle . A consistent feature of metabolism in the FL chickens is a glucose-insulin imbalance, where plasma glucose levels are lower and insulin levels are slightly elevated [11, 27]. The higher rate of lipogenesis observed in liver of FL chickens contributes to a greater accumulation of abdominal fat in this genotype [45, 46]. The search for major genes contributing to the divergence in adiposity between the FL and LL has primarily focused on the liver [3, 19, 29–32]. In the present study, the Del-Mar 14K Integrated Systems microarray was used to examine gene expression profiles of abdominal fat in juvenile FL and LL cockerels across six ages (1–11 wk). This analysis of time-course transcriptional profiles has provided the first panoramic view of the abdominal fat transcriptome in the FL and LL chickens and given functional insight into the 2.5-fold divergence of adiposity. In particular, we have discovered numerous DE genes that are involved in hemostasis (blood coagulation), adipokine signaling, thyroid hormone and retinol action, and lipogenesis in abdominal fat of FL and LL chickens during juvenile development. These findings are unlike previous microarray studies of adipose tissue in meat-type chickens.
Higher expression of blood coagulation factors in adipose tissue of LL chickens
A large number of genes involved in hemostasis were differentially expressed in adipose tissue of LL chickens (Table 3). Several coagulation factors identified in our transcriptional analysis of adipose tissue are either proteases (i.e., F2, F9, PLG, PROC, and CFB) or protease inhibitors (i.e., A2M, ANXA5, SERPINA1, and SERPIND1). We also found higher expression of carboxypeptidases [carboxypeptidase B2 (CPB2 or thrombin-activatable fibrinolysis inhibitor) and carboxypeptidase M (CPM)] in abdominal fat of LL chickens. Our discovery of a higher abundance of several genes encoding blood clotting factors in LL chickens is quite peculiar given that fattening, rather than leanness, in mammals is usually associated with the prothrombotic state [47–50]. In fact, obesity in humans is described as chronic low-grade inflammation where expression of hemostatic genes [e.g., serine peptidase inhibitor, clade E (nexin, plasminogen activator inhibitor type 1), member 1 (PAI-1), thrombin, fibrinogen and von Willebrand factor (VWF)] are positively associated with greater deposition of adipose tissue [47, 51]. The adipokine PAI-1 (SERPINE1) encodes a secreted regulator of fibrinolysis, which serves as a biomarker for metabolic syndrome in humans . Although PAI-1 has not been mapped to the chicken genome, we did find higher expression of the plasminogen activator inhibitor 1 RNA-binding protein (SERBP1) in abdominal fat of the FL chickens. The SERBP1 protein greatly increases the degradation of PAI-1 mRNA in rat hepatoma cells . In addition, SERBP1 functions as a partner with the progesterone receptor membrane component 1 (PGRMC1) in mediating the anti-apoptotic action of progesterone on the female reproductive tract of humans . Our identification of SERBP1 and its up-regulation in adipose tissue of the FL suggest that a functional homolog of PAI-1 does exist in the chicken. Another related member of the same clade as PAI-I, SERPINE2 was not differentially expressed in abdominal fat of FL and LL chickens according to microarray analysis. Since SERPINE2 was one of the most stably expressed genes in our qRT-PCR analysis, it was used as a housekeeping gene to normalize gene expression. Another hemostatic gene up regulated in abdominal fat of LL chickens was thrombospondin 2 (THBS2), which inhibits adipogenesis in mammals .
Hemostatic proteins have several functions, some of which include removal of signal peptides, activation of zymogens, transport of enzymes, or degradation of active enzymes. Given that many adipokines have a high functional and structural similarity to the classic coagulation factors and other hemostatic factors (e.g., ANGPTL4 contains a fibrinogen-like domain), it is reasonable to assume that these proteases act on pre-pro-adipokines or other secreted proteins expressed in adipose tissue. Little is known about the expression of blood coagulation genes in visceral fat or their role in the development of adiposity in chickens. Using K-means clustering (data not shown), we found that the expression profiles for most genes involved in coagulation were clustered with those of adipokines; this general trend was verified by qRT-PCR analysis (Figures 3 and 4). For example, secreted proteins ADIPOQ and ATRN have expression patterns that are similar to those of serine proteases (e.g., F2) and protease inhibitors (e.g., ANXA5). Further, the developmental profiles of RBP4 and ANGPTL4 were similar to that of A2M, a protease inhibitor and transporter of cytokines. The positive correlation of expression patterns between coagulation factors and adipokines is not surprising given that many adipokines are associated with hemostatic or inflammatory processes (e.g. RARRES2); and conversely, several genes involved in coagulation are considered as adipokines (e.g. PAI-1, A2M, F2 and FGA). Furthermore, a similar transcriptional analysis of liver from the same individual FL and LL birds failed to reveal differential expression of these blood coagulation factors . The lack of a parallel effect of genotype on hepatic expression of coagulation factors in the FL and LL chickens suggests that their ectopic expression in abdominal fat is specific and without consequence to systemic hemostasis.
Adipokines identified in abdominal fat of FL and LL chickens
A prime example of proteolytic processing of adipokines comes from our discovery of chemerin [or RARRES2 (retinoic acid receptor responder, tazarotene induced, 2)], which is expressed at higher levels in abdominal fat of LL chickens. Chemerin is a recently discovered adipokine that regulates adipogenesis; and chemerin can be transformed into a pro-inflammatory protein, a cell adhesion factor or an anti-inflammatory peptide, depending upon cleavage by specific proteases [55–57]. After removal of the N-terminal signal peptide, pro-chemerin is processed at the C-terminal end by serine proteases to generate an active pro-inflammatory adipokine, which can be cleaved further at its C-terminal end by cysteine proteases to generate an anti-inflammatory peptide . Active chemerin appears to exert its action by binding its extracellular receptor CMKLR1 on adipocytes and/or CCRL2 on activated macrophages, which then forms an adhesive bridge between these two resident cells in adipose tissue during the inflammatory response . Adipocyte-derived chemerin causes insulin resistance in skeletal muscle cells ; and as a secreted adipokine, chemerin regulates myogenesis by providing negative cross-talk between adipose tissue and skeletal muscle . Consequently, chemerin functions as a chemokine for leukocytes, an adipokine that regulates angiogenesis, and a biomarker of metabolic syndrome and obesity in humans [61–63]. However, chemerin appears to be associated with leanness in the chicken.
Retinol binding protein 4 (RBP4), the main transporter of retinol in blood, is another adipokine that was expressed higher in abdominal fat of LL chickens at 5 and 7 wk of age. Like chemerin, RBP4 serves as a biomarker of obesity-related diseases including insulin resistance, dyslipidemia, hypertension, and visceral obesity in both adult and adolescent humans . Similar to other genes involved in coagulation, chemerin and RBP4 are expressed higher in LL, rather than FL chickens. In contrast, several adipokines (ADIPOQ, ADIPOR1 and ATRN) found in abdominal fat of FL and LL chickens are regulated similar to mammals. For example, ADIPOQ is inversely related to fatness and it is associated with increased insulin sensitivity in mammals . From the qRT-PCR analysis, we observed late up-regulation of ADIPOQ (wk 7–11) and its receptor ADIPOR1 (wk 9) in LL chickens, which suggests that adipose tissue of FL chickens could be less sensitive to insulin at these ages. Attractin (ATRN) is a neuropeptide involved in melanocortin signaling and regulation of food intake, which suppresses diet-induced obesity . Our qRT-PCR analysis shows that the expression of ATRN markedly increases in abdominal fat after 5 wk of age in both the FL and LL; furthermore, the expression pattern of ATRN is strikingly similar to that of ADIPOQ and ADIPOR1. The adipokine ANGPTL4 was expressed higher in LL chickens at 1, 5, 7 and 11 wk of age, although this difference was not statistically significant by qRT-PCR analysis. Originally, ANGPTL4 was identified as a secreted “fasting-induced adipose factor (FIAF)” in the mouse that was sharply up regulated by fasting and a target gene of the transcription factor PPARA . In fact, ANGPTL4 is a potent irreversible inhibitor of lipoprotein lipase (LPL) activity, which leads to hypertriglyceridemia . Our qRT-PCR analysis shows a 3-fold increase in expression of LPL in the LL chicken at 7 wk. Of particular interest, ANGPTL4 promotes the cleavage of LPL, while the proteolytic cleavage of ANGPTL4 by proprotein convertase releases a more potent inhibitor of LPL activity—the N terminal domain . Thus, abdominal fat of chickens is enriched with adipokines, which can exert either local (autocrine/paracrine) or systemic (endocrine) actions after proteolytic processing and secretion into circulation (Table 3).
Our initial survey of global gene expression in abdominal fat of juvenile FL and LL chickens highlights another important feature of the avian endocrine system—the virtual absence of several important adipokines normally found in mammals. A few examples of adipokines not yet mapped to the current draft of the chicken genome (galGAL4), include leptin (LEP), omentin (ITLN1), resistin (RETN), tumor necrosis factor alpha (TNFA), and PAI-1. The existence of the chicken LEP gene remains a great and unresolved controversy [70–74], especially since extensive expressed sequence tag (EST)  and whole genome sequencing projects have failed to identify a bona fide LEP gene in the chicken. Furthermore, the LEP gene is absent from the genome of all birds sequenced so far (i.e., chicken, turkey, zebra finch, budgerigar and duck). However, the leptin receptor (LEPR) gene is expressed in several chicken tissues [75–79]; and chicken LEPR is capable of activating the JAK-STAT pathway in vitro[80, 81]. Similarly, components of TNF signaling are up regulated in the hypothalamus of LL chickens , although TNFA is yet to be identified in chickens. Despite the absence of several mammalian adipokines (i.e., LEP, TNFA, RETN, PAI-1, APOE, and ITLN1) and metabolic enzymes (i.e., LIPE), adipogenesis and lipid metabolism in the chicken are robustly regulated by mechanisms that are, for the most part, similar to those described in mammals.
Retinol metabolism and retinoic acid signaling in adipose tissue
Another remarkable observation from the present study was the over expression of 13 genes in abdominal fat of LL chickens that control metabolism of retinol, the precursor of retinoic acid (RA), which itself is a major chemical activator of multiple transcription factors controlling lipogenesis. The primary source of retinol is dietary plant-based β-carotene, which is symmetrically cleaved by the enzyme β-carotene monooxygenase 1 (BCMO1) into two molecules of retinal. Recently, we discovered mutations in the proximal promoter of BCMO1, which are responsible for variation in the color of breast meat in another F2 resource population of meat-type chickens . Another enzyme, β-carotene oxygenase 2 (BCO2), asymmetrically cleaves one molecule of β-carotene to generate one molecule of retinal and a by-product (e.g., β-apo-14′-carotenal), which acts downstream to block signaling of PPARG . The BCO2 gene in chickens was originally identified as the yellow skin gene, which controls the β-carotene content and thereby yellow pigmentation of the skin . Our qRT-PCR analysis of these two β-carotene degrading enzymes (Figure 8), showed only a main effect of age on expression of BCMO1, whereas the abundance of BCO2 was greater in abdominal fat of LL chickens, producing a main effect of genotype. Another study found increased expression of BCO2 in adipocytes from BCMO1 knockout mice and that dietary β-carotene reduces adiposity of mice—but only in the presence of a functional BCMO1 enzyme . This study also demonstrates the importance of BCMO1 in generating the precursor (retinal) for RA, which inhibits activation of PPARG and its lipogenic target genes that are mainly metabolic enzymes, adipokines and transport proteins (see Additional file 6). Our study shows higher expression of both BCMO1 and BCO2 in abdominal fat of the LL chickens after 5 wk of age, which presumably would lead to generation of more retinal and RA. This idea is supported by the differential expression of several genes involved in retinol metabolism and RA signaling in adipose tissue of FL and LL chickens. These genes are involved in transport of retinol (RBP4, TTR and RBP7), metabolism of retinol (RDH1, RETSAT, ADH1C, ADH5, and CYP2E1), and respond to RA (RARRES2, GPRC5C, and NOL7). In 3T3-L1 preadipocytes, RA inhibits adipogenesis via up-regulation of the transcriptional modulator SMAD3 . Interestingly, two members of the SMAD family (SMAD5 and SMAD6) were up-regulated (main effect of age) in adipose tissue of LL chickens (Additional file 3). The ligand (RA) activates its nuclear receptors (RAR and RXR), which can form heterodimers with other ligand-dependent transcription factors (e.g., LXR, PPARG and THR) to initiate transcription of numerous downstream target genes. Thus, RA seems to play an important role in reduction of adipogenesis and adiposity in the LL chickens.
Visceral adipose tissue as a major site of lipogenesis in chickens
Because the liver is widely considered as the primary site of de novo lipid synthesis in birds, most transcriptional studies of lipogenesis in the chicken have focused on liver rather than adipose tissue. A targeted low-density array enabled an initial transcriptional analysis of liver [at a single age (8 wk)] in the FL and LL chickens . This study showed up-regulation of several lipogenic genes (ACACA, FASN, SCD, APOA1, SREBF1, and MDH2) in the FL chickens. Examination of hepatic gene expression at three ages (1, 4, and 7 wk) in another population of chickens divergently selected on abdominal fatness revealed differential expression of several genes involved in lipid metabolism, including ACAT1, CEBPγ, FABP1, APOA1, MDH1, APOD and PPARG. A time-course (1–11 wk) transcriptional study of liver in juvenile FL and LL chickens revealed 1,805 differentially expressed (DE) genes, mostly between 7 and 11 wk . These functional genes identified in the liver of juvenile FL and LL birds were transcription factors, metabolic enzymes, transport proteins, differentiation factors, signaling molecules and adipokines.
In contrast, there have been only a few transcriptional studies of adipose tissue in the chicken. For example, a comparison of abdominal fat between meat-type (broiler) and egg-type chickens (layer) at a single physical age (10 wk), albeit at different physiological ages, focused attention on the up-regulation of LPL in broiler chickens and higher expression of APOA1 in layers . Another study using abdominal fat samples taken at 7 wk from a different population of divergently selected fat and lean chickens reported the differentially expression of 230 adipose genes [153 were up-regulated in the fat chickens, while 77 were up-regulated in the lean birds] . Their conclusion that TNFA plays a key role in lipid metabolism of the chicken is surprising, since this adipokine has not been mapped to the chicken genome sequence. A recent transcriptional study of chicken abdominal fat compared a short-term (5 hr) fasting response with acute insulin immunoneutralization . Surprisingly, fasting provoked larger changes in adipose gene transcription (1,780 DE genes) than did insulin neutralization with only 92 DE genes, which confirms the insensitivity of chicken abdominal fat to insulin . In contrast, more than a thousand genes were either differentially expressed in liver or leg muscle of the same birds following insulin immunoneutralization . Nonetheless, short-term fasting in the chicken depressed the expression of 40 genes in abdominal fat that are involved in the synthesis and storage of lipid, while a number of adipose genes that control lipolysis and oxidation of fatty acids were up-regulated by fasting or insulin neutralization .
The present study has identified a large number of lipogenic genes that are up regulated in abdominal fat of FL chickens (Table 3). A prime example of this lipogenic group is our clone for GH1 (GenBank accession BI390457) that corresponds to the short form of chicken GH (scGH), which lacks a signal peptide and is highly expressed in ocular tissue [90, 91], pituitary gland and heart of chick embryos . The short alternatively-spliced (16.5 kDa) isoform of full length GH (20 kDa) appears to function as an “intracrine” factor within the cell . Our discovery of higher expression of scGH in abdominal fat of the FL chicken supports the idea of a local lipogenic action of GH on adipose tissue, rather than the lipolytic response usually observed in mammals . In fact, our earlier work clearly established the lipogenic action of exogenous GH in juvenile chickens [93–96].
Up-regulation of two transcription factors (SREBF1 and THRSPA) and the histone deacetylase SIRT2 in abdominal fat of the FL were accompanied by higher expression of multiple genes involved in the generation and metabolism of lipids (DHCR7, FADS2, FASN, HMGCR, HMGCS2, LSS, MVD, SCD and SC5DL). The higher expression of the transcription factor SREBF1 and 12 lipogenic target genes in the FL strongly supports our idea that the divergence in abdominal fatness of FL and LL chickens could be related to differential expression of several lipogenic genes in abdominal fat of the FL. For example, FADS2, which catalyzes the rate limiting step in synthesis of highly unsaturated fatty acids, was highly up regulated in abdominal fat of FL chickens; binding sites for both SREBF1 and PPARA are found in the promoter region of FADS2. Likewise, SREBF1 regulates transcription of several genes that control synthesis of fatty acids, including acetyl-CoA carboxylase alpha (ACACA), which catalyzes the rate-limiting step in fatty acid synthesis [98, 99]. Furthermore, MID1 interacting protein 1 (MID1IP1) [or THRSP-like (THRSPL)] enhances ACACA polymerization and its enzymatic activity . Adipose tissue from FL chickens shows higher expression of THRSPA, a transcriptional regulator of several lipogenic genes in the chicken [3, 78, 96, 101]. Earlier, we discovered a 9-base pair deletion near the putative DNA-binding domain of chicken THRSPA and demonstrated association of this insertion/deletion polymorphism with abdominal fatness traits in multiple resource populations of chickens . Mutations in the THRSPA gene of chickens [101–103], ducks  and geese  are associated with fatness traits and are perhaps of potential use as molecular markers in poultry breeding programs. Furthermore, THRSP is a major regulator of adipogenesis in skeletal muscle of beef cattle [106, 107] and of lipogenesis in the lactating mammary gland of the dairy cow [108–110]. Interestingly, the THRSP-null mouse shows reduced lipogenesis in the mammary gland  and pups from the THRSP-null mouse exhibit reduced body weight gain due to diminished milk triglycerides . In humans, amplification of the THRSP locus is associated with lipogenic breast cancer ; and, as such, THRSP serves as a marker of aggressive breast cancer and a potential target of anti-cancer drugs . In humans, expression of THRSP in adipose tissue is depressed by transition from a lipogenic fed state to a lipolytic state induced by a 48 hr fast . These observations support the idea that THRSP is a transcriptional activator of several lipogenic enzymes (ACLY, FASN and ME) in the mouse . THRSP is activated in response to T3, glucose and insulin and inhibited by polyunsaturated fatty acids , cyclic AMP or glucagon . Recent work has shown that induction of THRSP increases expression of FASN in cultured hepatocyte cells and RNAi-mediated knock-down of THRSP depresses expression of FASN . Another study showed that FASN co-precipitates with THRSP in nuclear extracts from the mouse (referenced in ). The exact mechanism by which THRSP and MID1IP1 interact and work as regulators of gene transcription is currently unknown. These genes are highly expressed in fatty tissues of birds and mammals, where they regulate the expression and activity of multiple lipogenic enzymes. The proximal (4 kb) promoter region of THRSPA contains four putative binding sites for PPARG and four SREBF sites (L.A. Cogburn, unpublished observations). In the present study, we found higher expression of THRSPA in abdominal fat of FL chickens at all ages, except at 7 wk. In the rat, the far-upstream region of the THRSP promoter contains three T3-THR response elements (TREs) . Thus, THRSPA is responsive to metabolically active thyroid hormone (T3) generated by the activation enzyme DIO1, whereas the enzyme DIO3 is responsible for degradation of metabolically active T3 and conversion of the prohormone (T4) to metabolically inactive reverse T3 (rT3) . The up regulation of DIO3 in adipose tissue of juvenile LL chickens (1–11 wk) suggests that less T3 would be available to activate THRSPA transcription, which was observed in the LL. Thioredoxin interacting protein (TXNIP) is another important regulator of hepatic glucose metabolism  that also mediates hypothalamic control over energy utilization and adiposity in the mouse . The up-regulation of TXNIP in abdominal fat of the FL during the period of maximal fatness (3–11 wk) could contribute to their enhanced lipogenesis and adiposity. Likewise, we have discovered another putative sensor of glucose, the sweet taste receptor 1 (TAS1R1) gene, which is differentially expressed in the hypothalamus  and abdominal fat (Figure 2) of FL and LL chickens. Our observation of higher expression of TAS1R1 in the hypothalamus of the FL and abdominal fat of the LL suggest tissue specific regulation of this important tissue glucose sensor [124–126].
Increased lipolysis in abdominal fat of LL chickens
In contrast to the enhanced lipogenic state found in abdominal fat of FL chickens, the LL show higher expression of numerous genes involved in lipolysis (Table 3). Two cytochrome P450 family members (CYP27A1 and CYP2E1) were expressed at higher levels in abdominal fat of the LL when compared to FL chickens. CYP27A1 is involved in clearance of cholesterol via bile excretion, whereas CYP2E1 is strongly induced in white adipose tissue of the rat by prolonged fasting . The beta-subunit (HADHB) of mitochondrial tri-functional protein (MTP), a complex that catalyzes the final three steps of β-oxidation of long chain fatty acids, was also up regulated in adipose tissue of LL chickens. MTP knockout mice exhibit neonatal hypoglycemia and sudden neonatal death indicating its essential role in β-oxidation of long chain fatty acids . Two members of the hydroxysteroid (17-β) dehydrogenase family (17β-HSD, members 4 and 6) were also expressed higher in adipose tissue of the LL. The significance of HSD17B4 in β-oxidation of branched chain fatty acids was demonstrated in HSD17B4 knockout mice, which were unable to degrade phytanic and pristanic acids . Since the other 17β-HSD (HSD17B6) exhibits retinol dehydrogenase activity , its up-regulation in adipose tissue of the LL chicken suggests increased availability of all trans-retinoic acid. In addition, PDK4, which inhibits the pyruvate dehydrogenase complex and conversion of pyruvate to acetyl-CoA , was highly up-regulated in abdominal fat of the LL between 7 and 11 wk. The high expression of PDK4 in the liver of chicken embryos , whose energy is derived exclusively from yolk lipids, supports a similar action of PDK4 in adipose tissue of the LL chickens. Furthermore, the expression of PDK4 increased sharply in abdominal fat of two-week-old chickens by a 5-h fast or insulin immunoneutralization . The tumor suppressor protein p53 enhances lipid catabolism and induces expression of guanidinoacetate N-methyltransferase (GAMT), which enhances β-oxidation of fatty acids ; both of these genes were up regulated in visceral fat of the LL chickens. Another gene up regulated in the LL that protects against oxidative stress is SOD3, which is expressed at higher levels in the liver of low-growth (leaner) chickens . Thus, the present time-course transcriptional analysis of abdominal fat in juvenile FL and LL chickens provides compelling evidence for enhanced lipolysis in adipose tissue of the LL.
The present study adds a new dimension to function of visceral fat as a proteolytic processor of adipokines and other endocrine signals that control lipid metabolism. In particular, the LL cockerels exhibit high expression of several blood coagulation factors in adipose tissue, albeit not in liver. Some of these changes in the LL occur before the divergence in fatness. These hemostatic proteases and protease inhibitors could be involved in activation of adipokines, chemokines and other metabolic ligands that contribute to suppression of lipogenesis and adipogenesis in the LL. Furthermore, abdominal fat of the LL chickens has higher expression of genes involved in mobilization, utilization and export of lipids than does the FL. Several transcription factors have a larger number of target genes expressed higher in the LL that could also favor suppression of abdominal fat accretion. In contrast, abdominal fat of the FL chickens expresses a greater abundance of numerous target genes involved in lipogenesis and adipogenesis, which could contribute to their greater adiposity. The higher expression of these target genes in FL chickens appears after the onset of divergence in fatness. Therefore, abdominal (visceral) fat of the chicken could play a more significant role in lipogenesis and adiposity than previously considered. The assumption that the liver of birds serves as the major site of lipogenesis needs to be re-examined.
Availability of supporting data
The minimum information about microarray experiments (MIAME)-compliant microarray data described in this article are available in the NCBI Gene Expression Omnibus (GEO) under the accession number GSE37585. Additional file 3 provides annotated DE gene sets from statistical analysis of the microarray data; whereas Additional files 4, 5, 6 provide annotation, functional analysis and mapping of DE genes to biological functions and canonical pathways by IPA software.
Fat line chickens
Lean line chickens
Main effect of genotype
- A × G:
Age by genotype interaction
Main effect of age
Week of age
Expression quantitative trait loci
Gallus gallus chromosome 5
Ingenuity Pathway Analysis
Institut National de la Recherche Agronomique
Quantitative reverse transcriptase polymerase chain reaction
Linear models for statistical analysis of microarray data
False discovery rate
Statistical Analysis System
General linear models
Least Squares Means.
International Chicken Genome Sequencing Consortium: Sequence and comparative analysis of the chicken genome provide unique perspectives on vertebrate evolution. Nature. 2004, 432: 695-716. 10.1038/nature03154.
Burt DW: Emergence of the chicken as a model organism: implications for agriculture and biology. Poult Sci. 2007, 86: 1460-1471.
Cogburn LA, Porter TE, Duclos MJ, Simon J, Burgess SC, Zhu JJ, Cheng HH, Dodgson JB, Burnside J: Functional genomics of the chicken–a model organism. Poult Sci. 2007, 86: 2059-2094.
Dodgson JB: The chicken genome: some good news and some bad news. Poult Sci. 2007, 86: 1453-1459.
Stern CD: The chick: a great model system becomes even greater. Dev Cell. 2005, 8: 9-17.
Model Organisms for Biomedical Research.http://www.nih.gov/science/models,
Belo PS, Romsos DR, Leveille GA: Blood metabolites and glucose metabolism in the fed and fasted chicken. J Nutr. 1976, 106: 1135-1143.
Simon J: Chicken as a useful species for the comprehension of insulin action. Crit Rev Poult Biol. 1989, 2: 121-148.
Simon J: Insulin in birds: metabolic effects and possible implications in genetically fat and lean chickens. Leanness in domestic birds.--Genetic, metabolic and hormonal aspects. Edited by: Leclercq B, Whitehead CC. 1987, London: Butterworths, 253-268.
Simon J, Guillaumin S, Chevalier B, Derouet M, Guy G, Marche G, Ricard FH, Leclercq B: Plasma glucose-insulin relationship in chicken lines selected for high or low fasting glycaemia. Br Poult Sci. 2000, 41: 424-429. 10.1080/713654969.
Touchburn S, Simon J, Leclercq B: Evidence of a glucose-insulin imbalance and effect of dietary protein and energy level in chickens selected for high abdominal fat content. J Nutr. 1981, 325: 335-
Goodridge AG, Ball EG: Lipogenesis in the pigeon: in vivo studies. Amer J Physiol. 1967, 213: 245-249.
Rankinen T, Zuberi A, Chagnon YC, Weisnagel SJ, Argyropoulos G, Walts B, Perusse L, Bouchard C: The human obesity gene map: the 2005 update. Obesity (Silver Spring). 2006, 14: 529-644. 10.1038/oby.2006.71.
Ikeobi CON, Woolliams JA, Morrice DR, Windsor D, Burt DW, Hocking PM: Quantitative trait loci affecting fatness in the chicken. Anim Genet. 2002, 33: 428-435. 10.1046/j.1365-2052.2002.00911.x.
Jennen DGJ, Vereijken ALJ, Bovenhuis H, Crooijmans RPMA, Veenendaal A, van der Poel JJ, Groenen MAM: Detection and localization of quantitative trait loci affecting fatness in broilers. Poult Sci. 2004, 83: 295-301.
Abasht B, Pitel F, Lagarrigue S, Le Bihan-Duval E, Pascale LR, Demeure O, Vignoles F, Simon J, Cogburn L, Aggrey S, Vignal A, Douaire M: Fatness QTL on chicken chromosome 5 and interaction with sex. Genet Sel Evol. 2006, 38: 297-311. 10.1186/1297-9686-38-3-297.
Lagarrigue S, Pitel F, Carré W, Abasht B, Le Roy P, Neau A, Amigues Y, Sourdioux M, Simon J, Cogburn LA, Aggrey S, Leclercq B, Vignal A, Douaire M: Mapping quantitative trait loci affecting fatness and breast muscle weight in meat-type chicken lines divergently selected on abdominal fatness. Genet Sel Evol. 2006, 38: 85-97. 10.1186/1297-9686-38-1-85.
Mignon GL, Pitel F, Gilbert H, Bihan-Duval EL, Vignoles F, Demeure O, Lagarrigue S, Simon J, Cogburn LA, Aggrey SE, Douaire M, Roy PL: A comprehensive analysis of QTL for abdominal fat and breast muscle weights on chicken chromosome 5 using a multivariate approach. Anim Genet. 2008, 40: 157-164.
Le Mignon G, Desert C, Pitel F, Leroux S, Demeure OI, Guernec G, Abasht B, Douaire M, Le Roy P, Lagarrigue S: Using transcriptome profiling to characterize QTL regions on chicken chromosome 5. BMC Genomics. 2009, 10: 575-10.1186/1471-2164-10-575.
Simon J, Milenkovic D, Godet E, Cabau C, Collin A, Metayer-Coustard S, Rideau N, Tesseraud S, Derouet M, Crochet S, Cailleau-Audouin E, Hennequet-Antier C, Gespach C, Porter TE, Duclos MJ, Dupont J, Cogburn LA: Insulin immuno-neutralization in fed chickens: effects on liver and muscle transcriptome. Physiol Genomics. 2012, 44: 283-292. 10.1152/physiolgenomics.00057.2011.
Ji B, Ernest B, Gooding J, Das S, Saxton A, Simon J, Dupont J, Metayer-Coustard S, Campagna S, Voy B: Transcriptomic and metabolomic profiling of chicken adipose tissue in response to insulin neutralization and fasting. BMC Genomics. 2012, 13: 441-10.1186/1471-2164-13-441.
Dupont J, Métayer-Coustard S, Ji B, Ramé C, Gespach C, Voy B, Simon J: Characterization of major elements of insulin signaling cascade in chicken adipose tissue: apparent insulin refractoriness. Gen Comp Endocrinol. 2012, 176: 86-93. 10.1016/j.ygcen.2011.12.030.
Nilsson C, Raun K, Yan F, Larsen MO, Tang-Christensen M: Laboratory animals as surrogate models of human obesity. Acta Pharmacol Sin. 2012, 33: 173-181. 10.1038/aps.2011.203.
Leclercq B, Blum JC, Boyer JP: Selecting broilers for low or high abdominal fat: initial observations. Brit Poult Sci. 1980, 21: 107-113. 10.1080/00071668008416644.
Baéza E, Le Bihan-Duval E: Chicken lines divergent for low or high abdominal fat deposition: a relevant model to study the regulation of energy metabolism. Animal. 2013, 7: 965-973. 10.1017/S1751731113000153.
Leclercq B: Genetic selection of meat-type chickens for high or low abdominal fat content. Leanness in Domestic Birds. Edited by: Leclercq B, Whitehead CC. 1988, London: Butterworths, 25-40.
Simon J, Leclercq B: Longitudinal study of adiposity in chickens selected for high or low abdominal fat content; further evidence of a glucose-insulin imbalance in the fat line. J Nutr. 1982, 112: 1961-1973.
Hermier D, Quignard-Boulange A, Dugail I, Guy G, Salichon MR, Brigant L, Ardouin B, Leclercq B: Evidence of enhanced storage capacity in adipose tissue of genetically fat chickens. J Nutr. 1989, 119: 1369-1375.
Carré W, Bourneuf E, Douaire M, Diot C: Differential expression and genetic variation of hepatic messenger RNAs from genetically lean and fat chickens. Gene. 2002, 299: 235-243. 10.1016/S0378-1119(02)01077-6.
Daval S, Lagarrigue S, Douaire M: Messenger RNA levels and transcription rates of hepatic lipogenesis genes in genetically lean and fat chickens. Genet Sel Evol. 2000, 32: 521-531. 10.1186/1297-9686-32-5-521.
Assaf S, Lagarrigue S, Daval S, Sansom M, Leclercq B, Michel J, Pitel F, Alizadeh M, Vignal A, Douaire M: Genetic linkage and expression analysis of SREBP and lipogenic genes in fat and lean chicken. Comp Biochem Physiol B, Biochem Mol Biol. 2004, 137: 433-441. 10.1016/j.cbpc.2004.02.005.
Bourneuf E, Herault F, Chicault C, Carre W, Assaf S, Monnier A, Mottier S, Lagarrigue S, Douaire M, Mosser J, Diot C: Microarray analysis of differential gene expression in the liver of lean and fat chickens. Gene. 2006, 372: 162-170.
Le Mignon G, Pitel F, Gilbert H, LeBihan-Duval E, Vignoles F, Demeure O, Lagarrigue S, Simon J, Cogburn LA, Aggrey SE, Douarin NM, Le Roy P: A comprehensive analysis of QTL for abdominal fat and breast muscle weights on chicken chromosome 5 using a multivariate approach. Anim Genet. 2009, 40: 157-164. 10.1111/j.1365-2052.2008.01817.x.
Chirgwin JM, Przybla AE, MacDonald RJ, Rutter WJ: Isolation of biologically active ribonucleic acid from sources enriched in ribonuclease. Biochemistry. 1979, 18: 5294-5299. 10.1021/bi00591a005.
Smyth GK: Limma: Linear models for microarray data. 2005, New York: Springer, 397-420.
Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-
Benjamini Y, Hochberg Y: Controlling the false discovery rate: a practical and powerful approach to multiple testing. J Royal Stat Soc. 1995, 57B: 289-300.
Cotton EST Database.http://www.leonxie.com/referencegene.php,
Vandesompele J, Preter KD, Pattyn F, Poppe B, Roy NV, Paepe AD, Speleman F: Accurate normalization of real-time quantitative RT-PCR data by geometric averaging of multiple internal control genes. Genome Biol. 2002, 3: research0034-
Leclercq B, Hermier D, Guy G: Metabolism of very low density lipoproteins in genetically lean or fat lines of chicken. Reprod Nutr Dev. 1990, 30: 701-715. 10.1051/rnd:19900607.
Geraert PA, Macleod MG, Leclercq B: Energy metabolism in genetically fat and lean chickens: diet-and cold-induced thermogenesis. J Nutr. 1988, 118: 1232-1239.
Saadoun A, Leclercq B: In vivo lipogenesis in genetically fat and lean chickens of various ages. Comp Biochem Physiol. 1986, 83B: 607-611.
Saadoun A, Leclercq B: In vivo lipogenesis of genetically lean and fat chickens: effects of nutritional state and dietary fat. J Nutr. 1987, 117: 428-435.
Faber DR, De Groot PG, Visseren FLJ: Role of adipose tissue in haemostasis, coagulation and fibrinolysis. Obes Rev. 2009, 10: 554-563. 10.1111/j.1467-789X.2009.00593.x.
Allman-Farinelli MA: Obesity and venous thrombosis: a review. Semin Thromb Hemost. 2011, 37: 903-907. 10.1055/s-0031-1297369.
Schäfer K, Konstantinides S: Adipokines and thrombosis. Clin Exper Pharmacol and Physiol. 2011, 38: 864-871. 10.1111/j.1440-1681.2011.05589.x.
Lorenzet R, Napoleone E, Cutrone A, Donati MB: Thrombosis and obesity: cellular bases. Thromb Res. 2012, 129: 285-289. 10.1016/j.thromres.2011.10.021.
Palomo I, Alarcon M, Moore-Carrasco R, Argiles JM: Hemostasis alterations in metabolic syndrome (review). Int J Mol Med. 2006, 18: 969-974.
Heaton JH, Dlakic WM, Dlakic M, Gelehrter TD: Identification and cDNA cloning of a novel RNA-binding protein that interacts with the cyclic nucleotide-responsive sequence in the type-1 plasminogen activator inhibitor mRNA. J Biol Chem. 2001, 276: 3341-3347. 10.1074/jbc.M006538200.
Zhang L, Kanda Y, Roberts DJ, Ecker JL, Losel R, Wehling M, Peluso JJ, Pru JK: Expression of progesterone receptor membrane component 1 and its partner serpine 1 mRNA binding protein in uterine and placental tissues of the mouse and human. Mol Cell Endocrinol. 2008, 287: 81-89. 10.1016/j.mce.2008.02.012.
Shitaye HS, Terkhorn SP, Combs JA, Hankenson KD: Thrombospondin-2 is an endogenous adipocyte inhibitor. Matrix Biol. 2010, 29: 549-556. 10.1016/j.matbio.2010.05.006.
Ernst MC, Sinal CJ: Chemerin: at the crossroads of inflammation and obesity. Trends Endocrinol Metab. 2010, 21: 660-667. 10.1016/j.tem.2010.08.001.
Bondue B, Wittamer V, Parmentier M: Chemerin and its receptors in leukocyte trafficking, inflammation and metabolism. Cytokine Growth Factor Rev. 2011, 22: 331-338. 10.1016/j.cytogfr.2011.11.004.
Roman AA, Parlee SD, Sinal CJ: Chemerin: a potential endocrine link between obesity and type 2 diabetes. Endocrine. 2012, 42: 243-251. 10.1007/s12020-012-9698-8.
Yoshimura T, Oppenheim JJ: Chemerin reveals its chimeric nature. J Exp Med. 2008, 205: 2187-2190. 10.1084/jem.20081736.
Sell H, Laurencikiene J, Taube A, Eckardt K, Cramer A, Horrighs A, Arner P, Eckel J: Chemerin is a novel adipocyte-derived factor inducing insulin resistance in primary human skeletal muscle cells. Diabetes. 2009, 58: 2731-2740. 10.2337/db09-0277.
Yang H, Li F, Kong X, Yuan X, Wang W, Huang R, Li T, Geng M, Wu G, Yin Y: Chemerin regulates proliferation and differentiation of myoblast cells via ERK1/2 and mTOR signaling pathways. Cytokine. 2012, 60: 646-652. 10.1016/j.cyto.2012.07.033.
Goralski KB, McCarthy TC, Hanniman EA, Zabel BA, Butcher EC, Parlee SD, Muruganandan S, Sinal CJ: Chemerin, a novel adipokine that regulates adipogenesis and adipocyte metabolism. J Biol Chem. 2007, 282: 28175-28188. 10.1074/jbc.M700793200.
Bozaoglu K, Segal D, Shields KA, Cummings N, Curran JE, Comuzzie AG, Mahaney MC, Rainwater DL, VandeBerg JL, MacCluer JW, Collier G, Blangero J, Walder K, Jowett JBM: Chemerin Is associated with metabolic syndrome phenotypes in a Mexican-American population. J Clin Endocr Metab. 2009, 94: 3085-3088. 10.1210/jc.2008-1833.
Bozaoglu K, Curran JE, Stocker CJ, Zaibi MS, Segal D, Konstantopoulos N, Morrison S, Carless M, Dyer TD, Cole SA, Goring HHH, Moses EK, Walder K, Cawthorne MA, Blangero J, Jowett JBM: Chemerin, a novel adipokine in the regulation of angiogenesis. J Clin Endocrinol Metab. 2010, 95: 2476-2485. 10.1210/jc.2010-0042.
Conroy R, Espinal Y, Fennoy I, Accacha S, Boucher-Berry C, Carey D, Close S, DeSantis D, Gupta R, Hassoun A, Iazzetti L, Jacques F, Jean A, Michel L, Pavlovich K, Rapaport R, Rosenfeld W, Shamoon E, Shelov S, Speiser P, Ten S, Rosenbaum M: Retinol binding protein 4 is associated with adiposity-related co-morbidity risk factors in children. J Pediatr Endocrinol Metab. 2011, 24: 913-919.
Shehzad A, Igbal W, Shehzad O, Lee YS: Adiponectin: regulation of its production and its role in human diseases. Hormones (Athens). 2012, 11: 8-20.
Matarese G, La Cava A: The intricate interface between immune system and metabolism. Trends Immunol. 2004, 25: 193-200. 10.1016/j.it.2004.02.009.
Kersten S, Mandard S, Tan NS, Escher P, Metzger D, Chambon P, Gonzalez FJ, Desvergne B, Wahli W: Characterization of the fasting-induced adipose factor FIAF, a novel peroxisome proliferator-activated receptor target gene. J Biol Chem. 2000, 275: 28488-28493.
Shan L, Yu XC, Liu Z, Hu Y, Sturgis LT, Miranda ML, Liu Q: The angiopoietin-like proteins ANGPTL3 and ANGPTL4 inhibit lipoprotein lipase activity through distinct mechanisms. J Biol Chem. 2009, 284: 1419-1424.
Lei X, Shi F, Basu D, Huq A, Routhier S, Day R, Jin W: Proteolytic processing of angiopoietin-like protein 4 by proprotein convertases modulates its inhibitory effects on lipoprotein lipase activity. J Biol Chem. 2011, 286: 15747-15756. 10.1074/jbc.M110.217638.
Friedman-Einat M, Boswell T, Horev G, Girishvarma G, Dunn IC, Talbot RT, Sharp PJ: The chicken leptin gene: Has it been cloned?. Gen Comp Endocrinol. 1999, 115: 354-363. 10.1006/gcen.1999.7322.
Sharp PJ, Dunn IC, Waddington D: Chicken leptin. Gen Comp Endocrinol. 2008, 158: 2-4. 10.1016/j.ygcen.2008.05.018.
Simon J, Rideau N, Taouis M: Reply to viewpoints by PJ Sharp, IC Dunn, D Waddington and T Boswell [Chicken Leptin: General and Comparative Endocrinology, 158, 2–4 (2008)]. Gen Comp Endocrinol. 2009, 161: 159-
Pitel F, Faraut T, Bruneau G, Monget P: Is there a leptin gene in the chicken genome? Lessons from phylogenetics, bioinformatics and genomics. Gen Comp Endocrinol. 2010, 167: 1-5. 10.1016/j.ygcen.2009.10.006.
Yosefi S, Hen G, Rosenblum CI, Cerasale DJ, Beaulieu M, Criscuolo F, Friedman-Einat M: Lack of leptin activity in blood samples of Adélie penguin and bar-tailed godwit. J Endocrinol. 2010, 207: 113-122. 10.1677/JOE-10-0177.
Carré W, Wang X, Porter TE, Nys Y, Tang J-S, Bernberg E, Morgan R, Burnside J, Aggrey SE, Simon J, Cogburn LA: Chicken genomics resource: sequencing and annotation of 35,407 chicken ESTs from single and multiple tissue cDNA libraries and CAP3 assembly of a chicken gene index. Physiol Genomics. 2006, 25: 514-524. 10.1152/physiolgenomics.00207.2005.
Horev G, Einat P, Aharoni T, Eshdat Y, Friedman-Einat M: Molecular cloning and properties of the chicken leptin-receptor (CLEPR) gene. Mol Cell Enocrinol. 2000, 162: 95-106. 10.1016/S0303-7207(00)00205-7.
Ohkubo T, Tanaka M, Nakashima K: Structure and tissue distribution of chicken leptin receptor (cOb-R) mRNA. Biochem Biophys Acta. 2000, 1491: 303-308. 10.1016/S0167-4781(00)00046-4.
Cogburn LA, Wang X, Carré W, Rejto L, Porter TE, Aggrey SE, Simon J: Systems-wide chicken DNA microarrays, gene expression profiling and discovery of functional genes. Poult Sci. 2003, 82: 939-951.
Byerly MS, Simon J, Cogburn LA, Le Bihan-Duval E, Duclos MJ, Aggrey SE, Porter TE: Transcriptional profiling of hypothalamus during development of adiposity in genetically selected fat and lean chickens. Physiol Genomics. 2010, 42: 157-167. 10.1152/physiolgenomics.00029.2010.
Adachi H, Takemoto Y, Bungo T, Ohkubo T: Chicken leptin receptor is functional in activating JAK-STAT pathway in vitro. J Endocrnol. 2008, 197: 335-342. 10.1677/JOE-08-0098.
Hen G, Yosefi S, Ronin A, Einat P, Rosenblum CI, Denver RJ, Friedman-Einat M: Monitoring leptin activity using the chicken leptin receptor. J Endocrinol. 2008, 197: 325-333. 10.1677/JOE-08-0065.
Le Bihan-Duval E, Nadaf J, Berri C, Pitel F, Graulet B, Godet E, Leroux SY, Demeure O, Lagarrigue S, Duby C, Cogburn LA, Beaumont CM, Duclos MJ: Detection of a cis eQTL controlling BCMO1 gene expression leads to the identification of a QTG for chicken breast meat color. Plos One. 2011, 6: e14825-10.1371/journal.pone.0014825.
Tourniaire F, Gouranton E, von Lintig J, Keijer J, Luisa BM, Amengual J, Lietz G, Landrier JF: b-Carotene conversion products and their effects on adipose tissue. Genes Nutr. 2009, 4: 179-187. 10.1007/s12263-009-0128-3.
Eriksson J, Larson G, Gunnarsson U, Bed’hom B, Tixier-Boichard M, Strömstedt L, Wright D, Jungerius A, Vereijken A, Randi E, Jensen P, Andersson L: Identification of the yellow skin gene reveals a hybrid origin of the domestic chicken. PLoS Genet. 2008, 4: e1000010-10.1371/journal.pgen.1000010.
Amengual J, Gouranton E, van Helden YGJ, Hessel S, Ribot J, Kramer E, Kiec-Wilk B, Razny U, Lietz G, Wyss A, Dembinska-Kiec A, Palou A, Keijer J, Landrier JF, Bonet ML, von Lintig J: Beta-carotene reduces body adiposity of mice via BCMO1. Plos One. 2011, 6: e20644-10.1371/journal.pone.0020644.
Marchildon F, St-Louis C, Akter R, Roodman V, Wiper-Bergeron NL: Transcription factor Smad3 is required for the inhibition of adipogenesis by retinoic acid. J Biol Chem. 2010, 285: 13274-13284. 10.1074/jbc.M109.054536.
Wang HB, Wang QG, Zhang XY, Gu XF, Wang N, Wu SB, Li H: Microarray analysis of genes differentially expressed in the liver of lean and fat chickens. Animal. 2010, 4: 513-522. 10.1017/S1751731109991388.
Wang H, Li H, Wang Q, Wang Y, Han H, Shi H: Microarray analysis of adipose tissue gene expression profiles between two chicken breeds. J Biosci. 2006, 31: 565-573. 10.1007/BF02708408.
Wang HB, Li H, Wang QG, Zhang XY, Wang SZ, Wang YX, Wang XP: Profiling of chicken adipose tissue gene expression by genome array. BMC Genomics. 2007, 8: 193-10.1186/1471-2164-8-193.
Takeuchi S, Haneda M, Teshigawara K, Takahashi S: Identification of a novel GH isoform: A possible link between GH and melanocortin systems in the developing chicken eye. Endocrinology. 2001, 142: 5158-5166. 10.1210/en.142.12.5158.
Baudet ML, Martin B, Hassanali Z, Parker E, Sanders EJ, Harvey S: Expression, translation, and localization of a novel, small growth hormone variant. Endocrinology. 2007, 148: 103-115.
Nam SY, Lobie PE: The mechanism of effect of growth hormone on preadipocyte and adipocyte function. Obesity Rev. 2000, 1: 73-86. 10.1046/j.1467-789x.2000.00015.x.
Cogburn LA: Endocrine manipulation of body composition in broiler chickens. Crit Rev Poult Biol. 1991, 3: 283-305.
Moellers RF, Cogburn LA: Chronic intravenous infusion of chicken growth hormone increases body fat content of young broiler chickens. Comp Biochem Physiol. 1995, 110A: 47-56.
Cogburn LA, Tang J-S, Cui J, Sofer L, Leclercq B, Simon J, Burnside J: DNA microarray analysis of gene expression in liver of broiler chickens divergently selected for growth rate. Poult Sci. 2000, 79 (Suppl. 1): 72-
Wang X, Carré W, Saxton A, Cogburn LA: Manipulation of thyroid status and/or GH injection alters hepatic gene expression in the juvenile chicken. Cytogenet Genome Res. 2007, 117: 174-188. 10.1159/000103178.
Lattka E, Eggers S, Moeller G, Heim K, Weber M, Mehta D, Prokisch H, Illig T, Adamski J: A common FADS2 promoter polymorphism increases promoter activity and facilitates binding of transcription factor ELK1. J Lipid Res. 2010, 51: 182-191. 10.1194/jlr.M900289-JLR200.
Yin L, Zhang Y, Charron T, Hillgartner FB: Thyroid hormone, glucagon, and medium-chain fatty acids regulate transcription initiated from promoter 1 and promoter 2 of the acetyl-CoA carboxylase-[alpha] gene in chick embryo hepatocytes. Biochimica et Biophysica Acta (BBA) - Gene Structure and Expression. 2000, 1517: 91-99. 10.1016/S0167-4781(00)00267-0.
Yin L, Zhang Y, Hillgartner FB: Sterol regulatory element-binding protein-1 interacts with the nuclear thyroid hormone receptor to enhance acetyl-CoA aarboxylase-a transcription in hepatocytes. J Biol Chem. 2002, 277: 19554-19565. 10.1074/jbc.M111771200.
Kim CW, Moon YA, Park SW, Cheng D, Kwon HJ, Horton JD: Induced polymerization of mammalian acetyl-CoA carboxylase by MIG12 provides a tertiary level of regulation of fatty acid synthesis. Proc Natl Acad Sci USA. 2010, 107: 9626-9631. 10.1073/pnas.1001292107.
Wang X, Carré W, Zhou H, Lamont SJ, Cogburn LA: Duplicated Spot 14 genes in the chicken: characterization and identification of polymorphisms associated with abdominal fat traits. Gene. 2004, 332: 79-88.
D’Andre Hirwa C, Yan W, Wallace P, Nie Q, Luo C, Li H, Shen X, Sun L, Tang J, Li W, Zhu X, Yang G, Zhang X: Effects of the thyroid hormone responsive spot 14a gene on chicken growth and fat traits. Poult Sci. 2010, 89: 1981-1991. 10.3382/ps.2009-00582.
Cao ZP, Wang SZ, Wang QG, Wang YX, Li H: Association of Spot14a gene polymorphisms with body weight in the chicken. Poult Sci. 2007, 86: 1873-1880.
Zhan K, Hou ZC, Li HF, Xu GY, Zhao R, Yang N: Molecular cloning and expression of the duplicated thyroid hormone responsive Spot 14 (THRSP) genes in ducks. Poult Sci. 2006, 85: 1746-1754.
Su S, Zhu H, Li Q, Xie Z: Molecular cloning and sequence analysis of Spot 14 alpha in geese. Brit Poult Sci. 2009, 50: 459-466. 10.1080/00071660903110893.
Graugnard DE, Piantoni P, Bionaz M, Berger LL, Faulkner DB, Loor JJ: Adipogenic and energy metabolism gene networks in longissimus lumborum during rapid post-weining growth in Angus x Simmental cattle fed high-starch or low-starch diets. BMC Genomics. 2009, 10: 142-10.1186/1471-2164-10-142.
Graugnard DE, Berger LL, Faulkner DB, Loor JJ: High starch diets induce precocious adipogenic gene network up-regulation in longissimus lumborum of early-weaned Angus cattle. Brit J Nutr. 2010, 103: 953-963. 10.1017/S0007114509992789.
Harvatine KJ, Bauman DE: SREBP1 and thyroid hormone responsive Spot 14 (S14) are involved in the regulation of bovine mammary lipid synthesis during diet-induced milk fat depression and treatment with CLA. J Nutr. 2006, 136: 2468-2474.
Piantoni P, Bionaz M, Graugnard D, Daniels K, Everts R, Rodriguez-Zas S, Lewin H, Hurley H, Akers M, Loor J: Functional and gene network analyses of transcriptional signatures characterizing pre-weaned bovine mammary parenchyma or fat pad uncovered novel inter-tissue signaling networks during development. BMC Genomics. 2010, 11: 331-10.1186/1471-2164-11-331.
Bauman DE, Harvatine KJ, Lock AL: Nutrigenomics, rumen-derived bioactive fatty acids, and the regulation of milk fat synthesis. Ann Rev Nutr. 2011, 31: 299-319. 10.1146/annurev.nutr.012809.104648.
Zhu Q, Anderson GW, Mucha GT, Parks EJ, Metkowski JK, Mariash CN: The Spot 14 protein is required for de novo lipid synthesis in the lactating mammary gland. Endocrinology. 2005, 146: 3343-3350. 10.1210/en.2005-0204.
Anderson GW, Zhu Q, Metkowski J, Stack MJ, Gopinath S, Mariash CN: The Thrsp null mouse (Thrsptm1cnm) and diet-induced obesity. Mol Cell Endocrinol. 2009, 302: 99-107. 10.1016/j.mce.2009.01.005.
Moncur JT, Park JP, Memoli VA, Mohandas TK, Kinlaw WB: The “Spot 14” gene resides on the telomeric end of the 11q13 amplicon and is expressed in lipogenic breast cancers: Implications for control of tumor metabolism. Proc Natl Acad Sci. 1998, 95: 6989-6994. 10.1073/pnas.95.12.6989.
Kinlaw WB, Quinn JL, Wells WA, Roser-Jones C, Moncur JT: Spot 14: A marker of aggressive breast cancer and a potential therapeutic target. Endocrinology. 2006, 147: 4048-4055. 10.1210/en.2006-0463.
Kirschner LS, Mariash CN: Adipose S14 mRNA is abnormally regulated in obese subjects. Thyroid. 1999, 9: 143-148. 10.1089/thy.1999.9.143.
Kinlaw WB, Church JL, Harmon J, Mariash CN: Direct evidence for the role of the “Spot 14” protein in the regulation of lipid synthesis. J Biol Chem. 1995, 270: 16615-16618. 10.1074/jbc.270.28.16615.
Jump DB, Clarke SD, MacDougald O, Thelen A: Polyunsaturated fatty acids inhibit S14 gene transcription in rat liver and cultured hepatocytes. Proc Natl Acad Sci U S A. 1993, 90: 8454-8458. 10.1073/pnas.90.18.8454.
Kinlaw WB, Schwartz HL, HAMBLIN PS, Mariash CN, Oppenheimer JH: Triiodothyronine rapidly reverses inhibition of S14 gene transcription by glucagon. Endocrinology. 1988, 123: 2255-2260. 10.1210/endo-123-5-2255.
Breuker C, Moreau A, Lakhal L, Tamasi V, Parmentier Y, Meyer U, Maurel P, Lumbroso S, Vilarem MJ, Pascussi JM: Hepatic expression of thyroid hormone-responsive Spot 14 protein is regulated by constitutive androstane receptor (NR1I3). Endocrinology. 2010, 151: 1653-1661. 10.1210/en.2009-1435.
Liu H-C, Towle HC: Functional synergism between multiple thyroid hormone response elements regulates hepatic expression of the rat S14 gene. Mol Endocrinol. 1994, 8: 1021-1037. 10.1210/me.8.8.1021.
Darras VM, Van Herck SL: Iodothyronine deiodinase structure and function: from ascidians to humans. J Endocrinol. 2012, 215: 189-206. 10.1530/JOE-12-0204.
Chutkow WA, Patwari P, Yoshioka J, Lee RT: Thioredoxin-interacting protein (Txnip) is a critical regulator of hepatic glucose production. J Biol Chem. 2008, 283: 2397-2406.
Blouet C, Liu SM, Jo YH, Chua S, Schwartz GJ: TXNIP in Agrp neurons regulates adiposity, energy expenditure, and central leptin sensitivity. J Neurosci. 2012, 32: 9870-9877. 10.1523/JNEUROSCI.0353-12.2012.
Margolskee RF, Dyer J, Kokrashvili Z, Salmon KSH, Ilegems E, Daly K, Maillet EL, Ninomiya Y, Mosinger B, Shirazi-Beechey SP: T1R3 and gustducin in gut sense sugars to regulate expression of Na+−glucose cotransporter 1. Proc Natl Acad Sci. 2007, 104: 15075-15080. 10.1073/pnas.0706678104.
Ren X, Zhou L, Terwilliger R, Newton SS, de Araujo IE: Sweet taste signaling functions as a hypothalamic glucose sensor. Front Integr Neurosci. 2009, 3: 12-
Mace OJ, Lister N, Morgan E, Shepherd E, Affleck J, Helliwell P, Bronk JR, Kellett GL, Meredith D, Boyd R, Pieri M, Bailey PD, Pettcrew R, Foley D: An energy supply network of nutrient absorption coordinated by calcium and T1R taste receptors in rat small intestine. J Physiol. 2009, 587: 195-210. 10.1113/jphysiol.2008.159616.
Yoshinari K, Sato T, Okino N, Sugatani J, Miwa M: Expression and induction of cytochromes P450 in rat white adipose tissue. J Pharmacol Exp Ther. 2004, 311: 147-154. 10.1124/jpet.104.067066.
Ibdah JA, Paul H, Zhao Y, Binford S, Salleng K, Cline M, Matern D, Bennett MJ, Rinaldo P, Strauss AW: Lack of mitochondrial trifunctional protein in mice causes neonatal hypoglycemia and sudden death. J Clin Invest. 2001, 107: 1403-1409. 10.1172/JCI12590.
Baes M, Huyghe S, Carmeliet P, Declercq PE, Collen D, Mannaerts GP, Van Veldhoven PP: Inactivation of the peroxisomal multifunctional protein-2 in mice impedes the degradation of not only 2-methyl-branched fatty acids and bile acid intermediates but also of very long chain fatty acids. J Biol Chem. 2000, 275: 16329-16336. 10.1074/jbc.M001994200.
Baker ME: Evolution of 17ß-hydroxysteroid dehydrogenases and their role in androgen, estrogen and retinoid action. Mol Cell Endocrinol. 2001, 171: 211-215. 10.1016/S0303-7207(00)00414-7.
Holness MJ, Sugden MC: Regulation of pyruvate dehydrogenase complex activity by reversible phosphorylation. Biochem Soc Trans. 2003, 31: 1143-1151. 10.1042/BST0311143.
Goldstein I, Rotter V: Regulation of lipid metabolism by p53 – fighting two villains with one sword. Trends Endocrinol Metab. 2012, 23: 567-575. 10.1016/j.tem.2012.06.007.
This work was supported by a grant from the United States Department of Agriculture, Initiative for Future Agricultural and Food Systems, Animal Genome Program (USDA-IFAFS; Award # 00-52100-9614) to LAC, TEP, SEA, and JS. The longitudinal study and tissue sampling of the divergent lines required the contribution of a large number of scientists and technicians from the laboratories of INRA UR83 and breeding facilities at INRA UE1295 Pôle d’Expérimentation Avicole de Tours, F-37380 Nouzilly, France. And although not named, each individual’s essential contribution toward completion of this project is gratefully acknowledged. Portions of this work were also supported by a grant from the USDA Cooperative State Research, Education, and Extension Service (CSREES) (Grant # 2009-34562-20008) and the Avian Biosciences Center, University of Delaware.
The authors declare that they have no competing interests.
LAC, SEA, JS, ELBD, MJD and TP secured funding, designed and carried out the animal experiment. LAC, WC, XW and TEP designed, constructed and printed the Del-Mar 14 K Chicken Integrated Systems microarrays. WC isolated RNA from the 48 abdominal fat samples, labeled, hybridized and scanned the microarrays. CWR performed statistical analysis of microarray data, completed the Ingenuity Pathway Analysis (IPA) and qRT-PCR analyses, and wrote the first draft of the manuscript. All coauthors reviewed, revised, commented on and approved the final version of the manuscript.
Wilfrid Carré, Xiaofei Wang contributed equally to this work.
Electronic supplementary material
Additional file 2:Description of the primers used for the qRT-PCR analysis. A Microsoft Excel file containing a single work sheet “Primer information” provides the gene symbol, gene name, forward and reverse primer sequences, GenBank accession and amplicon size (bp) for each primer used for qRT-PCR analysis. (XLSX 17 KB)
Additional file 3:Differentially expressed (DE) gene lists. A Microsoft Excel file containing three work sheets (“Main Effect of Genotype”, “Age X Genotype Interaction” and “Main Effect of Age”) that provide information about DE genes identified as the main effect of genotype, age x genotype interaction, or the main effect of age, respectively. Each list provides the clone ID, gene symbol, gene name, Entrez protein ID, log2 ratio (FL/LL) and FDR adjusted P-value for each DE gene. (XLSX 278 KB)
Additional file 4:IPA summary of DE genes related to diseases and disorders. A Microsoft Excel file containing six work sheets. “Overall IPA Summary” presents a synopsis of functionally annotated diseases and disorders in abdominal fat of FL-LL chickens. This work sheet provides information on IPA functional annotation, P-value for over-representation of DE genes, a list of genes and the number of genes assigned to each biological function. Five worksheets provide annotation of DE genes listed in “Overall IPA Summary”. Each worksheet provides the Entrez protein ID, gene symbol and log2 ratio for functionally annotated groups of genes. (XLSX 23 KB)
Additional file 5:IPA canonical pathways in abdominal fat of FL and LL chickens. A Microsoft Excel file containing a single worksheet “Canonical Pathways”. The top of the worksheet provides the canonical pathway, P-value associated with over-representation and ratio (number of genes present in dataset/number of known genes assigned to that pathway in the Ingenuity® Knowledge Base. Below the summary, a table for each canonical pathway provides the gene symbol and log2 expression ratio (FL/LL). (XLSX 12 KB)
Additional file 6:IPA annotated DE genes in their molecular and cellular function. This Microsoft Excel file contains two worksheets: “IPA Summary” presents nine functional categories of DE genes found in abdominal fat of the FL and LL chickens that are related to “Lipid Metabolism”. This worksheet provides the biological function, P-values of over-representation, the number of genes, and a gene list for each functional category. Four of the functionally annotated groups (bold type) are expanded in the second worksheet “Lipid Metabolism” to include protein ID, gene symbol, and log2 expression ratio (FL/LL). (XLSX 15 KB)
Additional file 7:Verification of differential expression of transcription factors by qRT-PCR analysis. This PowerPoint slide file shows qRT-PCR analysis of six transcription factors. Each data point represents LSMEANS (n = 4 birds/genotype) of normalized expression values. A two-factor ANOVA was used to determine significance (P≤0.05). The shaded box in each panel indicates significant effects of age (A), genotype (G) and/or the A x G interaction; the parenthesis shows the common standard error (SE) of LSMEANS for that gene as determined by the GLM procedure in SAS. (PPTX 139 KB)
Additional file 8:qRT-PCR analysis of genes involved in thyroid hormone and retinol metabolism and signaling. The abundance of genes involved in signaling and metabolism of thyroid hormone (A., left side) and retinol (B., right side) was verified by quantitative reverse transcription PCR (qRT-PCR) analysis. Each data point represents LSMEANS (n = 4 birds/genotype) of normalized expression values. A two-factor ANOVA was used to determine significance (P≤0.05). The shaded box in each panel indicates significant effects of age (A), genotype (G) and/or the A x G interaction; the parenthesis shows the common standard error (SE) of LSMEANS for that gene as determined by the GLM procedure in SAS. (PPTX 164 KB)
Authors’ original submitted files for images
Below are the links to the authors’ original submitted files for images.
About this article
Cite this article
Resnyk, C.W., Carré, W., Wang, X. et al. Transcriptional analysis of abdominal fat in genetically fat and lean chickens reveals adipokines, lipogenic genes and a link between hemostasis and leanness. BMC Genomics 14, 557 (2013). https://doi.org/10.1186/1471-2164-14-557
- Transcriptional regulators
- Hemostatic genes
- Retinoic acid signaling
- Thyroid hormone action
- Polygenic trait
- Visceral obesity
- Gene interaction networks
- Canonical metabolic/regulatory pathways