牡蛎相关微病毒科噬菌体基因组的鉴定和进化分析

刘敏, 刘畅, 刘广锋, 朱鹏, 谢科明, 姜敬哲

刘敏, 刘畅, 刘广锋, 朱鹏, 谢科明, 姜敬哲. 牡蛎相关微病毒科噬菌体基因组的鉴定和进化分析[J]. 南方水产科学, 2023, 19(5): 134-142. DOI: 10.12131/20230038
引用本文: 刘敏, 刘畅, 刘广锋, 朱鹏, 谢科明, 姜敬哲. 牡蛎相关微病毒科噬菌体基因组的鉴定和进化分析[J]. 南方水产科学, 2023, 19(5): 134-142. DOI: 10.12131/20230038
LIU Min, LIU Chang, LIU Guangfeng, ZHU Peng, XIE Keming, JIANG Jingzhe. Identification and evolutionary analysis of genome of oyster-associated Microviridae[J]. South China Fisheries Science, 2023, 19(5): 134-142. DOI: 10.12131/20230038
Citation: LIU Min, LIU Chang, LIU Guangfeng, ZHU Peng, XIE Keming, JIANG Jingzhe. Identification and evolutionary analysis of genome of oyster-associated Microviridae[J]. South China Fisheries Science, 2023, 19(5): 134-142. DOI: 10.12131/20230038

牡蛎相关微病毒科噬菌体基因组的鉴定和进化分析

基金项目: 广东省重点领域研发计划项目 (2022B1111030001);国家自然科学基金面上项目 (31972847);中国水产科学研究院中央级公益性科研院所基本科研业务费专项资金资助 (2023TD44);中国水产科学研究院南海水产研究所中央级公益性科研院所基本科研业务费专项资金资助 (2021SD05, 2022XK02);广东省省级现代农业产业技术体系 (2019KJ141)
详细信息
    作者简介:

    刘 敏 (1999—),女,硕士研究生,研究方向为渔业生物病害防治。E-mail: 1099535669@qq.com

    通讯作者:

    姜敬哲  (1980—),男,研究员,博士,研究方向为渔业生物病害防治。E-mail: jingzhejiang@hotmail.com

  • 中图分类号: S 917.1

Identification and evolutionary analysis of genome of oyster-associated Microviridae

  • 摘要:

    具有滤食习性的牡蛎富集了水体中包括病毒在内的大量病原体,是一个极具价值的病毒库。为了对牡蛎相关的病毒库进行深入研究,在前期对华南沿海多地采集的香港牡蛎 (Crassostrea hongkongensis) 进行病毒组测序,并对测序数据进行质控、组装及物种注释后,挑选其中5条被鉴定为微病毒科 (Microviridae) 的基因组序列进行多维度分析,包括宿主预测、开放阅读框预测、主要衣壳蛋白系统发育与三维结构预测、主要衣壳蛋白与外部支架蛋白的进化关联以及病毒丰度分析等。结果显示,5个病毒的宿主均为埃希氏菌属 (Escherichia);其中病毒基因组序列HSd1-5344568聚类在Bullavirinae分支中,说明其为该亚科成员;其余4条未聚类到任何已知亚科中,应属于一个单独的未分类亚科;微病毒主要衣壳蛋白和外部支架蛋白的进化树之间的联系表明两个蛋白的进化规律不同。

    Abstract:

    Oysters with filter feeding habits enrich a large number of pathogens including viruses in water, like a very valuable virus bank. To conduct an in-depth research on oyster-related viruses, we sequeced the Hong Kong oyster (Crassostrea hongkongensis) collected along the South China coast. After the quality control, assembly and taxonomy annotation of the sequencing data, we selected five genomic sequences which were identified as Microviridae for multi-dimensional analyses, such as host prediction, open reading frame and gene function prediction, phylogeny and three-dimensional structure prediction of major capsid proteins, evolutionary association between major capsid proteins and external scaffold proteins, as well as virus abundance analysis. The results show that the hosts of the five viruses were all Escherichia; one of the virus genome sequences was clustered in Bullavirinae branch, which indicates that it is a member of the subfamily; the other four genome sequences were not clustered into any known subfamilies so they should belong to a single unclassified subfamily; the relationship between the evolutionary tree of the main capsid proteins and external scaffold proteins indicates that the evolution rules of the two proteins were different.

  • 氨氮是水产养殖中蛋白质的分解产物,主要由剩余食物和鱼粪分解产生[1],其主要在水中以氨 (NH3) 和离子态铵 (NH4 +) 的形式存在。NH3是一种脂溶性氨水,容易穿透细胞膜而对鱼类造成毒害[2]。大多数鱼类对环境中的氨氮浓度很敏感,是导致鱼病发生的主要环境因子之一。氨氮胁迫会损伤鱼类的鳃组织,渗透进入血液淋巴组织,机体抗氧化酶系统受损,鳃、消化器官等组织结构发生病变从而影响机体呼吸及消化功能,甚至引起发病死亡[3-4]。因此,研究氨氮胁迫对水产动物生理指标及组织结构的影响,可为胁迫因子对机体的影响监测提供基础指标并明确鱼类对氨氮的耐受性,对促进水产养殖的可持续发展具有重要意义。目前,氨氮对水生生物影响的研究已有很多,如黑鲷 (Acanthopagrus schlegeli)[5]、牙鲆 (Paralichthys olivaceus)[6]、黄颡鱼 (Pelteobagrus fulvidraco)[7]、克氏原螯虾 (Procambarus clarkia)[8]、奥尼罗非鱼 (Oreochromis niloticus×O. areus)[9]、方斑东风螺 (Babylonia areolata)[10]等。黄厚见等[11]发现,随着氨氮浓度的升高,梭鱼 (Liza haematocheila) 的胃蛋白酶活力显著下降。不同氨氮浓度胁迫下,福瑞鲤 (Cyprinus carpio)[12] 与刺参 (Stichopus japonicus) [13]的Na+/K+-ATP酶与呼吸代谢酶活力随氨氮浓度的升高而降低。研究表明高于0.60 mg·L−1的氨氮会使吉富罗非鱼幼鱼的组织结构受损,显著影响罗非鱼的消化酶水平,从而使其生长率明显低于对照组[14]。团头鲂 (Megalobrama amblycephala) 幼鱼在氨氮胁迫6 h后,鳃丝血管扩张,上皮组织增生,肝细胞肿胀,胁迫48 h,鳃小片卷曲,上皮细胞部分脱落;肝细胞部分溶解、血窦扩张,形成点状病灶,组织结构受到严重破坏[15]。不同的水生生物对氨氮的耐受能力不同,不同规格的同一品种鱼对氨氮的耐受能力也存在差异[16]

    翘嘴鳜 (Siniperca chuatsi) 为典型肉食性鱼类,喜食活饵。目前,翘嘴鳜的养殖由于过度追求数量,养殖密度高,养殖过程中投喂大量的饵料鱼,饵料残渣和鱼的代谢产物在水中快速积累易导致水体中的氨氮浓度快速升高[17-18]。生产实践表明翘嘴鳜对氨氮非常敏感,氨氮胁迫下其生长特性、消化功能均可能发生变化,进而影响了养殖效益。目前国内外针对翘嘴鳜幼鱼的研究主要集中在生长选育、食性改良等方面,而关于氨氮胁迫的研究却很少。本实验以翘嘴鳜幼鱼鳃、消化组织为研究对象,分析半致死浓度(LC50)氨氮胁迫下,不同胁迫时间,翘嘴鳜幼鱼鳃丝呼吸代谢酶 [ 乳酸脱氢酶 (LDH)、己糖激酶(HK)]与Na+/K+-ATP酶以及胃囊、肠道消化酶[ 淀粉酶 (AMS)、胃蛋白酶、脂肪酶 (LPS)]活力的变化,探究急性氨氮胁迫对翘嘴鳜幼鱼鳃与消化道生理功能的影响机制,进而实时监测养殖水体的氨氮浓度变化,预防病害,以期为翘嘴鳜的生态健康养殖提供依据。

    本研究所用翘嘴鳜幼鱼均选自南京农业大学无锡渔业学院南泉养殖基地,暂养于室内控温循环水养殖系统 (400 L·桶−1)。选取体表无损伤、规格整齐的翘嘴鳜幼鱼300尾。实验开始前,翘嘴鳜幼鱼于室内控温循环水养殖系统暂养14 d。暂养期间保持自然光照周期,暂养水体实验条件为溶解氧 (DO)质量浓度≥6 mg·L−1、氨氮质量浓度≤0.05 mg·L−1、水温 (22.0±0.5) ℃。暂养期间饵料鱼每2 d投放1次,数量为翘嘴鳜的2倍。

    室内循环水系统保持水温 (22.1±0.3) ℃,暂养14 d后选取规格整齐的翘嘴鳜幼鱼120尾[ 体质量 (13.55±0.69) g,体长 (9.26±0.35) cm]进行LC50测试。设定0、30、60和90 mg·L−1氯化铵 (NH4Cl) (试剂号75-09-2,国药集团化学试剂有限公司) 4个浓度梯度。每组浓度设置3个平行,每个平行放置10尾鱼。在96 h观测期中每隔8 h使用W-1型多参数水质分析仪 (杭州陆恒生物科技有限公司) 检测各水体氨氮浓度,将其调整为实验设计的氨氮浓度,每24 h换水1/4。观测翘嘴鳜幼鱼的活性变化,计数并捞出死亡个体 (鱼体侧翻且鳃盖停止扇动视为死亡)。

    暂养结束后,依据氨氮胁迫96 h LC50,设置对照组 (0 mg·L−1) 与实验组 (48.65 mg·L−1),实验组与对照组设置3个重复,选取体质量相近、活性良好的翘嘴鳜幼鱼180尾 [体质量 (13.31±0.49) g、体长 (9.73±0.46) cm] 随机平均放置于6个养殖桶中 (400 L),每桶30尾。对照组为完全曝气的自来水 (氨氮实测值为0.05 mg·L−1)。实验期间水温 (22.0±0.5) ℃, pH维持在7.70±0.10,DO质量浓度维持在 (6.13±0.12) mg·L−1, 氨氮胁迫时长为96 h,实验期间保持水体静止,停止进食,持续充氧,每隔8 h使用W-1型多参数水质分析仪 (杭州陆恒生物科技有限公司) 检测各水体氨氮浓度,将其调整为实验设计的氨氮浓度,每24 h换水1/4。

    在氨氮胁迫第0、第6、第12、第24、第48和第96小时,每桶分别随机选择3尾翘嘴鳜并用50 mg·L−1 MS-222轻度麻醉,采集翘嘴鳜幼鱼活体的第二鳃弓处鳃丝、胃囊、肠道,−80 ℃冻存用于后续酶活力测定。

    Na+/K+-ATP酶、LDH、HK、胃蛋白酶、LPS、AMS检测试剂盒购自南京建成生物研究所。取0.1 g样品加入9倍体积的生理盐水研磨后3 000 r·min−1 4 ℃离心10 min,取上清液−80 ℃保存备用。各种酶活力按试剂盒中的说明书进行操作,均采用比色法进行测定。

    采用定磷法测定。取同批10%组织匀浆加0.68%的生理盐水稀释至1%,用考马斯亮兰法测定组织中的蛋白含量。以每小时每毫克组织蛋白中ATP酶分解ATP产生1 μmol无机磷的量为1个ATP酶活力单位。

    每毫升上清37 ℃与基质作用15 min,在反应体系中产生1 μmol丙酮酸为1单位。通过测定样品孔与标准孔的吸光值来计算酶活力。

    在37 ℃、pH 7.6的条件下,每克组织蛋白在本反应体系中每分钟生成1 mmol·L−1的NADP定义为1个酶活力单位。在340 nm波长处,测定吸光度的增加值来计算酶活力。

    每毫克组织蛋白37 ℃下, 每分钟分解蛋白生成1 μg氨基酸相当于1个酶活力单位,通过测定660 nm处的吸光值进行计算。

    在37 ℃条件下,每升上清与底物反应1 min后,每消耗1 μmol底物为一个酶活力单位。通过测定420 nm处吸光度值的差值来计算酶活力。

    组织中每毫克蛋白在37 ℃与底物作用30 min,水解10 mg淀粉定义为1个淀粉酶活力单位,根据与碘生成的蓝色复合物来计算AMS活力。

    翘嘴鳜幼鱼96 h LC50分析采用直线内插法。实验结果分析使用SPSS 20.0 软件。酶活力检测结果以3个平行组数据 “平均值±标准差 ($\overline X \pm {\rm{SD}} $)”表示,实验相同采样时间实验组与对照组的比较采用独立样本t检验进行统计分析。对相同处理组不同时间点的数据进行单因素方差分析 (One-Way ANOVA),差异显著时,使用Duncan's检验法进行差异性分析 (P<0.05)。

    4种NH4Cl浓度处理96 h后,翘嘴鳜幼鱼呈现出不同的应激行为甚至死亡。起初躁动不安、四处游动,随着时间延长,反应更加剧烈。48 h后游动迟缓,尾鳍呈不同程度向内侧弓起,90 mg·L−1实验组鱼尾鳍与身体几乎成直角。在翘嘴鳜幼鱼氨氮胁迫96 h后,0、30、60、90 mg·L−1氨氮浓度组分别死亡0、8、18、29尾,因此计算得出96 h氨氮LC50为48.65 mg·L−1

    在96 h氨氮胁迫中各组翘嘴鳜幼鱼鳃Na+/K+-ATP酶活力随时间呈降低-升高-降低的变化趋势 (图1-a)。氨氮胁迫至第6小时,鳃Na+/K+-ATP酶活力显著降低 (P<0.05),随即呈升高趋势;至第48小时已显著高于对照组,达到最高值 (对照组的1.84倍);胁迫至第96小时,Na+/K+-ATP酶活力降低,但仍显著高于对照组 (P<0.05)。

    图  1  急性氨氮胁迫对翘嘴鳜幼鱼鳃呼吸代谢酶活力的影响
    不同小写字母表示同一处理不同时间点之间差异显著 (P<0.05);*. 实验组和对照组在胁迫后同一时间点差异显著 (P<0.05);实验组氨氮质量浓度为48.65 mg·L−1,对照组氨氮质量浓度为0 mg·L−1
    Figure  1.  Effect of acute ammonia stress on respiratory metabolism activity in gill of S. chuatsi
    Different lowercase superscripts indicate significant difference (P<0.05) in the same treatment at different time; *. Significant difference between ammonia stress group and the control group at the same time (P<0.05); concentrations of ammonia in the stress group and the control group are 48.65 mg·L−1 and 0 mg·L−1, respectively.

    急性氨氮胁迫96 h,LDH活力呈升高的变化趋势 (图1-b)。至第12小时,鳃LDH活力开始呈现升高的趋势;至第96小时显著高于对照组,达到最高值 (对照组的2.13倍)。各组翘嘴鳜幼鱼鳃HK活力呈逐步升高的变化趋势 (图1-c)。至第6小时,鳃HK活力显著升高 (P<0.05);至第48小时显著高于对照组,达到最高值 (对照组的1.70倍)。

    在96 h急性氨氮胁迫中,胃囊AMS活力随时间呈先降低后升高的变化趋势 (图2-a)。对照组胃囊AMS活力与各实验组之间差异显著 (P<0.05)。氨氮胁迫至第6小时,胃囊AMS活力显著降低 (P<0.05);胁迫至第12小时,胃囊AMS活力达到最低值 (对照组的0.37倍),随即呈现升高的趋势;至第96小时已显著高于对照组,达到最高值 (对照组的1.23倍)。实验组胃囊胃蛋白酶活力随时间呈升高-降低的变化趋势 (图2-b)。受氨氮胁迫的影响,胁迫至第6小时,胃囊胃蛋白酶活力显著降低 (P<0.05),随即呈升高的趋势;胁迫至第48小时,胃囊胃蛋白酶活力达到最高值 (对照组的1.44倍);胁迫至第96小时,实验组胃囊胃蛋白酶活力降低,仍显著高于对照组 (P<0.05)。实验组胃囊LPS随时间呈升高-降低的趋势 (图2-c)。氨氮胁迫至第12小时,胃囊LPS活力显著升高 (P<0.05);第48小时胃囊LPS活力达到最高值 (对照组的1.71倍);胁迫至第96小时,实验组胃囊LPS活力降低,仍显著高于对照组 (P<0.05)。

    图  2  急性氨氮胁迫对翘嘴鳜幼鱼胃、肠道消化酶活力的影响
    Figure  2.  Effect of acute ammonia stress on digestive enzyme activity in stomach and intestine of S. chuatsi

    在96 h急性氨氮胁迫中,各组翘嘴鳜幼鱼肠道AMS活力随时间呈先降低后升高的变化趋势 (图2-d)。受氨氮胁迫的影响,对照组肠道AMS活力与各实验组差异显著 (P<0.05)。氨氮胁迫6 h,AMS活力显著降低 (P<0.05);第24小时肠道AMS活力达到最低值 (对照组的0.49倍),随即呈升高趋势;至第48小时已显著高于对照组,达到最高值 (对照组的1.12倍)。各组翘嘴鳜幼鱼肠道胃蛋白酶活力随时间呈先降低后升高的变化趋势 (图2-e)。氨氮胁迫至第6小时,肠道胃蛋白酶活力降低,但与对照组无显著差异 (P<0.05);胁迫至第12小时,实验组肠道胃蛋白酶活力显著降低,达到最低值 (对照组的0.89倍),随即呈升高趋势;至第96小时已显著高于对照组,达到最高值 (对照组的1.15倍)。实验组肠道LPS随时间呈逐渐升高趋势 (图2-f)。胁迫12 h肠道LPS活力显著升高 (P<0.05);至第96小时肠道LPS活力达到最高值 (对照组的2.99倍)。

    Na+/K+-ATP酶对环境渗透压的变化十分敏感,在硬骨鱼类中具有调节渗透压功能的重要作用[19],当水体环境中的高浓度氨氮通过鱼类鳃部进入机体时,鱼类也能够通过广泛分布于鳃组织中的Na+/K+-ATP酶与Na+/K+ (NH4 +) 载体以跨膜运输的方式将氨排出体外[20-21]。攀鲈 (Anabas testudineus) 在氨氮胁迫下其鳃部Na+/K+-ATP酶活力在氨氮胁迫1和6 d后均显著升高,且与胁迫时间呈正相关[22],这表明水体环境中氨氮浓度的变化能够诱导或抑制Na+/K+-ATP酶活力。本实验发现在氨氮胁迫下翘嘴鳜幼鱼鳃部Na+/K+-ATP酶活力持续升高,说明随着胁迫时间的延长,机体诱导提高了Na+/K+-ATP酶活力,通过泌氯细胞膜两侧的转运蛋白与Na+/K+ (NH4 +) 载体将氨排出体外。胁迫至第96小时,Na+/K+-ATP酶仍显著高于对照组,但已呈下降趋势,这可能是由于高浓度的氨氮影响了Na+/K+-ATP酶的蛋白结构,使酶活力降低。鳃部组织细胞功能虽然受到氨氮毒性的影响,但其酶活力仍显著高于对照组,表明鳃主动渗透调节能力仍能通过排出氨氮维持鱼体渗透压平衡。

    研究发现随着氨氮胁迫浓度的升高,硬骨鱼的耗氧率与排氨率均受到抑制[23]。HK与LDH是糖酵解过程中的重要酶类,HK能够将葡萄糖磷酸化以进行下一步的代谢反应,是糖无氧酵解反应的限速酶。LDH能够在无氧条件下催化丙酮酸生成乳酸,是生物无氧代谢的标志酶[24-25]。本研究发现翘嘴鳜幼鱼鳃组织中HK与LDH活力随氨氮胁迫时间延长而升高,并一直维持在较高水平,这说明氨氮胁迫下机体有氧呼吸过程可能受到抑制,并通过提高无氧代谢的方式提供机体所需能量。呼吸代谢模式从有氧呼吸向无氧呼吸转变的趋势,在低氧胁迫研究中则更为明显,如低氧胁迫下花鲈(Lateolabrax japanicus)幼鱼肝脏中糖原含量显著降低,LDH活力显著升高[26]。于晓[13]也发现氨氮胁迫下刺参组织中HK活力在短时间内迅速升高。因此在环境中氨氮水平升高时翘嘴鳜幼鱼呼吸代谢功能受氨氮毒性的影响,从而通过增加无氧代谢水平以维持机体能量供应,并缓解部分供氧压力,有利于机体提高对氨氮环境的适应性。

    在氨氮胁迫下,随着外源性氨氮不断通过鳃、皮肤、肠道等途径进入鱼体,导致鱼类体内氨氮浓度升高,进而引发氨中毒,鱼类会产生代谢紊乱、生长迟缓、组织损伤病变等一系列毒性效应[27]。高氨氮环境下,鱼类为了维持体内氨氮积累-代谢的平衡,不仅能够活化抗氧化酶、细胞修复因子、免疫因子以修复组织氧化损伤,通过排泄系统、呼吸系统将过多的氨氮排出体外,还可以通过降低活动频率、消化代谢水平以减少体内氨的生成,从而降低氨氮毒性水平[28]。本实验中,氨氮胁迫至第6小时,翘嘴鳜幼鱼胃囊与肠道AMS活力、胃囊胃蛋白酶活力显著下降 (P<0.05)。这说明氨氮胁迫初期,鱼类一方面降解、排出外源性氨氮,同时降低机体消化代谢水平,进一步降低体内氨氮水平,以规避毒性效应的不利影响。在三疣梭子蟹 (Portunus trituberculatus) 的急性氨氮胁迫研究中也发现AMS、胃蛋白酶活力在胁迫初期呈下降趋势[29]。然而在持续的氨氮胁迫下,鱼类为维持体内氨氮积累-代谢的平衡,需要消耗的能量远远高于机体的正常代谢水平,这也刺激了鱼类相应提高消化酶活力以增加能量供给,增强机体对氨氮环境的适应性。实验中氨氮胁迫第12小时后,AMS、胃蛋白酶、LPS呈上升趋势。这说明在持续的高浓度氨氮胁迫下,机体诱导消化酶活力升高,分解体内糖类、蛋白质、脂肪供应机体适应高氨氮环境所消耗的能量。随着氨氮胁迫时间的延长,氨氮不断以NH3的形式经鳃、表皮和肠黏膜等进入鱼体,胃、肠道等消化器官中的血氨含量不断升高,器官表层黏膜持续胁迫于高氨氮环境中可导致消化器官损伤,消化机能下降[30]。本实验中鳜鱼幼鱼胁迫至第96小时,胃囊AMS与肠道胃蛋白酶、LPS活力再次升高,说明持续氨氮胁迫条件下,鱼类通过代谢、排出的方式排出过多的外源性氨氮,降低血氨含量,能量需求大量增加。这可能是由于消化酶活力受机体调控“代偿”性升高,以维持机体的能量供给[31]。本研究还发现,在96 h胁迫过程中AMS活力先降低后升高的变化程度显著高于胃蛋白酶与LPS,这可能与糖类的代谢机制有关。这与李波[32]、常志成等[26]的研究结果相似,即AMS活力急剧升高可能是由于鱼类消化代谢糖类物质供应能量时氧卡系数最高,在供氧不足时鱼类诱导AMS活力急剧升高,以消化代谢体内糖类物质来供应能量。

  • 图  1   主要衣壳蛋白的相似性网络聚类图

    注:图中的点表示主要衣壳蛋白序列,绿点代表ICTV分类的微病毒科 (n=15),橙点表示牡蛎样本中鉴定为微病毒科的序列 (n=220),蓝点 (n=741) 表示上述两个来源的序列在nr库中最相似的序列,灰线表示连接两个点的边 (即两点序列间blastp打分值)。

    Figure  1.   Identity network of major capsid proteins

    Note: The dots in the figure indicate the major capsid protein sequences, among which the green dots (n=15) represent the Microviridae classified by ICTV; the orange dots (n=220) indicate sequences identified as Microviridae in oyster samples; the blue dots (n=741)  indicate the sequences from the two sources mentioned above which are most similar in the NR library; the gray line indicates the edge connecting the two dots (Based on the score values of blastp).

    图  2   主要衣壳蛋白进化树及对应的基因组结构图

    注:MAFFT 比对主要衣壳蛋白的氨基酸序列,TrimAl 序列对齐,用iqtree建树,itol可视化,最大似然数;自展值:1 000,自动选择最佳替代模型,临界值70%。进化树上的红色星形表示来自牡蛎样品的微病毒科序列,序列号为:GX1-198598、­­T4S1-854210、T4S1-22425、ML1-11067和HSd1-5344568;ID的背景色表示用Cherry预测的宿主类型;ID右边一列表示ICTV分类的亚科和属,以及未分类病毒的样品来源;右侧为基因组结构示意图,不同颜色表示不同的标志基因。

    Figure  2.   Phylogenetic tree of major capsid proteins and related schematic maps of genome structure

    Note: The amino acid sequences of the major coat proteins were aligned by using MAFFT, trimmed by using trimAL sequences, phylogenetic tree was drawn by using iqtree, visualized by itol, Maximum Likelihood Estimation tree; bootstrap replication: 1 000, automatic selection of the best alternative model, and site coverage cutoff is 70%. The red ID on the phylogenetic tree indicates the sequence from the oyster sample. The background color of ID represents the host type predicted by Cherry; the column on the right side of ID shows the subfamilies and genera of ICTV classification, as well as the sample sources of unclassified viruses; and on the right side is a schematic map of genome structure, with different colors indicating different marker genes.

    图  3   微病毒科基因组序列间的两两平均核苷酸相似性

    Figure  3.   ANI value among genomic sequences

    图  4   微病毒科噬菌体主要衣壳蛋白三维结构预测图

    Figure  4.   Three-dimensional structure prediction of replication proteins of oyster-related Microviridae

    图  5   微病毒科噬菌体主要衣壳蛋白与外部支架蛋白的进化关系

    注:使用 MAFFT 进行基因组序列比对,TrimAl 对齐氨基酸序列,iqtree 建树,itol可视化,最大似然数;自展值:1 000,替代模型:自动选择。图中紫色序列ID为Gokushovirinae序列,红色序列ID为牡蛎相关微病毒科序列,粉色序列ID为Bullavirinae序列。左边为主要衣壳蛋白的进化树,右边为外部支架蛋白的进化树,同一个基因组的两个蛋白用直线连接。

    Figure  5.   Linkage between phylogenetic tree of major coat protein and that of external scaffolding protein of Microviridae phage

    Note: The amino acid sequences of the major coat proteins were aligned by using MAFFT, trimmed by using trimAL sequences, phylogenetic tree was drawn by using iqtree, visualized by itol, Maximum Likelihood Estimation; bootstrap replication: 1 000, automatic selection of the best alternative model, and site coverage cutoff value is 70%. The purple ID in the figure is Gokushovirinae; the red ID is the oyster-associated Microviridae sequence, and the pink ID is the Bullavirinae sequence. The phylogenetic tree of the major capsid protein is shown on the left, and the evolutionary tree of the external scaffold protein is shown on the right, with two proteins of the same genome connected by straight lines.

    表  1   牡蛎相关微病毒在各测序文库中的 TPM 丰度

    Table  1   TPM abundance of oyster-related microviruses in sequencing libraries

    样本 ID     
    Sample ID     
    病毒基因组 ID Genome of virus ID
    HSd1-5344568GX1-198598T4S1-22425T4S1-854210ML1-11067
    GX170519 0 0.0 012 0 0 0
    S1-DR 0.0 006 0 0 0 0
    S2-D-2 0.0 010 0 0 0 0
    S2-DR 0.0 005 0 0 0 0
    S3-D-2 0.0 002 0 0 0 0
    S3-D 0.0 055 0 0 0 0
    样本 ID     
    Sample ID     
    病毒基因组 ID Genome of virus ID
    HSd1-5344568 GX1-198598 T4S1-22425 T4S1-854210 ML1-11067
    S3-DR 0.0 023 0 0 0 0
    S4-D-2 0.0 004 0 0 0 0
    S5-D-2 0.0 017 0 0 0 0
    S6-D-2 0.0 006 0 0 0 0
    S7-D-2 0.0 003 0 0 0 0
    T4S170523 0 0 0.000 8 0.000 7 0
    T5S170523 0 0 0.001 3 0 0
    下载: 导出CSV
  • [1] 李辉尚, 李坚明, 秦小明, 等. 中国牡蛎产业发展现状、问题与对策:基于鲁、闽、粤、桂四省区的实证分析[J]. 海洋科学, 2017, 41(11): 125-129.
    [2] 李晨, 谢晓晨, 王博, 等. 牡蛎细菌病的研究进展[J]. 环境生态学, 2022, 4(4): 59-64.
    [3]

    CHANG R Y, WONG J, MATHAI A, et al. Production of highly stable spray dried phage formulations for treatment of Pseudomonas aeruginosa lung infection[J]. Eur J Pharm Biopharm, 2017, 12(1): 1-13.

    [4] 赵虹泽. 噬菌体对鸡白痢治疗效果评价及其对盲肠菌群的影响[D]. 武汉: 华中农业大学, 2022: 1-10.
    [5]

    DOSS J, CULBERTSON K, HAHN D, et al. A review of phage therapy against bacterial pathogens of aquatic and terrestrial organisms[J]. Viruses, 2017, 9(3): 50-60. doi: 10.3390/v9030050

    [6]

    PEDULLA M L, FORD M E, HOUTZ J M, et al. Origins of highly mosaic mycobacteriophage genomes[J]. Cell, 2003, 113(2): 171-182. doi: 10.1016/S0092-8674(03)00233-2

    [7]

    SALMOND G P, FINERAN P C. A century of the phage: past, present and future[J]. Nat Rev Microbiol, 2015, 13(12): 777-786. doi: 10.1038/nrmicro3564

    [8]

    GILDEA L, AYARIGA J A, ROBERTSON B K. Bacteriophages as biocontrol agents in livestock food production[J]. Microorganisms, 2022, 10(11): 2126-2145. doi: 10.3390/microorganisms10112126

    [9]

    DESNUES C, RODRIGUEZ-BRITO B, RAYHAWK S, et al. Biodiversity and biogeography of phages in modern stromatolites and thrombolites[J]. Nature, 2008, 452(7185): 340-343. doi: 10.1038/nature06735

    [10]

    ROSARIO K, DAYARAM A, MARINOV M, et al. Diverse circular ssDNA viruses discovered in dragonflies (Odonata: Epiprocta)[J]. J Gen Virol, 2012, 93(Pt 12): 2668-2681.

    [11]

    JIANG J Z, FANG Y F, WEI H Y, et al. A remarkably diverse and well-organized virus community in a filter-feeding oyster[J]. Microbiome, 2023, 11(1): 2-16. doi: 10.1186/s40168-022-01431-8

    [12]

    BRENTLINGER K L, HAFENSTEIN S, NOVAK C R, et al. Microviridae, a family divided: isolation, characterization, and genome sequence of phiMH2K, a bacteriophage of the obligate intracellular parasitic bacterium Bdellovibrio bacteriovorus[J]. J Bacteriol, 2002, 184(4): 1089-1094. doi: 10.1128/jb.184.4.1089-1094.2002

    [13]

    QUAISER A, DUFRESNE A, BALLAUD F, et al. Diversity and comparative genomics of Microviridae in Sphagnum-dominated peatlands[J]. Front Microbiol, 2015, 6: 375-385.

    [14]

    ZHANG T, BREITBART M, LEE W H, et al. RNA viral community in human feces: prevalence of plant pathogenic viruses[J]. PLoS Biol, 2006, 4(1): 3-12.

    [15]

    BREITBART M, HAYNES M, KELLEY S, et al. Viral diversity and dynamics in an infant gut[J]. Res Microbiol, 2008, 159(5): 367-373. doi: 10.1016/j.resmic.2008.04.006

    [16]

    BENBOW R M, HUTCHISON C A, FABRICANT J D, et al. Genetic map of bacteriophage φX174[J]. J Virol, 1971, 7(5): 549-558. doi: 10.1128/jvi.7.5.549-558.1971

    [17]

    CHEN S, ZHOU Y, CHEN Y, et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): i884-i890. doi: 10.1093/bioinformatics/bty560

    [18]

    LI D, LIU C M, LUO R, et al. MEGAHIT: an ultra-fast single-node solution for large and complex metagenomics assembly via succinct de Bruijn graph[J]. Bioinformatics, 2015, 31(10): 1674-1676. doi: 10.1093/bioinformatics/btv033

    [19]

    BUCHFINK B, XIE C, HUSON D H. Fast and sensitive protein alignment using DIAMOND[J]. Nat Methods, 2015, 12(1): 59-60. doi: 10.1038/nmeth.3176

    [20]

    HUSON D H, BEIER S, FLADE I, et al. MEGAN community edition-interactive exploration and analysis of large-scale microbiome sequencing data[J]. PLoS Comput Biol, 2016, 12(6): e1004957. doi: 10.1371/journal.pcbi.1004957

    [21]

    HYATT D, CHEN G L, LOCASCIO P F, et al. Prodigal: prokaryotic gene recognition and translation initiation site identification[J]. BMC Bioinform, 2010, 11(1): 119-130. doi: 10.1186/1471-2105-11-119

    [22]

    关迎晖, 向勇, 陈康. 基于Gephi的可视分析方法研究与应用[J]. 电信科学, 2013,2 9(S1): 112-119.

    [23]

    KATOH K, STANDLEY D M. MAFFT multiple sequence alignment software version 7: improvements in performance and usability[J]. Mol Biol Evol, 2013, 30(4): 772-780. doi: 10.1093/molbev/mst010

    [24]

    CAPELLA-GUTIÉRREZ S, SILLA-MARTíNEZ J M, GABALDÓN T. TrimAl: a tool for automated alignment trimming in large-scale phylogenetic analyses[J]. Bioinformatics, 2009, 25(15): 1972-1973. doi: 10.1093/bioinformatics/btp348

    [25]

    NGUYEN L T, SCHMIDT H A, von HAESELER A, et al. IQ-TREE: a fast and effective stochastic algorithm for estimating maximum-likelihood phylogenies[J]. Mol Biol Evol, 2014, 32(1): 268-274.

    [26]

    LEE I, OUK KIM Y, PARK S C, et al. OrthoANI: an improved algorithm and software for calculating average nucleotide identity[J]. Int J Syst Evol Microbiol, 2016, 66(2): 1100-1103. doi: 10.1099/ijsem.0.000760

    [27]

    TISZA M J, BELFORD A K, DOMÍNGUEZ-HUERTA G, et al. Cenote-Taker 2 democratizes virus discovery and sequence annotation[J]. Virus Evol, 2021, 7(1): veaa100. doi: 10.1093/ve/veaa100

    [28]

    SHANG J, SUN Y. CHERRY: a computational method for accurate prediction of virus-prokaryotic interactions using a graph encoder-decoder model[J]. Brief Bioinform, 2022, 23(5): 182-198. doi: 10.1093/bib/bbac182

    [29]

    PATRO R, DUGGAL G, LOVE M I, et al. Salmon provides fast and bias-aware quantification of transcript expression[J]. Nat Methods, 2017, 14(4): 417-419. doi: 10.1038/nmeth.4197

    [30]

    HOPKINS M, KAILASAN S, COHEN A, et al. Diversity of environmental single-stranded DNA phages revealed by PCR amplification of the partial major capsid protein[J]. Isme J, 2014, 8(10): 2093-2103. doi: 10.1038/ismej.2014.43

    [31] 李灏, 丁子元, 徐林通, 等. 论我国水产养殖病害控制技术现状与发展趋势[J]. 农业与技术, 2015, 35(24): 171.
    [32] 陈愿. 噬菌体在水产养殖业中的研究进展[J]. 水产学报, 2021, 45(9): 1605-1615.
    [33]

    ROHWER F, THURBER R V. Viruses manipulate the marine environment[J]. Nature, 2009, 459(7244): 207-212. doi: 10.1038/nature08060

    [34]

    ROUX S, ENAULT F, ROBIN A, et al. Assessing the diversity and specificity of two freshwater viral communities through metagenomics[J]. PLoS One, 2012, 7(3): e33641. doi: 10.1371/journal.pone.0033641

    [35]

    KRUPOVIC M, FORTERRE P. Microviridae goes temperate: microvirus-related proviruses reside in the Genomes of Bacteroidetes[J]. PLoS One, 2011,6(5): e19893.

    [36] 李振, 张建城, 曹振辉, 等. 噬菌体控制主要水产养殖类致病菌的研究进展[J]. 畜牧与兽医, 2015, 47(8): 138-143.
    [37]

    PENG Y, JIN Y, LIN H, et al. Application of the VPp1 bacteriophage combined with a coupled enzyme system in the rapid detection of Vibrio parahaemolyticus[J]. J Microbiol Methods, 2014, 98(3): 99-104.

    [38]

    GAUDU P, YAMAMOTO Y, JENSEN P R, et al. Genetics of Lactococci[J]. Microbiol Spectr, 2019, 7(4): 361-362.

    [39] 洪宝华, 马荣荣, 袁娜. 养殖梭鱼格氏乳球菌的分离鉴定及致病性研究[J]. 农业生物技术学报, 2020, 28(8): 1458-1470.
    [40]

    ZHU P, LIU G F, LIU C, et al. Novel RNA viruses in oysters revealed by virome[J]. iMeta, 2022, 1(4): e65.

  • 期刊类型引用(9)

    1. 王恒杰,戴梦杨,王倩,熊信宇,王灿莉,袁向阳. 环境胁迫因子对鱼类健康影响的研究进展. 中国农学通报. 2025(02): 157-164 . 百度学术
    2. 罗寒雨,刘笑,郑菲菲,毛天宇,冯柯尧,孙艳红,祝东梅. 翘嘴鳜肝脏组织原代细胞培养技术的建立. 水产科学. 2024(05): 808-816 . 百度学术
    3. 陈旭 ,左涛 ,周胜杰 ,杨蕊 ,于刚 ,秦传新 ,马振华 . 饲喂甘草后低盐胁迫对尖吻鲈相关酶活性的影响. 水产科学. 2023(03): 404-412 . 百度学术
    4. 张成硕 ,赵岩 ,王艳玲 ,曾萌冬 ,赵金良 . 酸碱胁迫对鳜存活率、组织结构及解毒酶活性的影响. 水产科学. 2023(04): 640-647 . 百度学术
    5. 卫明亮,张志伟,张志勇,林志杰,祝斐,贾超峰,孟乾,徐大凤,张曹进. 冷应激对黑鲷组织损伤及细胞凋亡基因表达的影响. 南方水产科学. 2022(05): 110-117 . 本站查看
    6. 杨慧,胡华蓉. 华东地区活鱼运输应激反应研究及预防建议. 黑龙江水产. 2022(06): 68-69 . 百度学术
    7. 韩朝婕,陈屹洋,贺振楠,张严匀,周文礼,高金伟,贾旭颖. 氨氮胁迫对水产动物生长、消化酶及免疫影响的研究进展. 河北渔业. 2021(05): 32-35 . 百度学术
    8. 张三珊,刘海粟,林妙华,廖绍安,王安利,付胜利. 大蒜和茯苓对草鱼幼鱼肠道结构、功能及肠道微生物群落的影响. 饲料工业. 2021(22): 6-13 . 百度学术
    9. 刘炜,周国勤,裴雪莹,茆健强,陈树桥,张雷鸣,石晓兰. 运输密度对兴凯湖翘嘴鲌水箱水质的影响. 中国农学通报. 2021(35): 124-130 . 百度学术

    其他类型引用(7)

图(5)  /  表(1)
计量
  • 文章访问数:  306
  • HTML全文浏览量:  94
  • PDF下载量:  61
  • 被引次数: 16
出版历程
  • 收稿日期:  2023-03-06
  • 修回日期:  2023-06-04
  • 录用日期:  2023-06-25
  • 网络出版日期:  2023-06-29
  • 刊出日期:  2023-10-04

目录

/

返回文章
返回