Ingestion, digestion and food selection of crimson snapper(Lutjanus erythopterus) larvae and juveniles
-
摘要: 为了探明仔、稚鱼阶段红鳍笛鲷(Lutjanus erythopterus)投喂轮虫的最适密度和颗粒饲料驯化时间，文章研究了不同轮虫投喂密度下，红鳍笛鲷仔鱼的生长、存活和食物选择，并采用生长、存活、RNA/DNA比率和消化道上皮细胞高度等参数作为评价指标，研究了颗粒饲料驯化时间 [13 dph (孵化后天数，W13)、16 dph (W16)、19 dph (W19)和22 dph (W22)] 对红鳍笛鲷仔、稚鱼的影响。结果表明，轮虫投喂密度显著影响红鳍笛鲷仔鱼的摄食、饵料选择、生长和存活，轮虫投喂密度为10~20 个·mL–1时，仔鱼的生长和成活率均无显著差异但显著高于1个·mL–1和30个·mL–1处理组。颗粒饲料驯化实验中，W19和W22处理组个体的生长和成活率显著高于其他两组，W13处理组个体的RNA/DNA比率最低。22 dph时，W13和W16处理组个体消化道上皮细胞高度明显低于其他两组。根据研究结果，轮虫和卤虫无节幼体混合喂养期间，轮虫密度对仔、稚鱼的食物选择性存在显著影响，建议在红鳍笛鲷初始投喂阶段使用10~20 个·mL–1的轮虫密度投喂仔鱼，红鳍笛鲷13 dph即可驯料，而最佳驯料期为16~22 dph。Abstract: To determine the optimal feeding density of rotifer and weaning time for Lutjanus erythopterus larvae and juveniles, we studied the impact of different rotifer feeding densities on the growth, survival and food selectivity for juvenile L. erythopterus. Then we investigated the effect of weaning time [13 dph (days post hatching, W13), 16 dph (W16), 19 dph (W19) and 22 dph (W22)] on the rearing performance of larval and juvenile L. erythopterus by using the evaluation indicators of growth, survival, RNA/DNA ratio and epithelial cell height of the digestive tract. The results show that the feeding density of rotifer affected the ingestion, food selection, growth and survival rate of L. erythopterus juveniles significantly. When the feeding density of rotifer was 10–20 ind·mL–1, there was no significant difference in the growth and survival rate of the juveniles, but significantly higher than treatment groups of 1 ind·mL–1 and 30 ind·mL–1. In the weaning experiment, the growth and survival rates of W19 and W22 groups were significantly higher than those of the other two groups, and the W13 treatment group had the lowest RNA/DNA ratio. At 22 dph, the height of the epithelial cells of digestive tract was significantly lower in the W13 and W16 groups than in the other two groups. It is indicated that the rotifer density had significant impact on the food selectivity of L. erythopterus during mixed feeding period of rotifers and Artemia nauplii. The rotifer density of 10–20 ind·mL–1 is recommended for the initial feeding period of L. erythopterus, and the weaning of L. erythopterus can be started on 13 dph, but the best timing was 16–22 dph.
-
长体圆鲹(Decapterus macrosoma)，又名长身圆鲹，隶属于鲈形目、鲹科、圆鲹属，主要分布于中国南海、印度尼西亚、澳洲和日本南部沿海等地[1]。长体圆鲹在中国南海分布较广，是南海灯光围网渔业的主要捕捞对象，具有较高的经济价值[2-5]。目前，国内外学者关于长体圆鲹的研究主要集中在生长繁殖[6]和资源评估[7-8]方面，与种群遗传和分子标记相关的研究报道较少。微卫星标记仅见翟云等[9]开发蓝圆鲹微卫星标记中获得5个跨物种通用标记可于用长体圆鲹，并无专门针对长体圆鲹开发的微卫星标记。种群遗传信息的匮乏，将大大影响对其资源的评估和长期有效的管理。
微卫星分子标记因是共显性标记，具有多态性高、变异性强、数据易统计等突出优点[10]，广泛应用于海洋生物遗传结构及遗传多样性分析[11-12]。但由于微卫星标记通用性较差，常常具有极强的种属特异性。鱼类微卫星标记开发中多以二核苷酸重复为主[13-15]，普遍认为它们具有较高的遗传变异[16]，但是也有部分学者研究认为三、四核苷酸重复位点较二核苷酸重复具有更高的筛选效率和多态性[17-19]。
本研究通过RAD-Seq高通量测序方法开发长体圆鲹二、三核苷酸微卫星分子标记，并对测试群体进行多样性分析，同时比较二、三核苷酸的筛选效率和多态性差异，旨在为长体圆鲹种群遗传结构及遗传多样性分析提供技术基础，并为该资源的评估和管理提供帮助。
1. 材料与方法
1.1 样品采集与基因组DNA提取
长体圆鲹样品采集于中国南海中沙群岛东部海域，共35尾。剪取部分肌肉样品加入无水乙醇保存。每个样品剪取少量肌肉组织，使用海洋动物组织基因组DNA提取试剂盒(天根，北京) 提取基因组DNA，0.8%的琼脂糖凝胶电泳检测DNA提取质量，−20 ℃保存备用。
1.2 高通量测序与引物合成
使用HiSeq2000高通量测序仪(Illumina，USA) 对长体圆鲹基因组DNA进行RAD-seq (测序服务由广州基迪奥生物科技有限公司提供)，经生物信息学搜索出微卫星位点[20]。使用Premier 5.0软件在重复单元侧翼序列上选择性设计出112条引物，主要参数为：G-C含量为40%~60%，引物长度为18~25 bp，退火温度为45~60 ℃，预期产物长度180~320 bp。送上海英潍捷基贸易有限公司合成引物。
1.3 引物筛选与分型检测
选取3个样本混合成的基因组DNA为模板，优化PCR反应条件，对引物进行首轮筛选，琼脂糖电泳检测是否能扩增出稳定且均一的目的片段。之后选取8尾个体的基因组DNA作为模板，使用三引物法[21]，利用M13荧光接头引物进行PCR扩增，扩增产物送华大基因公司经毛细管电泳进行等位基因分型，检测引物是否具有多态性。PCR反应体系为15 μL，其中包括10×PCR Buffer 1.5 μL，2.5 mmol·L–1 MgCL2 1.2 μL，2 mmol·L–1 dNTPs 2 μL，M13正向引物(10 μmol·L–1) 0.2 μL，M13反向引物(10 μmol·L–1) 0.6 μL，M13通用荧光引物(10 μmol·L–1) 0.5 μL，Taq酶(5 U·μL–1) 0.15 μL，DNA模版1 μL，加双蒸水至15 μL。PCR扩增程序为：94 ℃预变性5 min；94 ℃变性45 s，55~60 ℃退火45 s，72 ℃延伸45 s，35个循环；94 ℃变性45 s，53 ℃退火45 s，72 ℃延伸45 s，8个循环；72 ℃延伸30 min。
1.4 长体圆鲹群体遗传学评价
使用35尾长体圆鲹个体的基因组DNA为模板，对通过筛选的微卫星标记的种群遗传学特征进行评价。PCR反应体系和条件、等位基因分型方法如上。使用软件Genepop 4.0[22]对每个标记的种群遗传学特征值进行计算，包括等位基因数(Na)、表观杂合度(Ho)和期望杂合度(He)，进行“哈迪-温伯格”平衡(HWE)检验和连锁不平衡检测，并对P值进行Bonferroni校正。使用Cervus 3.0.7[23]软件计算多态信息含量(PIC)。
2. 结果
2.1 高通量测序结果与微卫星位点分析
RAD-seq高通量测序共获得长体圆鲹基因组原始数据2.06 G，GC含量为41.43%，Q30达93.05%。说明测序结果质量较好，可用于后续分析。搜索后共获得微卫星序列58 180条，一至六核苷酸重复微卫星位点70 508个，其中二核苷酸重复微卫星位点最多(37 646个)，占总数的53.39% (表1)，说明二核苷酸重复为主要的微卫星类型。二核苷酸重复微卫星位点共有4种重复类型，4种类型重复微卫星数量相差较大，AC/GT类有29 754个，占二核苷酸重复的68.4%；AG/CT类有6 487个，占17.2%；AT/TA有1 340个，占3.6%；GC/CG仅有65个，占0.17%。
表 1 长体圆鲹基因组中不同类型SSR统计Table 1. Different types of SSR statistics in D. macrosoma genome重复单元
repeat unit微卫星数量/个
number of microsatellite占比/%
ratio一核苷酸 mono-nucleotide 8 184 11.61 二核苷酸 di-nucleotide 37 646 53.39 三核苷酸 tri-nucleotide 13 960 19.80 四核苷酸 tetra-nucleotide 7 741 10.98 五核苷酸 penta-nucleotide 2 255 3.20 六核苷酸 hexa-nucleotide 722 1.02 合计 total 70 508 100.00 2.2 PCR引物设计和筛选
选取112条二、三核苷酸重复序列设计引物，其中二核苷酸重复为81对，三核苷酸重复为31对。经过筛选后，共有27对引物通过筛选(表2)，27对引物扩增的序列中18个位点为二核苷酸重复，重复次数为9~14次；8个位点为三核苷酸重复，重复次数为6~10次。二核苷酸重复位点检出效率为22.2%，三核苷酸重复位点检出效率为25.8%。
表 2 27对长体圆鲹微卫星引物信息Table 2. Information of 27 pairs of primers in D. macrosoma位点
locus引物序列 (5'−3')
primer sequence重复单元
repeat motif退火温度/℃
annealing temperature期望长度/bp
allele sizeDma03 F:CCACGCCTATTGAGTTACAGA (CA)9 60 186 R:GAGCCAGTGGATGAACAGAGT Dma07 F:GCCCCTGTGGGTGTGTGA (CA)9 60 225 R:GGGTGGTGGGTTCGGTTT Dma12 F:GAACCAGTGCCTACAATAGA (AC)9 60 243 R:CTGCTCACGGTAAGTCCA Dma15 F:ACAGGAAGGAACAGGACAG (TG)10 55 254 R:TATTGAAGTGAAAAAGCCG Dma22 F:CGCTGTTGAAATGAAGAAGA (GT)10 60 317 R:AGTGATGTCGCCTCATAAAT Dma23 F:AAACTGAGGGCGAGATAGAGG (AC)10 55 190 R:CCACAGGCTGAGTAAACCAAC Dma26 F:ATCCCATTCACCGACATAG (TG)10 58 258 R:CTGTGGTATCGTTCCCTGT Dma28 F:TGATTGGCTTCTACTCTGC (AC)10 55 281 R:AGTGGCTTGTTTGACTCTTAT Dma36 F:GGATGTAGTGAAGAGGGGAG (GT)11 55 239 R:CACAATCAGTGTTATGGCAG Dma38 F:GCCAATAAAGGCAAACAGT (CA)11 60 227 R:ATCCGAGACAAAGACATACAA Dma39 F:AGTGTGCTGACTTTTCTCTG (CA)11 55 241 R:TTATTGTTTGTTGTCTGGGT Dma45 F:CTCCTTTTTCTTCTTCCTCT (CA)11 60 281 R:CTACCTGCTCTTCAACTCAT Dma51 F:TGACAGCCTCCACTACTCC (GA)12 55 225 R:GCTAACCAGACACGCAAA Dma54 F:AAAGCCCATCTGTCTCGT (GT)12 60 202 R:TGTTTCAGTCCGTTCCTG Dma58 F:TCAAGAGGGAGTGGGAGC (AC)12 58 279 R:TCAAATGGGTGTTTAGCG Dma64 F:GCTCAGACTGCGTGGACA (TG)13 55 314 R:GCTGGTGAACAACAGGACA Dma72 F:TTCCGCAGGCATAAAAAC (CT)13 58 301 R:CCAAGGTCCGCTACACTA Dma76 F:TTCTCGCTGACCTGCTTG (TG)14 55 253 R:GCGTCCTCGTCGTCTTTC Dma81 F:GAGACACGGTCAGAAAACA (TGC)6 60 216 R:GGAAGTAGGACTCTAGGGG Dma82 F:CTGTCACTCCATTCCTATTCC (GTT)6 58 244 R:CCTACATTTGTGCTTTTGTTC Dma83 F:CTCTAAAGCCGACCTAACC (CTT)6 58 239 R:TGTCTCAACACAGCGAAAC Dma84 F:AAACTAACTCATCACCAG (TGT)6 55 283 R:AAACGACAGGAACTCAAT Dma85 F:CTCACTTTGACCCAACCAG (AGG)6 55 256 R:CCTTTCACCGAGACACCAG Dma131 F:TGCGGATGGGTGGTAGTGT (GGT)8 55 208 R:ATTGCTGGTAGTCGGTGGC Dma132 F:CCCAGTGAGACCAGAACCA (GCT)8 55 268 R:GACCCGTAGACAGGAGAGT Dma135 F:GTTGTTGTTTTTTTCCTT (GCA)9 55 301 R:CATCAGTCTGGCTTTATA Dma145 F:ACGATACAGCAGCCGAAG (TCA)10 60 197 R:AGTGATGTCGCCTCATAAAT 2.3 微卫星标记的种群遗传学评价
使用1个采集自南海东南部海域的长体圆鲹群体对筛选合格的微卫星标记进行种群遗传学评价。所有27个标记在测试群体中共检测到285个等位基因，等位基因数为5~17，Ho为0.342 9~0.857 1，平均为0.631 7；He为0.538 3~0.911 8，平均为0.7968。PIC为0.497~0.886，平均为0.780 9 (表3)，表明开发的微卫星位点具有较高的多态性。共有19个标记等位基因频率符合“哈迪-温伯格”平衡。连锁不平衡检测表明各位点间无连锁不平衡现象。
表 3 长体圆鲹微卫星标记的种群遗传学特征Table 3. Characteristics of microsatellite loci in D. macrosoma位点
locusN Na Ho He PHWE PIC Dma03 35 8 0.857 1 0.790 9 0.042 0 0.746 0 Dma07 35 9 0.771 4 0.855 9 0.553 4 0.825 0 Dma12 34 11 0.685 7 0.816 6 0.106 4 0.809 0 Dma15* 35 14 0.542 9 0.911 8 0.000 0 0.890 0 Dma22* 32 14 0.485 7 0.813 2 0.000 0 0.864 0 Dma23 34 17 0.714 3 0.864 8 0.003 3 0.868 0 Dma26* 34 14 0.428 6 0.869 1 0.000 0 0.870 0 Dma28 35 11 0.742 9 0.837 7 0.035 3 0.803 0 Dma36 35 9 0.685 7 0.786 3 0.054 5 0.743 0 Dma38* 34 12 0.628 6 0.851 2 0.000 0 0.849 0 Dma39* 33 12 0.514 3 0.782 0 0.000 0 0.793 0 Dma45 34 12 0.771 4 0.847 8 0.139 2 0.844 0 Dma51* 31 12 0.485 7 0.784 5 0.000 0 0.859 0 Dma54 35 15 0.771 4 0.900 6 0.139 0 0.878 0 Dma58 34 12 0.771 4 0.845 6 0.449 7 0.843 0 Dma64* 31 13 0.428 6 0.788 3 0.000 0 0.865 0 Dma72 35 8 0.628 6 0.713 9 0.155 7 0.659 0 Dma76 35 12 0.685 7 0.908 1 0.003 6 0.886 0 Dma81 35 8 0.628 6 0.795 0 0.045 7 0.752 0 Dma82 35 5 0.485 7 0.538 3 0.023 6 0.497 0 Dma83* 34 8 0.342 9 0.705 8 0.000 0 0.674 0 Dma84 35 7 0.542 9 0.704 8 0.005 9 0.649 0 Dma85 35 8 0.828 6 0.747 8 0.876 9 0.704 0 Dma131 35 9 0.628 6 0.717 6 0.646 6 0.661 0 Dma132 35 8 0.657 1 0.746 6 0.365 5 0.692 0 Dma135 34 10 0.742 9 0.840 5 0.241 8 0.835 0 Dma145 34 7 0.600 0 0.749 3 0.077 4 0.725 0 注：N. 有效样品数；Na. 等位基因数；Ho. 表观杂合度；He. 期望杂合度；PHWE. “哈迪-温伯格”平衡显著性检验P值；PIC. 多态信息含量；*. 经Bonferroni校正后显著背离“哈迪-温伯格”平衡 (校正P<0.001 85) Note: N. effective number of samples; Na. number of alleles; Ho. observed heterozygosity; He. expected heterozygosity; PHWE. Hardy–Weinberg probability test; PIC. polymorphism information content; *. significant deviation from HWE after Bonferroni's correction (adjusted P-value<0.001 85) 3. 讨论
3.1 高通量测序发掘微卫星序列的技术优势
传统微卫星标记开发方法耗时长、花费高、技术难度大。以磁珠富集法为例，标记开发过程中基因组DNA浓度、接头连接效率、富集过程中的杂交温度以及洗脱条件的控制等因素都会影响微卫星筛选的效率[24-25]，且最终获得的有效微卫星序列仅几百条[26-27]。相比较而言，高通量测序技术开发微卫星标记，省略了建库、克隆、筛选等，只需提取基因组DNA测序，利用生物信息学手段可直接获取微卫星序列，通常是传统方法获得微卫星序列数目的几百倍[15,28-29]，具有高效、便捷、准确的特点，能够满足短时间内大批量微卫星位点的开发需求，比如连锁图谱构建、QTL定位等[30-31]。
本次RAD-seq高通量测序共获得长体圆鲹基因组原始数据2.06 G，GC含量为41.43%，测序质量Q30达93.05%；共获得微卫星序列58 180条，一至六核苷酸重复微卫星位点70 508个。说明测序质量稳定高效，并获得了数量庞大、类型丰富的长体圆鲹微卫星序列，可用于后续长体圆鲹微卫星标记的大规模开发和相关遗传学研究。
3.2 不同核苷酸重复微卫星位点比较
本次高通量测序结果表明在长体圆鲹微卫星位点中二核苷酸重复为主要重复类型，AC/GT类重复数量最为丰富，GC/CG重复较为少见。此结果与大量水产动物微卫星位点研究结果相一致[32-34]，差异仅在于比例多少，以及除二核苷酸重复占主要地位外其他核苷酸重复的含量差异。熊良伟等[33]对中华鳑鲏(Rhodeus sinensis)微卫星的分析中，二核苷酸占总微卫星位点的53.59%，其次为单核苷酸重复，二核苷酸重复中AC/GT类占60.63%，GC/CG仅占0.32%。在裸体异鳔鳅鮀 (Xenophysogobio nudicorpa)中[32]，二核苷酸重复占总微卫星位点比例高达83.15%，AC/GT类重复占49.36%，GC/CG重复仅有4个。
多数鱼类开发的微卫星标记以二核苷酸重复为主，但研究表明，在人类基因组中三核苷酸重复序列与遗传疾病的发生有关，并且具有较高的多态性和遗传稳定性[35]。部分学者对三、四核苷酸重复微卫星标记的研究结果存在差异。房祖业等[28]对大刺鳅 (Mastacembelus armatus) 二、三、四核苷酸重复微卫星标记的筛选发现二核苷酸重复较三、四核苷酸重复具有更高的筛选效率和多态性；鲁翠云等[17]、谭照君等[18]、李文升等[19]的研究认为三、四核苷酸具有更高的多态性和分型效果。长体圆鲹二、三核苷酸的筛选效率分别为22.2%和29.0%，PIC分别为0.827 4和0.687 7 (表4)。就筛选效率而言，三核苷酸重复略高于二核苷酸重复，但二者相差不大。PIC为衡量种群遗传变异程度的重要指标[36]，二核苷酸重复多态性明显高于三核苷酸重复。本文中长体圆鲹二核苷酸重复筛选效率低于三核苷酸重复，但多态性二核苷酸重复明显高于三核苷酸重复。因此，筛选效率和多态性的差异可能由种属差异或其他多种因素导致。
表 4 长体圆鲹二、三核苷酸重复微卫星标记的比较Table 4. Comparison on di- and trinucleotide-repeated microsatellite loci in D. macrosoma序列
sequence引物数
primer number重复次数
repeat number筛选效率
efficiencyPIC 二核苷酸重复
di-nucleotide-repeated18 9~14 22.2% 0.827 4 三核苷酸重复
tri-nucleotide-repeated9 6~10 29.0% 0.687 7 通过筛选的27对引物中18个位点为二核苷酸重复，重复次数为9~14次不等；9个位点为三核苷酸重复，重复次数为6~10次不等，符合Ellegren[37]提出的真核生物微卫星位点重复大部分在30次重复以下。但Ellegren[37]认为二核苷酸重复以15~19次为主，本文中高通量测序获得的二核苷酸重复主要在6~15次。基于Weber[38]的研究结果，重复次数高的微卫星在种群中表现出的多态性较高，龚小玲等[39]对澳洲鳗鲡 (Anguilla australis) 进行标记开发时发现，微卫星重复序列的重复次数过高会影响PCR效果，应选择居中的重复次数为宜。长体圆鲹二核苷酸重复PIC为0.827 4，具有较高多态性，表明选择6~15次的二核苷酸重复是合适的。
3.3 微卫星标记的种群遗传学特征
群体杂合度的高低反映了群体在多个基因座上的遗传变异及群体遗传多样性丰富度[19]。本研究中长体圆鲹中沙群体的平均Ho为0.631 7，平均He为0.796 8，说明长体圆鲹该群体的遗传多样性较高。平均Ho和He存在差异，说明存在杂合子缺失或者纯合子过剩的情况。PIC也是衡量群体遗传多样性的重要指数，Botstein等[36]认为基因标记PIC>0.5为高度多态位点，0.25<PIC<0.5为中度多态位点，PIC<0.25为低度多态性位点，通常不作为遗传多样性分析。本文中长体圆鲹位点除1个为中度多态外，其他位点均为高度多态位点。表明开发所得的长体圆鲹微卫星标记在中沙群体中具有较好的遗传稳定性和丰富的遗传多样性。
在所有27个位点中有8个位点偏离了“哈迪-温伯格”平衡，这些位点不适合进一步的遗传分析。近亲杂交、无效等位基因、种群退化和自然选择等因素皆可能导致微卫星位点偏离HWE[15]。
-
![]()
图 1 相同驯料方法4种不同驯料时间实验设计
W13. 13~22 dph；W16. 16~25 dph；W19. 19~26 dph；W22. 22~33 dph；图7、图8同此
Figure 1. Experimental design of same weaning regime started on four different days post hatching
AN. artemia nauplii; MF.mixed feeding; WC. weanning complete; W13. 13−22 dph; W16. 16−25 dph; W19. 19−26 dph; W22. 22−33 dph; the same cases in Fig.7 and Fig.8.
![]()
图 2 不同密度轮虫投喂下9 dph红鳍笛鲷仔鱼的特定生长率和成活率
不同小写字母表示差异显著 (P<0.05)，后图同此
Figure 2. Specific growth rate and survival rate of 9 dph L. erythopterus juveniles at different rotifer densities
Different lowercase superscripts indicate significant difference (P<0.05). The same case in the following figures.
![]()
图 3 不同密度轮虫投喂下3 dph、5 dph、7 dph、9 dph红鳍笛鲷仔鱼摄食量
Figure 3. Food consumption of 3 dph, 5 dph, 7 dph and 9 dph L. erythopterus juvenles at different rotifer densities
![]()
图 4 不同密度轮虫投喂下7 dph红鳍笛鲷仔鱼日摄食量和轮虫利用率
Figure 4. Daily ingestion and utilization rate of rotifer for L. erythopterus juveniles on 7 dph at different rotifer densities
![]()
图 5 不同密度轮虫投喂下红鳍笛鲷仔鱼7 dph和9 dph的饵料选择系数
Figure 5. Food selection of 7 dph and 9 dph L. erythopterus juveniles at different rotifer densities
![]()
图 6 不同驯料时间红鳍笛鲷的特定生长率和成活率
Figure 6. Specific growth rate and survival rate of L. erythopterus juveniles at different weaning time
![]()
图 7 不同驯料时间红鳍笛鲷RNA/DNA比率
Figure 7. RNA/DNA ratio of L. erythopterus at different weaning time
-
[1] HU J, LIU Y, MA Z, et al. Feeding and development of warm water marine fish larvae in early life[M]//YUFERA M. Emerging issues in fish larvae research. Cham, Switzerland: Springer International Publishing AG, 2018: 275-296.
[2] 郭浩宇, 张秀梅, 张宗航, 等. 许氏平鲉仔、稚鱼的摄食特性及幼鱼胃排空率[J]. 水产学报, 2017, 41(2): 285-296. [3] 于欢欢, 李炎璐, 陈超, 等. 棕点石斑鱼(♀)×鞍带石斑鱼(♂)杂交F1仔、稚、幼鱼的摄食与生长特性分析[J]. 中国水产科学, 2015, 22(5): 968-977. [4] 刘利平, 刘登攀, 蒲金成, 等. 日本鳗鲡仔鱼的开口饵料和行为特征[J]. 水产学报, 2017, 41(5): 703-710. [5] 王晓龙, 温海深, 张美昭, 等. 花鲈初孵仔鱼饥饿不可逆点的确定及摄食节律研究[J]. 中国海洋大学学报(自然科学版), 2017, 47(5): 57-64. [6] 高小强, 洪磊, 刘志峰, 等. 美洲西鲱仔鱼不可逆点及仔、稚鱼摄食特性研究[J]. 水产学报, 2015, 39(3): 392-400. [7] 杨育凯, 虞为, 林黑着, 等. 豹纹鳃棘鲈仔鱼饥饿实验和不可逆点研究[J]. 南方水产科学, 2017, 13(6): 90-96. [8] 隋延鸣, 庄亚润, 周凯, 等. 轮虫、卤虫无节幼体及蝇蛆三种开口饵料对黄颡鱼仔鱼生长、存活及免疫酶活性的影响[J]. 水产学杂志, 2018, 31(2): 20-24. [9] MA Z, GUO H, ZHANG D, et al. Food ingestion, consumption and selectivity of pompano, Trachinotus ovatus (Linnaeus 1758) under different rotifer densities[J]. Aquacult Res, 2015, 46(11): 2593-2603.
[10] MA Z, QIN J G, HUTCHINSON W, et al. Food consumption and selectivity by larval yellowtail kingfish Seriola lalandi cultured at different live feed densities[J]. Aquacult Nutr, 2013, 19(4): 523-534.
[11] TANAKA Y, SATOH K, YAMADA H, et al. Assessment of the nutritional status of field-caught larval Pacific bluefin tuna by RNA/DNA ratio based on a starvation experiment of hatchery-reared fish[J]. J Exp Mar Bio Ecol, 2008, 354(1): 56-64.
[12] PILAR OLIVAR M, DIAZ M V, CHÍCHARO M A. Tissue effect on RNA:DNA ratios of marine fish larvae[J]. Scientia Marina, 2009, 73(S1): 171-182.
[13] CHEN B N, QIN J G, CARRAGHER J F, et al. Deleterious effects of food restrictions in yellowtail kingfish Seriola lalandi during early development[J]. Aquaculture, 2007, 271(1/2/3/4): 326-335.
[14] 刘皓, 张玉红, 罗杰, 等. 红鳍笛鲷(Lutjanus erythopterus)卵巢发育的组织学研究[J]. 海洋与湖沼, 2016, 47(1): 269-275. [15] 程大川, 马振华, 江世贵. 红鳍笛鲷仔、稚鱼异速生长[J]. 水生生物学报, 2017, 41(1): 206-213. [16] 李加儿, 区又君, 刘匆. 红笛鲷和卵形鲳鲹鳃的扫描电镜观察与功能探讨[J]. 海洋水产研究, 2007, 28(6): 45-50. [17] CUI K, CHENG D C, MA Z H, et al. Ontogenetic development of digestive enzymes in larval and juvenile crimson snapper Lutjanus erythopterus (Bloch 1790)[J]. Aquacult Res, 2017, 48(8): 4533-4544.
[18] CUI K, FU Z, CHENG D, et al. Development of immune functionality in larval and juvenile crimson snapper Lutjanus erythropterus (Bloch 1790)[J]. Aquacult Rep, 2018, 10: 1-7.
[19] HOPKINS K D. Reporting fish growth: a review of the basics[J]. J World Aquacult Soc, 1992, 23(3): 173-179.
[20] ELLIOTT J M, PERSSON L. The estimation of daily rates of food consumption for fish[J]. J Anim Ecol, 1978, 47(3): 977-991.
[21] IVLEV V S. Experimental ecology of the feeding of fishes[J]. Copeia, 1962, 1: 234-236.
[22] KLUMPP D W, von WESTERNHAGEN H. Nitrogen balance in marine fish larvae: influence of developmental stage and prey density[J]. Mar Biol, 1986, 93(2): 189-200.
[23] ZEHRA S, KHAN M A. Dietary histidine requirement of fingerling Catla catla (Hamilton) based on growth, protein gain, histidine gain, RNA/DNA ratio, haematological indices and carcass composition[J]. Aquacult Res, 2016, 47(4): 1028-1039.
[24] MA Z, QIN J G, NIE Z. Morphological changes of marine fish larvae and their nutrition need[M]//NIROOMAND P K. Larvae: morphology, biology and life cycle. New York: Nova Science Publisher, Inc., 2012: 1-20.
[25] WELLINGTON C G, MAYER C M, BOSSENBROEK J M, et al. Effects of turbidity and prey density on the foraging success of age 0 year yellow perch Perca flavescens[J]. J Fish Biol, 2010, 76(7): 1729-1741.
[26] RUZICKAJ J, GALLAGER S M. The importance of the cost of swimming to the foraging behavior and ecology of larval cod (Gadus morhua) on Georges Bank[J]. Deep-Sea Res II, 2006, 53(23/24): 2708-2734.
[27] THEILACKER G H, PORTER S M. Condition of larval walleye pollock, Theragra chalcogramma, in the western Gulf of Alaska assessed with histological and shrinkage indices[J]. Fish Bull, 1995, 93: 333-344.
[28] HAMZA N, MHETLI M, KESTEMONT P. Effects of weaning age and diets on ontogeny of digestive activities and structures of pikeperch (Sander lucioperca) larvae[J]. Fish Physiol Biochem, 2007, 33(2): 121-133.
[29] DIAZ M V, PAJARO M, OLIVAR M P. Nutritional condition of Argentine anchovy Engraulis anchoita larvae in connection with nursery ground properties[J]. Fish Res, 2011, 109(2/3): 330-341.
[30] MA Z, ZHENG P, GUO H, et al. Effect of weaning time on the performance of Trachinotus ovatus (Linnaeus 1758) larvae[J]. Aquacult Nutr, 2015, 21(5): 670-678.
[31] MA Z H, QIN J G, HUTCHINSON W, et al. Responses of digestive enzymes and body lipids to weaning times in yellowtail kingfish Seriola lalandi (Valenciennes, 1833) larvae[J]. Aquacult Res, 2014, 45(6): 973-982.
-
期刊类型引用(8)
1. 刘家光. 鳗鱼骨酶解及美拉德反应制备优质风味物工艺研究. 福建轻纺. 2024(03): 11-16 . 
百度学术
2. 王扬铎，苏永昌，王晓燕，施文正，刘智禹. 仿刺参不同部位ACE抑制活性分析及活性肽制备工艺优化. 食品工业科技. 2024(10): 187-197 . 百度学术
3. 郭星晨，李华鑫，张钰璇，马金璞，杨具田，李贞子，高丹丹. 三种描述符对食源性血管紧张素转化酶抑制二肽定量构效关系研究. 食品与发酵工业. 2024(10): 290-301 . 百度学术
4. 侯蕊，梁璋成，林晓婕，林晓姿，张秀红，何志刚. 红曲糟源酶解蛋白肽的功能性评价. 福建农业学报. 2024(03): 354-361 . 百度学术
5. 娄肖肖，马洪鹏，高姣姣，邵伟，郑楠，赵艳坤. 乳源生物活性肽生物学功能研究进展. 乳业科学与技术. 2024(04): 31-37 . 百度学术
6. 易一帆，王越，程宽，张敏，刘学平，丁翠翠. 胶原蛋白肽的制备、生物活性及应用研究进展. 食品工业科技. 2024(23): 396-404 . 百度学术
7. 王莹莹，侯喜龙，金桥，曲敏，李伟，佟长青. 大鲵骨胶原肽的研究及应用. 农产品加工. 2024(23): 73-77+84 . 百度学术
8. 赵欣，成晓瑜，张顺亮，王乐，岳宜静，刘博文，王辉. 不同脱苦处理对骨胶原蛋白肽品质及抗氧化特性的影响. 肉类研究. 2023(11): 13-20 . 百度学术
其他类型引用(5)
下载:
计量
- 文章访问数:
- HTML全文浏览量:
- PDF下载量:
- 被引次数: 13
粤公网安备 44010502001741号
