摘要:
为了分析苯系物(BTEXs)联合暴露后仿刺参(Apostichopus japonicus)管足转录组差异表达基因,采用Illumina HiSeqTM 2000测序技术,对1.0 mg·L-1苯系物(B)联合暴露12 h后和对照组(C)的仿刺参管足组织分别进行转录组测序。经Trinity软件进行de novo组装,获得了145 675条unigenes。利用公共数据库进行同源比对,共注释了35 330条unigenes。对比分析苯系物联合暴露组和对照组的转录组测序结果,获得了2 418个差异表达基因(DEGs) (|Log2 Fold changes| ≥ 1且FDR ≤ 0.001),其中,上调表达和下调表达基因分别为1 049和1 369个。GOseq分析结果显示,158个DEGs显著富集在149个GO terms中,包括103个生物学过程、17个细胞组分和29个分子功能(P < 0.05);KEGG代谢通路分析结果显示994个差异基因映射到268条代谢通路,这些差异表达基因参与的生理过程与其他生物的同源基因参与的信号传导、癌症、外源性化合物的生物降解代谢等过程相类似。上述结果为在转录组水平筛选苯系物的生物标志物,解析苯系物对仿刺参毒性作用的分子机制提供了科学参考。
Abstract:
In order to identify and analyze the differentially expressed genes (DEGs), transcriptome sequencing in the tube feet of sea cucumber Apostichopus japonicus exposed to 0.0 (C), 1.0 mg·L-1 BTEXs (B) for 12 h, was performed respectively using the Illumina HiSeqTM 2000 platform. The clean reads were then de novo assembled into 145 675 unigenes, and 35 330 unigenes were annotated by a similarity search against the public databases. By comparing B and C using the criteria |Log2 Fold changes| ≥ 1 and false discovery rate (FDR) ≤ 0.001, 2 418 DEGs were identified, among which 1 049 were up-regulated and 1 369 were down-regulated. GOseq analysis revealed that for B vs. C, 158 DEGs were highly enriched in 149 GO terms, including 103 biological processes, 17 cellular components and 29 molecular function terms (P value < 0.05). Pathways associated with BTEXs challenge were also mined. The result indicated that 994 DEGs were enriched in 268 pathways, including signal transduction, cancer, xenobiotics biodegradation and metabolism pathways. The results obtained in this study could be used for the screening of BTEXSs biomarkers on the transcriptome level, and provide the reference to understand the molecular mechanism of toxic effects of BTEXs on sea cucumber A. japonicus.