刺参响应灿烂弧菌侵染的基因组DNA甲基化水平和转录组差异及其关联分析  

Genomic DNA Methylation Levels and Transcriptome Differences of Apostichopus japonicus in Response to Vibrio splendidus Infection and Their Association Analysis

作　　者：李欣容 廖梅杰[2,3] 李彬 荣小军[2,3] 畅孟阳 王印庚[2,3] 于永翔 张正 范瑞用[4] 刘清兵 LI Xinrong;LIAO Meijie;LI Bin;RONG Xiaojun;CHANG Mengyang;WANG Yingeng;YU Yongxiang;ZHANG Zheng;FAN Ruiyong;LIU Qingbing(College of Fisheries and Life Science,Shanghai Ocean University,Shanghai 201306,China;Key Laboratory ofSustainable and Development of Marine Fisheries,Ministry of Agriculture and Rural Affairs,Yellow Sea Fisheries ResearchInstitute,Chinese Academy of Fishery Sciences,Qingdao,Shandong 266071,China;Pilot National Laboratory for MarineScience and Technology(Qingdao),Laboratory for Marine Fisheries Science and Food Production Processes,Qingdao,Shandong 266071,China;Qingdao Ruizi Company,Qingdao,Shandong 266408,China)

机构地区：[1]上海海洋大学水产与生命学院,上海201306 [2]中国水产科学研究院黄海水产研究所农业农村部海洋渔业可持续发展重点实验室,山东青岛266071 [3]青岛海洋科学与技术试点国家实验室海洋渔业科学与食物产出过程功能实验室,山东青岛266071 [4]青岛瑞滋集团有限公司,山东青岛266408

出　　处：《渔业科学进展》2022年第3期176-185,共10页Progress in Fishery Sciences

基　　金：国家重点研发计划(2018YFD0901603);山东省农业良种工程课题(2020LZGC015);中国水产科学研究院中央级公益性科研院所基本科研业务费专项资金(2020TD40;2021GH05)共同资助。

摘　　要：为探讨病原菌胁迫下刺参(Apostichopus japonicus)基因组DNA甲基化水平和转录组表达的差异,本研究采用人工攻毒侵染胁迫,获得刺参化皮体壁组织,并以未攻毒组的健康体壁组织为对照,对刺参2种体壁组织进行全基因组甲基化测序(WGBS)和转录组高通量测序,解析刺参体壁基因组DNA甲基化差异,筛选响应病原胁迫的差异甲基化区域和差异表达基因。同时,通过基因组甲基化和转录组联合分析,筛选负相关关联基因,为解析刺参响应灿烂弧菌(Vibrio splendidus)侵染的分子机理提供基础数据。结果显示,刺参对照组和侵染组基因组总甲基化水平分别为(3.59±0.04)%和(3.87±0.27)%,侵染组甲基化水平显著升高;在发生甲基化的位点中,攻毒组和对照组的m CpG占比分别为83.06%和81.91%,m CpG为最主要的甲基化形式。本研究共筛选出差异甲基化区域(DMRs)626677个,注释到23706个功能基因。转录组测序共检测到29290个基因,筛选到差异表达基因496个,其中,上调基因214个,下调基因282个。在基因组甲基化与转录组的关联分析中,筛选到180个负相关关联基因,其中,差异甲基化区域位于启动子区域的负相关基因为60个。对负相关关联基因的GO和KEGG富集分析,筛选到相关通路和LRR、hsp20和CARD等关键基因。本研究将为解析刺参抗病的表观遗传调控机制提供数据,也为刺参抗病性状选育提供科学参考。DNA methylation is an important epigenetic modification that plays a key role in gene expression regulation.In this study,two groups of sea cucumbers(Apostichopus japonicus)were prepared.One group had skin ulceration syndrome body wall(PT16S)under stress from Vibrio splendidus infection at a concentration of 1×106 CFU/mL(LD50);the other group had a healthy body wall(PT10H).Genomic DNA methylation levels and gene expression differences between the two groups were detected using whole-genome bisulfite sequencing(WGBS)and transcriptome sequencing.The key gene ontology(GO)terms and KEGG terms engaged in the immune response were selected using association analysis.The results showed that the total methylation levels of the A.japonicus genome of PT10H and PT16S were(3.59±0.04)%and(3.87±0.27)%,respectively.The methylation levels of the A.japonicus genome under pathogen challenge significantly increased;m CpG accounted for 83.06%and 81.91%of all the methylated sites in PT16S and PT10H,respectively,indicating that m CpG was the most important methylation form in the sea cucumber.A total of 626677 differentially methylated regions(DMRs)were screened and annotated into 23706 functional genes.A total of 496 differentially expressed genes were screened,including 214 up-regulated and 282 down-regulated genes.A total of 180 negatively correlated genes were isolated using association analysis between genomic methylation and transcriptome,of which 60 genes had DMRs located in promoter regions.Based on GO and KEGG enrichment of the 180negatively correlated genes,key genes such as LRR,Hsp20,and CARD were selected to play a critical role in the immune response.The results would provide primary data for the epigenetic regulatory mechanisms of A.japonicus and provide a theoretical reference for A.japonicus breeding.

关 键 词：刺参 DNA甲基化 转录组 关联分析 

分 类 号：S968.9[农业科学—水产养殖]