基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析

伞利择; 刘宝锁; 张楠; 郭梁; 郭华阳; 朱克诚; 张殿昌

doi:10.12131/20210347

基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析

伞利择^{1, 2,},
刘宝锁^{2, 3},
张楠^{2, 3},
郭梁^{2, 3},
郭华阳^{2, 3},
朱克诚^{2, 3},
张殿昌^{2, 3, ,}

1.
河北农业大学海洋学院，河北秦皇岛 066003
2.
中国水产科学研究院南海水产研究所/农业农村部南海渔业资源开发利用重点实验室，广东广州 510300
3.
广东省海洋生物种业工程技术研究中心，广东广州 510300

基金项目: 国家自然科学基金项目 (U20A2064)；财政部和农业农村部: 国家海水鱼产业技术体系资助 (CARS-47)；国家重点研发计划项目 (2018YFD0900301)；中国水产科学研究院中央级公益性科研院所基本科研业务费专项资金资助 (2020TD29)

详细信息

作者简介:
伞利择 (1996—)，男，硕士研究生，研究方向为水产动物遗传育种。E-mail: cdsanlz@163.com

通讯作者:
张殿昌 (1977—)，男，研究员，博士，从事水产种质资源与遗传育种研究。E-mail: zhangdch@scsfri.ac.cn

中图分类号: S 917.4
计量
- 文章访问数: 740
- HTML全文浏览量: 200
- PDF下载量: 44
出版历程
- 收稿日期: 2021-11-23
- 修回日期: 2022-02-05
- 录用日期: 2022-02-17
- 网络出版日期: 2022-03-02
- 刊出日期: 2022-10-04

Mining of InDel marker and association analysis of hypoxia tolerance traits in Trachinotus ovatus based on resequencing

SAN Lize^{1, 2,},
LIU Baosuo^{2, 3},
ZHANG Nan^{2, 3},
GUO Liang^{2, 3},
GUO Huayang^{2, 3},
ZHU Kecheng^{2, 3},
ZHANG Dianchang^{2, 3, ,}

1.
Ocean College, Hebei Agricultural University, Qinhuangdao 066003, China
2.
South China Sea Fisheries Research Institute, Chinese Academy of Fishery Sciences/Key Laboratory of South China Sea Fishery Resources Exploitation & Utilization, Ministry of Agriculture and Rural Areas, Guangzhou 510300, China
3.
Guangdong Provincial Engineer Technology Research Center of Marine Biological Seed Industry, Guangzhou 510300, China

摘要

摘要: 卵形鲳鲹 (Trachinotus ovatus) 是高耗氧率的鱼类，低氧极易导致其死亡。筛选其耐低氧性状InDel分子标记并挖掘影响耐低氧性状的功能基因，可为选育具有较强耐低氧能力的卵形鲳鲹提供有益指导。利用全基因组重测序技术分析了卵形鲳鲹基因InDel差异，挖掘卵形鲳鲹与耐低氧性状显著关联的InDel位点，并探讨了与耐低氧性状相关的候选基因。结果表明，测序共获得693.48 Gb数据，其中Q30的平均值为90.8%，注释分析发现了共2 574 178个InDel位点。敏感组50尾卵形鲳鲹所特有的InDel共249 395个，其中2 209个位于外显子上，在核苷酸切除修复信号通路和细胞黏着分子中存在变异基因。将InDel位点与耐低氧性状关联分析发现，3个InDel位点 (InDel 22883061、InDel 24919481和InDel 14451779) 接近显著性阈值，注释分析获得了9个候选基因。筛选到的InDel位点对后期分子标记选择育种的选择和鉴定有重要价值，注释到的候选基因为卵形鲳鲹耐低氧的机理研究提供了基础和依据。
- 卵形鲳鲹 /
- InDel /
- 重测序 /
- 耐低氧性状
Abstract: Trachinotus ovatus is a fish with high oxygen consumption rate, and hypoxia can easily lead to its death. Screening for InDel molecular markers for hypoxia tolerance traits in T. ovatus and discovering functional genes affecting hypoxia tolerance traits can provide guidances for selecting and breeding T. ovatus with strong hypoxia tolerance. In this study, we applied the whole genome resequencing technology to analyze the InDel differences of T. ovatus gene, excavated the InDel sites that were significantly associated with hypoxia tolerance traits in T. ovatus, and explored the candidate genes related to hypoxia tolerance traits. Altogether 693.48 Gb was obtained by sequencing, and the average value of Q30 was 90.8%. A total of 2 574 178 InDel markers were found by annotation analysis. A total of 249 395 InDel markers were found in the 50 individuals in the tolerance group, of which 2 209 were located in exons. It is found that there were mutated genes in the nucleotide excise repair signaling pathway and cell adhesion molecules. Three InDel markers (InDel 22883061, InDel 24919481 and InDel 14451779) were found to be close to the significance threshold by association analysis between InDel markers and hypoxia tolerance trait. Nine candidate genes were obtained by annotation analysis. The selected InDel markers were of great value for selecting and identifying molecular marker selection breeding. The annotated candidate genes provide a basis and references for the study of the hypoxia tolerance mechanism of T. ovatus.
- Trachinotus ovatus /
- InDel /
- Resequencing /
- Hypoxia tolerance trait

HTML全文

乌鳢 (Channa argus) 是一种备受欢迎的淡水鱼，常作为原料生产调理鱼片或冷冻鱼片，但其加工过程中产生的鱼头、鱼骨等下脚料利用率较低，造成较大的资源浪费，并产生较严重的环境问题。中国素有熬制鱼头汤的饮食习惯，但较为耗时。乌鳢鱼头富含矿物质元素、氨基酸和不饱和脂肪酸，是熬制鱼汤的优选原料。因此开发方便即食的乌鳢鱼头汤制品具有较好的市场前景。目前市场上存在少量的鱼汤制品多为浓缩型汤膏，体积大且含水量高，保藏条件较严格，其消费受到较大的限制。随着干燥技术在速溶汤上的应用，出现了一些速溶汤制品。如采用微波冷冻干燥方法制备的速溶蔬菜汤^[1]；通过微波预糊化与红外冷冻干燥制备的奶油蘑菇汤粉^[2]；Singh等^[3]用喷雾干燥法制成的蘑菇乳清汤粉可在30 ℃下保存8个月。但鲜见作为粉末型的方便鱼汤产品。为适应现代人快节奏的生活方式和对美味鱼汤产品的需求，本研究选取喷雾干燥法制备乌鳢鱼头汤粉，探索干燥参数对干燥效果和产品品质的影响，优化出最佳干燥条件，为冲调式鱼头汤粉的生产提供技术支撑，从而提高乌鳢加工副产物的利用价值。

1. 材料与方法

1.1 材料与试剂

新鲜乌鳢 [(1.0±0.2) kg] 购于江苏省镇江市吉麦隆超市，宰杀后取鱼头，用碎冰保藏于泡沫箱中运至实验室；麦芽糊精 (食品级) 购于河南新乡中信化学有限公司。

1.2 仪器与设备

JSM-7800F场发射扫描电子显微镜 (日本电子株式会社)；均质机 (上海申鹿均质机有限公司)；OM-1500A小型喷雾干燥器 (上海欧蒙有限公司)；Lab-1C-50真空冷冻干燥机 (北京博医康实验仪器有限公司)；3810 (PAL-1) Digital Pocket Refractometer糖度计 (日本Atago公司)；MB27快速水分测定仪 (奥豪斯仪器 (常州) 有限公司)；CR-400色差计 (日本Konica Minolta, Inc.)。

1.3 乌鳢鱼头汤的制备

将新鲜乌鳢鱼头清洗后沥干，用刀剁成小块，称质量，用料酒和少量姜块腌制30 min，沥干备用。不锈钢锅内放入约50 mL的大豆油加热，放入鱼头块 [(250±5) g] 煎炸40~60 s至变色，然后再转入汤锅中[ m(料)∶V(液)=1∶6]，大火 (98 ℃) 熬煮30 min，再转小火 (90 ℃) 熬制，在熬煮至第60分钟时添加1%的食盐，总熬煮时间为150 min。当汤体温度降至40 ℃时，用勺子除去上层的油层，选取下层的鱼汤利用棉纱布进行过滤，得率约为最初添加水体积的40%，然后进行后续干燥。

1.4 喷雾干燥工艺参数筛选

1.4.1 麦芽糊精的添加量

物料黏度是喷雾干燥过程中的一个重要因素^[4]。选择常用的麦芽糊精作为助干剂控制黏度，添加质量浓度为0、5、10、15和20 g·L⁻¹。喷雾干燥进料速度设为10 mL·min⁻¹，进风温度为180 ℃，结束后收集样品并计算集粉率，根据粉末的水分质量分数、色泽、溶解性和堆积密度等指标及干燥过程中的粘壁情况确定麦芽糊精的最佳添加量。

1.4.2 进风温度对喷雾干燥效果的影响

在1.4.1确定麦芽糊精用量的基础上，将进风温度设定为60、170、180和190 ℃，进料速度固定为10 mL·min⁻¹，进行喷雾干燥实验。测定所得产品的集粉率、水分质量分数、色泽、溶解性和堆积密度等，筛选出最优的进风温度。

1.4.3 进料速度对喷雾干燥效果的影响

在1.4.1确定麦芽糊精用量的基础上，将进料速度设定为6、9、12和15 mL·min⁻¹，进口温度设为1.4.2得到的最优温度，进行喷雾干燥实验。测定产品集粉率，以及粉末的水分质量分数、色泽、溶解性和堆积密度等，确定最优的进料速度。

1.5 喷雾干燥鱼头汤粉特性评价

1.5.1 粉末颗粒形貌与粒径分析

利用扫描电镜观察所得粉末的形貌和表面特征。将粉末样品粘在贴有双面胶的圆柱形铝底座上，然后进行离子溅射喷金，在不同放大倍数下以15 kV的加速电压进行检测。选取适当放大倍数 (1 000×、5 000×和10 000×) 作为每个样品的代表。采用激光衍射粒度仪分析颗粒的粒度。选用通用的样品池，将2.0 g样品悬浮于20 mL蒸馏水中超声30 s，折光系数为1.50，利用光衍射法进行样品分析。平均粒径用体积平均直径D[3,4]表示，粒径分布用跨度 (Span) 表征。跨度指数计算为：

$$ I_{\rm {span }}=\left(D_{90}-D_{10}\right) / D_{50} $$

(1)

式中：I_span为跨度指数；D₉₀为90%颗粒的平均粒径 (μm)；D₅₀为50%颗粒的平均粒径 (μm)；D₁₀为10%颗粒的平均粒径 (μm)。

1.5.2 固形物质量分数测定

使用糖度计测定鱼头汤固形物的质量分数。

1.5.3 集粉率计算

参考郑唯等^[5]的方法，计算公式为：

$$ P=m_{0} /\left(m_{1}+m_{2}\right) \times 100 {\text{%}} $$

(2)

式中：P为集粉率 (%)；m₀为粉末质量 (g)；m₁为液料中的干物质质量 (g)；m₂为助干剂麦芽糊精添加质量 (g)。

1.5.4 堆积密度测定

参考Teo等^[6]的方法，即粉末的质量与在量筒中固定体积的比值。称量10 mL量筒中可堆积的样品质量，重复3次，取平均值，计算公式为：

$$ \rho=m / V $$

(3)

式中：ρ为堆积密度 (g·mL⁻¹)；m为粉末样品质量 (g)；V为粉末样品体积 (mL)。

1.5.5 色度测定

利用色差仪对所得粉末颜色进行测定分析，得到样品的明度 (L^*)、红绿色度 (a^*)、蓝黄色度 (b^*) 并计算白度。

1.5.6 水分质量分数测定

采用快速水分测定仪进行分析。

1.5.7 冲调时间

参考李秀军等^[7]的方法，用量筒量取 50~55 ℃的蒸馏水 8 mL 放入 10 mL试管中冲调1 g样品，并用玻璃棒搅拌，观察并记录样品完全溶解所需时间。

1.5.8 吸湿性

参考Cai和Corke^[8]的方法并适当修改。取1 g样品置于含饱和氯化钠 (NaCl) 溶液 (相对湿度75%) 的干燥器中，在25 ℃条件下放置1周后称质量，吸湿性为每100 g样吸附的水分质量 [g·(100 g)⁻¹]：

$$ A=\left(W_{1}+W_{2}\right)/W_1 \times 100 {\text{%}} $$

(4)

式中：A为吸湿性 (%)；W₁为样品初始质量 (g)；W₂为1周后的样品质量 (g)。

1.5.9 溶解度

取0.5 g样品加入50 mL蒸馏水中，搅拌器搅拌30 min，离心5~10 min (3 000~4 000 r·min⁻¹)，取上清液10 mL转移至玻璃皿中，烘干至恒质量：

$$ S=m / M \times 100 {\text{%}} $$

(5)

式中：S为溶解度 [g·(100 mL)⁻¹]；m为上清液中固形物的质量 (g)；M为样品溶液体积 (mL)。

1.6 数据分析

所有实验重复3次，结果取平均值。显著性分析用SPSS 20.0软件进行，以P<0.05 作为标准。采用Origin 8.0软件绘图。

2. 结果与分析

2.1 乌鳢鱼头汤喷雾干燥参数优化

2.1.1 助干剂对喷雾干燥产品的影响

助干剂 (麦芽糊精) 对鱼头汤喷雾干燥产品品质的影响见表1。未添加麦芽糊精时，鱼头汤经喷雾干燥后所得到的粉末积粉率和堆积密度均最低，同时冲调时间最长并存在严重的粘壁现象。而添加麦芽糊精可显著提高鱼头汤粉末的集粉率，当其质量浓度高于10 g·L⁻¹时，集粉率增至39.2%。但麦芽糊精质量浓度进一步增加对集粉率无显著影响。助干剂质量浓度较低时，料液中固形物的包埋效果差，影响玻璃化温度，导致雾化不充分^[9]，易出现粘壁现象，造成集粉率低。而随着添加量的增加，助干剂能在物料表面形成一层膜，从而改善粘壁现象，集粉率变大^[10-11]。此外，随着麦芽糊精质量浓度的增加，产品的水分质量分数逐渐下降且趋于稳定，当大于15 g·L⁻¹时水分质量分数降至5%以下，有利于长期存放。结果表明麦芽糊精的添加可改变物料间的粘黏程度，使颗粒间隔增大，利于水分的脱除，干燥加快。因此，麦芽糊精不仅可提高集粉率还有助于干燥过程中的脱水^[12-13]。

表 1 助干剂对产品喷雾干燥的影响

Table 1. Effect of carrier on spray drying of product

麦芽糊精质量浓度 Addition of maltdextrin/(g·L⁻¹)	集粉率 Powder collection rate/%	水分质量分数 Moisture mass fraction/%	堆积密度 Bulk density/ (g·mL⁻¹)	冲调时间 t/s	亮度 L^*	红绿色度 a^*	蓝黄色度 b^*	白度 Brightness	粘壁情况 Wall sticking
0	13.5^a	4.84^b	0.365^a	117^d	50.13^b	−0.14^d	7.01^c	49.64^b	粘壁严重
5	19.9^b	6.07^c	0.456^d	104^c	48.00^a	−0.31^a	3.39^b	47.89^a	粘壁严重
10	38.5^c	5.08^b	0.441^c	84^b	48.63^a	−0.26^b	2.55^a	48.57^a	部分粘壁
15	37.7^c	4.17^a	0.434^b	77^a	51.49^b	−0.35^a	2.39^a	51.43^b	轻微粘壁
20	39.2^c	4.03^a	0.429^b	74^a	50.60^b	−0.22^c	2.13^a	50.55^b	几乎不粘壁
注：同一列标注不同小写字母表示组间存在显著差异(P<0.05)；下表同此。 Note: Different lowercase letters within the same column indicate significant difference between groups (P<0.05). The same case in the following tables.

下载: 导出CSV

| 显示表格

此外，麦芽糊精的添加还会影响产品的堆积密度。未添加麦芽糊精时所得到的粉末堆积密度仅为0.365 g·mL⁻¹，当添加质量浓度为5 和15 g·L⁻¹时，堆积密度增大；而当添加质量浓度超过15 g·L⁻¹时，堆积密度则无明显差异，可能是由于粉末水分质量分数不同造成。麦芽糊精的添加还导致冲调时间的缩短，表明产品的溶解度随着麦芽糊精的增多而提高。麦芽糊精影响了产品的色泽，未添加麦芽糊精或添加量较低时，b^*较高，粉末偏黄；随着添加量的增加，L^*增大，b^*显著降低并趋于稳定，白度升高，但添加量过大会有浓郁的麦芽糊精味，影响产品原有风味。因此，综合各参数，确定麦芽糊精的添加质量浓度为15 g·L⁻¹。该条件下干燥粉末的水分质量分数低于5%，溶解性较好且色泽佳。

2.1.2 进风温度对喷雾干燥产品的影响

进风温度对喷雾干燥过程中能量传递及利用率会产生较大影响，是影响喷雾干燥效果和产品品质的重要参数。进风温度对鱼头汤喷雾干燥的影响见表2。随着温度的上升，粉末的集粉率呈现先增大后减小的趋势，在170 ℃时达到最大。温度较低时，液体无法迅速干燥，半干的颗粒会粘壁，导致集粉率下降；温度较高时，热量能够充分传递给喷出的雾滴，干燥室内的气流分子无规则运动加剧，水分蒸发速率加快，物料在气化后能迅速干燥。提高干燥速度，粘壁现象得到改善，粉末不易结块。但若温度过高，粉体表面会快速形成较硬的外壳，导致颗粒内部水分不易向外扩散，造成水分含量又有所升高。此外，在较高温度下蛋白质和麦芽糊精粉末则容易糊化，不仅影响产品最终的口感，还会影响其色泽，造成品质下降。黄绍天和王步江^[14]报道进风温度过高会导致喷雾干燥后期料液温度大于黏流温度，糖类等物质发生轻微焦化现象，产生轻微的热粘壁现象。

表 2 进风温度对产品喷雾干燥的影响

Table 2. Effect of inlet air temperature on spray drying of product

温度 Temperature/ ℃	集粉率 Powder collection rate/%	水分质量分数 Moisture mass fraction/%	堆积密度 Bulk density/ (g·mL⁻¹)	冲调时间 t/min	亮度 L^*	红绿色度 a^*	蓝黄色度 b^*	白度 Brightness	粘壁情况 Wall sticking
160	41.86^c	6.11^c	0.390^b	1.40^b	54.59^b	−0.38^b	4.29^a	54.39^b	部分粘壁
170	43.02^d	4.77^a	0.362^b	1.32^a	56.29^c	−0.48^a	4.66^b	56.04^b	轻微粘壁
180	37.70^b	6.09^c	0.401^c	1.29^a	51.49^a	−0.32^b	4.30^a	51.29^a	轻微粘壁
190	35.41^a	4.97^b	0.304^a	1.46^c	51.50^a	−0.47^a	5.34^c	51.20^a	少量粘壁

下载: 导出CSV

| 显示表格

粉末的水分质量分数的变化规律与集粉率基本一致，在170 ℃时最低 (4.77%)，产品得率最高 (表2)。此外，干燥粉末的堆积密度在180 ℃时达到最大，而粉末的冲调时间在170和180 ℃时较短且无显著性差异 (P>0.05)。但当进风温度过高时，粉末的冲调时间则会变长，可能是因为温度低时物料干燥不充分导致溶解性降低，而温度过高会形成坚硬的外壳，则需要更多的时间才能被溶解。此外，进风温度对产品色泽也产生了显著影响，当升高至190 ℃时，b^*明显增大，可能是由于温度过高使物料糊化从而导致粉末颜色偏黄。粘壁情况除了在160 ℃下不理想外，其他组均呈现较好的效果。综上，本实验选择最佳的进风温度为170 ℃。

2.1.3 进料速度对产品喷雾干燥的影响

进料速度也是影响喷雾干燥过程中产品品质的重要参数。当进料速度过低时，鱼头汤的喷雾干燥时间较长，不仅增加能耗，还会影响生产效率。进料速度对干燥产品的影响见表3。集粉率随进料速度的增加总体上呈下降趋势，在较低进料速度范围内集粉率较高，但当进料速度增加到一定程度时，集粉率则显著降低 (P<0.05)。当进料速度较低时，水分能够充分蒸发，物料可在雾化室得到较彻底的干燥，但当进料速度过快时，物料雾化不充分，水分蒸发不彻底，会导致干燥不完全，产生大量粘壁，出现结块，同时还易造成聚集堵塞影响喷雾。结果显示，当鱼头汤进料速度为6 mL·min⁻¹时，水分质量分数最低 (3.91%)，堆积密度与其他组无明显差异，并且在该条件下溶解性较好，冲调时间小于1 min，且轻微粘壁。综上，本实验选择进料速度为6 mL·min⁻¹。

表 3 进料速度对产品喷雾干燥的影响

Table 3. Effect of feeding speed on spray drying of product

进料速度 Processing velocity/ (mL·min⁻¹)	集粉率 Powder collection rate/%	水分质量分数 Moisture mass fraction/%	堆积密度 Bulk density/ (g·mL⁻¹)	冲调时间 t/s	亮度 L^*	红绿色度 a^*	蓝黄色度 b^*	白度 Brightness	粘壁情况 Wall sticking
6	44.55^c	3.91^a	0.301^b	58^b	52.75^a	−0.30^b	4.01^b	52.58^a	轻微粘壁
9	40.64^b	4.54^b	0.297^b	47^a	55.80^b	−0.39^a	3.84^a	55.63^b	部分粘壁
12	42.29^b	4.93^c	0.296^b	75^d	53.53^a	−0.36^ab	4.36^c	53.32^a	轻微粘壁
15	34.99^a	5.41^d	0.287^a	65^c	54.89^b	−0.43^a	3.87^ab	54.72^b	严重粘壁

下载: 导出CSV

| 显示表格

2.2 喷雾干燥鱼头粉品质评价

2.2.1 产品宏观形态

喷雾干燥获得的是较细的白色粉状产品 (图1)。蒋丽施等^[15]利用白乌鱼 (Opniocepnalus argus)熬煮鱼汤，通过喷雾干燥所得到的鱼汤粉的白度也呈现白色，且麦芽糊精添加量越多白度越高。

图 1 喷雾干燥鱼头汤粉外观

Figure 1. Product by spray drying

下载: 全尺寸图片幻灯片

2.2.2 产品含水量

含水量是检验干燥产品的一个重要指标，与产品存放条件和时间有关^[16]。喷雾干燥过程中，物料与热空气充分接触并进行热量交换，水分从颗粒内部快速扩散，物料在雾化后形成细小的颗粒，呈均匀多孔的球状，有利于水分的扩散。本实验显示鱼头汤经过喷雾干燥得到的粉末水分质量分数仅为4.03%，有利于产品的后期存放。蒋丽施等^[15]同样采用最佳进风温度为170 ℃对白乌鱼鱼汤进行喷雾干燥，所得鱼汤粉的水分质量分数为10.26%。这可能与不同原料来源的鱼汤中组分种类和含量不同有关。

2.2.3 产品堆积密度

堆积密度是一个反映粉体质构的参数，与微观结构有关^[17-18]。堆积密度大可以减少包装和运输成本^[19]。表4显示了喷雾干燥得到的产品堆积密度。喷雾干燥过程中，液滴会被雾化，在这个过程中会受重力、阻力和碰撞等影响产生形变，产生不同的颗粒结构，导致粉末堆积密度增大。进料温度、进料速度和鱼汤的固形物含量是影响鱼头汤粉的重要因素^[20]，利用本研究优化的喷雾干燥条件所得鱼头汤粉末的堆积密度为0.542 g·mL⁻¹。朱琳芳^[20]利用喷雾干燥制备的鳙 (Aristichthys nobilis) 汤粉的堆积密度约为0.25 g·mL⁻¹。堆积密度与喷雾过程中所形成的粉末粒径有关，粒径越小，堆积密度越大。

表 4 喷雾干燥的乌鳢鱼头汤粉品质指标

Table 4. Quality index of C. argus head soup powder made by spray drying

水分质量分数 Moisture mass fraction/%	堆积密度 Bulk density/(g·mL⁻¹)	溶解度 Solubility/[g·(100 mL)⁻¹]	吸湿性 Hygroscopicity/[g·(100 g)⁻¹]
4.03±0.15	0.54±0.09	99.13±1.98	24.47±1.65

下载: 导出CSV

| 显示表格

2.2.4 产品溶解度

本实验中喷雾干燥产品的溶解度较高 (表4)，可能是因为喷雾干燥过程中所形成的颗粒直径小，且较为均匀，空隙更多，与水分接触的通道多，导致粉末溶解度较高^[21]，这与付露莹等^[22]的研究结果一致。鳙鱼汤的喷雾干燥粉也具有相当高的溶解度，但随进料浓度和进料温度的升高而降低，随进料速度的增加而增大^[20]。

2.2.5 产品吸湿性

喷雾干燥制得的粉末吸湿性较低，主要与添加的助干剂有关。在较低浓度下，麦芽糊精的添加提高了产量，降低了粉体的吸湿性^[23-24]，因为其具有较少的亲水性基团，显示出较低的吸湿性。通过本研究优化的喷雾干燥条件获得鱼头汤粉的吸湿性约为24.47 g·(100 g)⁻¹ (表4)。Ferrari等^[25]在对黑莓 (Rubus fruticosus) 果肉进行喷雾干燥时加入5%~10%的麦芽糊精，产品的吸湿性随麦芽糊精添加量的升高而降低，但过高时会导致粉末的吸湿性增大^[26-27]，因此，适量的助干剂才有利于降低产品的吸湿性。

2.2.6 粉体粒径分布

通过喷雾干燥制备的鱼头汤粉粉末粒径及其分布见表5。D[4,3]所代表的平均粒径显示喷雾干燥样品的平均粒径较小。根据D₉₀、D₅₀和D₁₀计算得到跨度指数接近1，表明粉末粒径分布较窄^[28]，进一步说明喷雾干燥获得的粉体大小更均匀。在进行喷雾干燥时，雾化的过程会将液体先分散为小雾滴，受热后水分蒸发，形成细小的粉末状颗粒，得到的粉体粒径小且分布相对均匀。

表 5 粉末粒径和跨度指数

Table 5. Particle sizes and span index of powder

粉末粒径 Particle size/μm				跨度指数 Span
D[4,3]	D₉₀	D₅₀	D₁₀	跨度指数 Span
2.84±1.15	5.43±0.57	2.80±0.07	0.87±0.01	1.63±0.15

下载: 导出CSV

| 显示表格

2.2.7 颗粒形貌

粉体的微观形态及结构与干燥产品的流动性、速溶性、堆积密度等指标密切相关。喷雾干燥得到的鱼头汤粉末的微观形态扫描电镜图见图2。如前所述，在喷雾干燥过程中经过雾化器液体被均匀地分散成雾滴，使得制成的颗粒形状基本相同。放大1 000倍大多呈现球型颗粒状或扁球状，而放大5 000和10 000倍喷雾干燥的颗粒形状较圆、非常规则，表面光滑、大小均匀且多孔，分散性好，颗粒间的空隙大。该结果进一步解释了喷雾干燥的粉末含水量低，同时具较好的流动性、速溶性和堆积密度的原因。一些喷雾干燥的颗粒会融合在一起，这可能与吸湿性、物料浓度、添加的辅料黏度以及成膜能力有关^[29]。

图 2 扫描电镜结果

Figure 2. SEM of the spray drying samples

下载: 全尺寸图片幻灯片

3. 结论

麦芽糊精可作为乌鳢鱼头汤喷雾干燥的助干剂，能有效改善喷雾干燥过程中的粉末粘壁现象，当添加质量浓度为15 g·L⁻¹、喷雾干燥进风温度为170 ℃、进料速度为6 mL·min⁻¹时，乌鳢鱼头汤喷雾干燥效果最佳，制得的鱼头汤粉末集粉率最高，水分质量分数低于5%，堆积密度较大，色泽佳，速溶性较好。此外，产品颗粒呈球状，粒径小且分布较均匀，表面光滑，颗粒间空隙较大，吸湿性较低，利于存放。总之，喷雾干燥方法可应用于方便快捷的鱼头汤粉末产品生产，但规模化生产工艺还有待进一步研究。

图 1 卵形鲳鲹基因组 GC含量和测序质量分布图

Figure 1. GC content and sequence quality in T.ovatus

下载: 全尺寸图片幻灯片

图 2 InDel在基因组上的分布位置

Figure 2. Region of InDel on genome

下载: 全尺寸图片幻灯片

图 3 InDel在基因组上的分布密度

Figure 3. Distribution density of InDel on genome

下载: 全尺寸图片幻灯片

图 4 KEGG富集分析

Figure 4. KEGG enrichment scatter plot

下载: 全尺寸图片幻灯片

图 5 100尾卵形鲳鲹亲缘关系热图 (a)、主成分分析 (b)、耐低氧性状全基因组关联分析的QQ图 (c) 和曼哈顿图 (d)

Figure 5. Heat map of genetic relationship (a), principal component analysis (b), Q-Q plot (c) and Manhattan plot (d) of genome-wide association analysis for low oxygen tolerance of 100 individuals of T. ovatus

下载: 全尺寸图片幻灯片

表 1 鉴定到的不同长度InDel的数量

Table 1 Number of InDels of different sizes identified

InDel类型 InDel Type	InDel大小 InDel Size/bp	数量 Number	百分比 Percentage/%	InDel类型 InDel Type	InDel大小 InDel Size/bp	数量 Number	百分比 Percentage/%
插入 Insertion	1	496002	44.94	缺失 Deletion	1	532256	37.7
	2	189179	17.14		2	292212	20.7
	3	58625	5.31		3	107471	7.6
	4	69620	6.31		4	106766	7.6
	5	23055	2.09		5	43919	3.1
	6	43192	3.91		6	62130	4.4
	7	13832	1.25		7	27167	1.9
	8	28729	2.60		8	39844	2.8
	9	11719	1.06		9	20564	1.5
	10	20236	1.83		10	27837	2.0
	≥11	149421	13.54		≥11	152545	10.8
总计 Total　　　　1103610				总计 Total　　　　1412711

下载: 导出CSV

表 2 变异基因GO功能分类注释

Table 2 GO classification of mutated genes

类别 Classification	条目 Term	变异基因数量 Number of variant genes
细胞组分 Cellular component	染色体	7
	顶体泡	6
	线粒体外膜	4
	电压门控钾离子通道复合物	3
	肌间盘	3
分子功能 Molecular function	核酸结合	11
	GTP结合	8
	电压门控离子通道活性	5
	镁离子结合	5
	转录辅抑制活性	5
生物过程 Biological process	细胞对DNA损伤刺激的反应	7
	Wnt信号通路	6
	细胞器组织	6
	骨骼肌组织发育	5
	DNA整合	5

下载: 导出CSV

参考文献(50)

[1]	PAERL H W, OTTEN T G. Harmful cyanobacterial blooms: causes, consequences, and controls[J]. Microb Ecol, 2013, 65(4): 995-1010. doi: 10.1007/s00248-012-0159-y
[2]	YIN F, GONG H, KE Q, et al. Stress, antioxidant defence and mucosal immune responses of the large yellow croaker Pseudosciaena crocea challenged with Cryptocaryon irritans[J]. Fish Shellfish Immunol, 2015, 47(1): 344-351. doi: 10.1016/j.fsi.2015.09.013
[3]	MOREIRA M, SCHRAMA D, SOARES F, et al. Physiological responses of reared sea bream (Sparus aurata Linnaeus, 1758) to an Amyloodinium ocellatum outbreak[J]. J Fish Dis, 2017, 40(11): 1545-1560. doi: 10.1111/jfd.12623
[4]	ONUKWUFOR J O, WOOD C M. The osmorespiratory compromise in rainbow trout (Oncorhynchus mykiss): the effects of fish size, hypoxia, temperature and strenuous exercise on gill diffusive water fluxes and sodium net loss rates[J]. Comp Biochem Physiol A, 2018, 219: 10-8.
[5]	OBIRIKORANG K A, ACHEAMPONG J N, DUODU C P, et al. Growth, metabolism and respiration in Nile tilapia (Oreochromis niloticus) exposed to chronic or periodic hypoxia[J]. Comp Biochem Physiol A, 2020, 248: 110768. doi: 10.1016/j.cbpa.2020.110768
[6]	POULSEN S B, JENSEN L F, NIELSEN K S, et al. Behaviour of rainbow trout Oncorhynchus mykiss presented with a choice of normoxia and stepwise progressive hypoxia[J]. J Fish Biol, 2011, 79(4): 969-979. doi: 10.1111/j.1095-8649.2011.03069.x
[7]	WOOD A T, CLARK T D, ELLIOTT N G, et al. The effects of constant and cyclical hypoxia on the survival, growth and metabolic physiology of incubating Atlantic salmon (Salmo salar)[J]. Aquaculture, 2020, 527: 735449. doi: 10.1016/j.aquaculture.2020.735449
[8]	SOLLID J, de ANGELIS P, GUNDERSEN K, et al. Hypoxia induces adaptive and reversible gross morphological changes in crucian carp gills[J]. J Exp Biol, 2003, 206(20): 3667-3673. doi: 10.1242/jeb.00594
[9]	MITROVIC D, DYMOWSKA A, NILSSON G E, et al. Physiological consequences of gill remodeling in goldfish (Carassius auratus) during exposure to long-term hypoxia[J]. Am J Physiol Regul Integr Comp Physiol, 2009, 297(1): R224-234. doi: 10.1152/ajpregu.00189.2009
[10]	WU C B, LIU Z Y, LI F G, et al. Gill remodeling in response to hypoxia and temperature occurs in the hypoxia sensitive blunt snout bream (Megalobrama amblycephala)[J]. Aquaculture, 2017, 479: 479-486. doi: 10.1016/j.aquaculture.2017.06.020
[11]	MUNEER P M A, SIVANANDAN R, GOPALAKRISHNAN A, et al. Development and characterization of RAPD and microsatellite markers for genetic variation analysis in the critically endangered yellow catfish Horabagrus nigricollaris (Teleostei: Horabagridae)[J]. Biochem Genet, 2011, 49(1/2): 83-95.
[12]	WATANABE T, YOSHIDA M, NAKAJIMA M, et al. Linkage mapping of AFLP and microsatellite DNA markers with the body color-and sex-determining loci in the guppy (Poecilia reticulata)[J]. Zool Sci, 2005, 22(8): 883-889. doi: 10.2108/zsj.22.883
[13]	POOMPUANG S, NA-NAKORN U. A preliminary genetic map of walking catfish (Clarias macrocephalus)[J]. Aquaculture, 2004, 232(1/2/3/4): 195-203.
[14]	LEI D J, ZHAO G, XIE P, et al. Analysis of genetic diversity of Leuciscus leuciscus baicalensis using novel microsatellite markers with cross-species transferability[J]. Genet Mol Res, 2017, 16(2): 16029376.
[15]	TSAI H Y, HAMILTON A, TINCH A E, et al. Genome wide association and genomic prediction for growth traits in juvenile farmed Atlantic salmon using a high density SNP array[J]. BMC Genom, 2015, 16: 969. doi: 10.1186/s12864-015-2117-9
[16]	TSAI H Y, HAMILTON A, TINCH A E, et al. Genomic prediction of host resistance to sea lice in farmed Atlantic salmon populations[J]. Genet Sel Evol, 2016, 48(1): 47. doi: 10.1186/s12711-016-0226-9
[17]	JIN Y, ZHOU T, GENG X, et al. A genome-wide association study of heat stress-associated SNPs in catfish[J]. Anim Genet, 2017, 48(2): 233-236. doi: 10.1111/age.12482
[18]	CHEN F, LAI F L, LUO M J, et al. The genome-wide landscape of small insertion and deletion mutations in Monopterus albus[J]. J Genet Genom, 2019, 46(2): 75-86. doi: 10.1016/j.jgg.2019.02.002
[19]	VASEMAGI A, GROSS R, PALM D, et al. Discovery and application of insertion-deletion (INDEL) polymorphisms for QTL mapping of early life-history traits in Atlantic salmon[J]. BMC Genom, 2010, 11: 156. doi: 10.1186/1471-2164-11-156
[20]	BRITTEN R J, ROWEN L, WILLIAMS J, et al. Majority of divergence between closely related DNA samples is due to InDels[J]. Proc Nat Acad Sci USA, 2003, 100(8): 4661-4665. doi: 10.1073/pnas.0330964100
[21]	ROCKAH-SHMUEL L, TOTH-PETROCZY A, SELA A, et al. Correlated occurrence and bypass of frame-shifting insertion-deletions (InDels) to give functional proteins[J]. PLOS Genet, 2013, 9(10): e1003882. doi: 10.1371/journal.pgen.1003882
[22]	罗林. 奶牛乳脂性状候选基因SNP及InDel位点筛选及其与产奶性状关联分析[D]. 大庆: 黑龙江八一农垦大学, 2021: 4-8.
[23]	岳晓鹏. 基于甘蓝型油菜基因组重测序开发InDel标记[D]. 武汉: 华中农业大学, 2014: 6.
[24]	陈静, 何吉祥, 樊佳佳, 等. 草鱼MyoD基因SNP和InDel标记的筛选及其与生长性状的关联分析[J]. 江苏农业学报, 2018, 34(3): 612-616. doi: 10.3969/j.issn.1000-4440.2018.03.019
[25]	CHEN S, ZHOU Y, CHEN Y, et al. Fastp: an ultra-fast all-in-one FASTQ preprocessor[J]. Bioinformatics, 2018, 34(17): 884-890. doi: 10.1093/bioinformatics/bty560
[26]	MCKENNA A, HANNA M, BANKS E, et al. The genome analysis toolkit: a MapReduce framework for analyzing next-generation DNA sequencing data[J]. Genome Res, 2010, 20(9): 1297-1303. doi: 10.1101/gr.107524.110
[27]	CINGOLANI P, PATEL V M, COON M, et al. Using Drosophila melanogaster as a model for genotoxic chemical mutational studies with a new program, SnpSift[J]. Front Genet, 2012, 3: 35.
[28]	WANG K, LI M Y, HAKONARSON H. ANNOVAR: functional annotation of genetic variants from high-throughput sequencing data[J]. Nucleic Acids Res, 2010, 38(16): e164. doi: 10.1093/nar/gkq603
[29]	PURCELL S, NEALE B, TODD-BROWN K, et al. PLINK: a tool set for whole-genome association and population-based linkage analyses[J]. Am J Human Genet, 2007, 81(3): 559-575. doi: 10.1086/519795
[30]	YANG J, LEE S H, GODDARD M E, et al. GCTA: a tool for genome-wide complex trait analysis[J]. Am J Human Genet, 2011, 88(1): 76-82. doi: 10.1016/j.ajhg.2010.11.011
[31]	ZHOU X, STEPHENS M. Genome-wide efficient mixed-model analysis for association studies[J]. Nature Genet, 2012, 44(7): 821-824. doi: 10.1038/ng.2310
[32]	ZHONG X, WANG X, ZHOU T, et al. Genome-wide association study reveals multiple novel QTL associated with low oxygen tolerance in hybrid catfish[J]. Mar Biotechnol (NY), 2017, 19(4): 379-390. doi: 10.1007/s10126-017-9757-5
[33]	WANG X, LIU S, JIANG C, et al. Multiple across-strain and within-strain QTLs suggest highly complex genetic architecture for hypoxia tolerance in channel catfish[J]. Mol Genet Genom, 2017, 292(1): 63-76. doi: 10.1007/s00438-016-1256-2
[34]	杨洁, 赫佳, 王丹碧, 等. InDel标记的研究和应用进展[J]. 生物多样性, 2016, 24(2): 237-243. doi: 10.17520/biods.2015205
[35]	SAN L Z, LIU B S, LIU B, et al. Genome-wide association study reveals multiple novel SNPs and putative candidate genes associated with low oxygen tolerance in golden pompano Trachinotus ovatus (Linnaeus 1758)[J]. Aquaculture, 2021, 544: 737098. doi: 10.1016/j.aquaculture.2021.737098
[36]	王慧芳, 周光现, 孙永峰, 等. 基于全基因组重测序技术的狮头鹅InDel标记分析[J]. 畜牧兽医学报, 2021, 52(3): 662-675. doi: 10.11843/j.issn.0366-6964.2021.03.010
[37]	吴迷, 汪念, 沈超, 等. 基于重测序的陆地棉InDel标记开发与评价[J]. 作物学报, 2019, 45(2): 196-203.
[38]	常玉晓. 水稻T-DNA插入突变体侧翼序列的分离和水稻DNA复制蛋白RPA1a及RPA2-3基因的功能研究[D]. 武汉: 华中农业大学, 2010: 16-20.
[39]	黄芳. 斑马鱼DNA聚合酶Delta四亚基和相关因子的制备鉴定及赤点石斑鱼IκBα基因的克隆与功能分析[D]. 镇江: 江苏大学, 2016: 3-8.
[40]	ZEBDA N, DUBROVSKYI O, BIRUKOV K G. Focal adhesion kinase regulation of mechanotransduction and its impact on endothelial cell functions[J]. Microvasc Res, 2012, 83(1): 71-81. doi: 10.1016/j.mvr.2011.06.007
[41]	黄智康, 江世贵, 周发林, 等. 基于InDel标记的斑节对虾早期性别鉴定方法的建立[J]. 南方水产科学, 2020, 16(3): 113-118. doi: 10.12131/20190222
[42]	牛姣姣. 芝麻枯萎病抗性关联分析[D]. 南京: 南京农业大学, 2014: 24-43.
[43]	薛蕾. 山羊InDels筛选鉴定及其与生产性状关联分析[D]. 重庆: 西南大学, 2018: 30-33.
[44]	YANG Z L, ZOU L Q, SUN T T, et al. Genome-wide association study using whole-genome sequencing identifies a genomic region on chromosome 6 associated with comb traits in Nandan-Yao chicken[J]. Front Genet, 2021, 12: 682501. doi: 10.3389/fgene.2021.682501
[45]	ELLIS J M, WONG W, WOLFGANG M J. Acyl coenzyme a thioesterase 7 regulates neuronal fatty acid metabolism to prevent neurotoxicity[J]. Mol Cell Biol, 2013, 33(9): 1869-1882. doi: 10.1128/MCB.01548-12
[46]	MAEKAWA R, MUTO H, HATAYAMA M, et al. Dysregulation of erythropoiesis and altered erythroblastic NMDA receptor-mediated calcium influx in Lrfn2-deficient mice[J]. PLOS ONE, 2021, 16(1): e0245624. doi: 10.1371/journal.pone.0245624
[47]	FORSTER J R, LOCHNIT G, STOHR H. Proteomic analysis of the membrane palmitoylated protein-4 (MPP4)-associated protein complex in the retina[J]. Exp Eye Res, 2009, 88(1): 39-46. doi: 10.1016/j.exer.2008.09.016
[48]	RUIZ-HERNANDEZ A, ROMERO-NAVA R, HUANG F Y, et al. Altered function and expression of the orphan GPR135 at the cardiovascular level in diabetic Wistar rats[J]. J Recept Sig Transd, 2018, 38(5/6): 484-491.
[49]	GROS-LOUIS F, KRIZ J, KABASHI E, et al. Als2 mRNA splicing variants detected in KO mice rescue severe motor dysfunction phenotype in Als2 knock-down zebrafish[J]. Human Mol Genet, 2008, 17(17): 2691-2702. doi: 10.1093/hmg/ddn171
[50]	HUANG L J, SZYMANSKA K, JENSEN V L, et al. TMEM237 is mutated in individuals with a joubert syndrome related disorder and expands the role of the TMEM family at the ciliary transition zone[J]. Am J Human Genet, 2011, 89(6): 713-730. doi: 10.1016/j.ajhg.2011.11.005

施引文献(3)

期刊类型引用(3)

1.	焦晓磊，李清，樊威，苏建，罗煜，吴俊，邓语，张崟. 白乌鱼头骨和脊骨理化成分分析及营养价值比较. 农产品加工. 2025(01): 52-56 . 百度学术
2.	王益娟，周芳，焦文娟，赵甜甜，林以琳，张业辉，南海军，黄文. 基于非靶向代谢组学的药膳鱼汤代谢产物差异分析. 中国调味品. 2024(05): 148-155 . 百度学术
3.	陈胜军，路美明，相欢，薛勇，黄卉，李来好，胡晓. 乌鳢营养评价与加工保鲜技术研究进展. 肉类研究. 2023(02): 40-45 . 百度学术

其他类型引用(0)

资源附件(0)

图(5) / 表(2)

计量

文章访问数: 740
HTML全文浏览量: 200
PDF下载量: 44
被引次数: 3

1. 材料与方法
1.1 材料与试剂
1.2 仪器与设备
1.3 乌鳢鱼头汤的制备
1.4 喷雾干燥工艺参数筛选
1.4.1 麦芽糊精的添加量
1.4.2 进风温度对喷雾干燥效果的影响
1.4.3 进料速度对喷雾干燥效果的影响
1.5 喷雾干燥鱼头汤粉特性评价
1.5.1 粉末颗粒形貌与粒径分析
1.5.2 固形物质量分数测定
1.5.3 集粉率计算
1.5.4 堆积密度测定
1.5.5 色度测定
1.5.6 水分质量分数测定
1.5.7 冲调时间
1.5.8 吸湿性
1.5.9 溶解度
1.6 数据分析
2. 结果与分析
2.1 乌鳢鱼头汤喷雾干燥参数优化
2.1.1 助干剂对喷雾干燥产品的影响
2.1.2 进风温度对喷雾干燥产品的影响
2.1.3 进料速度对产品喷雾干燥的影响
2.2 喷雾干燥鱼头粉品质评价
2.2.1 产品宏观形态
2.2.2 产品含水量
2.2.3 产品堆积密度
2.2.4 产品溶解度
2.2.5 产品吸湿性
2.2.6 粉体粒径分布
2.2.7 颗粒形貌
3. 结论

1. 材料与方法
1.1 材料与试剂
1.2 仪器与设备
1.3 乌鳢鱼头汤的制备
1.4 喷雾干燥工艺参数筛选
1.4.1 麦芽糊精的添加量
1.4.2 进风温度对喷雾干燥效果的影响
1.4.3 进料速度对喷雾干燥效果的影响
1.5 喷雾干燥鱼头汤粉特性评价
1.5.1 粉末颗粒形貌与粒径分析
1.5.2 固形物质量分数测定
1.5.3 集粉率计算
1.5.4 堆积密度测定
1.5.5 色度测定
1.5.6 水分质量分数测定
1.5.7 冲调时间
1.5.8 吸湿性
1.5.9 溶解度
1.6 数据分析
2. 结果与分析
2.1 乌鳢鱼头汤喷雾干燥参数优化
2.1.1 助干剂对喷雾干燥产品的影响
2.1.2 进风温度对喷雾干燥产品的影响
2.1.3 进料速度对产品喷雾干燥的影响
2.2 喷雾干燥鱼头粉品质评价
2.2.1 产品宏观形态
2.2.2 产品含水量
2.2.3 产品堆积密度
2.2.4 产品溶解度
2.2.5 产品吸湿性
2.2.6 粉体粒径分布
2.2.7 颗粒形貌
3. 结论

参考文献(50)

施引文献(3)

资源附件(0)

基于卵形鲳鲹基因组重测序的InDel标记挖掘与耐低氧性状关联分析

作者简介: 伞利择 (1996—)，男，硕士研究生，研究方向为水产动物遗传育种。E-mail: cdsanlz@163.com

通讯作者: 张殿昌 (1977—)，男，研究员，博士，从事水产种质资源与遗传育种研究。E-mail: zhangdch@scsfri.ac.cn

计量

出版历程