- 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
- Christopher W Resnyk1,
- Wilfrid Carré†1, 2,
- Xiaofei Wang†1, 3,
- Tom E Porter4,
- Jean Simon5,
- Elisabeth Le Bihan-Duval5,
- Michael J Duclos5,
- Sam E Aggrey6 and
- Larry A Cogburn1Email author
https://doi.org/10.1186/1471-2164-14-557
© Resnyk et al.; licensee BioMed Central Ltd. 2013
- Received: 29 April 2013
- Accepted: 12 August 2013
- Published: 16 August 2013
Abstract
Background
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).
Results
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.
Conclusions
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.
Keywords
- Adipogenesis
- Transcriptional regulators
- Hemostatic genes
- Lipogenesis
- Adipokines
- Retinoic acid signaling
- Thyroid hormone action
- Polygenic trait
- Visceral obesity
- Gene interaction networks
- Canonical metabolic/regulatory pathways
Background
The chicken was the first avian species and domestic animal selected for complete genome sequencing and assembly [1]. Subsequently, the chicken has emerged as a premier model in animal agriculture [2–4] and developmental biology [5]. Although now recognized as a model organism for biomedical research [6], 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 [12] and, like humans [13], 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 [20]. 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 [21]. 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 [22]. 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 [23]. 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 [26]. 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 [29], quantitative RT-PCR [30, 31] or targeted low-density array [32]. Our preliminary analysis of the liver transcriptome in the FL and LL chickens during juvenile development revealed 1,805 differentially expressed (DE) genes [3]. 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 [19].
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.
Methods
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.
Microarray analysis
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 [34], 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 [35].
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 [36] package in R (version 2.11.1) [37]. 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 [38] 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 [39] 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 [40]. 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.
Results
Phenotypic measurements
Phenotypic measurements from juvenile FL and LL cockerels
Age (wk) | ||||||
|---|---|---|---|---|---|---|
1 | 3 | 5 | 7 | 9 | 11 | |
Body weight (kg) | ||||||
FL | 0.115 | 0.544 | 1.297 | 1.983 | 2.693 | 3.222 |
LL | 0.123 | 0.551 | 1.204 | 1.964 | 2.787 | 3.281 |
Abdominal fat (g) | ||||||
FL | 0.5 | 13* | 38* | 88* | 124* | 150* |
LL | 0.4 | 5 * | 15* | 31* | 54* | 59* |
FL/LL ratio | 1.2 | 2.6 | 2.5 | 2.8 | 2.3 | 2.5 |
Abdominal fat (%BW) | ||||||
FL | 0.4 | 2.3* | 2.9* | 4.4* | 4.6* | 4.6* |
LL | 0.3 | 1.0* | 1.2* | 1.6* | 1.9* | 1.8* |
FL/LL ratio | 1.3 | 2.3 | 2.4 | 2.8 | 2.4 | 2.6 |
Abdominal fat gene expression
Venn diagram showing unique and shared genes among main effect of age (A) or genotype (G), and their interaction (A × G). This diagram shows the number of differentially expressed (DE) genes that are common across contrasts and those that are unique to G (P≤0.05), A (P≤0.001), or the A × G interaction (P≤0.05).
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 [41], 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 [42] 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.
Top biological functions of DE genes in abdominal fat of juvenile FL and LL chickens*
Diseases and disorders | P-value | # Genes | |
|---|---|---|---|
Developmental disorder | 2.76E-07 | 33 | |
Hereditary disorder | 3.01E-06 | 71 | |
Inflammatory disease | 7.14E-06 | 11 | |
Organismal injury and abnormalities | 4.51E-05 | 31 | |
Metabolic disease | 4.77E-05 | 41 | |
Molecular and cellular functions | |||
Lipid metabolism | 6.06E-05 | 46 | |
Small Molecule biochemistry | 6.06E-05 | 43 | |
Physiological system development and function | |||
Hematological system development and function | 1.87E-05 | 34 | |
Organ morphology | 2.56E-05 | 7 | |
Renal system development/function | 2.56E-05 | 10 | |
Embryonic development | 1.06E-04 | 23 | |
Cardiovascular system function | 1.75E-04 | 8 | |
Top canonical pathways | P -Value | Genes † | Ratio † |
Coagulation system | 2.56E-08 | (7/38) | 0.184 |
Intrinsic prothrombin activation pathway | 1.75E-04 | (6/34) | 0.176 |
Extrinsic prothrombin activation pathway | 4.65E-04 | (3/20) | 0.15 |
Acute phase response signaling | 5.11E-08 | (15/178) | 0.08 |
Gene interaction network in abdominal fat of LL chickens associated with hemostasis. Functional gene interactions networks were identified by Ingenuity Pathway Analysis (IPA®) software. This network shows direct gene interactions mainly in abdominal fat of LL chickens related to “Hematological System Development and Function” (A). The IPA® Upstream Regulator Analysis identified transcription factors with direct actions on differentially expressed target genes in abdominal fat of FL and LL chickens. This analysis of upstream regulators (based on expected responses from literature and observed responses in the data set) predicts inhibition (blue color) of hepatic nuclear factor 1A (HNF1A) (B) and peroxisome proliferator-activated receptor gamma (PPARG) (C), which would lead to inhibition (blue edges or lines) of target gene expression. Red gene symbols indicate higher expression in the FL and green gene symbols indicate higher expression in the LL.
Verification of differential expression of hemostatic genes by qRT-PCR analysis. The abundance of six genes associated with blood coagulation was determined by quantitative reverse transcription PCR (qRT-PCR) analysis. Data points represent Least Squares Means (LSMEANS; n = 4 birds/genotype) of normalized expression values generated by the general linear models (GLM) procedure in Statistical Analysis System (SAS) software. A two-factor (genotye and age) analysis of variance (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.
Verification of differential expression of adipokines by qRT-PCR analysis. The abundance of eight adipokines was determined by quantitative reverse transcription PCR (qRT-PCR) analysis. Data points represent 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 × G interaction; the parenthesis shows the common standard error (SE) of LSMEANS for that gene determined by the GLM procedure in SAS.
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
Transcriptional regulation of gene interaction network in abdominal fat of FL and LL chickens controlling lipogenesis. Functional gene interactions and up-stream regulators were identified by IPA (gene symbols and color schemes as described in Figure 2). Direct interactions (solid lines) were found among transcription regulators [JUN, SREBF1, SIRT2, MID1IP1, NR1H2 (LXRB) and THRSP] and lipogenic genes (A). ✝THRSP and its paralog MID1IP1 (or THRSP-like, THRSPL) were not con sidered as transcription regulators by IPA software. THRSPA was added to this network based on microarray and qRT-PCR analysis (Additional files 3 and 7) and known involvement of THRSP in regulating expression of lipogenic enzymes across multiple species of birds and mammals. This analysis of upstream regulators predicts activation of JUN (B) and sterol response element binding factor 1 (SREBF1) (C), which would lead to inhibition or activation [orange edges (lines)] expression of DE target genes. Gene symbol color indicates higher expression in the FL (red) or higher expression (green) in the LL.
Expression of genes associated with lipid metabolism by qRT-PCR analysis. mRNA expressions of eight genes involved in lipid metabolism were determined by quantitative reverse transcription PCR (qRT-PCR). 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 × G interaction; the parenthesis shows the common standard error (SE) of LSMEANS for that gene determined by the GLM procedure in SAS.
Transcriptional regulators of DE genes controlling lipogenesis in abdominal fat of FL and LL chickens. A large number of DE lipogenic genes interact with two transcriptional regulators, SIRT2 and PPARD (A). The IPA Upstream Regulator Analysis (B) predicts that the up-regulation of SIRT2 leads to activation of five lipogeneic genes (orange-edged arrows), whereas, the predicted inhibition of PPARD would lead to down regulation (blue-edged arrows) of seven DE target genes in the FL [or up-regulation (green gene symbols) in the LL]. The predicted activation of major lipogenic genes (ALDH2, CCL13, FASN and SCD) would be blocked (blunt orange arrow) by PPARD.
Ligand activated nuclear receptors and other transcription factors
Functional categories of DE and prior candidate genes expressed in abdominal fat and the average fold change (FL/LL) as determined by microarray and/or qRT-PCR analyses
Functional category | Symbol | Gene name | Microarray (FC)* | qRT-PCR (FC)* |
|---|---|---|---|---|
Hemostasis | A2M | Alpha-2-macroglobulin | −1.89 | −1.10 |
AGT | Angiotensinogen | 1.20 | - | |
ANG | Angiogenin | −2.51 | - | |
CFB | Complement factor B | −1.49 | - | |
CPB2 | Carboxypeptidase B2 | −1.43 | - | |
CPM | Carboxypeptidase M | −1.32 | - | |
F2 | Thrombin | −1.85 | −1.35 | |
F9 | Christmas factor | −1.51 | −4.04 | |
FGA | Fibrinogen alpha | −2.61 | - | |
PLG | Plasminogen | −1.79 | - | |
PROC | Protein C | −1.39 | −3.54 | |
SERPINA1 | Antitrypsin | −1.75 | - | |
SERPIND1 | Heparin cofactor | −2.00 | −1.50 | |
THBS2 | Thrombospondin 2 | −1.17 | - | |
Adipokines | ADIPOQ | Adiponectin | 1.03 | −1.48 |
ANGPTL4 | Angiopoietin-like 4 | 1.01 | −1.58 | |
ATRN | Attractin | −1.12 | −1.22 | |
CFD | Adipsin | 1.24 | - | |
LPL | Lipoprotein lipase | - | −1.41 | |
NAMPT | Visfatin | - | −1.20 | |
RARRES2 | Chemerin | −1.32 | −1.54 | |
RBP4 | Retinol binding protein 4 | −2.33 | −1.11 | |
Lipogenesis | DHCR7 | 7-Dehydrocholesterol reductase | 1.11 | - |
FADS2 | Fatty acid desaturase 2 | 1.21 | - | |
FASN | Fatty acid synthase | 1.36 | 1.60 | |
G6PC | Glucose-6-phosphatase, catalytic subunit | 1.46 | - | |
scGH | Growth hormone, chicken, short form | 1.15 | - | |
HMGCR | 3-Hydroxy-3-methylglutaryl-CoA reductase | 1.09 | - | |
INSIG2 | Insulin induced gene 2 | 1.74 | - | |
LCAT | Lecithin-cholesterol acyltransferase | 1.32 | - | |
MVD | Mevalonate (diphospho) decarboxylase | 1.20 | - | |
SCD | Stearoyl-CoA desaturase | 1.48 | 1.88 | |
SREBF1 | Sterol regulatory element binding transcription factor 1 | 1.12 | 1.32 | |
THRSPA | Thyroid hormone responsive spot 14 A | - | 1.64 | |
TXNIP | Thioredoxin interacting protein | - | 1.80 | |
Lipolysis | ACAT1 | Acetyl-CoA acetyltransferase 1 | −3.18 | - |
ADH1C | Alcohol dehydrogenase 1C (class I), gamma polypeptide | −1.81 | - | |
APOA1 | Apolipoprotein A-I | −1.16 | - | |
APP | Amyloid beta (A4) precursor protein | −1.15 | - | |
BCMO1 | beta-carotene 15,15′-monooxygenase | - | −1.13 | |
BCO2 | beta-carotene oxygenase 2 | −1.15 | −1.48 | |
CYP27A1 | Cytochrome P450, family 27, subfamily A, polypeptide 1 | −1.14 | - | |
CYP2E1 | Cytochrome P450, family 2, subfamily E, polypeptide 1 | −1.80 | - | |
EHHADH | Enoyl-CoA, hydratase/3-hydroxyacyl CoA dehydrogenase | −1.09 | - | |
GAMT | Guanidinoacetate N-methyltransferase | −1.21 | - | |
HADHB | Hydroxyacyl-CoA dehydrogenase (trifunctional protein) | −1.10 | - | |
HSD17B4 | Hydroxysteroid (17-beta) dehydrogenase 4 | −1.92 | - | |
HSD17B6 | Hydroxysteroid (17-beta) dehydrogenase 6 | −1.19 | - | |
IRS1 | Insulin receptor substrate 1 | −1.59 | - | |
PDK4 | pyruvate dehydrogenase kinase, isozyme 4 | - | −1.99 | |
PHYH | Phytanoyl-CoA 2-hydroxylase | −1.56 | - | |
SLC2A2 | Facilitated glucose transporter 2 (GLUT2) | −2.23 | - | |
SOD3 | superoxide dismutase 3, extracellular | −1.10 | −1.20 | |
TP53 | Tumor protein p53 | −1.29 | - | |
UCP3 | Uncoupling protein 3 (mitochondrial, proton carrier) | −1.21 | - |
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).
Gene interaction network of nuclear receptors, co-activators and regulators of gene transcription in abdominal fat of juvenile FL and LL chickens. This gene network shows direct interactions of seven transcriptional regulators [CEBPZ, RXRG, NR1H4 or farnesoid X receptor (FXR), NCOA1 or steroid receptor coactivator 1 (SRC-1), THRA, THRSP and MID1IP1 (or THRSP-like, THRSPL)] and their target genes. Gene symbol color indicates higher expression in the FL (red) or higher expression (green) in the LL.
Upstream regulators of gene transcription in abdominal fat of juvenile FL and LL chickens. Ingenuity® Upstream Regulator Analysis revealed a large number of transcriptional regulators (see Table 4) controlling lipid metabolism genes in abdominal fat (A). This IPA analysis shows “up-stream regulators” and their downstream targets found among DE fatty acid metabolism genes identified in abdominal fat of the FL and LL chickens. Differentially expressed gene targets regulated by six additional transcription factors are shown (B). The IPA prediction of activation (orange lines and symbols) or (blue lines and symbols) inhibition states is based on prior knowledge accrued by Ingenuity® Knowledge Base and expression values of differentially expressed genes identified by microarray analysis of abdominal fat in juvenile FL and LL chickens. Gene symbol color indicates higher expression in the FL (red) or higher expression (green) in the LL.
Transcriptional regulators of genes that control the divergence of abdominal fatness in FL and LL chickens
Symbol | NCBI Entrez gene name | Activation z-score | P-value of overlap | # Target molecules |
|---|---|---|---|---|
CEBPA | CCAAT/enhancer binding protein (C/EBP), alpha | 0.379 | 1.61E-15 | 34 |
CEBPB | CCAAT/enhancer binding protein (C/EBP), beta | 1.935 | 1.75E-10 | 25 |
CREB1 | cAMP responsive element binding protein 1 | 0.527 | 1.64E-05 | 15 |
NR0B2 | nuclear receptor subfamily 0, group B, member 2 (SHP) | −0.84 | 3.74E-08 | 10 |
NR1H2 | nuclear receptor subfamily 1, group H, member 2 (LXRB) | −1.512 | 2.64E-10 | 12 |
NR1H3 | nuclear receptor subfamily 1, group H, member 3 (LXRA) | −0.2 | 4.69E-11 | 14 |
NR1H4 | nuclear receptor subfamily 1, group H, member 4 (FXR) | 1.076 | 1.40E-06 | 11 |
NR5A2 | nuclear receptor subfamily 5, group A, member 2 (LRH1) | −1.412 | 3.57E-04 | 7 |
PPARA | peroxisome proliferator-activated receptor alpha | −1.339 | 1.78E-32 | 54 |
PPARD | peroxisome proliferator-activated receptor delta | −0.767 | 3.13E-13 | 21 |
PPARG | peroxisome proliferator-activated receptor gamma | −1.629 | 1.38E-27 | 48 |
PPARGC1B | peroxisome proliferator-activated receptor gamma, coactivator 1 beta | 1.488 | 1.37E-12 | 12 |
RXRA | retinoid X receptor, alpha | −0.932 | 3.32E-18 | 31 |
SREBF1 | sterol regulatory element binding transcription factor 1 | 0.511 | 2.91E-20 | 29 |
SREBF2 | sterol regulatory element binding transcription factor 2 | 1.171 | 2.78E-15 | 17 |
THRA | thyroid hormone receptor, alpha | −0.246 | 2.08E-08 | 12 |
THRB | thyroid hormone receptor, beta | 1.135 | 7.18E-11 | 19 |
Discussion
The divergent FL and LL chickens were originally developed as experimental models to study genetic and endocrine mechanisms controlling excessive abdominal fatness [43]. 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 [44]. The divergence of adiposity between the FL and LL chickens occurs at 3 wk of age [27]; hyperplasia of adipocytes was found as early as 2 wk of age in the FL [28], followed by marked hypertrophy of adipocytes by 9 wk of age [27]. 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 [44]. 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 [51]. 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 [52]. 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 [53]. 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 [54].
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 [3]. 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 [58]. 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 [58]. Adipocyte-derived chemerin causes insulin resistance in skeletal muscle cells [59]; and as a secreted adipokine, chemerin regulates myogenesis by providing negative cross-talk between adipose tissue and skeletal muscle [60]. 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 [64]. 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 [65]. 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 [66]. 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 [67]. In fact, ANGPTL4 is a potent irreversible inhibitor of lipoprotein lipase (LPL) activity, which leads to hypertriglyceridemia [68]. 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 [69]. 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) [75] 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 [79], 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 [82]. 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 [83]. 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 [84]. 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 [85]. 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 [86]. 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 [32]. 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[87]. 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 [3]. 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 [88]. 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] [89]. 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 [21]. 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 [22]. In contrast, more than a thousand genes were either differentially expressed in liver or leg muscle of the same birds following insulin immunoneutralization [20]. 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 [21].
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 [91]. The short alternatively-spliced (16.5 kDa) isoform of full length GH (20 kDa) appears to function as an “intracrine” factor within the cell [91]. 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 [92]. 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[97]. 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 [100]. 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 [101]. Mutations in the THRSPA gene of chickens [101–103], ducks [104] and geese [105] 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 [111] and pups from the THRSP-null mouse exhibit reduced body weight gain due to diminished milk triglycerides [112]. In humans, amplification of the THRSP locus is associated with lipogenic breast cancer [113]; and, as such, THRSP serves as a marker of aggressive breast cancer and a potential target of anti-cancer drugs [114]. 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 [115]. These observations support the idea that THRSP is a transcriptional activator of several lipogenic enzymes (ACLY, FASN and ME) in the mouse [116]. THRSP is activated in response to T3, glucose and insulin and inhibited by polyunsaturated fatty acids [117], cyclic AMP or glucagon [118]. 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 [119]. Another study showed that FASN co-precipitates with THRSP in nuclear extracts from the mouse (referenced in [110]). 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) [120]. 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) [121]. 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 [122] that also mediates hypothalamic control over energy utilization and adiposity in the mouse [123]. 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 [79] 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 [127]. 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 [128]. 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 [129]. Since the other 17β-HSD (HSD17B6) exhibits retinol dehydrogenase activity [130], 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 [131], 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 [3], 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 [21]. The tumor suppressor protein p53 enhances lipid catabolism and induces expression of guanidinoacetate N-methyltransferase (GAMT), which enhances β-oxidation of fatty acids [132]; 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 [3]. 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.
Conclusions
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.
Notes
Abbreviations
- FL:
-
Fat line chickens
- LL:
-
Lean line chickens
- DE:
-
Differentially expressed
- G:
-
Main effect of genotype
- A × G:
-
Age by genotype interaction
- A:
-
Main effect of age
- wk:
-
Week of age
- eQTL:
-
Expression quantitative trait loci
- GGA5:
-
Gallus gallus chromosome 5
- IPA:
-
Ingenuity Pathway Analysis
- INRA:
-
Institut National de la Recherche Agronomique
- qRT-PCR:
-
Quantitative reverse transcriptase polymerase chain reaction
- h:
-
Hour(s)
- L:
-
Light
- D:
-
Dark
- LIMMA:
-
Linear models for statistical analysis of microarray data
- FDR:
-
False discovery rate
- SAS:
-
Statistical Analysis System
- GLM:
-
General linear models
- LSMEANS:
-
Least Squares Means.
Declarations
Acknowledgements
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.
Authors’ Affiliations
References
- 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.Google Scholar
- Burt DW: Emergence of the chicken as a model organism: implications for agriculture and biology. Poult Sci. 2007, 86: 1460-1471.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Dodgson JB: The chicken genome: some good news and some bad news. Poult Sci. 2007, 86: 1453-1459.PubMedGoogle Scholar
- Stern CD: The chick: a great model system becomes even greater. Dev Cell. 2005, 8: 9-17.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Simon J: Chicken as a useful species for the comprehension of insulin action. Crit Rev Poult Biol. 1989, 2: 121-148.Google Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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-Google Scholar
- Goodridge AG, Ball EG: Lipogenesis in the pigeon: in vivo studies. Amer J Physiol. 1967, 213: 245-249.PubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Life Technologies.http://www.invitrogen.com/,
- Smyth GK: Limma: Linear models for microarray data. 2005, New York: Springer, 397-420.Google Scholar
- Smyth GK: Linear models and empirical bayes methods for assessing differential expression in microarray experiments. Stat Appl Genet Mol Biol. 2004, 3: Article3-PubMedGoogle Scholar
- 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.Google Scholar
- 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-PubMed CentralPubMedGoogle Scholar
- Cogburn Laboratory.http://cogburn.dbi.udel.edu/,
- Ingenuity Systems.http://www.ingenuity.com/,
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Saadoun A, Leclercq B: In vivo lipogenesis in genetically fat and lean chickens of various ages. Comp Biochem Physiol. 1986, 83B: 607-611.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Allman-Farinelli MA: Obesity and venous thrombosis: a review. Semin Thromb Hemost. 2011, 37: 903-907. 10.1055/s-0031-1297369.PubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- Palomo I, Alarcon M, Moore-Carrasco R, Argiles JM: Hemostasis alterations in metabolic syndrome (review). Int J Mol Med. 2006, 18: 969-974.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Yoshimura T, Oppenheim JJ: Chemerin reveals its chimeric nature. J Exp Med. 2008, 205: 2187-2190. 10.1084/jem.20081736.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Sharp PJ, Dunn IC, Waddington D: Chicken leptin. Gen Comp Endocrinol. 2008, 158: 2-4. 10.1016/j.ygcen.2008.05.018.PubMedGoogle Scholar
- 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-PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- Cogburn LA: Endocrine manipulation of body composition in broiler chickens. Crit Rev Poult Biol. 1991, 3: 283-305.Google Scholar
- 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.Google Scholar
- 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-Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.Google Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Kirschner LS, Mariash CN: Adipose S14 mRNA is abnormally regulated in obese subjects. Thyroid. 1999, 9: 143-148. 10.1089/thy.1999.9.143.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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-PubMed CentralPubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMed CentralPubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
- Holness MJ, Sugden MC: Regulation of pyruvate dehydrogenase complex activity by reversible phosphorylation. Biochem Soc Trans. 2003, 31: 1143-1151. 10.1042/BST0311143.PubMedGoogle Scholar
- 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.PubMedGoogle Scholar
Copyright
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
