- Research article
- Open Access
Identification of differentially expressed genes and pathways between intramuscular and abdominal fat-derived preadipocyte differentiation of chickens in vitro
BMC Genomics volume 20, Article number: 743 (2019)
The distribution and deposition of fat tissue in different parts of the body are the key factors affecting the carcass quality and meat flavour of chickens. Intramuscular fat (IMF) content is an important factor associated with meat quality, while abdominal fat (AbF) is regarded as one of the main factors affecting poultry slaughter efficiency. To investigate the differentially expressed genes (DEGs) and molecular regulatory mechanisms related to adipogenic differentiation between IMF- and AbF-derived preadipocytes, we analysed the mRNA expression profiles in preadipocytes (0d, Pre-) and adipocytes (10d, Ad-) from IMF and AbF of Gushi chickens.
AbF-derived preadipocytes exhibited a higher adipogenic differentiation ability (96.4% + 0.6) than IMF-derived preadipocytes (86.0% + 0.4) (p < 0.01). By Ribo-Zero RNA sequencing, we obtained 4403 (2055 upregulated and 2348 downregulated) and 4693 (2797 upregulated and 1896 downregulated) DEGs between preadipocytes and adipocytes in the IMF and Ad groups, respectively. For IMF-derived preadipocyte differentiation, pathways related to the PPAR signalling pathway, ECM-receptor interaction and focal adhesion pathway were significantly enriched. For AbF-derived preadipocyte differentiation, the steroid biosynthesis pathways, calcium signaling pathway and ECM-receptor interaction pathway were significantly enriched. A large number of DEGs related to lipid metabolism, fatty acid metabolism and preadipocyte differentiation, such as PPARG, ACSBG2, FABP4, FASN, APOA1 and INSIG1, were identified in our study.
This study revealed large transcriptomic differences between IMF- and AbF-derived preadipocyte differentiation. A large number of DEGs and transcription factors that were closely related to fatty acid metabolism, lipid metabolism and preadipocyte differentiation were identified in the present study. Additionally, the microenvironment of IMF- and AbF-derived preadipocyte may play a significant role in adipogenic differentiation. This study provides valuable evidence to understand the molecular mechanisms underlying adipogenesis and fat deposition in chickens.
Chicken is generally accepted as one of the main protein sources worldwide. In the last several decades, meat quality has decreased as a result of genetic selection for growth rate and feed conversion. Intramuscular fat (IMF) content, an important factor influencing meat quality, contributes to multiple meat quality characteristics, such as flavour, tenderness and juiciness [1,2,3]. Abdominal fat (AbF) is an important carcass trait in chickens. A higher growth rate induces larger fiber diameters and lower IMF deposition, which severely deteriorates the quality of meat [4, 5]. However, the overemphasis on selection for a rapid growth rate leads to excessive fat accumulation, especially AbF accumulation . Excessive fat is often discarded as waste [6,7,8]. Reducing the levels of AbF and increasing the levels of IMF have therefore become a major breeding goals in the chickens industry [6, 9].
Previous studies have indicated that adipose tissues from different locations display unique physiological and biochemical characteristics [10,11,12]. In addition, glucose utilization and lipid metabolism mechanisms and hormone sensitivities are different among tissues from different locations [13,14,15,16]. IMF has specific biological features compared with fat from other locations. Previous studies have suggested that AbF has higher triglyceride (TG) levels than IMF tissue [17,18,19]. Hrdinka C et al. demonstrated that the fatty acid composition of AbF differs significantly from that of IMF . Zhou et al. found that dietary supplementation with 3% conjugated linoleic acid (CLA) decreased AbF accumulation but increased IMF content . Leng et al. indicated that a desirable broiler line with high IMF content but low AbF content could be obtained by genetic selection . This may be because AbF deposition and IMF deposition are subject to different regulatory mechanisms. However, the mechanisms underlying regional differences in chicken adipogenesis remain unknown.
In the current study, Ribo-Zero RNA-Seq was used to systematically identify differentially expressed genes (DEGs), mRNAs (DEMs) and novel genes (DENGs) and different pathways between preadipocytes and adipocytes of IMF and AbF. These data may contribute to a more thorough understanding of tissue-specific adipogenic differentiation and poultry meat quality.
Ribo-zero RNA-Seq of different chicken different adipose tissue-derived preadipocytes and adipocytes
Intramuscular and abdominal preadipocytes were isolated and cultured in growth medium until they reached 80–90% confluence (Fig. 1a). Microscopy showed that the IMF preadipocytes shared the same fibroblast-like morphology as the AbF preadipocytes (Fig. 1b). To construct intramuscular and abdominal adipogenic differentiation models, MDI medium supplemented with oleic acid was used for adipogenic differentiation. After induction with adipogenic agents for 10 days, chicken preadipocytes readily differentiated into mature adipocytes, and lipid droplets were visible under a microscope after 10 days of induction (Fig. 1b). AbF-derived preadipocytes exhibited a higher adipogenic differentiation ability (96.4% + 0.6) than IMF-derived preadipocytes (86.0% + 0.4) (p < 0.01). The expression level of the adipogenic marker genes peroxisome proliferator-activated receptor γ (PPARγ, PPARG) and fatty acid binding protein 4 (FABP4, ap2) significantly increased with adipogenic differentiation (p < 0.01) (Fig. 1c). Pearson correlation analysis showed that the gene expression correlation coefficient within each group was noticeably higher than that between the groups, reflecting a good linear correlation between the independent samples of preadipocytes or adipocytes in the IMF and AbF groups (Pearson correlation coefficient, r = 0.98) (Fig. 1d). Principal component analysis (PCA) also showed global differences among the preadipocyte, adipocyte, IMF and AbF groups (Fig. 1e). All evidence suggested that our data was repeatability and reproducibility.
Global analysis of gene expression patterns in chicken adipocytes
As shown in Table 1, 963,374,122 raw reads were produced from 8 cDNA libraries. We identified a considerable number of genes in preadipocytes and adipocytes derived from chicken breast muscle and abdominal tissues. The percentage of clean reads in each library ranged from 94.21 to 96.09%. The mapping percentage of the 8 samples ranged from 88.82 to 93.42%. We found that the genomic loci of the genes were widely distributed across chromosomes (Fig. 2a). The sequencing depth was saturated at 68 M reads for each library (Fig. 2b). The mapping percentages of the different samples on different regions of the genome are displayed in Fig. 2c. More than 75% of the reads were mapped to gene regions. In preadipocytes and adipocytes derived from breast muscle and abdominal tissues, less than 0.1% of the reads were mapped to splicing sites of the genome.
DEG, DEM and DENG profiles between different chicken adipose tissue-derived preadipocytes and adipocytes
To identify potential candidate genes related to adipogenic differentiation, we examined the expression level of genes in preadipocytes and adipocytes. A total of 2039 genes were found to be differentially expressed between preadipocytes and adipocytes in both the IMF and AbF groups (Fig. 3a). As expected, we noticed that large numbers of genes or transcription factors related to adipogenic differentiation and lipid metabolism were differentially expressed, including peroxisome proliferator-activated receptor γ (PPARG), bone morphogenetic protein 4 (BMP4), fatty acid synthase (FASN), adiponectin, C1Q and collagen domain-containing (ADIPOQ), perilipin 2 (PLIN2), lipin 1 (LPIN1), and carnitine palmitoyltransferase 1A (CPT1A) (Fig. 3b). The ten most abundant DEGs between preadipocytes and adipocytes in the IMF and AbF groups are presented in Table 2. Furthermore, we determined the global gene expression profiles in preadipocytes and adipocytes between different groups. The number of differentially expressed genes (DEGs) (Additional file 1: Table S1), mRNAs (DEMs) (Additional file 2: Table S2) and novel genes (DENGs) (Additional file 3: Table S3) between the different groups are shown in Fig. 4a. As illustrated in Fig. 4b, the expression levels of the genes of the two different samples exhibited a significant positive correlation (R2 > 0.81, p < 0.01) (Fig. 4b). We found that 2742 DEGs were upregulated and 1705 DEGs were downregulated in the AbAd group compared with the IMAd group, that 3437 DEGs were upregulated and 1310 DEGs were downregulated in the AbPre group compared with the IMPre group, that 2797 DEGs were upregulated and 1896 DEGs were downregulated in the AbPre group compared with the AbAd group, and that 2055 DEGs were upregulated and 2348 DEGs were downregulated in the IMPre group compared with the IMAd group (Fig. 4c).
Functional enrichment analysis of DEGs involved in the adipogenic differentiation of chicken preadipocytes
To investigate the functions of the DEGs in chicken adipogenic differentiation, GO (Additional file 4: Table S4) and KEGG pathway (Additional file 5: Table S5) analyses were performed in the present study. Our results suggested that the DEGs in the AbPre vs AbAd and IMPre vs IMAd comparisons were significantly enriched in ECM-receptor interaction, the PPAR signalling pathway, and focal adhesion (Fig. 5). Interestingly, we found that ABC transporters, glutathione metabolism and fatty acid biosynthesis were significantly enriched for the IMF groups, while steroid biosynthesis and the p53 signalling pathway were significantly enriched for the AbF groups. The DEGs in the IMPre vs AbPre and IMAd vs AbAd comparisons were significantly enriched in the focal adhesion, ECM-receptor interaction, and fatty acid metabolism pathways, among others (Fig. 5). To identify the gene expression patterns associated with adipogenic differentiation in both IMF- and AbF-derived adipocytes, the DEGs shared between the IMF group and the AbF group were analysed in the present study. Our results suggested that the shared DEGs were enriched for the metabolism, cellular processes and translation terms (Fig. 6a) and for pathways including ECM-receptor interaction, DNA replication, the cell cycle and the PPAR signalling pathway (Fig. 6b). Furthermore, we found that the shared downregulated DEGs (number: 750) were significantly enriched in the cell cycle, focal adhesion and purine metabolism pathways, while the upregulated DEGs (number: 792) were significantly enriched in the PPAR signalling pathway, the adipocytokine signalling pathway and the retinol metabolism pathway (Fig. 7).
Differentially expressed transcription factors and their potential interacting genes involved in adipogenic differentiation
To identify the differential gene expression patterns of transcription factors in the context of adipogenic differentiation, we compared the expression levels of transcription factors (TFs). Interestingly, our results suggested that most TFs shared the same gene expression patterns in both the IMF and AbF groups, while KLF9 and MYOG showed significantly different expression patterns between the groups (Fig. 8). Furthermore, as shown in Fig. 9, the TFs and genes related to adipogenic differentiation were significantly positively correlated (r > 0.86, p < 0.01).
Genes expressed in a fat depot-specific manner and the validation of the DEGs by qRT-PCR
To identify whether gene expression is fat depot-specific (Fig. 10), the expression levels of genes related to adipogenic differentiation at different differentiation stages and in different groups (IMF and AbF) were compared. Our results suggested that most genes have different gene expression patterns during the adipogenic differentiation process. Most genes reached the highest expression levels after 6 days of induced differentiation in the AbF group and after 8 days in the IMF group (Fig. 11). To confirm the results from RNA-Seq, qRT-PCR was performed in the present study. Our qRT-PCR results were consistent with the RNA-Seq results (Additional file 7: Figure S1).
Integrated analysis of DEG-pathway network between different fat-derived chicken preadipocytes and adipocytes
To further understand the adipogenic differentiation-regulated mechanisms of different fat-derived chicken preadipocytes, we visualized the integrated DEG-pathway networks for the IMF (Fig. 12a) and AbF (Fig. 12b) groups. Our results showed that a large number of DEGs in the IMF group were mainly enriched in the PPAR signaling pathway, fatty acid biosynthesis, focal adhesion and ECM-receptor interaction (Fig. 12a). For the AbF group, we noticed that the DEGs were mainly enriched in the PPAR signaling pathway, focal adhesion, ECM-receptor interaction, steroid biosynthesis and the p53 signaling pathway (Fig. 12b).
Fat deposition is mainly dependent on the proliferation, differentiation and maturation of preadipocytes . IMF content is an important factor that contributes to the tenderness, juiciness and flavour of meat and thus affects meat quality. High AbF content causes low slaughter efficiency [5, 7]. It is known that IMF has high genetic correlations with abdominal fat weight (AFW) and moderate correlations with AF percentage (AFP) in chickens . Previous studies have suggested that the proliferation and differentiation abilities of genes involved in lipid metabolism are dramatically lower in intramuscular adipocytes (IMAs) than in subcutaneous adipocytes (SAs) [11, 14, 24,25,26]. In the present study, lipids accumulated in chicken IMF and AbF preadipocytes at the late stage of differentiation. The oil red O staining results and the expression levels of two well-known adipogenic markers demonstrated that the model of adipocyte differentiation was successfully established. We noticed that the accumulation of lipids in chicken AbF adipocytes was higher than that in IMF adipocytes, consistent with the findings of previous studies in pigs [11, 27, 28].
Ribo-Zero RNA-Seq has been applied as an efficient method to explore transcriptional characteristics because it can capture both poly(A) + and poly(−) transcripts [9, 29, 30]. However, to our knowledge, none of the previous studies were carried out to compare gene expression profiles between IMF and AbF adipocytes, especially in chickens. The process by which preadipocytes differentiate into mature adipocytes is complex and is regulated by various transcription factors . Previous studies have identified a large number of transcription factors, including CCAAT/enhancer-binding protein (C/EBP), peroxisome proliferator-activated receptors (PPARs) and sterol regulatory element-binding protein (SREBP) [32,33,34,35].
In this study, we detected global gene expression profiles in preadipocytes and adipocytes, providing large amounts of information for further studies on the regulatory mechanisms underlying poultry adipogenic differentiation and tissue-specific fat deposition. Our results revealed that many genes related to lipid metabolism, such as matrix metallopeptidase 2 (MMP2), extracellular fatty acid-binding protein (EXFABP), CD36, prostaglandin D2 synthase (PTGDS), chondroadherin (CHAD), laminin alpha 2 (LAMA2), bone morphogenetic protein 4 (BMP4) and collagen type VI alpha 1 chain (COL6A1), were predominantly expressed in IMF adipocytes, whereas apolipoprotein A1 (APOA1), fatty acid-binding protein 4 (FABP4), perilipin 2 (PLIN2), fatty acid synthetase (FASN), pyruvate dehydrogenase kinase 4 (PDK4), collagen type IV alpha 1 chain (COL4A1) and glycerol-3-phosphate dehydrogenase 1-like 2 (GPD1L2) were highly expressed in AbF adipocytes; these findings suggest that these genes might be involved in tissue-specific fat deposition in chickens. Previous studies have suggested that transcription factors play important roles in the regulation of adipocyte differentiation [36,37,38,39]. Transcription factors, such as Kruppel-like factors (KLFs) (KLF9, KLF6 and KLF15), CCAAT enhancer-binding protein alpha (CEBPA), forkhead box O3 (FOXO3), myogenin (MYOG), sterol regulatory element-binding transcription factor 1 (SREBF1), nuclear receptor subfamily 3 group C member 2 (NR3C2), GATA-binding protein 2 (GATA2) and myogenic differentiation 1 (MYOD1), were also differentially expressed between preadipocytes and adipocytes in the different adipose tissues. KLFs and GATAs are transiently induced to control the preadipocyte-to-adipocyte transition [38, 40,41,42]. Furthermore, Jiang et al. suggested that NR3C1/4 may participate in intramuscular adipogenic differentiation by binding to glucocorticoid response elements in the promoters of glucocorticoid-responsive genes to activate their transcription and by regulating other transcription factors [43,44,45].
In the present study, functional annotation analysis of the DEGs revealed that these genes play important roles in some lipid metabolism- and adipogenic differentiation-related pathways, such as the PPAR signalling pathway, ECM-receptor interaction and fatty acid metabolism [34, 43, 46]. PPARG, the most adipocyte-specific and adipogenic member of the PPAR family, is mainly expressed in adipose tissue and plays an important role in lipid metabolism and adipocyte differentiation [47, 48]. The extracellular matrix (ECM) plays an important role in the regulation of proliferation, adipogenic differentiation and migration of preadipocytes . Adipocyte differentiation can also be affected by fatty acid metabolism via the regulation of transcription factors [34, 49].
We noticed that APOA1, CD36, LAMA2, CHAD, MMP1, RRM2, VWF and PIK3R5 were the most upregulated genes during the adipogenic differentiation processes between the IMF and AbF groups. Therefore, these genes may be involved in the positively regulating position specificity fat deposition. In addition, the differences in secretory functions and hormone sensitivities between IMF and AbF might be caused by the position-specific regulation of adipose tissue in poultry. Further studies are necessary for to elucidate the associated cell microenvironments and paracrine signaling pathways and the effects of muscle-specific regulation on adipose tissue in poultry.
In conclusion, our current study showed that abdominal fat (AbF) preadipocytes accumulate more lipids than intramuscular fat (IMF) preadipocytes. This study presents the first analysis of gene expression during the differentiation of intramuscular and abdominal preadipocytes in chickens. A total of 2039 DEGs were identified by a pairwise comparison of preadipocytes at different stages of differentiation. The DEGs were found to be involved in the PPAR signalling pathway, fatty acid biosynthesis, ECM-receptor interaction and focal adhesion, consistent with previous reports on preadipocyte differentiation in chickens. These DEGs and pathways might play significant roles in intramuscular preadipocyte differentiation in chickens. Our findings provide a solid foundation for future studies on the molecular mechanisms underlying tissue-specific fat deposition and on strategies for the improvement of meat quality in poultry.
All animal experiments were performed according to the guidelines of Henan Agricultural University (Institutional Animal Care and Use Committee (IACUC), Permit No. 11–0085, Date: 06–2011). All efforts were made to decrease animal suffering. All birds were euthanized by intraperitoneal injection of pentobarbital (Sigma, St. Louis, MO, USA) (1.0 mg/mL in methanol) at a dose of 45 mg/kg of body weight. The Gushi chickens were provided by the Animal Center of Henan Agricultural University. All birds were raised in the same environmental conditions with ad libitum water and food.
Primary preadipocyte isolation and culture in vitro
Primary preadipocytes were isolated from the breast muscles and abdominal fat of two-week-old Gushi chickens according our previously described method . In brief, breast muscle and abdominal adipose tissues were separated from the body of each chicken under sterile conditions. The tissues were washed using phosphate-buffered saline (PBS) supplemented with penicillin (100 units/mL) and streptomycin (100 μg/mL). The washed tissue was cut into 1-mm3 pieces and then digested with collagenase type II (1 mg/mL, Solarbio, Beijing, China) at 37 °C for 90 min. The digested cell suspension was filtered using 200- and 500-mesh screens to separate the stromal-vascular fraction from undigested tissue and mature adipocytes, and the fraction was then centrifuged at 1000 x g for 5 min. Preadipocytes were plated onto a 6-well culture plate at a density of 1 × 105 cells/mL and maintained in Dulbecco’s modified Eagle’s medium/Ham’s nutrient mixture F-12 (DMEM/F12) supplemented with 10% FBS (Gibco, Beijing, China) with penicillin (100 units/mL) and streptomycin (100 μg/mL) in a humidified atmosphere with 5% (v/v) CO2 at 37 °C. For intramuscular preadipocytes, the differential adherence method was used to separate them from other cells. The basal medium was replaced with fresh medium after 2 hours.
Preadipocyte adipogenic differentiation assay
Preadipocytes were cultured in 6-well plates until they reached 80–90% confluence. For adipogenic differentiation, upon reaching confluence (0d), the cells were exposed to differentiation medium consisting of basal medium (DMEM/F12 with 10% FBS) supplemented with 50 nM insulin, 1 μM dexamethasone, 0.5 mM 3-isobutyl-1-methylxanthine (DMI) and 300 μM oleate (dissolved in DMSO). The cells were collected at day 0, 2, 4, 6, 8, and 10 after induction. Each stage included three biological replicates (n = 3, three wells of cells were collected for every 2 days). The cell samples were stored at − 80 °C until use.
RNA isolation, library preparation and sequencing
Total RNA was extracted using TRIzol reagent (TaKaRa, Dalian, China). RNA was quantified with a Qubit 2.0 (Thermo Fisher Scientific, Waltham, MA, USA) and a Nanodrop ND-2000 spectrophotometer (Thermo Fisher Scientific). Qualified total RNA was further purified with an RNA Clean XP Kit (Beckman Coulter, Inc., Kraemer Boulevard, Brea, CA, USA) and an RNase-Free DNase Set (QIAGEN, GmBH, Germany). RNA purity was assessed using an Agilent Bioanalyzer 2100 (Agilent Technologies, Santa Clara, CA, USA) with a threshold RNA integrity number > 8. The total RNA was stored at − 80 °C until use. After RNA samples were selected for library construction and deep sequencing, the rRNA was removed, and the mRNA was enriched using magnetic beads with oligo (dT) primers. RNA sequencing libraries were generated with a VAHTS™ Total RNA-seq (H/M/R) Library Prep Kit for Illumina® following the manufacturer’s instructions. The quality of all libraries was confirmed using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). The libraries were then analysed by using one lane of a 150 + 150-nt paired-end Illumina HiSeq 2500 run with Illumina sequencing primers. The quality of the raw data was examined using FastQC.
Read mapping and transcriptome assembly
The splice-mapping algorithm of HISAT2 (2.0.4)  was used to perform genome mapping of the pre-processed reads. The clean data were mapped to the Gallus gallus reference genome (GGA5) with HISAT2 (2.0.4), and the default parameters were used.
Analysis of differential gene expression patterns
To determine the differentially expressed genes, Cuffdiff  was used to calculate the expected number of Fragments Per Kilobase of exon model per Million mapped reads (FPKM) for each gene. StringTie (1.3.0) [53, 54] was used for quantitative analysis of the transcripts to obtain the count numbers and FPKM values of the transcripts in each sample. EdgeR  was used for differential gene or mRNAs analysis between samples. For biological replicates, transcripts or genes with a Qvalue (adjusted p value) < 0.05 and a |fold change| ≥ 2 were defined as differentially expressed genes (DEGs) or mRNAs (DEMs) between the two groups.
Functional enrichment analysis of the differentially expressed genes (DEGs)
Gene Ontology (GO) enrichment analysis was performed with the DAVID database (https://david.ncifcrf.gov/). For Kyoto Encyclopedia of Genes and Genomes (KEGG) analysis, the differentially expressed genes or transcripts were analyzed on the KEGG online website (http://www.genome.jp/kegg/). The gene-pathway interaction networks for the DEGs were visualized with Cytoscape 3.4.0 (http://www.cytoscape.org/) .
Oil red O staining
An Oil Red O staining assay was performed according to the methods in our previous study . Briefly, adipocytes were gently washed with cold PBS three times and fixed with 4% paraformaldehyde for 30 min. Then, the fixed cells were gently washed with cold PBS three times, incubated with 60% filtered Oil Red O solution for 40 min and then observed under a phase-contrast microscope to check for positive cells appearing red. The cells were then washed three times with deionized cold PBS and photographed using an Olympus CKX41-F32FL microscope (Olympus, Tokyo, Japan). Subsequently, Oil Red O was eluted from the stained cells with 100% isopropanol (v/v) and quantified with a microplate reader (Thermo Fisher Scientific) at 500 nm. ImageJ software (National Institutes of Health, Bethesda, MD) was used to estimate the adipogenic differentiation ability. The cells with lipid droplets were regarded as differentiated cells: Adipogenic differentiation ratio = (differentiated cells count / total cells count) × 100%.
RNA isolation and real-time quantitative PCR (qRT-PCR)
Primers for the DEGs were designed using Primer3Plus online software (http://www.primer3plus.com/cgi-bin/dev/primer3plus.cgi) (Additional file 6: Table S6). qPCR was performed using SYBR® Green PCR Master Mix (TaKaRa, Dalian, China). The PCR mixture contained 5 μL of SYBR® Premix Ex Taq II (2×), 0.5 μL of forward primer (10 μM), 0.5 μL of reverse primer (10 μM), and 200 ng of cDNA, with RNA-free water added to 10 μL. The qRT-PCR was conducted in a LightCycler 96 system (Roche) by the SYBR Green method. The program included an initial step of 95 °C for 5 min; 40 cycles of 95 °C for 30 s, 60 °C for 30 s, and 72 °C for 20 s; and melting curve generation, which was performed as follows: 95 °C for 10 s, annealing at 65 °C for 20 s, and heating through a continuous temperature gradient to 97 °C with 5 acquisitions/s. Chicken GAPDH was selected as an internal control gene. All samples were examined in triplicate. All data were analysed using the 2-ΔΔCt method. All data are shown as fold changes in gene expression compared with the gene expression in the 0d group.
Statistical significance between two experimental groups was evaluated with a T-test for comparisons in SPSS 20.0 statistical software (IBM, Chicago, IL, USA). Statistical significance among three or more experimental groups was evaluated by one-way ANOVA followed by Dunnett’s test for multiple comparisons in SPSS 22.0. Graphics were drawn using GraphPad Prism 7 (GraphPad Software, San Diego, CA, USA) and RStudio 1.1.453 software. All data are expressed as the mean ± standard error (SE). A P-value ≤0.05 was considered statistically significant, and a P-value ≤0.01 was considered extreme significant.
Availability of data and materials
The full data sets have been submitted to NCBI with bioproject asscession number: PRJNA429489 and Sequence Read Archive (SRA) accession number: SRR6459505, SRR6459506, SRR6459503, SRR6459504, SRR6459509, SRR6459510, SRR6459507 and SRR6459508.
The datasets supporting the conclusions of this article are included within the article and its additional files.
Bone morphogenetic protein 4
CCAAT enhancer binding protein alpha
CCAAT/enhancer-binding protein alpha
Collagen type IV alpha 1 chain
Collagen type VI alpha 1 chain
Differentially expressed genes
Differentially expressed mRNAs
Differentially expressed novel genes
Dubecco’s modified Eagle medium
Extracellular fatty acid-binding protein
Fatty acid binding protein 4
Fatty acid synthetase
Forkhead box O3
GATA binding protein 2
Glycerol-3-phosphate dehydrogenase 1-like 2
Insulin induced gene 1/2
Kyoto Encyclopedia of Genes and Genomes
- Kruppel-like factors:
(KLFs) (KLF9, KLF6 and KLF15)
Laminin, alpha 2
Matrix metallopeptidase 1
Myogenic differentiation 1
Nuclear receptor subfamily 3 group C member 2
Pyruvate dehydrogenase kinase 4
Peroxisome proliferator-activated receptors
Prostaglandin D2 synthase
Quantitative real time PCR
Sterol regulatory element binding transcription factor 1
Sterol regulatory element-binding protein
Nishimura T, ., Hattori A, ., Takahashi K.: Structural changes in intramuscular connective tissue during the fattening of Japanese black cattle: effect of marbling on beef tenderization. J Anim Sci. 1999, 77(1):93–104.
Chartrin P, ., Méteau K, ., Juin H, ., Bernadet MD, Guy G, ., Larzul C, ., Rémignon H, ., Mourot J, ., Duclos MJ, Baéza E: Effects of intramuscular fat levels on sensory characteristics of duck breast meat. J Poult Sci 2006, 85(5):914–922.
Fernandez X, ., Monin G, ., Talmant A, ., Mourot J, ., Lebret B: Influence of intramuscular fat content on the quality of pig meat - 2. Consumer acceptability of m longissimus lumborum. J Meat Sci 1999, 53(1):67–72.
Dransfield E, Sosnicki AA. Relationship between muscle growth and poultry meat quality. J Poult Sci. 1999;78(5):743–6.
Ismail I, Joo S. Poultry Meat Quality in Relation to Muscle Growth and Muscle Fiber Characteristics. Korean J Food Sci Anim Resour. 2017;37(6):873–83.
Leng L, Zhang H, Dong JQ, Wang ZP, Zhang XY, Wang SZ, Cao ZP, Li YM, Li H. Selection against abdominal fat percentage may increase intramuscular fat content in broilers. J Anim Breed Genet. 2016;133(5):422–8.
Emmerson DA. Commercial approaches to genetic selection for growth and feed conversion in domestic poultry. Poult Sci. 1997;76(8):1121.
Jiang M, Fan WL, Xing SY, Wang J, Li P, Liu RR, Li QH, Zheng MQ, Cui HX, Wen JJ. Effects of balanced selection for intramuscular fat and abdominal fat percentage and estimates of genetic parameters. Poult Sci. 2017;96(2):282–7.
Qiu F, Xie L, Ma JE, Luo W, Zhang L, Chao Z, Chen S, Nie Q, Lin Z, Zhang XJ. Lower Expression of SLC27A1 Enhances Intramuscular Fat Deposition in Chicken via Down-Regulated Fatty Acid Oxidation Mediated by CPT1A. Front Physiol. 2017;8:449.
Sun WX, Wang HH, Jiang BC, Zhao YY, Xie ZR, Xiong K, Chen J. Global comparison of gene expression between subcutaneous and intramuscular adipose tissue of mature Erhualian pig. Genet Mol Res. 2013;12(4):5085–101.
Zhou G, Wang S, Wang Z, Zhu X, Shu G, Liao W, Yu K, Gao P, Xi Q, Wang XJ. Global comparison of gene expression profiles between intramuscular and subcutaneous adipocytes of neonatal landrace pig using microarray. Meat Sci. 2010;86(2):440–50.
Cui H, Zheng M, Zhao G, Liu R, Wen J. Identification of differentially expressed genes and pathways for intramuscular fat metabolism between breast and thigh tissues of chickens. BMC Genomics. 2018;19(1):55.
Smith SB, Crouse JD. Relative contributions of acetate, lactate and glucose to lipogenesis in bovine intramuscular and subcutaneous adipose tissue. J Nutr. 1984;114(4):792–800.
Chu W, Wei W, Han H, Gao Y, Liu K, Tian Y, Jiang Z, Zhang L, Chen J. Muscle-specific downregulation of GR levels inhibits adipogenesis in porcine intramuscular adipocyte tissue. Sci Rep. 2017;7(1):510.
Songbo W, Guixuan Z, Gang S, Lina W, Xiaotong Z, Ping G, Qianyun X, Yongliang Z, Li Y, Qingyan J. Glucose utilization, lipid metabolism and BMP-Smad signaling pathway of porcine intramuscular preadipocytes compared with subcutaneous preadipocytes. Cell Physiol Biochem. 2013;31(6):981–96.
Han H, Wei W, Chu W, Liu K, Tian Y, Jiang Z, Chen J. Muscle Conditional Medium Reduces Intramuscular Adipocyte Differentiation and Lipid Accumulation through Regulating Insulin Signaling. Int J Mol Sci. 2017;18(8):1799.
Hood RL, Allen CE. Cellularity of bovine adipose tissue. J Lipid Res. 1973;14(6):605–10.
Hood RL, Allen CE. Lipogenesis in isolated intramuscular adipose tissue from four bovine muscles. J Anim Sci. 1978;46(6):1626.
Lu L, Cui H, Zheng M, Zhao G, Jie WJ. Comparative analysis of differentially expressed genes related to triglyceride metabolism between intramuscular fat and abdominal fat in broilers. Br Poult Sci. 2018;59(5):514–20. https://doi.org/10.1080/00071668.2018.1483573.
Hrdinka C, ., Zollitsch W, Knaus W, ., Lettner F: Effects of dietary fatty acid pattern on melting point and composition of adipose tissues and intramuscular fat of broiler carcasses. J Poult Sci 1996, 75(2):208–215.
Zhou J, Shi SY, Wang JH, Pan L, Ke-Zong QI JCJoVS: Effects of conjugated linoleic acid on different adipose tissue deposits and immune function of yellow broiler chickens. 2008.
Cristancho AG, Lazar MA. Forming functional fat: a growing understanding of adipocyte differentiation. Nat Rev Mol Cell Biol. 2011;12(11):722.
Chen JL, Zhao GP, Zheng MQ, Wen J, Yang N. Estimation of genetic parameters for contents of intramuscular fat and inosine-5′-monophosphate and carcass traits in Chinese Beijing-You chickens. Poult Sci. 2008;87(6):1098.
Gardan D, Gondret F, Louveau I. Lipid metabolism and secretory function of porcine intramuscular adipocytes compared with subcutaneous and perirenal adipocytes. Am J Physiol Endocrinol Metab. 2006;291(2):E372–80.
Gondret F, ., Guitton N, ., Guillerm-Regost C, ., Louveau I, . Regional differences in porcine adipocytes isolated from skeletal muscle and adipose tissues as identified by a proteomic approach. J Anim Sci 2008, 86(9):2115–2125.
Hua ZG, Xiong LJ, Yan C, Qing ZY, Hui GP, Tian YJ, Xin ZR. Comparison of the adipogenesis in intramuscular and subcutaneous adipocytes from Bamei and landrace pigs. Biochemistry and cell biology = Biochimie et biologie cellulaire. 2014;92(4):259–67.
Maryline K, Pierre S. A review of the factors influencing the development of intermuscular adipose tissue in the growing pig. Meat Sci. 2011;88(2):213–20.
Kouba M, ., Bonneau M, . Compared development of intermuscular and subcutaneous fat in carcass and primal cuts of growing pigs from 30 to 140kg body weight. Meat Sci 2009, 81(1):270–274.
Sun X, ., Li M, ., Sun Y, ., Cai H, ., Li R, ., Wei X, ., Lan X, ., Huang Y, ., Lei C, ., Chen H, . The developmental transcriptome landscape of bovine skeletal muscle defined by Ribo-zero ribonucleic acid sequencing. J Anim Sci 2015, 93(12):5648–5658.
Zhao W, He X, Hoadley KA, Parker JS, Hayes DN, Perou CM. Comparison of RNA-Seq by poly (a) capture, ribosomal RNA depletion, and DNA microarray for expression profiling. BMC Genomics. 2014;15(1):1–11.
Lefterova M, Lazar M. New developments in adipogenesis. Trends Endocrinol Metab. 2009;20(3):107–14.
Rosen ED, Hsu CH, Wang X, Sakai S, Freeman MW, Gonzalez FJ, Spiegelman BM. C/EBPα induces adipogenesis through PPARγ: a unified pathway. Genes Dev. 2002;16(1):22.
Frédéric P, Martin K, Namjin C, Acharawan TN, Thanaset S, Rita MDO, Mark L, Mcburney MW, Leonard G. Sirt1 promotes fat mobilization in white adipocytes by repressing PPAR-gamma. Nature. 2004;429(6993):771–6.
Kim JB, Spiegelman BM. ADD1/SREBP1 promotes adipocyte differentiation and gene expression linked to fatty acid metabolism. Genes Dev. 1996;10(9):1096–107.
Cowherd RM, Lyle RE, Jr MG. Molecular regulation of adipocyte differentiation ☆. Semin Cell Dev Biol. 1999;10(1):3.
Nakae J, Kitamura T, Kitamura Y, Arden KC, Accili DJ. The forkhead transcription factor Foxo1 regulates adipocyte differentiation. Dev Cell. 2003;4(1):119–29.
Tontonoz P, Kim JB, Graves RA, Spiegelman BMJM, Biology C. ADD1: a novel helix-loop-helix transcription factor associated with adipocyte determination and differentiation. Mol Cell Biol. 1993;13(8):4753–9.
Tong Q, Dalgin G, Xu H, Ting CN, Leiden JM, Hotamisligil GS. Function of GATA Transcription Factors in Preadipocyte-Adipocyte Transition. Science. 2000;290(5489):134–8.
Tong Q, Tsai J, Tan G, Dalgin G, Hotamisligil GS. Interaction between GATA and the C/EBP family of transcription factors is critical in GATA-mediated suppression of adipocyte differentiation. Mol Cell Biol. 2005;25(2):706–15.
Jiang S, Wei H, Song T, Yang Y, Zhang F, Zhou Y, Peng J, Jiang SJ. KLF13 promotes porcine adipocyte differentiation through PPARγ activation. Cell Biosci. 2015;5(1):28.
Yamamoto K, Sakaguchi M, Medina RJ, Niida A, Sakaguchi Y, Miyazaki M, Kataoka K, Huh NH. Transcriptional regulation of a brown adipocyte-specific gene, UCP1, by KLF11 and KLF15. Biochem Biophys Res Commun. 2010;400(1):175–80.
Pei H, ., Yao Y, ., Yang Y, ., Liao K, ., J-R WJ: Krüppel-like factor KLF9 regulates PPARγ transactivation at the middle stage of adipogenesis. Cell Death Differ. 2011, 18(2):315–327.
Shuzhong J, Hongkui W, Tongxing S, Yang Y, Jian P, Siwen J. Transcriptome comparison between porcine subcutaneous and intramuscular stromal vascular cells during adipogenic differentiation. PLoS One. 2013;8(10):e77094.
Mi-Jeong L, Da-Wei G, Burkey BF, Fried SK. Pathways regulated by glucocorticoids in omental and subcutaneous human adipose tissues: a microarray study. Am J Physiol Endocrinol Metab. 2011;300(3):E571.
Chi-Yi Y, Oleg M, Lee JV, Joanna T, Charlie H, Speed TP, Jen-Chywan WJ. Genome-wide analysis of glucocorticoid receptor binding regions in adipocytes reveal gene network involved in triglyceride homeostasis. PLoS One. 2010;5(12):e15188.
Cai H, Li M, Sun X, Plath M, Li C, Lan X, Lei C, Huang Y, Bai Y, Qi X, et al. Global Transcriptome Analysis During Adipogenic Differentiation and Involvement of Transthyretin Gene in Adipogenesis in Cattle. Front Genet. 2018;9:463.
Michael L, Lazar MA. The many faces of PPARgamma. Cell. 2005;123(6):993–9.
Wang Y, ., Mu Y, ., Li H, ., Ding N, ., Wang Q, ., Wang Y, ., Wang S, ., Wang N, . Peroxisome proliferator-activated receptor-gamma gene: a key regulator of adipocyte differentiation in chickens. Poult Sci 2008, 87(2):226.
Schaffer JE, Lodish HF. Expression cloning and characterization of a novel adipocyte long chain fatty acid transport protein. Cell. 1994;79(3):427–36.
Zhang M, Li DH, Li F, Sun JW, Jiang RR, Li ZJ, Han RL, Li GX, Liu XJ, Kang XT. Integrated analysis of MiRNA and genes associated with meat quality reveals that Gga-MiR-140-5p affects intramuscular fat deposition in chickens. Cell Physiol Biochem. 2018;46(6):2421.
Kim D, Langmead B, Salzberg SL. HISAT: a fast spliced aligner with low memory requirements. Nat Methods. 2015;12(4):357–60.
Trapnell C, Hendrickson DG, Sauvageau M, Goff L, Rinn JL, Pachter L. Differential analysis of gene regulation at transcript resolution with RNA-seq. Nat Biotechnol. 2013;31(1):46–53.
Pertea M, Pertea GM, Antonescu CM, Chang TC, Mendell JT, Salzberg SL. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads. Nat Biotechnol. 2015;33(3):290–5.
Pertea M, Kim D, Pertea GM, Leek JT, Salzberg SL. Transcript-level expression analysis of RNA-seq experiments with HISAT, StringTie and Ballgown. Nat Protoc. 2016;11(9):1650–67.
Robinson MD, DJ MC, Smyth GK. edgeR: a Bioconductor package for differential expression analysis of digital gene expression data. Bioinformatics (Oxford, England). 2010;26(1):139–40.
Shannon P, Markiel A, Ozier O, Baliga NS, Wang JT, Ramage D, Amin N, Schwikowski B, Ideker T. Cytoscape: a software environment for integrated models of biomolecular interaction networks. Genome Res. 2003;13(11):2498–504.
We are grateful to the editor and reviewers for their comments and suggestions. The authors would like to thank Shanghai SHBIO Biotech Co., Ltd. for analyzing the data.
This study was supported by Program for Innovation Research Team of Ministry of Education (IRT16R23), Key Science and Technology Research Project of Henan Province (151100110800), National Natural Science Foundation of China (31572356), National Natural Science Foundation of China (31501948), Research on the Foundation and Frontier Technology of Henan Science and Technology Department (162300410162), the China Agriculture Research System (CARS-40-K04). The funding agencies had not involved in the experimental design, analysis and interpretation of the data or writing of the manuscript.
Ethics approval and consent to participate
All animal experiments were carried out in accordance with the Guidelines for Experimental Animals established by the Ministry of Agricultural of China (Beijing, China). This study was approved by the Animal Welfare Committee of College of Animal Science and Veterinary Medicine, Henan Agricultural University (Zhengzhou, Henan, China) (Permit Number:11–0085; Date: 06–2011).
Consent for publication
The authors declare that they have no competing interests.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
Additional file 1: Table S1.
DEGs among four groups
Additional file 2: Table S2.
DEMs among four groups
Additional file 3: Table S3.
DENGs among four groups
Additional file 4: Table S4.
GO analysis associated with DEGs among four
Additional file 5: Table S5.
KEGG Pathways associated with DEGs among four
Additional file 6: Table S6.
Primers used for qRT-PCR in this study
Additional file 7: Figure S1.
The correlation analysis between RNA-Seq data and qRT-PCR results
Rights and permissions
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
About this article
Cite this article
Zhang, M., Li, F., Ma, Xf. et al. Identification of differentially expressed genes and pathways between intramuscular and abdominal fat-derived preadipocyte differentiation of chickens in vitro. BMC Genomics 20, 743 (2019). https://doi.org/10.1186/s12864-019-6116-0
- Abdominal fat
- Intramuscular fat
- Adipocyte differentiation