基于转录组测序筛选宁蒗高原鸡胫长性状的信号通路和关键基因  

Screening Signal Pathways and Key Genes for Shank Length Trait of Ninglang Plateau Chickens Based on Transcriptome Sequencing

作　　者：李金燕 冷堂健 熊宝 郭盘江[1] 吴培福[1] 陈粉粉[1] 周杰珑[1] LI Jinyan;LENG Tangjian;XIONG Bao;GUO Panjiang;WU Peifu;CHEN Fenfen;ZHOU Jielong(College of Biological Science and Food Engineering,Southwest Forestry University,Kunming,650224;Labai Township Ninglang Plateau Chickens Breeding Farm of Ninglang Yi Autonomous County,Lijiang City,Yunnan Province,Ninglang,674309)

机构地区：[1]西南林业大学生物与食品工程学院,昆明650224 [2]云南省丽江市宁蒗彝族自治县拉伯乡宁蒗高原鸡原种保种场,宁蒗674309

出　　处：《基因组学与应用生物学》2024年第11期1794-1812,共19页Genomics and Applied Biology

基　　金：云南省农业基础研究联合专项面上项目(202101BD070001-070)资助。

摘　　要：为挖掘影响宁蒗高原鸡(Gallus gallus domestica)胫发育的信号通路和关键基因,本研究采集第6周龄(6W组)、12周龄(12W组)和18周龄(18W组)健康母鸡跖骨生长板组织进行转录组测序,筛选差异表达基因(differentially expressedgene,DEGs)并随机选取6个基因进行实时荧光定量PCR(real time quantitative PCR,RT-qPCR)验证。通过KEGG通路分类、KEGG和GO富集分析、蛋白质-蛋白质相互作用(protein-protein interaction,PPI)网络可视化分析,以及cytoHubba插件中的最大权团中心度(maximal clique centrality,MCC)、最大邻域分量密度(density of maximum neighborhood component,DMNC)和最大邻域分量(maximum neighborhood component,MNC)等3种算法,筛选出参与宁蒗高原鸡胫发育的候选信号通路和关键基因。结果表明,12W vs 6W组得到66个DEGs,上调DEGs 15个,下调DEGs 51个;18W vs 6W组得到1814个DEGs,上调DEGs 768个,下调DEGs 1046个;18W vs 12W组得到2044个DEGs,上调DEGs 1049个,下调DEGs 995个。6个基因的mRNA相对表达量变化趋势与转录组测序结果一致,表明测序结果可信。KEGG通路分类筛选出312个与生长发育相关的DEGs,KEGG富集分析进一步得到105条信号通路(P<0.05),GO富集分析表明信号转导、ATP结合、细胞外间隙和细胞外区域等条目被显著富集(P<0.05)。综合KEGG和GO富集分析结果,得到15条信号通路和164个DEGs,PPI进一步分析筛选出可能参与宁蒗高原鸡胫发育的信号通路,根据表达量聚类和cytoHubba插件算法,得到调控胫发育的关键基因。本研究筛选到Wnt、TGF-β、PI3K-Aκt、ECM-受体相互作用、细胞周期、黏着斑和刺猬共7条信号通路,以及IHH、COL9A1、COL6A2、FGF7、HGF和WNT4等6个基因为影响宁蒗高原鸡胫发育的候选信号通路和关键基因;并且IHH、COL9A1和WNT4基因在胫发育前期发挥作用,FGF7、HGF和COL6A2基因主要在胫发育后期发挥作用,IHH基因为影响宁蒗高原鸡胫发育的核心基因。该研究为探In order to explore the signaling pathways and key genes affecting the development of shank in Ninglang Plateau Chickens(Gallus gallus domestica),the metatarsal growth plate tissues of healthy hens at 6 weeks of age(6W group),12 weeks of age(12W group)and 18 weeks of age(18W group)were collected for transcriptome sequencing.The differentilly expressed genes(DEGs)were screened and 6 genes were randomly selected for real time quantitative PCR(RT-qPCR)verification.Through KEGG pathway classification,KEGG and GO enrichment analysis,protein-protein interaction(PPI)network visualization analysis and 3 algorithms of maximal clique centrality(MCC),density of maximum neighborhood component(DMNC),and maximum neighborhood component(MNC)in cytoHubba plug-in,candidate signaling pathways and key genes involved in shank development of Ninglang Plateau chickens were screened out.The results showed that 66 DEGs were obtained in the 12W vs 6W group,15 DEGs were up-regulated and 51 DEGs were down-regulated.In the 18W us 6W group,1814 DEGs were obtained,768 DEGs were up-regulated and 1046 DEGs were down-regulated.In the 18W us 12W group,2044 DEGs were obtained,1049 DEGs were up-regulated and 995 DECs were down-regulated.The change trend of mRNA relative expression of 6 genes was consistent with the results of transcriptome sequencing,indicating that the sequencing results were credible.KECG pathway classification screened 312 DEGs related to growth and development.KEGG enrichment analysis further obtained 105 signaling pathways(P<0.05),and signal transduction,ATP binding,extracellular space,and extracellular region were significantly enriched in GO enrichment analysis(P<0.05).Combining KEGG and CO enrichment analysis results,15 signaling pathways and 164 DEGs were obtained.PPI further analyzed and screened out signaling pathways that may be involved in the development of shank in Ninglang Plateau chickens.According to the expression clustering and cytoHubba plug-in algorithm,the key genes regulating the development of shank were obt

关 键 词：宁蒗高原鸡 胫长性状 转录组 信号通路 关键基因 

分 类 号：S831[农业科学—畜牧学]