Analysis of nutritional composition and quality of muscles in Hemibagrus guttatus, H. macropterus and H. pluriradiatus
-
摘要:
为探究斑鳠 (Hemibagrus guttatus)、大鳍鳠 (H. macropterus) 和越鳠 (H. pluriradiatus) 的肌肉营养成分和品质差异，采用生化分析方法对3种鱼背部肌肉中的一般营养成分含量以及氨基酸、脂肪酸和矿物元素的组成与含量进行测定和评价。结果显示，斑鳠的水分含量显著高于大鳍鳠和越鳠 (p<0.05)；越鳠和大鳍鳠的蛋白质含量显著高于斑鳠 (p<0.05)；粗脂肪和总糖含量从高到低依次为越鳠、大鳍鳠、斑鳠，差异显著 (p<0.05)；灰分含量则无显著性差异 (p>0.05)。3种鱼类均测得16种常见氨基酸，其中以谷氨酸含量最高，必需氨基酸 (Essential amino acids, EAA) 中以赖氨酸含量最高；总氨基酸含量 (Total amino acids, TAA) 和必需氨基酸指数 (Essential amino acid index, EAAI) 分别为斑鳠最低，越鳠最高。氨基酸评分 (Amino acid score, AAS) 与化学评分 (Chemical score, CS) 显示，3种鱼类的第一限制氨基酸均为甲硫氨酸+半胱氨酸。3种鱼类均测得24种脂肪酸，比较多不饱和脂肪酸含量，班鳠含量最高，且3种鱼类均以n-6系列的C18:2n6c (亚油酸) 为主。斑鳠的致动脉粥样硬化指数 (Atherogenic index, AI=0.47) 和血栓形成指数 (Thrombosis index, TI=1.21) 最低。3种鱼类肌肉中富含多种无机元素，常量元素钾 (K) 含量最高，呈“高钾低钠 (Na)”的特征；微量元素中铁 (Fe)、锌 (Zn)、硼 (B) 含量较高；硒 (Se) 仅在越鳠中检出，锰 (Mn) 仅在大鳍鳠中检出。研究表明，越鳠在蛋白质、粗脂肪、总糖、氨基酸总量及必需氨基酸指数方面均高于斑鳠和大鳍鳠。斑鳠在水分和多不饱和脂肪酸含量方面高于越鳠和大鳍鳠。
Abstract:To explore the differences in the muscle nutritional components and quality among Hemibagrus guttatus, H. macropterus, and H. pluriradiatus, we applied the biochemical analysis methods to determine and evaluate the contents of general nutritional components, as well as the composition and contents of amino acids, fatty acids, and mineral elements in the dorsal muscles of the three fish species. The results show that the water content of H. guttatus was significantly higher than that of H. macropterus and H. pluriradiatus (p<0.05); the protein contents of H. pluriradiatus and H. macropterus were significantly higher than that of H. guttatus (p<0.05). The contents of crude fat and total sugar decreased following a descending order of H. pluriradiatus>H. macropterus>H. guttatus, with significant differences (p<0.05), while there was no significant differences in the ash content (p>0.05). Sixteen common amino acids had been detected in the three fish species. Glutamic acid had the highest content among them, and lysine had the highest content among essential amino acids (EAA). The total amino acid content (TAA) and essential amino acid index (EAAI) were the lowest in H. guttatus but the highest in H. pluriradiatus. Amino acid score (AAS) and chemical score (CS) indicate that the first limiting amino acid of the three fish species was methionine+cysteine. Twenty-four fatty acids had been detected in the three fish species. In terms of the contents of polyunsaturated fatty acids, H. guttatus had the highest content, and C18:2n6c (linoleic acid) of the n-6 series was the main one in the three fish species. The atherogenic index (AI=0.47) and thrombogenic index (TI=1.21) of H. guttatus were the lowest. The muscles of the three fish species were rich in various inorganic elements. The content of the major element potassium (K) was the highest, showing the characteristic of "high potassium and low sodium". The contents of trace elements iron (Fe), zinc (Zn), and boron (B) were relatively high. Selenium (Se) was only detected in H. pluriradiatus, and manganese (Mn) was only detected in H. macropterus. The study indicates that the protein, crude fat, total sugar, total amino acid content, and essential amino acid index of H. pluriradiatus are higher than H. guttatus and H. macropterus. And the water content and polyunsaturated fatty acid content of H. guttatus are higher than H. pluriradiatus and H. macropterus.
-
Keywords:
- Hemibagrus guttatus /
- H. macropterus /
- H. pluriradiatus /
- Muscle /
- Nutritional component
-
牡蛎 (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 斑鳠、大鳍鳠、越鳠肌肉的常规营养成分 (鲜质量)
Table 1 Basic components in muscle of H. guttatus, H. macropterus and H. pluriradiatus (Wet mass)
g·100 g−1 物种
Species水分
Moisture灰分
Ash粗蛋白
Crude Protein粗脂肪
Crude Lipid总糖
Total sugar斑鳠 H. guttatus 80.75±0.35a 1.35±0.07a 16.90±0.00a 0.70±0.00a 0.26±0.01a 大鳍鳠 H. macropterus 76.60±0.14b 1.35±0.07a 19.10±0.42b 1.95±0.07b 0.28±0.02a 越鳠 H. pluriradiatus 76.50±0.00b 1.50±0.00a 19.50±0.00b 2.40±0.00c 0.32±0.01b 注：同一列数据上字母不同表示差异显著 (p<0.05)，字母相同表示差异不显著 (p>0.05)。 Note: Different letters within the same column represent significant differences (p<0.05), while the same letters represent insignificant differences (p>0.05). 
下载: 导出CSV
表 2 斑鳠、大鳍鳠、越鳠肌肉中的脂肪酸组成及含量
Table 2 Composition and content of fatty acids in muscle of H. guttatus, H. macropterus and H. pluriradiatus
% 脂肪酸
Fatty
acids斑鳠
H.
guttatus大鳍鳠
H.
macropterus越鳠
H.
pluriradiatus月桂酸 (C12:0) 0.02±0.00a 0.05±0.00b 0.06±0.00c 豆蔻酸 (C14:0) 1.22±0.03a 1.91±0.04b 2.13±0.00c 十五烷酸 (C15:0) 0.14±0.00a 0.20±0.00b 0.24±0.00c 棕榈酸 (C16:0) 21.85±0.21a 21.35±0.21a 23.00±0.14b 十七烷酸 (C17:0) 0.19±0.00a 0.23±0.00a 0.27±0.00b 硬脂酸 (C18:0) 9.45±0.16a 7.59±0.24a 7.56±0.01a 花生酸 (C20:0) 0.24±0.01a 0.32±0.00b 0.35±0.00c 山嵛酸 (C22:0) 0.3±0.01a 0.15±0.00a 0.16±0.00b 木焦油酸 (C24:0) 0.30±0.01a 0.07±0.00a 0.06±0.00b 饱和脂肪酸 ΣSFA 33.71±0.07a 31.87±0.02b 33.82±0.14a 棕榈油酸 (C16:1n-7) 1.28±0.01a 1.96±0.02b 1.91±0.00c 豆蔻一烯酸 (C14:1n5) 0.09±0.00a 0.06±0.00b 0.06±0.00c 油酸 (C18:1n-9c) 22.80±0.14a 31.20±0.28b 29.55±0.07c 二十碳一烯酸 (C20:1) 1.20±0.02a 1.94±0.02b 1.85±0.00c 芥酸 (C22:1n-9) 0.34±0.01a 0.29±0.00b 0.24±0.01c 二十四碳一烯酸
(C24:1n-9)0.43±0.01a 0.26±0.01b 0.22±0.00c 单不饱和脂肪酸ΣMUFA 26.13±0.11a 35.7±0.27b 33.84±0.07c 亚油酸 (C18:2n-6c) 18.35±0.21a 19.85±0.07b 21.15±0.07c 亚麻酸 (C18:3n-3) 1.50±0.04a 1.70±0.02b 1.85±0.00c 二十碳二烯酸 (C20:2) 1.75±0.06a 1.31±0.04b 1.19±0.01c 二十碳三烯酸 (C20:3n-3) 4.51±0.14a 0.89±0.03a 0.8±0.01b 二十碳四烯酸 (C20:4n-6) 4.19±0.06a 1.43±0.02b 1.26±0.01c 二十碳五烯酸EPA
(C20:5n-3)0.34±0.00a 0.64±0.01b 0.61±0.00c 二十二碳六烯酸
DHA (C22:6n-3)5.09±0.11a 5.41±0.17b 4.39±0.05b 二十二碳四烯酸 (C22:4n-6) 1.14±0.04a 0.15±0.00a 0.14±0.00b 二十二碳五烯酸 (C22:5n-3) 3.31±0.07a 1.01±0.01a 0.90±0.00b 多不饱和脂肪酸ΣPUFA 40.18±0.22a 32.37±0.34b 32.28±0.14b 不饱和脂肪酸ΣUFA 66.31±0.10a 68.07±0.08b 66.12±0.21a ΣPUFA n-3 14.75±0.28a 9.64±0.21b 8.55±0.06c ΣPUFA n-6 23.68±0.12a 21.42±0.10b 22.55±0.06c n-3/n-6 0.62 0.45 0.38 ΣUFA/ΣSFA 1.97 2.14 1.96 ΣPUFA/ΣSFA 1.19 1.02 0.95 高胆固醇血症指数 HI 23.07 23.26 25.13 致动脉粥样硬化指数 AI 0.47 0.43 0.46 致血栓形成指数TI 1.21 1.03 1.11 注：HI表示高胆固醇血症指数 (C14:0+C16:0)；AI表示致动脉粥样硬化指数; TI表示致血栓形成指数。不同的上标字母表示组间差异显著 (p<0.05)。 Note: HI. Hypercholesterolemia index (C14:0+C16:0); AI. Atherogenic index; TI. Thrombogenic index. Different superscript letters represent significant differences among the groups (p<0.05).
下载: 导出CSV
表 3 斑鳠、大鳍鳠、越鳠肌肉中氨基酸的组成与含量
Table 3 Composition and content of hydrolyzed amino acids in muscle of H. guttatus, H. macropterus and H. pluriradiatus
g·100 g−1 氨基酸
Amino acids斑鳠
H.
guttatus大鳍鳠
H.
macropterus越鳠
H.
pluriradiatus缬氨酸 Val* 0.77±0.01a 0.89±0.04b 0.92±0.01b 亮氨酸Leu* 1.33±0.03a 1.52±0.09b 1.59±0.01b 苏氨酸 Thr* 0.77±0.02a 0.88±0.03b 0.91±0.00b 赖氨酸 Lys* 1.61±0.04a 1.80±0.11b 1.89±0.02b 异亮氨酸Ile* 0.73±0.01a 0.84±0.05b 0.88±0.01b 苯丙氨酸 Phe* 0.67±0.02a 0.77±0.04b 0.81±0.01b 甲硫氨酸 Met* 0.37±0.03a 0.38±0.01a 0.55±0.01b 天冬氨酸 Asp*** 1.71±0.04a 1.98±0.08b 2.06±0.01b 丝氨酸 Ser 0.67±0.02a 0.80±0.00b 0.81±0.00b 谷氨酸Glu*** 2.75±0.06a 3.10±0.13b 3.19±0.01b 甘氨酸 Gly*** 0.75±0.06a 0.85±0.00b 0.93±0.04b 丙氨酸 Ala*** 0.97±0.04a 1.14±0.01b 1.16±0.01b 酪氨酸 Tyr** 0.51±0.03a 0.61±0.04b 0.66±0.02b 脯氨酸 Pro 0.55±0.03a 0.61±0.00b 0.62±0.05b 组氨酸 His** 0.25±0.01a 0.29±0.02b 0.30±0.00b 精氨酸 Arg** 0.98±0.03a 1.15±0.03b 1.16±0.01b 氨基酸总量TAA 15.39±0.41a 17.61±0.59b 18.44±0.00b 必需氨基酸EAA 6.25±0.16a 7.08±0.38b 7.55±0.07b 半必需氨基酸HEAA 1.74±0.01a 2.05±0.03b 2.12±0.01b 鲜味氨基酸DAA 6.18±0.19a 7.07±0.19b 7.34±0.04b 非必需氨基酸NEAA 9.14±0.21a 10.53±0.23b 10.89±0.07b EAA/TAA% 40.61 40.20 40.94 EAA/NEAA% 68.38 67.24 69.33 DAA/TAA% 40.16 40.15 39.80 注：*必需氨基酸；**半必需氨基酸 ；***鲜味氨基酸。不同的上标字母表示组间差异显著 (p<0.05)。 Note: *. Essential amino acids; **. Half essential amino acids; ***. Delicious amino acids. Different superscripts represent significant differences among the groups (p<0.05).
下载: 导出CSV
表 4 斑鳠、大鳍鳠、越鳠必需氨基酸组成的评价
Table 4 Evaluation of essential amino acid composition in H. guttatus, H. macropterus and H. pluriradiatus
氨基酸
Amino acidsFAO评分模式
FAO score model/
(mg·g−1)鸡蛋蛋白
Egg protein/
(mg·g−1)斑鳠
H. guttatus大鳍鳠
H. macropterus越鳠
H. pluriradiatusAAS CS AAS CS AAS CS 苏氨酸 Thr 250 292 1.14 0.97 1.16 0.99 1.17 1.00 缬氨酸 Val 310 410 0.92 0.69 0.94 0.71 0.95 0.72 甲硫氨酸+半胱氨酸Met+Cys 220 386 0.62 0.35 0.56 0.32 0.80 0.46 异亮氨酸 Ile 250 331 1.08 0.82 1.10 0.83 1.12 0.85 亮氨酸 Leu 440 534 1.12 0.92 1.13 0.93 1.16 0.95 苯丙氨酸+酪氨酸Phe+Tyr 380 565 1.14 0.77 1.19 0.80 1.24 0.83 赖氨酸 Lys 340 441 1.75 1.44 1.73 1.34 1.78 1.77 必需氨基酸指数EAAI 78.52 78.70 84.40
下载: 导出CSV
表 5 斑鳠、大鳍鳠、越鳠重要矿质元素组成及质量分数
Table 5 Composition and mass fraction of mineral elements in H. guttatus, H. macropterus and H. pluriradiatus
mg·kg−1 矿质元素
Mineral element班鳠
H.
guttatus大鳍鳠
H.
macropterus越鳠
H.
pluriradiatus常量元素
Macroelement钾 (K) 3 285 3 782 3 328 镁 (Mg) 336 295 345 钠 (Na) 444 279 284 钙 (Ca) <3 14 31 微量元素
Microelement铁 (Fe) 1.9 2.3 1.4 铜 (Cu) 2.4 0.21 0.15 锌 (Zn) 4.6 2.9 2.5 硒 (Se) <0.1 <0.1 0.90 锰 (Mn) <0.1 0.11 <0.1 硼 (B) 2.2 5.3 11 注：常量元素<3表示未检出，微量元素<0.1表示未检出。 Note: Macroelement<3 represents undetected, and microelement<0.1 represents undetected.
下载: 导出CSV
-
[1] 陈继位, 麦克, 王仕刚. 斑鳠人工繁殖技术[J]. 农业工程技术, 2023, 43(28): 111-112. [2] 李晓莉, 朱永久, 杨德国, 等. 大鳍鳠人工繁殖研究进展及前景展望[J]. 水产科技情报, 2023, 50(2): 128-132. [3] 杜孝田. 元江越鳠的生物学研究[D]. 武汉: 华中农业大学, 2010: 1-6 [4] 黄飞鹤, 陈琴, 罗永巨, 等. 斑鳠的含肉率及肌肉营养评价[J]. 福建水产, 2001(2): 8-14. [5] 吴小明, 周立斌, 彭飞影. 斑鳠的含肉率及肌肉营养成分分析[J]. 水产科学, 2010, 29(7): 425-428. doi: 10.3969/j.issn.1003-1111.2010.07.010 [6] 陈定福, 何学福, 周启贵. 长吻(鱼危)与大鰭鳠的含肉率及鱼肉营养成分的比较研究[J]. 淡水渔业, 1988(5): 21-23,13. [7] 向枭, 叶元土, 周兴华, 等. 大鳍鳠的消化能力与营养价值[J]. 水产学报, 2003(4): 371-376. [8] 吴雪萍, 叶红霞, 姚又菊, 等. 3个地区可口革囊星虫营养成分及重金属含量的分析与评价[J]. 热带海洋学报, 2024, 43(2): 92-107. [9] PELLETT P L, YOUNG V R. Nutritional evaluation of protein foods[M]. Tokyo: The united national university publishing company, 1980: 26-29.
[10] 王光亚. 新编《食物成分表》[J]. 营养学报, 1991(2): 165. [11] 杨育凯, 黄小林, 舒琥, 等. 不同生境下黄斑篮子鱼肌肉营养成分比较分析[J]. 南方水产科学, 2023, 19(1): 128-135. doi: 10.12131/20210345 [12] 王林娜, 田永胜, 李振通, 等. 棕点石斑鱼、杂交“金虎斑”及珍珠龙胆的肌肉营养成分分析及品质评价[J]. 水产学报, 2023, 47(9): 113-121. [13] LIN S H, ZHANG M, LUO Q, et al. Effects of flow velocity on the muscle nutrient content and flesh quality of the Pacific abalone (Haliotis discus hannai)[J]. Aquaculture, 2024, 582: 740492. doi: 10.1016/j.aquaculture.2023.740492
[14] 杨学坤, 韩雪, 王清峰, 等. 2种饲养模式条件下三穗鸭的肌肉营养检测分析与评价[J]. 贵州畜牧兽医, 2024, 48(2): 11-14. doi: 10.3969/j.issn.1007-1474.2024.02.005 [15] 王恒志, 张新党, 薛绍伟, 等. 雌雄丝尾鳠肌肉营养成分分析与评价[J]. 科学养鱼, 2019(5): 64-66. [16] 王鑫鑫, 杜辉, 张志华, 等. 大黄堡湿地主要经济鱼类营养成分分析与评价[J]. 饲料研究, 2024, 47(13): 99-104. [17] 贺伟平, 吴敬东, 孙建林, 等. 野生长胸异鮡肌肉营养成分分析与评价[J]. 水产科技情报, 2023, 50(5): 314-319. [18] 张志华, 张达娟, 李泽利, 等. 于桥水库4种经济鱼类肌肉营养成分分析与评价[J]. 饲料工业, 2024, 45(16): 92-98. [19] 喻亚丽, 鲁晓蓉, 董立学, 等. 低盐和淡水养殖模式下草鱼肌肉营养和风味品质评价[J]. 大连海洋大学学报, 2024, 39(3): 470-479. [20] 贾成霞, 杨慕, 胡庆杰, 等. 6种野生鱼肌肉营养成分分析与评价[J]. 现代食品科技, 2023, 39(6): 27-36. [21] 秦志清, 梁萍, 林建斌, 等. 闽江马口鱼肌肉营养成分分析与评价[J]. 渔业研究, 2022, 44(5): 508-515. [22] 李绍明, 梁志强, 伍远安, 等. 华鳈肌肉营养成分分析及评价[J]. 安徽农业科学, 2011(31): 19233-19235. doi: 10.3969/j.issn.0517-6611.2011.31.078 [23] 胡玉婷, 江河, 凌俊, 等. 野生与池养麦穗鱼营养成分比较[J]. 安徽农业科学, 2017, 45(4): 95-98. doi: 10.3969/j.issn.0517-6611.2017.04.032 [24] 钟文武, 周睿, 鲁祥兴, 等. 澜沧江云南吻孔鲃肌肉营养成分分析与评价[J]. 饲料研究, 2024, 47(10): 118-124. [25] 黄辉, 储忝江, 王宇希, 等. 钱塘江华鳈肌肉营养成分分析[J]. 渔业研究, 2024, 46(3): 271-278. [26] 中国康复医学会心血管病专业委员会, 中国营养学会临床营养分会, 中华预防医学会慢性病预防与控制分会, 等. 心血管疾病营养处方专家共识[J]. 中华内科杂志, 2014, 53(2): 151-158. doi: 10.3760/cma.j.issn.0578-1426.2014.02.021 [27] 唐迪, 邹烨, 仰榴青. GC-MS分析木瓜籽油中的脂肪酸组成[J]. 江苏农业科学, 2012, 40(10): 301-302. doi: 10.3969/j.issn.1002-1302.2012.10.112 [28] 姚远, 张燕萍, 章海鑫, 等. 鄱阳湖流域野生和养殖翘嘴鳜肌肉营养成分评价[J]. 饲料研究, 2024, 47(14): 113-117. [29] 文亚洲, 韩登武, 王喜军, 等. 高山美利奴羊肉品质分析及脂肪酸特征研究[J]. 饲料研究, 2024(20): 110-114. [30] 宋长虹, 唐生, 郝克非, 等. 中国居民日常食物中脂肪酸含量的分析[J]. 食品与机械, 2014, 30(5): 61-63. [31] 郑华梁, 郑琳. 鱼油补充剂改善心血管亚健康人群血脂的疗效与安全性的Meta分析[J]. 现代养生, 2024, 24(8): 568-574. doi: 10.3969/j.issn.1671-0223(s).2024.08.003 [32] 黄和, 任波, 孙小曼. ω-3多不饱和脂肪酸健康机制及应用研究进展[J]. 食品科学技术学报, 2024, 42(5): 1-12. doi: 10.12301/spxb202400679 [33] 刘玉鹏, 夏佳, 秦炯. Omega-3多不饱和脂肪酸对脑发育及脑功能的价值[J]. 中国实用儿科杂志, 2023(10): 741-745. [34] DAAN K, ALESSANDRO M, ADALBERTA A F, et al. Comparative ecologic relationships of saturated fat, sucrose, food groups, and a Mediterranean food pattern score to 50-year coronary heart disease mortality rates among 16 cohorts of the Seven Countries Study[J]. Eur J Clin Nutr, 2018, 72(8): 1103-1110. doi: 10.1038/s41430-018-0183-1
[35] ULBRICHT T L V, SOUTHGATE D A T. Coronary heart disease: seven dietary factors[J]. Lancet, 1991, 338(8773): 985-992. doi: 10.1016/0140-6736(91)91846-M
[36] 张圣杰, 周康奇, 潘贤辉, 等. 多棱角螺腹足和内脏团营养成分及风味物质分析[J]. 淡水渔业, 2024, 54(6): 87-95. doi: 10.3969/j.issn.1000-6907.2024.06.010 [37] 孙若鹏, 杨宾兰, 刘康, 等. 广西北部湾海域不同地区红树蚬营养成分分析与品质评价[J]. 广东海洋大学学报, 2024, 44(4): 63-70. [38] 李元跃, 李容伟, 陈融斌, 等. 红树林下人工笼养中华乌塘鳢 (Bostrychus sinensis) 肌肉营养成分的分析与评价[J]. 海洋与湖沼, 2022, 53(5): 1170-1179. doi: 10.11693/hyhz20211200330 [39] 汤雨. 甘氨酸对凡纳滨对虾生长、生化指标及体成分的影响[D]. 天津: 天津农学院, 2021: 35-36 [40] 白海锋, 杨公社, 袁永锋, 等. 黄河滩莲菜地套养河蟹肌肉营养成分分析[J]. 江西水产科技, 2019(3): 5-7. doi: 10.3969/j.issn.1006-3188.2019.03.003 [41] 吴聪聪, 姜群, 张晓君, 等. 中华鳖5个养殖群体肌肉营养成分分析与评价[J]. 广东海洋大学学报, 2024, 44(2): 54-61. doi: 10.3969/j.issn.1673-9159.2024.02.008 [42] 姚丽丽, 董玉珊, 韩蒙蒙, 等. 肉制品中滋味物质释放、感知与呈味机制的研究进展[J]. 食品与发酵工业, 2024, 1-12. [43] 史晋绒, 刘点, 王永明, 等. 美丽华沙鳅肌肉营养成分分析与评价[J]. 中国农学通报, 2023, 39(17): 124-130. doi: 10.11924/j.issn.1000-6850.casb2022-0502 [44] 寇红岩, 周萌, 黄燕华, 等. 钙、磷和铁矿物质元素对中华鳖生长和免疫的影响[J]. 饲料研究, 2020, 43(4): 134-137. [45] LU S Y, SIWAPORN P K, CHEN C S, et al. Effects of black soldier fly larvae (Hermetia illucens L. ) as feed supplements on muscle nutrient composition, meat quality, and antioxidant capacity in Qianbei goat[J]. Anim Biosci, 2024(37): 2167-2177.
[46] WANG Z Y, HUANG P Y, LI X Y, et al. Comparison of growth performance, nutritional composition, and muscle transcriptome between two cultured varieties of the Chinese mitten crab (Eriocheir sinensis )[J]. Fishes, 2024, 9(4): 132. doi: 10.3390/fishes9040132
[47] WAN L T, PENG Y Y, YU H H, et al. Comparing the muscle nutritional quality of eight common wild-caught economic shrimp species from the East China Sea[J]. J Aquat Food Prod Technol, 2022, 31(6): 549-564. doi: 10.1080/10498850.2022.2081062
[48] 赵亭亭, 张岩, 陈超, 等. 3种养殖石斑鱼的肌肉营养成分分析与品质评价[J]. 渔业科学进展, 2018, 39(6): 89-96. [49] 战楠, 焦杏春. 微量元素与健康[J]. 地球, 2024(1): 27-30. [50] 许安阳, 王安利. 锌铜铁对水产动物免疫机能的影响[J]. 水产科学, 2004(9): 36-39. doi: 10.3969/j.issn.1003-1111.2004.09.012
计量
- 文章访问数: 0
- HTML全文浏览量: 0
- PDF下载量: 0
粤公网安备 44010502001741号
