Effects of nitrite and microplastic stress on immune, detoxification metabolism and osmoregulation-related indicators in gills of Litopenaeus vannamei
-
摘要:
鳃是对虾重要的呼吸器官，是亚硝酸盐毒性效应的主要靶器官，也是微塑料主要富集的部位之一。鳃组织参与了对虾渗透调节、氮排泄、免疫功能等生理过程，对对虾维持机体健康具有重要意义。为研究亚硝酸盐和微塑料单因素及复合胁迫对凡纳滨对虾 (Litopenaeus vannamei) 鳃组织生理功能的影响，将对虾分为对照组、20 mg·L−1亚硝酸盐胁迫组 (NIT)、10 µg·L−1微塑料胁迫组 (MP)、20 mg·L−1亚硝酸盐和10 µg·L−1微塑料复合胁迫组 (NM)，于第14天测定对虾鳃中相关指标变化。结果显示：1) 氧化应激指标丙二醛 (MDA)、过氧化氢 (H2O2) 浓度及抗氧化酶超氧化物歧化酶 (SOD)、过氧化氢酶 (CAT)、谷胱甘肽过氧化物酶 (GPx) 活性在胁迫后发生不同程度的改变；2) 解毒代谢指标细胞色素P450基因 (CYP450) 和谷胱甘肽S-转移酶基因 (GST) 相对表达水平在胁迫后发生不同程度的紊乱；3) 渗透调节指标如离子转运酶液泡型ATP酶 (VATP)、Na+/H+交换体7 (NHE7)、Na+/K+-ATP酶亚基α (NKA-α)、Na+/K+-ATP酶亚基β (NKA-β)、碳酸酐酶 (CA) 和水通道蛋白TIP4-1 (TIP4)、钙通道蛋白1 (CCP)、氯通道蛋白2 (CLC)、水通道蛋白 (AQP) 基因的相对表达水平在胁迫后发生不同程度的紊乱；4) 细胞凋亡因子基因 (CASP-3) 相对表达水平在3个胁迫组均显著降低 (P<0.05)。由此推断，亚硝酸盐和微塑料胁迫会诱导凡纳滨对虾鳃中免疫、解毒代谢和渗透调节相关指标变化，进而影响其正常的生理功能。
Abstract:Gill is an important respiratory organ of prawns, the main target organ of nitrite toxicity effect, and also one of the main enrichment sites of microplastics. The gill tissue participates in the physiological processes of prawns such as osmoregulation, nitrogen excretion, immune function, etc., which is important for maintaining the prawns' health. In order to investigate the effects of single and combined stress of nitrite and microplastics on the physiological functions of gill tissues of Litopenaeus vannamei, we designed the control group, 20 mg·L−1 nitrite stress group (NIT), 10 μg·L−1 microplastic stress group (MP), 20 mg·L−1 nitrite and 10 μg·L−1 microplastic composite stress group (NM), and then measured the changes of immune and osmotic regulation in gills of shrimps on the 14th day. The results show that: 1) The oxidative stress indicators such as the contents of malondialdehyde (MDA) and hydrogen peroxide (H2O2), and the activities of superoxide dismutase (SOD), catalase (CAT) and glutathione peroxidase (GPx) changed at different degrees after the stress. 2) The detoxification metabolic indicators such as the relative expression levels of cytochrome P450 (CYP450), glutathione S-transferase (GST) and apoptosis factor (CASP-3) were disturbed at different degrees after the stress. 3) The relative gene expression levels of the osmoregulation indexes such as ion transporters (VATP, NHE7, NKA-α, NKA-β and CA) and channel proteins (TIP4, CCP, CLC and AQP) occurred at different degrees of disorder after the stress. 4) The apoptosis index sucha as the relative expression level of apoptosis factor (CASP-3) gene decreased significantly in the three groups (P<0.05). Thus, it is inferred that nitrite and microplastic stress can induce the changes of immunity, detoxification metabolism and osmotic regulation in the gills of L. vannamei, affecting its normal physiological functions.
-
Keywords:
- Litopenaeus vannamei /
- Gill /
- Nitrite /
- Microplastics /
- Detoxification metabolism /
- Osmotic regulation
-
牡蛎 (Crassostrea) 是一种双壳类软体动物，广泛分布于世界各地沿岸潮间带，其肉质鲜美、营养丰富，是世界第一大养殖贝类，也是我国支柱性海水养殖种类之一[1]。由于养殖环境恶化、养殖密度过高和种质退化等原因，病害问题已经成为制约牡蛎产业发展的重要因素。其中细菌性病原的危害尤其严重[2]。细菌病害的发生与养殖密度和水质环境密切相关，但由于牡蛎大多在开放海域养殖，直接施用消毒剂、抗生素等药物会被海水快速稀释，不但起不到防治作用，还会对海洋生态环境造成破坏。
噬菌体是专性侵染细菌的病毒，能高效、精准地攻击宿主细菌，在海洋生态系统中扮演着重要角色。使用噬菌体治疗细菌性疾病已在医疗[3]、畜禽[4]和水产[5]领域得到了广泛实践，是减抗替抗的重要选择。虽然细菌变异很快，但噬菌体拥有高度多样化的基因库和极高的变异速度[6]，能快速获得新功能从而适应宿主及环境变化[7-8]，实现控制宿主细菌丰度、微生物群落组成，并影响动物宿主健康，甚至对海洋生物地球化学循环产生影响[9]。迄今为止，对具有双链DNA (double strands DNA, dsDNA) 基因组的噬菌体的研究较为充分，已经有大量关于它们的多样性、进化及在细菌性病害治疗方面的研究报道。然而，针对单链DNA (single strand DNA, ssDNA) 噬菌体的研究仍然非常缺乏。与dsDNA噬菌体相比，ssDNA噬菌体的基因组较小，且在各种环境中广泛存在[10]，便于直接人工合成，通常不携带毒力或抗性基因等有害基因，是进行功能噬菌体开发的理想类群。
本课题组的前期研究发现，养殖牡蛎体内蕴藏着异常丰富、多样且全新的病毒和噬菌体类群，其中大量与微病毒科 (Microviridae) 和圆环病毒科 (Circoviridae) 等ssDNA病毒相关[11]。微病毒科噬菌体是一类具有环状ssDNA基因组的噬菌体[12]，在水环境、动物和人类内脏的内容物中均有检出[13-15]。微病毒科包含两个经由国际病毒分类委员会 (International Committee of Taxonomy of Viruses, ICTV) 批准的亚科：Bullavirinae 和Gokushovirinae，两个亚科之间的区别主要在于基因组组成和主要衣壳蛋白的结构不同。Bullavirinae以微病毒PhiX174为代表，包括Alphatrevirus、Gequatrovirus和Sinsheimervirus 3个属，据报道，该亚科可感染肠杆菌属 (Enterobacter) 细菌，并已通过其代表病毒株PhiX174进行了深入的分子学研究[16]。Gokushovirinae亚科因其可感染细胞内寄生生物衣原体属 (Chlamydia)、螺原体属(Spiroplasma) 和蛭弧菌属 (Bdellovibrio)[12]，通常被称为“寄生虫的寄生虫”。
近年来，在多种宏基因组样本中也发现了大量未分类的微病毒科噬菌体，或许会进一步扩充微病毒的宿主范围[13-15]。本研究对华南沿海多地养殖牡蛎样本的宏病毒组测序数据进行深入挖掘，发现并鉴定到5个感染埃希氏菌属 (Escherichia) 的新型微病毒科噬菌体，相关成果可为推动微病毒科噬菌体多样性与分类研究、研发具有良好抑菌作用的新型ssDNA噬菌体制剂提供参考。
1. 材料与方法
1.1 数据来源
本课题组已完成华南多地养殖牡蛎样品的宏病毒组测序工作，共获得了约25亿条测序原始序列 (reads)[11]。使用 Fastp (V0.20.0)[17] 对原始序列进行质控，包括去除低质量和接头序列，然后使用Megahit (V1.2.9)[18]将过滤后的reads组装成重叠群 (Contig)，进一步使用Diamond (V0.9.24.125)[19]比对美国国家生物技术信息中心 (National Center for Biotechnology Information, NCBI) 的非冗余蛋白质数据库 (NR)，从而对Contig序列进行注释。使用Megan 6[20]对注释结果进行进一步分类鉴定，选择其中注释为微病毒科的病毒序列，经病毒基因组完整性鉴定后共获得了227个几乎完整的疑似微病毒基因组。
1.2 微病毒序列筛选
选用微病毒科中经ICTV分类的8个属 (分属于两个亚科) 的15条序列作为参考株，与新获得的227条微病毒序列一并分析。使用Prodigal (V2.6.3)[21]预测开放阅读框 (Open reading frame, ORF)，run type选择normal mode，mode选取meta。使用Diamond (V0.9.24.125) 对预测的ORF进行Blastp比对，比对数据库为NCBI非冗余蛋白质数据库。筛选出其中比对结果为主要衣壳蛋白且相似度最高的前10条序列，除去重复序列后，共得到741条相似度较高的主要衣壳蛋白序列。使用gephi (V0.9.2)[22]绘制主要衣壳蛋白的网络关系图，选取与Bullavirinae聚类的5条序列进行深入分析 (Contig ID: HSd1-5344568、GX1-198598、T4S1-22425、T4S1-854210、ML1-11067)。
1.3 主要衣壳蛋白进化树分析
使用MAFFT (V7.515)[23]将所选取的序列进行比对，并使用TrimAl[24]对结果进行裁剪。采用iqtree (V2.2.0.3)[25]进行最大似然系统发育分析，设置为自动选择最优模型，Bootstrap为1 000，最后使用itol (https://itol.embl.de/) 可视化树文件。
1.4 基因组结构分析
使用Orthologous Average Nucleotide Identity Tool (OAT) (V0.93.1)[26]计算5条病毒序列与其在进化树上的相邻序列之间的平均核苷酸相似性 (Average Nucleotide Identity, ANI)。使用Cenote Taker2 (V2.1.3)[27]对有代表性的若干基因组的ORF进行功能注释，并使用SnapGene (V4.3.6, www.snapgene.com) 用于打开Cenote-Taker2输出文件，将基因组结构可视化，展示在进化树中对应的序列之后。
1.5 宿主预测
使用CHERRY (V0.9.24.125)[28]对进化树上的组装病毒序列进行宿主预测。CHERRY集成了多种基于蛋白质和 DNA的序列特征 (如噬菌体之间的蛋白质信息、噬菌体与原核生物之间的序列相似性、CRISPR 信号以及k-mer频率信息)，利用这些特征计算配对病毒和原核生物之间形成感染的概率，选择概率最高的为预测的细菌宿主。
1.6 病毒丰度分析
为了计算每个病毒的相对丰度值，首先用salmon (V0.13.1) index[29] 命令对合并后的5条病毒基因组序列构建索引 (Index)，作为salmon比对的参考基因组数据集；然后再利用salmon quant[29] 命令将牡蛎病毒组文库全部的clean reads与参考基因组进行映射 (Mapping)，统计能映射到这5条基因组的reads数量，根据调整后的TPM (Transcripts per million) 计算公式计算这5条病毒序列的相对丰度值。
2. 结果
2.1 牡蛎微病毒科相关病毒的发现与聚类
本研究在前期工作的基础上，以中国华南沿海多地养殖牡蛎为研究对象，进行宏病毒组测序后，去除低质量序列并组装，得到3 347 421条非冗余的Contig序列，并与NCBI非冗余蛋白序列库进行blastp比对和分类鉴定，最终获得728 784条疑似病毒来源的Contig。由图1可见，牡蛎样本中 (橙色点) 发现了大量的疑似微病毒科噬菌体，其中绝大多数为全新的分类群，与NCBI中未分类的相关噬菌体聚类，并与ICTV参考毒株 (绿色点) 差异明显。由于未分类病毒类群的分析难度较高，本研究进一步筛选其中与Bullavirinae参考株序列聚类在一起的5条疑似微小病毒的基因组序列 (图1，橙色圈) 进一步分析。
![]()
图 1 主要衣壳蛋白的相似性网络聚类图注：图中的点表示主要衣壳蛋白序列，绿点代表ICTV分类的微病毒科 (n=15)，橙点表示牡蛎样本中鉴定为微病毒科的序列 (n=220)，蓝点 (n=741) 表示上述两个来源的序列在nr库中最相似的序列，灰线表示连接两个点的边 (即两点序列间blastp打分值)。Figure 1. Identity network of major capsid proteinsNote: 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.2 微病毒科成员的系统发育、宿主来源与基因组结构
主要衣壳蛋白是微病毒科分类的经典系统发育标志基因[30]。为了解这5条病毒序列的进化关系和分类界定，对5条序列、与其相关的近缘病毒及经ICTV认定的微病毒科参考株序列的主要衣壳蛋白序列汇总构建进化树 (图2)。进化树分成3个大的分支，其中HSd1-5344568与其他4条序列聚类到不同的分支中，与已被报道的序列QBQ83335.1亲缘关系最近，聚集在亚科Bullavirinae的大分支。其余4条序列聚类在Bullavirinae和Gokushovirinae两个亚科分支的中间分支中，其中T4S1-22425、T4S1-854210和 GX1-198598这3条序列聚类在同一个小分支中，且没有其他已知的相似序列。利用Cherry软件进行宿主预测的结果显示，Gokushovirinae分支的宿主主要为螺原体、蛭弧菌和衣原体属。Bullavirinae分支以及新发现的分支中大多数噬菌体的宿主为埃希氏菌属、肠杆菌属和乳球菌属 (Lactococcus) ，符合已报道的宿主范围，并且本研究的5条噬菌体基因组均来自埃希氏菌属 (棕色标注)。
![]()
图 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 structureNote: 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.由样品来源看 (图2，中间注释信息)，未知病毒数据主要来源于动物、废水和沉积物共3类样品。但在5个牡蛎病毒各自的进化分支上，同一分支的相邻序列也均来自于动物宏基因组样本，说明这些微病毒的宿主可能是在动物宿主体内专性共生的细菌。
根据基因组功能基因注释 (图2) 发现：首先，微病毒科噬菌体具有比较保守的基因组结构。主要衣壳蛋白 (黄色) 和复制相关蛋白 (粉色) 是其最具标志性的两个保守基因，两个基因并不相邻，中间会夹杂一些其他基因。其次，微病毒科的编码基因均为同一方向，即均为正义链编码，与圆环病毒等其他环状ssDNA病毒不同。最后，微病毒科3个主要亚科之间存在差异明显的基因组结构。位于主要衣壳蛋白和复制相关蛋白之间的次要衣壳蛋白 (橙色) 是Gokushovirinae的特征基因；两种刺突蛋白 (绿色和深蓝色) 和外部支架蛋白 (浅蓝色) 是Bullavirinae的特征基因；而作为Bullavirinae姊妹分支的中间分支同样拥有外部支架蛋白，但刺突蛋白似乎并不是其特征基因。值得一提的是，T4S1-22425基因组中未能预测到复制相关蛋白的基因，这或许说明该病毒拥有非常独特的复制相关蛋白，当然也无法排除因测序或组装导致的序列错误。本研究还发现，归属于Bullavirinae的HSd1-5344568与相似序列QBQ83335.1的基因组结构较为一致，但不同于ICTV标准Bullavirinae毒株，它们都缺少次要刺突蛋白。
2.3 平均核苷酸一致性
由图2可知，GX1-198598、T4S1-22425和T4S1-854210 3条序列的主要衣壳蛋白聚类在一起，HSd1-5344568与QBQ83335.1、ML1-11067与AXL15490.1分别在同一个聚类中。为了进一步比较几个近缘病毒之间的相似性，进一步分析了这些病毒基因组间的两两平均核苷酸相似性 (ANI) (图3)。GX1-198598、T4S1-22425和T4S1-854210 3个基因组之间的ANI为0。HSd1-5344568和ML1-11067与各自最相近序列的ANI分别为66.23%和62.09%。综合进化树的聚类和ANI的结果，推测5条序列中HSd1-5344568属于Bullavirinae亚科的一个新属，其余4条序列则属于一个新的亚科。ICTV的分类标准判断，这5条序列与任一参考微病毒科病毒序列的平均核苷酸相似性均低于70%，因此至少应该分别代表了5个新的属。
2.4 主要衣壳蛋白结构预测
通过结构预测结果，可以辅助判断证实蛋白的功能及活性。从本研究对5个新型病毒的主要衣壳蛋白结构域三维结构预测结果可见 (图4)，5条序列的主要衣壳蛋白的空间结构存在较高的相似性，所有病毒的主要衣壳蛋白均具有保守的8股β桶核心 (也称为病毒果冻卷)。其中图4-a、4-b和4-c为进化树上相邻的3个分支，这3个序列的三维结构也非常相似。HSd1-5344568 (图4-d)与其相似序列QBQ83335.1 (图4-e)、ML1-11067 (图4-f)与其相似序列AXL15490.1 (图4-g) 之间的空间结构相似性均较高。
![]()
图 4 微病毒科噬菌体主要衣壳蛋白三维结构预测图Figure 4. Three-dimensional structure prediction of replication proteins of oyster-related Microviridae2.5 主要衣壳蛋白进化树与外部支架蛋白的进化规律与重组
为了更好地了解微病毒中各基因的进化规律，将主要衣壳蛋白进化树与外部支架蛋白的进化树进行了关联分析。由图5所示，无论是外部支架蛋白还是主要衣壳蛋白，Gokushovirinae序列都聚集在一个独立的进化分支中，说明Gokushovirinae亚科的进化位置或生态位相对独立，未发现与其他亚科病毒的重组现象。另一方面，虽然外部支架蛋白进化树的总体结构与主要衣壳蛋白相似，即Gokushovirinae、Bullavirinae和中间分支相互独立 (图5中间黑色和蓝色连线)，但发现其中有两个病毒 (NC001330和KF044309) 的基因可能发生了重组 (图5中间红色连线)，即重组了属于Bullavirinae的主要衣壳蛋白和属于中间分支的外部支架蛋白。值得注意的是，KF044309与其姊妹分支KX266637的宿主均为肠杆菌属 (图2)，与临近的其他病毒宿主 (埃希氏菌属) 明显不同。该结果提示微病毒的中间分支类群与Bullavirinae类群可能在一定程度上存在宿主或生态位重叠，从而为两种病毒基因组的重组提供了机会。
![]()
图 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 phageNote: 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.2.6 牡蛎相关微病毒的丰度
分析病毒序列在测序数据中的丰度可为了解病毒的生态功能提供线索，噬菌体的丰度与其细菌宿主的丰度密切相关。为了计算每个噬菌体序列的相对丰度，统计了Salmon软件映射到5条基因组上的reads数量，其中只有除ML1-11067外的4条序列有映射结果。由TPM丰度表 (表1) 可见，4条序列在文库中的丰度均比较低，说明这些病毒的宿主在牡蛎体内的丰度并不高，对宿主的影响较小。HSd1-5344568在样品中的丰度最高、分布最广。其中ML1-11067在各个文库中都没有映射到reads，可能是组装软件Meghit与映射软件Salmon的算法不同所导致。
表 1 牡蛎相关微病毒在各测序文库中的 TPM 丰度Table 1. TPM abundance of oyster-related microviruses in sequencing libraries样本 ID
Sample ID病毒基因组 ID Genome of virus ID HSd1-5344568 GX1-198598 T4S1-22425 T4S1-854210 ML1-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 3. 讨论
迄今为止，大多数水生生物的微生物组研究都集中在细菌群落上，其病毒组直到最近才受到关注[31]。噬菌体与微生物群落和动物宿主都有很强的相关性[32]，在肠道等复杂的生态系统中，多种因素可能会影响噬菌体的丰度、多样性和生态作用。然而，由于噬菌体基因组中缺乏通用的基因标记，因此很难鉴定噬菌体的群落组成[33]。病毒宏基因组学对给定样本中的全体病毒核酸进行高通量测序，是探索新型病毒、挖掘病毒多样性的一种高效方法[15]。
最近的研究表明，人们已在多种环境中发现了众多微病毒科的成员，提示它们可能在不同的环境中发挥作用。如在不同自然栖息地或动物粪便和组织的宏基因组数据中，已经组装出数千个完整的微病毒基因组。这些环境微病毒都揭示了Gokushovirinae的巨大多样性，为将新发现的大量病毒进行合理分类，研究人员提出了一些假定亚科，其中包括从泥炭地生态系统中发现的 Aravirinae和Stokavirinae[13]、从白蚁肠道中发现的Sukshmavirinae[34]、蜻蜓相关的D 类群[10]以及从人类肠道中发现的Alpavirinae[35]。本研究结果显示，5条新型微病毒序列中，仅HSd1-5344568被分类为Bullavirinae，其余4条候选序列 (GX1-198598、T4S1-22425、T4S1-854210、ML1-11067) 虽然与ICTV认定的两个标准亚科没有聚类到一起，但明显属于介于这两个标准亚科的中间类群，与之前研究报道的其他假定类群不同。
近年来，养殖贝类的细菌性疾病连年暴发，严重阻碍了贝类养殖业的发展。目前的主流抗菌手段是使用抗生素，但其大规模且不规范的使用也带来了水产品药物残留、细菌耐药性和环境污染等突出问题，因此亟需寻求安全、天然无残留的抗生素替代品。噬菌体作为一种特异性高且生态友好的抗生素替代品引起了广泛关注，有实验表明噬菌体可以有效预防和控制水产养殖中的病原菌[32,36]。然而由于噬菌体具有宿主专一性，导致其抗菌谱较狭窄，但这一问题可以通过使用噬菌体鸡尾酒法来解决[37]，即使用3种或以上的噬菌体联合应用来杀灭细菌。所以在水生环境和水生动物样品中挖掘噬菌体的多样性对研发噬菌体抗菌制剂十分重要。
通过宿主预测软件的分析发现，这5个噬菌体的宿主可能为埃希氏菌属细菌，与先前报道的Bullavirinae噬菌体的宿主范围一致。埃希氏菌属属于肠杆菌科，包括大肠埃希氏菌 (E. coli)、费格森埃希氏菌 (E. fergusonii) 等5个种，其中致病性大肠埃希氏菌被认为是21世纪对人类健康产生重大影响的病原微生物之一[12]。Bullavirinae亚科的噬菌体还可能感染乳球菌属 (Lactococcus) 和肠杆菌属。乳球菌属包括格氏乳球菌 (L. garvieae)、乳酸乳球菌 (L. lactis) 等种或亚种，其中格氏乳球菌和鱼乳球菌 (L. piscium) 已被报道可引起鱼类、中华鳖 (Trionyx sinensis)、罗氏沼虾 (Macrobrachium rosenbergii) 等水产动物发生腹水性肝脾病变、腹腔和内脏出血以及体表溃疡等病症，是危害淡水养殖的常见病原菌[38-39]。肠杆菌属也是肠杆菌科中的一员，广泛分布于自然环境中，目前已发现该属的14个种，属于重要的机会性致病菌。根据中国细菌耐药网提供的数据，2015—2017年肠杆菌属检出率均排在前十位，尤其是该菌属下的阴沟肠杆菌检出率均占所检出革兰氏阴性菌第五位。
综上所述，本研究围绕5条与ICTV标准病毒株近缘的新型微病毒的全基因组序列，通过对包括宿主来源、基因组结构、主要衣壳蛋白、外部支架蛋白等多方面的分析，明确了5个微病毒的宿主来源及其分类位置，并发现了微病毒主要衣壳蛋白和外部支架蛋白存在的重组现象。本研究可为拓展海洋病毒多样性[40]、理解微病毒科病毒的分类进化规律、深入挖掘海洋噬菌体资源提供参考。
-
![]()
图 1 亚硝酸盐和微塑料胁迫对凡纳滨对虾鳃中氧化损伤指标的影响
注：数据上标不同字母表示组间差异显著 (P<0.05)；后图同此。
Figure 1. Effects of nitrite and microplastics stress on oxidative damage index in gills of L. vannamei
Note: Different letters on the bar indicate significant differences among the groups (P<0.05). The same case in the following figures.
![]()
图 2 亚硝酸盐和微塑料胁迫对凡纳滨对虾鳃中抗氧化酶活性的影响
Figure 2. Effects of nitrite and microplastics stress on antioxidant enzymes activity in gills of L. vannamei
![]()
图 3 亚硝酸盐和微塑料胁迫对凡纳滨对虾鳃组织免疫基因表达的影响
Figure 3. Effects of nitrite and microplastics stress on immune gene expression in gills of L. vannamei
![]()
图 4 亚硝酸盐和微塑料胁迫对凡纳滨对虾鳃组织渗透调节基因表达的影响
Figure 4. Effects of nitrite and microplastics stress on osmotic regulation gene expression in gills of L. vannamei
表 1 本研究所用引物序列
Table 1 Primer sequences used in this study
基因名称
Gene name正向引物 (5'—3')
Forward primer (5'–3')反向引物 (5'—3')
Reverse primer (5'–3')CASP-3 ATTATAGAGTGCCTGCGTTGCTACC TGCGATATCATGCCTAGATGCTTGG CYP450 CGCCGCCAAGGAAGTGAAGG TGGGAGGCGATGGTGGTGTC GST ACGAACACTACGAACAGAAGGATGC GCCAGGAAGTCGATGTAGGTTAGC VATP GCCAAGTATGAGCACCTAGCAGAC TCGTCCACCACCTTACTGAAGAGAG TIP4 GTCTCTCTGAACCCTTGCCAAACTC GTTGCCTATGATTACCTGCGTCCTC CCP CCCTGAGACCTCTGTTTGCTGTG TGGCAGTGTTCTTCGCATGTTCC NHE7 CACCGCCATCCTCACCAAGTTC ACAGAAGAGCACAGCCACAATACC CLC GCCGATAGGAGGTGTGTTGTTCAG CGAAGAAGCCACGCCAGTAGTTC NKA-α TGAAATCGTGTTTGCCCGTACCTC ACCATCACCAGTTACAGCCACAATG NKA-β GAACCCAGCCGACGAAGAATACG CAGCAACAATAGGTGGCAGGTAGC AQP GCAGCCATCTTGAAGGGAGTGAC ACGAGGACGAAGGTGATGAGGAG CA CCTATTCTGGCTCCCTCACTACCC GTTCATCCTCTGGGCAACACTCG β-actin AGCTCATTGTAGAAGGTGTGATGCC TCCTGACCCTGAAGTACCCCATTG 
下载: 导出CSV
-
[1] DUAN Y F, ZHANG Y, DONG H B, et al. Effect of dietary Clostridium butyricum on growth, intestine health status and resistance to ammonia stress in Pacific white shrimp Litopenaeus vannamei[J]. Fish Shellfish Immunol, 2017, 65: 25-33. doi: 10.1016/j.fsi.2017.03.048
[2] HUANG W T, YIN H, YANG Y Y, et al. Influence of the co-exposure of microplastics and tetrabromobisphenol A on human gut: simulation in vitro with human cell Caco-2 and gut microbiota[J]. Sci Total Environ, 2021, 778: 146264. doi: 10.1016/j.scitotenv.2021.146264
[3] WANG W N, WANG A L, ZHANG Y J, et al. Effects of nitrite on lethal and immune response of Macrobrachium nipponense[J]. Aquaculture, 2004, 232(1/2/3/4): 679-686.
[4] DUAN Y, LIU Q S, WANG Y, et al. Impairment of the intestine barrier function in Litopenaeus vannamei exposed to ammonia and nitrite stress[J]. Fish Shellfish Immunol, 2018, 78: 279-288. doi: 10.1016/j.fsi.2018.04.050
[5] HSIEH S L, WU Y C, XU R Q, et al. Effect of polyethylene microplastics on oxidative stress and histopathology damages in Litopenaeus vannamei[J]. Environ Pollut, 2021, 288: 117800. doi: 10.1016/j.envpol.2021.117800
[6] SABOROWSKI R, KOREZ P, RIESBECK S, et al. Shrimp and microplastics: a case study with the Atlantic ditch shrimp Palaemon varians[J]. Ecotox Environ Safe, 2022, 234: 113394. doi: 10.1016/j.ecoenv.2022.113394
[7] RAGUSA A, SVELATO A, SANTACROCE C, et al. Plasticenta: first evidence of microplastics in human placenta[J]. Environ Int, 2021, 146: 106274. doi: 10.1016/j.envint.2020.106274
[8] WANG D L, DI Z, WANG L M, et al. Effects of white spot syndrome virus infection on immuno-enzyme activities and ultrastructure in gills of Cherax quadricarinatus[J]. Fish Shellfish Immunol, 2012, 32(5): 645-650. doi: 10.1016/j.fsi.2012.01.005
[9] FOGUESATTO K, BOYLE R T, ROVANI M T, et al. Aquaporin in different moult stages of a freshwater decapod crustacean: expression and participation in muscle hydration control[J]. Comp Biochem Physiol A, 2017, 208: 61-69. doi: 10.1016/j.cbpa.2017.03.003
[10] LINGREL J B. Na, K-ATPase: isoform structure, function, and expression[J]. J Bioenerg Biomembr, 1992, 24(3): 263-270.
[11] 陈亭君. 低盐、高亚硝酸盐胁迫日本囊对虾转录组分析及MjTPS基因克隆和功能研究[D]. 湛江: 广东海洋大学, 2022: 64-67. [12] DUAN Y F, XIONG D L, WANG Y, et al. Toxicological effects of microplastics in Litopenaeus vannamei as indicated by an integrated microbiome, proteomic and metabolomic approach[J]. Sci Total Environ, 2020, 761(2): 143311.
[13] YAN M, LI W, CHEN X F, et al. A preliminary study of the association between colonization of microorganism on microplastics and intestinal microbiota in shrimp under natural conditions[J]. J Hazard Mater, 2021, 408: 124882. doi: 10.1016/j.jhazmat.2020.124882
[14] 于萍. 微塑料对中华绒螯蟹毒性效应的初步研究[D]. 上海: 华东师范大学, 2019: 8-22. [15] XIAO J, LUO S S, DU J H, et al. Transcriptomic analysis of gills in nitrite-tolerant and -sensitive families of Litopenaeus vannamei[J]. Comp Biochem Physiol C, 2022, 253: 109212.
[16] 周泽妍. 白洋淀入淀河流微塑料分布及其重金属吸附和风险评估[D]. 北京: 北京林业大学, 2021: 54-57. [17] GRAY A D, WEINSTEIN J E. Size- and shape-dependent effects of microplastic particles on adult daggerblade grass shrimp (Palaemonetes pugio)[J]. Environ Toxicol Chem, 2017, 36(11): 3074-3080. doi: 10.1002/etc.3881
[18] 林艳. 水体亚硝酸盐氮胁迫对鳙生理及肝鳃GSH/GSTs、HIF-1α通路关键因子表达的影响[D]. 南京: 南京农业大学, 2016: 16-22. [19] KOONGOLLA J B, LIN L, PAN Y F, et al. Occurrence of microplastics in gastrointestinal tracts and gills of fish from Beibu Gulf, South China Sea[J]. Environ Pollut, 2020, 258: 113734. doi: 10.1016/j.envpol.2019.113734
[20] 熊大林, 段亚飞, 陈成勋, 等. 高温与氨氮复合胁迫对凡纳滨对虾渗透调节的影响[J]. 水产科学, 2021, 40(4): 475-482. doi: 10.16378/j.cnki.1003-1111.19281 [21] GOPALAKRISHNAN S, CHEN F Y, THILAGAM H, et al. Modulation of immune-associated parameters and antioxidant responses in the crab (Scylla serrata) exposed to mercury[J]. Chemosphere, 2013, 90(3): 917-928. doi: 10.1016/j.chemosphere.2012.06.031
[22] 苏家齐, 祝华萍, 朱长波, 等. 盐度和钠离子/钾离子对凡纳滨对虾幼虾存活与组织结构的影响[J]. 南方水产科学, 2021, 17(5): 45-53. doi: 10.12131/20210011 [23] 王佳雯, 卢洁, 姚托, 等. 环境因子和免疫刺激对香港牡蛎免疫指标的影响[J]. 南方水产科学, 2021, 17(4): 18-26. doi: 10.12131/20210051 [24] LIU J, QU W, KADIISKA M B. Role of oxidative stress in cadmium toxicity and carcinogenesis[J]. Toxicol Appl Pharm, 2009, 238(3): 209-214. doi: 10.1016/j.taap.2009.01.029
[25] ARUN S, SUBRAMANIAN P. Antioxidant enzymes in freshwater prawn Macrobrachium malcolmsonii during embryonic and larval development[J]. Comp Biochem Phys B, 1998, 121(3): 273-277.
[26] YU P, LIU Z Q, WU D L, et al. Accumulation of polystyrene microplastics in juvenile Eriocheir sinensis and oxidative stress effects in the liver[J]. Aquat Toxicol, 2018, 200: 28-36. doi: 10.1016/j.aquatox.2018.04.015
[27] BURKINA V, ZAMARATSKAIA G, RANDAK T, et al. Verapamil does not modify catalytic activity of CYP450 in rainbow trout after long-term exposure[J]. Ecotox Environ Safe, 2012, 79(9): 148-152.
[28] GUO H, XIAN J A, WANG A L, et al. Analysis of digital gene expression profiling in hemocytes of white shrimp Litopenaeus vannamei under nitrite stress[J]. Fish Shellfish Immunol, 2016, 56: 1-11. doi: 10.1016/j.fsi.2016.06.059
[29] ZHANG W Y, LI J J, CHEN Y M, et al. Exposure time relevance of response to nitrite exposure: insight from transcriptional responses of immune and antioxidant defense in the crayfish, Procambarus clarkii[J]. Aquat Toxicol, 2019, 214: 105262. doi: 10.1016/j.aquatox.2019.105262
[30] HANSLIK L, HUPPERTSBERG S, KÄMMER N, et al. Rethinking the relevance of microplastics as vector for anthropogenic contaminants: adsorption of toxicants to microplastics during exposure in a highly polluted stream - analytical quantification and assessment of toxic effects in zebrafish (Danio rerio)[J]. Sci Total Environ, 2022, 816: 151640. doi: 10.1016/j.scitotenv.2021.151640
[31] MARCO G D, CONTI G O, GIANNETTO A, et al. Embryotoxicity of polystyrene microplastics in zebrafish Danio rerio[J]. Environ Res, 2022, 208: 112552. doi: 10.1016/j.envres.2021.112552
[32] CHANG C C, YEH M S, LIN H K, et al. The effect of Vibrio alginolyticus infection on caspase-3 expression and activity in white shrimp Litopenaeus vannamei[J]. Fish Shellfish Immunl, 2008, 25(5): 672-678. doi: 10.1016/j.fsi.2008.09.004
[33] SANTOS D, LUZIO A, FÉLIX L, et al. Oxidative stress, apoptosis and serotonergic system changes in zebrafish (Danio rerio) gills after long-term exposure to microplastics and copper[J]. Comp Biochem Physiol C, 2022, 258: 109363.
[34] FALEIROS R O, GOLDMAN M H S, FURRIEL R P M, et al. Differential adjustment in gill Na+/K+- and V-ATPase activities and transporter mRNA expression during osmoregulatory acclimation in the cinnamon shrimp Macrobrachium amazonicum (Decapoda, Palaemonidae)[J]. J Exp Biol, 2010, 213(22): 3894-3905. doi: 10.1242/jeb.046870
[35] SÁEZ A G, LOZANO E, ZALDÍVAR R A. Evolutionary history of Na, K-ATPases and their osmoregulatory role[J]. Genetica, 2009, 136: 479-490. doi: 10.1007/s10709-009-9356-0
[36] 李丽. 花鳗鲡幼鳗VHA、NKA及FXYD11基因适应盐度变化的分子调节机制[D]. 南京: 南京师范大学, 2015: 36-41. [37] ONKEN H, PUTZENLECHNER M. A V-ATPase drives active, electrogenic and Na+-independent Cl– absorption across the gills of Eriocheir sinensis[J]. J Exp Biol, 1995, 198(3): 767-774. doi: 10.1242/jeb.198.3.767
[38] 邢永泽, 乔莹, 米铁柱, 等. 海月水母螅状体横裂生殖过程中渗透压调节基因表达及酶活性变化分析[J]. 中国海洋大学学报(自然科学版), 2021, 51(S1): 19-23, 127. doi: 10.16441/j.cnki.hdxb.20200353 [39] 郭雯翡. 盐碱胁迫下青海湖裸鲤鳃抑制消减cDNA文库的构建及相关基因的克隆[D]. 上海: 上海海洋大学, 2012: 4-7. [40] GIFFARD-MENA I, BOULO V, AUJOULAT F, et al. Aquaporin molecular characterization in the sea-bass (Dicentrarchus labrax): the effect of salinity on AQP1 and AQP3 expression[J]. Comp Biochem Physoil A, 2007, 148(2): 430-444. doi: 10.1016/j.cbpa.2007.06.002
[41] 李蓓. 基于c-Myc/CLIC4通路探究脂多糖致大鼠海马损伤的致病机制[D]. 哈尔滨: 东北农业大学, 2019: 28-32.
计量
- 文章访问数: 491
- HTML全文浏览量: 185
- PDF下载量: 70
粤公网安备 44010502001741号
