基于转录组学的安格斯牛和中国西门塔尔牛脂肪沉积相关差异表达基因研究

doi:10.7668/hbnxb.20190684

摘要/Abstract

摘要： 为探明安格斯牛和中国西门塔尔牛脂肪沉积的相关差异表达基因，以其背膘组织为研究对象，提取总RNA，建立测序文库后利用华大BGISEQ-500进行测序，经生物信息学分析筛选出差异表达基因，并对差异表达基因进行GO功能富集和KEGG通路分析，最后采用荧光定量qPCR（qRT-PCR）验证转录组测序数据的准确性。中国西门塔尔牛对安格斯牛转录组数据结果显示有1 296个差异表达基因，其中上调基因802个，下调基因494个；对转录组数据的GO功能富集分析结果显示在生物过程、分子功能、细胞组分分别富集到2 205，890，3 149个条目，KEGG通路分析结果表明，差异基因显著富集在6个分支上。对其中与脂肪沉积相关的cAMP通路进行分析，筛选出与脂肪沉积相关的差异表达基因3个，分别为DKKL1、ACACB和PTGS2，并通过qRT-PCR对基因DKKL1、ACACB和PTGS2进行了验证，初步确定其为脂肪沉积的候选基因。通过对中国西门塔尔牛和安格斯牛的背膘组织测序获得的大量差异表达基因，为后续对脂肪沉积相关基因的筛选奠定了基础。

关键词: 牛, 脂肪沉积, DKKL1, ACACB, PTGS2, qRT-PCR

Abstract: To identify differentially expressed genes related to fat deposition in Angus and Chinese simmental cattle and analyze differentially expressed gene functions and pathways, the total RNA was extracted from the backfat, then sequenced by BGISEQ-500. The differentially expressed genes were screened by bioinformatics analysis, GO and KEGG enrichment analysis were used to study the function of differentially expressed genes and the candidate genes related to fat deposition were screened. The accuracy of RNA-seq results was verified by qRT-PCR. The results showed that there were 1 296 significant differentially expressed genes in Angus cattle compared with Chinese simental cattle, including 802 up genes and 494 down genes. Differential transcripts of GO function enrichment analysis revealed that molecular functions, cell components, and biological processes were enriched to 890 entries, 3 149 entries,and 2 205 entries. KEGG analysis showed that the differential transcripts were significantly enriched in 6 pathways. The cAMP pathway associated with fat deposition was analyzed, and the differentially expressed genes associated with fat deposition were screened out, the genes DKKL1, ACACB and PTGS2 were verified by qRT-PCR, and they were preliminarily identified as candidate genes for fat deposition. At the same time, a large number of differentially expressed genes were obtained by sequencing the backfat tissues of Chnese simmental cattle and Angus cattle, which laid a foundation for the subsequent screening of genes related to fat deposition.

Key words: Cattle, Fat deposition, DKKL1, ACACB, PTGS2, qRT-PCR

中图分类号:

Q78
S832

韩信嵘, 吴慧光, 吴江鸿, 王国富, 高树新. 基于转录组学的安格斯牛和中国西门塔尔牛脂肪沉积相关差异表达基因研究[J]. 华北农学报, 2020, 35(6): 208-216. doi: 10.7668/hbnxb.20190684.

HAN Xinrong, WU Huiguang, WU Jianghong, WANG Guofu, GAO Shuxin. Study on Differentially Expressed Genes Related to Fat Deposition Gene in Angus and Chinese Simmental Cattle Base on Transcriptome[J]. ACTA AGRICULTURAE BOREALI-SINICA, 2020, 35(6): 208-216. doi: 10.7668/hbnxb.20190684.

http://www.hbnxb.net/CN/Y2020/V35/I6/208

参考文献

[1] 王国富, 吴慧光, 赵新海, 季守财, 吴红君, 王东升, 刘胜敏, 高树新. 安格斯牛、海福特牛和中国西门塔尔牛的部分胴体性状比较分析[J].内蒙古民族大学学报(自然科学版), 2010,25(5):535-537.doi:10.3969/j.issn.1671-0185.2010.05.016.
Wang G F, Wu H G, Zhao X H,Ji S C, Wu H J, Wang D S, Liu S M, Gao S X. Comparative study on some carcass traits in Angus,Hereford and Chinese Simmental cattles[J]. Journal of Inner Mongolia University for Nationalities(Natural Sciences), 2010,25(5):535-537.
[2] 王小梅. 内蒙古地区不同品种肉牛生产性能和肉品质及脂肪代谢相关基因mRNA表达量的比较研究[D].呼和浩特:内蒙古农业大学, 2012.
Wang X M. Comparative study on different breed beef cattle's production,meat quality and mRNA expression of fatty metabolism related factors in Inner Mongolia[D].Hohhot:Inner Mongolia Agricultural University, 2012.
[3] 张国梁, 刘基伟, 胡成华, 吴健, 张嘉保. 西门塔尔、夏洛来、安格斯杂交牛不同营养水平的育肥效果分析[J].中国畜牧兽医, 2007,34(7):10-12.doi:10.3969/j.issn.1671-7236.2007.07.002.
Zhang G L, Liu J W, Hu C H, Wu J, Zhang J B. Fattening effect analysis of Simmental, Charolais and Angus crossbred cattle at different nutritional levels[J]. China Animal Husbandry & Veterinary Medicine, 2007,34(7):10-12.
[4] 史新平, 陈燕, 徐玲, 金生云, 宝金山, 牛红, 王泽昭, 李俊雅, 孙少华.肉用西门塔尔牛与和牛杂交群体的肉品质分析[J].中国畜牧兽医, 2018,45(4):953-960.doi:10.16431/j.cnki.1671-7236.2018.04.05.
Shi X P,Chen Y,Xu L,Jin S Y,Bao J S,Niu H,Wang Z Z,Li J Y,Sun S H.Analysis of beef Simmental and Wagyu hybrid cattle on meat quality[J]. Chinese Animal Husbandry and Veterinary Medicine,2018,45(4):953-960.
[5] Mazzucco J P, Goszczynski D E, Ripoli M V, Melucci L M, Pardo A M, Colatto E, Rogberg-Mu n oz A, Mezzadra C A, Depetris G J, Giovambattista G, Villarreal E L. Growth, carcass and meat quality traits in beef from Angus, Hereford and cross-breed grazing steers, and their association with SNPs in genes related to fat deposition metabolism[J]. Meat Science, 2016,114:121-129.doi:10.1016/j.meatsci.2015.12.018.
[6] Pertea M, Pertea G M, Antonescu C M, Chang T C, Mendell J T, Salzberg S L. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads[J]. Nature Biotechnology,2015,33(3):290-295.doi:10.1038/nbt.3122.
[7] Trapnell C, Roberts A, Goff L, Pertea G, Kim D, Kelley D R, Pimentel H, Salzberg S L, Rinn J L, Pachter L. Differential gene and transcript expression analysis of RNA-seq experiments with TopHat and Cufflinks[J]. Nature Protocols,2012,7(3):562-578.doi:10.1038/nprot.2012.016.
[8] Kong L, Zhang Y, Ye Z Q, Liu X Q, Zhao S Q, Wei L P, Gao G. CPC:assess the protein-coding potential of transcripts using sequence features and support vector machine[J]. Nucleic Acids Research, 2007,35(2):345-349.doi:10.1093/nar/gkm391.
[9] McKenna A, Hanna M, Banks E, Sivachenko A, Cibulskis K, Kernytsky A, Garimella K, Altshuler D, Gabriel S, Daly M, DePristo M A. The Genome Analysis Toolkit:a MapReduce framework for analyzing next-generation DNA sequencing data[J]. Genome Research, 2010,20(9):1297-1303.doi:10.1101/gr.107524.110.
[10] Shen S H, Park J W, Lu Z X, Lin L, Henry M D, Wu Y N, Zhou Q, Xing Y. rMATS:robust and flexible detection of differential alternative splicing from replicate RNA-Seq data[J]. Proceedings of the National Academy of Sciences of the United States of America, 2014,111(51):5593-5601.doi:10.1073/pnas.1419161111.
[11] Wang L K, Feng Z X, Wang X, Wang X W, Zhang X G. DEGseq:an R package for identifying differentially expressed genes from RNA-seq data[J]. Bioinformatics, 2010,26(1):136-138.doi:10.1093/bioinformatics/btp612.
[12] Kanehisa M, Goto S, Hattori M, Aoki-Kinoshita K F, Itoh M, Kawashima S, Katayama T, Araki M, Hirakawa M. From genomics to chemical genomics:new developments in KEGG[J]. Nucleic Acids Research,2006,34(S1):354-357.doi:10.1093/nar/gkj102.
[13] Kanehisa M, Araki M, Goto S, Hattori M, Hirakawa M, Itoh M, Katayama T, Kawashima S, Okuda S, Tokimatsu T, Yamanishi Y. KEGG for linking genomes to life and the environment[J]. Nucleic Acids Research,2008,36(S1):480-484.doi:10.1093/nar/gkm882.
[14] Kim D, Langmead B, Salzberg S L. HISAT:a fast spliced aligner with low memory requirements[J]. Nature Methods,2015,12(4):357-360.doi:10.1038/nmeth.3317.
[15] Okeudo N J, Moss B W. Interrelationships amongst carcass and meat quality characteristics of sheep[J]. Meat Science, 2005,69(1):1-8.doi:10.1016/j.meatsci.2004.04.011.
[16] Ravnskjaer K, Madiraju A, Montminy M. Role of the cAMP pathway in glucose and lipid metabolism[J]. Handbook of Experimental Pharmacology,2016,233:29-49.doi:10.1007/164_2015_32.
[17] 李常青. 外源性cAMP对秦川犊牛生长发育及调控机理的研究[D].杨凌:西北农林科技大学, 2010.
Li C Q.The effect on the growth and controt mechanism in Qinchuan calves by exogenous camp[D]. Yangling:Northwest A&F University, 2010.
[18] Acin-Perez R, Salazar E, Kamenetsky M, Buck J, Levin L R, Manfredi G. Cyclic AMP produced inside mitochondria regulates oxidative phosphorylation[J]. Cell Metabolism,2009,9(3):265-276.doi:10.1016/j.cmet.2009.01.012.
[19] Okuda H, Morimoto C, Tsujita T. Effect of substrates on the cyclic AMP-dependent lipolytic reaction of hormone-sensitive lipase[J]. Journal of Lipid Research, 1994,35(7):1267-1273.
[20] McNeel R L, Mersmann H J. β-adrenergic receptor subtype transcripts in porcine adipose tissue[J]. Journal of Animal Science, 1995,73(7):1962-1971.doi:10.2527/1995.7371962x.
[21] Cruciat C M, Niehrs C. Secreted and transmembrane wnt inhibitors and activators[J]. Cold Spring Harbor Perspectives in Biology, 2013,5(3):a015081.doi:10.1101/cshperspect.a015081.
[22] Mao B Y, Wu W, Davidson G, Marhold J, Li M F, Mechler B M, Delius H, Hoppe D, Stannek P, Walter C, Glinka A, Niehrs C. Kremen proteins are Dickkopf receptors that regulate Wnt/beta-catenin signalling[J]. Nature, 2002,417(6889):664-667.doi:10.1038/nature756.
[23] Krupnik V E, Sharp J D, Jiang C, Robison K, Chickering T W, Amaravadi L, Brown D E, Guyot D, Mays G, Leiby K, Chang B, Duong T, Goodearl A D J, Gearing D P, Sokol S Y, McCarthy S A. Functional and structural diversity of the human Dickkopf gene family[J]. Gene,1999,238(2):301-313.doi:10.1016/s0378-1119(99)00365-0.
[24] Dakhova O, O'Day D, Kinet N, Yucer N, Wiese M, Shetty G, Ducy P. Dickkopf-like1 regulates postpubertalspermatocyte apoptosis and testosterone production[J]. Endocrinology, 2009,150(1):404-412.doi:10.1210/en.2008-0673.
[25] 李亮, 程彦伟. 乙酰辅酶A羧化酶在治疗肥胖中的潜在作用[J].生命的化学, 2007,27(2):180-182.doi:10.3969/j.issn.1000-1336.2007.02.030.
Li L, Cheng Y W. The potential effects of acetyl-CoA carboxylase in curing obesity[J]. Chemistry of Life, 2007,27(2):180-182.
[26] Oh W, Abu-Elheiga L, Kordari P, Gu Z W, Shaikenov T, Chirala S S, Wakil S J. Glucose and fat metabolism in adipose tissue of acetyl-CoA carboxylase 2 knockout mice[J]. Proceedings of the National Academy of Sciences of the United States of America,2005,102(5):1384-1389.doi:10.1073/pnas.0409451102.
[27] Abu-Elheiga L, Oh W, Kordari P, Wakil S J. Acetyl-CoA carboxylase 2 mutant mice are protected against obesity and diabetes induced by high-fat/high-carbohydrate diets[J]. Proceedings of the National Academy of Sciences of the United States of America, 2003,100(18):10207-10212.doi:10.1073/pnas.1733877100.
[28] 高丽南, 王思珍, 曹颖霞, 朱晓岩, 刘国文. 营养与品种对肉牛血清中AMPK活性与血液生化指标关系的研究[J].中国畜牧杂志. 2012,48(9):29-33.
doi,10.3969/j.issn.0258-7033.2012.09.008.
Gao L N, Wang S Z, Cao Y X, Zhu X Y, Liu G W.The effects of nutrition and bread on AMPK activity in serum and relation ship between AMPK activity and blood indicators[J]. Chinese Journal of Animal Science, 2012,48(9):29-33.
[29] Warner T D, Mitchell J A. Cyclooxygenase-3(COX-3):filling in the gaps toward a COX continuum[J]. Proceedings of the National Academy of Sciences of the United States of America, 2002,99(21):13371-13373.doi:10.1073/pnas.222543099.
[30] Hu E, Tontonoz P, Spiegelman B M. Transdifferentiation of myoblasts by the adipogenic transcription factors PPAR gamma and C/EBP alpha[J]. Proceedings of the National Academy of Sciences of the United States of America,1995,92(21):9856-9860.doi:10.1073/pnas.92.21.9856.
[31] Siebert A, Goren I, Pfeilschifter J, Frank S. Anti-Inflammatory effects of rosiglitazone in obesity-impaired wound healing depend on adipocyte differentiation[J]. PLoS One, 2016,11(12):e0168562.doi:10.1371/journal.pone.0168562.
[32] Belda B J, Thompson J T, Eser P O, Vanden Heuvel J P. 10e12z CLA alters adipocyte differentiation and adipocyte cytokine expression and induces macrophage proliferation[J]. The Journal of Nutritional Biochemistry, 2012,23(5):510-518.doi:10.1016/j.jnutbio.2011.02.009.
[33] Kennedy A, Overman A, Lapoint K, Hopkins R, West T, Chuang C C, Martinez K, Bell D, McIntosh M. Conjugated linoleic acid-mediated inflammation and insulin resistance in human adipocytes are attenuated by resveratrol[J]. Journal of Lipid Research,2009,50(2):225-232.doi:10.1194/jlr.M800258-JLR200.
[34] Yang Y, Wei H K, Song T X, Cai A L, Zhou Y F, Peng J, Jiang S W, Peng J. E4BP4 mediates glucocorticoid-regulated adipogenesis through COX2[J]. Molecular and Cellular Endocrinology, 2017,450:43-53.doi:10.1016/j.mce.2017.04.015.
[35] Fain J N, Ballou L R, Bahouth S W. Obesity is induced in mice heterozygous for cyclooxygenase-2[J]. Prostaglandins & Other Lipid Mediators, 2001,65(4):199-209.doi:10.1016/S0090-6980(01)00136-8.

选择文件类型/文献管理软件名称

选择包含的内容