不同体重鹅肝脏转录组的功能分析

doi:10.11843/j.issn.0366-6964.2020.05.010

摘要/Abstract

摘要： 旨在探讨肝调节营养物质吸收转化影响鹅平均日增重的分子机制。本研究依据平均日增重将70日龄的四川白鹅进行分组（高：H；中：M；低：L），然后采用RNA-Seq技术分析各组肝RNA表达，筛选出差异表达基因，并对其进行功能富集及聚类关联分析。结果显示，H、M和L组的平均日增重差异显著（P<0.05），分别为（52.24±0.41）、（44.44±0.42）和（37.80±0.47）g·d^-1。从测序数据分别鉴定出14 632、14 651和14 819个表达基因，其中有13 776个基因共表达。Top 100的基因表达量占所有基因表达量的65.38%~70.37%，其中，有85个共表达基因，17个特异性表达基因。鉴定出435个差异表达基因，主要富集于细胞凋亡、蛋白重折叠、胰岛素分泌、营养反应、细胞增长、氨基酸转运、cAMP信号通路等生物学过程。聚类构建的5个差异表达基因模块与平均日增重和脂质代谢指标显著相关（|r|≥ 0.75，P<0.05）。总体而言，本研究鉴定了体重相关的潜在候选基因和代谢通路，初步探讨了肝调节营养物质吸收转化影响鹅体重的分子机制，以使人们更深入地了解肝代谢如何影响鹅的体重及有助于设计新的选择策略来改善鹅的生产。

关键词: 鹅, 体重, 肝脏, 转录组, 功能

Abstract: This study aimed to explore the molecular mechanism by which the liver regulates the absorption and conversion of nutrients and affects the body weight of goose. RNA-Seq technology was used to analyze the RNA expression in livers of Sichuan White goose with different average daily gain (high:H; middle:M; low:L) at 70-days-old. The differentially expressed genes were screened, their functional enrichment and cluster correlation analysis were performed. The results showed that the average daily gain of the H, M, and L groups were significantly different (P<0.05), which were (52.24±0.41), (44.44±0.42), and (37.80±0.47) g·d^-1, respectively. There were 14 632, 14 651 and 14 819 expression genes identified, of which 13 776 genes were co-expressed. The expression level of the top 100 genes accounted for 65.38%-70.37% of all gene expression, including 85 co-expression genes and 17 specific expression genes. Four hundred thirty five differentially expressed genes were identified, and mainly enriched in apoptosis, protein refolding, insulin secretion, nutritional response, cell growth, amino acid transport, and cAMP signaling pathway. The 5 gene modules constructed from differentially expressed genes were significantly associated with average daily gain and lipid metabolism index (|r|≥ 0.75, P<0.05). Overall, body weight related potential candidate genes and metabolic pathways were identified. The molecular mechanism by which the liver regulates nutrient absorption and conversion and affects the weight of goose was explored, which make us better understand how liver metabolism affects the weight of goose and can help us to design new selection strategies to improve goose production.

Key words: goose, body weight, liver, transcriptome, function

中图分类号:

S835.2

章杰, 熊子标, 李龙娇, 田旭, 揭晓蝶, 程雅婷, 万堃, 刘安芳, 向邦全, 何航. 不同体重鹅肝脏转录组的功能分析[J]. 畜牧兽医学报, 2020, 51(5): 987-996.

ZHANG Jie, XIONG Zibiao, LI Longjiao, TIAN Xu, JIE Xiaodie, CHENG Yating, WAN Kun, LIU Anfang, XIANG Bangquan, HE Hang. Functional Analysis of Liver Transcriptomes of Goose with Different Body Weight[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(5): 987-996.

参考文献 49

[1]	李慧锋,张臻,朱文进,等.白来航蛋鸡产蛋前后肝转录组功能分析[J].畜牧兽医学报,2017,48(9):1624-1634.LI H F,ZHANG Z,ZHU W J,et al.Functional analysis of liver transcriptome before and after white leghorn hen begin laying eggs[J].Acta Veterinaria et Zootechnica Sinica,2017,48(9):1624-1634.(in Chinese)
[2]	FAM B C,JOANNIDES C N,ANDRIKOPOULOS S.The liver:key in regulating appetite and body weight[J]. Adipocyte, 2012,1(4):259-264.
[3]	MUOZ R,ESTANY J,TOR M,et al.Hepatic lipogenic enzyme expression in pigs is affected by selection for decreased backfat thickness at constant intramuscular fat content[J].Meat Sci,2013,93(3):746-751.
[4]	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.
[5]	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.
[6]	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].J Anim Sci Biotechnol,2014,5(1):32.
[7]	YE M H,ZHOU B,WEI S S,et al.Transcriptomic analysis identifies candidate genes related to intramuscular fat deposition and fatty acid composition in the breast muscle of Squabs (Columba)[J].G3,2016,6(7):2081-2090.
[8]	WANG X,ZHANG Y Q,ZHANG X Z,et al.The comprehensive liver transcriptome of two cattle breeds with different intramuscular fat content[J].Biochem Biophys Res Commun,2017,490(3):1018-1025.
[9]	GUNEWARDENA S S,YOO B,PENG L,et al.Deciphering the developmental dynamics of the mouse liver transcriptome[J]. PLoS One,2015,10(10):e0141220.
[10]	MIRANDA A,LÓPEZ-CARDONA A P,LAGUNA-BARRAZA R,et al.Transcriptome profiling of liver of non-genetic low birth weight and long term health consequences[J].BMC Genomics,2014,15(1):327.
[11]	PINGEL H.Waterfowl production for food security[J].Lohmann Inf,2011,46(2):32-42.
[12]	LU L Z,CHEN Y,WANG Z,et al.The goose genome sequence leads to insights into the evolution of waterfowl and susceptibility to fatty liver[J].Genome Biol,2015,16(1):89.
[13]	NRC.Nutrient requirements of poultry[M].8th ed.Washington,DC:National Academy Press,1994.
[14]	KIM D,LANGMEAD B,SALZBERG S L.HISAT:a fast spliced aligner with low memory requirements[J].Nat Methods,2015,12(4):357-360.
[15]	PERTEA M,KIM D,PERTEA G M,et al.Transcript-level expression analysis of RNA-seq experiments with HISAT,StringTie and Ballgown[J].Nat Protoc,2016,11(9):1650-1667.
[16]	DILLIES M A,RAU A,AUBERT J,et al.A comprehensive evaluation of normalization methods for Illumina high-throughput RNA sequencing data analysis[J].Brief Bioinform,2013,14(6):671-683.
[17]	ROBINSON M D,MCCARTHY D J,SMYTH G K.edgeR:a Bioconductor package for differential expression analysis of digital gene expression data[J].Bioinformatics,2010,26(1):139-140.
[18]	BENJAMINI Y,YEKUTIELI D.The control of the false discovery rate in multiple testing under dependency[J].Ann Stat,2001,29(4):1165-1188.
[19]	DA W H,SHERMAN B T,LEMPICKI R A.Systematic and integrative analysis of large gene lists using DAVID bioinformatics resources[J].Nat Protoc,2009,4(1):44-57.
[20]	WOLFS M G,RENSEN S S,BRUIN-VAN DIJK E J,et al.Co-expressed immune and metabolic genes in visceral and subcutaneous adipose tissue from severely obese individuals are associated with plasma HDL and glucose levels:a microarray study[J].BMC Med Genomics,2010,3(1):34.
[21]	LEE J N,DUTTA R K,MAHARJAN Y,et al.Catalase inhibition induces pexophagy through ROS accumulation[J].Biochem Biophys Res Commun,2018,501(3):696-702.
[22]	ALLALI-HASSANI A,PAN P W,DOMBROVSKI L,et al.Structural and chemical profiling of the human cytosolic sulfotransferases[J]. PLoS Biol,2007,5(5):e97.
[23]	HARDWICK J P.Cytochrome P450 omega hydroxylase (CYP4) function in fatty acid metabolism and metabolic diseases[J]. Biochem Pharmacol,2008,75(12):2263-2275.
[24]	JOSEPH R,POSCHMANN J,SUKARIEH R,et al.ACSL1 is associated with fetal programming of insulin sensitivity and cellular lipid content[J]. Mol Endocrinol,2015,29(6):909-920.
[25]	TAKENOUCHI T,SHIMADA H,UEHARA T,et al.A paradoxical thrombogenic mutation in factor II at the target site of arthropod bleeding toxin[J].Eur J Med Genet,2019,62(2):93-95.
[26]	WINGES A,GARCIA T B,PRAGER P,et al.Osmotic expression of aldose reductase in retinal pigment epithelial cells:involvement of NFAT5[J].Graefes Arch Clin ExpOphthalmol,2016,254(12):2387-2400.
[27]	HILLIAN A D,SEBASTIAN B M,NAGY L E.Complement and non-alcoholic liver disease:Role of C4 and C1q in diet-induced obesity in mice[J].Mol Immunol,2010,47(13):2282.
[28]	FANG B A,KOVAČEVIČ Ž,PARK K C,et al.Molecular functions of the iron-regulated metastasis suppressor,NDRG1,and its potential as a molecular target for cancer therapy[J].Biochim Biophys Acta,2014,1845(1):1-19.
[29]	NAKANO I,DOUGHERTY J D,KIM K,et al.Phosphoserine phosphatase is expressed in the neural stem cell niche and regulates neural stem and progenitor cell proliferation[J].Stem Cells,2007,25(8):1975-1984.
[30]	LU W W,ZHANG H,LI Y M,et al.Gastric cancer-derived heat shock protein-gp96 peptide complex enhances dendritic cell activation[J].World J Gastroenterol,2017,23(24):4390-4398.
[31]	LLANOS S,SERRANO M.Depletion of ribosomal protein L37 occurs in response to DNA damage and activates p53 through the L11/MDM2 pathway[J].Cell Cycle,2010,9(19):4005-4012.
[32]	CHEN J,GUO K X,KASTAN M B.Interactions of nucleolin and ribosomal protein L26(RPL26) in translational control of human p53 mRNA[J].J Biol Chem,2012,287(20):16467-16476.
[33]	KOBAYASHI T,SASAKI Y,OSHIMA Y,et al.Activation of the ribosomal protein L13 gene in human gastrointestinal cancer[J].Int J Mol Med,2006,18(1):161-170.
[34]	HARRIS T A,YAMAKUCHI M,FERLITO M,et al.MicroRNA-126 regulates endothelial expression of vascular cell adhesion molecule[J].Proc Natl Acad Sci U S A,2008,105(5):1516-1521.
[35]	WANG Z Q,MILNE K,WEBB J R,et al.CD74 and intratumoral immune response in breast cancer[J]. Oncotarget, 2016,8(8):12664-12674.
[36]	LI G Y,ZHANG L,LIU J Z,et al.shRNA-mediated RPS15A silencing inhibits U937 acute myeloid leukemia cell proliferation and enhances apoptosis[J].Mol Med Rep,2016,13(5):4400-4406.
[37]	FREUNDT J K,FROMMEYER G,W TZEL F,et al.The transcription factor ATF4 promotes expression of cell stress genes and cardiomyocyte death in a cellular model of atrial fibrillation[J].Biomed Res Int,2018,2018:3694362.
[38]	CUI H X,ZHENG M Q,LIU R R,et al.Liver dominant expression of fatty acid synthase (FAS) gene in two chicken breeds during intramuscular-fat development[J].Mol Biol Rep,2012,39(4):3479-3484.
[39]	MARIONI J C,MASON C E,MANE S M,et al.RNA-seq:an assessment of technical reproducibility and comparison with gene expression arrays[J].Genome Res,2008,18(9):1509-1517.
[40]	GAO G L,LI Q,ZHAO X Z,et al.Transcriptome profiling of the hypothalamus during prelaying and laying periods in Sichuan white geese (Anser cygnoides)[J].Anim Sci J,2015,86(8):800-805.
[41]	NING D,HAN Q,LI Q,et al.Comprehensive analysis of Sichuan white geese (Anser cygnoides) transcriptome[J].Anim Sci J,2014,85(6):650-659.
[42]	XU Q,ZHAO W M,CHEN Y,et al.Transcriptome profiling of the goose (Anser cygnoides) ovaries identify laying and broodiness phenotypes[J].PLoS One,2013,8(2):e55496.
[43]	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.
[44]	YU Y,PING J,CHEN H,et al.A comparative analysis of liver transcriptome suggests divergent liver function among human,mouse and rat[J].Genomics,2010,96(5):281-289.
[45]	MUKHERJEE R,YUN J W.Long chain acyl CoA synthetase 1 and gelsolin are oppositely regulated in adipogenesis and lipogenesis[J].Biochem Biophys Res Commun,2012,420(3):588-593.
[46]	HAN C C,WANG J W,XU H Y,et al.Effect of overfeeding on plasma parameters and mRNA expression of genes associated with hepatic lipogenesis in geese[J].Asian Austral J Anim,2008,21(4):590-595.
[47]	LI X R,ZHANG P Y,WANG B B,et al.NDRG1 negatively regulates proliferation and milk bio-synthesis of bovine epithelial cells via the mTOR signaling pathway[J].Res Vet Sci,2019,124:158-165.
[48]	CAI K,EL-MERAHBI R,LOEFFLER M,et al.Ndrg1 promotes adipocyte differentiation and sustains their function[J].Sci Rep,2017,7(1):7191.
[49]	COSTA C B,LUNARDELLI P A,FONTES P K,et al.Influence of cAMP modulator supplementation of in vitro culture medium on Bos taurus indicus embryos[J].Theriogenology,2020,141:134-141.