贝莱斯芽孢杆菌LG37全基因组测序分析及无机氮代谢相关候选基因的筛选

刘广鑫, 董晏君, 赵丽娟, 邓益琴, 程长洪, 马红玲, 江建军, 冯娟, 郭志勋, 林蠡

刘广鑫, 董晏君, 赵丽娟, 邓益琴, 程长洪, 马红玲, 江建军, 冯娟, 郭志勋, 林蠡. 贝莱斯芽孢杆菌LG37全基因组测序分析及无机氮代谢相关候选基因的筛选[J]. 南方水产科学, 2022, 18(3): 57-67. DOI: 10.12131/20210149
引用本文: 刘广鑫, 董晏君, 赵丽娟, 邓益琴, 程长洪, 马红玲, 江建军, 冯娟, 郭志勋, 林蠡. 贝莱斯芽孢杆菌LG37全基因组测序分析及无机氮代谢相关候选基因的筛选[J]. 南方水产科学, 2022, 18(3): 57-67. DOI: 10.12131/20210149
LIU Guangxin, DONG Yanjun, ZHAO Lijuan, DENG Yiqin, CHENG Changhong, MA Hongling, JIANG Jianjun, FENG Juan, GUO Zhixun, LIN Li. Sequencing of whole genome of Bacillus velezensis LG37 and screening of inorganic nitrogen metabolism candidate genes[J]. South China Fisheries Science, 2022, 18(3): 57-67. DOI: 10.12131/20210149
Citation: LIU Guangxin, DONG Yanjun, ZHAO Lijuan, DENG Yiqin, CHENG Changhong, MA Hongling, JIANG Jianjun, FENG Juan, GUO Zhixun, LIN Li. Sequencing of whole genome of Bacillus velezensis LG37 and screening of inorganic nitrogen metabolism candidate genes[J]. South China Fisheries Science, 2022, 18(3): 57-67. DOI: 10.12131/20210149

贝莱斯芽孢杆菌LG37全基因组测序分析及无机氮代谢相关候选基因的筛选

基金项目: 中国水产科学研究院南海水产研究所中央级公益性科研院所基本科研业务费专项资金资助 (2020TS04);广东省重点领域研发计划项目 (2019B020215001);中国-东盟海上合作基金 (CAMC-2018F);广东省海洋与渔业厅项目 (D21822202);财政部和农业农村部国家现代农业产业技术体系资助 (CARS-48)
详细信息
    作者简介:

    刘广鑫 (1988—),男,助理研究员,博士,从事渔业生物病害防治研究。E-mail: Guangxin_liu1988@163.com

    通讯作者:

    郭志勋 (1970—),男,研究员,博士,从事渔业生物病害防治研究。E-mail: guozhixun1@163.com

    林 蠡 (1970—),男,教授,博士,从事渔业生物病害防治研究。E-mail: linli@zhku.edu.cn

  • 中图分类号: S 949

Sequencing of whole genome of Bacillus velezensis LG37 and screening of inorganic nitrogen metabolism candidate genes

  • 摘要: 前期研究发现贝莱斯芽孢杆菌 (Bacillus velezensis) LG37可高效同化无机氮,但其机理尚不清楚。为解读其高效同化无机氮的机理,结合三代PacBio RS II 和二代Illumina HiSeq 2000 测序技术对贝莱斯芽孢杆菌LG37进行全基因组测序,在此基础上利用NR、KEGG、eggNOG、GO和CARD数据库进行序列注释、分析,并通过本地Blast+对无机氮代谢相关基因进行挖掘。测序结果表明:1) 贝莱斯芽孢杆菌LG37的基因组为3 929 697 bp的环状染色体,GC含量为46.5%,包含3 854个蛋白质编码基因、86个tRNA 基因和27个rRNA基因。2) 共筛选出无机氮代谢相关候选基因94个,主要涉及编码感应蛋白、转录调控因子、转运蛋白、氧化还原酶和同化酶等,并对这些基因的GO功能进行了注释分析。综上,LG37 全基因组测序及无机氮代谢相关基因的分析为芽孢杆菌降低养殖水体中无机氮的研究提供了基因水平数据,为芽孢杆菌微生态制剂降低水体中无机氮的应用研究提供了理论依据。
    Abstract: It has been found that Bacillus velezensis can assimilate inorganic nitrogen efficiently. However, the underlying mechanism of inorganic nitrogen assimilation remains enigmatic. In order to elucidate the mechanism, we sequenced the complete genome of LG37 by PacBio RS II and Illumina HiSeq 2000, and then annotated and analyzed the sequence by the database of NR, KEGG, eggNOG, GO and CARD. Finally, we screened the genes related to inorganic nitrogen metabolism by local Blast+. The results show that: 1) The genome contained one circular chromosomal with a size of 3 929 697 bp and a GC-content of 46.5%. Gene prediction and annotation was performed to acquire a total of 3 854 protein-coding genes, 86 tRNA genes and 27 rRNA genes. 2) A total of 94 inorganic nitrogen metabolism candidate genes were screened by local Blast+. These genes were involved into coding sensing protein, transcriptional regulator, transporter, oxidoreductase and assimilator, etc.. In conclusion, the whole genome sequencing and data analysis of LG37 provide data at gene level and theoretical basis for functional study and application of Bacillus in reducing inorganic nitrogen in aquaculture water.
  • 中国对虾 (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  不同物质对中国对虾酚氧化酶活性的抑制作用
    Fig. 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  复合保鲜剂对中国对虾冷藏期间防黑效果评价
    Fig. 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值和菌落总数变化
    Fig. 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的样品,以此类推。
    Fig. 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   LG37 环状基因组图谱

    Figure  1.   Circular genome map of B. velezensis LG37 strain

    图  2   贝莱斯芽孢杆菌基于16S rDNA 序列构建的系统进化树

    Figure  2.   Neighbor-joining tree of B. velezensis LG37 based on 16S rDNA sequences

    图  3   COG功能分类

    Figure  3.   COG classification

    图  4   LG37基因功能注释KEGG代谢通路

    Figure  4.   Gene KEGG pathway classification map of LG37

    图  5   氮代谢通路 (ko00910)

    Figure  5.   Nitrogen metabolism (ko00910)

    表  1   贝莱斯芽孢杆菌LG37 基因组特性

    Table  1   Genome features of B. velezensis LG37 strain

    特性
    Feature
    数值
    Value
    基因组大小 Genome size/bp 3 929 697
    GC-含量 GC-content 46.5%
    质粒数量 Plasmid number 0
    总基因 Total genes 3 967
    蛋白编码基因 Protein-coding genes 3 854
    转运 RNA tRNA 86
    核糖体 RNA rRNA 27
    编码区域大小 Coding region size/bp 3 495 864
    编码区域GC-含量 GC-content of coding region 47.3%
    编码区域/全基因组 Coding region/Genome length 89.0%
    间隔区域大小 Intergenic region size/bp 433 833
    间隔区域占比 Ratio of intergenic region 11.0%
    下载: 导出CSV

    表  2   LG37 基因组氮代谢通路及其相关基因

    Table  2   Related genes of nitrogen metabolism pathways of LG37

    通路
    Pathway ID
    KEGG 描述
    KEGG description
    基因
    Gene
    ko00910 氮代谢 Nitrogen metabolism orf00490, orf00539, orf00544, orf00545, orf00546, orf00547, orf00817, orf01213, orf01254, orf01968, orf02226, orf02227, orf02368, orf03807, orf03808, orf03809, orf03810, orf03811
    M00531 同化硝酸盐还原 Assimilatory nitrate reduction orf00539, orf03807, orf03809
    M00530 异化硝酸盐还原 Dissimilatory nitrate reduction orf00544, orf00545, orf00547, orf03808, orf03809, orf03810, orf03811
    M00529 反硝化 Denitrification orf00544, orf00545, orf00547, orf03809
    M00804 完全硝化 Complete nitrification orf03809
    下载: 导出CSV

    表  3   无机氮代谢候选基因

    Table  3   Candidate genes of inorganic nitrogen metabolism

    基因
    Gene
    大小
    Size/bp
    蛋白
    Protein
    GO-分子功能
    GO-Molecular function
    orf00084 741 Type III pantothenate kinase YacB 泛酸激酶活性
    orf00106 246 Putative septation protein SpoVG 分子功能的负调控
    orf00148 885 Pyridoxal 5'-phosphate synthase subunit PdxS 谷氨酰胺水解活性
    orf00173 2 004 Nitrate reductase YyaE 硝酸还原酶活性
    orf00246 1 416 Arginine utilization regulatory protein RocR 转录因子结合
    orf00250 1 206 Ornithine aminotransferase RocD 鸟氨酸氧酸转氨酶活性
    orf00316 963 Iron(3+)-hydroxamate-binding protein YxeB 无机离子转运与代谢
    orf00379 1 017 Respiratory nitrate reductase NarI 硝酸还原酶活性
    orf00490 1 287 Glutamate dehydrogenase RocG 谷氨酸脱氢酶 (NAD+) 活性
    orf00535 1 377 Cytochrome cd1-nitrite reductase-like YwhL 亚硝酸盐还原酶活性
    orf00539 1 131 Nitrate transporter NarT 跨膜转运
    orf00544 3 687 Nitrate reductase alpha chain NarG 硝酸还原酶活性
    orf00545 1 464 Nitrate reductase beta chain NarH 硝酸还原酶活性
    orf00546 558 Nitrate reductase NarJ 未折叠蛋白结合
    orf00547 672 Nitrate reductase gamma chain NarI 硝酸还原酶活性
    orf00608 318 Urease subunit gamma 氨基酸转运与代谢
    orf00609 375 Urease subunit beta 氨基酸转运与代谢
    orf00610 1 710 Urease subunit alpha 氨基酸转运与代谢
    orf00623 351 Nitrogen regulatory protein P-II GlnB 酶调节活性
    orf00624 1 212 Ammonium transporter NrgA 铵跨膜转运蛋白活性
    orf00630 282 Stage III sporulation protein D SpoIIID DNA结合转录因子活性
    orf00783 726 Glucosamine-6-phosphate deaminase NagB 葡萄糖胺-6-磷酸脱氨酶活性
    orf00817 582 YvdA 无机离子转运与代谢
    orf00959 2 112 YvgW 阳离子转运ATP酶活性
    orf00966 1 716 Sulfite reductase [NADPH] CysI 亚硫酸盐还原酶 (NADPH) 活性
    orf00998 1 047 ABC transporter permease protein YvrB 转运体活性
    orf01042 357 Uncharacterized protein YusI 氧化还原酶活性
    orf01058 1 398 ABC transporter ATP-binding protein ATP结合
    orf01105 237 Nitrogen-fixing NifU domain-containing protein 铁硫簇结
    orf01175 2 403 Cation:proton antiporter 单价无机阳离子跨膜转运蛋白活性
    orf01213 1 032 Nitronate monooxygenase Ncd2 硝酸单加氧酶活性
    orf01254 549 YtiB 无机离子转运与代谢
    orf01262 813 Nitrate transport system permease protein YtlD 跨膜转运蛋白活性
    orf01263 783 Nitrate ABC transporter permease YtlC 跨膜转运
    orf01264 1 005 Nitrate ABC transporter periplasmic protein YtlA ATP酶活性
    orf01290 753 Quaternary-amine-transporting ATPase 传输ATP酶活性的季铵盐化合物
    orf01357 1 191 Nitric oxide dioxygenase 一氧化氮双加氧酶活性
    orf01435 1 962 Threonine--tRNA ligase 1 ThrS ATP结合
    orf01491 225 Spore germination protein GerE DNA结合
    orf01591 657 GlnP 氨基酸转运与代谢
    orf01592 651 GlnM 氨基酸转运与代谢
    orf01593 828 GlnH 氨基酸转运与代谢
    orf01617 801 Formate/nitrite transporter 跨膜转运蛋白活性
    orf01659 726 RNA polymerase sigma factor DNA结合转录因子活性
    orf01864 855 Nitrogen assimilation regulatory protein nac DNA结合转录因子活性
    orf01911 450 Ferric uptake regulation protein DNA结合转录因子活性
    orf01968 1 275 Glutamate dehydrogenase RocG 谷氨酸脱氢酶硝酸还原酶活性
    orf02215 1 983 Nitrate reductase 硝酸还原酶活性
    orf02225 903 HTH-type transcriptional regulator GltC DNA结合转录因子活性
    orf02226 903 Glutamate synthase (NADPH/NADH) GltB 谷氨酸合酶 (NADPH) 活性
    orf02227 4 560 Glutamate synthase (NADPH/ NADH) GltD 谷氨酸合酶 (NADPH) 活性
    orf02368 1 335 Glutamine synthetase GlnA 谷氨酸氨连接酶活性
    orf02369 405 HTH-type transcriptional regulator GlnR DNA结合
    orf02384 315 Ammonium compound efflux SMR transporter 膜的组成部分
    orf02385 354 Ammonium compound efflux SMR transporter 膜的组成部分
    orf02417 261 Stage V sporulation protein S 核酸结合
    orf02418 795 2',3'-cyclic-nucleotide 2'-phosphodiesterase YmdB 2',3'-环核苷酸2'-磷酸二酯酶活性
    orf02465 723 Uridylate kinase PyrH ATP结合
    orf02483 363 Chemotaxis protein CheY 磷脂酶信号转导系统
    orf02588 783 RNA polymerase sigma factor DNA结合转录因子活性
    orf02638 930 Glutaminase 氨基酸转运与代谢
    orf02680 666 Potassium uptake protein KtrA 阳离子跨膜转运蛋白活性
    orf02747 1 914 YkvW 阳离子转运ATP酶活性
    orf02756 342 Putative transcriptional regulator 转录调控,DNA模板
    orf02787 1 353 YkrM 阳离子跨膜转运蛋白活性
    orf02794 783 Uncharacterized membrane protein YkoY 膜的组成部分
    orf02804 742 HTH-type transcriptional regulator TnrA 核心启动子结合
    orf02828 315 Ammonium compound efflux SMR transporter 膜的组成部分
    orf02830 339 Ammonium compound efflux SMR transporter 膜的组成部分
    orf02854 999 Anion permease 无机磷酸盐跨膜转运蛋白活性
    orf03008 396 ArsC family transcriptional regulator SpxA 电子转移活性
    orf03161 1 215 Cation/H (+) antiporter YhaU 溶质:质子逆向转运活性
    orf03186 396 Putative fluoride ion transporter CrcB 无机阴离子跨膜转运蛋白活性
    orf03187 360 Putative fluoride ion transporter CrcB 无机阴离子跨膜转运蛋白活性
    orf03209 435 HTH-type transcriptional regulator NsrR DNA结合
    orf03262 834 ABC-type nitrate transport system 离子跨膜转运
    orf03263 990 Putative binding protein SsuA ATP酶活性
    orf03264 768 Aliphatic sulfonates import protein SsuB 阴离子跨膜转运蛋白活性
    orf03448 720 Probable transcriptional regulatory protein 转录调控,DNA模板
    orf03546 330 Ammonium compound efflux SMR transporter 膜的组成部分
    orf03547 315 Ammonium compound efflux SMR transporter 膜的组成部分
    orf03597 963 Arsenic resistance protein 无机阴离子跨膜转运蛋白活性
    orf03677 822 Probable manganese catalase YdbD 无机离子转运与代谢
    orf03682 366 Ammonium compound efflux SMR transporter 膜的组成部分
    orf03807 1 206 Assimilatory nitrate reductase NasA 硝酸铁氧还蛋白还原酶活性
    orf03808 2 328 Nitrite reductase large subunit NasB 亚硝酸盐还原酶 [NAD(P)H] 活性
    orf03809 2 133 Assimilatory nitrate reductase NasC 硝酸还原酶活性
    orf03810 2 418 Nitrite reductase [NAD(P)H] NasD 亚硝酸盐还原酶 [NAD(P)H] 活性
    orf03811 321 Assimilatory nitrite reductase [NAD(P)H] NasE 亚硝酸盐还原酶 [NAD(P)H] 活性
    orf03812 1 440 NasF 辅酶代谢
    orf03840 1 257 Transport system atp-binding protein opuaa 传输ATP酶活性的季铵盐化合物
    orf03889 984 Glutaminase 1 GlsA1 谷氨酰胺酶活性
    orf03890 1 437 GlnT 假定的钠/谷氨酰胺转运体
    orf03963 1 434 Sodium-independent anion transporter 次级活性硫酸盐跨膜转运蛋白活性
    下载: 导出CSV
  • [1] 胡晓娟, 文国樑, 田雅洁, 等. 不同培养条件下菌株NB5 对氨氮的去除效果研究[J]. 南方水产科学, 2020, 16(6): 89-96.
    [2] 方金龙, 王元, 房文红, 等. 氨氮胁迫下白斑综合征病毒对凡纳滨对虾的致病性[J]. 南方水产科学, 2017, 13(4): 52-58. doi: 10.3969/j.issn.2095-0780.2017.04.007
    [3] 肖炜, 李大宇, 徐杨, 等. 慢性氨氮胁迫对吉富罗非鱼幼鱼生长、免疫及代谢的影响[J]. 南方水产科学, 2015, 11(4): 81-87. doi: 10.3969/j.issn.2095-0780.2015.04.012
    [4] 郭国强, 孙红飞, 张永耀. 分子氨对草鱼鱼种红细胞渗透脆性的影响[J]. 水产科学, 2010, 29(8): 489-491. doi: 10.3969/j.issn.1003-1111.2010.08.011
    [5] 孙侦龙, 朱永祥, 刘大勇, 等. 非离子氨氮和亚硝酸盐氮对暗纹东方鲀稚鱼的急性毒性试验[J]. 水产科学, 2015, 34(3): 135-139.
    [6]

    KUEBUTORNYE F K A, ABARIKE E D, LU Y. A review on the application of Bacillus as probiotics in aquaculture[J]. Fish Shellfish Immunol, 2019, 87: 820-828. doi: 10.1016/j.fsi.2019.02.010

    [7]

    LIU G X, VIJAYARAMAN S B, DONG Y J, et al. Bacillus velezensis LG37: transcriptome profiling and functional verification of GlnK and MnrA in ammonia assimilation[J]. BMC Genom, 2020, 21(1): 215. doi: 10.1186/s12864-020-6621-1

    [8]

    GOODWIN S, MCPHERSON J D, MCCOMBIE W R. Coming of age: ten years of next-generation sequencing technologies[J]. Nat Rev Genet, 2016, 17(6): 333-51.

    [9]

    GALPERIN M Y, KRISTENSEN D M, MAKAROVA K S, et al. Microbial genome analysis: the COG approach[J]. Brief Bioinform, 2019, 20(4): 1063-1070. doi: 10.1093/bib/bbx117

    [10]

    CASTAÑEDA C D, GAMBLE J N, WAMSLEY K G S, et al. In ovo administration of Bacillus subtilis serotypes effect hatchability, 21-day performance, and intestinal microflora[J]. Poult Sci, 2021, 100(6): 101125. doi: 10.1016/j.psj.2021.101125

    [11]

    TAKAMI H, TOYODA A, UCHIYAMA I, et al. Complete genome sequence and expression profile of the commercial lytic enzyme producer Lysobacter enzymogenes M497-1[J]. DNA Res, 2017, 24(2): 169-177.

    [12]

    MOLTON J S, LEE I R, BERTRAND D, et al. Stool metagenome analysis of patients with Klebsiella pneumoniae liver abscess and their domestic partners[J]. Int J Infect Dis, 2021, 107: 1-4. doi: 10.1016/j.ijid.2021.04.012

    [13] 朱玉玲, 彭晶, 唐诗哲, 等. 宏基因组技术在极端环境酯酶挖掘中的应用[J]. 生命科学研究, 2021, 25(2): 169-175.
    [14] 魏大鹏, 单洪伟, 马甡, 等. 混料设计优化复合菌剂比例的研究[J]. 南方水产科学, 2014, 10(1): 86-91. doi: 10.3969/j.issn.2095-0780.2014.01.013
    [15]

    TAMURA K, STECHER G, PETERSON D, et al. MEGA6: molecular evolutionary genetics analysis version 6.0[J]. Mol Biol Evol, 2013, 30(12): 2725-2729. doi: 10.1093/molbev/mst197

    [16]

    NAKANO K, SHIROMA A, SHIMOJI M, et al. Advantages of genome sequencing by long-read sequencer using SMRT technology in medical area[J]. Hum Cell, 2017, 30(3): 149-161. doi: 10.1007/s13577-017-0168-8

    [17]

    GAN Y Q, ZHANG T, GAN Y Q, et al. Complete genome sequences of two Enterococcus faecium strains and comparative genomic analysis[J]. Exp Ther Med, 2020, 19(3): 2019-2028.

    [18]

    BENNETT J P Jr, KEENEY P M, BROHAWN D G. RNA sequencing reveals small and variable contributions of infectious agents to transcriptomes of postmortem nervous tissues from amyotrophic lateral sclerosis, Alzheimer's disease and Parkinson's disease subjects, and increased expression of genes from disease-activated microglia[J]. Front Neurosci, 2019, 13: 235. doi: 10.3389/fnins.2019.00235

    [19]

    EMSLEY C, KING S, NYULASI I, et al. A GLIMmer of insight into lung transplant nutrition: enhanced detection of malnutrition in lung transplant patients using the GLIM criteria[J]. Clin Nutr, 2021, 40(5): 2521-2526. doi: 10.1016/j.clnu.2021.02.047

    [20]

    CHAN P P, LOWE T M. tRNAscan-SE: searching for tRNA genes in genomic sequences[J]. Methods Mol Biol, 2019, 1962: 1-14.

    [21]

    LAGESEN K, HALLIN P, RØDLAND E A, et al. RNAmmer: consistent and rapid annotation of ribosomal RNA genes[J]. Nucleic Acids Res, 2007, 35(9): 3100-8. doi: 10.1093/nar/gkm160

    [22]

    CAMACHO C, COULOURIS G, AVAGYAN V, et al. BLAST+: architecture and applications[J]. BMC Bioinformatics, 2009, 10: 421. doi: 10.1186/1471-2105-10-421

    [23]

    NAYAK V S. Nonelectrolytic production of caustic soda and hydrochloric acid from sodium chloride[J]. Ind Eng Chem Res, 1996, 35(10): 3808-3811. doi: 10.1021/ie960043h

    [24] 张静, 高婷婷, 李勇, 等. 蛋白营养对工业养殖大菱鲆 (Scophthatmus maximus L.)幼鱼生长、氨氮排泄及肌肉氨基酸的效应[J]. 渔业科学进展, 2016, 37(6): 34-41. doi: 10.11758/yykxjz.20151124001
    [25]

    EBELING J M, TIMMONS M B, BISOGNI J J. Engineering analysis of the stoichiometry of photoautotrophic, autotrophic, and heterotrophic removal of ammonia-nitrogen in aquaculture systems[J]. Aquaculture, 2006, 257(1/2/3/4): 346-358.

    [26] 雷阳, 张倩, 陈钰, 等. 对虾养殖高效降解氨氮微生态菌的筛选与鉴定[J]. 福建农业科技, 2019(10): 16-20.
    [27]

    ZHU X H, ZHANG S M, ZHOU L Y, et al. Probiotic potential of Bacillus velezensis: antimicrobial activity against non-O1 Vibrio cholerae and immune enhancement effects on Macrobrachium nipponense[J]. Aquaculture, 2021, 541: 736817. doi: 10.1016/j.aquaculture.2021.736817

    [28] 张德锋, 高艳侠, 可小丽, 等. 贝莱斯芽孢杆菌LF01基因组序列分析及其代谢产物的生防作用[J]. 水产学报, 2022, 46(2): 196-206.
    [29]

    FAN B, BLOM J, KLENK H P, et al. Bacillus amyloliquefaciens, Bacillus velezensis, and Bacillus siamensis form an "operational group B. amyloliquefaciens" within the B. subtilis species complex[J]. Front Microbiol, 2017, 8: 22.

    [30]

    MURUKE M S H, CAMP H J O D, SEMESI A K, et al. The level of enzymes involved in the allantoin metabolism of Bacillus fastidiosus grown under different conditions[J]. Curr Microbiol, 1995, 30(1): 45-7. doi: 10.1007/BF00294523

    [31]

    NIU T, LV X, LIU Z, et al. Synergetic engineering of central carbon and nitrogen metabolism for the production of N-acetylglucosamine in Bacillus subtilis[J]. Biotechnol Appl Biochem, 2020, 67(1): 123-132. doi: 10.1002/bab.1845

    [32]

    WANG B, ZHANG D, CHU S, et al. Genomic analysis of Bacillus megaterium NCT-2 reveals its genetic basis for the bioremediation of secondary salinization soil[J]. Int J Genomics, 2020, 2020: 4109186.

    [33]

    MASSAD R S, LOUBET B, TUZET A, et al. Relationship between ammonia stomatal compensation point and nitrogen metabolism in arable crops: current status of knowledge and potential modelling approaches[J]. Environ Pollut, 2008, 154(3): 390-403.

    [34]

    TANG Y, YE Z, WEI Y, et al. Vertebrate paralogous CRMPs in nervous system: evolutionary, structural, and functional interplay[J]. J Mol Neurosci, 2015, 55(2): 324-34. doi: 10.1007/s12031-014-0327-2

    [35] 顾志良, 耿拓宇. 鸡重要性状主效基因和QTL 的研究进展[J]. 中国家禽, 2003(S1): 130-134.
    [36]

    AWADA H, THAPA B, VISCONTE V. The genomics of myelodysplastic syndromes: origins of disease evolution, biological pathways, and prognostic implications[J]. Cells, 2020, 9(11): 2512. doi: 10.3390/cells9112512

    [37] 丁忠涛, 张锐, 郭三堆. 棉花抗逆相关基因GhDr1的克隆及生物信息学分析[J]. 生物技术通报, 2011(1): 99-106.
  • 期刊类型引用(1)

    1. 程思忠,谢岩黎,孙淑敏,马卫宾,李倩,杨玉辉. 黄曲霉毒素B_1降解菌的筛选鉴定及降解酶挖掘. 食品安全质量检测学报. 2023(04): 1-7 . 百度学术

    其他类型引用(3)

推荐阅读
Lh原油和0# 柴油乳化液对凡纳滨对虾肝胰腺抗氧化酶活性及相关功能基因表达的影响
沈楚焰 et al., 南方水产科学, 2025
尖翅燕鱼染色体水平基因组与特征分析
欧阳焱 et al., 南方水产科学, 2024
基于cbam-unet的大菱鲆鱼卵识别计数方法
钱程 et al., 南方水产科学, 2024
草鱼subfatin分子鉴定及表达特性分析
杨博雅 et al., 南方水产科学, 2024
构树bzip基因家族鉴定及硒代谢相关基因筛选
BU Xian-chen et al., JOURNAL OF SOUTHERN AGRICULTURE, 2025
一株多重耐药准肺炎克雷伯菌的基因组和耐药基因分析
崔晶花 et al., 遵义医科大学学报, 2024
Pheiges: all-cell-free phage synthesis and selection from engineered genomes
Levrier, Antoine et al., NATURE COMMUNICATIONS, 2024
Genomic and transcriptomic analysis of checkpoint blockade response in advanced non-small cell lung cancer
Ravi, Arvind et al., NATURE GENETICS, 2023
Genom7gnet: an efficient n7-methylguanosine site prediction approach based on a nucleotide language model
IEEE-ACM TRANSACTIONS ON COMPUTATIONAL BIOLOGY AND BIOINFORMATICS, 2024
Identification of transcription factors associated with the disease-free survival of triple-negative breast cancer through weighted gene co-expression network analysis
CYTOJOURNAL, 2024
Powered by
图(5)  /  表(3)
计量
  • 文章访问数:  916
  • HTML全文浏览量:  201
  • PDF下载量:  74
  • 被引次数: 4
出版历程
  • 收稿日期:  2021-05-19
  • 修回日期:  2021-07-06
  • 录用日期:  2021-09-28
  • 网络出版日期:  2021-10-18
  • 刊出日期:  2022-06-04

目录

/

返回文章
返回