黄鳝生长差异基因发掘及其调控机制初步研究

霍欢欢, 刘瑜, 周秋白, 郭枫, 隗黎丽, 彭墨, 张燕萍, 陈文静

霍欢欢, 刘瑜, 周秋白, 郭枫, 隗黎丽, 彭墨, 张燕萍, 陈文静. 黄鳝生长差异基因发掘及其调控机制初步研究[J]. 南方水产科学, 2020, 16(1): 1-8. DOI: 10.12131/20190176
引用本文: 霍欢欢, 刘瑜, 周秋白, 郭枫, 隗黎丽, 彭墨, 张燕萍, 陈文静. 黄鳝生长差异基因发掘及其调控机制初步研究[J]. 南方水产科学, 2020, 16(1): 1-8. DOI: 10.12131/20190176
HUO Huanhuan, LIU Yu, ZHOU Qiubai, GUO Feng, WEI Lili, PENG Mo, ZHANG Yanping, CHEN Wenjing. Primary study on differentially expressed genes screening of Monopterus albus and their regulation mechanism[J]. South China Fisheries Science, 2020, 16(1): 1-8. DOI: 10.12131/20190176
Citation: HUO Huanhuan, LIU Yu, ZHOU Qiubai, GUO Feng, WEI Lili, PENG Mo, ZHANG Yanping, CHEN Wenjing. Primary study on differentially expressed genes screening of Monopterus albus and their regulation mechanism[J]. South China Fisheries Science, 2020, 16(1): 1-8. DOI: 10.12131/20190176

黄鳝生长差异基因发掘及其调控机制初步研究

基金项目: 国家特色淡水鱼产业体系建设专项资金 (CARS-46);国家自然科学基金项目 (31360641)
详细信息
    作者简介:

    霍欢欢 (1986—),男,博士,讲师,从事鱼类营养及免疫研究。E-mail: hhh16430@163.com

    通讯作者:

    周秋白 (1965—),男,博士,教授,从事鱼类营养研究。E-mail: zhouqiubai@163.com

  • 中图分类号: S 917.4

Primary study on differentially expressed genes screening of Monopterus albus and their regulation mechanism

  • 摘要: 为揭示黄鳝 (Monopterus albus) 生长相关基因的调控机制,对相同亲本具有显著性生长差异的黄鳝肝脏进行了转录组测序分析。结果显示,转录组测序共得到19 149个基因,其中差异基因598个,差异基因中有303个基因显著上调,295个显著下调。KEGG通路分析发现,598个差异基因分属262条通路中,其中有38条通路显著富集。GO功能注释发现,与生长相关的差异基因有7个,分别为col1α1、nkx6.1、nnosplexina4、igfbp1、pcgf1和h3.3,这些基因表达水平的变化可能对黄鳝神经、内分泌和消化系统的发育及生理活动产生了调节作用从而影响了黄鳝的生长。结合KEGG通路分析发现,col1α1所在的利什曼病通路和nnos所在的精氨酸和脯氨酸代谢通路富集显著,说明其对黄鳝生长具有重要调节作用。
    Abstract: In order to reveal the regulatory mechanism of growth related genes in eels (Monopterus albus), we carried out a transcriptome sequencing analysis on the liver of the eels that had significant growth difference with the same parent. A total of 19 149 genes were obtained by transcriptome sequencing, among which 598 were differentially expressed (303 and 295 were significantly up-regulated and down-regulated, respectively). KEGG pathway analysis shows that 598 differentially expressed genes belonged to 262 pathways, among which 38 pathways are significantly enriched. The GO functional annotation reveals that the seven differentially related genes were col1α1, nkx6.1, nnos, plexina4, igfbp1, pcgf1and h3.3, and variation in the expression levels of these genes had regulated the development and physiological activities of the nervous system, endocrine system and digestive system of the eels, which then affected their growth. Combined with KEGG pathway analysis, it is found that the leishmaniasis pathway, in which col1α1 is located, and the arginine and proline metabolic pathway, in which nnos is located, are significantly enriched, revealing their important influence on the growth of eels.
  • 黄鳝 (Monopteru albus),俗名鳝鱼,广泛分布于我国各湖泊、水库和稻田等淡水水域。其肉质鲜美,含有丰富的必需氨基酸和脂肪酸及其他特殊营养素,具有很高的药用和滋补功能,一直深受消费者青睐[1]。黄鳝是目前我国大力推广养殖的重要名特优水产品之一,2016年我国人工养殖黄鳝产量达到38.6万吨。从市场角度看,规格一致的鱼苗或成鱼有利于养殖品种的商品化,但在鱼类养殖过程中,总有部分鱼生长缓慢,这不仅浪费养殖空间,也导致饲料浪费,严重影响了养殖经济效益。黄鳝养殖在这方面尤为突出,即使是同一亲本所产、在相同环境中生长的黄鳝,这种生长差异依然存在且显著。影响鱼类生长的因素很多,如种质差异、营养水平、温度、密度、水质等[2],但目前对黄鳝生长差异产生的机制尚不清楚,亟待研究探讨。

    近年来,转录组学被广泛应用于生物学研究,对于功能基因发掘、转录调控机制、分子标记开发和信号通路等各方面的研究具有重要作用[3]。目前,转录组学已广泛应用于鱼类发育[4]、进化[5]、免疫[6]和抗病机理[7]研究中。本文通过转录组测序发掘黄鳝生长差异基因及其调控通路,初步阐明造成黄鳝生长差异的基因调控机制,以期为促进黄鳝工业化养殖提供理论依据。

    实验鳝苗取自江西农业大学水产基地,均为饲养在相同水箱中、同一亲本黄鳝当年产卵孵化的鳝苗;黄鳝养殖用的配合饲料大宗原料购买于南昌大佑农生物技术有限公司。养殖1年后,从同一水箱中取个体差异显著的健康黄鳝,个体较大的为实验组 (AEG),个体较小的为对照组 (ACG),具体数据见表1,取其肝脏液氮冷冻后放−80 ℃冰箱备用。

    表  1  样品信息
    Table  1.  Information of samples
    样品Sample体质量Body mass/g全长Total length/cm
    实验组Treatment group AEG1 42.27 35.8
    AEG2 44.10 36.5
    AEG3 35.78 34.6
    对照组Control group ACG1 6.30 19.8
    ACG2 6.61 18.9
    ACG3 6.12 20.2
    下载: 导出CSV 
    | 显示表格

    取50~100 mg肝脏组织,加入一定比例的TRIzol,在低温下迅速匀浆,室温放置5 min,使其充分裂解,经氯仿抽提、异丙醇沉淀、75%乙醇洗涤后,室温开盖晾约5 min,用适量RNase-free水溶解RNA后于−80 ℃保存备用。使用1%的琼脂糖凝胶电泳检测其完整性。

    总RNA提取以后,使用美国纽英伦生物技术公司试剂盒 (#E7530) 进行cDNA建库。首先用带有Oligo (dT) 磁珠富集真核生物的mRNA加入fragmentation buffer将mRNA打断成短片段并用六碱基随机引物合成cDNA第一条链。然后加入缓冲液、DNA聚合酶Ⅰ、dNTP、RNase和缓冲溶液合成cDNA第二条链。使用快速PCR抽提试剂盒,尾端修复,纯化双链cDNA片段,并引入单碱基“A”使其与IIIumina测序接头链接。最后进行PCR扩增,通过琼脂糖凝胶电泳分离连接产物,PCR扩增富集目标片段。文库构建完成后,使用Agilent 2100 Bioanalyzer对文库进行检测以及使用ABI StepOnePlus Real-Time PCR System对文库浓度进行定量检测,合格后用IIIumina HiSeq TM 4000对cDNA文库进行测序。测序服务由北京博云华康基因科技有限公司提供。

    转录组测序后原始数据Raw reads含有低质量的reads,经过筛选过滤得到高质量的Clean reads。使用HISAT程序将得到的Clean reads比对到参考基因组上。利用DEGseq差异分析软件包进行差异基因筛选,首先计算差异倍数 (Fold change,FC) log2值和P,只有同时符合log2绝对值大于2和P的绝对值小于0.001的基因才被确定为显著差异基因 (Differentially expressed gene, DEG),然后对差异基因进行KEGG通路和GO功能富集性分析。

    对于筛选到的7个与生长相关的差异表达基因进行qPCR验证。引物设计见表2,18 S作为内参基因。荧光定量PCR反应条件为95预变性90 s,95 ℃变性5 s,60 ℃退火15 s,72 ℃延伸20 s,共40个循环。每个样品重复3次,实验数据按照Livak的2−△△Ct方法处理。

    表  2  荧光定量PCR引物信息
    Table  2.  Information of primers used for qPCR
    基因GeneqPCR引物序列 (5'–3')Primer sequence of qPCR退火温度Annealing temperature/℃片段大小Fragment size/bp
    col1α1 F:AGTTGTTTGCGGACCGAGAT 60.0 110
    R:GCAATCTGGCATTTCCTCACA 59.2
    nkx6.1 F:GGACAAAGATGGGAAACGAAA 56.7 96
    R:GCCAGGTATTTGGTCTGTTCA 58.2
    nnos F:CTATCAGTCTGGATGCCACAAC 58.8 115
    R:CAGAGCCCAACAGAAACATTAG 57.3
    plexina4 F:TGCTGAGAACCCTGAGTGGATA 60.6 159
    R:TAGCATTTGCGGTTGTCTTCAT 58.9
    pcgf1 F:CAGCCCTTACTCAACCTCAAA 57.9 167
    R:GCATCTGGCACAGCATCTACG 61.7
    igfbp1 F:CAGAGAGCCTTGGAAAAGATTG 57.3 171
    R:CTTGCCGTTCCAGGAGTGT 59.9
    h3.3 F:ATTTTGAGTTGCGGCGATTA 56.4 181
    R:GTAACGATGGGGCTTCTTCAC 59.0
    18S F:GTGGAGCGATTTGTCTGGTTA 57.8 162
    R:CGGACATCTAAGGGCATCAC 57.7
    下载: 导出CSV 
    | 显示表格

    本实验共测6个样本,平均每个样本产出约41 000 000 Clean reads,约合6.0 Gb数据量 (表3)。用Q20 (单个碱基的测序错误率低于1%) 和Q30 (单个碱基的测序错误率低于0.1%) 对这些数据的质量进行检测,Q20 和Q30值分别为98.02%和94.75%以上,说明测序质量很好。将测序结果与参考基因组进行比对,匹配率在74.32%以上。将匹配到参考基因的Reads组装成基因,共得到19 149个基因。以长势差的ACG组为对照,使用RSEM计算基因与转录表达水平,分析发现差异表达基因有598个,其中有303个基因上调,295个下调。差异表达基因的火山图和统计图分别见图1-a图1-b

    表  3  测序结果统计
    Table  3.  Statistics of sequencing results
    样品Sample过滤后ReadsFiltered ReadsQ20/%Q30/%GC含量GC content/%匹配率Matching ratio/%
    对照组Control groupACG14120116898.0394.7546.7875.48
    ACG24104770498.0494.7746.6175.48
    ACG34092654098.0294.7848.1275.98
    实验组Treatment groupAEG14100432498.0994.9646.8774.32
    AEG24098404698.1295.0047.0974.44
    AEG34099936698.0794.8847.3375.21
    下载: 导出CSV 
    | 显示表格
    图  1  差异基因火山和统计图
    红色代表上调,蓝色代表下调,黑色代表无差异,每个点代表一个基因
    Figure  1.  Volcano-plots (a) and statistics (b) of differentially expressed genes
    Up-regulated genes are in red; down-regulated ones are in blue; black represents no difference; each dot represents a gene.

    将598个差异表达基因进行GO分析 (图2),这些基因分属于生物过程、细胞组分和分子功能三大类下的42个分支。生物过程组中含有19个分支,其中单生物体过程、细胞过程和代谢过程占比最高分别为19.57%、19.23%和15.72%;细胞组分组中有15个分支,其中细胞、细胞部分和膜分占比最高,分别为11.87%、11.87%和10.87%;在分子功能分组中有8个分支,其中结合、催化活性和运输活性占比最高,分别为16.22%、13.38%和3.51%。

    图  2  差异表达基因GO功能分类图
    1. 单生物体过程;2. 细胞过程;3. 代谢过程;4. 生物调节;5. 生物过程调节;6. 定位;7. 刺激应答;8. 发展过程;9. 多细胞生物体过程;10. 生物体细胞组成或起源;11. 信号传导;12. 正调节生物过程;13. 负调节生物过程;14;生长;15. 运动;16. 免疫系统过程;17. 行为;18多生物体过程;19. 节律过程;20. 细胞;21. 细胞部分;22. 膜;23. 膜部分;24. 细胞器;25. 复杂大分子;26. 细胞器部分;27. 细胞外区域;28. 细胞外区域部分;29. 膜封闭腔;30. 细胞连接;31. 细胞外模型;32. 超分子纤维;33. 突触;34. 突触部分;35. 结合;36. 催化活性;37. 运输活性;38. 分子功能调节;39. 信号传感活性;40. 结构分子活性;41. 分子传感活性;42. 核苷酸结合转录因子活性
    Figure  2.  GO functional classification map of differentially expressed genes
    1. Single-organism process; 2. Cellular process; 3. Metabolic process; 4. Biological regulation; 5. Regulation of biological regulation; 6. Localization; 7. Response to stimulus; 8. Developmental process; 9. Multicellular organismal process; 10. Cell compent organization or biogenesis; 11. Signaling; 12. Positive regulation of biological process; 13. Negative regulation of biological process; 14. Gowth; 15. Locomotion; 16. Immune system process; 17. Behavior; 18. Multi-organism process; 19. Rhythmic process; 20. Cell; 21. Cell part; 22. Membrane; 23. Membrane part; 24. Organelle; 25. Macromolecular complex; 26. Organelle part; 27. Extracellular region; 28. Extracellular region part; 29. Membrane-enclosed lumen; 30. Cell junction; 31. Extracellular matrix; 32. Supramolecular fiber; 33. Synapse; 34. Synapse part; 35. Binding; 36. Catalytic activity; 37. Transporter activity; 38. Molecular function regulation; 39. Signal transducer activity; 40. Structural molecule activity; 41. Molecular transducer activity; 42. Nucleic acid binding transcription factor activity

    将差异表达基因序列与KEGG数据库中的数据进行BlastX比对注释,结果显示598个差异表达基因分布在262条KEGG通路中,其中显著富集的通路有38条 (P<0.05),在这38条通路中有11条代谢相关通路富集效果极显著 (P<0.01)。

    GO功能注释分析发现,与生长相关的差异表达基因有7个,分别为Ⅰ型胶原α1 (collagen typeⅠalpha 1,col1α1)、转录因子nkx6.1 (NK6 homebox 1)、神经性一氧化氮合酶 (Neuronal nitric oxide synthasennos)、神经丛蛋白家族A4 (plexina4)、类胰岛素生长因子结合蛋白-1 (Insulin-like growth factors binding protein 1, igfbp1)、多梳环指蛋白1 (Polycomb group RING finger protein 1, pcgf1) 和组蛋白3.3 (histone 3.3,h3.3)。这7个基因中,除col1α1和nkx6.1基因显著下调外,其余5个基因明显上调。这7个基因分属不同的通路 (表4),这些通路分属于4种类型,其中col1α1、nkx6.1和h3.3基因所在通路属于人类疾病,nnos基因所在通路属于新陈代谢,plexina4基因所在通路属于环境信息加工,igfbp1和pcgf1基因所在通路属于细胞过程。这些通路中只有col1α1所在的利什曼病通路和nnos所在的精氨酸和脯氨酸代谢通路是显著富集的 (P<0.05),其他5个差异基因所在通路富集效果不显著 (P>0.05)。

    表  4  生长相关差异表达基因
    Table  4.  Growth-related differentially expressed genes
    基因名称Gene name基因IDGene ID差异倍数Fold change所属通路Belonged pathway通路IDPathway ID
    col1α1 109951101 −2.5 利什曼病 Ko05140
    nkx6.1 109952563 −3.3 青少年成熟性糖尿病 Ko04950
    nnos 109959109 4.0 精氨酸和脯氨酸代谢 Ko00330
    plexina4 109959815 3.4 细胞黏附分子 Ko04514
    igfbp1 109961253 2.3 p53信号通路 Ko04115
    pcgf1 109972045 3.3 干细胞潜能调节通路 Ko04550
    h3.3 109974907 2.2 癌症的转录失调 Ko05202
    下载: 导出CSV 
    | 显示表格

    荧光定量PCR结果表明 (图3),所选的7个差异表达基因与转录组测序结果表达趋势一致。col1α1和nkx6.1这两个基因表达量均为对照组高于实验组 (实验组为1,对照组分别为3.807和1.725),nnosplexina4、igfbp1、pcgf1和h3.3这5个基因表达量则是实验组高于对照组 (实验组为1,对照组分别为0.497、0.511、0.012、0.872和0.168)。通过荧光定量PCR验证,表明转录组测序结果可靠。

    图  3  差异表达基因qPCR验证图
    Figure  3.  Validation of differentially expressed genes by qPCR

    生长分化是鱼类生长过程中的普遍现象。有观点认为,摄食不足会导致生长分化的产生[8],但杨帆等[9]研究表明,饱食投喂且增加投喂频率并不能改善黄鳝的生长分化现象。为了解黄鳝生长分化的调节机制,本研究进行了转录组测序。本次测序共得到19 149个基因,其中差异表达基因有598个,303个基因上调,295个下调。对这些差异基因进行KEGG通路分析发现,富集的KEGG通路多数与代谢相关。本研究发现,差异表达基因显著富集的38条通路中有15条与代谢有关,其中与脂肪代谢相关的有8条。在前20条富集极显著的通路中,更是有11条与代谢相关,其中影响最显著的代谢通路为类固醇生物合成。

    GO功能注释分析发现,与生长相关的差异表达基因有7个,分别为col1α1、nkx6.1、nnosplexina4、igfbp1、pcgf1和h3.3。col1α1基因显著下调,说明在ACG组的黄鳝中col1α1表达水平显著高于AEG组。Ⅰ型胶原 (COL1) 是由α1和α2两条肽链组成的,col1α1基因的过度表达使α1和α2链的比例发生变化,当比例超过2∶1时会导致骨骼密度、骨骼结构和骨骼质量发生变化从而引起骨质疏松性骨折进而影响生物体正常发育[10]。转录因子nkx6.1基因最早在胰腺β细胞中发现,与胰腺的发育相关[11]。研究发现nkx6.1基因参与调控胰腺β细胞的二次分化,胰腺β细胞能够分泌胰岛素,nkx6.1基因的突变会导致胰腺β细胞无法形成,影响胰岛素的分泌。对鸡胚胎发育的研究发现,过量表达的nkx6.1基因在胚胎发育早期能够促进脊髓Olig2的表达,但在晚期反而会抑制Olig2的表达,这说明nkx6.1在不同时期功能会发生变化。本实验中,ACG组的nkx6.1表达量明显高于AEG组,说明ACG组黄鳝过量表达nkx6.1可能会抑制胰腺β细胞的二次分化,减少胰岛素的分泌,从而影响黄鳝的生长发育。nnos存在于神经元和神经纤维中,其主要的功能就是在细胞间传递信息,具有传递和调节的作用[12]。正常状态下nnos可以产生少量NO维持细胞的生理活动。NO参与多种生理活动调节,对于生长因子增殖、T细胞活化、神经发育和神经再生都有促进作用。除此之外,对美国红鱼 (Sciaenops ocellatus) 添加维生素C的实验表明,nnos基因可以像inos一样被诱导表达参与机体的免疫应答[13]。本实验中nnos表达差异倍数最大 (4倍),NO含量水平的降低会导致胰岛素抗体的出现,使胰岛素不能正常发挥作用,最终导致黄鳝生长发育受阻。plexina4是神经丛蛋白a家族 (a1—a4) 最晚被发现的一个,plexina4对于中枢神经和外周神经的修复和再生具有重要作用[14]plexina4在脑信号蛋白3的信号传递中,起到重要的转导作用[15]。作为膜结合脑信号蛋白6A和6B的受体和信号转导因子,plexina4对于皮质脊髓束和海马组织中苔藓纤维的形成和发育至关重要[16]。有研究发现,在视觉神经和运动神经的发育信号通路中plexina4起到了指导作用,另外对于通路的维持和再生也具有重要作用[17]。本实验中,ACG组plexina4的表达量降低,有可能造成黄鳝神经系统发育不完全从而导致生长受阻。

    类胰岛素生长因子结合蛋白 (Insulin-like growth factors binding protein, Igfbps) 负责保护Igf防止其降解并调节它的生物活性[18]。Igfbps家族目前已知含有6种,从Igfbp1至Igfbp6,对硬骨鱼类的研究发现,不同组织至少含有1种Igfbp而在鱼类血液中至少含有3种主要的Igfbps[19]。在鱼类胚胎发育的各个时期都能检测到igfbp1的表达,孵化后主要在肝脏中表达。无论是成鱼还是胚胎,低氧诱导都可以使igfbp1的表达量显著增加。Igfbp1的主要功能是运输类胰岛素生长因子1 (Igf1)。Igfbp1通过与Igf1结合,调控Igf1的生物学功能:例如,介导Igf与受体之间的亲和力;控制Igf1的运输与代谢;延长Igf1的半衰期;决定Igf1的细胞通透率;调节Igf1的作用位点等。除此之外,Igfbp1自身也具有其他不依赖Igf1的生物学功能,包括控制细胞增殖、抑制机体代谢、参与肿瘤抑制、诱导细胞凋亡和促进血糖升高等。相关研究还发现许多内分泌疾病 (甲亢、糖尿病和性腺发育障碍等) 与Igfbp1有着重要联系[20],虽然对其作用机制尚不了解,但研究发现类固醇激素 (胰岛素、促肾上腺激素和雌激素等) 对igfbp1的表达有调控作用。本实验中igfbp1在ACG组中表达量显著降低,这可能造成Igf1不能完全发挥其生物学功能,从而导致黄鳝发育受阻。Pcgf1是多梳蛋白家族的一种,属于多梳抑制复合体 (Polycomb repression complex, Pcr) Pcr1。Pcgf1在神经系统高度表达因此又称为神经系统多梳蛋白1 (Nervous system polycomb 1, Nspc1)。Pcgf1蛋白在哺乳动物中可以对细胞周期进行调控,对造血干细胞的增殖和分化具有重要作用[21]。Pcr1环指部分的亲水性表面具有E3泛素连接酶活性,它可以通过改变染色质的状态来抑制基因表达[22]。有研究表明pcgf1基因的敲除会使细胞增殖能力下降[23],本实验中ACG组的pcgf1表达量降低可能导致了造血干细胞的增殖下降,从而影响黄鳝的生长发育。组蛋白H3的变体H3.3是一种重要的母源因子,能在受精后替换精子中的鱼精蛋白,参与雄性原核的重编程[24]。母源H3.3会重新激活细胞核的多潜能基因oct4,敲除h3.3后关键的多潜能基因转录水平降低,体细胞核不能被完全重编程导致胚胎不能正常发育[25-26]。注入外源h3.3 mRNA可以弥补这种缺陷。h3.3被认为是转录活性的标志[27]h3.3能够促进基因的转录表达,维持基因组稳定,保证rDNA的转录,促进rDNA的表达,重编码供体细胞核使其成为具有全能性的胚胎[28]。本实验中,ACG组h3.3表达量的降低可能导致了一些关键生长基因 (例如igfbp1、nnosplexina4和pcgf1等) 的转录表达降低,从而影响了黄鳝的生长发育。

    本文通过对生长差异显著的黄鳝进行转录组测序分析,找到了7个与黄鳝生长相关的差异表达基因,这些基因与黄鳝的神经系统、代谢系统和内分泌系统等有密切关系。结合KEGG通路分析,发现col1α1所在的利什曼病通路和nnos所在的精氨酸和脯氨酸代谢通路富集显著,说明其对黄鳝生长具有重要影响,但具体调节机制还需要进一步研究。在7个生长相关的差异表达基因中h3.3与胚胎发育相关,它既能维持基因组稳定,又能保证DNA的正确转录,还能激活多潜能基因重新编程体细胞;由于h3.3具有调节基因转录表达的功能,所以这些差异表达基因的出现是否与h3.3基因 (ACG组) 的表达下调相关还需进一步研究。

  • 图  1   差异基因火山和统计图

    红色代表上调,蓝色代表下调,黑色代表无差异,每个点代表一个基因

    Figure  1.   Volcano-plots (a) and statistics (b) of differentially expressed genes

    Up-regulated genes are in red; down-regulated ones are in blue; black represents no difference; each dot represents a gene.

    图  2   差异表达基因GO功能分类图

    1. 单生物体过程;2. 细胞过程;3. 代谢过程;4. 生物调节;5. 生物过程调节;6. 定位;7. 刺激应答;8. 发展过程;9. 多细胞生物体过程;10. 生物体细胞组成或起源;11. 信号传导;12. 正调节生物过程;13. 负调节生物过程;14;生长;15. 运动;16. 免疫系统过程;17. 行为;18多生物体过程;19. 节律过程;20. 细胞;21. 细胞部分;22. 膜;23. 膜部分;24. 细胞器;25. 复杂大分子;26. 细胞器部分;27. 细胞外区域;28. 细胞外区域部分;29. 膜封闭腔;30. 细胞连接;31. 细胞外模型;32. 超分子纤维;33. 突触;34. 突触部分;35. 结合;36. 催化活性;37. 运输活性;38. 分子功能调节;39. 信号传感活性;40. 结构分子活性;41. 分子传感活性;42. 核苷酸结合转录因子活性

    Figure  2.   GO functional classification map of differentially expressed genes

    1. Single-organism process; 2. Cellular process; 3. Metabolic process; 4. Biological regulation; 5. Regulation of biological regulation; 6. Localization; 7. Response to stimulus; 8. Developmental process; 9. Multicellular organismal process; 10. Cell compent organization or biogenesis; 11. Signaling; 12. Positive regulation of biological process; 13. Negative regulation of biological process; 14. Gowth; 15. Locomotion; 16. Immune system process; 17. Behavior; 18. Multi-organism process; 19. Rhythmic process; 20. Cell; 21. Cell part; 22. Membrane; 23. Membrane part; 24. Organelle; 25. Macromolecular complex; 26. Organelle part; 27. Extracellular region; 28. Extracellular region part; 29. Membrane-enclosed lumen; 30. Cell junction; 31. Extracellular matrix; 32. Supramolecular fiber; 33. Synapse; 34. Synapse part; 35. Binding; 36. Catalytic activity; 37. Transporter activity; 38. Molecular function regulation; 39. Signal transducer activity; 40. Structural molecule activity; 41. Molecular transducer activity; 42. Nucleic acid binding transcription factor activity

    图  3   差异表达基因qPCR验证图

    Figure  3.   Validation of differentially expressed genes by qPCR

    表  1   样品信息

    Table  1   Information of samples

    样品Sample体质量Body mass/g全长Total length/cm
    实验组Treatment group AEG1 42.27 35.8
    AEG2 44.10 36.5
    AEG3 35.78 34.6
    对照组Control group ACG1 6.30 19.8
    ACG2 6.61 18.9
    ACG3 6.12 20.2
    下载: 导出CSV

    表  2   荧光定量PCR引物信息

    Table  2   Information of primers used for qPCR

    基因GeneqPCR引物序列 (5'–3')Primer sequence of qPCR退火温度Annealing temperature/℃片段大小Fragment size/bp
    col1α1 F:AGTTGTTTGCGGACCGAGAT 60.0 110
    R:GCAATCTGGCATTTCCTCACA 59.2
    nkx6.1 F:GGACAAAGATGGGAAACGAAA 56.7 96
    R:GCCAGGTATTTGGTCTGTTCA 58.2
    nnos F:CTATCAGTCTGGATGCCACAAC 58.8 115
    R:CAGAGCCCAACAGAAACATTAG 57.3
    plexina4 F:TGCTGAGAACCCTGAGTGGATA 60.6 159
    R:TAGCATTTGCGGTTGTCTTCAT 58.9
    pcgf1 F:CAGCCCTTACTCAACCTCAAA 57.9 167
    R:GCATCTGGCACAGCATCTACG 61.7
    igfbp1 F:CAGAGAGCCTTGGAAAAGATTG 57.3 171
    R:CTTGCCGTTCCAGGAGTGT 59.9
    h3.3 F:ATTTTGAGTTGCGGCGATTA 56.4 181
    R:GTAACGATGGGGCTTCTTCAC 59.0
    18S F:GTGGAGCGATTTGTCTGGTTA 57.8 162
    R:CGGACATCTAAGGGCATCAC 57.7
    下载: 导出CSV

    表  3   测序结果统计

    Table  3   Statistics of sequencing results

    样品Sample过滤后ReadsFiltered ReadsQ20/%Q30/%GC含量GC content/%匹配率Matching ratio/%
    对照组Control groupACG14120116898.0394.7546.7875.48
    ACG24104770498.0494.7746.6175.48
    ACG34092654098.0294.7848.1275.98
    实验组Treatment groupAEG14100432498.0994.9646.8774.32
    AEG24098404698.1295.0047.0974.44
    AEG34099936698.0794.8847.3375.21
    下载: 导出CSV

    表  4   生长相关差异表达基因

    Table  4   Growth-related differentially expressed genes

    基因名称Gene name基因IDGene ID差异倍数Fold change所属通路Belonged pathway通路IDPathway ID
    col1α1 109951101 −2.5 利什曼病 Ko05140
    nkx6.1 109952563 −3.3 青少年成熟性糖尿病 Ko04950
    nnos 109959109 4.0 精氨酸和脯氨酸代谢 Ko00330
    plexina4 109959815 3.4 细胞黏附分子 Ko04514
    igfbp1 109961253 2.3 p53信号通路 Ko04115
    pcgf1 109972045 3.3 干细胞潜能调节通路 Ko04550
    h3.3 109974907 2.2 癌症的转录失调 Ko05202
    下载: 导出CSV
  • [1] 周秋白, 朱长生, 吴华东, 等. 饲料中不同脂肪源对黄鳝生长和组织中脂肪酸含量的影响[J]. 水生生物学报, 2011, 35(2): 246-255.
    [2]

    VRTÍLEK M, ŽÁK J, POLAČIK M, et al. Rapid growth and large body size in annual fish populations are compromised by density-dependent regulation[J]. J Fish Biol, 2019, 95(2): 673-678. doi: 10.1111/jfb.14052

    [3]

    GAO Z, LUO W, LIU H, et al. Transcriptome analysis and SSR/SNP markers information of the blunt snout bream (Megalobrama amblycephala)[J]. PLoS One, 2012, 7(8): 1-10.

    [4]

    RIBAS L, ROBLEDO D, GÓMEZTATO A, et al. Comprehensive transcriptomic analysis of the process of gonadal sex differentiation in the turbot (Scophthalmus maximus)[J]. Mol Cell Endocrinol, 2016, 15(422): 132-149.

    [5]

    ROBLEDO D, HERMIDA M, RUBIOLO J A, et al. Integrating genomic resources of flatfish (Pleuronectiformes) to boost aquaculture production[J]. Comp Biochem Phys D, 2016, 3(21): 41-55.

    [6]

    YIN F, GAO Q, TANG B, et al. Transcriptome and analysis on the complement and coagulation cascades pathway of large yellow croaker (Larimichthys crocea) to ciliate ectoparasite Cryptocaryon irritans infection[J]. Fish Shellfish Immun, 2016, 50: 127-141.

    [7]

    DU X, LI Y, LI D, et al. Transcriptome profiling of spleen provides insights into the antiviral mechanism in Schizothorax prenanti, after poly (I: C) challenge[J]. Fish Shellfish Immun, 2017, 62: 13-23.

    [8]

    MCCARTHY I D, CARTER C G, HOULIHAN D F. The effect of feeding hierarchy on individual variability in daily feeding of rainbow trout, Oncorhynchus mykiss (Walbaum)[J]. J Fish Biol, 1992, 41(2): 257-263. doi: 10.1111/j.1095-8649.1992.tb02655.x

    [9] 杨帆, 张世萍, 韩凯佳, 等. 投喂频率对黄鳝幼鱼摄食、生长及饵料利用效率的影响[J]. 淡水渔业, 2011, 41(3): 50-54. doi: 10.3969/j.issn.1000-6907.2011.03.009
    [10]

    XIAO S, HUANG Q, KUNG A W C. Genetics of osteoporosis in Chinese[J]. Int J Rheum Dis, 2010, 11(4): 359-365.

    [11]

    GUPTA S K, WESOLOWSKA-ANDERSEN A, RINGGAARD A K, et al. NKX6.1, induced pluripotent stem cell reporter lines for isolation and analysis of functionally relevant neuronal and pancreas populations[J]. Stem Cell Res, 2018, 29: 220-231. doi: 10.1016/j.scr.2018.04.010

    [12]

    YU S Y, ZHANG M, LUO J, et al. Curcumin ameliorates memory deficits via neuronal nitric oxide synthase in aged mice[J]. Prog Neuro-Psychoph, 2013, 45(4): 47-53.

    [13] 周立斌, 白茹, 马细兰, 等. 饲料维生素C对美国红鱼 (Sciaenops ocellatus) 组织NO含量、NOS活力和nNOS基因表达的影响[J]. 海洋与湖沼, 2014, 45(5): 1122-1126. doi: 10.11693/hyhz20140700193
    [14]

    GUTEKUNST C A, STEWART E N, FRANZ C K, et al. Plexina4 distribution in the adult rat spinal cord and dorsal root ganglia[J]. J Chem Neuroanat, 2012, 44(1): 1-13. doi: 10.1016/j.jchemneu.2012.03.002

    [15]

    SMOLKIN T, NIR-ZVI I, DUVSHANI N, et al. Complexes of plexin-A4 and plexin-D1 convey semaphorin-3C signals to induce cytoskeletal collapse in the absence of neuropilins[J]. J Cell Sci, 2018, 131(9): 5733-5783.

    [16]

    TAN C, LU N N, WANG C K, et al. Endothelium-derived semaphorin 3G regulates hippocampal synaptic structure and plasticity via neuropilin-2/PlexinA4[J]. Neuron, 2019, 101(5): 920-937. doi: 10.1016/j.neuron.2018.12.036

    [17]

    EMERSON S E, LIGHT S E, EBERT A M. Neuronal expression patterns of the PlexinA family during zebrafish development[J]. Gene Expr Patterns, 2018, 27: 56-66. doi: 10.1016/j.gep.2017.10.007

    [18]

    WANG J, LI Y C, DENG M, et al. Serum insulin-like growth factor 1 and its binding protein 3 as prognostic factors for the incidence, progression, and outcome of hepatocellular carcinoma: a systematic review and meta-analysis[J]. Oncotarget, 2017, 8(46): 81098-81108.

    [19]

    SHIMIZU M, DICKHOFF W W. Circulating insulin-like growth factor binding proteins in fish: their identities and physiological regulation[J]. Gen Comp Endocr, 2017, 252: 150-161. doi: 10.1016/j.ygcen.2017.08.002

    [20]

    FENG X, LIN J, XING S, et al. Higher IGFBP-1 to IGF-1 serum ratio predicts unfavourable survival in patients with nasopharyngeal carcinoma[J]. BMC Cancer, 2017, 17(1): 90-100. doi: 10.1186/s12885-017-3068-0

    [21]

    YUKIYA Y, AKIHIRO A, NOBUHIKO E. Clarifying the impact of polycomb complex component disruption in human cancers[J]. Mol Cancer Res, 2014, 12(4): 479-484. doi: 10.1158/1541-7786.MCR-13-0596

    [22]

    SHUKLA S, ZHAO Q, YING W, et al. Small molecule inhibitors of ring1B-Bmi1 E3 ligase target polycomb repressive complex 1 activity and regulate cell proliferation[J]. Cancer Res, 2015, 75(15): 3520-3520.

    [23] 关珊丽, 王健, 谢忆, 等. 多梳蛋白CBX家族在哺乳动物中结构与功能多样性的研究进展[J]. 中国科学:生命科学, 2017, 47(8): 847-862.
    [24]

    LIN C J, KOH F, WONG P, et al. Hira-mediated H3.3 incorporation is required for DNA replication and ribosomal RNA transcription in the mouse zygote[J]. Dev Cell, 2014, 30(3): 268-279. doi: 10.1016/j.devcel.2014.06.022

    [25]

    WEN D C, BANASZYNSKI L A, ROSENWAKS Z, et al. H3.3 replacement facilitates epigenetic reprogramming of donor nuclei in somatic cell nuclear transfer embryos[J]. Nucleus, 2014, 5(5): 369-375. doi: 10.4161/nucl.36231

    [26]

    WEN D C, BANASZYNSKI L A, LIU Y, et al. Histone variant H3.3 is an essential maternal factor for oocyte reprogramming[J]. P Natl Acad Sci USA, 2014, 111(20): 7325-7330. doi: 10.1073/pnas.1406389111

    [27]

    DELBARRE E, IVANAUSKIENE K, KUNTZIGER T, et al. DAXX-dependent supply of soluble (H3.3 –H4) dimers to PML bodies pending deposition into chromatin[J]. Genome Res, 2013, 23(3): 440-451. doi: 10.1101/gr.142703.112

    [28] 黄星卫, 程香荣, 王楠, 等. 组蛋白H3变体H3.3及其在细胞重编程中的作用[J]. 遗传, 2018, 40(3): 186-196.
  • 期刊类型引用(2)

    1. 林彬彬,袁泉,田志新,潘显斌,周文宗,徐震. 基于SSA- LSTM模型的黄鳝池溶氧预测研究. 渔业现代化. 2023(01): 71-79 . 百度学术
    2. 曹晓莉,李昭林,胡毅. 低鱼粉饲料中添加牛磺酸对黄鳝生长、消化率及肠道酶活性的影响. 南方水产科学. 2021(05): 64-70 . 本站查看

    其他类型引用(2)

图(3)  /  表(4)
计量
  • 文章访问数:  5097
  • HTML全文浏览量:  2764
  • PDF下载量:  82
  • 被引次数: 4
出版历程
  • 收稿日期:  2019-09-08
  • 修回日期:  2019-11-02
  • 网络出版日期:  2019-11-29
  • 刊出日期:  2020-02-04

目录

/

返回文章
返回