Transcriptome profile of liver at different physiological stages reveals potential mode for lipid metabolism in laying hens
- Hong Li†1,
- Taian Wang†1,
- Chunlin Xu1,
- Dandan Wang1,
- Junxiao Ren1,
- Yanmin Li1,
- Yadong Tian1, 2,
- Yanbin Wang1, 2,
- Yuping Jiao1, 4,
- Xiangtao Kang1, 2, 3Email author and
- Xiaojun Liu1, 2, 3Email author
© Li et al. 2015
Received: 27 May 2015
Accepted: 21 September 2015
Published: 9 October 2015
Liver is an important metabolic organ that plays a critical role in lipid synthesis, degradation, and transport; however, the molecular regulatory mechanisms of lipid metabolism remain unclear in chicken. In this study, RNA-Seq technology was used to investigate differences in expression profiles of hepatic lipid metabolism-related genes and associated pathways between juvenile and laying hens. The study aimed to broaden the understanding of liver lipid metabolism in chicken, and thereby to help improve laying performance in the poultry industry.
RNA-Seq analysis was carried out on total RNA harvested from the liver of juvenile (n = 3) and laying (n = 3) hens. Compared with juvenile hens, 2567 differentially expressed genes (1082 up-regulated and 1485 down-regulated) with P ≤ 0.05 were obtained in laying hens, and 960 of these genes were significantly differentially expressed (SDE) at a false discovery rate (FDR) of ≤0.05 and fold-change ≥2 or ≤0.5. In addition, most of the 198 SDE novel genes (91 up-regulated and 107 down-regulated) were discovered highly expressed, and 332 SDE isoforms were identified. Gene ontology (GO) enrichment and KEGG (Kyoto Encyclopedia of Genes and Genomes) pathway analysis showed that the SDE genes were most enrichment in steroid biosynthesis, PPAR signaling pathway, biosynthesis of unsaturated fatty acids, glycerophospholipid metabolism, three amino acid pathways, and pyruvate metabolism (P ≤ 0.05). The top significantly enriched GO terms among the SDE genes included lipid biosynthesis, cholesterol and sterol metabolic, and oxidation reduction, indicating that principal lipogenesis occurred in the liver of laying hens.
This study suggests that the majority of changes at the transcriptome level in laying hen liver were closely related to fat metabolism. Some of the SDE uncharacterized novel genes and alternative splicing isoforms that were detected might also take part in lipid metabolism, although this needs further investigation. This study provides valuable information about the expression profiles of mRNAs from chicken liver, and in-depth functional investigations of these mRNAs could provide new insights into the molecular networks of lipid metabolism in chicken liver.
Liver is the main metabolic organ where more than 70 % of the de novo synthesis of fatty acids occurs in chicken [1–3]. Especially, the liver plays an important role in lipid synthesis, degradation, and transport processes. During the hen laying cycle, hydrophobic lipids including triacylglycerols, cholesteryl esters, cholestery esters, and free fatty acids are synthesized in the liver and assembled to form egg-yolk precursors such as very-low density lipoprotein (VLDL) and vitellogenin particles. The particles are then secreted into the circulation and transferred to the developing oocyte to meet the requirements for embryo growth and development [4–7]. The molecular regulatory mechanisms of these crucial physiological processes have been investigated extensively since the 1970s, and thus are reasonably well understood [8–11].
It is generally agreed that the physiological processes for lipids synthesis, secretion, and transfer in the liver of laying hens are regulated tightly by estrogen. Estrogen regulates the transcription of target genes containing consensus estrogen response elements through the estrogen receptors (ERs) ERα, ERβ, and G protein-coupled receptor (GPR30) [12–14]. The different receptors play distinct roles in gene regulation . Previous studies have revealed that estrogen physiological functions could be mediated by different receptors in various species [15, 16]. However, the specific nuclear receptor subtype that mediates the production of yolk precursors in chicken liver is unclear .
Although most of the genes and their products involved in hepatic lipid metabolism are highly similar in poultry and mammalian species, the functions of some of these genes and their products are considered to be slightly different in poultry compared with their functions in mammals [4, 18–20]. For instance, a recent study on lysophosphatidylglycerol acyltransferase 1 (LPGAT1) indicated that LPGAT1 may play an important role in lipid synthesis in mice  rather than in poultry. Moreover, it has been suggested that poultry species may have lost some of the genes related to lipid metabolism during the evolutionary process . Therefore, the range of genes and their products involved in hepatic lipid metabolism in laying hen remains to be fully elucidated .
How VLDL particles are assembled and secreted in chicken liver is still not fully understood. In mammals, it has been well documented that microsomal triglyceride transfer protein (MTTP) assists in lipoprotein assembly to form low-density lipoprotein [21, 24–28]. The formation of VLDL particles in avian species is tightly regulated by estrogen, and a previous study demonstrated that the up-regulation of MTTP in liver was not required for increased VLDL assembly during the laying period in chicken . Therefore, understanding the synthesis, formation, and transport of yolk precursors in the liver of laying hens is important.
In recent years, the study of gene regulation and interactions has broadened considerably because of advances in genomics, epigenomics, and bioinformatics, as well as with the development of next generation sequencing. RNA-Seq is a novel gene expression profiling technology based on high-throughput sequencing . Compared with other large-scale gene expression profiling methods, RNA-Seq is superior in detecting mRNA expression in different tissues or at different development stages in a single assay, which can help reveal novel genes and splice variants .
In this study, expression profiles of hepatic lipid metabolism-related genes and associated pathways were investigated between juvenile and laying hens (two different physiological stages) using RNA-Seq technology. Because lipogenesis is known to be highly stimulated in the liver of sexually mature hens and to eliminate genes that may be unrelated to lipid metabolism, liver expression profiles were compared between juvenile hens and laying hens. Bioinformatics tools were used to analyze the major differentially expressed genes and pathways. The present study provides an overview of the genes related to lipid metabolism that play a significant role during embryonic development by synthesizing components of the egg yolk.
Animals and liver tissue samples preparation
All animal experiments were performed in accordance with the protocol approved by the Institutional Animal Care and Use Committee (IACUC) of Henan Agricultural University. The experimental animals used in this study were one strain of the Chinese domestic breed laying hens (Lushi green shell chicken).
All the chickens were raised in cages under the same environment with ad libitum conditions. Six hens were selected randomly from two different physiological stages, juvenile hens and laying hens. The three juvenile hens were slaughtered when they were 20 weeks old (L20), and the three laying hens were slaughtered when they were 30 weeks old (L30). Liver tissue samples were harvested immediately. The collected samples were immediately snap-frozen in liquid nitrogen and stored at −80 °C for further use.
Total RNA was extracted from the chicken liver tissues using TRIzol® reagents following the manufacturer’s manual (Invitrogen, Carlsbad, CA). Degradation and contamination of the total RNA was detected on 1 % agarose gels. The purity of the total RNA was assessed using a NanoPhotometer® spectrophotometer (IMPLEN, CA). The integrity was estimated using a RNA Nano 6000 Assay Kit with the Agilent Bioanalyzer 2100 system (Agilent technologies, Santa Clara, CA). The RNA concentration was checked with a Qubit® RNA Assay Kit in a Qubit® 2.0 flurometer (Life Technologies, CA). The 28S/18S ratio of the qualified RNA ranged from 1.8 to 2.0 and the RNA integrity values ranged from 8.0 to 10.0. RNA samples were stored at −80 °C for further analysis.
RNA-Seq library construction and sequencing
Six mRNA libraries were constructed, one for each of the samples (L20-1, L20-2, L20-3 and L30-1, L30-2, L30-3). A total of 3 μg RNA per sample was prepared for mRNA sequencing using the TruSeq RNA Sample Prep Kit v2 (Illumina) according to the manufacturer’s protocol. Briefly, the mRNA was isolated from the total RNA using oligo (dT) beads with two rounds of oligo-dT purification. Following the rRNA depletion step, the purified RNA was fragmented with the Ribo-Zero rRNA Removal Kits (Epicentre). First-strand cDNA synthesis was performed using the Invitrogen random hexamer primers and Superscript II reverse transcriptase (Invitrogen). The second-strand was synthesized using Invitrogen DNA polymerase 1 (Invitrogen). End repair and poly-adenylation were performed, and the mRNAs were ligated to adapters before PCR amplification. The enriched cDNA templates that were 100 nucleotides (nt) long were purified and used for further analysis. The libraries were qualified using a Qubit® 2.0 Fluorometer (Invitrogen) and Qubit dsDNA HS Assay Kit (Invitrogen). The purity and size of the libraries were checked on an Agilent 2100 Bioanalyzer (Agilent Technologies). The adapter-ligated cDNA fragment libraries were run on an Illumina GAIIx analyzer to complete the cluster generation and primer hybridization. Then the Illumina PE flow cell (v3-HS) carrying clusters were sequenced with paired-end 2 × 100 nt multiplex on an Illumina HiSeq 2500 platform following the manufacturer’s instruction (Illumina).
Transcriptome sequencing data processing and annotation
After the sequencing was completed, image data was outputted and transformed into raw reads and stored with a FASTQ format. The obtained raw reads were cleaned using the FASTX-Toolkit (version: 0.0.13) . Reads with adapter, low quality at 3′ end, containing fuzzy N bases, rRNA, sequences shorter than 20 nt and low quality with Q <20 were removed. The resultant clean reads from each sample library were used for the downstream analyses. The clean reads were mapped to the chicken genome assembly (galGal4), which we downloaded from Ensembl , using the spliced mapping algorithm in TopHat2 (version: 2.0.9) . We used TopHat2 as the mapping tool because it can generate a database of splice junctions based on Ensemble annotations of galGal4 and thus can produce a better mapping result than other non-splice mapping tools.
Transcript identification and alternative splicing analysis
We used the reference annotation-based transcript (RABT) assembly method in Cufflink (version: 2.1.1)  to construct and identify both known and novel transcripts from the TopHat2 alignment results. The AStalavist software (version: 3.2) [36, 37] can characterize alternative splicing (AS) for whole transcriptome data from reference annotated transcripts. We used AStalavist to estimate AS events within and between groups. The differentially expressed isoforms were estimated by Cufflink.
Quantification of differential mRNA expression levels
The Cuffdiff was used to analyze the differential expression genes. In our study, the false discovery rate (FDR) was used to determine the threshold of the P-value in multiple tests and analyses. Genes were identified as differentially expressed (DE) genes when P ≤ 0.05. DE genes with fold changes ≥2 or ≤0.5 (FDR ≤0.05) were identified as significantly differentially expressed (SDE) genes .
Quantitative real time PCR (qRT-PCR)
To confirm the repeatability and accuracy of the RNA-Seq gene expression data obtained from the chicken liver libraries, qRT-PCR was carried out on 12 randomly selected DE genes that were prepared from the total RNA. The PrimeScript™ RT Reagent kit with gDNA Eraser (TaKaRa, Dalian, China) was used to synthesize the first-strand cDNA. The qRT-PCRs were performed on a LightCycler® 96 Real-Time PCR system (Roche Applied Science) in a 20-μl reaction volume containing 2 μL cDNA, 10 μL 2 × SYBR®Premix Ex Taq™ II (TliRNaseH Plus) (TaKaRa), 0.5 μL each of forward and reverse primers (10 μM), and 7 μL deionized water. The β-actin gene was used as the reference gene, and all the qRT-PCR gene-specific primers were designed using the Oligo 6.0 software . The primer sequences are presented in Additional file 1: Table S1. The qPCR amplification procedure was as follows: 95 °C for 3 min, 40 cycles of 95 °C for 12 s, 61 °C for 40 s, 72 °C for 30 s, and an extension for 10 min at 72 °C. All the reactions were run with three replicates, and the relative gene expression levels were analyzed using the comparative C T method (also referred to as the 2-△△CT method) . In this study, the Wilcox rank sum test was used. The statistical analyses were performed with R for windows version 3.2.0 , with the test conducted as a one-sided tail test and a significance level of P ≤ 0.05. The values are presented as mean ± standard error.
Functional annotation analyses
Identification of expressed transcripts in the chicken liver transcriptome
Characteristics of the reads from six chicken liver libraries
Raw bases (Gb)
Q20 value (%)
GC content (%)
Mapped unique readsb
Mapping ratio (%)c
The correlation of transcript expression levels between samples is a crucial indicator for the reliability of the experimental results and the rationality of sampling. Generally, the Pearson correlation coefficient shall be no less than 0.92 (r2 ≥ 0.92) . We performed correlation analyses among the six samples to determine whether differential gene expression was observed between the L20 and L30 groups. The Pearson correlation coefficient demonstrated that the expression levels of the three biological replicates in each group (L20 and L30) were similar based on the normalized FPKM values (i.e., all r2 ≥ 0.93).
Identification of differentially expressed genes and isoforms between the two physiological development stages
In this study, we identified a total of 13,532 genes in the chicken liver libraries; 1767 of them were novel genes and 198 of these novel genes showed significant changes in expression (91 up-regulated and 107 down-regulated) between the L20 and L30 groups (FDR ≤0.05) (Additional file 2: Table S2). Among the annotated genes, we identified 2567 DE genes (1082 up-regulated and 1485 down-regulated) in L30 compared with L20 with P ≤ 0.05; 960 of these were SDE genes (473 up-regulated and 487 down-regulated) with a fold-change ≥2 or ≤0.5 (FDR ≤0.05) (Additional file 3: Table S3).
Real-time PCR validation of differential genes expression
Expression patterns of the 12 mRNAs selected for qRT-PCR validation
Functional analysis of differentially expressed genes
To better understand the regulation network of lipid synthesis and transport during egg production, we analyzed the functional distribution of the DE genes in the liver of laying hens liver compared with the liver of juvenile hens.
Top gene ontology clusters of SDE genes between liver samples from juvenile and laying hens
Sterol metabolic process
APOB, APOA1, HMGCR, CYP7A1, LCAT, INSIG1, FDPS, LSS, SC4MOL, VLDLR, DHCR24
Cholesterol metabolic process
APOB, APOA1, HMGCR, CYP7A1, LCAT, INSIG1, FDPS, LSS, VLDLR, DHCR24
FAR1, APOB, APOA1, LCAT, INSIG1, ACSBG2, AACS, VLDLR
Steroid metabolic process
OSBPL3, HMGCR, FDPS, LSS, SC4MOL, APOB, APOA1, LCAT, CYP7A1, INSIG1, OSBPL10, VLDLR, DHCR24
APOB, APOA1, LCAT, INSIG1, VLDLR
APOB, TPRXL, APOA1, PPARG, ANXA1, ATP11A, CETP, VTG2, ATP8B3, VLDLR
Sterol biosynthetic process
HMGCR, FDPS, LSS, SC4MOL, DHCR24
APOB, APOA1, LCAT, INSIG1, VLDLR
APOB, APOA1, CETP, APOVLDLII, VLDLR
Plasma lipoprotein particle
APOB, APOA1, CETP, APOVLDLII, VLDLR
APOB, TPRXL, APOA1, PPARG, ANXA1, ATP11A, CETP, VTG2, ATP8B3, VLDLR
Cholesterol biosynthetic process
HMGCR, FDPS, LSS, DHCR24
APOB, APOVLDLII, VLDLR
Very-low-density lipoprotein particle
APOB, APOVLDLII, VLDLR
Triglyceride-rich lipoprotein particle
APOB, APOVLDLII, VLDLR
Steroid biosynthetic process
HMGCR, FDPS, LSS, SC4MOL, DHCR24
APOB, APOA1, VLDLR
KEGG pathways associated with SDE genes between liver samples from juvenile and laying hens
SOAT1, CYP51A1, DHCR7, LSS, HSD17B7, SC4MOL, NSDHL, DHCR24, FDFT1
PPAR signaling pathway
SCD, PPARG, FADS2, DBI, APOA1, ACSL1, CYP7A1, FABP3, FABP1, FABP5, ACSL5, ACAA1, ANGPTL4
Arginine and proline metabolism
ALDH7A1, ASL2, P4HA2, GATM, ASS1, GLUD1, MAOB, GAMT, AGMAT, PRODH
Biosynthesis of unsaturated fatty acids
PECR, SCD, FADS1, ELOVL2, FADS2, ELOVL6, ACAA1
Glycine, serine and threonine metabolism
GATM, BHMT, MAOB, GCAT, GAMT, PSAT1, GLDC
GPD2, GPD1, PLA2G4A, DGKQ, PLA2G12A, GNPAT, ETNK2, PISD, GPAM, AGPAT2
ME1, LDHB, ALDH7A1, AKR1B1, PDHA1, ACSS2, PDHX
Alanine, aspartate and glutamate metabolism
ASL2, ASS1, GLUD1, ABAT, GPT2, PPAT
KYNU, ALDH7A1, CYP1A1, MAOB, ACMSD, HAAO, AFMID
Glyoxylate and dicarboxylate metabolism
HAO1, PGP, ACO1, AFMID
GALNT2, GALNT1, GALNT6, WBSCR17, GALNTL1, ST6GALNAC1
Valine, leucine and isoleucine degradation
BCAT1, ALDH6A1, ALDH7A1, BCKDHB, ABAT, HIBADH, ACAA1
ABCG8, ABCG5, TAP2, ABCB6, ABCA3, ABCG1
Lipid synthesis and transfer is a dynamic and complex process, and previous studies have suggested that the enzymes involved in this process could play different roles in mammal and chicken or other avian species . Indeed, recent studies have shown that some mammalian genes related to lipid metabolism have been lost in chicken . As a consequence, the regulation in gene expression of lipid metabolism in chicken liver is yet to be fully understood.
DE genes are considered to be important regulatory factors of lipid synthesis and transport in liver during the laying stage of chicken. In this study, we obtained a total of 2567 DE genes between juvenile hens and laying hens livers using RNA-Seq technology. Some of these may participate in lipid biological synthesis, assemble, and transfer at the two different physiological stages. For example, SCD-1 (Stearoyl-CoA desaturase) together with FADS2 (previously named Δ6 desaturase) were up-regulated in the lipogenesis of the PPAR signaling pathway in the liver of laying hen in this study. SCD, which is regulated by a hormone, is a rate limiting enzyme of monounsaturated fatty acid synthesis in liver, and the mRNA expression and activity of SCD-1 have been shown to be triggered by insulin to promote fat synthesis . FADS2, which catalyzes the initial desaturation step to synthesize the long chain polyunsaturated fatty acid (LC-PUFAs), was found to occur mainly in the liver of laying hens , suggesting that FADS2 may contribute to yolk formation. In the liver of severe negative energy balance cows, the down-regulation of FADS2 was shown to suppress the synthesis of LC-PUFAs, arachidonic acid, and eicosapentaenoic acid . Taken together, these results indicated that the RNA-Seq data generated by this study was sufficiently representative of the chicken liver transcriptome.
VLDL assembly occurs in the endoplasmic reticulum of hepatic cells. To assemble a lipoprotein particle that is competent for transport through the secretory pathway, ApoB has to interact with triglycerides, cholesteryl esters, free cholesterol, and phospholipids. This is a highly regulated process that requires the activity of MTTP [54, 55]. It was reported that chicken MTTP contained functionally important domains that are commonly found in the large lipid transport protein family . However, another study showed that MTTP did not respond to the increase of VLDL induced by estrogen either in vivo or in vitro . In this study, MTTP was not significantly differentially expressed in the liver of laying hens compared to juvenile hens. This study is consistent with a previous finding in the coordinated up-regulation of protein components, such as ApoB (8-fold change) and ApoVLDL-II (320-fold change), along with the up-regulation of lipid synthesis led to increased production of VLDL during an egg-laying cycle in avian [57, 58]. The female-specific yolk precursor proteins VTG1 and VTG2, which are synthesized in the liver and depend on estrogen stimulation , were found to be abundantly expressed in the liver of laying hens, and the expression of VTG2 was the main subtype . Estrogen works via the estrogen receptors (ERα and ERβ) and GPR30 that regulate the transcription of target genes, which contain estrogen response elements. Previous research has suggested that in vitro, ERα rather than ERβ and GPR30 could mediate estrogen’s effects on stimulating vitellogenin and ApoVLDL production, while in vivo, ERβ was up-regulated in liver of laying hen in comparison to pullet . Our finding is consistent with previous report that the expression of ERβ but not ERα was significantly up-regulated in vivo. Since so many genes are directly or indirectly involved in the complex physiological process in liver during egg-laying stage, it is hard to believe that estrogen mediates these genes only via ERβ. Therefore, these results may still need to be further validated in future studies.
FABPs have been reported to have multiple biological functions, including roles in hepatic fatty acid oxidation [61, 62], intracellular fatty acid transport , storage, and export, as well as in cholesterol and phospholipid metabolism [64–66]. In this study, FABP1 and FABP3 were both significantly up-regulated in the liver of laying hens compared with juvenile hens, which suggests that they may promote lipid metabolism in the PPAR signaling pathway to meet the requirements of laying eggs. Acyl-CoA binding protein (ACBP, also known as DBI) was reported to act as an endogenous modulator to regulate the levels of gonadal hormones in vivo . The transcriptional factor sterol regulatory element binding protein (SREBP-1) and fatty acid synthase (FASN) genes were both found to be elevated coordinately in laying chicken liver that could synthesize fatty acids de novo , which was consistent with a previous report . Peroxisomal proliferator-activated receptor α (PPARα, a transcriptional factor) controls the expression of fatty acid oxidative metabolism by modulating the expression of peroxisomal acyl-CoA oxidase and mitochondrial carnitine palmitoyltransferase , and it has been reported to be highly expressed in rodent liver  and swine adipose tissue . However, in the present study, PPARα was suppressed in the laying hen liver, suggesting that adipose tissue may oxidize sizeable quantities of fatty acids in avian species, and perhaps also in other mammalian species.
LPGAT1 belongs to a large group of acyltransferases and is a member of the lysophosphatidic acid acyltransferase family. LPGAT1 promotes hepatic lipogenesis in mice  and also may be involved in triacylglycerol synthesis and secretion in liver . However, in this study we found that LPGAT1 was down-regulated in laying hens liver. In addition, in a related study we showed that down-regulated LPGAT1 was induced by estrogen both in vivo and in vitro (data not shown). All these results suggested that LPGAT1 may have different expression patterns in mammals and avian related to specific functions in regulating fatty acid synthesis. Furthermore, LPGAT1 may have multiple subcellular localizations, and could therefore potentially have multiple functions in different cells or within the same cells .
In mammals, lipogenesis is known to occur in liver, adipose tissue, and mammary gland, whereas, in avian species, it occurs mainly in avian liver . During the egg laying stage, fat synthesis in chicken liver is especially active . The GO annotation cluster analyses (Additional file 5: Table S5) showed that the SDE genes were involved mainly in lipid biosynthesis, transport and localization, sterol and cholesterol metabolism, as well as in immune response and some other processes. In poultry, the ovary cannot synthesize lipids; therefore, liver lipoproteins are transferred in the plasma and deposited into the oocytes to form the egg yolk in laying hens. Therefore, lipid synthesis in chicken liver and lipoprotein transfer plays a crucial role on the egg production performance of hens. Some of the SDE genes that are not be involved in lipid metabolism may instead contribute to liver homeostasis in response to the dramatic increase in lipogenesis and protein biosynthesis in the liver of hens at the laying stage.
Alternative splicing of pre-mRNA plays an important role in regulating gene expression in higher eukaryotes. A previous report indicated that 40–60 % of human genes have alternative splicing isoforms, although some variants exist only in relatively low abundance . It has been shown that proteins with different functions can be produced by a diverse array of mRNAs derived from a single pre-mRNA, suggesting that alternative splicing is a crucial mechanism for regulating life . The three alternative splicing isoforms α, δ, and γ of the PPAR gene were detected in our transcriptomic data, isoform γ was significantly down-regulated in L30 compared with L20, while isoforms α and δ were not differentially expressed. It has been shown that PPARα and PPARγ may play significant roles in glucose and lipid metabolism in the early life stage of mouse . Moreover, the DE novel genes detected in this study may provide important information about liver lipid metabolism in chicken. For example, a significantly up-regulated novel gene ENSGALG00000014190 with four alternative splicing isoforms was observed in our transcriptome sequencing data. This gene was predicated to encode a protein of 357 amino acids that could take part in the lipid metabolic process (UniProt: F1NXW6), which requires confirmation. Another up-regulated novel gene ENSGALG00000023444 with three isoforms was also observed, but its isoforms and function need to be investigated further.
Lipid metabolism is controlled by multiple pathways and influenced by multiple genes. These pathways include the PPAR signaling pathway, steroid biosynthesis, steroid hormone biosynthesis, and biosynthesis of unsaturated fatty acids . In our KEGG analysis, the PPAR signaling pathway, which is essential for lipid metabolism, showed one of the most significant associations with the SDE genes in the livers of laying hens. Eighteen DE genes involved in the PPAR pathway (Additional file 6: Figure S1); 11 were up-regulated and seven were down-regulated. In the PPAR pathway, a cytochrome P450 (CYP7A1) catalyzes the rate limiting step of conversion of cholesterol into bile acids. CYP7A1 is also involved in the KEGG Bile secretion pathway, and was reported to be up-regulated in severe negative energy balance cows . The altered expression patterns of hepatic genes in the PPAR signaling pathway could play a role in regulating the lipid metabolism. In addition, a total nine DE genes (Additional file 7: Figure S2) which all were SDE ones were found to be involved in the steroid biosynthesis pathway and all of them were up-regulated except sterol O-acyltransferase (SOAT1, esterification to fatty acids), which suggests that this pathway was quite active in steroid hormone synthesis. SOATs (SOAT1 and SOAT2) are known to synthesize cholesterol fatty acid esters using fatty acids released from membrane phospholipids .
During the laying stage, gene expression is highly stimulated in liver to support the metabolic changes associated with the development of the reproductive organs. In the present study, we identified 960 SDE genes with a fold change ≥2 or ≤0.5 (FDR ≤0.05) in the livers of laying hens compared with juvenile hens. Although species-specific differences should be considered when comparing chicken with mammalian systems, the current findings appear to be consistent with conservation of lipid metabolism and adipogenesis processes in chicken and mammal. The chicken liver transcriptome reported here could greatly broaden our understanding of the regulation and networks of gene expression related to liver lipid metabolism in hens at different physiological stages. Our results will serve as important resource for revealing the mechanism of lipid metabolism during egg-laying stage.
This study generated transcriptomic data using RNA-Seq technology that will help to expand our understanding of the molecular repertoire of lipid metabolism-related genes at different physiological stages in chicken. Differences in expressed genes were found between the juvenile and egg laying stages, including highly expressed novel genes, splice isoforms, and pathways. These findings will be a valuable resource for biological investigations of liver lipid metabolism-related genes in chicken, and may also provide clues for understanding the molecular mechanisms in other poultry and mammalian species.
The raw sequencing data, mapped data, and data for visualization of the RNA-Seq analyses of the chicken liver transcriptome data at different physiological stages have been deposited in the Gene Expression Omnibus (GEO) at the National Center for Biotechnology Information (NCBI) under accession number GSE70010.
This research was supported by the Earmarked Fund for Modern Agro-Industry Technology Research System (no. CARS-41-K04), the Program for Innovation Research Team of Ministry of Education (IRT1236), and the International Cooperation Project of Henan, China (144300510044).
Open AccessThis 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.
- O’hea E, Leveille G. Lipogenesis in isolated adipose tissue of the domestic chick (Gallus domesticus). Comp Biochem Physiol. 1968;26(1):111–20.View ArticlePubMedGoogle Scholar
- Leveille GA, O’Hea EK, Chakrabarty K. In vivo lipogenesis in the domestic chicken. Exp Biol Med. 1968;128(2):398–401.View ArticleGoogle Scholar
- Wang X, Yang L, Wang H, Shao F, Yu J, Jiang H, et al. Growth Hormone-Regulated mRNAs and miRNAs in Chicken Hepatocytes. PLoS One. 2014;9(11):e112896.PubMed CentralView ArticlePubMedGoogle Scholar
- Wiskocil R, Bensky P, Dower W, Goldberger RF, Gordon JI, Deeley RG. Coordinate regulation of two estrogen-dependent genes in avian liver. Proc Natl Acad Sci U S A. 1980;77(8):4474–8.PubMed CentralView ArticlePubMedGoogle Scholar
- Walzem RL, Hansen RJ, Williams DL, Hamilton RL. Estrogen induction of VLDLy assembly in egg-laying hens. J Nutr. 1999;129(2S Suppl):467S–72S.PubMedGoogle Scholar
- Brady L, Romsos DR, Leveille GA. In vivo estimation of fatty acid synthesis in the chicken (Gallus domesticus) utilizing 3H2O. Comp Biochem Physiol B. 1976;54(3):403–7.PubMedGoogle Scholar
- Nikolay B, Plieschnig JA, Subik D, Schneider JD, Schneider WJ, Hermann M. A novel estrogen-regulated avian apolipoprotein. Biochimie. 2013;95(12):2445–53.PubMed CentralView ArticlePubMedGoogle Scholar
- Schneider WJ. Yolk precursor transport in the laying hen. Curr Opin Lipidol. 1995;6(2):92–6.View ArticlePubMedGoogle Scholar
- Bujo H, Hermann M, Lindstedt KA, Nimpf J, Schneider WJ. Low density lipoprotein receptor gene family members mediate yolk deposition. J Nutr. 1997;127 Suppl 5:801S–4S.PubMedGoogle Scholar
- Schneider WJ. Low density lipoprotein receptor relatives in chicken ovarian follicle and oocyte development. Cytogenet Genome Res. 2007;117(1–4):248–55.View ArticlePubMedGoogle Scholar
- Schneider WJ. Receptor-mediated mechanisms in ovarian follicle and oocyte development. Gen Comp Endocrinol. 2009;163(1–2):18–23.View ArticlePubMedGoogle Scholar
- Kumar V, Chambon P. The estrogen receptor binds tightly to its responsive element as a ligand-induced homodimer. Cell. 1988;55(1):145–56.View ArticlePubMedGoogle Scholar
- Paech K, Webb P, Kuiper GG, Nilsson S, Gustafsson J, Kushner PJ, et al. Differential ligand activation of estrogen receptors ERalpha and ERbeta at AP1 sites. Science. 1997;277(5331):1508–10.View ArticlePubMedGoogle Scholar
- Pache JC, Burton DW, Deftos LJ, Hastings RH. A carboxyl leucine-rich region of parathyroid hormone-related protein is critical for nuclear export. Endocrinology. 2006;147(2):990–8.View ArticlePubMedGoogle Scholar
- Menuet A, Anglade I, Flouriot G, Pakdel F, Kah O. Tissue-specific expression of two structurally different estrogen receptor alpha isoforms along the female reproductive axis of an oviparous species, the rainbow trout. Biol Reprod. 2001;65(5):1548–57.View ArticlePubMedGoogle Scholar
- Leanos-Castaneda O, Van Der Kraak G. Functional characterization of estrogen receptor subtypes, ERalpha and ERbeta, mediating vitellogenin production in the liver of rainbow trout. Toxicol Appl Pharmacol. 2007;224(2):116–25.View ArticlePubMedGoogle Scholar
- Li J, Leghari IH, He B, Zeng W, Mi Y, Zhang C. Estrogen stimulates expression of chicken hepatic vitellogenin II and very low-density apolipoprotein II through ER-alpha. Theriogenology. 2014;82(3):517–24.View ArticlePubMedGoogle Scholar
- Kirchgessner TG, Heinzmann C, Svenson KL, Gordon DA, Nicosia M, Lebherz HG, et al. Regulation of chicken apolipoprotein B: cloning, tissue distribution, and estrogen induction of mRNA. Gene. 1987;59(2–3):241–51.View ArticlePubMedGoogle Scholar
- Hermier D, Catheline D, Legrand P. Relationship between hepatic fatty acid desaturation and lipid secretion in the estrogenized chicken. Comp Biochem Physiol A Physiol. 1996;115(3):259–64.View ArticlePubMedGoogle Scholar
- Mason TM. The role of factors that regulate the synthesis and secretion of very-low-density lipoprotein by hepatocytes. Crit Rev Clin Lab Sci. 1998;35(6):461–87.View ArticlePubMedGoogle Scholar
- Soh J, Iqbal J, Queiroz J, Fernandez-Hernando C, Hussain MM. MicroRNA-30c reduces hyperlipidemia and atherosclerosis in mice by decreasing lipid synthesis and lipoprotein secretion. Nat Med. 2013;19(7):892–900.PubMed CentralView ArticlePubMedGoogle Scholar
- Dakovic N, Terezol M, Pitel F, Maillard V, Elis S, Leroux S, et al. The loss of adipokine genes in the chicken genome and implications for insulin metabolism. Mol Biol Evol. 2014;31(10):2637–46.View ArticlePubMedGoogle Scholar
- Riegler B, Besenboeck C, Bauer R, Nimpf J, Schneider WJ. Enzymes involved in hepatic acylglycerol metabolism in the chicken. Biochem Biophys Res Commun. 2011;406(2):257–61.View ArticlePubMedGoogle Scholar
- Hussain MM, Rava P, Pan X, Dai K, Dougan SK, Iqbal J, et al. Microsomal triglyceride transfer protein in plasma and cellular lipid metabolism. Curr Opin Lipidol. 2008;19(3):277–84.View ArticlePubMedGoogle Scholar
- Hussain MM, Iqbal J, Anwar K, Rava P, Dai K. Microsomal triglyceride transfer protein. a multifunctional protein. Front Biosci. 2003;8:s500–506.View ArticlePubMedGoogle Scholar
- Hussain MM, Rava P, Walsh M, Rana M, Iqbal J. Multiple functions of microsomal triglyceride transfer protein. Nutr Metab (Lond). 2012;9:14.View ArticleGoogle Scholar
- Berriot-Varoqueaux N, Aggerbeck LP, Samson-Bouma M, Wetterau JR. The role of the microsomal triglygeride transfer protein in abetalipoproteinemia. Annu Rev Nutr. 2000;20:663–97.View ArticlePubMedGoogle Scholar
- Rustaeus S, Lindberg K, Stillemark P, Claesson C, Asp L, Larsson T, et al. Assembly of very low density lipoprotein: a two-step process of apolipoprotein B core lipidation. J Nutr. 1999;129(Suppl 2S):463S–6S.PubMedGoogle Scholar
- Ivessa NE, Rehberg E, Kienzle B, Seif F, Hermann R, Hermann M, et al. Molecular cloning, expression, and hormonal regulation of the chicken microsomal triglyceride transfer protein. Gene. 2013;523(1):1–9.View ArticlePubMedGoogle Scholar
- Wang Z, Gerstein M, Snyder M. RNA-Seq: a revolutionary tool for transcriptomics. Nat Rev Genet. 2009;10(1):57–63.PubMed CentralView ArticlePubMedGoogle Scholar
- Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley DR, et al. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks. Nat Protoc. 2012;7(3):562–78.PubMed CentralView ArticlePubMedGoogle Scholar
- FASTX-Toolkit T: http://hannonlab.cshl.edu/fastx_toolkit/index.html. Accessed 26 Jun 2015.
- Ensemble Galgal4:ftp://ftp.ensembl.org/pub/release-75/fasta/gallus_gallus/dna/Gallus_gallus.Galgal4.75.dna.toplevel.fa.gz. Accessed 26 Jun 2015.
- Kim D, Pertea G, Trapnell C, Pimentel H, Kelley R, Salzberg SL. TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions. Genome Biol. 2013;14(4):R36.PubMed CentralView ArticlePubMedGoogle Scholar
- Kong L, Zhang Y, Ye ZQ, Liu XQ, Zhao SQ, Wei L, et al. CPC: assess the protein-coding potential of transcripts using sequence features and support vector machine. Nucleic Acids Res. 2007;35(Web Server issue):W345–349.PubMed CentralView ArticlePubMedGoogle Scholar
- Foissac S, Sammeth M. ASTALAVISTA: dynamic and flexible analysis of alternative splicing events in custom gene datasets. Nucleic Acids Res. 2007;35(Web Server issue):W297–299.PubMed CentralView ArticlePubMedGoogle Scholar
- Sammeth M. Complete alternative splicing events are bubbles in splicing graphs. J Comput Biol. 2009;16(8):1117–40.View ArticlePubMedGoogle Scholar
- Anders S, Pyl PT, Huber W. HTSeq—a Python framework to work with high-throughput sequencing data. Bioinformatics. 2015;31(2):166–9.PubMed CentralView ArticlePubMedGoogle Scholar
- Chen C, Deng B, Qiao M, Zheng R, Chai J, Ding Y, et al. Solexa sequencing identification of conserved and novel microRNAs in backfat of Large White and Chinese Meishan pigs. PLoS One. 2012;7(2):e31426.PubMed CentralView ArticlePubMedGoogle Scholar
- Rchlik W. OLIGO 7 Primer Analysis Software. Methods Mol Biol. 2007;402:35–60.View ArticleGoogle Scholar
- Schmittgen TD, Livak KJ. Analyzing real-time PCR data by the comparative CT method. Nat Protoc. 2008;3(6):1101–8.View ArticlePubMedGoogle Scholar
- Gentleman RR. In: Team RDC, editor. A language and environment for statistical computing. Vienna, Austria: R Foundation for Statistical Computing; 2008.Google Scholar
- Dennis Jr G, Sherman BT, Hosack DA, Yang J, Gao W, Lane HC, et al. DAVID: database for annotation, visualization, and integrated discovery. Genome Biol. 2003;4(5):3.View ArticleGoogle Scholar
- Huang DW, Sherman BT, Lempicki RA. Bioinformatics enrichment tools: paths toward the comprehensive functional analysis of large gene lists. Nucleic Acids Res. 2009;37(1):1–13.PubMed CentralView ArticleGoogle Scholar
- Chen R, Zhang F, Song L, Shu Y, Lin Y, Dong L, et al. Transcriptome profiling reveals that the SM22α-regulated molecular pathways contribute to vascular pathology. J Mol Cell Cardiol. 2014;72:263–72.View ArticlePubMedGoogle Scholar
- Han W, Zou J, Wang K, Su Y, Zhu Y, Song C, et al. High-Throughput Sequencing Reveals Hypothalamic MicroRNAs as Novel Partners Involved in Timing the Rapid Development of Chicken (Gallus gallus) Gonads. PLoS One. 2015;10(6):e0129738.PubMed CentralView ArticlePubMedGoogle Scholar
- Kang X, Liu G, Liu Y, Xu Q, Zhang M, Fang M. Transcriptome profile at different physiological stages reveals potential mode for curly fleece in Chinese tan sheep. PLoS One. 2013;8(8):e71763.PubMed CentralView ArticlePubMedGoogle Scholar
- Lander ES, Linton LM, Birren B, Nusbaum C, Zody MC, Baldwin J, et al. Initial sequencing and analysis of the human genome. Nature. 2001;409(6822):860–921.View ArticlePubMedGoogle Scholar
- Johnson JM, Castle J, Garrett-Engele P, Kan Z, Loerch PM, Armour CD, et al. Genome-wide survey of human alternative pre-mRNA splicing with exon junction microarrays. Science. 2003;302(5653):2141–4.View ArticlePubMedGoogle Scholar
- Wang ET, Sandberg R, Luo S, Khrebtukova I, Zhang L, Mayr C, et al. Alternative isoform regulation in human tissue transcriptomes. Nature. 2008;456(7221):470–6.PubMed CentralView ArticlePubMedGoogle Scholar
- Lefevre P, Diot C, Legrand P, Douaire M. Hormonal regulation of stearoyl coenzyme-A desaturase 1 activity and gene expression in primary cultures of chicken hepatocytes. Arch Biochem Biophys. 1999;368(2):329–37.View ArticlePubMedGoogle Scholar
- Jacobi SK, Lin X, Corl BA, Hess HA, Harrell RJ, Odle J. Dietary arachidonate differentially alters desaturase-elongase pathway flux and gene expression in liver and intestine of suckling pigs. J Nutr. 2011;141(4):548–53.View ArticlePubMedGoogle Scholar
- McCabe M, Waters S, Morris D, Kenny D, Lynn D, Creevey C. RNA-seq analysis of differential gene expression in liver from lactating dairy cows divergent in negative energy balance. BMC Genomics. 2012;13(1):193.PubMed CentralView ArticlePubMedGoogle Scholar
- Gordon DA, Jamil H. Progress towards understanding the role of microsomal triglyceride transfer protein in apolipoprotein-B lipoprotein assembly. Biochim Biophys Acta. 2000;1486(1):72–83.View ArticlePubMedGoogle Scholar
- Gordon DA, Wetterau JR, Gregg RE. Microsomal triglyceride transfer protein: a protein complex required for the assembly of lipoprotein particles. Trends Cell Biol. 1995;5(8):317–21.View ArticlePubMedGoogle Scholar
- Rava P, Hussain MM. Acquisition of triacylglycerol transfer activity by microsomal triglyceride transfer protein during evolution. Biochemistry. 2007;46(43):12263–74.PubMed CentralView ArticlePubMedGoogle Scholar
- Capony F, Williams DL. Apolipoprotein B of avian very low density lipoprotein: characteristics of its regulation in nonstimulated and estrogen-stimulated rooster. Biochemistry. 1980;19(10):2219–26.View ArticlePubMedGoogle Scholar
- Chan L, Jackson R, O’malley B, Means A. Synthesis of very low density lipoproteins in the cockerel. Effects of estrogen. J Clin Invest. 1976;58(2):368–79.PubMed CentralView ArticlePubMedGoogle Scholar
- Deeley RG, Gordon JI, Burns AT, Mullinix KP, Binastein M, Goldberger RF. Primary activation of the vitellogenin gene in the rooster. J Biol Chem. 1977;252(22):8310–9.PubMedGoogle Scholar
- Evans MI, Silva R, Burch J. Isolation of chicken vitellogenin I and III cDNAs and the developmental regulation of five estrogen-responsive genes in the embryonic liver. Genes Dev. 1988;2(1):116–24.View ArticlePubMedGoogle Scholar
- Veerkamp J. Fatty acid-binding protein and its relation to fatty acid oxidation. Mol Cell Biol. 1993;123:101–6.Google Scholar
- Kaikaus RM, Sui Z, Lysenko N, Wu NY, de Montellano PO, Ockner R, et al. Regulation of pathways of extramitochondrial fatty acid oxidation and liver fatty acid-binding protein by long-chain monocarboxylic fatty acids in hepatocytes. Effect of inhibition of carnitine palmitoyltransferase I. J Biol Chem. 1993;268(36):26866–71.PubMedGoogle Scholar
- Woodford J, Behnke W, Schroeder F. Liver fatty acid binding protein enhances sterol transfer by membrane interaction. Mol Cell Biochem. 1995;152(1):51–62.PubMedGoogle Scholar
- Jefferson J, Powell D, Rymaszewski Z, Kukowska-Latallo J, Lowe J, Schroeder F. Altered membrane structure in transfected mouse L-cell fibroblasts expressing rat liver fatty acid-binding protein. J Biol Chem. 1990;265(19):11062–8.PubMedGoogle Scholar
- Martin GG, Danneberg H, Kumar LS, Atshaves BP, Erol E, Bader M, et al. Decreased liver fatty acid binding capacity and altered liver lipid distribution in mice lacking the liver fatty acid-binding protein gene. J Biol Chem. 2003;278(24):21429–38.View ArticlePubMedGoogle Scholar
- Atshaves BP, McIntosh AM, Lyuksyutova OI, Zipfel W, Webb WW, Schroeder F. Liver fatty acid-binding protein gene ablation inhibits branched-chain fatty acid metabolism in cultured primary hepatocytes. J Biol Chem. 2004;279(30):30954–65.View ArticlePubMedGoogle Scholar
- Dong E, Matsumoto K, Watanabe H. Diazepam binding inhibitor (DBI) reduces testosterone and estradiol levels in vivo. Life Sci. 2002;70(11):1317–23.View ArticlePubMedGoogle Scholar
- Gondret F, Ferré P, Dugail I. ADD-1/SREBP-1 is a major determinant of tissue differential lipogenic capacity in mammalian and avian species. J Lipid Res. 2001;42(1):106–13.PubMedGoogle Scholar
- Hillgartner FB, Salati LM, Goodridge AG. Physiological and molecular mechanisms involved in nutritional regulation of fatty acid synthesis. Physiol Rev. 1995;75(1):47–76.PubMedGoogle Scholar
- Ferré P. The biology of peroxisome proliferator-activated receptors relationship with lipid metabolism and insulin sensitivity. Diabetes. 2004;53 suppl 1:S43–50.View ArticlePubMedGoogle Scholar
- Bergen WG, Mersmann HJ. Comparative aspects of lipid metabolism: impact on contemporary research and use of animal models. J Nutr. 2005;135(11):2499–502.PubMedGoogle Scholar
- Ding S, Schinckel A, Weber T, Mersmann H. Expression of porcine transcription factors and genes related to fatty acid metabolism in different tissues and genetic populations. J Anim Sci. 2000;78(8):2127–34.PubMedGoogle Scholar
- Hiramine Y, Emoto H, Takasuga S, Hiramatsu R. Novel acyl-coenzyme A: monoacylglycerol acyltransferase plays an important role in hepatic triacylglycerol secretion. J Lipid Res. 2010;51(6):1424–31.PubMed CentralView ArticlePubMedGoogle Scholar
- Traurig MT, Orczewska JI, Ortiz DJ, Bian L, Marinelarena AM, Kobes S, et al. Evidence for a role of LPGAT1 in influencing BMI and percent body fat in Native Americans. Obesity. 2013;21(1):193–202.PubMed CentralView ArticlePubMedGoogle Scholar
- Klasing K. Nutritional modulation of resistance to infectious diseases. Poult Sci. 1998;77(8):1119–25.View ArticlePubMedGoogle Scholar
- Wang N-Y, Lai H-C, Au L-C. Methods for enrichment of a mRNA isoform with specific alternative splicing. Anal Biochem. 2006;359(2):189–93.View ArticlePubMedGoogle Scholar
- Schwerk C, Schulze-Osthoff K. Regulation of apoptosis by alternative pre-mRNA splicing. Mol Cell. 2005;19(1):1–13.View ArticlePubMedGoogle Scholar
- Zheng J, Xiao X, Zhang Q, Yu M, Xu J, Wang Z. Maternal protein restriction induces early-onset glucose intolerance and alters hepatic genes expression in the peroxisome proliferator-activated receptor pathway in offspring. J Diabetes Investig. 2015;6(3):269–79.PubMed CentralView ArticlePubMedGoogle Scholar
- Fatima A, Waters S, O’Boyle P, Seoighe C, Morris DG. Alterations in hepatic miRNA expression during negative energy balance in postpartum dairy cattle. BMC Genomics. 2014;15(1):28.PubMed CentralView ArticlePubMedGoogle Scholar
- Cases S, Novak S, Zheng Y-W, Myers HM, Lear SR, Sande E, et al. ACAT-2, a second mammalian acyl-CoA: cholesterol acyltransferase its cloning, expression, and characterization. J Biol Chem. 1998;273(41):26755–64.View ArticlePubMedGoogle Scholar