基于RNA-Seq高通量测序技术的大口黑鲈转录组分析

黄勇, 龚望宝, 陈海刚, 熊建利, 孙西红

黄勇, 龚望宝, 陈海刚, 熊建利, 孙西红. 基于RNA-Seq高通量测序技术的大口黑鲈转录组分析[J]. 南方水产科学, 2019, 15(1): 106-112. DOI: 10.12131/20180066
引用本文: 黄勇, 龚望宝, 陈海刚, 熊建利, 孙西红. 基于RNA-Seq高通量测序技术的大口黑鲈转录组分析[J]. 南方水产科学, 2019, 15(1): 106-112. DOI: 10.12131/20180066
HUANG Yong, GONG Wangbao, CHEN Haigang, XIONG Jianli, SUN Xihong. Sequencing and bioinformatic analysis for transcriptome of Micropterus salmoides based on RNA-seq[J]. South China Fisheries Science, 2019, 15(1): 106-112. DOI: 10.12131/20180066
Citation: HUANG Yong, GONG Wangbao, CHEN Haigang, XIONG Jianli, SUN Xihong. Sequencing and bioinformatic analysis for transcriptome of Micropterus salmoides based on RNA-seq[J]. South China Fisheries Science, 2019, 15(1): 106-112. DOI: 10.12131/20180066

基于RNA-Seq高通量测序技术的大口黑鲈转录组分析

基金项目: 国家科技支撑计划项目 (2012BAD25B04);河南省自然科学基金项目 (162300410069);河南省高等学校重点科研项目 (16A240002);广东省渔业生态环境重点实验室开放基金 (LFE-2016-13)
详细信息
    作者简介:

    黄 勇(1979—),男,博士,副教授,从事水产经济动物非编码RNA调控研究。E-mail: huangyong1979111@126.com

  • 中图分类号: S 917.4

Sequencing and bioinformatic analysis for transcriptome of Micropterus salmoides based on RNA-seq

  • 摘要:

    文章以大口黑鲈(Micropterus salmoides)组织作为研究对象,利用RNA-seq技术进行转录本测序和数据分析,经拼接组装,最终获得35 659条unigenes,序列平均长度738 bp,序列长度中位数 (N50)为1 052 bp。另外从长度分布与GC含量等方面对unigenes进行评估,数据显示测序质量好、可信度高。使用6大数据库 (KOG、Nr、Pfam、Swiss-Prot、GO和KEGG) 注释大口黑鲈转录组unigenes,分别对应有15 832、21 279、14 524、16 973、15 024和11 185条unigenes获得注释。其中,5 617条unigenes在以上所有数据库中同时注释成功,17 253条unigenes至少被一个数据库注释。KEGG分析结果显示,获得注释的11 185条unigenes被划分到267个代谢通路中,参与信号转导通路的unigenes数量最多,共有1 349条(12.06%)。另外还检测到4 030个微卫星 (SSR)位点。通过对大口黑鲈转录组测序,获得了大量的转录组信息,为大口黑鲈的功能基因克隆、基因组学、遗传多样性分析、分子标记开发及遗传改良等研究奠定了基础。

    Abstract:

    The transcripts of Micropterus salmoides were obtained by RNA-seq technology to conduct a transcriptomic analysis. A total of 35 659 unigenes were generated by de novo assembly with an average length of 738 bp and N50 of 1 052 bp. The unigenes were assessed for length distribution and GC content. The sequencing data were of high quality and reliability. A total of 15 832, 21 279, 14 524, 16 973, 15 024 and 11 185 unigenes were annotated from the KOG, Nr, Pfam, Swiss-Prot, GO and KEGG databases, respectively, among which 5 617 were annotated in all the databases and 17 253 were annotated in at least one database. Blasted with KEGG pathway, 11 185 unigenes were annotated, belonging to 267 categories. Of all the pathways, the number of unigenes joining in the signal transduction was the most (1 349, 12.06%). Finally, 4 030 SSRs were identified. The results provide rich data to understand transcriptome information of M. salmoides and lay the foundation for further research on functional gene cloning, genomics, genetic diversity analysis, molecular marker exploitation and genetic improvement in M. salmoides.

  • 中国对虾 (Fenneropenaeus chinensis),又名中国明对虾、明虾,体型大,价格昂贵,是黄、渤海区重要的经济虾类之一。其蛋白质含量高,富含镁 (Mg)、铁 (Fe) 等多种微量元素,营养价值高[1]。中国对虾在贮运流通过程中易发生黑变,严重影响商品价值。研究表明,虾体中的酪氨酸在有氧条件下经酚氧化酶 (Phenoloxidase, PO) 催化生成L-3,4-二羟基苯丙氨酸 (Levodopa, L-DOPA),进一步生成醌类化合物,而醌类物质极不稳定,可自发聚合生成黑色素,或者与蛋白质官能团交联生成黑色素,导致黑变发生[2]。虾的防黑保鲜一直是国内外研究的热点,目前虾类保鲜除控制温度外,主要采用化学方法。生产中常用亚硫酸盐类、4-己基间苯二酚 (4-hexylresorcinol, 4-HR) 等化学试剂来延缓虾的黑变。GB 2760—2014《食品添加剂使用标准》明确规定了焦亚硫酸钠和4-HR的残留限量,即虾体内二氧化硫 (SO2) 残留量不得超过0.1 g·kg−1,4-HR残留量不得超过1 mg·kg−1。国内外学者对亚硫酸盐类[3-5]和4-HR[6-7]在虾类防黑保鲜中的应用进行了系列研究,发现在限量范围内,亚硫酸盐和4-HR的防黑效果有限,且过量使用则可能产生安全问题[8]。近年来的相关研究多聚于一些天然植物提取物,如植物多酚[9]、石榴皮提取物[10]、绿茶提取物[11]等,在抑制虾类黑变方面表现出较大的应用潜力。将这些天然植物提取物与焦亚硫酸钠、4-HR等复配使用,有助于解决单一成分防黑效果不理想的问题[12]

    本研究在GB 2760限定范围内对具有潜在防黑作用的物质进行筛选,将常用的焦亚硫酸钠、4-HR与筛选确定的天然植物提取物进行复配,通过正交试验确定了复合防黑保鲜剂的最优配比,并对其在中国对虾防黑保鲜方面的应用效果进行了验证。研究结果为虾类绿色高效防黑保鲜剂的开发与应用提供参考。

    鲜活中国对虾购于青岛市崂山区沙子口码头,体长 (19.7±1.25) cm,体质量 (62.5±10.0) g,充氧保活条件下2 h内运至实验室,经暂养后流水冲洗,然后加冰猝死,置于无菌样品袋中备用。

    石榴皮提取物 (纯度为90%)、葡萄籽提取物 (纯度为95%) (均为食品级) 购自陕西昊辰生物科技有限公司;竹叶抗氧化物 (食品级) 购自河南豫中生物科技有限公司;三氯乙酸 (TCA)、抗坏血酸、焦亚硫酸钠、硫代巴比妥酸 (TBA)、4-己基间苯二酚 (均为分析纯) 购自国药集团化学试剂有限公司。

    METTLER TOLEDO (梅特勒托利多) 分析天平 (苏州赛恩斯仪器有限公司);T18高速分散机 (德国IKA公司);SHA-B恒温振荡器 (常州国语仪器制造有限公司);5PX-25085H-II生化培养箱 (上海新苗医疗器械有限公司)。

    参照Zhang等[13]的方法提取酚氧化酶粗酶,并略作调整。将虾的头胸部经液氮冷冻后研磨成粉,按照1∶3 (mv) 的比例将头胸部粉末与缓冲液 [0.05 mol·L−1、pH 7.2磷酸盐缓冲液,含0.2% (w) Brij-35和1 mol·L–1氯化钠 (NaCl)] 混合,4 ℃下搅拌30 min,过滤,将上清液于4 ℃、8 000 r·min–1离心30 min,过滤,得到的上清液即为粗酶液。

    参照Sae-Leaw和Benjakul[14]的方法对酚氧化酶粗酶进行纯化。称取一定质量的硫酸铵固体加入到酚氧化酶粗酶液中,使其质量分数达50%,静置30 min后于8 000 r·min–1条件下离心30 min,去除上清液,收集沉淀于透析袋中透析24 h,期间更换3次透析液,透析结束后进行浓缩,8 000 r·min–1离心30 min,除去不溶物。经硫酸铵沉淀后的粗酶液上样于弱阴离子交换柱DEAE Sepharose Fast Flow (1.6 cm×25 cm),用0.05 mol∙L–1的磷酸盐 (pH 7.2) 清洗柱至A280 nm小于0.05,分别使用含0.1、0.2、0.3、0.4、0.5 mol·L–1 NaCl的磷酸盐缓冲液 (0.05 mol·L–1, pH 7.2) 进行洗脱,流速1 mL·min–1,测定各管的酚氧化酶活性,将有酚氧化酶活性和蛋白峰的管分别合并,利用超滤管对其除盐浓缩后置于−80 ℃中保存备用。

    参照Nirmal等[15]的方法并略有改动。将不同质量浓度 (1.0、3.0、5.0、7.0、10.0 g·L–1) 的抗坏血酸、竹叶抗氧化物、石榴皮提取物、葡萄籽提取物与酚氧化酶溶液按1∶1 (vv)的比例混合,室温下反应30 min后与底物L-DOPA溶液 (15 mmol·L–1,去离子水溶解) 混合,45 ℃水浴中孵育3 min,测定其475 nm处的动力学曲线。对照组以去离子水代替。试验结果以相对于酚氧化酶最大酶活的百分比表示。1个酶活力单位(U)定义为,在475 nm波长下,单位体积酶溶液 (mL) 在单位时间 (min) 内使酶促反应吸光度值增加0.001,用U·mL–1酶液表示。

    预试验发现,焦亚硫酸钠溶液质量浓度不超过15 g·L–1时,处理对虾体内的SO2残留量低于0.1 g·kg–1;4-HR溶液质量浓度不超过50 mg·L–1时,4-HR残留量不高于1 mg·L–1。基于预试验结果,并参照文献[12],选取抗坏血酸、竹叶抗氧化物、焦亚硫酸钠、4-HR,分别设置质量浓度为1.0~5.0 g·L–1、1.0~5.0 g·L–1、13.0~15.0 g·L−1、30.0~50.0 mg·L–1。采用4因素、3水平L9 (34) 正交试验 (表1) 对复合保鲜剂的配比进行优化。以中国对虾酚氧化酶活性为指标,对照组以去离子水代替,每组取3次平行样。

    表  1  复合保鲜剂配方正交试验设计
    Table  1.  Orthogonal experiment design of compound anti-blackening preservative
    水平
    Level
    A:抗坏血酸
    Ascorbic acid/(g·L–1)
    B:竹叶抗氧化物
    Bamboo leaf antioxidant/(g·L–1)
    C:焦亚硫酸钠
    Sodium metabisulfite/(g·L–1)
    D:4-HR
    4-hexylresorcinol/(mg·L–1)
    11.01.013.030.0
    23.03.014.040.0
    35.05.015.050.0
    下载: 导出CSV 
    | 显示表格

    按1.3.3得到的复合保鲜剂的最佳配比配制成溶液,虾与保鲜剂的比例为1∶4 (mv),浸泡时间为1 min,浸泡后沥水2~3 min。

    按照1.3.4的方法处理后,于4 ℃条件下贮藏,分别在第0、第24、第48、第72、第96、第120、第144 小时取样,同时对照组以去离子水浸泡后于同一条件下取样,进行黑变评分并测定相关指标。

    参考Montero等[6]的方法,由6名经培训的感官品评人员采用表2的标准对样品的黑变情况进行评分。

    表  2  中国对虾黑变评分标准
    Table  2.  Evaluation criteria for blakening of F. chinensis
    分值
    Score
    黑变情况
    Condition of blackening
    0 无黑变,即体表完全没有黑点
    0.1~1.0 轻微黑变,即体表有少量黑点或整体黑变面积较小
    1.1~2.0 中等黑变,即体表有较多的黑点或整体黑变面积稍大
    2.1~3.0 严重黑变,即体表有大量黑点或整体黑变面积大
    下载: 导出CSV 
    | 显示表格

    TBARS值的测定参照GB 5009.181—2016 《食品安全国家标准 食品中丙二醛的测定》中第二法 分光光度法,结果以mg·kg−1表示。

    TVB-N值参照GB 5009.228—2016《食品中挥发性盐基氮的测定》中的微量扩散法测定,结果以mg·100 g–1表示。

    菌落总数参照GB 4789.2—2022《食品微生物学检验 菌落总数测定》,采用平板计数,计数单位为lg CFU·g–1

    将经复合保鲜剂处理后的中国对虾于4 ℃条件下贮藏,分别在第0、第72、第144 小时取样,分析其微生物多样性。采用HTAB (Hexadecyltrimethy ammonium bromide) 法[16]从样本中提取基因组DNA,对16S rDNA的V3–V4区进行扩增,电泳检测PCR扩增产物的纯度,用QuantiFluorTM荧光计定量,将纯化的扩增产物等量混合,连接测序接头,构建测序文库,Illumina PE250上机测序 (由广州基迪奥生物科技有限公司协助完成)[17]。参照Huptas等[18]的方法筛选数据后进行聚类分析,将一致性超过95%的序列聚类为操作分类单位 (Operational taxonomic unit, OTU) 去除聚类比对过程中的嵌合体序列,得到最终有效数据,并进行OTU丰度统计分析。

    试验重复2次,每次取3次平行样,结果以“平均值±标准差 ($\bar x $±s)”表示。采用t-检验进行组间差异显著性分析,P<0.05为差异显著,P<0.01为差异极显著。使用Origin Pro 2022b软件作图。

    抗坏血酸、葡萄籽提取物、石榴皮提取物、竹叶抗氧化物对酚氧化酶的抑制作用见图1。当抗坏血酸的质量浓度为3.0 g·L–1时,可抑制90%以上的酚氧化酶活性。抗坏血酸可以通过还原黑变过程中的中间产物-醌类物质来阻断黑变进程[19],抗坏血酸还可以作为竞争性抑制剂与酚氧化酶的双铜活性中心结合,减少酶与底物的接触,进而抑制黑变[20]。在试验选定的浓度范围内,葡萄籽提取物对酚氧化酶的抑制作用与其质量浓度几乎呈线性相关,但过高浓度的葡萄籽提取物溶液颜色过深,对虾的感官品质有不良影响,因此不适宜作为复合保鲜剂的配方。石榴皮提取物对酚氧化酶的抑制作用有限,当质量浓度提高至7.0 g·L–1时,也仅能将相对酶活控制在70%左右。竹叶抗氧化物质量浓度为3.0 g·L–1时几乎可以完全抑制酚氧化酶活性,竹叶提取物由黄酮、酚类、酸类等物质构成,具有良好的抑菌和抗氧化作用[21]。因此,后续试验选择抗坏血酸和竹叶抗氧化物作为复配组分。

    图  1  不同物质对中国对虾酚氧化酶活性的抑制作用
    Figure  1.  Inhibitory effect of different substances on phenoloxidase activity of F. chinensis

    将筛选出的抗坏血酸、竹叶抗氧化物与焦亚硫酸钠、4-HR复配,以酚氧化酶活性为指标进行正交试验 (表3)。结果显示极差:RC> RD> RA> RB,说明在选定的浓度范围内,对酚氧化酶活性影响大小依次为焦亚硫酸钠>4-HR>抗坏血酸>竹叶抗氧化物。较优方案为A2B3C3D3,即复合保鲜剂的最优配比为:抗坏血酸质量浓度3.0 g·L–1、竹叶抗氧化物质量浓度5.0 g·L–1、焦亚硫酸钠质量浓度15.0 g·L–1、4-HR质量浓度50.0 mg·L–1

    表  3  L9 (34) 正交试验结果
    Table  3.  L9 (34) orthogonal test results
    序号
    No.
    A:抗坏血酸
    Ascorbic acid
    B:竹叶抗氧化物
    Bamboo leaf antioxidant
    C:焦亚硫酸钠
    Sodium metabisulfite
    D:4-HR
    4-hexylresorcinol
    酚氧化酶活性
    PO activity/(U·L−1)
    11111115.00
    2123246.67
    3132385.00
    4213324.67
    52221103.33
    6231281.67
    7312296.67
    83213110.00
    9333151.67
    K1246.67236.33306.67270.00
    K2209.67260.00285.00225.00
    K3258.34218.33123.00219.67
    k182.2278.78102.2290.00
    k269.8986.6795.0075.00
    k386.1172.7841.0073.22
    R48.6723.67162.0050.33
    下载: 导出CSV 
    | 显示表格

    中国对虾经复合保鲜剂处理后的黑变情况见图2。对照组在冷藏期间黑变评分迅速上升,贮藏第24小时,头胸和躯干部出现黑斑;第48小时黑变加重;第72小时黑斑遍布对虾全身,这与郭芳等[22]的观察结果基本一致。经复合保鲜剂处理后,中国对虾在整个贮藏期间的黑变评分均显著低于对照组 (P<0.05),尤其是前96 h仅发生了轻微黑变,即使在贮藏144 h后虾的感官品质也在可接受范围内,表明复合保鲜剂有效延缓了中国对虾黑变进程,这与其含有的抑菌成分和抗氧化成分有效抑制黑变过程中腐败菌的生长以及酚氧化酶氧化反应有关[23]。此外,徐扬等[24]和López-Caballero等[25]研究发现焦亚硫酸钠对虾类的黑变控制效果有限,仅在贮藏前期起到控制黑变的效果,贮藏后期 (4~5 d) 与未处理之间无显著性差异。然而,本研究中复合保鲜剂处理组在贮藏期内与对照组差异显著 (P<0.05),说明研制的复合保鲜剂的抑菌和抗氧化活性更强,从而抑制虾类黑变的效果更优,可有效延长虾类的货架期。

    图  2  复合保鲜剂对中国对虾冷藏期间防黑效果评价
    Figure  2.  Evaluation of anti-blackening effect of compound preservative on F. chinensis

    TBARS值可用于评估水产品的脂肪氧化程度,氧化程度越高,TBARS值越大[26]。中国对虾的初始TBARS值为0.067 mg·kg–1。冷藏过程中,两组对虾TBARS值的增长趋势相似 (图3-a)。相较对照组,处理组增长趋势较为缓慢,且在整个贮藏期间TBARS值均显著低于对照组 (P<0.05),表明复合保鲜剂能够有效延缓中国对虾脂质氧化,这与其中的抗坏血酸、竹叶抗氧化物等成分都具有良好的抗氧化作用有关[27-28],有效抑制了脂质氧化引起的对虾在贮藏过程中的品质劣变。杨峰等[29]通过对南极磷虾 (Euphausia superba) 的保鲜研究同样发现复配后的保鲜剂能够显著抑制脂质氧化,与单一抑黑剂相比效果更优。

    图  3  中国对虾冷藏期间TBARS值、TVB-N值和菌落总数变化
    Figure  3.  Changes in TBARS value, TVB-N value, and total bacterial counts during storage of F. chinensis

    TVB-N值是评价水产品鲜度的重要指标[30]。中国对虾冷藏期间TVB-N值变化见图3-b。在冷藏过程中,两组样品的TVB-N值呈增加趋势。冷藏24 h以及后续的取样点,处理组的TVB-N值均显著低于对照组 (P<0.05)。根据GB 2733—2015,对虾的TVB-N的限量值为30.0 mg·100 g–1。冷藏第144 小时,对照组TVB-N值达32.65 mg·100 g–1,而此时此处理组TVB-N值仅为24.83 mg·100 g–1。复合保鲜剂中的成分对微生物的生长有抑制作用,减少了微生物对蛋白质的降解,导致TVB-N值增加缓慢,这与Gokoglu等[31]观察不同抑黑剂处理红虾 (Aristaeomorpha foliacea) 后的TVB-N值变化结果一致。此外,Senapati等[32]研究发现经4-HR处理后凡纳滨对虾 (Litopenaeus vannamei) 的TVB-N值低于焦亚硫酸钠处理组。加之,竹叶抗氧化物中的黄酮类化合物成分以及抗坏血酸低pH的特性使其均具有抑制微生物生长的作用[33-34],表明研制的复合保鲜剂可以更有效地控制中国对虾在贮藏过程中的品质变化。

    菌落总数是评定水产品鲜度的重要指标,水产品腐败与微生物的生长繁殖密切相关[35]。如图3-c所示,对照组、处理组的初始菌落总数分别为2.69和2.45 lg CFU·g–1。在冷藏前48 h,两组对虾的菌落总数增长迅速,48 h后增长趋缓。冷藏过程中对照组的菌落总数均显著高于处理组 (P<0.05)。冷藏第120小时,对照组的菌落总数已达5.11 lg CFU·g–1,此时处理组的菌落总数仅为4.65 lg CFU·g–1,仍在一级鲜度范围内[36],说明复合保鲜剂能够在一定程度上抑制微生物的繁殖,这与复合保鲜剂中的4-HR[37]、焦亚硫酸钠[38]、竹叶抗氧化物[39]、抗坏血酸[40]4种组分均具备一定的抑菌作用有关。因此,研制的复合保鲜剂可有效延缓中国对虾的菌落总数升高,一定程度上延长了对虾的货架期。

    中国对虾经复合保鲜剂处理后,在冷藏过程中的菌群结构变化情况见图4。第0天的对照组与处理组的初始菌群略有差别,对照组中别弧菌属 (Aliivibrio)、黏着杆菌属 (Tenacibaculum)、交替假单胞菌属 (Pseudoalteromonas) 比例相对较高,而处理组则是别弧菌属、黄杆菌属 (Flavobacterium)、弧菌属 (Vibrio) 比例占优。

    图  4  中国对虾冷藏过程属水平上的细菌分布堆叠图
    注:C. 对照组;T. 处理组;C72. 对照组贮藏72 h的样品,以此类推。
    Figure  4.  Stacked map of bacteria species distribution of F. chinensis at genus level during storage
    Note: C. Control group; T. Treatment group; C72. Samples stored for 72 h in the control group, and so on.

    对照组中的别弧菌属在贮藏前期比例迅速增加,第72小时达54%,而贮藏后期比例减少;发光杆菌属 (Photobacterium) 随冷藏时间的延长占比逐渐升高,第144小时占比增至45%,成为中国对虾的优势菌属。Kuuliala等[41]在大西洋银鳕 (Gadus morhua) 中也观察到发光杆菌属在冷藏过程中比例持续增长的现象,说明发光杆菌属可能是导致水产品腐败的重要菌属,其生长代谢产物是中国对虾冷藏后期不良气味形成的主要来源。处理组中的优势菌属为致腐能力相对较弱的别弧菌属,其比例随着贮藏时间的延长持续增加,而致腐能力较强的发光杆菌属在冷藏后期明显减少,表明复合保鲜剂可有效抑制发光杆菌属的繁殖。因此,复合保鲜剂的抑菌活性可达到延缓水产品在贮藏过程中快速腐败的目的。Peng等[42]研究表明一种新型的羟基吡喃酮-氨基硫脲衍生物对虾体内的副溶血性弧菌 (V. parahaemolyticus) 表现出强烈的抑菌活性,有效延长了虾类的货架期,与本研究结果相似。

    本研究经对选取的多种具有抑黑作用的物质进行单因素试验和正交试验,最终确定复合防黑保鲜剂较优配比为:抗坏血酸质量浓度3.0 g·L–1、竹叶抗氧化物质量浓度5.0 g·L–1、焦亚硫酸钠质量浓度15.0 g·L–1、4-HR质量浓度50.0 mg·L–1。经复合保鲜剂处理后,中国对虾黑变被明显抑制,同时该保鲜剂具有抑制细菌生长、延缓脂质氧化的作用。微生物多样性分析结果表明,发光杆菌属和别弧菌属是中国对虾冷藏过程的优势菌属,而复合保鲜剂可以有效抑制发光杆菌属的生长繁殖。因此,研制的复合保鲜剂具有良好的抑黑保鲜效果,具有较好的应用前景。

  • 图  1   Unigenes GC含量分布

    Figure  1.   GC content distribution of unigenes

    图  2   转录组transcript和unigenes长度分布

    A. transcript的长度分布;B. unigenes的长度分布

    Figure  2.   Length distribution of assembled transcripts and unigenes of transcriptome

    A. size distribution of transcripts; B. size distribution of unigenes

    图  3   Unigenes GO功能注释

    1. 转录、DNA依赖性;2. 转录调控、DNA依赖性;3. 蛋白质转运;4. 凋亡;5. 细胞分化;6. 细胞黏附;7. 蛋白质水解;8. 多细胞器官发育;9. 细胞周期;10. 细胞分化;11. 胞内信号转导;12. mRNA加工;13. 信号转导;14. 小G蛋白介导信号转导;15. 转运;16. 转录正调控;17. 有丝分裂;18. 染色质修饰;19. 转录负调控;20. RNA剪切;21. 胞内蛋白质运输;22. 翻译;23. DNA修复;24. 转录;25. 内噬作用;26. 细胞核;27. 细胞浆;28. 必须膜;29. 胞液;30. 细胞质膜;31. 内质网膜;32. 核仁;33. 线粒体;34. 胞外区;35. 核质;36. 隔膜;37. 高尔基体膜;38. 细胞骨架;39. 高尔基体;40. 细胞核周区;41. ATP结合;42. 锌离子结合;43. 蛋白结合;44. DNA结合;45. 金属离子结合;46. RNA结合;47. 钙离子结合;48. 特异序列DNA结合;49. 结合;50. 蛋白质丝氨酸特异酶结合

    Figure  3.   GO functional annotation of unigenes

    1. transcription, DNA-dependent; 2. regulation of transcription, DNA-dependent; 3. protein transport; 4. apoptosis; 5. cell division; 6. cell adhesion; 7. proteolysis; 8. multicellular organismal development; 9. cell cycle; 10. cell differentiation; 11. intracellular signal transduction; 12. mRNA procession; 13. signal transduction; 14. small GTPase mediated signal transduction; 15. transport; 16. negative regulation of transcription; 17. mitosis; 18. chromatin modification; 19. positive regulation of transcription; 20. RNA splicing; 21. intracellular protein; 22. translation; 23. DNA repair; 24. transcription; 25. endocytosis; 26. nucleus; 27. cytoplasm; 28. integral to membrane; 29. cytosol; 30. plasma membrane; 31. endoplasmic reticulum membrane; 32. nuclelous; 33. mitochondrion; 34. extracellular region; 35. nucleoplasm; 36. membrane; 37. golgi membrane; 38. cytoskeleton; 39. golgi apparatus; 40. perinuclear region of cytoplasm; 41. ATP binding; 42. zinc ion binding; 43. protein binding; 44. DNA binding; 45. metal ion binding; 46. RNA binding; 47. calcium ion binding; 48. specific DNA sequence binding; 49. binding; 50. protein serine/threonine-specific kinase

    图  4   Unigenes的KOG注释

    A. RNA加工和修饰;B. 染色体结构和动力学;C. 能源生产和转换;D. 细胞周期调控-细胞分裂-染色体分离;E. 氨基酸转运和代谢;F. 氨基酸转运和代谢;G. 碳水化合物转运和代谢;H. 辅酶转运和代谢;I. 脂质转运和代谢;J. 翻译-核糖体结构-生物合成;K. 转录;L. 复制-重组-修复;M. 细胞壁-细胞膜合成;N. 细胞运动;O. 翻译后修饰-蛋白质周转-分子伴侣;P. 无机离子转运与代谢;Q. 次生代谢产物生物合成、转运和分解代谢;R. 一般功能预测;S. 功能未知;T. 信号转导机制;U. 细胞内分泌和囊泡运输;V. 防御机制;W. 胞外结构;Y. 核结构;Z. 细胞骨架

    Figure  4.   KOG annotation of unigenes

    A. RNA processing and modification; B. chromatin structure and dynamics; C. energy production and conversion; D. cell cycle control, cell division, chromosome partitioning; E. amino acid transport and metabolism; F. nucleotide transport and metabolism; G. carbohydrate transport and metabolism; H. coenzyme transport and metabolism; I. lipid transport and metabolism; J. translation, ribosomal structure and biogenesis; K. transcription; L. replication, recombination and repair; M. cell wall/membrane/envelope biogenesis; N. cell motility; O. posttranslational modification, protein turnover, chaperones; P. inorganic ion transport and metabolism; Q. secondary metabolites biosynthesis, transport and catabolism; R. general function prediction only; S. function unknown; T. signal transduction mechanisms; U. intracellular trafficking, secretion, and vesicular transport; V. defense mechanisms; W. extracellular structure; Y. nuclear structure; Z. cytoskeleton

    图  5   Unigenes的KEGG注释

    Figure  5.   KEGG annotation of unigenes

    表  1   转录组数据拼接结果统计

    Table  1   Assembly result of transcriptome data

    序列种类
    category
    总数
    total number
    N50长度/bp
    N50 length
    总长度/bp
    total length
    最大长度/bp
    maximum length
    中等长度/bp
    median length
    最小长度/bp
    minimum length
    平均长度/bp
    average length
    transcript 42 264 1 331 36 549 208 35 516 522 201 864
    unigene 35 659 1 052 26 340 493 35 516 472 201 738
     注:N50表示将转录本从长到短排序,依次累加碱基数,当累计碱基数达到转录本总碱基数的50%时的转录本的长度  Note: N50 of transcript or unigenes was calculated by ordering all sequences, then adding the lengths from longest to shortest until the summed length exceeded 50% of the total length of all sequences.
    下载: 导出CSV
  • [1]

    LI Y, ZHOU Z, TIAN M, et al. Exploring single nucleotide polymorphism (SNP), microsatellite (SSR) and differentially expressed genes in the jellyfish (Rhopilema esculentum) by transcriptome sequencing[J]. Mar Genom, 2017, 34: 31-37. doi: 10.1016/j.margen.2017.01.007

    [2]

    CARRUTHERS M, YURCHENKO A A, AUGLEY J J, et al. De novo transcriptome assembly, annotation and comparison of four ecological and evolutionary model salmonid fish species[J]. BMC Genom, 2018, 19(1): 32. doi: 10.1186/s12864-017-4379-x

    [3]

    JORGE P H, MASTROCHIRCO F A, HATA M E, et al. Genetic characterization of the fish Piaractus brachypomus by microsatellites derived from transcriptome sequencing[J]. Front Genet, 2018, 9: 46. doi: 10.3389/fgene.2018.00046

    [4] 罗辉, 叶华, 肖世俊, 等. 转录组学技术在水产动物研究中的运用[J]. 水产学报, 2015, 39(4): 598-607.
    [5]

    CUI J, XU J, ZHANG S, et al. Transcriptional profiling reveals differential gene expression of Amur ide (Leuciscus waleckii) during spawning migration[J]. Int J Mol Sci, 2015, 16(6): 13959-13972.

    [6] 邱忠营, 刘菲, 张克瑶, 等. 疣蝗转录组分析[J]. 基因组学与应用生物学, 2016, 35(8): 1989-1998.
    [7]

    SUDHAGAR A, KUMAR G, EL-MATBOULI M. Transcriptome analysis based on RNA-Seq in understanding pathogenic mechanisms of diseases and the immune system of fish: a comprehensive review[J]. Int J Mol Sci, 2018, 19(1): 1-19.

    [8]

    CAO S, ZHU L, NIE H, et al. De novo assembly, gene annotation, and marker development using Illumina paired-end transcriptome sequencing in the Crassadoma gigantea[J]. Gene, 2018, 658: 54-62. doi: 10.1016/j.gene.2018.03.019

    [9]

    JIA Z, WANG Q, WU K, et al. De novo transcriptome sequencing and comparative analysis to discover genes involved in ovarian maturity in Strongylocentrotus nudus[J]. Comp Biochem Physiol D, 2017, 23: 27-38.

    [10]

    LU T, SUN Y, MA Q, et al. De novo transcriptomic analysis and development of EST-SSR markers in the Siberian tiger (Panthera tigris altaica)[J]. Mol Genet Genom, 2016, 291(6): 2145-2157. doi: 10.1007/s00438-016-1246-4

    [11]

    RAI R, CHAUHAN S K, SINGH V V, et al. RNA-seq analysis reveals unique transcriptome signatures in systemic lupus erythematosus patients with distinct autoantibody specificities[J]. PLoS One, 2016, 11(11): e0166312. doi: 10.1371/journal.pone.0166312

    [12]

    LI S, SHEN L, SUN L, et al. Small RNA-Seq analysis reveals microRNA-regulation of the Imd pathway during Escherichia coli infection in Drosophila[J]. Dev Comp Immunol, 2017, 70: 80-87. doi: 10.1016/j.dci.2017.01.008

    [13]

    SALEM M, REXROAD C E, WANG J N, et al. Characterization of the rainbow trout transcriptome using Sanger and 454-pyrosequencing approaches[J]. BMC Genomics, 2010, 11: 564. doi: 10.1186/1471-2164-11-564

    [14]

    COPPE A, PUJOLAR J M, MAES G E, et al. Sequencing, de novo annotation and analysis of the first Anguilla anguilla transcriptome: EeelBase opens new perspectives for the study of the critically endangered european eel[J]. BMC Genomics, 2010, 11: 635. doi: 10.1186/1471-2164-11-635

    [15]

    PEREIRO P, BALSEIRO P, ROMERO A, et al. High-throughput sequence analysis of turbot (Scophthalmus maximus) transcriptome using 454-pyrosequencing for the discovery of antiviral immune genes[J]. PLoS One, 2012, 7(5): e35369. doi: 10.1371/journal.pone.0035369

    [16] 赵文, 高峰英, 石振广. 达氏鳇鱼肌肉组织转录组测序与功能分析[J]. 水产学报, 2014, 38(9): 1255-1262.
    [17]

    MU Y N, DING F, CUI P, et al. Transcriptome and expression profiling analysis revealed changes of multiple signaling pathways involved in immunity in the large yellow croaker during Aeromonas hydrophila infection[J]. BMC Genom, 2010, 11: 506. doi: 10.1186/1471-2164-11-506

    [18]

    WANG R J, SUN L Y, BAO L S, et al. Bulk segregant RNA-seq reveals expression and positional candidate genes and allele-specific expression for disease resistance against enteric septicemia of catfish[J]. BMC Genom, 2013, 14: 929. doi: 10.1186/1471-2164-14-929

    [19]

    SUTHERLAND B J, KOCZKA K W, YASUIKE M, et al. Comparative transcriptomics of Atlantic Salmo salar, chum Oncorhynchus keta and pink salmon O. gorbuscha during infections with salmon lice Lepeophtheirus salmonis[J]. BMC Genom, 2014, 15(1): 200. doi: 10.1186/1471-2164-15-200

    [20]

    DANG Y, XU X, SHEN Y, et al. Transcriptome analysis of the innate immunity-related complement system in spleen tissue of Ctenopharyngodon idella infected with Aeromonas hydrophila[J]. PLoS One, 2016, 11(7): e0157413. doi: 10.1371/journal.pone.0157413

    [21] 赵刚, 龚全, 刘亚, 等. 基于Illumina高通量测序的岩原鲤转录组分析[J]. 西南农业学报, 2016, 29(7): 1743-1749.
    [22] 许建, 赵建, 徐礼鸣, 等. 基于RNA-Seq技术的鲮转录组分析[J]. 大连海洋大学学报, 2014, 29(6): 556-560. doi: 10.3969/J.ISSN.2095-1388.2014.06.003
    [23] 全迎春, 马冬梅, 白俊杰, 等. 大口黑鲈转录组SNPs筛选及其与生长的关联分析[J]. 水生生物学报, 2016, 11(40): 294-300.
图(5)  /  表(1)
计量
  • 文章访问数:  4750
  • HTML全文浏览量:  2049
  • PDF下载量:  89
  • 被引次数: 0
出版历程
  • 收稿日期:  2018-04-09
  • 修回日期:  2018-06-14
  • 录用日期:  2018-07-24
  • 网络出版日期:  2018-12-04
  • 刊出日期:  2019-02-04

目录

/

返回文章
返回