基于RNA-Seq技术的塔里木马鹿毛色相关差异表达基因筛选及初步分析

doi:10.11843/j.issn.0366-6964.2020.10.009

Abstract

Abstract: The aim of this study was to screen and analyze the hair color related genes of Tarim red deer based on RNA-Seq technology. Using Illumina Hi Seq TM2000 sequencing platform, the skin tissue of Tarim red deer and Tianshan red deer were sequenced by transcriptome. After quality control and assembly, the obtained sequences were compared with the annotations in NR, Swiss-Prot, COG, KOG, KEGG, GO and Pfam databases, and the differentially expressed genes were screened, functional annotated and enriched. The results showed that 25 038 Unigenes with annotated information were obtained by sequencing. The comparative analysis showed that there were 922 differentially expressed genes between Tarim red deer and Tianshan red deer, of which 495 were up-regulated and 427 were down-regulated. The results of GO functional enrichment analysis showed that 568 differentially expressed genes were enriched in 61 GO terms and participated in biological processes, cell components and molecular functions, respectively. The enrichment analysis of KEGG metabolic pathway showed that the most significant metabolic pathway enriched by differentially expressed genes was the ECM-receptor interaction. Real-time fluorescence quantitative PCR (qRT-PCR) was used to analyze the transcriptional level changes of 7 candidate genes related to hair color to verify the accuracy and reliability of transcriptome sequencing results. The expression trend of these genes was consistent with the transcriptome sequencing results. ECM-receptor interaction, protein digestion and absorption, PI3K-Akt signaling pathway and tyrosine related to melanin synthesis may be related to regulation of hair color of Tarim red deer; Candidate genes MITF, Ggt1, VDR, PTPRF, CⅡTA, ARPC5L, POMC may play an important role in hair color formation of Tarim red deer. The results of this study will provide abundant experimental data for studying the molecular regulatory mechanism of hair color related genes and exploring the potential new genes in Tarim red deer.

Key words: Cervus elaphus yarkandensis, RNA-Seq technology, hair color, differentially expressed genes, qRT-PCR

CLC Number:

Q953.3

Subinur ELI, Tayerjan MAMAT, Buweihailiqiemu ABABAIKERI, Mahmut HALIK. Screening and Preliminary Analysis of Tarim Red Deer (Cervus elaphus yarkandensis) Hair Color Differentially Expressed Genes Based on RNA-Seq[J]. Acta Veterinaria et Zootechnica Sinica, 2020, 51(10): 2413-2424.

References 40

[1]	吴宇婷.哺乳动物毛色形成机制与影响因素[J].四川动物,2011,30(6):1003-1007.WU Y T.The formation mechanism of mammal coat color and its possible influencing factors[J].Sichuan Journal of Zoology, 2011,30(6):1003-1007.(in Chinese)
[2]	日沙来提·吐尔地,古力巴哈尔·吾斯曼,艾斯卡尔·买买提,等.新疆部分野生哺乳动物被毛的扫描电镜分析[J].电子显微学报, 2011,30(3):280-284.TURDI R,OSMAN G,MAMAT A,et al.Scanning electron microscope analysis of part of mammalian wildlife hair in Xinjiang[J]. Journal of Chinese Electron Microscopy Society,2011,30(3):280-284.(in Chinese)
[3]	MONTOLIU L,OETTING W S,BENNETT D C.Color Genes[EB/OL].2016.http://www.espcr.org/micemut.
[4]	ANSORGE W J.Next-generation DNA sequencing techniques[J].New Biotechnol,2009,25(4):195-203.
[5]	MOLIN A D,BARUZZO G,CAMILLO B D.Single-cell RNA-sequencing:assessment of differential expression analysis methods[J]. Front Genet,2017,8:62.
[6]	吴良涛,郑敏,华敏,等.基于RNA-Seq筛选鸭肠炎病毒感染鸭脾差异表达免疫相关基因[J].畜牧兽医学报, 2017, 48(2):297-306.WU L T,ZHENG M,HUA M,et al.Screening of differentially expressed immune-related genes from duck spleen with duck enteritis virus infection based on RNA-Seq technology[J].Acta Veterinaria et Zootechnica Sinica,2017,48(2):297-306.(in Chinese)
[7]	付言峰,李兰,赵为民,等.转录组测序分析猪胚胎附植的中期和后期卵巢差异表达基因[J].畜牧兽医学报, 2018,49(9):1879-1888.FU Y F,LI L,ZHAO W M,et al.Differentially expressed genes analysis between mid-implantation and post-implantation in porcine ovary by RNA sequencing[J].Acta Veterinaria et Zootechnica Sinica,2018,49(9):1879-1888.(in Chinese)
[8]	杨峰,车天宇,米璐,等.转录组测序技术在生物学研究中的现状及展望[J].畜牧与兽医,2019,51(3):133-138.YANG F,CHE T Y,MI L,et al.Current status andprospect of biological research on transcriptome sequencing[J].Animal Husbandry & Veterinary Medicine,2019,51(3):133-138.(in Chinese)
[9]	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.
[10]	TARIQ M A,KIM H J,JEJELOWO O,et al.Whole-transcriptome RNAseq analysis from minute amount of total RNA[J].Nucleic Acids Res,2011,39(18):e120.
[11]	RICHARD H,SCHULZ M H,SULTAN M,et al.Prediction of alternative isoforms from exon expression levels in RNA-Seq experiments[J].Nucleic Acids Res,2010,38(10):e112.
[12]	侯扶江,常生华,于应文.我国鹿况简介[J].草业科学,2003,20(11):47-50.HOU F J,CHANG S H,YU Y W.Brief introduction to deer in China[J].Pratacultural Science,2003,20(11):47-50.(in Chinese)
[13]	马合木提·哈力克,艾尼瓦尔·吐米尔,大泰司纪之,等.新疆塔里木马鹿[M].乌鲁木齐:新疆大学出版社,2012.HALIK M,TOMUR A,OHTAISHI N,et al. Tarim Red Deer of XinJiang in China[M]. Urumqi:Xinjiang University Press,2012. (in Chinese)
[14]	汪松.中国濒危动物红皮书[M].北京:科学出版社,1998:417.WANG S. China red data book of endangered animals[M]. Beijing:Science Press,1998:417. (in Chinese)
[15]	JAMES M,DOLEN JR.A deer of many lands-aguide to the subspecies of the red deer Cervus elaphus L[J].Zoonooz, 1988, 62(10):4-34.
[16]	苗德成,李健梅.塔里木马鹿的基本特性[J].特种经济动植物,2008,11(10):6-7.MIAO D C,LI J M. Basic features of the Tarim red deer[J]. Special Economic Animal and Plant,2008,11(10):6-7. (in Chinese)
[17]	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 Res,1997,25(17):3389-3402.
[18]	XIE C,MAO X Z,HUANG J J,et al.KOBAS 2.0:a web server for annotation and identification of enriched pathways and diseases[J].Nucleic Acids Res,2011,39(S2):W316-W322.
[19]	EDDY S R.Profile hidden Markov models[J].Bioinformatics,1998,14(9):755-763.
[20]	LANGMEAD B,TRAPNELL C,POP M,et al.Ultrafast and memory-efficient alignment of short DNA sequences to the human genome[J].Genome Biol,2009,10(3):R25.
[21]	LI B,DEWEY C N.RSEM:accurate transcript quantification from RNA-Seq data with or without a reference genome[J].BMC Bioinformatics,2011,12(1):323.
[22]	ANDERS S,HUBER W.Differential expression analysis for sequence count data[J].Genome Biol,2010,11(10):R106.
[23]	LENG N,DAWSON J A,THOMSON J A,et al.EBSeq:an empirical bayes hierarchical model for inference in RNA-seq experiments[J].Bioinformatics,2013,29(8):1035-1043.
[24]	HARRINGTON N P,SURUJBALLI O P,WATERS W R,et al.Development and evaluation of a real-time reverse transcription-PCR assay for quantification of gamma interferon mRNA to diagnose tuberculosis in multiple animal species[J].Clin Vaccine Immunol,2007,14(12):1563-1571.
[25]	GAUDET P,ŠKUNCA N,HU J C,et al.Primer on the gene ontology[J].Methods Mol Biol,2017,1446:25-37.
[26]	KANEHISA M,FURUMICHI M,TANABE M,et al.KEGG:new perspectives on genomes,pathways,diseases and drugs[J]. Nucleic Acids Res,2017,45(D1):D353-D361.
[27]	MORTAZAVI A,WILLIAMS B A,MCCUE K,et al.Mapping and quantifying mammalian transcriptomes by RNA-Seq[J].Nat Meth,2008,5(7):621-628.
[28]	朱志明,陈红萍,林如龙,等.山麻鸭开产期和产蛋高峰期卵巢组织转录组分析[J].中国农业科学,2016,49(5):998-1007.ZHU Z M,CHEN H P,LIN R L,et al.Transcriptome analysis of ovary tissue in early laying period and egg laying peak period of Shanma ducks[J].Scientia Agricultura Sinica,2016,49(5):998-1007.(in Chinese)
[29]	王文龙,冯陈晨,红梅,等.不同发育阶段斯氏副柔线虫比较转录组学分析[J].中国农业科学,2017,50(23):4644-4655.WANG W L,FENG C C,HONG M,et al.The comparative transcriptome analysis of Parabronema skrjabini at different developmental stages[J].Scientia Agricultura Sinica,2017,50(23):4644-4655.(in Chinese)
[30]	成志敏,张宁芳,王媛媛,等.基于RNA-Seq技术筛选影响猪肌纤维性状的候选基因[J].畜牧兽医学报,2019,50(5):918-929.CHENG Z M,ZHANG N F,WANG Y Y,et al.Screening of candidate genes for muscle fiber characteristics in pig using RNA-Seq[J].Acta Veterinaria et Zootechnica Sinica,2019,50(5):918-929.(in Chinese)
[31]	MIAO X Y,LUO Q M.Genome-wide transcriptome analysis between Small-tail Han sheep and the Surabaya fur sheep using high-throughput RNA sequencing[J].Reproduction,2013,145(6):587-596.
[32]	BENDER D A.Tricarboxylic acid cycle[M]//CABALLERO B.Encyclopedia of Food Sciences and Nutrition.2nd ed.Amsterdam:Academic Press,2003:5851-5856.
[33]	BRILLANTES A M B,ONDRIAS K,SCOTT A,et al.Stabilization of calcium release channel (ryanodine receptor) function by FK506-binding protein[J].Cell,1994,77(4):513-523.
[34]	JAYARAMAN T,MARKS A R.Rapamycin-FKBP12 blocks proliferation,induces differentiation,and inhibits cdc2 kinase activity in a myogenic cell line[J].J Biol Chem,1993,268(34):25385-25388.
[35]	MARIMAN E C M,WANG P.Adipocyte extra-cellular matrix composition,dynamics and role in obesity[J].Cell Mol Life Sci,2011,67(8):1277-1292.
[36]	马淑慧,李亚楠,赫晓燕,等.POMC在不同毛色绵羊皮肤中的表达与定位[J].中国农业科学,2015,48(9):1818-1824.MA S H,LI Y N,HAO X Y,et al.POMC expression and localization in sheep skin with different coat color[J].Scientia Agricultura Sinica,2015,48(9):1818-1824.(in Chinese)
[37]	ASCHAFFENBURG R,DREWRY J.Genetics of the β-lactoglobulins of cow's milk[J].Nature,1957,180(4582):376-378.
[38]	KUSS A W,GOGOL J,GELDERMANN H.Associations of a polymorphic AP-2 binding site in the 5'-flanking region of the bovine β-lactoglobulin gene with milk proteins[J].J Dairy Sci,2003,86(6):2213-2218.
[39]	杨广礼.动物毛色形成的调控机制研究[J].黑龙江畜牧兽医,2014(3):45-48.YANG G L. Regulation mechanism of animal coat color formation[J]. Heilongjiang Animal Science and Veterinary Medicine,2014(3):45-48. (in Chinese)
[40]	刘小辉.坝上长尾鸡资源保护及羽色多样性分子遗传基础[D].保定:河北农业大学,2018.LIU X H.Resources protection and molecular genetic basis of plumage color diversity of Bashang Long-tail chickens[D]. Baoding:Hebei Agricultural University,2018.(in Chinese)

Screening and Preliminary Analysis of Tarim Red Deer (Cervus elaphus yarkandensis) Hair Color Differentially Expressed Genes Based on RNA-Seq

PDF (PC)

Knowledge

Abstract

Cite this article

share this article

References 40

Related Articles 15

Recommended Articles

Metrics

Comments

[1]	GAO Yawei, PENG Di, SUN Zhaoyang, YAN Ziyue, CUI Kai, MA Zefang. Mining the Molecular Mechanism of Exogenous Melatonin Affecting the Development of Mink Ovary Based on Transcriptome Data [J]. Acta Veterinaria et Zootechnica Sinica, 2024, 55(2): 607-618.
[2]	WANG Jiandong, TANG Yulin, WANG Min, ZHANG Baosuo, YANG Fuqiang, GAO Haihui, YU Yang, GUO Yansheng. The Mechanism of Lycium barbarum Polysaccharide against Immunosuppression Induced by Cyclophosphamide in Chicks Based on RNA-Seq Technique [J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(8): 3519-3532.
[3]	WANG Zixuan, WANG Qiao, ZHANG Jin, Astrid Lissette Barreto Sánchez, ZHENG Maiqing, LI Qinghe, CUI Huanxian, AN Bingxing, ZHAO Guiping, WEN Jie, LI Hegang. Transcriptome Based Screening of Functional Genes Related to Heat Stress Resistance in Beijing You Chickens and Guangming Broilers [J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(5): 1905-1914.
[4]	OU Zhengmiao, ZHOU Jiawen, LIU Lili, WU Yun, CHEN Fenfen. Screening and Expression Analysis of Genes Related to Lipid Metabolism in Liver Tissue of Wuliangshan Sooty Chicken Based on RNA-Seq [J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(3): 976-988.
[5]	HAN Haoyuan, LI Shikai, YANG Ruiqiao, LI Manman, LI Jun, HA Si, ZHAO Jinyan, WEI Hongfang, QUAN Kai. Mining Key Candidate Genes for High Reproduction Performance of Huai Goats Based on Transcriptome Sequencing [J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(12): 5077-5090.
[6]	CAI Jiawei, ZHANG Chen, JIN Rongshuai, BAO Zhiyuan, ZHANG Xiyu, WANG Fan, ZHAI Pin, ZHAO Bohao, CHEN Yang, TANG Xianwei, WU Xinsheng. Analysis of Testicular Tissue Morphology and Semen Transcriptome of Male Rabbits under Heat Stress [J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(11): 4653-4663.
[7]	ZHANG Yinliang, YUE Qiaoxian, HUANG Chenxuan, CHEN Hui, WANG Dehe, ZHOU Rongyan. Construction of Tibia Transcript Profile of Laying Hens at the Early and Late Laying Stages and Analysis of Genes Related to Bone Metabolism [J]. Acta Veterinaria et Zootechnica Sinica, 2023, 54(10): 4164-4173.
[8]	MAO Yanni, CHANG Jiawei, LI Na, WANG Xin, KANG Xinyun, MA Qiang, MA Liang, WANG Guiqin. Transcriptome Differential Expression Analysis of Staphylococcus aureus in Biofilm State and Planktonic State [J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(8): 2697-2707.
[9]	WANG Chuanchuan, MU Tong, FENG Xiaofang, YU Baojun, ZHANG Juan, GU Yaling. Identification of Key Candidate Genes for Milk Fat Metabolism in Dairy Cows Based on Transcriptome Sequencing [J]. Acta Veterinaria et Zootechnica Sinica, 2022, 53(10): 3421-3433.
[10]	ZHANG Zhidong, RYAN Eoin, LI Yanmin. Quantitative Analysis of FMDV Receptor Integrin mRNA Levels in Ovine Tissues at Different Developmental Stages [J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(9): 2626-2632.
[11]	MA Xin, YANG Guang, LI Xiangyu, LIU Yueyue, CHEN Yong. Effects of Intravenous Infusion of Atrial Natriuretic Peptide on Blood Lipid Metabolism Hormones and Transcriptome Analysis of Tail Adipose Tissue of Altay Sheep [J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(11): 3076-3088.
[12]	FANG Yan, ZHOU Junwen, GUAN Weijun, JIANG Lin, PU Yabin, ZHAO Qianjun, HE Xiaohong, MA Yuehui. Exploring the Molecular Mechanism of Bactrian Camel's Desert Adaptation Based on Rumen Transcriptome [J]. Acta Veterinaria et Zootechnica Sinica, 2021, 52(1): 77-87.
[13]	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.
[14]	SHI Huibin, WANG Yuqin, ZHANG Xiaohui, WU Zhibo, YANG Fang. 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.
[15]	YU Yan, GU Yu, WANG Chengdong, HU Jing, YANG Aining, YANG Qiao, LI Desheng, ZHANG Hemin, DENG Linhua, LING Shanshan, ZUO Zhicai, PENG Guangneng, ZHONG Zhijun, MA Xiaoping. Dynamic Changes of Dectin-1, TLR-2,TLR-4 and Cytokines in Mice Infected with Microsporum gypseum [J]. ACTA VETERINARIA ET ZOOTECHNICA SINICA, 2019, 50(7): 1449-1457.