中国农业科学 ›› 2019, Vol. 52 ›› Issue (19): 3471-3484.doi: 10.3864/j.issn.0578-1752.2019.19.016

• 畜牧·兽医·资源昆虫 • 上一篇    下一篇

苏博美利奴羊毛囊发育相关lncRNA与mRNA共表达网络的构建

陈华峰1,2,田可川2(),黄锡霞1(),阿布来提·苏来曼1,何军敏3,田月珍2,徐新明2,付雪峰2,赵冰茹4,朱桦1,哈尼克孜·吐拉甫2   

  1. 1 新疆农业大学动物科学学院,乌鲁木齐 830052
    2 新疆畜牧科学院畜牧研究所,乌鲁木齐 830011
    3 甘肃农业大学动物科学技术学院,兰州 730070
    4 中国农业大学动物科技学院,北京 100083
  • 收稿日期:2018-12-19 接受日期:2019-07-03 出版日期:2019-10-01 发布日期:2019-10-11
  • 通讯作者: 田可川,黄锡霞
  • 作者简介:陈华峰,E-mail:1504159731@qq.com。
  • 基金资助:
    国家自然科学基金(31460593);国家自然科学基金(31360543);国家现代农业产业技术体系建设项目(CARS-40);新疆绒毛用羊遗传育种与繁殖自治区重点实验室项目(XJYS1105);新疆农业大学研究生科研创新项目(XJAUGRI-2017013)

Construction of Co-expression Network of lncRNA and mRNA Related to Hair Follicle Development of Subo Merino Sheep

CHEN HuaFeng1,2,TIAN KeChuan2(),HUANG XiXia1(),Ablat Sulayman1,HE JunMin3,TIAN YueZhen2,XU XinMing2,FU XueFeng2,ZHAO BingRu4,ZHU Hua1,Hanikezi Tulafu2   

  1. 1 College of Animal Science, Xinjiang Agricultural University, Urumqi 830052
    2 Institute of Animal Science, Xinjiang Academy of Animal Science, Urumqi 830011
    3 College of Animal Science and Technology, Gansu Agricultural University, Lanzhou 730070
    4 College of Animal Science and Technology, China Agricultural University, Beijing 100083
  • Received:2018-12-19 Accepted:2019-07-03 Online:2019-10-01 Published:2019-10-11
  • Contact: KeChuan TIAN,XiXia HUANG

摘要:

【目的】 随着高通量测序技术的发展和完善,海量的转录组测序数据随之涌现,基因网络的方法也越来越多被用于研究中;细毛羊毛囊发育受多个基因及lncRNA共同调控,单独研究某一分子并不足以发现其调控机制,本研究旨在构建毛囊发育相关lncRNA和mRNA的共表达网络,挖掘细毛羊毛囊发育的潜在候选基因。【方法】 以苏博美利奴羊胎龄65 d(D65)、85 d(D85)、105 d(D105)、135 d(D135)的胎儿以及出生后7 d(A7)、30 d(A30)的羔羊肩胛部皮肤组织,每个时期有3个生物学重复,共18个样本,进行转录组测序以获得6个不同发育时期的lncRNA与mRNA表达谱数据,筛选出相邻时期的差异表达lncRNA与mRNA,运用加权基因共表达网络分析(weighted gene co-expression network analysis,WGCNA)方法构建共表达模块,使用DAVID在线工具进行GO(gene ontology)和 KEGG pathway富集分析以找到与毛囊发育相关的模块,最后从目标模块筛选高连通度的lncRNA与mRNA使用Cytoscape软件进行网络可视化。【结果】 从表达谱数据中筛选得到9 070个差异表达的lncRNA与mRNA,运用WGCNA方法得到11个模块,使用DAVID富集分析发现honeydew1、paleturquoise、skyblue2模块中的基因富集在皮肤发育、毛囊发育、毛囊形态发生、Wnt信号通路负调控、细胞粘附等生物过程,以及Wnt信号通路、TGF-β信号通路、Hedgehog信号通路、紧密连接、MAPK信号通路、ECM-受体相互作用等毛囊发育相关的信号通路,并筛选这3个模块中连通度高的lncRNA和mRNA构建了子网络,得到包括TGFB2CTSBSFNSPINT1FAM83GGSDMAMPZL3VIMCRABP1在内的毛囊发育相关基因,预测出ENSOART00000029117、TCONS_00489976、TCONS_00376759等24个lncRNA的可能靶基因。【结论】 首次运用WGCNA方法构建了苏博美利奴羊毛囊发育相关lncRNA与mRNA的共表达网络,鉴别出3个毛囊发育相关的共表达模块,发现多个与毛囊发育相关的潜在候选基因,并预测出24个lncRNA可能的靶基因。

关键词: 苏博美利奴羊, WGCNA, mRNA, lncRNA, 毛囊发育

Abstract:

【Objective】 With the development and improvement of high-throughput sequencing technology, massive transcriptome sequencing data has emerged, and more and more methods of gene network have been used in research; the hair follicles development of fine wool sheep is regulated by multiple genes and lncRNAs. It is not enough to study a certain molecule alone to discover its regulatory mechanism. The aim of this study was to construct a co-expression network of lncRNAs and mRNAs related to hair follicle development, and to explore potential candidate genes for hair follicle development. 【Method】 The study was performed on fetal skin tissues of the 65th, 85th, 105th, and 135th days and lambs skin tissues of the 7th and 30th days from Subo Merino sheep, with 3 biological replicates in each period, a total of 18 samples were sequenced by transcriptome to obtain lncRNAs and mRNAs expression profiles of 6 different developmental stages, and the differential expressed lncRNAs and mRNAs in adjacent periods was screened to construct a co-expression module by Weighted Gene Co-expression Network Analysis (WGCNA) method, and GO (Gene Ontology) and KEGG pathway enrichment analysis of DAVID online tool was used to find hair follicle related modules. Finally, the high-interconnect lncRNAs and mRNAs were screened from the target module for network visualization using Cytoscape software.【Result】 From the expression profile data, 9070 differentially expressed lncRNAs and mRNAs were screened, and 11 modules were obtained by WGCNA method. The DAVID enrichment analysis revealed that the genes in the honeydew1 module, paleturquoise module, and skyblue2 module were involved in biological processes such as skin development, hair follicle development, hair follicle morphogenesis, negative regulation of Wnt signaling pathway, and cell adhesion, et al; and the signaling pathways involved in hair follicle development such as Wnt signaling pathway, TGF-β signaling pathway, Hedgehog signaling pathway, tight junction, MAPK signaling pathway, ECM-receptor interaction, et al. The sub-networks of lncRNAs and mRNAs with high connectivity were screened in these 3 modules, and the hair follicle development related genes including TGFB2, CTSB, SFN, SPINT1, FAM83G, GSDMA, MPZL3, VIM and CRABP1 were obtained, and possible target genes of 24 lncRNAs including ENSOART00000029117, TCONS_00489976, TCONS_00376759 were predicted.【Conclusion】 In this study, we first constructed a co-expression network of lncRNAs and mRNAs related to the hair follicles development of Subo merino sheep using WGCNA method, identified 3 co-expression modules related to hair follicle development, and found several potential candidate genes related to hair follicle development, and predicted possible target genes of 24 lncRNAs.

Key words: Subo Merino sheep, WGCNA, mRNA, lncRNA, hair follicle development