分子营养 Molecular Nutrition

不同脂肪源饲粮育成猪肝脏转录组差异分析

  • 李建平 ,
  • 秦贵信 ,
  • 赵志辉 ,
  • 张巧灵 ,
  • 王大力 ,
  • 孙博兴 ,
  • 姜海龙
展开
  • 1. 吉林农业大学动物科学技术学院, 长春 130118;
    2. 河南牧业经济学院畜牧工程系, 郑州 450011;
    3. 吉林大学农业试验基地, 长春 130062;
    4. 吉林大学动物医学学院, 长春 130062;

收稿日期: 2015-01-06

  网络出版日期: 2015-07-07

基金资助

吉林省安全优质肉猪现代生产关键技术示范(20140301013ny);国家科技支撑计划课题生猪健康养殖模式构建与示范(2012bad39b03)

Differentially Analysis of Liver Transcriptomes of Finishing Pigs Fed Different Dietary Fats Sources

  • LI Jianping ,
  • QIN Guixin ,
  • ZHAO Zhihui ,
  • ZHANG Qiaoling ,
  • WANG Dali ,
  • SUN Boxing ,
  • JIANG Hailong
Expand
  • 1. College of Animal Science and Technology, Jilin Agricultural University, Changchun 130118, China;
    2. Department of Animal Husbandry, Henan University of Animal Husbandry and Economy, Zhengzhou 450011, China;
    3. Agriculture Test Station, Jilin University, Changchun 130062, China;
    4. College of Animal Science and Veterinary Medicine, Jilin University, Changchun 130062, China

Received date: 2015-01-06

  Online published: 2015-07-07

摘要

本试验旨在肥育猪饲粮中分别添加红花籽油和椰子油,采集育成猪肝脏组织,进行高通量转录组测序,找出2种处理间的差异表达基因。利用Illumina HiSeqTM 2500高通量RNA-seq测序技术对2种处理育成猪肝脏进行转录组测序,使用TopHat2软件将测序得到的Reads序列与猪参考基因组(Sscrofa 10.2)序列比对,找出差异表达基因,并在Nr、GO和KEGG数据库中进行功能注释、富集分析和聚类分析。结果显示:红花籽油组和椰子油组肝脏差异表达基因共有938个,与椰子油组相比,红花籽油组表达上调基因有479个,下调基因有459个;GO功能分类注释到细胞组成、生物学过程和分子功能数据库中差异表达基因数分别有773、768和729个;注释到KEGG通路中差异表达基因数346个,显著富集通路为酮体生成与降解通路和类萜骨架生物合成通路(P < 0.05)。

本文引用格式

李建平 , 秦贵信 , 赵志辉 , 张巧灵 , 王大力 , 孙博兴 , 姜海龙 . 不同脂肪源饲粮育成猪肝脏转录组差异分析[J]. 动物营养学报, 2015 , 27(7) : 2128 -2139 . DOI: 10.3969/j.issn.1006-267x.2015.07.018

Abstract

The study was to detect differentially expressed genes in liver of finishing pigs which were divided into two groups and separately fed safflower seed oil and coconut oil. The chosen pigs were slaughtered and liver tissue were isolated and their mRNA were respectively extracted. The liver transcriptomes were analyzed by Illumina HiSeq TM 2500 high-throughput RNA sequencing system. All RNA-Seq reads were mapped on the reference pig genome (Sscrofa10.2) using TopHat2 software. All of differentially expressed genes were annotated using Nr, GO and KEGG databases. The results showed that there were 938 differentially expressed genes in liver between two groups. 479 genes were up-regulated and 459 genes were down-regulated in safflower seed oil group vs. coconut oil group. There were 773 differentially expressed genes annotated into database of cellular component, 768 into biological process, 729 into molecular function in GO databases. There were 346 differentially expressed genes annotated into KEGG databases, and clustered pathways of differentially expressed genes included synthesis and degradation of ketone bodies and terpenoid backbone biosynthesis(P < 0.05).

参考文献

[1] ROSSI R,PASTORELLI G,CANNATA S,et al.Recent advances in the use of fatty acids as supplements in pig diets:a review[J].Animal Feed Science and Technology,2010,162(1/2):1-11.

[2] RAMÍREZ M,AMATE L,GIL A.Absorption and distribution of dietary fatty acids from different sources[J].Early Human Development,2001,65(Suppl.2):S95-S101.

[3] SUN H,JIANG T,WANG S B,et al.The effect of LXRα,ChREBP and Elovl6 in liver and white adipose tissue on medium- and long-chain fatty acid diet-induced insulin resistance[J].Diabetes Research and Clinical Practice,2013,102(3):183-192.  

[4] VALLIM T,SALTER A M.Regulation of hepatic gene expression by saturated fatty acids[J].Prostaglandins,Leukotrienes and Essential FattyAcids,2010,82(4/5/6):211-218.

[5] TAKEUCHI H,MATSUO T,TOKUYAMA K,et al.Diet-induced thermogenesis is lower in rats fed a lard diet than in those fed a high oleic acid safflower oil diet,a safflower oil diet or a linseed oil diet[J].The Journal of Nutrition,1995,125(4):920-925.

[6] SHIMOMURA Y,TAMURA T,SUZUKI M.Less body fat accumulation in rats fed a safflower oil diet than in rats fed a beef tallow diet[J].The Journal of Nutrition,1990,120(11):1291-1296 .

[7] IYER M N H,SARMAH B C,TAMULI M K,et al.Effect of dietary sunflower oil and coconut oil on adipose tissue gene expression,fatty acid composition and serum lipid profile of grower pigs[J].Archives of Animal Nutrition,2012,66(4):271-282.  

[8] JIA L H,LIU Y,LI Y Z.Determination of wholesome elements and heavy metals in safflower (Carthamus tinctorius L.) from Xinjiang and Henan by ICP-MS/ICP-AES[J].Journal of Pharmaceutical Analysis,2011,1(2):100-103.  

[9] ASGARPANAH J,KAZEMIVASH N.Phytochemistry,Pharmacology and medicinal properties of Carthamus tinctorius L.[J].Chinese Journal of Integrative Medicine,2013,19(2):153-159.  

[10] HIGA R,KURTZ M,MAZZUCCO M B,et al.Folic acid and safflower oil supplementation interacts and protects embryos from maternal diabetes-induced damage[J].Molecular Human Reproduction,2012,18(5):253-264.  

[11] APPAIAH P,SUNIL L,PRASANTH KUMAR P K,et al.Composition of coconut testa,coconut kernel and its oil[J].Journal of the American Oil Chemists Society,2014,91(6):917-924.  

[12] DEBMANDAL M,MANDAL S.Coconut (Cocos nucifera L.:Arecaceae):in health promotion and disease prevention[J].Asian Pacific Journal of Tropical Medicine,2011,4(3):241-247.  

[13] KIM D,PERTEA G,TRAPNELL C,et al.TopHat2:accurate alignment of transcriptomes in the presence of insertions,deletions and gene fusions[J].Genome Biology,2013,14(4):R36.

[14] TRAPNELL C,WILLIAMS B A,PERTEA G,et al.Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation[J].Nature Biotechnology,2010,28(5):511-515.  

[15] ALTSCHUL S F,MADDEN T L,SCHÄFFER A A,et al.Gapped BLAST and PSI-BLAST:a new generation of protein database search programs[J].Nucleic Acids Research,1997,25(17):3389-3402.  

[16] DENG Y Y,LI J Q,WU S F,et al.Integrated nr database in protein annotation system and its localization[J].Computer Engineering,2006,32(5):71-73,76.

[17] ASHBURNER M,BALL C A,BLAKE J A,et al.Gene ontology:tool for the unification of biology[J].Nature Genetics,2000,25(1):25-29.  

[18] KANEHISA M,GOTO S,KAWASHIMA S,et al.The KEGG resource for deciphering the genome[J].Nucleic Acids Research,2004,32(Suppl.1):D277-D280.

[19] CHEN C Y,AI H S,REN J,et al.A global view of porcine transcriptome in three tissues from a full-sib pair with extreme phenotypes in growth and fat deposition by paired-end RNA sequencing[J].BMC Genomics,2011,12(1):448.

[20] COROMINAS J,RAMAYO-CALDAS Y,PUIG-OLIVERAS A,et al.Analysis of porcine adipose tissue transcriptome reveals differences in de novo fatty acid synthesis in pigs with divergent muscle fatty acid composition[J].BMC Genomics,2013,14(1):843.

[21] SEVANE N,BIALADE F,VELASCO S,et al.Dietary inulin supplementation modifies significantly the liver transcriptomic profile of broiler chickens[J].PLoS One,2014,9(6):e98942.

[22] LI R W,RINALDI M,CAPUCO A V.Characterization of the abomasal transcriptome for mechanisms of resistance to gastrointestinal nematodes in cattle[J].Veterinary Research,2011,42(1):114.

[23] FAN R W,XIE J S,BAI J M,et al.Skin transcriptome profiles associated with coat color in sheep[J].BMC Genomics,2013,14(1):389.

[24] MORTON N M,NELSON Y B,MICHAILIDOU Z,et al.A stratified transcriptomics analysis of polygenic fat and lean mouse adipose tissues identifies novel candidate obesity genes[J].PLoS One,2011,6(9):e23944

[25] PATEL O V,CASEY T,DOVER H,et al.Homeorhetic adaptation to lactation:comparative transcriptome analysis of mammary,liver,and adipose tissue during the transition from pregnancy to lactation in rats[J].Functional & Integrative Genomics,2011,11(1):193-202.  

[26] XING K,ZHU F,ZHAI L W,et al.The liver transcriptome of two full-sibling Songliao black pigs with extreme differences in backfat thickness[J].Journal of Animal Science and Biotechnology,2014,5(1):32.

[27] ESTEVE-CODINA A,KOFLER R,PALMIERI N,et al.Exploring the gonad transcriptome of two extreme male pigs with RNA-seq[J].BMC Genomics,2011,12(1):552.

[28] RAMAYO-CALDAS Y,MACH N,ESTEVE-CODINA A,et al.Liver transcriptome profile in pigs with extreme phenotypes of intramuscular fatty acid composition[J].BMC Genomics,2012,13(1):547.

[29] CHAKRABORTY S,CAI Y,TARR M A.In vitro oxidative footprinting provides insight into apolipoprotein B-100 structure in low-density lipoprotein[J].Proteomics,2014,14(21/22):2614-2622.

[30] LI H N,YAN J R,OU W J,et al.Construction of a biotinylated cameloid-like antibody for lable-free detection of apolipoprotein B-100[J].Biosensors and Bioelectronics,2015,64:111-118.

[31] SU Q Z,TSAI J,XU E,et al.Apolipoprotein B100 acts as a molecular link between lipid-induced endoplasmic reticulum stress and hepatic insulin resistance[J].Hepatology,2009,50(1):77-84.  

[32] CAVIGLIA J M,GAYET C,OTA T,et al.Different fatty acids inhibit apoB100 secretion by different pathways:unique roles for ER stress,ceramide,and autophagy[J].Journal of Lipid Research,2011,52(9):1636-1651.  

[33] BHASKAR,KUMARI N,GOYAL N.Cloning,characterization and sub-cellular localization of gamma subunit of T-complex protein-1 (chaperonin) from Leishmania donovani[J].Biochemical and Biophysical Research Communications,2012,429(1/2):70-74.

[34] FOLGUEIRA C,REQUENA J M.A postgenomic view of the heat shock proteins in kinetoplastids[J].FEMS Microbiology Reviews,2007,31(4):359-377.  

[35] SHONHAI A,MAIER A G,PRZYBORSKI J M,et al.Intracellular protozoan parasites of humans:the role of molecular chaperones in development and pathogenesis[J].Protein & Peptide Letters,2011,18(2):143-157.  

[36] CYRNE L,GUERREIRO P,CARDOSO A C,et al.The Tetrahymena chaperonin subunit CCTη gene is coexpressed with CCTγ gene during cilia biogenesis and cell sexual reproduction[J].FEBS Letters,1996,383(3):277-283.  

[37] SEIXAS C,CRUTO T,TAVARES A,et al.CCTα and CCTδ chaperonin subunits are essential and required for cilia assembly and maintenance in Tetrahymena[J].PLoS One,2010,5(5):e10704.

[38] GRANTHAM J,RUDDOCK L W,ROOBOL A,et al.Eukaryotic chaperonin containing T-complex polypeptide 1 interacts with filamentous actin and reduces the initial rate of actin polymerization in vitro[J].Cell Stress & Chaperones,2002,7(3):235-242.  

[39] FRYDMAN J,NIMMESGERN E,ERDJUMENT-BROMAGE H,et al.Function in protein folding of TRiC,a cytosolic ring complex containing TCP-1 and structurally related subunits[J].EMBO Journal,1992,11(13):4767-4778.

[40] BOURKE G J,ELALAMI W,WILSON S J,et al.Slow axonal transport of the cytosolic chaperonin CCT with Hsc73 and actin in motor neurons[J].Journal of Neuroscience Research,2002,68(1):29-35.  

[41] DAVIS T A,LOOS B,ENGELBRECHT A-M.AHNAK:the giant jack of all trades[J].Cellular Signalling,2014,26(12):2683-2693.  

[42] HAASE H.Ahnak,a new player in β-adrenergic regulation of the cardiac L-type Ca2+ channel[J].Cardiovascular Research,2007,73(1):19-25.  

[43] MARG A,HAASE H,NEUMANN T,et al.AHNAK1 and AHNAK2 are costameric proteins:AHNAK1 affects transverse skeletal muscle fiber stiffness[J].Biochemical and Biophysical Research Communications,2010,401(1):143-148.  

[44] DE MORRÉE A,DROOG M,GRAND MOURSEL L,et al.Self-regulated alternative splicing at the AHNAK locus[J].The FASEB Journal,2012,26(1):93-103.  

[45] HUANG Y C,LAVAL S H,VAN REMOORTERE A,et al.AHNAK,a novel component of the dysferlin protein complex,redistributes to the cytoplasm with dysferlin during skeletal muscle regeneration[J].The FASEB Journal,2007,21(3):732-742.  

[46] ALVAREZ J L,PETZHOLD D,PANKONIEN I,et al.Ahnak1 modulates L-type Ca2+ channel inactivation of rodent cardiomyocytes[J].Pflügers Archiv-European Journal of Physiology,2010,460(4):719-730.

[47] NEBERT D W,WIKVALL K,MILLER W L.Human cytochromes P450 in health and disease[J].Philosophical Transactions of the Royal Society of London Series B-Biological Sciences,2013,368(1612):20120431.

[48] COOPER D Y,ESTABROOK R W,ROSENTHAL O.The stoichiometry of C21 hydroxylation of steroids by adrenocortical microsomes[J].Journal of Biological Chemistry,1963,238:1320-1323.

[49] ORRENIUS S.Some aspects on the hydroxylation of drugs,steroid hormones and fatty acids (omega-oxidation) in rat liver microsomes[J].Hoppe-Seylers Zeitschrift für Physiologische Chemie,1968,349(11):1619-1621.  

[50] BAILÓN-SOTO C E,GALAVIZ-HERMNDEZ C,LAZALDE-RAMOS B P,et al.Influence of CYP1A1*2C on high triglyceride levels in female Mexican indigenous tarahumaras[J].Archives of Medical Research,2014,45(5):409-416.  

[51] RUSSELL D W.The enzymes,regulation,and genetics of bile acid synthesis[J].Annual Review of Biochemistry,2003,72(1):137-174.  

[52] NORLIN M,WIKVALL K.Enzymes in the conversion of cholesterol into bile acids[J].Current Molecular Medicine,2007,7(2):199-218.  

[53] MAST N,GRAHAM S E,ANDERSSON U,et al.Cholesterol binding to cytochrome P450 7A1,a key enzyme in bile acid biosynthesis[J].Biochemistry,2005,44(9):3259-3271.  

[54] ABRAHAMSSON A,KRAPIVNER S,GUSTAFSSON U,et al.Common polymorphisms in the CYP7A1 gene do not contribute to variation in rates of bile acid synthesis and plasma LDL cholesterol concentration[J].Atherosclerosis,2005,182(1):37-45.  

[55] BARCELOS A L V,CHIES R,ALMEIDA S E M,et al.Association of CYP7A1-278A>C polymorphism and the response of plasma triglyceride after dietary intervention in dyslipidemic patients[J].Brazilian Journal of Medical and Biological Research,2009,42(6):487-493.  

[56] 王镜岩,朱圣庚,徐长法.生物化学[M].3版.北京:高等教育出版社,2002.

[57] GOLDHABER S B.Trace element risk assessment:essentiality vs. toxicity[J].Regulatory Toxicology and Pharmacology,2003,38(2):232-242.  
文章导航

/