Morphological development of embryos and larvae of Mystus macropterus
-
摘要:
大鳍鳠 (Mystus macropterus) 是一种亟待开发的优质养殖品种。探明养殖大鳍鳠早期发育阶段的特征和规律，可为其人工催产和苗种培育提供理论依据。采用显微镜观察法，对其胚胎及仔稚鱼发育各时期的形态特征和发育特点进行了研究，并对其形态学指标进行了测量，推算出仔稚鱼阶段生长特性的回归方程。结果表明：大鳍鳠受精卵呈椭圆形、黄褐色，卵径为 (3.46±0.3) mm；平均水温26.2 ℃下，大鳍鳠胚胎发育历时79 h，积温为2065.26 h·℃，根据发育时胚胎的典型外部形态特征，将其分为受精卵、胚盘期、卵裂期、囊胚期、原肠期、神经胚期、器官分化期和出膜期共8个阶段28个时期；初孵仔鱼全长 (7.31±0.52) mm，仔鱼期历时10 d，仔鱼前期卵黄囊大部分被吸收，全长特定生长率 (SGRL) 为9.79%，体质量特定生长率 (SGRm) 为12.31%，出膜后7~8 d仔鱼大部分开口，9 d卵黄囊消失，仔鱼后期SGRL为4.15%，SGRm为12.25%。稚鱼期历时11 d，各鳍基本长成，体型和体色接近成鱼，SGRL为4.30%，SGRm为13.17%，全长22~24 mm。统计模型显示，大鳍鳠仔稚鱼的生长属于异速生长类型。
Abstract:Mystus macropterus is a high-quality aquaculture variety that urgently needs to be developed. Exploring the characteristics and patterns of early developmental stages of cultured M. macropterus can provide basic data for its seedling cultivation and large-scale efficient seed production. We studied the morphological characteristics and developmental features of its embryos and juveniles by microscopic observation, and measured the morphological indicators to derive regression equations for the growth characteristics of juveniles. The results show that the fertilized eggs of the fish were oval and yellowish brown, with a diameter of (3.46±0.3) mm. At incubation temperature of 26.2 ℃, the embryos took 79 h to develop for eight stages with 28 periods, including fertilization, cleavage, blastocyst, gastrula, neuroembryo, organ differentiation and hatching, and the total accumulative temperature was 2065.26 h·℃. The total length of the larvae was (7.31±0.52) mm, the larvae stage lasted for 10 d, and most of the yolk sac was absorbed at early stage of larvae. The full-length specific growth rate (SGRL) was 9.79%, and the weight-specific growth rate (SGRm) was 12.31%. Most of the larvae opened 7–8 d after membrane emergence, and the yolk sac disappeared on 9th day. The SGRL and SGRm were 4.15% and 12.25% at later stage of larvae, respectively. The juvenile stage lasted for 11 d. The fins were basically mature, and the body size and body color were close to those of adult fish. SGRL was 4.30%, SGRm was 13.17%, and the total length was 22–24 mm. The statistical model shows that the growth of the larvae and juveniles is of allometric type.
-
黄鳝 (Monopterus albus) 属于硬骨鱼纲、合鳃鱼目、合鳃鱼科，为淡水肉食性鱼类[1]，主要分布于中国、日本、印度，以及泰国、印度尼西亚、缅甸等一些东南亚国家[2]。在我国，黄鳝分布在除青藏高原以外的广泛区域，是我国重要的名特优淡水经济鱼类之一，其肉质鲜美，具有很高的营养价值[3]。根据《中国渔业统计年鉴》数据，湖北、江西、安徽、湖南、四川是我国黄鳝养殖主产区，在2016—2020年间年均产量超30万吨，其中湖北省的产量最高，其年产量占全国的40%以上[4-8]。黄鳝通常在完成第一次产卵后会由雌性转为雄性，导致适龄繁育的雌性亲本数量稀缺，极大限制了其产业发展[9]。目前，我国黄鳝养殖主要依赖于野生苗种捕捞，质量参差不齐，这种高强度捕捞对黄鳝养殖产业的可持续性及种质资源造成了极大破坏。在育种方面，目前尚无国家认定的黄鳝养殖品种，因此开展黄鳝品种选育、筛选优良品系、进行家系管理，是当下黄鳝育种工作中的迫切需求。笔者课题组已实现黄鳝全人工繁育技术的突破，已应用于规模化育苗，并开展了家系选育工作，为下一步黄鳝良种的选育奠定了基础。
在黄鳝后续良种选育工作中，亲子鉴定及系谱管理技术是避免其近交衰退的关键技术。传统家系选育工作存在繁琐及环境差异等问题，传统标记手段也存在幼体标记困难及造成应激、损伤、脱落等问题。开发高效亲子鉴定技术可为黄鳝良种选育和种质资源保护及利用提供技术支持。亲子鉴定又称亲权鉴定，可以通过遗传特征相互验证亲本和子代之间是否存在亲缘关系。在育种生产中开展亲子鉴定可有效避免近交导致的种质退化问题。微卫星序列 (Simple sequence repeats, SSR) 具有高度的多态性、稳定性、呈共显性遗传、数量丰富且遵循孟德尔遗传定律的优点，在水生动物亲子鉴定中应用广泛。文萍等[10]使用2组微卫星多重PCR体系对黄喉拟水龟 (Mauremys mutica) 进行亲子鉴定，在父本信息未知时，母本的鉴定准确率为87%；苗贵东等[11]通过2组微卫星多重PCR体系对大菱鲆 (Scophthalmus maximus) 7个家系开展了亲子鉴定，准确率达到98.26%；谢敏敏等[12]通过建立2组微卫星多重PCR体系，对鼋 (Pelochelys cantorii) 家系进行亲子鉴定分析，得到10个位点的双亲排除率为98%；Dong等[13]为了评估中国明对虾 (Fenneropenaeus chinensis) 养殖群体亲缘鉴定的可行性，使用5个微卫星位点，获得了92.9%的鉴定成功率；辛苗苗等[14]使用7个微卫星位点对多倍体中华鲟 (Acipenser sinensis) 进行亲子鉴定分析，得到双亲累计排除率为99.99%；Yang等[15]利用10个微卫星位点对鳜 (Siniperca chuatsi) 5个全同胞家系进行了亲子鉴定，得到的实际鉴定率为95%，且随机选择了69个后代进行双盲试验，结果表明95%的子代可以正确聚类；Zhu等[16]筛选出了14个用于评价斑节对虾 (Penaeus monodon) 亲子关系的微卫星标记，并用其中多态性最高的6个微卫星标记获得了99%的鉴定准确率。
目前，黄鳝已开发出较多的微卫星标记，如刘臻等[17]通过富集文库法开发了20个微卫星标记，设计了8对4碱基重复高度多态性的微卫星引物，分析了3种不同性别表型黄鳝群体的遗传多样性；Li等[18]使用富集文库法分离了16对微卫星引物，分析了两个黄鳝自然群体的遗传多样性和群体结构；Zhang等[19]通过RAD (Restriction-site associated DNA tags，简化基因组) 测序在基因组中鉴定出了9 897个微卫星，并筛选出28对多态性引物。但上述研究开发的微卫星标记仅用于遗传多样性分析，目前尚未见利用微卫星多重PCR体系对黄鳝进行亲子鉴定的相关研究报道。本研究利用笔者课题组前期发表的高度多态性微卫星标记[20]和基于黄鳝基因组数据开发的新的微卫星标记[21]，结合繁育的黄鳝全同胞家系，建立了基于微卫星多重PCR的亲子鉴定技术，为后续黄鳝育种以及种质资源保护提供基础。
1. 材料与方法
1.1 实验材料
黄鳝样本取自中国水产科学研究院长江水产研究所窑湾基地。雌、雄亲本按照1∶1比例得到11个全同胞家系，雌、雄亲本完成繁殖后取背部肌肉于−80 ℃保存，每个家系在孵化出膜至30日龄随机选取10尾于−80 ℃保存。随机挑选4个黄鳝家系子代各10尾和随机候选亲本4尾，同样−80 ℃保存。
1.2 基因组DNA的提取
基因组DNA采用苯酚-氯仿法[22]提取，用1% (w) 琼脂糖凝胶电泳以及分光光度计检测DNA浓度及完整性。将合格的DNA用双蒸水 (ddH2O) 稀释至100 ng·μL−1，于−20 ℃保存。
1.3 多态性微卫星引物的筛选
黄鳝全基因组数据来源于中国科学院国家基因研究中心数据库 (CRA003062)。使用MISA对黄鳝全基因组序列中的完美型SSR进行检索[23]，检索标准与课题组前期筛选微卫星的标准一致。挑取了148个SSR位点，同时提取SSR位点前后200 bp的序列，使用Primer premier 3.0软件设计引物。随机选取6尾黄鳝的混合DNA为模板，进行微卫星引物的扩增效果及初步的多态性验证，共获得35对扩增效果理想的引物。另挑选课题组前期发表的具有高度多态性的15对微卫星引物，共50对微卫星引物，以参与繁殖的22尾黄鳝亲本DNA为模板进行筛选，选择扩增等位基因数较多、出峰质量较好的位点。
1.4 多重PCR体系的构建及优化
对筛选获得的16对微卫星引物进行荧光修饰 (FAM、HEX、ROX、TAMRA) ，由武汉天一辉远科技有限公司合成，以亲本DNA为模板进行多重PCR构建。构建过程中避免产生引物二聚体、发卡结构和错配等现象，优化模板DNA浓度、退火温度、延伸时间等反应条件。根据片段大小将引物分为两组，两组引物等比混合开展预实验，根据出峰结果调整引物浓度，最终确定最佳的多重PCR体系。
PCR反应体系10 μL，包含2×PCR Mix 5 μL，10 pmol·μL−1上下游引物各0.5 μL，基因组DNA 1 μL以及灭菌超纯水3 μL。扩增程序为：95 ℃预变性5 min；95 ℃变性30 s，62~53 ℃退火30 s，72 ℃延伸30 s，共10个循环，其中退火温度每次循环下降1 ℃；95 ℃变性30 s，52 ℃退火30 s，72 ℃延伸30 s，共25个循环；72 ℃延伸20 min。使用1% (w) 的琼脂糖凝胶电泳检测扩增产物的完整性，然后由ABI3730XL测序仪进行毛细管电泳检测 (由武汉天一辉远生物科技有限公司完成)，使用GeneMarker[24]读取各位点等位基因大小。
1.5 数据统计及分析
将整理的亲本与子代基因型结果导入Cervus 3.0[25]进行等位基因频率分析、模拟分析和亲子鉴定分析。模拟分析设置模拟子代数为10 000个，候选亲本数为父母本各11对，对100%的候选亲本进行抽样模拟分析，允许1%的分型错误，置信度为95%。在性别已知和未知的情况下，基于16个微卫星位点模拟分析不同候选亲本数情况下的鉴定率，设置模拟子代数目为10 000个，候选亲本数分别为20、50、100、150和200对，抽样比例为100%，允许1%的分型错误，置信度为95%。根据多态信息含量大小将16个微卫星位点排序，分析微卫星位点数与鉴定率之间的关系。统计所有子代样本在某一位点出现的所有基因型，单个样本在该位点存在对应基因型则记为1，否则记为0，构建基因型0/1矩阵，并根据该矩阵利用NTSYS软件[26]对110尾家系子代进行聚类分析。利用NTSYS软件对4个家系共40尾子代进行聚类分析，验证微卫星多重PCR体系的准确率。
其中单个位点的非亲排除率 (NEP) 由Cervus 3.0软件计算得到，累积非亲排除率 (CEP) 计算公式如下：
$$ R_{\rm{CEP}}=1-{\prod} _{i=1}^{n}{(1-R}_{\mathrm{N}\mathrm{E}\mathrm{P},i}) $$ (1) 式中：RCEP为累积非亲排除率；RNEP, i为第i个微卫星位点的NEP值；n为微卫星位点排除概率的数目。CE-1P对应用NE-1P计算，CE-2P对应用NE-2P计算，CE-PP对应用NE-PP计算。
2. 结果
2.1 微卫星的筛选与多重PCR建立
从黄鳝基因组预测的148个微卫星位点，筛选获得了扩增效果较好的微卫星引物35对。进而以30尾亲本为模板，从50对微卫星引物中筛选得到16对多态性高、扩增效果好的微卫星引物。
将筛选的16对微卫星引物进行不同组合并优化反应程序、调整引物浓度比例，最终确定2 组多重PCR (表1)。2组多重PCR的扩增效率都较好，分型结果符合预期。A组多重PCR体系各位点如见图1。
表 1 黄鳝 2 组微卫星多重PCR的引物信息Table 1. Primer information of two multiplex PCR sets of microsatellites in M. albus分组
Group位点
Locus重复单元
Reapet motif退火温度
Tm/℃引物序列
Primer sequence荧光标记
Fluorescence label片段长度
Size/bp比例
RatioA Mta027 (GCA)7 52 F: CTGCGGTAACAAGCGTATCA
R: CTGGGGATCCCAGTCAAACTFAM 147~154 3 Mta141 (GT)17 52 F: CAGAGATGGTCGACTGGTGA
R: AGAGCCATGGGAAGCACTTAFAM 185~312 2 Mta031 (ATG)8 52 F: AGAGGCAGCTCATGGACACT
R: TTCATGGCTGAGCTGACTTGFAM 242~265 3 Ma28 (ATC)7 52 F: GGGATTGTCTGGAATGCTGT
R: GTGAACCCTGAACAGACGGTHEX 144~168 3 Ma39 (GAT)8 52 F: AGGTGAAGGGGAACACACTG
R: TTGCCCTGTCCATTTTTCTCHEX 193~224 3 Ma71 (TCA)13 52 F: TCACATTGCCCAGAGAACAG
R: GGCAGCATCAAGAGACCTTCHEX 258~269 2 Ma24 (CTG)9 52 F: TCGTCTGGCTCAGAGGAGAT
4R: AGCTGCAGCAGGAGGAATACTAMRA 125~143 7 Mta026 (CAG)8 52 F5: TGAACATCCCTTCTTCCACC
R: GTCATGGGTTTGTGTTGCTGTAMRA 171~182 4 Ma85 (TGC)9 52 F: GTGGAAGAAGCTGGATGAGC
R: CACATGGGGTGTTCTCACTGTAMRA 208~236 5 Mta013 (AGG)9 52 F: GAGCGTCTTTTCCATCCTTG
R: CCTACTGCTGCTTCTGGTCCROX 131~150 4 Mta025 (GCT)7 52 F: GCCAAGCAAACGTGTGAGTA
R: AACGAGCTGCGTGTTAAGGTROX 172~184 3 Ma78 (TAAGA)5 52 F: CTGACGGTTGTGTGTAACGG
R: CACGTACGCCGGTTAAACTTROX 214~235 3 B Mta033 (TGT)7 52 F: TGGAGGCAGGAGAAGGAGTA
R: TGAAAGACTCCTCGCAACCTFAM 158~170 1 Mta021 (TGT)7 52 F: AGAGTCTTGCTCGGTTTCCA
R: AGAGTCTTGCTCGGTTTCCAFAM 201~210 1 Ma5 (TCTAC)7 52 F: TTCAGGTTCGGGACTTTGAC
R: AGACGGTTGTGTCAGGAAGGFAM 231~256 1 Mta148 (GT)15 52 F: TTTTCCCAGCAGCTGATTTC
R: CATGACAACAGGACGCAAACHEX 158~174 1 2.2 亲子鉴定分析
16个微卫星位点在黄鳝22尾亲本和110尾子代中的遗传参数如表2所示。等位基因数 (Na) 为3~9，平均等位基因数为5.562；观测杂合度和期望杂合度分别为0.421~0.771和0.382~0.797，平均观测杂合度和期望杂合度分别为0.627和0.619；多态信息含量 (PIC) 为0.347~0.766，平均多态信息含量为0.564，其中Ma71、Ma78、Mta021、Mta027和Mta033表现为中度多态性，其余11个位点均表现为高度多态性。16个位点均未偏离Hardy-Weinberg平衡。
表 2 16 个微卫星位点的遗传参数Table 2. Genetic parameters of 16 SSR loci位点
Locus等位基因数
Na观测杂合度
Ho期望杂合度
He多态信息含量
PICNE-1P NE-2P NE-PP 无效等位基因频率
FNullMa5 5 0.758 0.704 0.647 0.723 0.555 0.380 −0.037 6 Ma24 6 0.550 0.546 0.513 0.833 0.662 0.472 −0.018 2 Ma28 9 0.679 0.660 0.613 0.747 0.574 0.384 −0.017 3 Ma39 9 0.771 0.797 0.766 0.570 0.393 0.206 0.004 0 Ma71 5 0.485 0.494 0.444 0.873 0.732 0.579 0.005 4 Ma78 4 0.489 0.490 0.401 0.880 0.780 0.662 −0.001 5 Ma85 8 0.583 0.717 0.676 0.682 0.504 0.308 0.106 6 Mta013 4 0.583 0.647 0.586 0.775 0.615 0.444 0.046 6 Mta021 4 0.667 0.578 0.485 0.830 0.714 0.574 −0.076 9 Mta025 4 0.652 0.588 0.522 0.822 0.676 0.517 −0.055 4 Mta026 4 0.636 0.621 0.545 0.799 0.657 0.498 −0.020 5 Mta027 3 0.432 0.382 0.347 0.928 0.803 0.675 −0.086 5 Mta031 8 0.756 0.746 0.714 0.640 0.457 0.259 −0.019 3 Mta033 4 0.606 0.509 0.430 0.870 0.756 0.625 −0.100 9 Mta141 8 0.735 0.764 0.732 0.620 0.438 0.244 0.018 8 Mta148 5 0.644 0.666 0.603 0.760 0.600 0.430 0.008 3 双亲未知时单个亲本的累积非亲排除率 Combined non-exclusion probability (First parent): 0.013 807 65 已知单亲时另一个亲本的累积非亲排除率 Combined non-exclusion probability (Second parent): 0.000 339 45 双亲未知时父母本组合的累积非亲排除率 Combined non-exclusion probability (Parent pair): 0.000 001 26 注：NE-1P代表双亲基因型未知时，单个位点的单亲排除率；NE-2P表示已知单亲基因型时，另一亲本单个位点的排除率；NE-PP表示双亲基因型未知时，单个位点的双亲排除率。 Note: NE-1P represents exclusion probability and of the first parent given only the genotype; NE-2P represents exclusion probability for one candidate parent given the genotype of a known parent of the opposite sex; NE-PP represents exclusion probability probability of a parent pair given only the genotype. 将16个微卫星位点根据多态信息含量由大到小排序探究微卫星数量与鉴定率之间的关系。结果如图2所示，随着微卫星位点数的增加，亲子鉴定的模拟鉴定率和实际鉴定率逐渐上升。当位点数为8时，亲子鉴定的模拟鉴定率和实际鉴定率均达到95%。当位点数继续增加，模拟鉴定率无限接近100%，实际鉴定率则保持在95%左右。实际亲子鉴定结果显示，110尾子代中仅有6尾子代的三联LOD≤0，无法准确匹配到真实的父母本，运用16个微卫星位点的实际鉴定成功率为95%。
利用构建的2组多重PCR体系计算黄鳝11个家系单个位点的排除概率和累计排除概率，结果见表2。NE-1P、NE-2P、NE-PP分别为0.570~0.928、0.393~0.803、0.143~0.675，CE-1P、CE-2P、CE-PP分别为0.999 999 99、0.999 999 91、0.999 964 76，模拟亲本已知的情况准确率为99.96%，置信度为95%。说明在理想情况下，使用这16个微卫星标记进行黄鳝的亲子鉴定可以在置信度95%的情况下达到99.96%的准确率。
此外，分别开展了在候选亲本性别已知和未知的情况下进行不同候选亲本数 (20、50、100、150、200对) 的模拟分析。结果如图3所示，当候选亲本数为20对，在候选亲本性别已知或未知的情况下，父本、母本和父母本联合都能够得到99%以上的鉴定率。当候选亲本数为50对时，无论父母本性别是否已知，父母本联合的鉴定率也都在95%左右，单亲本的鉴定率也在93%左右。当候选亲本数为100对时，无论父母本性别是否已知，父母本联合的鉴定率保持在95%以上，但单亲本的鉴定率下降到了90%以下。当候选亲本数为150对时，在候选亲本性别已知和未知的情况下，父母本联合仍然能够得到95%以上的鉴定率。当候选亲本数继续增加至200对时，鉴定率在父母本联合且性别已知的情况下依然有较高的鉴定率 (94.88%)，而在性别未知的情况下，鉴定率明显下降 (89.56%)。单亲本的鉴定率在亲本数为200对时下降至80%左右。
![]()
图 3 黄鳝16个微卫星标记基于不同亲本数的模拟分析Figure 3. Simulation analysis of 16 pairs microsatellite markers of M. albus based on different candidate parent2.3 家系聚类分析
NTSYS聚类分析的结果如图4所示，结果表明11个家系的子代仅有2尾无法正确聚类，分别为个体1-07和101-07，其余子代个体均可以正确聚类，鉴定准确率为98.18%。
![]()
图 4 110个黄鳝家系聚类分析图注：右侧数字代表家系，其中1-07和101-07为聚类错误的子代编号。Figure 4. Cluster analysis diagram of 110 families of M. albusNote: The numbers on the right represent families, and 1-07 and 101-07 are clustered incorrectly.2.4 群体遗传多样性分析
利用两组多重PCR体系对黄鳝11个子代家系群体和全部候选亲本群体独立进行遗传多样性分析，结果见表3。11个子代家系群体的平均等位基因数为0.489，平均观测杂合度和期望杂合度分别为0.632和0.489，平均多态信息含量为0.394。比较候选亲本群体和11个子代家系群体的遗传多样性，发现候选亲本群体的平均等位基因数、观测杂合度和多态信息含量均要高于11个子代家系群体。
表 3 黄鳝子代和亲本群体的遗传多样性Table 3. Genetic diversities of offsprings and parents of M. albus群体
Population平均等位基因数
Mean Na平均观测杂合度
Mean Ho平均期望杂合度
Mean He平均多态信息含量
Mean PIC子代-1 Offspring-1 2.750 0.619 0.495 0.408 子代-41 Offspring-41 2.813 0.688 0.554 0.453 子代-68 Offspring-68 2.375 0.669 0.475 0.375 子代-70 Offspring-70 2.063 0.488 0.386 0.301 子代-77 Offspring-77 2.563 0.719 0.536 0.436 子代-79 Offspring-79 2.125 0.575 0.430 0.334 子代-81 Offspring-81 2.563 0.806 0.574 0.457 子代-85 Offspring-85 2.125 0.494 0.384 0.311 子代-88 Offspring-88 2.375 0.556 0.466 0.373 子代-96 Offspring-96 2.438 0.639 0.506 0.411 子代-101 Offspring-101 3.188 0.700 0.574 0.474 亲本 Parents 5.625 0.556 0.636 0.570 2.5 亲子鉴定的应用
利用另外4个家系的子代和候选亲本44尾个体进行聚类分析。结果如图5，表明仅有个体10-04和23-08无法正确聚类，子代聚类准确率为95%，并且所有4个候选亲本均能正确聚类，说明本实验建立的微卫星多重PCR体系可以用于黄鳝家系的亲子鉴定。
![]()
图 5 亲子鉴定应用聚类分析注：“-”前的数字代表家系、后面的数字代表子代标号，M代表家系对应的母本，F代表家系对应的父本，其中10-04为聚类错误的子代编号。Figure 5. Cluster analysis of parentage test applicationNote: The number before "-" represents the family; the number after "-" represents the offspring label; M represents the female parent corresponding to the family; F represents the male parent corresponding to the family, and 10-04 was clustered incorrectly.3. 讨论
3.1 微卫星多重PCR体系
自20世纪80年代Chamberlain等[27]成功建立多重PCR (Multiplex polymerase chain reaction, MPCR) 体系以来，以其高效、敏感、大幅降低成本、提高检测效率的优势，在基因分型、遗传分析等领域得到广泛应用[28-31]。建立多重PCR反应体系，需要保证同种荧光标记引物的扩增产物不会重叠，以避免基因型误读，从而实现多重优势。本实验筛选获得16个微卫星位点的平均He、平均Ho和平均PIC分别为0.627、0.619和0.564。根据Botstein的标准[32]，11个位点为高度多态性位点，5个位点为中度多态性位点，这一结果与在兰州鲇 (Silurus lanzhouensis) [33]和圆口铜鱼 (Coreius guichenoti) [34]中用于亲子鉴定的微卫星位点多态性水平相近，表明这16个位点具有丰富的遗传多样性，能够用于后续的亲子鉴定。张毅等[35]提出影响多重PCR的关键因素是体系中引物之间是否兼容以及引物浓度的比例。本实验将16对微卫星引物依据扩增产物的片段大小及初步实验，在退火温度确定的情况下根据电泳结果不断地调整引物浓度的比例来保证每对引物的扩增效果，最终构建了2组多重PCR体系，包括A组的十二重PCR体系和B组的四重PCR体系。其中本实验A组十二重PCR体系经过不断调整最终确定了各引物浓度的比例，而B组四重PCR体系仅实验一次就获得了良好的扩增效果，说明当多重PCR体系中引物数量多的时候，引物浓度的比例尤为重要。
3.2 微卫星标记在黄鳝亲子鉴定中的应用
近年来，微卫星亲子鉴定被广泛应用于种质资源保护、新品种选育及种群遗传管理等方面，尤其在品种选育工作中能避免近亲交配并指导生产[36-38]。在微卫星亲子鉴定中，主要根据累计排除概率 (CE-PP) 判定。赵桐茂[39]研究表明，如果CE-PP<80%，则不能确定存在亲子关系；如果95%<CEP≤99%，则可能存在亲子关系；如果CE-PP≥99.73%，则可以认为存在亲子关系。本实验中，16个微卫星位点的CE-PP达到99.99%，说明这批微卫星可以用于黄鳝的亲子鉴定；16个微卫星位点的模拟鉴定率高达99.96%，实际鉴定率达95%，这一结果与哲罗鱼 (Hucho taimen)[40]和黄河鲤 (Cyprinus carpio haematoperus)[41]的结果相近。在实际生产进行亲子鉴定时，往往是在更多的候选亲本中鉴定亲本和子代的关系。本实验进行了不同亲本数的模拟分析，结果表明16个微卫星标记在亲本性别未知的150对亲本和亲本性别已知的200对亲本的情况下，均可达到95%的准确率。利用另外4个黄鳝家系对建立的微卫星多重PCR体系进行验证，结果显示仅有2尾子代聚类错误，说明本实验建立的多重PCR体系可用于黄鳝亲子鉴定。尽管随着家系数量继续增加，亲子鉴定的准确率会逐渐降低，但本实验建立的亲子鉴定技术可以满足一定规模黄鳝育种对亲本数的要求，可以为开展黄鳝新品种选育工作、选育家系规模确定提供参考。此外，根据子代聚类分析的结果，可以将不同家系的子代区分开来，实际生产中在亲本缺失的情况下，可以对混养的不同家系黄鳝进行区分，以有效防止近亲繁殖。
在亲子鉴定中出现实际鉴定率低于模拟鉴定率的情形，原因可能是基因型误读及无效等位基因的存在[42]。根据O'reilly等[43]的研究，在基因分型过程中存在2%~3%的分型错误；此外，二、三碱基重复的微卫星标记在读取基因型大小时可能会因为碱基重复单元较短而出现基因型的误读，使得基因分型的结果存在一定误差。本实验中，实际鉴定率 (95%) 略低于模拟鉴定率 (99.96%) ，筛选的16个微卫星位点中有且仅有一个位点 (Ma85) 的无效等位基因频率大于10%，这可能是造成实际鉴定率略低的原因。
3.3 群体的遗传多样性分析
群体的遗传多样性主要通过有效等位基因数、杂合度 (Ho/He) 和多态信息含量表示，遗传参数越高，多态性越高，说明群体的遗传多样性越高。通过分析发现，各个位点的PIC与NE-PP大致呈反比关系，与韩叶等[44]对大麻哈鱼 (Oncorhynchus keta) 的研究结果一致。比较各子代家系与亲本群体之间的遗传多样性，发现亲本群体的平均期望杂合度和平均PIC均高于子代各群体，说明亲本群体的遗传多样性高于各子代家系群体的遗传多样性，在各子代家系群体中，子代家系81的杂合度最高，子代家系101和子代家系77次之。比较各子代家系群体之间的多态信息含量，发现子代家系101的多态信息含量最高，子代家系81和41次之。综合杂合度和多态信息含量的比较，得出子代家系81和101具有较高的遗传多样性，说明其具有较高的育种潜力，可以作为家系选育的基础群体。11个子代家系群体观测杂合度为0.631，期望杂合度为0.618，多态信息含量为0.562，整体而言杂合度和多态信息含量均低于我国长江中下游湖北、安徽、江苏、浙江4省的野生黄鳝群体[45-46]和江西1个野生黄鳝群体，但高于天津3个野生黄鳝群体[47]，说明选育家系子代群体的遗传多样性相比野生群体有所下降。
4. 结论
本研究筛选得到了16个高度多态性的微卫星标记，建立了2组可用于黄鳝亲子鉴定的微卫星多重PCR体系，在黄鳝110尾子代及22尾亲本中达到了95%的鉴定率，研究结果有助于黄鳝的家系选育与管理，并可为其种质资源保护和良种选育提供理论依据和技术手段。
-
![]()
图 1 大鳍鳠胚胎发育各时期
注：a. 受精卵；b. 胚盘期；c. 2细胞期；d. 4细胞期；e. 8细胞期；f. 16细胞期；g. 32细胞期；h. 多细胞期；i. 囊胚早期；j. 囊胚中期；k. 囊胚晚期；l. 原肠早期；m. 原肠中期；n. 原肠晚期；o. 神经胚期；p. 胚孔封闭期；q. 体节出现期；r. 原肾出现期；s. 耳囊形成期；t. 尾芽分离期；u. 肌肉效应期；v. 心脏跳动期；w. 眼囊形成期；x. 耳石出现期；y. 第二对须原基出现期；z. 胸鳍原基出现期；aa. 眼色素出现期；ab. 出膜仔鱼；标尺=500 μm。
Figure 1. Embryonic development of M. macropterus
Note: a. Fertilized egg; b. Blastodisc stage; c. 2-cell stage; d. 4-cell stage; e. 8-cell stage; f. 16-cell stage; g. 32-cell stage; h. Multicellular stage; i. Early blastula stage; j. Middle blastula stage; k. Late blastula stage; l. Early gastrula stage; m. Middle gastrula stage; n. Late gastrula stage; o. Neurula stage; p. Blastopore closed stage; q. Appearance of myomere; r. Proto-renal stage; s. Appearance of otic capsule; t. Tail bud separation period; u. Muscular effect stage; v. Heart pulsation; w. Optic capsule stage; x. Otolith formation stage; y. The second pair of whisker primordia; z.Thoracic fin primordium; aa. Occupation of ocular pigment; ab. Hatching stage; Bar=500 μm.
![]()
图 2 大鳍鳠仔稚鱼发育
注：a. 出膜1 d；b. 出膜2 d；c. 出膜3 d；d. 出膜4 d；e. 出膜5 d；f. 出膜6 d；g. 出膜7d；h. 出膜8 d；i. 出膜9 d；j. 出膜10 d；k. 出膜 11 d；l. 出膜12 d；m. 出膜13 d；n. 出膜14 d；o. 出膜16 d；p. 出膜17 d；q. 出膜18 d；r. 出膜19 d；s. 出膜20 d；t. 出膜21 d；标尺=1 mm。
Figure 2. Development of M. macropterus larva and juvenile
Note: a. 1 d post hatching; b. 2 d post hatching; c. 3 d post hatching; d. 4 d post hatching; e. 5 d post hatching; f. 6 d post hatching; g. 7 d post hatching; h. 8 d post hatching; i. 9 d post hatching; j. 10 d post hatching; k. 11 d post hatching; l. 12 d post hatching; m. 13 d post hatching; n. 14 d post hatching; o. 16 d post hatching; P. 17 d post hatching; q. 18 d post hatching; r. 19 d post hatching; s. 20 d post hatching; t. 21 d post hatching; Bar=1 mm.
表 1 大鳍鳠胚胎发育时序
Table 1 Process of embryonic development of M. macropterus
发育时期
Developmental stage受精后时间
Time after
fertilization/
h持续时间
Duration/
h平均卵径
Average egg
diameter/
mm平均水温
Average
temperature/
℃积温
Accumulative
temperature/
(h·℃)图1
Fig. 1受精卵
Fertilized egg受精卵
Fertilized egg0 2.00 3.46±0.30 26.1 52.20 a 胚盘期
Blastodisc stage胚盘期
Blastodisc stage2.00 1.35 3.44±0.16 26.2 35.37 b 卵裂期
Cleavage stage2 细胞期
2-cell stage3.35 3.15 3.52±0.12 26.3 82.85 c 4 细胞期
4-cell staged 8 细胞期
8-cell stagee 16 细胞期
16-cell stagef 32 细胞期
32-cell stageg 多细胞期
Multicellular stageh 囊胚期
Blastula stage囊胚早期
Early blastula stage6.50 12.65 3.53±0.17 26.3 332.70 i 囊胚中期
Middle blastula stagej 囊胚晚期
Late blastula stagek 原肠期
Gastrula stage原肠早期
Early gastrula stage19.15 7.05 3.42±0.18 26.3 185.42 l 原肠中期
Middle gastrula stagem 原肠晚期
Late gastrula stagen 神经胚期
Neurula stage神经胚期
Neurula stage26.20 13.47 3.40±0.05 26.0 350.22 o 胚孔封闭期
Blastopore closed stagep 器官形成期
Organogenesis stage体节出现期
Appearance of myomere39.67 32.25 3.57±0.28 26.1 841.73 q 原肾出现期
Proto-renal stager 耳囊形成期
Appearance of otic capsules 尾芽分离期
Tail bud separation periodt 肌肉效应期
Muscular effect stageu 心脏跳动期
Heart pulsationv 眼囊形成期
Optic capsule stagew 耳石出现期
Otolith formation stagex 第二对须原基出现
The second pair of whisker primordiay 胸鳍原基出现期
Thoracic fin primordiumz 眼色素出现期
Occupation of ocular pigmentaa 出膜期
Hatching stage出膜仔鱼
Hatching stage71.90 7.08 26.1 184.79 ab 
下载: 导出CSV
表 2 大鳍鳠不同日龄仔、稚鱼生长变化情况
Table 2 Growth variation of larval and juvenile M. macropterus at different time
日龄
Days post
hatching/d全长
Total
length/mm体质量
Body
mass/g全长特定
增长率
SGRL/%体质量特定
增长率
SGRm/%1 7.31±0.52 0.011±0.005 9.79 12.31 2 10.01±0.96 0.013±0.001 3 10.19±0.42 0.016±0.001 4 10.85±1.76 0.018±0.003 5 11.58±0.22 0.018±0.003 6 11.93±1.19 0.020±0.001 7 12.96±1.80 0.021±0.001 4.15 12.25 8 13.74±3.45 0.021±0.002 9 13.92±0.45 0.024±0.003 10 14.68±1.08 0.030±0.001 11 15.81±2.16 0.033±0.002 4.30 13.17 12 17.00±3.87 0.034±0.013 13 18.26±2.95 0.037±0.011 14 18.38±3.72 0.055±0.003 15 19.19±2.60 0.056±0.013 16 21.65±2.25 0.061±0.001 17 22.65±3.79 0.062±0.007 18 23.06±2.93 0.089±0.004 19 23.11±4.08 0.100±0.011 20 24.11±0.21 0.115±0.032 21 24.32±0.05 0.122±0.019 22 25.00±0.23 0.125±0.001 2.28 12.13 23 25.32±1.23 0.127±0.007 24 26.02±0.78 0.130±0.004 25 27.28±2.55 0.140±0.011 26 27.98±0.56 0.163±0.032 27 28.36±1.78 0.189±0.019 28 29.50±2.21 0.227±0.005 29 29.72±0.66 0.227±0.019 30 30.02±1.68 0.230±0.035
下载: 导出CSV
表 3 大鳍鳠与其他鲶形目鱼类胚胎发育比较
Table 3 Comparison of embryonic developmental characteristics between M. macropterus and other Siluriformes fish
种类
Specis卵径
Egg diameter/
mm吸水后卵径
Egg diameter
after water
absorption/
mm出膜方式
Hatching
pattern出膜长度
Hatching
length/mm胚胎发育水温
Temperature/
℃胚胎发育历时
Embryonic
development文献
Reference大鳍鳠
M. macropterus2.72±0.2 3.46±0.30 尾部破膜 7.31±0.52 26 79 h 本文 黄颡鱼
Pelteobagrus fulvidraco1.86~2.26 尾部破膜 4.8~5.5 24~27.5 62 h 50 min [18] 长吻鮠
Leiocassis longirostris2.18~2.35 2.51~2.58 尾部破膜 5.33~6.62 21~27 59 h 5 min [19] 瓦氏黄颡鱼
Pelteobagrus vachelli1.67 2.04 尾部破膜 4.2~4.35 20~24.5 49 h 15 min [20] 24~28 64 h 乌苏里鮠
Pseudobagrus ussuriensis2.0 2.5 尾部破膜 7 25.5~26 49 h 26 min [21] 斑鳠
Mystus guttatus3.1~3.8 4.0~4.2 尾部破膜 27.5±0.5 55 h [22] 兰州鲶
Silurus lanzhouensis2.07~2.4 3.44~3.46 尾部破膜 4.5 21.5~22.5 53 h 10 min [23] 24~25 39 h 37 min 27~28 30 h 20 min 南方大口鲶
S. soldatovi meridionalis2.089 3.055 尾部破膜 5.0~7.05 16.5~18.5 53 h 25 min [24] 鲶
S. asotus1.4~1.7 4.05~5.7 尾部破膜 4.64~5.32 27.5~31 29.5 h [25] 胡子鲶
Claris fuscus1.7~1.9 1.9~2.1 腹部卵黄囊破膜 4.8~5.1 28.5~31 28 h 25 min [26] 革胡子鲶
C. gariepinus1.23~1.36 1.34~1.43 尾部破膜 3.12~3.62 26.4~29.2 23 h 45 min [27]
下载: 导出CSV
-
[1] 葛玲瑞, 刘金凤, 刘科均, 等. 大鳍鳠流水养殖技术要点[J]. 水产养殖, 2022, 43(11): 51-52. doi: 10.3969/j.issn.1004-2091.2022.11.013 [2] 杨德国, 周剑光, 吴国犀, 等. 长江大鳍鳠的人工繁殖[J]. 中国水产科学, 1998, 5(2): 27-31. [3] 朱成科, 郑永华, 高元彧, 等. 大鳍鳠的生物学特性和人工繁养殖前景研究[J]. 水利渔业, 2007, 27(5): 86-89. doi: 10.3969/j.issn.1003-1278.2007.05.037 [4] 许倩, 刘良国, 黄文轩, 等. 沅水五强溪水库大鳍鳠的遗传多样性分析[J]. 四川农业大学学报, 2019, 37(2): 253-258. [5] 杨德国, 吴国犀, 周剑光, 等. 大鳍鳠亲鱼池塘驯养和培育技术[J]. 淡水渔业, 1996, 26(2): 3-6. [6] 马跃岗, 朱杰, 陈元坤, 等. 大鳍鳠的人工繁殖与苗种培育研究[J]. 重庆水产, 2012(2): 20-22. [7] 李晓莉, 朱永久, 杨德国, 等. 大鳍鳠人工繁殖研究进展及前景展望[J]. 水产科技情报, 2023, 50(2): 128-132. [8] 殷名称. 鱼类早期生活史研究与其进展[J]. 水产学报, 1991, 15(4): 348-358. [9] PUVANESWARI S, MARIMUTHU K, KARUPPASAMY R, et al. Early embryonic and larval development of Indian catfish, Heteropneustes fossilis[J]. Eurasian J Biosic, 2009, 3: 84-96.
[10] 王德寿, 杨松林. 大鳍鳠的可量性状、长重关系和肥满度[J]. 西南师范大学学报(自然科学版), 1992, 17(4): 510-515. doi: 10.13718/j.cnki.xsxb.1992.04.018 [11] CHATAKONDI N G, TORRANS E L. The influence of ovarian fluid pH of stripped unfertilized channel catfish, Ictalurus punctatus, eggs on the hatching success of channel catfish [female]×blue catfish, Ictalurus furcatus[male], hybrid catfish eggs[J]. J World Aquac Soc, 2012, 43(4): 585-593. doi: 10.1111/j.1749-7345.2012.00577.x
[12] NAHIDUZZAMANA M, HASSAN M M, ROY P K, et al. Sperm cryopreservation of the Indian major carp, Labeo calbasu: effects of cryoprotectants, cooling rates and thawing rates on egg fertilization[J]. Anim Report Sci, 2012, 136(1/2): 133-138.
[13] 林红军, 刘莎, 曾玲. 黄颡鱼人工授精及胚胎发育研究[J]. 黑龙江畜牧兽医, 2019(12): 138-140, 170. doi: 10.13881/j.cnki.hljxmsy.2018.09.0099 [14] 韩军军, 贺疆滔, 陈朋, 等. 扁吻鱼胚胎发育及仔稚鱼形态与生长观察研究[J]. 南方水产科学, 2021, 17(1): 59-66. [15] 许静, 谢从新, 邵俭, 等. 雅鲁藏布江尖裸鲤胚胎和仔稚鱼发育研究[J]. 水生态学杂志, 2011, 32(2): 86-95. doi: 10.15928/j.1674-3075.2011.02.019 [16] 王亚利, 刘玥, 田佳佳, 等. 温度对四川华鳊仔、稚鱼生长发育及存活率的影响[J]. 渔业科学进展, 2021, 42(5): 47-54. doi: 10.19663/j.issn2095-9869.20200308002 [17] 钟全福, 樊海平, 薛凌展, 等. 大刺鳅仔、稚鱼的生长及摄食节律[J]. 广东海洋大学学报, 2020, 40(2): 28-34. [18] 王令玲, 仇潜如, 邹世平, 等. 黄颡鱼胚胎和胚后发育的观察研究[J]. 淡水渔业, 1989(5): 9-12. [19] 苏良栋, 何学福, 张耀光, 等. 长吻鮠 Leiocassis longirostris Gunther胚胎发育的初步观察[J]. 淡水渔业, 1985, 15(4): 2-4. [20] 袁伦强. 温度对瓦氏黄颡鱼 (Pelteobagrus vachelli)胚胎和仔鱼发育的影响[D]. 重庆: 西南师范大学, 2005: 8-9. [21] 崔宽宽, 李贺密, 苗建春, 等. 乌苏里鮠胚胎发育观察[J]. 河北渔业, 2012(2): 32-36, 44. doi: 10.3969/j.issn.1004-6755.2012.02.013 [22] 焦宗垚, 陈赛, 陈昆慈, 等. 珠江斑鳠的早期胚胎发育观察[J]. 广东海洋大学学报, 2007, 27(3): 20-23. doi: 10.3969/j.issn.1673-9159.2007.03.005 [23] 吴旭东, 李力, 张锋, 等. 兰州鲇胚胎发育初步研究[J]. 当代水产, 2011, 36(1): 65-67. doi: 10.3969/j.issn.1674-9049.2011.01.021 [24] 谢小军. 南方大口鲶的胚胎发育[J]. 西南师范大学学报 (自然科学版), 1986, 13(3): 72-78. [25] 魏刚, 罗学成. 鲶胚胎和幼鱼发育的研究[J]. 四川师范学院学报 (自然科学版), 1994, 15(4): 350-355. [26] 潘炯华, 郑文彪. 胡子鲶的胚胎和幼鱼发育的研究[J]. 水生生物学集刊, 1982, 7(4): 437-444. [27] 潘炯华, 郑文彪. 革胡子鲇的胚胎和仔、稚鱼发育的研究[J]. 华南师范大学学报 (自然科学版), 1987(1): 19-28. [28] 张耀光, 王德寿, 罗泉笙. 大鳍鳠的胚胎发育[J]. 西南师范大学学报 (自然科学版), 1991, 16(2): 223-229. [29] 牛建功, 张涛, 沙文军, 等. 斑重唇鱼胚胎及早期仔鱼发育的初步观察[J]. 河北渔业, 2021(11): 7-12. doi: 10.3969/j.issn.1004-6755.2021.11.003 [30] HANSE T K, FALK-PETERSEN I B. The influence of rearing temperature on early development and growth of spotted wolffish Anarhichas minor (Olafsen)[J]. Aquac Res, 2002, 32(5): 369-378.
[31] THORN M W, MORBEY Y E. Egg size and the adaptive capacity of early life history traits in chinook salmon (Oncorhynchus tshawytscha)[J]. Evol Appl, 2018, 11(2): 205-219. doi: 10.1111/eva.12531
[32] HAN M Y, LUO M, YANG R, et al. Impact of temperature on survival and spinal development of golden pompano Trachinotus ovatus (Linnaeus 1758)[J]. Aquac Rep, 2020, 18: 100556. doi: 10.1016/j.aqrep.2020.100556
[33] TSUJI M, ABE H, HANYUU K, et al. Effect of temperature on survival, growth and malformation of cultured larvae and juveniles of the seven-band grouper Epinephelus septemfasciatus[J]. Fish Sci, 2014, 80(1): 69-81. doi: 10.1007/s12562-013-0680-z
[34] 王志强, 庞守忠. 黄颡鱼仔稚鱼发育和摄食习性研究[J]. 江苏农业科学, 2009(6): 311-313. doi: 10.3969/j.issn.1002-1302.2009.06.131 [35] 刘明鉴, 郭华阳, 高杰, 等. 尖翅燕鱼早期胚胎发育及仔稚鱼形态观察[J]. 南方水产科学, 2022, 18(4): 103-111. [36] 乔晔. 长江鱼类早期形态发育与种类鉴别[D]. 武汉: 中国科学院研究生院 (水生生物研究所), 2005: 13-17. [37] 徐滨, 朱祥云, 魏开金, 等. 拉萨裂腹鱼的胚胎及仔稚鱼发育特征[J]. 大连海洋大学学报, 2020, 35(5): 663-670. doi: 10.16535/j.cnki.dlhyxb.2019-215 [38] KAJIWARA K, NAKAYA M, SUZUKI K, et al. Effect of egg size on the growth rate and survival of wild walleye pollock Gadus chalcogrammus larvae[J]. Fish Oceanogr, 2022, 31(3): 238-254. doi: 10.1111/fog.12575
[39] RIDEOUT R M, TRIPPEL E A, LITVAK M K. Effects of egg size, food supply and spawning time on early life historysuccess of haddock Melanogrammus aeglefinus[J]. Mar Ecol Prog Ser, 2005, 285: 169-180. doi: 10.3354/meps285169
[40] 杨建, 耿龙武, 王雨, 等. 拟赤梢鱼的胚胎发育和仔稚鱼生长特性观察[J]. 水生生物学报, 2021, 45(3): 643-651. doi: 10.7541/2021.2020.017 [41] 张呈祥, 徐钢春, 徐跑, 等. 美洲鲥仔、稚、幼鱼的形态发育与生长特征[J]. 中国水产科学, 2010, 17(6): 1227-1235. [42] 赵红雪, 吴旭东, 李力, 等. 兰州鲇仔稚鱼发育的初步研究[J]. 当代水产, 2010, 35(12): 59-62. doi: 10.3969/j.issn.1674-9049.2010.12.014 [43] ÇELIK P, CIRILK Ş. Embryonic and larval development of serpae tetra Hyphessobrycon eques (Steindachner, 1882)[J]. Aquac Res, 2020, 51(1): 292-306. doi: 10.1111/are.14375
[44] 李晓莉, 朱永久, 杨德国, 等. 大规格大鳍鳠对隐蔽所的选择及其行为特征[J]. 中国水产科学, 2022, 29(12): 1768-1777. doi: 10.12264/JFSC2022-0270 [45] 张耀光, 何学福. 长吻鮠幼鱼发育的研究[J]. 水生生物学报, 1991, 15(2): 153-160. [46] 李萍, 庹云. 岩原鲤早期行为习性的初步观察[J]. 安徽农业科学, 2008, 36(2): 565-566. doi: 10.3969/j.issn.0517-6611.2008.02.078 [47] 黄小林, 杨育凯, 李涛, 等. 黄斑篮子鱼仔、稚、幼鱼形态观察与生长研究[J]. 南方水产科学, 2018, 14(5): 88-94. doi: 10.3969/j.issn.2095-0780.2018.05.011 [48] 施兆鸿, 彭士明, 罗海忠, 等. 人工育苗条件下舟山牙鲆仔稚鱼的摄食与生长[J]. 上海水产大学学报, 2007, 16(6): 566-557. [49] 王亚利, 刘玥, 田佳佳, 等. 四川华鳊仔稚鱼生长与形态发育研究[J]. 渔业科学进展, 2020, 41(4): 58-67. doi: 10.19663/j.issn2095-9869.20190627001 [50] 王雨福, 肖志忠, 刘清华, 等. 斑石鲷早期发育的异速生长模式[J]. 海洋科学, 2016, 40(5): 43-48. doi: 10.11759//hykx20141216005 -
期刊类型引用(1)
1. 韩豪祥，廉杰，金洪宇，马波. 额尔齐斯河黑鲫胚胎及早期仔鱼发育特征. 中国水产科学. 2024(06): 684-692 . 
百度学术
其他类型引用(0)
计量
- 文章访问数: 206
- HTML全文浏览量: 48
- PDF下载量: 41
- 被引次数: 1
粤公网安备 44010502001741号
