简州大耳羊肌内脂肪细胞成脂分化差异表达基因的筛选与鉴定

doi:10.11843/j.issn.0366-6964.2020.07.006

摘要/Abstract

摘要： 旨在通过转录组测序技术获得简州大耳羊肌内前体脂肪细胞成脂分化前后的差异表达基因，并经生物信息学分析获得相关信号通路及可能发挥作用的关键功能候选基因。本研究以7日龄的健康简州大耳羊公羊为试验动物（n=3），采用胶原酶消化法分离获得其肌内前体脂肪细胞；利用Illumina平台对肌内前体脂肪细胞和诱导分化5 d的肌内脂肪细胞cDNA样品（n=3）进行高通量测序；以|log₂fold change|>0和P<0.05为阈值筛选获得差异表达基因，并利用clusterProfiler R包对其进行GO功能和KEGG通路富集；最后利用实时荧光定量PCR（quantitative real-time PCR，qRT-PCR）技术检测功能候选基因在细胞分化前后的表达量变化。结果显示，共获得差异表达基因7 916个，其中4 143个为表达上调基因，主要富集到氧化磷酸化、核糖体合成和三羧酸循环通路等305条通路；3 773个为表达下调基因，且主要富集到甲状腺激素信号通路等303条通路；差异基因GO功能注释中52.8%为生物过程、13.4%是细胞组成和33.8%是分子功能。qRT-PCR结果表明，UCP3、ACACB、ACOT11、ACOX3、APOA1和WISP2在分化前后的山羊肌内脂肪细胞中的表达趋势与RNA-seq结果一致，提示这些基因适合作为下一步研究的功能候选基因。本研究筛选得到简州大耳羊肌内脂肪细胞成脂分化的差异基因，并确认了6个基因可作为功能候选基因。研究结果为阐明肉用山羊肌内脂肪细胞成脂分化的分子调控网络提供基础数据和系统资料。

关键词: 简州大耳羊, 肌内前体脂肪细胞, 转录组测序, 差异表达基因, 诱导分化

Abstract: This experiment aimed to obtain the differentially expressed genes before and after the adipogenic differentiation of Jianzhou Da'er goat intramuscular preadipocytes by the RNA-seq technology, and acquire the related signaling pathways and key candidate functional genes by bioinformatics analysis. The 7-day-old healthy male Jianzhou Da'er goat was used (n=3) and their intramuscular preadipocytes were isolated by collagenase digestion. Illumina platform was used to carry out the high-throughput sequencing to two groups of cDNA samples (n=3) of Jianzhou Da'er goat intramuscular preadipocytes and intramuscular adipocytes induced for 5 day, respectively. The differentially expressed genes were obtained according to the thresholds of|log₂fold change|>0 and P<0.05. The differentially expressed genes were enriched to GO terms and KEGG pathways using the clusterProfiler R. Ultimately, the qRT-PCR technology was performed to detect the change of expression levels of candidate functional genes before and after differentiation of preadipocytes. The results showed that there were 7 916 differentially expressed genes in the RNA-seq results during the adipogenic differentiation of Jianzhou Da'er goat intramuscular preadipocytes, including 4 143 up-regulated genes that were enriched to 305 KEGG pathways, including oxidative phosphorylation, ribosome synthesis, and citrate cycle pathways, etc., and 3 773 down-regulated genes that were enriched to 303 KEGG pathways, including thyroxine signaling pathways, and so on; In the results of GO function enrichment analysis, the biological process was 52.8% of GO terms, the cellular component was 13.4%, and the molecular function was 33.8%. The qRT-PCR results showed that the trends of UCP3, ACACB, ACOT11, ACOX3, APOA1, and WISP2 expression during the adipogenic differentiation of intramuscular preadipocytes were consistent with the results of RNA-seq, which indicated that these genes were the suitable candidate functional genes for further research. This study obtained the differentially expressed genes during the adipogenic differentiation of Jianzhou Da'er goat intramuscular adipocytes and verified 6 candidate functional genes. The results of this study provided the basic data and systematic information for elucidating the molecular regulatory networks during the adipogenic differentiation of meat goat intramuscular adipocytes.

Key words: Jianzhou Da'er goat, intramuscular preadipocytes, RNA-seq technology, differentially expressed genes, induced differentiation

中图分类号:

S827.2

谢光杰, 王永, 许晴, 朱江江, 林亚秋. 简州大耳羊肌内脂肪细胞成脂分化差异表达基因的筛选与鉴定[J]. 畜牧兽医学报, 2020, 51(7): 1525-1536.

XIE Guangjie, WANG Yong, XU Qing, ZHU Jiangjiang, LIN Yaqiu. Selection and Validation of the Differentially Expressed Genes during the Adipogenic Differentiation of Jianzhou Da'er Goat Intramuscular Adipocytes[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(7): 1525-1536.

参考文献 43

[1]	俄木曲者,熊朝瑞,范景胜,等.简州大耳羊产肉性能世代选育进展及肉品质研究[J].中国草食动物科学,2014,34(3):9-12.EMU Q Z,XIONG Z R,FAN J S,et al.Jianzhou Da'er goat generation of progress in breeding of sheep meat production performance and meat quality research[J].China Herbivore Science,2014,34(3):9-12.(in Chinese)
[2]	HOCQUETTE J F,GONDRET F,BAÉZA E,et al.Intramuscular fat content in meat-producing animals:development,genetic and nutritional control,and identification of putative markers[J].Animal,2010,4(2):303-319.
[3]	HUNT M R,LEGAKO J F,DINH T T N,et al.Assessment of volatile compounds,neutral and polar lipid fatty acids of four beef muscles from USDA Choice and Select graded carcasses and their relationships with consumer palatability scores and intramuscular fat content[J].Meat Sci,2016,116:91-101.
[4]	GONDRET F,LEBRET B.Feeding intensity and dietary protein level affect adipocyte cellularity and lipogenic capacity of muscle homogenates in growing pigs,without modification of the expression of sterol regulatory element binding protein[J].J Anim Sci,2002,80(12):3184-3193.
[5]	GONDRET F,MOUROT J,BONNEAU M.Comparison of intramuscular adipose tissue cellularity in muscles differing in their lipid content and fibre type composition during rabbit growth[J].Livest Prod Sci,1998,54(1):1-10.
[6]	LOWE C E,O'RAHILLY S,ROCHFORD J J.Adipogenesis at a glance[J].J Cell Sci,2011,124(16):2681-2686.
[7]	CHEN F F,WANG Y Q,TANG G R,et al.Differences between porcine Longissimus thoracis and semitendinosus intramuscular fat content and the regulation of their preadipocytes during adipogenic differentiation[J].Meat Sci,2019,147:116-126.
[8]	WANG Z,GERSTEIN M,SNYDER M.RNA-Seq:a revolutionary tool for transcriptomics[J].Nat Rev Genet,2009,10(1):57-63.
[9]	MARGUERAT S,BHLER J.RNA-seq:from tech-nology to biology[J].Cell Mol Life Sci,2010,67(4):569-579.
[10]	XU T Y,WU J,HAN P,et al.Circular RNA expre-ssion profiles and features in human tissues:a study using RNA-seq data[J].BMC Genomics,2017,18(6):680.
[11]	HIRBEC H,MARMAI C,DUROUX-RICHARD I,et al.The microglial reaction signature revealed by RNAseq from individual mice[J].Glia,2018,66(5):971-986.
[12]	YUN J Y,JIN H G,CAO Y,et al.RNA-Seq analysis reveals a positive role of HTR2A in Adipogenesis in Yan yellow cattle[J].Int J Mol Sci,2018,19(6):1760.
[13]	杜琛,付绍印,高鸿雁,等.绒山羊肌内脂肪细胞成熟前后比较转录组分析[J].畜牧兽医学报,2014,45(5):714-721.DU C,FU S Y,GAO H Y,et al.Transcriptome analysis of intramuscular preadipocytes and matureadipocyte in cashmere goats[J]. Acta Veterinaria et Zootechnica Sinica,2014,45(5):714-721.(in Chinese)
[14]	E G X,ZHAO Y J,MA Y H,et al.Identification of differentially expressed long noncoding RNAs and mRNAs involved with dominant follicle selection in goats using RNA-seq[J].Pak J Zool,2018,50(1):47-56.
[15]	施会彬,王玉琴,张小辉,等.基于转录组测序数据的山羊下丘脑组织繁殖相关信号通路及候选基因的筛选[J].畜牧兽医学报,2019,50(7):1380-1395.SHI H B,WANG Y Q,ZHANG X H,et al.Selection of related signaling pathways and candidate genes related to hypothalamus tissue reproduction in goats based on Transcriptome sequencing data[J].Acta Veterinaria et Zootechnica Sinica,2019,50(7):1380-1395.(in Chinese)
[16]	XU Q,LIN S,WANG Y,et al.Fibroblast growth factor 10(FGF10) promotes the adipogenesis of intramuscular preadipocytes in goat[J].Mol Biol Rep,2018,45(6):1881-1888.
[17]	RIO D C,ARES M JR,HANNON G J,et al.Purification of RNA Using TRIzol (TRI Reagent)[J].Cold Spring Harb Protoc,2010,2010(6):pdb.prot5439.
[18]	PARKHOMCHUK D,BORODINA T,AMSTISLAVSKIY V,et al.Transcriptome analysis by strand-specific sequencing of complementary DNA[J]. Nucleic Acids Res,2009,37(18):e123.
[19]	MORTAZAVI A,WILLIAMS B A,MCCUE K,et al.Mapping and quantifying mammalian transcriptomes by RNA-Seq[J].Nat Methods, 2008,5(7):621-628.
[20]	BRAY N L,PIMENTEL H,MELSTED P,et al.Near-optimal probabilistic RNA-seq quantification[J].Nat Biotechnol,2016, 34(5):525-527.
[21]	许晴,林森,朱江江,等.山羊肌内前体脂肪细胞诱导分化过程中内参基因的表达稳定性分析[J].畜牧兽医学报, 2018, 49(5):907-918.XU Q,LIN S,ZHU J J,et al.The expression stability analysis of reference genes in the process of goat intramuscular preadipocytes differentiation in Goat[J].Acta Veterinaria et Zootechnica Sinica,2018,49(5):907-918.(in Chinese)
[22]	LOVE M I,HUBER W,ANDERS S.Moderated estimation of fold change and dispersion for RNA-seqdata with DESeq2[J]. Genome Biol,2014,15(12):550.
[23]	王建蒙,付绍印,丽春,等.褪黑激素对内蒙古绒山羊皮肤中let-7家族及其靶基因的影响[J].畜牧兽医学报, 2019,50(7):1367-1379.WANG J M,FU S Y,LI C,et al.Effects of melatonin on let-7 family and its target genes in skin of inner mongolian cashmere goats[J].Acta Veterinaria et Zootechnica Sinica,2019,50(7):1367-1379.(in Chinese)
[24]	赵濛,李婕,奈日乐,等.角蛋白关联蛋白基因3.1在绒山羊不同时期皮肤中的表达[J].中国畜牧杂志,2016,52(9):10-13.ZHAO M,LI J,NAI R L,et al.Keratin associated protein 3.1 gene expression in different periods of cashmere goat skin[J].Chinese Journal of Animal Science,2016,52(9):10-13.(in Chinese)
[25]	KIM J B,SPIEGELMAN B M.ADD1/SREBP1 promotes adipocyte differentiation and gene expression linked to fatty acid metabolism[J].Genes Dev,1996,10(9):1096-1107.
[26]	BRENNERS,JACOB F,MESELSON M.An unstable intermediate carrying information from genes to ribosomes for protein synthesis[J].Nature,1961,190(4776):576-581.
[27]	SPINELLI J B,HAIGIS M C.The multifaceted contributions of mitochondria to cellular metabolism[J].Nat Cell Biol, 2018, 20(7):745-754.
[28]	SARASTE M.Oxidative Phosphorylation at the fin de siècle[J].Science,1999,283(5407):1488-1493.
[29]	FERNIE A R,CARRARI F,SWEETLOVE L J.Respiratory metabolism:glycolysis,the TCA cycle and mitochondrial electron transport[J].Curr Opin Plant Biol,2004,7(3):254-261.
[30]	VIGUERIE N,MILLET L,AVIZOU S,et al.Regulation of human adipocyte gene expression by thyroid hormone[J].J Clin Endocrinol Metab,2002,87(2):630-634.
[31]	ROSEN E D,WALKEY C J,PUIGSERVER P,et al.Transcriptional regulation of adipogenesis[J].Genes Dev,2000,14(11):1293-1307.
[32]	VIDAL-PUIG A,SOLANES G,GRUJIC D,et al.UCP3:an uncoupling protein homologue expressed preferentially and abundantly in skeletal muscle and brown adipose tissue[J].Biochem Biophys Res Commun,1997,235(1):79-82.
[33]	BEZAIRE V,SPRIET L L,CAMPBELL S,et al.Constitutive UCP3 overexpression at physiological levels increases mouse skeletal muscle capacity for fatty acid transport and oxidation[J].FASEB J,2005,19(8):977-979.
[34]	乔永.湖羊羔羊不同部位肌肉肌内脂肪沉积相关基因表达的发育性变化研究[D].南京:南京农业大学,2007.QIAO Y.Research on the developmental changes of MRNAS expression levels of the genes related to intramuscular fat in different part of muscles in Hu lamb[D].Nanjing:Nanjing Agricultural University,2007.(in Chinese)
[35]	ABU-ELHEIGA L,MATZUK M M,ABO-HASHEMA K A H,et al.Continuous fatty acid oxidation and reduced fat storage in mice lacking Acetyl-CoA Carboxylase 2[J].Science,2001,291(5513):2613-2616.
[36]	RIANCHO J A,VÁZQUEZ L,GARCÍA-PÉREZ M A,et al.Association of ACACB polymorphisms with obesity and diabetes[J]. Mol Genet Metab,2011,104(4):670-676.
[37]	TILLANDER V,ALEXSON S E H,COHEN D E.Deactivating fatty acids:acyl-CoA thioesterase-mediatedcontrol of lipid metabolism[J]. Trends Endocrinol Metab,2017,28(7):473-484.
[38]	HUNG J Y,CHIANG S R,LIU K T,et al.Over-expression and proliferation dependence of acyl-CoA thioesterase 11 and 13 in lung adenocarcinoma[J].Oncol Lett,2017,14(3):3647-3656.
[39]	MADUREIRA T V,CASTRO L F C,ROCHA E.Acyl-coenzyme A oxidases 1 and 3 in brown trout (Salmo trutta f.fario):can peroxisomal fatty acid β-oxidation be regulated by estrogen signaling?[J].Fish Physiol Biochem,2016,42(1):389-401.
[40]	DELGADO-LISTA J,PEREZ-JIMENEZ F,RUANO J,et al.Effects of variations in the APOA1/C3/A4/A5 gene cluster on different parameters of postprandial lipidmetabolism in healthy young men[J].J Lipid Res,2010,51(1):63-73.
[41]	HAMMARSTEDT A,HEDJAZIFAR S,JENNDAHL L,et al.WISP2 regulates preadipocyte commitment and PPARγactivation by BMP4[J].Proc Natl Acad Sci U S A,2013,110(7):2563-2568.
[42]	ZHANG Y,YU M,DAI M J,et al.miR-450a-5p within rat adipose tissue exosome-like vesicles promotes adipogenic differentiation by targeting WISP2[J].J Cell Sci,2017,130(6):1158-1168.
[43]	曹阳.杜寒杂交羊与小尾寒羊背最长肌全基因组甲基化与转录组比较分析[D].长春:吉林大学,2017.CAO Y.Comparative analysis on genome-wide methylation and RNA-seq in longissimus dorsi muscle between dorper×small tailed Han crossbred and small tailed Han sheep[D].Changchun:Jilin University,2017.(in Chinese)