豹纹鳃棘鲈差异流速下肝脏转录组分析

高进; 王永波; 刘金叶; 郭一兰; 符书源

doi:10.12131/20210125

豹纹鳃棘鲈差异流速下肝脏转录组分析

高进^{1, 2,},
王永波^{1, 2},
刘金叶^{2, 3},
郭一兰^{1, 2},
符书源^{1, 2, 3, ,}

1.
海南热带海洋学院/热带海洋生物资源利用与保护教育部重点实验室，海南三亚 572022
2.
海南省海洋与渔业科学院，海南海口 571126
3.
海南省热带海水养殖工程技术研究中心，海南海口 571126

基金项目: 海南省重点研发计划项目 (ZDYF2019074, ZDYF2020093)；2021海南省科研院所开发专项

详细信息

作者简介:
高　进 (1991—)，男，助理研究员，博士，从事水产动物基因组与数量遗传育种研究。E-mial: gaojin@hntou.edu.cn

通讯作者:
符书源 (1981—)，男，副研究员，硕士，从事海水鱼类繁育与健康养殖技术研究。E-mail: fulank23@sohu.com

中图分类号: S 917.4
计量
- 文章访问数: 651
- HTML全文浏览量: 240
- PDF下载量: 63
出版历程
- 收稿日期: 2021-04-21
- 修回日期: 2021-07-18
- 录用日期: 2021-07-28
- 网络出版日期: 2021-08-10
- 刊出日期: 2022-02-04

Transcriptome analysis of Plectropomus leopardus liver under different flow velocity

GAO Jin^{1, 2,},
WANG Yongbo^{1, 2},
LIU Jinye^{2, 3},
GUO Yilan^{1, 2},
FU Shuyuan^{1, 2, 3, ,}

1.
Hainan Tropical Ocean University/Key Laboratory of Utilization and Conservation for Tropical Marine Bioresources, Ministry of Education, Sanya 572022, China
2.
Hainan Academy of Ocean and Fisheries Sciences, Haikou 571126, China
3.
Hainan Provincial Engineering Research Center for Tropical Sea-Farming, Haikou 571126, China

摘要

摘要: 水流速度是影响鱼类生长的重要环境因子之一。为探究豹纹鳃棘鲈 (Plectropomus leopardus) 在不同流速条件下相关基因的功能和表达情况，利用RNA-Seq技术对差异流速下的豹纹鳃棘鲈肝脏组织进行了转录组分析。挑选相同繁育批次中规格一致的豹纹鳃棘鲈幼苗，分别在正常流速 (0.1 m·s⁻¹, Low flow velocity, LFV) 和高流速 (0.4 m·s⁻¹, High flow velocity, HFV) 实验组中养殖150 d后进行肝脏转录组测序分析，探究豹纹鳃棘鲈在差异流速下个体间的基因表达模式差异。结果显示，经筛选后共获得1 977个LFV-HFV显著性差异表达基因 (Differentially expressed genes, DGE)，其中上调基因999个，下调基因978个。GO功能注释分类发现，共有1 124个DEGs被GO数据库注释到并归属为56个功能类别。KEGG富集分析结果显示，573个DEGs参与了154条KEGG通路，其中富集最显著的为PPAR信号通路。LFV和HFV两组鱼的肝脏组织学观察结果显示，两者之间脂肪含量差异明显，LFV组肝脏中的脂肪含量显著高于HFV组 (P<0.05)。通过转录组分析，筛选了大量豹纹鳃棘鲈差异流速下的DGEs，为深入探讨豹纹鳃棘鲈对流速变化适应性的分子调控机制提供了技术支撑。
- 豹纹鳃棘鲈 /
- 差异流速 /
- 转录组分析
Abstract: Water flow velocity is one of important eco-environment factors which affects the fish growth. In order to explore the function and expression of related genes of Plectropomus leopardus under different flow velocity, we conducted a transcriptome analysis of liver tissue of P. leopardus under different flow velocity by RNA-seq technology. We selected the fish fry of P. leopardus with identical size from the same breeding batch and cultured them for 150 d with water flow velocity of 0.1 m·s⁻¹ (Low flow velocity, LFV) and 0.4 m·s⁻¹ (High flow velocity, HFV). Then, we conducted a transcriptome analysis on the liver so as to investigate the difference of gene expression patterns with different flow velocity. We had obtained a total of 1 977 differentially expressed genes (DEGs) by transcriptome analysis (999 up-regulated and 978 down-regulated for LFV-HFV, respectively). The GO functional annotation reveals that 1124 DEGs were annotated in Gene Ontology Consortium and assigned to 56 functional terms. KEGG pathway analysis shows that 573 DEGs belonged to 154 pathways, and PPAR signaling pathway was most significantly enriched. Histological observation of livers of the tested fish indicates that the difference in fat contents between LFV and HFV group was significant, and the fat content was obviously higher in LFV than in HFV (P<0.05). According to the transcriptome analysis, we excavated many DEGs under different flow velocity, which provids technical support for further research on the molecular regulation mechanism of adaptability to change in flow velocity of P. leopardus.
- Plectropomus leopardus /
- Different flow velocity /
- Transcriptome analysis

HTML全文

稻渔综合种养模式是在我国传统稻田养鱼的基础上，利用互利共生原则建立起来的多物种共栖、多层次利用的一种新型农业模式^[1-2]。这种养殖模式可以减少农药化肥的使用，提高水稻和鱼类品质，实现“一田两用”。金背鲤 (Cyprinus carpio var. Jinbei，暂定名) 因具有独特的“金背金额”外观特征而得名，是湘西地区在长期稻田养鱼过程中形成的鲤鱼地域性品种，其生长速度快于本地鲤，而且肉质厚实、口感鲜美，深受消费者的喜爱。

稻田养殖实现了多物种共栖，既可为养殖鱼类提供天然饵料，且其低养殖密度也可改善水质，促进鱼类生长，并提升鱼类的品质^[3-4]。有研究报道，鱼类肠道微生物是动态变化的，其群落结构受养殖模式和水体环境等因素的影响^[5]。良好的养殖模式会提升水产动物的营养与风味品质^[6]。如张艳霞等^[7]研究发现普通网箱养殖的大黄鱼 (Larimichthys crocea) 鱼腥味、哈喇味及青草味较重，而大围网箱养殖的则蘑菇味、脂香味、果香味较强。水产品的风味由气味和滋味两部分组成，它们可以通过复杂的气味协同、累加、掩盖、融合作用以及滋味的变调、对比、消杀、相乘作用相互影响，形成独特的风味^[8]。目前有关金背鲤的研究报道非常有限，探讨不同养殖模式下金背鲤肠道菌群和鱼肉风味品质的差异，对下一步金背鲤的繁殖优化具有重要意义。

本文采用高通量测序技术，对稻田和池塘养殖模式下金背鲤肠道微生物进行测序，同时利用液相色谱 (LC) 测定滋味物质、顶空固相微萃取-气相色谱-离子迁移谱联用技术 (HS-SPME-GC-IMS) 测定挥发性风味物质，结合感觉阈值与呈味、气味特征，计算滋味活性值和相对气味活度值，分析不同养殖模式下金背鲤的肠道菌群结构和菌群的差异及主要风味物质特征差异性，以为金背鲤下一步的选育繁殖提供研究方向，并为完善和推广稻田养殖模式提供参考数据。

1. 材料与方法

1.1 试剂

核苷酸标准物质 (纯度≥99%，美国药典 USP 标准品)；24种氨基酸标准品 (17种氨基酸为混合标准品 2.5 μmol·mL⁻¹，其他为分析纯，Sigma)；DNA凝胶回收试剂盒 (Qiagen公司) ；TruSeq文库构建试剂盒、MiSeq测序试剂盒 (美国Illumina公司)；溶菌酶 (美国BBI公司)。

1.2 主要仪器设备

SynergyHTX 多功能酶标仪 (美国 BioTe Tek公司)；Illumina HiSeq测序仪 (美国 Illumina 公司)；1200 高效液相色谱仪 (美国安捷伦公司)；1100 液相色谱仪 (美国安捷伦公司)；FlavourSpec^®风味分析仪 (山东海能科学仪器有限公司)。

1.3 方法

1.3.1 样品采集与处理

样品采集辰溪县凤凰山农业公司池塘养殖 (FGP) 组和稻田基地养殖 (FGF) 组，同一批次放养、体型大小相似的金背鲤，用冲氧袋运送至实验室，在FGP和FGF组随机各选5尾，在无菌状态下取出整个肠道，采集肠道内容物置于无菌的离心管中，液氮速冻，−80 ℃保存，用于提取总DNA。另外鱼放入塑料桶中静养 12 h后，各取5尾，剖杀后去鳞、内脏、皮等，手工取全鱼肌肉，搅碎机绞碎成均匀肉糜后分装，−40 ℃冻藏备用。

1.3.2 肠道微生物高通量测序分析

取肠道样本 (–80 ℃ 冰箱冻存) 放置在冰块上解冻，放入研钵后加入液氮研磨至粉末，使用预先冷冻处理的2 mL离心管，称量约100 mg样本进行总DNA提取。然后用乙醇沉淀纯化和1% (质量分数) 琼脂糖凝胶电泳检测所提取的DNA浓度和纯度，送百迈客生物公司使用Illumina Hiseq2500高通量测序平台进行测序分析。采用软件USEARCH version 10.0^[9]进行OTU分析，Mothurversionv.1.30^[10] 进行Alpha指数分析。通过RDP Classifier^[11] (置信度阈值为0.8) 进行物种注释，并分别在各个分类水平统计各样本的群落组成。使用LefSe进行组间差异显著性分析。

1.3.3 核苷酸分析

采用高效液相色谱仪分析核苷酸。样品采用10% (质量分数) 高氯酸溶液提取，提取液调节pH为 6.8，过0.22 µm滤膜上机测定，记录峰面积，响应值均应在仪器检测的线性范围之内，根据标准品的保留时间定性，外标法定量。

1.3.4 游离氨基酸分析

采用Agilent1100 液相色谱仪进行分析。称取0.50 g 样品于10 mL离心管中，加0.01 mol·L⁻¹ 盐酸 5 mL ，混匀，沸水浴30 min，离心取上清液，沉淀再加2 mL 0.01 mol·L⁻¹ 盐酸悬浮超声5 min后离心，合并上清液，定容至10 mL，过膜上机，通过安捷伦公司自动在线衍生化方法，比较游离氨基酸标准品的保留时间和峰面积来计算样品中游离氨基酸含量。

1.3.5 滋味活性值 (Taste activity value, TAV) 的计算

TAV 表示各呈味物质的浓度与该物质的阈值之比。参考徐永霞等^[12]的方法，计算不同养殖模式下金背鲤肌肉中各游离氨基酸和核苷酸的 TAV。

1.3.6 挥发性化合物分析

称取2.0 g 金背鲤肌肉样品置于一个20 mL顶空进样瓶中，在60 ℃下孵育15 min，随即用配置CTC多功能自动进样器和注射器的GC-IMS (Flavour Spec^®-G.A.S.，德国) 在不分流模式下进行后续分析。

测试条件：采用FS-SE-54-C4-CB-1 (15 m×0.53 mm) 色谱柱，柱温60 ℃，N₂为载气，流速150 mL·min⁻¹，IMS温度45 ℃；2 mL·min⁻¹维持至20 min，100 mL·min⁻¹维持至30 min，总运行30 min。采用LAV软件2.2.1 (G.A.S., 德国) 采集和处理GC-IMS数据，NIST和IMS数据库来对挥发物进行定性。

每个样品中加入10 μL 4-甲基-2戊醇 (20×10⁻⁶) 内标物计算各挥发性成分的相对含量。

1.3.7 相对气味活度值(Relative odor activityvalue, ROAV)法评价

采用ROAV法^[13]的计算公式如下：

$$ \mathrm{R}\mathrm{O}\mathrm{A}\mathrm{V}_i\approx \frac{{{{C}}_{\rm{ri}}}}{{{C}}_{\rm{rstan}}}\times \frac{T_{\rm{stan}}}{T_{{\rm{i}}}} \times 100{\text{%}} $$

(1)

式中：C_ri 为各挥发性组分的相对百分含量 (%，质量比)；C_rstan为对样品总体风味贡献最大组分的相对含量 (%，质量比)；T_stan 为对样品总体风味贡献最大组分相对应的感觉阈值 (μg·kg⁻¹)；T_i 为各挥发性组分相对应的感觉阈值 (μg·kg⁻¹)。

1.4 数据处理

通过 Excel 2010 和 SPSS 22.0 软件对数据进行统计分析，分析结果以“平均值±标准误 ($ \overline X $±SE)”表示。采用SPSS Statistics 23.0 软件的方差分析 (ANO-VA, P<0.05)；使用Origin 2018软件作微生物菌群箱线图和风味柱形图；微生物菌群与风味物质的相关性采用斯皮尔曼 (Spearman) 相关性系数计算。

2. 结果

2.1 不同养殖模式下金背鲤肠道菌群组成分析

2.1.1 OTU聚类及多样性分析

采用韦恩图对不同养殖模式相关特征进行分析，结果见图1 。两组样本共获得787个OTUs，稻田放养FGF组和池塘养殖FGP组共享的特征OTU有469个，占总数的59.6%。在特有的OTU中，FGF组196个(占比24.9%)，高于FGP组 (122个，占比15.5%)，主要分布于变形杆菌门、厚壁菌门和拟杆菌门。其Alpha多样性及差异如图2所示，FGF组的ACE指数、Chao1指数、香农-威纳指数 (Shannon-Wiener) 和辛普森指数 (Simpson) 均高于FGP组，但差异不显著 (P>0.05)。说明稻田放养组在物种丰富度和均匀度上优于池塘养殖组，但无显著性差异。

图 1 两种养殖模式的金背鲤肠道菌的特征韦恩图

Figure 1. Venn diagram of features of intestinal bacteria in C. carpio var. Jinbei of FGF and FGP groups

下载: 全尺寸图片幻灯片

图 2 两种养殖模式的金背鲤肠道微生物Alpha多样性指数

Figure 2. Alpha diversity of intestinal microorganisms in C. carpio var. Jinbei of FGF and FGP groups

下载: 全尺寸图片幻灯片

2.1.2 菌群结构组成差异性分析

对样本中相对丰度大于1%的OTU进行相对丰度和差异分析。如图3 所示，在生物学分类门水平上，两组养殖模式检出4个门水平，分别为变形菌门、梭杆菌门、厚壁菌门和拟杆菌门，相对丰度大于96%，这与黄颡鱼 (Pseudobagrus vachell)^[14]、草鱼 (Ctenopharyngodon idella)^[15]和罗非鱼 (Oreochromis mossambicus)^[16]肠道菌的结果一致。虽然养殖模式未改变优势菌门的种类，但影响了其相对丰度^[14]。根据图3 的单因素方差分析可知，除厚壁菌门无显著性差异 (P>0.05) 外，其变形菌门 (45.7% vs. 69.6%)、梭杆菌门 (42.8% vs. 19.5%)和拟杆菌门 (2.27% vs. 0.98%) 的相对丰度均存在显著性差异 (P<0.05)。

图 3 门水平上的物种相对丰度及差异分析

注：*. 差异显著 (P<0.05)。图4同此。

Figure 3. Relative abundance and differences of species at phylum level

Note: *. Significant differences (P<0.05). The same case in Fig. 4.

下载: 全尺寸图片幻灯片

属水平上的生物学分类如图4所示。优势菌属结构不变，其丰度值有差异，相对丰度大于1%的属有6个：嗜水气单胞菌 (Aeromonas)、鲸杆菌属 (Cetobacterium)、肠杆菌 (Enterobacteriaceae)、弧菌属 (Vibrio)、乳球菌属 (Lactococcus) 和克雷伯菌属 (Klebsiella)；根据图4的单因素方差分析，除乳球菌属外，其他5个属之间存在显著性差异 (P<0.05)。

图 4 属水平上的物种相对丰度及差异分析

Figure 4. Relative abundance and differences of species at genus level

下载: 全尺寸图片幻灯片

为进一步了解不同养殖模式下金背鲤肠道细菌群落结构的差异性，采用组间LEfSe分析，在不同组间寻找具有统计学差异的Biomarker，更清晰地筛选出了两组中有显著性差异的微生物 (LDA阈值为3)。由图5 可知，在FGF组主要为弧菌、拟杆菌 (Bacteroides)、交替单胞菌 (Alteromonadales) 、希瓦氏菌 (Shewanella)、嗜冷假单胞菌 (Pseudomonas psychrophila) 和Brevinema；FGP组主要为莫拉克斯氏菌 (Moraxella) 和克雷伯菌属。

图 5 两种养殖模式的金背鲤肠道微生物菌群的差异性分析

Figure 5. Difference analysis of microbial flora in intestine of C. carpio var. Jinbei of FGF and FGP groups

下载: 全尺寸图片幻灯片

2.2 不同养殖模式下金背鲤鱼肌肉中滋味组成分析

2.2.1 核苷酸及其关联产物分析

采用高效液相色谱分析两种养殖模式金背鲤肌肉中的核苷酸及其关联化合物，结果见表1。可以看出，两种养殖模式中ATP及其关联物的含量差异显著 (P<0.05)，FGF组中鲜味肌苷酸IMP含量 (1.676 g·kg⁻¹) 远远高于FGP组 (0.246 g·kg⁻¹)；而呈苦味的次黄嘌呤核苷HXR、次黄嘌呤HX却低于FGP组。其呈味物质含量高但贡献率不一定大，贡献率通常用TAV表示，其值等于样品中非挥发性滋味活性物质的浓度除以该化合物的阈值。当TAV>l时，表示该化合物对整体滋味存在贡献；当TAV<l时，表示该滋味物质对于整体滋味贡献不大^[16]。FGF组中肌苷酸的TAV高达 6.705 ，FGP组仅有 0.985；FGF组和FGP组的腺苷酸TVA值分别为0.349和0.291；此外，稻田放养FGF组的K值 (8.805%) 小于池塘养殖FGP组 (21.305%)。由此推测，稻田放养组的金背鲤肉质鲜味和鲜度优于池塘养殖组。

表 1 两种养殖模式的金背鲤肌肉中核苷酸及其关联产物含量及K值 (N=3)

Table 1. Nucleotides contents, their related products and K value in muscle of C. carpio var. Jinbei of FGF and FGP groups (N=3)

核苷酸种类 Nucleotide	阈值^[17] Threshold/ (g·kg⁻¹)	呈味特征 Taste characteristics	核苷酸含量 Nucleotide content/(g·kg⁻¹)		滋味活性值 TAV
核苷酸种类 Nucleotide	阈值^[17] Threshold/ (g·kg⁻¹)	呈味特征 Taste characteristics	稻田放养 FGF	池塘养殖 FGP	稻田放养 FGF	池塘养殖 FGP
腺苷三磷酸 ATP	—	—	0.051±0.0066^a	0.07686±0.0012^b
腺苷二磷酸 ADP	—	—	0.885±0.0271^b	1.935±0.0163^a
腺苷酸 AMP	0.50	鲜甜味	0.1747±0.0173^a	0.146±0.00026^b	0.349	0.291
肌苷酸 IMP	0.25	鲜味	1.676±0.1109^a	0.246±0.0055^b	6.705	0.985
次黄嘌呤核苷 HXR	—	苦味	0.0162±0.0021^b	0.0342±0.0017^a
次黄嘌呤 HX	—	苦味	0.0826±0.0062^b	0.193±0.0035^a
Σ呈味核苷酸 ΣFlavoring nucleotide			2.885±0.17^a	2.63±0.028^b
K值 K value/%			8.805	21.305
注：“—”表示未查到相关阈值和呈味特征；同行中不同字母间存在显著性差异 (P<0.05)；下同。 Note: "—" indicates that the relevant threshold and flavor characteristics are not found; values with different letters within the same line are signifi-cantly different (P<0.05). The same case in the following tables.

下载: 导出CSV

| 显示表格

2.2.2 游离氨基酸组成分析

从表2可以看出，两种模式共检出23种游离氨基酸，包括8种必需氨基酸、2种半必需氨基酸和13种非必需氨基酸。除羟脯氨酸、肌氨酸、天冬酰胺含量无显著性差异外 (P<0.05)，其余氨基酸均有显著性差异 (P>0.05)。在游离氨基酸、必需氨基酸、半必需氨基酸和非必需氨基酸总量上，两组养殖模式有显著性差异 (P<0.05)。在两组养殖模式中，苦味氨基酸总量占比较高，分别占游离氨基酸总量的50.19%和53.74%，两组间具有显著性差异 (P<0.05)；其次是甜味氨基酸，两组间差异不显著 (P>0.05)，占比分别为32.26%和28.46%；鲜味和芳香族氨基酸比例较低，但差异显著 (P<0.05)，尤其FGF组鲜味氨基酸高于FGP组，分别为2.24%和1.53%。

表 2 两种养殖模式下金背鲤肌肉中游离氨基酸与滋味活性值

Table 2. Amino acid and taste activity value in muscle of C. carpio var. Jinbei of FGF and FGP groups

氨基酸 Amino acid	阈值^[18-19] Threshold/ (g·kg⁻¹)	呈味特征^[20-22] Taste characteristics	游离氨基酸含量 Free amino acid content/(g·kg⁻¹)		滋味活性值 TAV
氨基酸 Amino acid	阈值^[18-19] Threshold/ (g·kg⁻¹)	呈味特征^[20-22] Taste characteristics	稻田放养 FGF	池塘养殖 FGP	稻田放养 FGF	池塘养殖 FGP
苏氨酸 Thr^★	2.60	甜 (+)	0.312±0.0042^b	0.362±0.0041^a	0.119	0.138
赖氨酸 Lys^★	0.50	甜/苦 (−)	0.254±0.0081^b	0.495±0.011^a	0.495	0.991
缬氨酸 Val^▲	0.40	甜/苦 (−)	0.297±0.0031^b	0.357±0.0023^a	0.737	0.888
蛋氨酸 Met^▲	0.30	苦/甜/硫 (−)	0.031±0.00002^b	0.0459±0.0002^a	0.103	0.153
色氨酸 Trp^◆	0.90	苦/芳香 (−)	0.0284±0.0004^a	0.0256±0.0002^b	0.035	0.031
苯丙氨酸 Phe^◆	0.90	苦/芳香 (−)	0.091±0.0004^b	0.119±0.0027^a	0.101	0.13
异亮氨酸 Ile^▲	0.90	苦 (−)	0.0501±0.0006^b	0.819±0.0003^a	0.055	0.091
亮氨酸 Leu^▲	1.90	苦 (−)	0.0982±0.0026^b	0.137±0.00007^a	0.051	0.072
Σ必需氨基酸 ΣEAA			1.162±0.019^b	1.623±0.021^a	1.69	2.49
组氨酸 His^▲	0.20	酸/苦 (−)	1.825±0.002^b	2.139±0.016^a	9.117	10.64
精氨酸 Arg^▲	0.50	甜/苦 (−)	0.0816±0.0018^b	0.181±0.001^a	0.161	0.36
Σ半必需氨基酸 ΣSEAA			1.907±0.0038^b	2.32±0.017^a	9.28	11.00
天冬氨酸 Asp^●	0.03	鲜/酸 (+)	0.0281±0.0005^a	0.0143±0.0003^b	0.923	0.483
谷氨酸 Glu^●	0.05	鲜/酸 (+)	0.114±0.0014^a	0.0949±0.004^b	2.268	1.954
丙氨酸 Ala^★	0.60	鲜/甜 (+)	0.315±0.015^a	0.251±0.014^b	0.542	0.401
丝氨酸 Ser^★	1.50	甜 (+)	0.131±0.0027^b	0.159±0.0001^a	0.086	0.106
甘氨酸 Gly^★	1.30	甜 (+)	1.088±0.017^a	0.905±0.0018^b	0.827	0.697
脯氨酸 Pro^★	3.00	甜/苦 (+)	0.124±0.0016^b	0.174±0.012^a	0.041	0.055
肌氨酸 Sar	—	—	0.125±0.012^a	0.122±0.001^a
酪氨酸 Tyr^◆	0.91	苦/芳香 (−)	0.0693±0.0015^b	0.0869±0.001^a	0.075	0.095
半胱氨酸 Cys^▲	—	苦/甜/硫 (−)	0.639±0.027^a	0.550±0.007^b
Σ非必需氨基酸 ΣNEAA			2.689±0.083^a	2.357±0.041^b	4.279	2.868
羟脯氨酸 Hyp	—	—	0.575±0.01^a	0.604±0.02^a
谷氨酰胺 Gln	—	—	0.08±0.0019^b	0.102±0.0009^a
天冬酰胺 Asn	—	—	0.0031±0.0001^a	0.0034±0.0004^a
瓜氨酸 Cit	—	—	0.0357±0.0006^b	0.0654±0.0001^a
游离氨基酸总量 Total free amino acid			6.361±0.12^b	7.137±0.098^a
^● 鲜味氨基酸总量/占比 Total delicious amino acid/Proportion			0.1425±0.0019^a/2.24%		0.109±0.0042^b/1.53%
^★ 甜味氨基酸总量/占比 Total Sweet amino acid/Proportion			2.052±0.013^a/32.26%		2.031±0.033^a/28.46%
^▲ 苦味氨基酸总量/占比 Total bitter amino acid/Proportion			3.193±0.01^b/50.19%		3.836±0.0071^a/53.74%
^◆ 芳香氨基酸总量/占比 Total aromatic amino acids/Proportion			0.189±0.0021^b/2.96%		0.232±0.0039^a/3.25%

下载: 导出CSV

| 显示表格

呈味氨基酸含量高时，是否其贡献率也高，通常还要看滋味活度值的大小，因呈味氨基酸的滋味活度值及其含量高低直接影响鱼肉的滋味^[18-21]，FGF组，谷氨酸TVA大于1，高于FGP组，对鲜味贡献大。组氨酸是苦味氨基酸，两组TAV均大于1，但有研究表明，呈味氨基酸对风味的贡献，还需看游离氨基酸与其他风味成分的相互作用^[23-24]。所以，结合呈味核苷酸的分析结果，在FGF组中高含量的IMP 能与谷氨酸和天冬氨酸产生协同效应，使产品更具鲜味^[25-26] 。

2.3 不同养殖模式下金背鲤鱼肌肉中挥发性化合物组成分析

本研究首次采用HS-GC-IMS技术，在两种养殖模式的金背鲤肌肉中检测到43个峰，40种化合物，其中定性出37种物质，包括10种醛、14种醇、6种酮、3种酸、2种酯、1种芳香烃、1种醚 (表3、图6) 。两种养殖模式的金背鲤肌肉中主要挥发性化合物为醛类 (9.54%~21.42%)、醇类(12.79%~19.19%)、酮类 (55.46%~58.01%) ，与Lu等^[27]检测的新鲜鲤肌肉中主要挥发性化合物的结果一致。在10种醛类化合物中，两种养殖模式的金背鲤除丁醛无显著性差异外 (P>0.05)，其余的醛类化合物均存在显著性差异 (P<0.05)；在14种醇类化合物中，1-丁醇、乙醇、芳樟醇、2-甲基-1-丁醇无显著性差异 (P>0.05)，其他的存在显著性差异 (P<0.05)；在6种酮类化合物中均存在显著性差异 (P<0.05)。而且从图4的指纹图谱可以明显看出，两组金背鲤肌肉中挥发性化合物的差异性，FGF组以壬醛、辛醛、庚醛、己醛、戊醛、E-2-己烯醛、苯甲醛、1-庚醇、1-戊烯-3-醇、芳樟醇、3-羟基-2-丁酮、3-辛酮、1-戊烯-3-酮、2-戊酮、甲酸丁酯、乙酸乙酯、β-罗勒烯、二甲基硫醚等化合物为主，FGP组以2-丙酮、2-甲基丙酸、3-甲基丁酸、丙酸、1-辛烯-3-醇、Z-3-己烯醇、1-己醇、1-戊醇、2-甲基丁醇、3-甲基丁醇、1-丁醇、异丁醇、1-丙醇、2-丙醇、乙醇、2-丁酮、丙酮等为主。

表 3 两组养殖模式的金背鲤肌肉中挥发性化合物相对含量及其ROVA值

Table 3. Volatile compounds and relative odor activity value in muscle of C. carpio var. Jinbei of FGF and FGP groups　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　　μg·kg⁻¹

化合物 Compound	阈值^[28-29] Threshold/ (μg·kg⁻¹)	气味特征^[28] Odour characteristics	稻田放养 FGF		池塘养殖 FGP
化合物 Compound	阈值^[28-29] Threshold/ (μg·kg⁻¹)	气味特征^[28] Odour characteristics	相对含量 Relative content	相对气味活度值 ROAV	相对含量 Relative content	相对气味活度值 ROAV
醛类 Aldehydes			21.42%		9.54%
壬醛 Nonanal	1	脂香、青草味	12.07±2.45^a	100	4.384±0.63^b	100
辛醛 Octanal	0.7	果香叶	6.32±1.1^a	74.80	3.633±0.49^b	43.01
庚醛 Heptanal	3	鱼腥、烤鱼、哈喇味	9.67±1.43^a	26.71	5.641±1.34^b	15.58
E-2-己烯醛 (E)-2-Hexenal	17	青香、脂香、果香味	3.42±1.48^a	1.67	2.54±0.26^b	1.24
己醛-M Hexanal-M	4.5	鱼腥、青草味	44.35±3.09^a	81.56	30.11±4.64^b	55.44
己醛-D Hexanal-D	4.5	鱼腥、青草味	31.66±3.85^a	58.29	11.60±4.72^b	21.36
戊醛 Pentanal	20	果香味	13.1±0.77^a	5.43	6.623±0.85^b	2.74
2-甲基丁醛2-Methylbutanal	1	果香味	4.08±0.29^b	33.8	6.80±1.2^a	56.34
3-甲基丁醛 3-Methylbutanal	1.1	巧克力味、腐臭味	2.71±0.14^b	20.41	4.84±1.11^a	36.45
丁醛 Butanal	1.3	花香、水果香味	5.07±0.72^a	32.31	5.66±0.53^a	36.07
苯甲醛 Benzaldehyde	350	苦杏仁味	2.6±0.44^a	0.062	1.63±0.42^b	0.038
醇类 Alcohols			12.79%		19.19%
1-庚醇 1-Heptanol	330	脂味、酒香	2.07±0.33^a	0.052	1.35±0.31^b	0.034
1-辛烯-3-醇 1-Octen-3-ol	10	蘑菇味、泥土味	3.64±0.54^b	3.16	5.80±2.38^a	4.81
1-己醇 1-Hexanol	5.6	青草味	3.38±0.54^b	5.00	6.22±2.01^a	9.20
正戊醇 1-Pentanol	150.2	面包酒香、果香	3.29±0.19^b	0.18	5.33±0.3^a	0.29
1-戊烯-3-醇 1-Penten-3-ol	358.1	肉香味、鱼腥味	10.68±1.09^a	0.25	9.03±2.05^b	0.021
1-丁醇 1-Butanol	5 000	温和的杂醇油气息	1.21±0.08^a	<0.01	1.04±0.07^a	<0.01
乙醇 Ethanol	100 000	酒味	29.21±6.39^a	<0.01	26.75±2.58^a	<0.01
异丙醇 2-Propanol	—	—	2.85±0.09^b		21.92±3.37^a
3-甲基-1-丁醇 3-Methyl-1-butanol	250	发酵味、油脂味、	1.26±0.12^b	0.04	3.34±0.69^a	0.11
正丙醇 1-Propanol	8 505	—	18.04±1.36^b	0.018	80.69±4.78^a	0.079
异丁醇 Isobutanol	7 000	—	1.26±0.09^b	<0.01	2.16±0.27^a	<0.01
Z-3-己烯-1-醇 (Z)-3-Hexen-1-ol	250	—	0.39±0.06^b	0.013	1.15±0.04^a	0.038
芳樟醇 Linalool	30 000	热带香辛料味	1.96±0.14^a	<0.01	1.46±0.36^a	<0.01
2-甲基-1-丁醇 2-Methyl-1-butanol	1 200	—	1.4±0.04^a	<0.01	1.65±0.25^a	0.01
酮类 Ketones			55.46%		58.01%
3-辛酮 3-Octanone	28	酮香、青香	2.35±0.15^a	0.70	1.76±0.31^b	0.52
3-羟基-2-丁酮 3-Hydroxy-2-butanone	800	脂香味	24.67±1.85^a	0.26	16.73±0.34^b	0.17
1-戊烯-3-酮 1-Penten-3-one	1	蘑菇味、烤焦味	4.54±1.32^a	37.61	1.72±0.16^b	14.25
2-戊酮 2-Pentanone	13 800	果味	4.7±0.2^a	<0.01	4.27±0.42^b	<0.01
2-丁酮-D 2-Butanone-D	50 000	醚香、果香、樟脑香	18.42±1.83^b	<0.01	72.65±3.42^a	0.012
2-丁酮-M 2-Butanone-M	50 000	醚香、果香、樟脑香	35.87±0.82^b	<0.01	40.29±1.69^a	<0.01
2-丙酮 2-Propanone	—	青草味	259.14±7.81^b		370.14±8.67^a
酯类 Esters			2.06%		0.94%
甲酸丁酯 Butyl formate	—		2.95±0.45^a		1.841±0.11^b
乙酸乙酯-M Ethyl Acetate-M	5	醚香、果香	7.65±0.9^a	12.68	4.70±0.37^b	7.79
乙酸乙酯-D Ethyl Acetate-D	5	醚香、果香	2.36±0.5^a	3.91	1.72±0.29^a	2.85
酸类 Acids			1.25%		3.36%
2-甲基丙酸 2-Methylpropanoic acid	—	—	2±0.12^a		2.14±0.15^a
3-甲基丁酸 3-Methylbutanoic acid	—	—	0.57±0.07^b		1.44±0.12^a
丙酸 Propanoic acid	—	—	5.29±0.75^b		25.86±0.47^a
烃类 Hydrocarbon			0.77%		0.19%
β-罗勒烯 beta-Ocimene	34	脂香、甜味	4.85±0.69^a	1.18	1.65±0.77^b	0.40
醚类 Ethers			5.47%		1.86%
二甲基硫醚 Dimethyl sulfide	300	果香、青香味	34.46±2.78^a	0.95	16.28±0.45^b	0.045
其他类 Others			0.78%		6.89%
未鉴定物质1 Unidentified substance 1	—	—	2.49±0.34^b		7.76±0.33^a
未鉴定物质2 Unidentified substance 2	—	—	2.04±0.04^b		50.80±4.06^a
未鉴定物质3 Unidentified substance 3	—	—	0.41±0.03^b		1.71±0.06^a

下载: 导出CSV

| 显示表格

图 6 两种养殖模式的金背鲤肌肉中挥发性化合物GC-IMS 指纹图谱

Figure 6. GC-IMS gallrey plot of volatile compounds in muscle of C. carpio var. Jinbei of FGF and FGP groups

下载: 全尺寸图片幻灯片

由于挥发性化合物对总体风味的贡献大小取决于挥发性组分在风味体系中的浓度和感觉阈值。根据挥发性成分的相对百分含量和感觉阈值数值，计算得出挥发性风味化合物的ROAV，ROAV≥1的化合物为样品的关键气味成分，0.1≤ROAV<1的化合物对整体风味具有修饰作用。如表3中所示，FGF组中关键气味化合物 (ROAV≥1) 有14种：壬醛、辛醛、庚醛、E-2-己烯醛、己醛、戊醛、2-甲基丁醛、3-甲基丁醛、丁醛、1-辛烯-3-醇、1-己醇、1-戊烯-3-酮、乙酸乙酯和β-罗勒烯；对整体风味具有修饰作用的化合物 (0.1≤ROAV<1) 有正戊醇、1-戊烯-3-醇、3-辛酮、3-羟基-2-丁酮和二甲基硫醚。FGP组中关键气味化合物 (ROAV≥1) 与FGF组比，缺β-罗勒烯，共13种；而对整体风味具有修饰作用的化合物 (0.1≤ROAV<1) 有正戊醇、3-甲基1-丁醇、3-辛酮、3-羟基-2-丁酮和β-罗勒烯。

2.4 不同养殖模式下金背鲤肠道细菌与肌肉中风味化合物的相关性分析

对主要菌群和风味物质进行相关性分析，结果如图7 所示，肌苷酸、腺苷酸和甜味氨基酸与弧菌属、拟杆菌属、鲸杆菌属呈正相关，与狭义梭状芽孢杆菌1 (Clostridium_sensu_stricto 1)、嗜水气单胞菌、摩根菌属 (Morganella) 和克雷伯菌属呈负相关。芳香氨基酸与弧菌属、拟杆菌属存在极显著负相关性 (0.001≤P<0.01)，与摩根菌属和克雷伯菌属存在极显著正相关性 (0.001≤P<0.01)。苦味和鲜味氨基酸与拟杆菌属呈显著负相关性 (P<0.05)，与克雷伯菌属呈显著正相关性 (P<0.05)。酯类化合物 (COO)、醛类化合物 (CHO)、芳烃类化合物 (AH)、二甲基硫醚 (DS) 与拟杆菌属和弧菌属呈极显著正相关性 (P<0.01)，与克雷伯菌属和摩根菌属呈显著负相关性 (P<0.05)。酮类化合物 (CO)、醇类化合物 (OH)、酸类化合物 (COOH) 与克雷伯菌属和摩根菌属呈显著正相关性 (P<0.05)，与拟杆菌属和弧菌属呈极显著负相关性 (P<0.01)。

图 7 肠道菌群结构与风味物质的相关性热图分析

Figure 7. Correlation heat map analysis of intestinal flora structure and flavor substances

下载: 全尺寸图片幻灯片

3. 讨论

鱼类肠道微生物起源于水体环境和饵料，其原始菌的建立则取决于水体及鱼卵表面的微生物种类^[30]，水体和饵料中的微生物通过鱼类的生活和摄食过程进入体内，经过适应后定殖下来，与宿主形成整体。多项研究发现鱼类等水生动物的遗传背景、饲养环境、饲料类型等因素显著影响了其肠道菌群结构^[31-32] 。肠道菌群对鱼类的营养循环具有强烈的促进作用，如变形菌门对多糖、蛋白质和其他有机物的降解和发酵极其重要^[33]。厚壁菌门中的乳酸菌可以降低肠道的pH，增强水生动物的抗菌能力^[34]。气单胞菌可以阻止病原菌在肠道和养殖环境中过多繁殖，从而改善水体环境并提高宿主免疫力 ^[35]。本研究对两种不同养殖模式下金背鲤的肠道菌群进行分析，发现两组样本虽然在物种丰富度和均匀度上无显著性差异，但其特征OTU值和菌落结构组成存在一定的差异，稻田放养组主要为弧菌、拟杆菌、交替单胞菌、希瓦氏菌和嗜冷假单胞菌，而池塘养殖组主要为莫拉克斯氏菌和克雷伯菌属。微生物与风味物质的相关性分析结果显示，弧菌、拟杆菌、克雷白氏杆菌和摩根菌属与不同的风味物质呈显著相关性 (0.01≤P<0.05 ＆ 0.001≤P<0.01 ＆ P<0.001)，从而引起鱼肉风味的差异，这可能是由于稻田环境给金背鲤提供了更加丰富多样的饵料，以及水质环境和微生物等综合作用的结果。

核苷酸和游离氨基酸是水产品中重要的滋味物质^[36]，鲜味肌苷酸IMP具有较强的鲜味，腺苷酸AMP具有鲜甜味^[37-38]，且两种核苷酸还能与谷氨酸和天冬氨酸产生协同效应，使产品更具鲜味^[23-24]，其典型的滋味组分类型差异和共有组分含量差异会导致水产品不同的特征滋味^[39]。本研究中两组养殖模式金背鲤的呈味核苷酸组成差异明显，稻田放养组呈味核苷酸IMP高达1.676 g·kg⁻¹，TAV为 6.705，远超池塘养殖组；游离氨基酸组成也存在显著性差异 (P<0.05)，稻田放养组的谷氨酸、甘氨酸、丙氨酸等关键氨基酸明显高于池塘养殖组。这可能是因为与单一的池塘养殖相比，稻田放养的鱼类活动频率增加，为鱼类提供了更多、更丰富的浮游生物和昆虫类天然饵料，增加了其高蛋白质营养补给，从而提升了鱼肉的鲜味。有大量研究数据验证了这一说法，如吕敏等^[40]报道稻田养殖黄颡鱼肌肉中的风味氨基酸含量显著高于池塘养殖组，段青源等^[41]报道野生大黄鱼的呈味氨基酸优于人工养殖的。

水产品中挥发性化合物组成是一个极其复杂的体系，由各种风味前体物质经过酶解、自氧化和其他许多复杂的化学反应而生成^[42]。通常认为醛类化合物主要来自于氨基酸的降解、多不饱和脂肪酸的氧化以及微生物的作用，其阈值较低，且具有气味加和作用^[43] 。两组养殖模式的醛类化合物的ROAV除苯甲醛小于0.1外，其余醛类的相对气味活度值均大于1，对整体风味具有重要贡献。如己醛、庚醛、壬醛、戊醛等化合物已被证明与鲜鱼的香味相关联^[44]，而且在新鲜鱼中，己醛通常产生一种原生味、鲜香的特征香味^[45]；E-2-己烯醛呈现出果香、青香和明显的脂香气^[42]；2-甲基丁醛、3-甲基丁醛在鳕 (Gadus) 中呈现特殊气味^[46]。醇类化合物是由脂肪酸的氧化分解或羰基化合物还原产生的^[47] ，ROAV≥1的1-辛烯-3-醇类似蘑菇气味，是水生动物肌肉中重要的风味物质^[48]。酮类化合物是由氨基酸降解、不饱和脂肪酸热氧化降解和微生物作用产生^[49]，有研究报道酮类对消减腥味有一定的贡献^[50]。如1-戊烯-3-酮呈蘑菇味、烤焦味，ROAV远远大于1，对整体风味贡献较大。酯类化合物中的乙酸乙酯具有果香味。因为水产品风味相当复杂，并且能够与许多其他物质产生重叠的风味效应^[51]。因此，人们闻到的鱼味有可能是几类化合物协同作用的结果。

据报道，鱼类肠道中某些微生物具有营养功能，可以生成维生素和挥发性脂肪酸，如Sugita等^[52]发现淡水鱼类鲤、金鱼 (Carassias auratas)、罗非鱼等肠道菌群中有产生维生素 B₁₂ 的菌株；李学梅^[53]发现鲤体内脂肪积累与肠道中的拟杆菌呈负相关，与厚壁菌呈正相关；Yang等^[54]研究发现草鱼肠道菌群显著富集与脂质代谢相关的基因 (甘油脂代谢、脂肪酸代谢、脂质生物合成)。部分鱼类肠道微生物的新陈代谢产物包括胞外酶，可以对蛋白质和脂肪等起促进作用，如Asfie等^[55]发现在牙鲆 (Paralichthys olivaceus) 仔稚鱼消化道中分离出的弧菌有生产蛋白酶的能力。

本研究结果显示，不同养殖模式下，弧菌、拟杆菌、克雷白氏杆菌和摩根菌属与呈味核苷酸、游离氨基酸和挥发性化合物呈极显著相关性 (P<0.01)，推断养殖模式会影响鱼的肠道微生物及其风味品质，且稻田放养组的金背鲤品质特征优于池塘养殖组。

图 1 不同流速下豹纹鳃棘鲈肝脏组织 (H.E染色)

a. HFV低倍镜观察；b. LFV低倍镜观察；c. HFV高倍镜观察；d. LFV高倍镜观察。

Figure 1. Liver tissues (H.E staining) of P. leopardus with different flow velocities

a. Observation of sections with low magnification for HFV; b. Observation of sections with low magnification for LFV; 　　　　　　c. Observation of sections with high magnification for HFV; d. Observation of sections with high magnification for LFV.

下载: 全尺寸图片幻灯片

图 2 不同流速下豹纹鳃棘鲈肝脏组织 (油红O染色)

a. HFV低倍镜观察；b. LFV低倍镜观察；c. HFV高倍镜观察；d. LFV高倍镜观察。

Figure 2. Liver tissues (Oil red O staining) of P. leopardus with different flow velocities

下载: 全尺寸图片幻灯片

图 3 豹纹鳃棘鲈基因GO分类图

Figure 3. Gene ontology (GO) assignment of genes of P. leopardus

下载: 全尺寸图片幻灯片

图 4 豹纹鳃棘鲈差异表达基因KEGG富集图

Figure 4. KEGG enrichment of differentially expressed genes of P. leopardus

下载: 全尺寸图片幻灯片

图 5 RT-PCR验证RNA-Seq结果

Figure 5. Validation of RNA-Seq data by using RT-qPCR

下载: 全尺寸图片幻灯片

图 6 SNP突变类型统计分布图

A. 腺嘌呤；T. 胸腺嘧啶；C. 胞嘧啶；G. 鸟嘌呤。

Figure 6. Statistical distribution of type of SNP mutations

A. Adenosine; T. Thymine; C. Cytosine; G. Guanine.

下载: 全尺寸图片幻灯片

表 1 本研究所用引物及其序列

Table 1 Sequences of primers used in this study

基因编号 Gene ID	基因名称 Gene name	引物序列 Primer sequence
utg000043l-1.624	cyp7b1	F: ACTTCATCGCCCTCTACCCTC
utg000043l-1.624	cyp7b1	R: TGAGCCTCTGACCGTCTTTG
utg000067l-2.111	nr1d2	F: CCTCTGGTTTCCATTACGGG
utg000067l-2.111	nr1d2	R: ATCAGGCATTTCTTGAAGCG
utg000067l-3.316	ciart	F: TTCAGTGAGAGCGAGCACAC
utg000067l-3.316	ciart	R: TTGGTTTCTTCAGGGCAGTG
utg000134l-0.116	cpt1a	F: AGCACCTGACTGACCGTAAGC
utg000134l-0.116	cpt1a	R: GCATCTCAAGTTCACTGGGTAAG
utg000150l-1.275	irs2	F: TGACATCAGCGACCCTTGTG
utg000150l-1.275	irs2	R: CGCCACTACTCTCTGTTGACG
utg000157l-1.277	scm2a	F: GCAGCAAAGACTGGAGCAAG
utg000157l-1.277	scm2a	R: TCGGTGAACTCATCTGGCAC
utg000241l-0.137	lat3	F: CAGGAAGAGATGTTGAACCTTGG
utg000241l-0.137	lat3	R: GCGAGGAAGATGAGACCAGAC
utg000253l-0.17	wsb2	F: GCCACAGGTTTGGAGAACAG
utg000253l-0.17	wsb2	R: AACACCAGGTCCCTCACTACAC
utg001131l-0.55	per1	F: CTCAAATACGCACTTCAATGTGTC
utg001131l-0.55	per1	R: TGAGGGTGTACTCTGAGGTGATG
utg000003l-2.243	acod	F: AGCAATGTTCTCCCTGAGGC
utg000003l-2.243	acod	R: CCAAAGCAAGGTCAAAGGATG
utg000129l-0.32	cyp2j2	F: GGCAACTTATTCTCTGTGGATTTC
utg000129l-0.32	cyp2j2	R: GCTGTCTCCCTGATTTACCAGTG
utg000134l-1.177	chka	F: CTTTCACATCACCATCATAAGGG
utg000134l-1.177	chka	R: CTTTGTTTGACTGTCGGGAATC
utg000537l-0.141	acsbg2	F: GCAGCAGAAGAGCCTGACCTAC
utg000537l-0.141	acsbg2	R: TAGATGCCAACAGCAAACCC
utg000714l-0.135	fasn	F: TGAGAAACACTCTGACACGAATG
utg000714l-0.135	fasn	R: TCAAAGCGTAGCCTCGGTAG

下载: 导出CSV

表 2 豹纹鳃棘鲈肝脏文库测序数据量统计

Table 2 Output statistics of P. leopardus liver library sequencing

样品编号 Sample ID	有效读长 Clean reads	GC含量 GC content/%	Q30的碱基数比例 Base ratio of Q30/%	比对效率 Mapping efficiency/%
SP10-1	26861 072	51.42	95.30	93.78
SP10-2	35829 110	51.51	95.05	93.53
SP10-3	23503 733	50.44	94.98	93.65
SP40-1	28594 219	51.54	95.23	93.16
SP40-2	28147 723	51.04	95.20	92.82
SP40-3	24705 424	51.00	94.80	92.77

下载: 导出CSV

表 3 豹纹鳃棘鲈不同流速显著差异通路中部分候选基因表达模式

Table 3 Expression pattern of partial candidate genes of P. leopardus in significant pathways with different flow velocity

基因ID Gene ID	log₂ (差异倍数) log₂ (Fold change)	FDR	登录号 Accession No.	基因名 Gene name
utg000134l-0.116	−5.33	1.87×10⁻³¹	XP_008289244.1	cpt1a
utg000043l-1.624	−3.55	9.54×10⁻¹⁶⁶	XP_008291411.1	cyp7b1
utg000150l-1.275	−2.30	1.64×10⁻¹³²	XP_018557299.1	irs2
utg000761l-0.82	−1.78	6.66×10⁻¹²⁵	XP_018522996.1	angptl4
utg000248l-0.27	−1.64	1.10×10⁻¹²	XP_018560228.1	socs3
utg000172l-0.359	−1.59	4.89×10⁻¹⁹⁵	XP_019209058.1	cyp46a1
utg000527l-0.63	−1.42	2.67×10⁻¹³	XP_018558522.1	acsl1
utg000455l_0.134	−1.41	9.36×10⁻³¹	XP_010731808.1	aqp7
utg000250l-2.242	−1.39	6.13×10⁻⁴⁶	XP_019119573.1	adipor2
utg000129l-0.162	−1.21	1.22×10⁻³⁷	XP_018522970.1	rxrgb
utg000097l-1.90	−1.08	1.01×10⁻⁷	XP_019120297.1	traf2
utg000349l-0.213	−1.03	1.05×10⁻⁵	XP_018534038.1	dgke
utg000794l-0.76	−1.02	1.54×10⁻⁶	XP_019126544.1	prkag2
utg000088l-2.212	−1.01	3.89×10⁻⁴⁷	XP_020505581.1	akr1a
utg000841l-0.8	1.01	4.13×10⁻⁸	XP_022072468.1	tradd
utg000329l-0.218	1.03	8.07×10⁻³²	XP_018527840.1	abhd16a
utg001296l-0.11	1.07	7.82×10⁻⁵	XP_018525976.1	slc27a4
utg000225l-0.155	1.11	6.08×10⁻³³	XP_018531738.1	lpin2
utg000681l-0.96	1.12	2.09×10⁻⁶	XP_018530981.1	acsl4
utg000325l-0.118	1.16	5.97×10⁻⁴	XP_005472382.1	hxk2
utg000841l_0.93	1.51	1.23×10⁻²⁴	XP_018522449.1	adpgk
utg000519l-0.406	1.82	1.04×10⁻²⁹	XP_018527726.1	kmt5c
utg000454l-0.43	1.98	3.21×10⁻¹²	XP_018528735.1	psat1
utg000537l-0.141	2.39	2.05×10⁻³³	XP_020511829.1	acsbg2
utg000393l-0.311	10.65	3.25×10⁻¹²⁸	XP_018529751.1	crhbp

下载: 导出CSV

表 4 SNP位点统计表

Table 4 Output statistics of SNP loci

样品编号 Sample ID	SNP位点数 Number of SNP	基因区SNP Genic SNP	基因间区SNP Intergenic SNP	转换类型占比 Transition/%	颠换类型占比 Transversion/%
SP10-1	161 388	121 218	40 170	30.88	93.78
SP10-2	178 270	133 923	44 347	31.22	93.53
SP10-3	171 685	129 389	42 296	31.37	93.65
SP40-1	168 801	126 800	42 001	31.69	93.16
SP40-2	171 772	128 821	42 951	31.79	92.82
SP40-3	165 186	123 814	41 372	31.79	92.77

下载: 导出CSV

参考文献(56)

[1]	FRISCH A, ANDERSON T. Physiological stress responses of two species of coral trout (Plectropomus leopardus and Plectropomus maculaus)[J]. Comp Biochem Physiol A, 2005, 140(3): 317-327. doi: 10.1016/j.cbpb.2005.01.014
[2]	关献涛, 吴洪喜, 马建忠, 等. 饵料中添加自然产物对豹纹鳃棘鲈生长和体色的影响[J]. 海洋学研究, 2018, 36(2): 80-91. doi: 10.3969/j.issn.1001-909X.2018.02.011
[3]	殷艳慧, 蒋万胜, 潘晓赋, 等. 水产养殖鱼类生长性状研究进展[J]. 中国水产科学, 2020, 27(4): 463-484.
[4]	王永波, 陈国华, 林彬, 等. 豹纹鳃棘鲈胚胎发育的初步观察[J]. 海洋科学, 2009, 33(3): 21-26.
[5]	BURGESS A I, CALLAN C K, TOUSE R, et al. Increasing survival and growth in larval leopard coral grouper (Plectropomus leopardus) using intensively cultured Parvocalanus crassirostris nauplii[J]. J World Aquacult Soc, 2020, 51(1): 171-182. doi: 10.1111/jwas.12635
[6]	陈超, 吴雷明, 李炎璐, 等. 豹纹鳃棘鲈 (Plectropomus leopardus) 早期形态与色素变化及添加剂对其体色的影响[J]. 渔业科学进展, 2014, 35(5): 83-90. doi: 10.11758/yykxjz.20140512
[7]	徐晓丽, 邵蓬, 李灏, 等. 豹纹鳃棘鲈致病性哈维氏弧菌的分离鉴定与系统发育分析[J]. 华中农业大学学报, 2014, 33(4): 112-118.
[8]	KENZO Y, KAZUHISA Y, KIMIO A, et al. Influence of light intensity on feeding, growth, and early survival of leopard local grouper (Plectropomus leopardus) larvae under mass-scale rearing conditions[J]. Aquaculture, 2008, 279(1/2/3/4): 55-62.
[9]	涂志刚, 蒋玉峰, 邱名毅. 三亚崖州区豹纹鳃棘鲈养殖现状与发展建议[J]. 中国水产, 2019(12): 54-55.
[10]	HOCKLEY F A, WILSON C A M E, BREW A, et al. Fish responses to flow velocity and turbulence in relation to size, sex and parasite load[J]. J R Soc Interface, 2014, 8(14): 1-11.
[11]	OGATA H Y, OKU H. Effects of water velocity on growth performance of juvenile Japanese flounder Paralichthys olivaceus[J]. J World Aquac Soc, 2000, 31(2): 225-231. doi: 10.1111/j.1749-7345.2000.tb00357.x
[12]	MERINO G E, PIEDRAHITA R H, CONKLIN D E. Effect of water velocity on the growth of California halibut (Paralichthys californicus) juveniles[J]. Aquaculture, 2007, 271(1/2/3/4): 206-215.
[13]	SUN G, LI M, WANG J, et al. Effects of flow rate on growth performance and welfare of juvenile turbot (Scophthalmus maximus L.) in recirculating aquaculture systems[J]. Aquacult Res, 2016, 47: 1341-1352. doi: 10.1111/are.12597
[14]	黄苑媚, 刘志刚, 谢恩义, 等. 水流速率对全缘马尾藻幼孢子体生长和生理活性的影响[J]. 广东海洋大学学报, 2014, 34(6): 45-50. doi: 10.3969/j.issn.1673-9159.2014.06.008
[15]	林明德, 陈刚, 马骞, 等. 杂交石斑鱼和母本褐点石斑鱼转录组测序及差异表达基因分析[J]. 广东海洋大学学报, 2019, 39(3): 15-23. doi: 10.3969/j.issn.1673-9159.2019.03.003
[16]	JULIANA C S, DOUGLAS D, MARTINS L F, et al. RNA-Seq differential expression analysis: an extended review and a software tool[J]. PLoS One, 2017, 12(12): e0190152. doi: 10.1371/journal.pone.0190152
[17]	WANG Z, GERSTEIN M, SNYDER M. RNA-Seq: a revolutionary tool for transcriptomics[J]. Nat Rev Genet, 2009, 10(1): 57-63. doi: 10.1038/nrg2484
[18]	QIAN X, BA Y, ZHUANG Q, et al. RNA-Seq technology and its application in fish transcriptomics[J]. Omics, 2014, 18(2): 98-110. doi: 10.1089/omi.2013.0110
[19]	罗辉, 叶华, 肖世俊, 等. 转录组学技术在水产动物研究中的运用[J]. 水产学报, 2015, 39(4): 598-607.
[20]	张旭, 周毅, 罗永巨, 等. 光周期影响罗非鱼脑组织转录组基因表达分析[J]. 西南农业学报, 2019, 32(11): 2704-2711.
[21]	赵超平. 卵形鲳鲹盐度适应调控机制研究[D]. 上海: 上海海洋大学, 2018: 9-34.
[22]	潘霞, 徐永健, 宁燕, 等. 温度胁迫对幼体大海马基因转录表达的影响[J]. 核农学报, 2020, 34(7): 1421-1431. doi: 10.11869/j.issn.100-8551.2020.07.1421
[23]	WANG L, YU C P, GUO L, et al. In silico comparative transcriptome analysis of two color morphs of the common coral trout (Plectropomus Leopardus)[J]. PLoS One, 2015(12): e0145868.
[24]	MEKUCHI M, SAKATA K, YAMAGUCHI T, et al. Trans-omics approaches used to characterise fish nutritional biorhythms in leopard coral grouper (Plectropomus leopardus)[J]. Sci Rep, 2017, 7(1): 9372. doi: 10.1038/s41598-017-09531-4
[25]	KIM D, LANGMEAD B, SALZBERG S L. HISAT: a fast spliced aligner with low memory requirements[J]. Nat Methods, 2015, 12(4): 357-360. doi: 10.1038/nmeth.3317
[26]	PERTEA M, PERTEA G M, ANTONESCU C M, et al. StringTie enables improved reconstruction of a transcriptome from RNA-seq reads[J]. Nat Biotechnol, 2015, 33(3): 290-295. doi: 10.1038/nbt.3122
[27]	ALTSCHUL S F, MADDEN T L, ZHANG J, et al. Gapped BLAST and PSI BLAST: a new generation of protein database search programs[J]. Nucleic Acids Res, 1997, 25(17): 3389-3402. doi: 10.1093/nar/25.17.3389
[28]	EDDY S R. Profile hidden Markov models[J]. Bioinformatics, 1998, 14(9): 755-763. doi: 10.1093/bioinformatics/14.9.755
[29]	YOUNG M D, WAKEFIELD M J, SMYTH G K, et al. Gene ontology analysis for RNA-seq: accounting for selection bias[J]. Genome Biol, 2014, 11(2): R14.
[30]	XIE C, MAO X, HUANG J, et al. KOBAS 2.0: a web server for annotation and identification of enriched pathways and diseases[J]. Nucleic Acids Res, 2011, 39(Suppl 2): W316-W322.
[31]	FLOREA L, SONG L, SALZBERG S L. Thousands of exon skipping events differentiate among splicing patterns in sixteen human tissues[J]. F1000 Res, 2013, 2: 188. doi: 10.12688/f1000research.2-188.v1
[32]	LOVE M I, HUBER W, ANDERS S. Moderated estimation of fold change and dispersion for RNA-seq data with DESeq2[J]. Genome Biol, 2014, 15(12): 550. doi: 10.1186/s13059-014-0550-8
[33]	LIVAK K J, SCHMITTGEN T D. Analysis of relative gene expression data using real-time quantitative PCR and the 2^{−∆∆ C_T} method[J]. Methods, 2001, 25(4): 402-408. doi: 10.1006/meth.2001.1262
[34]	许亚琴, 吴立新, 陈炜, 等. 水流对鱼类生理生态学影响的研究进展[J]. 现代农业科技, 2020, 4: 199-200. doi: 10.3969/j.issn.1007-5739.2020.07.122
[35]	娄宇栋, 冯建, 何娇娇, 等. 流速胁迫对美国红鱼的转录特性研究[J]. 浙江海洋大学学报(自然科学版), 2019, 38(1): 13-22.
[36]	许亚琴. 流速对拉氏鱥幼鱼生长、非特异性免疫能力及脂肪酸组成的影响[D]. 大连: 大连海洋大学, 2020: 8-19.
[37]	林浩然. 鱼类生长和生长激素分泌活动的调节[J]. 动物学报, 1996, 42(1): 69-79.
[38]	FAM B C, JOANNIDES C N, ANDRIKOPOULOS S. The liver: key in regulating appetite and body weight[J]. Adipocyte, 2012, 1(4): 259-264. doi: 10.4161/adip.21448
[39]	MUNOZ R, ESTANY J, TOR M, et al. Hepatic lipogenic enzyme expression in pigs is affected by selection for decreased backfat thickness at constant intramuscular fat content[J]. Meat Sci, 2013, 93(3): 746-751. doi: 10.1016/j.meatsci.2012.11.045
[40]	董忠典, 黎学友, 廖健, 等. 雌、雄弓背青鳉 (Oryzias curvinotus) 肝脏转录组比较分析[J]. 海洋与湖沼, 2020, 51(5): 1203-1213.
[41]	GAHR S A, VALLEJO R L, WEBER G M, et al. Effects of short-term growth hormone treatment on liver and muscle transcriptomes in rainbow trout (Oncorhynchus mykiss)[J]. Physiol Genomics, 2008, 32(3): 380-392. doi: 10.1152/physiolgenomics.00142.2007
[42]	刘凯, 谢楠, 冯晓宇, 等. 基于RNA-Seq技术对乌鳢和斑鳢肝脏的转录组分析[J]. 经济动物学报, 2015, 19(4): 213-219.
[43]	RISAU W, SARIOLA H, ZERWES H G, et al. Vasculogenesis and angiogenesis in embryonic-stem-cell-derived embryoid bodies[J]. Development, 1988, 102(3): 471-478. doi: 10.1242/dev.102.3.471
[44]	李建农, 蒋建东. 微管的生物学特性与药物研究[J]. 药学学报, 2003, 38(4): 311-315. doi: 10.3321/j.issn:0513-4870.2003.04.018
[45]	BECKWITH E J, YANOVSKY M J. Circadian regulation of gene expression: at the crossroads of transcriptional and post-transcriptional regulatory networks[J]. Curr Opin Genet Dev, 2014, 27: 35-42. doi: 10.1016/j.gde.2014.03.007
[46]	李优磊. PPAR信号通路在调控猪皮下脂肪与肌内脂肪差异沉积中的作用及机制研究[D]. 西安: 西北农林科技大学, 2018: 32-40.
[47]	杜滢. 氨基酸调控脂肪代谢的机制研究[D]. 北京: 中国科学院大学, 2013: 1-5.
[48]	吕怡航. 精氨酸对动物能量代谢的影响[J]. 饲料工业, 2014, 35(24): 40-46.
[49]	纪晨光. 胆汁酸对梗阻性黄疸肠黏膜屏障的保护作用及其机制的研究[D]. 石家庄: 河北医科大学, 2015: 1-14.
[50]	JIA W, WEI M L, RAJANI C, et al. Targeting the alternative bile acid synthetic pathway for metabolic diseases[J]. Protein Cell, 2021, 12(5): 411-425.
[51]	HAO F Q, TIAN M M, ZHANG X B, et al. Butyrate enhances CPT1A activity to promote fatty acid oxidation and iTreg differentiation[J]. PNAS, 2021, 118(22): e2014681118. doi: 10.1073/pnas.2014681118
[52]	梁计峻, 林亚秋, 俞雨阳, 等. 山羊CPT1A基因的克隆表达及肌内脂肪含量的相关性分析[J]. 华北农学报, 2019, 34(5): 231-238. doi: 10.7668/hbnxb.201751529
[53]	THOMAS Q, VALLIM D A, ELIZABETH J T. Pleiotropic roles of bile acids in metabolism[J]. Cell Metab, 2013, 17(5): 657-669. doi: 10.1016/j.cmet.2013.03.013
[54]	DONG X, PARK S, LIN X, et al. Irs1 and Irs2 signaling is essential for glucose homeostastis and systemic growth[J]. J Clin Invest, 2006, 116(1): 101-114. doi: 10.1172/JCI25735
[55]	ZHAO Q W, ZHANG Z, RONG W Q, et al. KMT5c modulates adipocyte thermogenesis by regulating Trp53 expression[J]. PNAS, 2020, 117(36): 22413-22422. doi: 10.1073/pnas.1922548117
[56]	D'ANDRE H C, PAUL W, SHEN X, et al. Identification and characterization of genes that control fat deposition in chickens[J]. J Anim Sci Biotechno, 2013, 4(1): 43. doi: 10.1186/2049-1891-4-43

施引文献(1)

期刊类型引用(1)

刘尊雷，杨林林，金艳，袁兴伟，张翼，张辉，许敏，程家骅，严利平. 基于CMSY和BSM的东海区重要渔业种类资源评估. 中国水产科学. 2023(06): 735-752 .

百度学术

其他类型引用(0)

资源附件(0)

图(6) / 表(4)

计量

文章访问数: 651
HTML全文浏览量: 240
PDF下载量: 63
被引次数: 1

1. 材料与方法
1.1 试剂
1.2 主要仪器设备
1.3 方法
1.3.1 样品采集与处理
1.3.2 肠道微生物高通量测序分析
1.3.3 核苷酸分析
1.3.4 游离氨基酸分析
1.3.5 滋味活性值 (Taste activity value, TAV) 的计算
1.3.6 挥发性化合物分析
1.3.7 相对气味活度值(Relative odor activityvalue, ROAV)法评价
1.4 数据处理
2. 结果
2.1 不同养殖模式下金背鲤肠道菌群组成分析
2.1.1 OTU聚类及多样性分析
2.1.2 菌群结构组成差异性分析
2.2 不同养殖模式下金背鲤鱼肌肉中滋味组成分析
2.2.1 核苷酸及其关联产物分析
2.2.2 游离氨基酸组成分析
2.3 不同养殖模式下金背鲤鱼肌肉中挥发性化合物组成分析
2.4 不同养殖模式下金背鲤肠道细菌与肌肉中风味化合物的相关性分析
3. 讨论

1. 材料与方法
1.1 试剂
1.2 主要仪器设备
1.3 方法
1.3.1 样品采集与处理
1.3.2 肠道微生物高通量测序分析
1.3.3 核苷酸分析
1.3.4 游离氨基酸分析
1.3.5 滋味活性值 (Taste activity value, TAV) 的计算
1.3.6 挥发性化合物分析
1.3.7 相对气味活度值(Relative odor activityvalue, ROAV)法评价
1.4 数据处理
2. 结果
2.1 不同养殖模式下金背鲤肠道菌群组成分析
2.1.1 OTU聚类及多样性分析
2.1.2 菌群结构组成差异性分析
2.2 不同养殖模式下金背鲤鱼肌肉中滋味组成分析
2.2.1 核苷酸及其关联产物分析
2.2.2 游离氨基酸组成分析
2.3 不同养殖模式下金背鲤鱼肌肉中挥发性化合物组成分析
2.4 不同养殖模式下金背鲤肠道细菌与肌肉中风味化合物的相关性分析
3. 讨论

参考文献(56)

施引文献(1)

资源附件(0)

豹纹鳃棘鲈差异流速下肝脏转录组分析

作者简介: 高 进 (1991—)，男，助理研究员，博士，从事水产动物基因组与数量遗传育种研究。E-mial: gaojin@hntou.edu.cn

通讯作者: 符书源 (1981—)，男，副研究员，硕士，从事海水鱼类繁育与健康养殖技术研究。E-mail: fulank23@sohu.com

计量

出版历程