Analysis of reproductive biological characteristics of Silurus asotus
-
摘要:
对沅江干流水域鲇 (Silurus asotus) 的生长与繁殖特征进行研究，为其种质资源的保护及合理利用提供科学依据。2018年5月至2020年4月在湖南省常德市沅江凌津滩库区采集到鲇738 尾，对其体长体质量关系及繁殖特性进行分析。结果显示，鲇体长和体质量呈幂函数关系，关系式为W=0.056L2.57(R=0.97)，幂指数b值为2.57，与理论匀速生长值3之间差异显著 (P<0.05) ，群体生长为异速型，群体雌雄性比为1.21∶1。鲇群体肥满度平均为0.88，雌性平均0.94，雄性平均0.81，雌性极显著高于雄性 (P<0.01)。性成熟系数雌性极显著高于雄性 (P<0.01)，雌性平均7.51，雄性平均1.85，均为5月最高、9月最低，雌、雄性成熟系数的周年变化趋势基本同步。卵径5月最高，达1 626.5 μm；5—8月逐步下降，其变化与实际产卵高峰期一致。个体平均绝对怀卵量为 (29 677±18 174) 粒，个体平均相对怀卵量为 (107±65) 粒·g−1。绝对怀卵量与体质量和体长均成直线相关，关系式分别为：F=437.36W+945.09 (R2=0.93, P<0.01) 和F=1 205.38L−21 285.27 (R2=0.87, P<0.01)，从性成熟系数的周年变化推测，鲇的产卵高峰期为5月中旬至8月上旬。
Abstract:To study the growth and reproductive characteristics of catfish (Silurus asotus) in the main stream of the Yuanjiang River, and to provide theoretical support for its protection and utilization, we had collected 738 individuals in the Yuanjiang Lingjintan Reservoir area of Changde City, Hunan Province, from May 2018 to April 2020. Then we analyzed the relationship between body length/body mass and reproductive characteristics. The results show a power function relationship between body length and body mass for Yuanjiang S. asotus, with the relationship of W=0.056L2.57 (R=0.97) and a significant difference between b-value and theoretical uniform growth value of 3 (P<0.05). The population growth was heteroscedastic, with a population sex ratio of 1.21:1. The average fertility of Yuanjiang S. asotus population was 0.88 (Female: 0.94, male: 0.81, female significantly higher than male, P<0.01). The sexual maturity coefficients of females were higher than those of males (P<0.01) (Female: averagely 7.51, male: averagely 1.85), both highest in May and lowest in September. The annual trends in sexual maturity coefficients for females and males were basically synchronized. The egg diameter reached a maximum of 1626.5 μm in May, gradually decreasing from May to August, which was consistent with the actual peak spawning period. The mean absolute egg mass was 29 677±18 174 and relative egg mass was (107±65) per·g−1. The absolute egg mass was linearly correlated with body mass and length as F=437.36W+945.09 (R2=0.93, P<0.01) and F=1 205.38L−21 285.27 (R2=0.87, P<0.01), respectively. The peak spawning period for Yuanjiang S. asotus was inferred from the annual variation in sexual maturity coefficients to be from mid-May to early August.
-
疣吻沙蚕 (Tylorrhynchus heterochaetus) 又名疵吻沙蚕，属环节动物门、多毛纲、叶须虫目、沙蚕科，俗称禾虫、流蜞，在中国东南沿海河口地区的泥沙质浅滩或稻田中广泛分布[1-2]。作为一种经济价值极高的多毛类，其味道鲜美且营养丰富，是中国广东、广西、福建、香港、澳门及东南亚各地独具特色的水产品，素有“水中冬虫夏草”之美誉[2-3]。更重要的是，疣吻沙蚕与水稻具有天然的共生关系，发展疣吻沙蚕稻田综合种养模式可显著增加水稻种植的经济效益[4]。“水稻+疣吻沙蚕”生态综合种养技术正日趋成熟，自2021年起连续3年入选广东省农业主推技术，在助力乡村振兴、保障国家粮食安全等方面应用潜力巨大。然而，疣吻沙蚕具有生境狭窄、种群地理隔离、群体恢复力低等物种特性，受栖息地破坏、过度捕捞、环境污染等因素的影响其自然种群呈逐年衰退的趋势，在有些地区甚至已经绝迹[4-5]。在突破疣吻沙蚕全人工繁殖技术的基础上，近年已实现其苗种的规模化培育，增养殖业的发展以及技术的革新培育了大量的人工繁育群体，但因遗传育种研究的滞后，疣吻沙蚕种质开始出现生长性能及抗逆性下降等退化现象[4]。为建立疣吻沙蚕资源管理策略并促进其可持续开发利用，特别是避免常年的产业化增养殖引起遗传结构单一化，防止人工繁育群体污染自然种群，有必要采用分子手段科学地分析及评估其种质现状，揭示种群遗传多样性水平与遗传结构，并以此指导新品种的培育与遗传改良研究。
目前，国内外针对疣吻沙蚕的研究主要集中在形态学、生活史、繁殖生物学、增养殖技术、营养与活性成分分析等方面[3,5]。Chen等[5-6]完成了疣吻沙蚕线粒体基因组测序，并基于线粒体COI序列分析了7个地理群体的遗传结构，除此之外，遗传学相关的研究鲜有报道；分子遗传信息的匮乏制约了疣吻沙蚕资源保护利用研究的深入开展。微卫星标记具有多态性丰富、共显性遗传、基因组覆盖广泛等优点[7]，该技术操作简便、检测周期短、稳定性高、重复性好，已广泛应用于谱系鉴定、亲缘关系检测、遗传多样性分析、种质资源评价、遗传图谱构建及基因定位等领域[8]。基因组survey是利用高通量测序技术实施小片段、低深度测序，然后基于Lander-Waterman模型进行k-mer分析，并根据k-mer频率和深度的统计结果评估物种基因组的大小与复杂程度等关键特征[9]。基因组survey不仅能为全基因组测序及高质量组装提供科学依据，而且对无参考基因组序列的物种而言，利用基因组survey数据大规模开发微卫星标记是当前最高效的策略之一，较传统方法具有效率高、周期短、成本低的优势[9]。基因组微卫星鉴定研究在水产经济动物中已广泛开展，筛选的标记获得了良好的遗传分析效果。例如，张永德等[9]在卵形鲳鲹 (Trachinotus ovatus) 基因组survey数据中检测到190 121个微卫星位点并成功开发了多态分子标记，为全基因组测序与组装、渔业资源保护利用及良种选育奠定了基础；上官清等[10]分析了斑鳢 (Channa maculata) 基因组微卫星特征，并筛选到20个多态性位点用于群体遗传多样性和遗传结构分析，为该物种的遗传监测、亲缘关系鉴定、种质资源养护及管理提供了技术支持。
本研究对疣吻沙蚕基因组进行低深度高通量测序，通过k-mer分析预测基因组大小、杂合度和重复比例等信息，对测序数据初步组装后搜索组装序列中的微卫星位点，分析其特征与分布规律，初步验证标记的有效性和多态性，以期指导基因组精细图谱的绘制，并为群体遗传学研究提供可靠的标记资源。
1. 材料与方法
1.1 实验材料
疣吻沙蚕采集自阳江市广东阳海农业技术发展有限公司疣吻沙蚕增养殖试验基地(111°55'23"E、 21°49'16"N) 保种的越南海防群体。用于基因组测序的疣吻沙蚕幼体用适量无菌水洗涤3次，解剖后剪取体壁肌肉组织装入2 mL冻存管，液氮速冻后置于−80 ℃保存。微卫星标记多态性验证群体的30尾个体经无菌水洗涤后，解剖剪取体壁肌肉组织，置于体积分数为95%的乙醇中在−20 ℃下保存。
1.2 基因组测序及生物信息学分析
运用苯酚-氯仿法提取疣吻沙蚕基因组DNA，经1.0% (w) 琼脂糖凝胶电泳检测完整度后，用NanoDrop 2000超微量分光光度计 (ThermoFisher，美国) 检测浓度和纯度。利用Covaris超声波破碎仪将基因组DNA片段化，筛选合适长度的DNA片段，经末端修复、加A尾、添加测序接头、纯化及PCR扩增等步骤建立350 bp小片段文库。测序文库经Qubit 2.0 (ThermoFisher，美国) 和Agilent 2100 Bioanalyzer (Agilent，美国) 检测浓度和插入片段大小后，用Illumina HiseqTM X Ten平台进行双末端测序。基因组测序工作委托北京诺禾致源科技股份有限公司完成。原始数据经质控和过滤后，用GCE 1.0.0[11]软件对有效数据进行k-mer分析。运用SOAPdenovo 2.01[12]软件，选择k-mer=41将有效数据组装至contig和scaffold级别，并统计GC含量和覆盖深度等信息。
1.3 微卫星位点搜索与引物设计
使用MISA (Microsatellite identification tool) 软件(http://pgrc.ipk-gatersleben.de/misa/misa.html)在长度大于500 bp的序列中搜索微卫星位点。运行参数为：重复基序长度1~6 bp，单碱基重复次数≥12次，二碱基重复次数≥6次，三碱基、四碱基重复次数≥5次，五碱基、六碱基重复次数≥4次；将间隔区域长度小于100 bp的相邻微卫星归为1个复合型位点。基于微卫星侧翼序列，用Primer 5[13]批量设计引物，主要参数为：引物长度18~27 bp，扩增产物100~300 bp，退火温度(Tm) 55~65 ℃，GC含量40%~60%，正、反向引物退火温差≤5 ℃；尽量避免出现发卡结构、二聚体、错配和引物二聚体；每个微卫星位点生成3~5对候选引物。随机选取50对微卫星引物委托生工生物工程(上海)股份有限公司合成。
1.4 多态微卫星标记的验证与筛选
使用天根生化科技(北京)有限公司的海洋动物组织基因组DNA提取试剂盒从疣吻沙蚕体壁肌肉组织中提取基因组DNA。PCR反应体系为20 μL，正、反向引物 (10 μmol·L−1) 各0.5 μL，DNA模板 (15 ng·μL−1) 1.0 μL，2×PCR Mix 10 μL，用超纯水补至20 μL。扩增反应由Bio-Rad My Cycler Thermal Cycler (Bio-Rad，美国) 完成，运行程序为：95 ℃预变性5 min；95 ℃变性30 s，59~60 ℃退火30 s，72 ℃延伸30 s，25个循环；72 ℃延伸5 min。PCR产物先进行1.5% (w) 琼脂糖凝胶电泳分析，参照预期片段大小筛选出有效扩增引物。有效扩增引物加荧光接头后 (正向引物5'端添加FAM荧光素)，在30尾疣吻沙蚕个体的基因组DNA中扩增以验证其多态性。PCR产物送至上海翼禾应用生物技术有限公司用3730XL测序分析仪 (Applied Biosystems，美国) 进行毛细管电泳分析，使用GeneMapper 3.2软件 (Applied Biosystems，美国) 进行基因分型。
1.5 数据分析
使用Excel 2016软件完成微卫星位点分布特征信息的统计分析和图表绘制。微卫星发生频率=含微卫星的序列总数/序列总数×100%；微卫星出现频率=微卫星总数/序列总数×100%；微卫星丰度 (个·Mb−1)=微卫星总数/序列总长度[14]。用GenAlEx 6.5软件[15]计算等位基因数 (Na)、有效等位基因数 (Ne)、观测杂合度 (Ho)、期望杂合度 (He)、多态信息含量 (PIC)。用Genepop在线软件 (https://genepop.curtin.edu.au/) 检验位点间的连锁不平衡及群体的哈迪-温伯格平衡 (Hardy-Weinberg Equilibrium, HWE)，并用Bonferroni法对显著性阈值进行校正。
2. 结果
2.1 疣吻沙蚕的基因组survey测序
低深度高通量测序产生的原始数据经质控后共获得57.48 Gb有效数据，碱基错误率为0.05%，Q20和Q30分别为95.56%和89.70%，表明基因组测序质量较高 (表1)。对有效数据进行k-mer分析 (k=17)，结果显示在深度为57时出现主峰值 (图1)，总k-mer为44 257 233 158，排除错误k-mer的误差影响后得到修正的基因组大小为759.53 Mb，杂合率为1.41%，重复序列比例为45.92%。选择k-mer=41将有效读段初步组装至contig和scaffold水平，最终获得的contig总长度为821 637 022 bp，最大长度为84 713 bp，N50为548 bp；scaffold总长度为840 375 821 bp，最大长度为89 326 bp，N50为662 bp。
表 1 疣吻沙蚕基因组 survey 测序数据统计Table 1. Statistics of genomic survey sequencing data of T. heterochaetus测序文库
Sequencing library原始数据量
Raw base/Gb有效数据比
Effective rate/%有效数据量
Clean base/Gb碱基错误率
Error rate/%Q20/% Q30/% GC 含量
GC content/%L1 33.12 99.68 33.01 0.04 96.23 90.77 39.13 L2 24.53 99.75 24.47 0.05 94.88 88.63 39.00 总计 Total 57.65 — 57.48 — — — — 均值 Mean — 99.72 — 0.05 95.56 89.70 39.07 ![]()
图 1 疣吻沙蚕基因组 k-mer 种类频率分布Figure 1. Frequency distribution of k-mer species in genome of T. heterochaetus2.2 疣吻沙蚕基因组微卫星位点特征
用MISA软件在109 881条组装序列中检测到 130 216个微卫星位点，长度共计2 341 179 bp；微卫星发生频率5.04%，出现频率5.97%，分布丰度为154.9个·Mb−1。疣吻沙蚕微卫星位点以单碱基和二碱基重复最为丰富，分别有45 582和42 298条，各占35.00%和32.48%；其次是三碱基重复 (18 782条)，占14.42%；六碱基数量最少，仅占2.44% (图2)。微卫星序列的重复数范围为4~56拷贝，主要集中在4~18拷贝，重复19次及以上的有3 703条，仅占3.50%；单碱基重复以12~16次最常见 (40 706条，89.30%)，二碱基重复集中在6~18次 (42 035条，99.38%)，三碱基重复以5~12次为主 (18 526条，98.64%)，四碱基重复以5~10次为主 (12 173条，98.91%)，而五碱基和六碱基重复主要为4~8次 (共11 176条，99.37%) (图2—图3)。
![]()
图 2 疣吻沙蚕基因组 6 种类型微卫星的数量与比例Figure 2. Number and proportion of six motif types of microsatellite loci in genome of T. heterochaetus![]()
图 3 疣吻沙蚕基因组微卫星重复数分布特征Figure 3. Distribution pattern of microsatellite repeat number in genome of T. heterochaetus疣吻沙蚕基因组微卫星共包含320种重复基序类型：单碱基2种、二碱基4种、三碱基10种、四碱基31种、五碱基91种、六碱基182种。单碱基重复基序拷贝数集中在12~15次，以C/G为主，占比58.02%；二碱基重复拷贝数多为6~10次，以AT/AT最为丰富，占比62.38%，其次是AC/GT (27.37%)、AG/CT (10.22%)，而CG/CG数量稀少 (0.02%)；在三碱基重复中，拷贝数多为5~12次，占比最高的是AAT/ATT，有6253条 (33.29%)，其次是ATC/GAT (32.15%)，而CCG/CGG数量最少，仅占0.06%；四碱基拷贝数多为5~8次，优势重复基序是AAAT/ATTT，有3567条 (28.98%)，其次是AATC/GATT (14.08%)、ACTC/GAGT (9.60%)、ATCC/GGAT (7.38%)；五碱基、六碱基拷贝数集中在4~7次，优势重复基序分别是AAAAT/ATTTT (17.07%) 和AACCCT/AGGGTT (11.73%) (图4)。
![]()
图 4 疣吻沙蚕基因组微卫星重复基序类型分布特征Figure 4. Distribution pattern of microsatellite motif types in genome of T. heterochaetus2.3 疣吻沙蚕基因组微卫星位点长度特征
根据重复序列长度可将微卫星位点分为两类：一类是长度达20 bp及以上的高度多态I型，另一类则是长度介于12~19 bp的中度多态II型[16]。疣吻沙蚕基因组微卫星位点的长度分布区间为12~336 bp，I型微卫星位点有43 240条，占33.21%，剩余的II型位点占66.79%；绝大部分I型微卫星位点的长度为20~39 bp，超过50 bp的位点仅占0.51% (图5-a)。进一步分析发现，在I型微卫星位点中，二、三、四、五碱基是最主要的重复类型，占85.92%，在后续多态标记的筛选中有较高的开发价值(图5-b)。
![]()
图 5 疣吻沙蚕基因组微卫星长度分布特征注：a. 不同长度区间微卫星数量及比例；b. 不同类型微卫星长度分布特征。Figure 5. Distribution pattern of length of microsatellite loci genome of T. heterochaetusNote: a. Number and percentage of microsatellite loci at different length intervals; b. Length distribution of the six motif types of microsatellite loci.2.4 疣吻沙蚕基因组微卫星多态标记验证筛选
使用Primer 5成功对37 370个微卫星位点设计了引物。随机选取50个位点合成引物并进行PCR验证，共获得41对 (82%) 有效扩增引物，表明MISA软件鉴定的微卫星位点具有较高有效性。多态性筛选结果显示，15个 (30%) 微卫星位点的引物表现出稳定且可重复的多态性 (表2，图6)。在30尾疣吻沙蚕中，15对引物共检测到87个等位基因，平均等位基因数5.800，等位基因频率为0.060~0.400，其中ThGM021位点的等位基因数最少 (Na=2.000)，ThGM004位点的等位基因数最多 (Na=12.000)。Ne为1.164~6.713，平均值为3.328；Ho为0.050~0.879，平均值为0.487；He为0.141~0.789，平均值为0.561；PIC为0.136~0.776，平均值为0.511 (表3)。在15个位点中，8个属高度多态性位点 (PIC>0.5)，5个属中度多态性位点 (0.25<PIC≤0.5)，2个属低度多态性位点 (PIC≤0.25)，表明该群体的遗传多样性较丰富。经Bonferroni校正后，有3个位点 (ThGM006、ThGM011、ThGM040) 偏离HWE，其中ThGM040属于高度多态性位点。各位点间无连锁不平衡现象。
表 2 疣吻沙蚕 15 对多态微卫星引物信息Table 2. Information of 15 polymorphic microsatellite loci in genome of T. heterochaetus位点
Locus引物序列 (5'—3')
Primer sequence (5'–3')重复单元
Repeat unit产物大小
Size/bp退火温度
Annealing temperature/℃ThGM004 F: TGCTGCTACTGCTACAGCTACTATG (TAC)18 289 60.0 R: CTGACAAAGTTTGGTGGCTG ThGM006 F: TGAAAATTAGTGTGATTTTGTCCC (CA)11 260 59.0 R: AGCCAACCAGAACATGAACA ThGM011 F: AACTTGGACTAAGGCTATCAAAAA (AG)17 220 59.0 R: CTTGGGGTTCATGCATCATT ThGM015 F: TTGGTTGTTATCCATGCACC (TAT)12 279 59.5 R: AGACAGCAGTGAAATAGCACCA ThGM017 F: ATTCGATAAGCATTCCACCG (ATGG)8 215 60.0 R: CTTGGTAGCTGGCCTGTCTC ThGM021 F: TGCGAAATGAGAAGTGAGCA (TA)10 277 60.0 R: TGCCTGTGTGGAATACCAAG ThGM024 F: ACCTGTCCACCCGTCATTTA (TAT)14 294 59.5 R: CCTTTAGGGGATGGCTACAA ThGM029 F: GAGCAAAATATTCAAGTTGGCA (ATT)12 243 59.0 R: TTGTTTGTCATATCTTCTAAAGAGCA ThGM033 F: GGAGTGGGGAGGATTTTAGC (TG)18 277 60.0 R: CCATGTACAGCATTCAGCCA ThGM035 F: GTAAGGGCAAGGGTTGTGAA (AG)13 226 60.0 R: ACCGTTACCCTAACCCCAAC ThGM038 F: TTACCCTGCCATCCTACCAG (TG)20 157 60.0 R: CTATTCTGCCAGTGGTCGCT ThGM040 F: GGATCCAGAAGGGGTAAAGC (TTA)11 239 59.5 R: GTTGGTCATGTTCCTGTTGC ThGM041 F: ACCAGCTGCTAGAGGCAGAC (ATG)7 260 60.0 R: TTAGGTCCTCACCCAGGGAT ThGM043 F: AAAAGCAAGTGGTAACACAAAATG (TCAT)11 272 59.5 R: CATTGGGCTCTGGGAATAAA ThGM047 F: CGACCTGCGGATTTAATTTG (TGG)12 148 60.0 R: ATATCTTGGCGGCGGATAG 注：F. 正向引物；R. 反向引物。 Note: F. Forward primer; R. Reverse primer. ![]()
图 6 部分多态微卫星标记的毛细管电泳分型结果Figure 6. A set of polymorphic microsatellite loci visualized by high-resolution capillary electrophoresis表 3 15 个多态微卫星位点在疣吻沙蚕群体中的遗传特征Table 3. Genetic characteristics of 15 polymorphic microsatellite loci in a T. heterochaetus population位点 Locus 等位基因数 Na 有效等位基因数 Ne 观测杂合度 Ho 期望杂合度 He 多态信息含量 PIC 哈迪-温伯格平衡的P 值 PHWE ThGM004 12 6.672 0.697 0.775 0.726 0.275 ThGM006 4 1.642 0.367 0.388 0.372 0.001* ThGM011 5 1.608 0.257 0.349 0.377 0.026* ThGM015 8 4.933 0.438 0.789 0.776 0.225 ThGM017 3 1.521 0.066 0.271 0.245 0.148 ThGM021 2 1.593 0.367 0.508 0.375 1.000 ThGM024 11 6.713 0.879 0.742 0.682 0.541 ThGM029 10 5.647 0.697 0.658 0.599 0.140 ThGM033 5 3.102 0.167 0.772 0.720 0.069 ThGM035 3 2.164 0.050 0.141 0.136 0.086 ThGM038 8 4.878 0.576 0.545 0.489 0.221 ThGM040 4 2.441 0.733 0.718 0.652 0.008* ThGM041 3 2.155 0.417 0.431 0.336 0.503 ThGM043 5 2.727 0.724 0.682 0.617 0.148 ThGM047 4 2.224 0.867 0.642 0.569 0.267 均值 Mean 5.800 3.328 0.487 0.561 0.511 — 注：*. Bonferroni法校正后显著偏离哈迪-温伯格平衡(P<0.05)；n=30。 Note: *. Significant departure from Hardy-Weinberg equilibrium after Bonferroni's correction (P<0.05); n=30. 3. 讨论
3.1 疣吻沙蚕基因组基本特征
不同物种间基因组的大小和复杂程度有明显差异，会直接影响到测序策略的选择及基因组的组装效果。因此，进行全基因组测序前须先评估物种基因组的基本特征。疣吻沙蚕基因组survey测序及k-mer分析估计其基因组为759.53 Mb，远大于水蛭 (Helobdella robusta)[17]和宽体金线蛭 (Whitmania Pigra)[18]等所有已报道的蛭纲物种，也较多毛纲的海蠕虫 (Capitella teleta)[17]、欧文虫 (Owenia fusiformis)[19]、巨型管虫 (Riftia pachyptila)[20]、Lamellibrachia luymesi[21]的大；与寡毛纲通俗腔蚓 (Metaphire vulgaris)[22]及多毛纲搓稚虫 (Streblospio benedicti)[23]的大小 (约0.7 Gb) 相当；但明显小于安德爱胜蚓 (Eisenia andrei)[24]、赤子爱胜蚓 (E. fetida)[25]，以及多毛纲的旋鳃虫(Spirobranchus lamarcki)[26]和深海管虫 (Paraescarpia echinospica)[27] (1.0~1.3 Gb)。物种间基因组大小的差异性在一定程度上与重复序列比例的高低有关。疣吻沙蚕基因组的重复序列比例 (45.92%) 低于旋鳃虫[26]和深海管虫[27]，但显著高于海蠕虫[17]、巨型管虫[20]，以及所有已公布的蛭纲动物 (重复比例均低于34%)[15-18]；其基因组杂合率 (1.41%) 与安德爱胜蚓[24]和赤子爱胜蚓[25]相当，远远高于L. luymesi (0.60%)[21]、搓稚虫 (0.29%)[23]与深海管虫 (0.63%)[27]，是目前已报道的杂合率最高的多毛类。此外，组装的contig和scalffold总长分别为821 637 022和840 375 821 bp，N50分别为548和662 bp，提示过高的杂合率会造成组装序列长度大于预估基因组大小，并导致显著偏低的N50指标[28]。综上所述，疣吻沙蚕基因组属复杂基因组类型，绘制染色体级别高质量全基因组图谱应优先考虑“PacBio+Illumina+Hi-C”策略。该研究数据为后续全基因组测序与组装提供了基础资料。
3.2 疣吻沙蚕基因组微卫星分布特征
水生动物基因组中单碱基重复微卫星占优势的现象较少，已报道的有中华绒螯蟹 (Eriocheir sinensis)[29]、鲤 (Cyprinus carpio)[30]、脊尾白虾 (Exopalaemon carinicauda)[31]和胡鲶 (Clarias batrachus)[32]。在疣吻沙蚕基因组中检测到 130 216个微卫星位点，位点数量随重复基序碱基数的增加而迅速减少，单碱基重复微卫星最丰富 (35.00%)，其次为二碱基重复 (32.48%)、三碱基重复 (14.42%)，这与以二碱基或三碱基为优势类型的鱼[9,33-35]、虾[36-37]、贝类[38]明显不同，特别是与亲缘关系较近的蛭类、寡毛类相比也表现出较大差异。宽体金线蛭[39]、天锡杜拉蚓 (Drawida gisti)[40]基因组中三碱基重复占绝对优势，主要基序类型分别是AAT、ATA和ATT、AAT。王斌等[41]对4种蛭类的转录组数据进行分析，同样发现三碱基重复微卫星占优势，在宽体金线蛭中甚至高达68.00%。作为最特殊的微卫星重复类型，三碱基重复微卫星可形成复杂的环-折叠构型来稳定DNA结构，从而更有利于转录过程中的解旋和蛋白质识别[33]。此外，疣吻沙蚕转录组中二碱基AT/TA重复基序频率 (32.19%) 最高，单碱基A/T重复类型(21.17%)其次，之后是三碱基AGC/GCT (29.06%) 基序 (未发表数据)，表明微卫星类型特征在转录组和基因组水平上存在差异，类似的差异性在圆鳍鱼 (Cyclopterus lumpus)[42]、凡纳滨对虾[36,43]、宽体金线蛭[39,41]中也被证实。不同物种表现出不同的微卫星优势类型分布规律，可能与其进化水平有关。
疣吻沙蚕基因组微卫星核心区集中在4~18拷贝，随着重复次数的增加，微卫星数量呈显著降低趋势，这可能与重复单元长度的不断增加使得微卫星稳定性降低或基序高频次重复导致更高的突变率有关[44]。进一步分析发现，二碱基至六碱基微卫星重复单元中A/T含量显著高于G/C含量，即微卫星序列表现出明显的A/T碱基优势，这与许多水生动物如脊尾白虾[31]、凡纳滨对虾[36]、虾夷扇贝 (Mizuhopecten yessoensis)[38]、宽体金线蛭[39]及旧金山湾卤虫 (Artemia franciscana)[45]中的研究结果一致。早期观点认为，微卫星富含A/T的原因可能是由于CpG甲基化后的C易脱氨基转变为T，导致G/C比例不断缩小，而突变的A/T碱基类型相应增多[33]；后来却发现这也可能与微卫星位点的产生方式，即与DNA的复制滑动存在一定关系[29,46]。微卫星富含A/T的原因可能是A/T含量高则Tm值降低，其序列容易发生DNA解链，并通过复制滑动机制和重组机制产生高A/T含量的重复类型的概率更高[38]。
3.3 疣吻沙蚕微卫星标记开发及其遗传分析效果
微卫星是重要的分子遗传学研究工具，开发多态标记是其广泛应用的关键。与传统方法相比，基于高通量测序技术开发微卫星标记更具优势，产生的基因组或转录组数据是标记开发的重要资源。然而与转录组微卫星相比，基因组微卫星的多态性往往更高且分布更广泛，可获得更优的基因组覆盖率[42]。Temnykh等[16]认为长度≥20 bp的微卫星位点多态性较高，长度10~20 bp的为中等多态性，小于10 bp的多态性低。在疣吻沙蚕基因组中共筛选到43 240条序列长度≥20 bp的I型微卫星位点，占比超过33%，远高于其转录组序列中20%的I型位点比例 (未发表数据)。挑选50个微卫星位点进行有效性及多态性验证，有41个 (82%) 位点的引物可扩增出特异性条带，其中15个 (30%) 表现出稳定且可重复的多态性。疣吻沙蚕基因组多态微卫星标记的筛选成功率与卵形鲳鲹[9]、凡纳滨对虾[36]、宽体金线蛭[39]中的结果相当。由此可见，疣吻沙蚕微卫星标记筛选成功率较高，获得的微卫星位点具有良好的开发潜力，是丰富且可靠的标记资源。
用微卫星标记开展群体遗传分析时，一般认为有效等位基因数越接近观测等位基因数，群体等位基因分布就越均匀；然而在实际分析中，通常会将全部条带视为有效等位基因，无效等位基因过剩便会造成等位基因分布不均[47]。从疣吻沙蚕基因组筛选的15个多态微卫星位点中，仅3个位点 (ThGM021、ThGM035、ThGM041) 的等位基因数接近有效等位基因数，提示大部分位的等位基因分布不均，这可能是因为验证群体的样本容量偏小引起了主效等位基因的缺失[9]。有研究表明，群体的最小样本容量一般与遗传分析选用的参数有关。利用He、PIC以及香农指数等评估遗传多样性时，群体样本容量达到27便能使分析结果接近总体水平的95%；但选择Na时，样本容量则需达到52以上[48]。疣吻沙蚕的群体样本量为30，故而选择He、PIC可更好地评价微卫星位点的多态性。15个位点的平均Ho (0.487)、平均He (0.561)和平均PIC (0.511)共同表明，疣吻沙蚕验证群体具有较丰富的遗传多样性，这与Chen等[5]利用COI标记的分析结果并不完全一致，可能是因为微卫星属于核标记，多态性更高且受选择性作用更小，其揭示遗传变异的灵敏度高于线粒体标记。通常认为，PIC越接近1，则群体杂合个体的比例越大、多态性越高：PIC<0.25为低度多态性，0.25≤PIC<0.5为中度多态性，PIC≥0.5为高度多态性[47]。在15个多态微卫星位点中，仅2个为低度多态性，其余13个 (86.7%) 均为高度或中度多态性。因此，基于疣吻沙蚕基因组序列筛选的微卫星标记多态性高且遗传信息丰富，可为群体遗传分析、种质资源评价、分子育种研究提供优质工具。
4. 结论
疣吻沙蚕基因组属于高杂合、高重复的复杂基因组，其微卫星位点的类型丰富且具备较好的多态性潜能，可作为有效资源用于微卫星标记的大规模开发，对种质资源评价与保护利用、种群遗传学以及分子育种研究具有实际价值。
-
![]()
图 3 鲇卵巢发育分期的月份变化
Figure 3. Variation of ovaries developmental stages of S. asotus from Yuanjiang River in different months
![]()
图 4 鲇雌性性成熟比例与体长关系
Figure 4. Relationship between body length and sex mature frequency of female S. asotus in Yuanjiang River
![]()
图 5 鲇 IV 期卵巢的卵径分布
Figure 5. Oocyte diameter-frequency distribution of S. asotus ovary at stage IV
表 1 鲇 1—12 月的性比
Table 1 Sex ratio of S. asotus from January to December
月份
Month雌 Female/尾 雄 Male/尾 雌∶雄 Female∶Male 1 35 31 1.13c 2 23 39 0.59b 3 40 20 2.00e 4 38 22 1.73d 5 23 35 0.66b 6 16 49 0.33a 7 44 16 2.75f 8 46 18 2.56f 9 33 27 1.22c 10 40 16 2.50f 11 25 43 0.58b 12 41 18 2.27e 注：同列不同上标字母示差异显著(P<0.05)。表2、表5同此。 Note: Different superscripts within the same column indicate significant differences (P<0.05). The same case in Table 2 and Table 5. 
下载: 导出CSV
表 2 鲇体长、体高、体质量的周年变化
Table 2 Annual changes in body length, body height and body mass of S. asotus
月份
Month体长
Body length/cm体高
Body height/cm体质量
Body mass/g空壳质量
Eviscerated mass/g范围
Range均值
Average范围
Range均值
Average范围
Range均值
Average范围
Range均值
Average1 23.50~38.50 30.02±4.45b 3.20~5.30 4.17±0.67a 92.00~499.00 229.21±96.08b 87.00~446.00 205.30±86.67b 2 26.40~39.20 32.29±3.45b 3.30~6.15 4.69±0.52b 128.00~600.00 295.17±97.14c 121.00~529.50 272.00±85.28c 3 14.80~43.20 34.38±3.74b 3.32~6.23 4.73±0.65b 55.00~691.00 306.97±116.53c 46.00~610.50 277.97±107.46c 4 25.40~44.00 35.62±3.74b 3.46~7.25 5.03±0.73c 135.00~710.00 375.46±98.72e 119.50~628.50 337.23±88.92e 5 14.80~39.30 26.71±5.62a 2.80~5.30 4.09±0.66a 55.00~583.00 249.63±124.11b 46.50~449.50 210.10±105.23b 6 14.50~35.58 30.58±4.32b 3.40~5.60 4.46±0.49b 154.00~524.50 327.00±97.73d 142.00~473.50 290.84±84.29d 7 27.10~39.40 33.14±3.01b 3.80~5.60 4.63±0.56b 168.00~449.50 288.70±88.70c 112.00~342.50 254.03±76.60c 8 30.70~40.11 33.95±2.91b 3.60~5.75 4.77±0.56b 209.00~644.00 340.78±105.34e 195.50~564.50 300.72±90.08d 9 21.52~36.70 26.77±4.09a 1.27~6.30 4.15±0.63a 91.50~583.00 177.59±101.51a 81.50~549.00 158.14±94.64a 10 25.40~41.20 32.57±3.95b 3.47~8.60 4.93±1.21c 152.50~772.00 327.50±160.20d 174.00~661.50 296.07±139.70d 11 21.50~49.73 39.52±4.81c 3.70~7.50 5.88±0.76d 184.00~1042.00 381.60±158.40e 124.00~917.50 341.43±137.50e 12 26.80~47.32 35.07±5.34b 3.35~6.38 4.92±0.82c 149.50~925.00 370.90±182.30e 134.00~872.50 338.06±162.20e
下载: 导出CSV
表 3 鲇性腺发育状况
Table 3 Gonadal development of S. asotus
月份
Month雄 Male 雌 Female 发育分期
Developmental stage颜色
Color发育分期
Developmental stage颜色
Color5 IV, IV—V, V 白色，乳白色 IV, V 草绿色 6 V, V—VI 白色，乳白色 IV, V 草绿色，黄绿色 7 V, V—VI, VI 白色，淡红色 IV, V, VI 草绿色，黄绿色 8 V—VI, VI 白色，淡红色 V—VI, VI 草绿色，黄绿色 9 VI, VI—II 淡红色 VI, VI—II 浅绿色，黄绿色 10 II 淡红色，浅肉红色 II 浅绿色，红黄色 11 II, II—III 浅肉红色 II, II—III 红黄色，浅黄色 12 II, II—III 肉红色 II, II—III 黄色，黄绿色 1 III 肉红色 III, III—IV 黄绿色，浅绿色 2 III—IV 深肉红色，白色 III, IV 黄绿色，浅绿色 3 III, IV 白色，乳白色 III, IV 黄绿色，浅绿色 4 III, IV 白色，乳白色 III, IV 黄绿色，浅绿色
下载: 导出CSV
表 4 鲇 1—12 月份的性成熟系数
Table 4 Gonadosomatic indexes of S. asotus from January to December
月份
Month体质量
Body mass/g空壳质量
Eviscerated mass/g性腺质量
Gonad mass/g性成熟系数
Gonadosomatic index肥满度
Fullness雌
Female雄
Male雌
Female雄
Male雌
Female雄
Male雌
Female雄
Male雌
Female雄
Male1 264.97 193.43 238.56 171.50 16.05 1.82 6.73 1.06 0.79 0.79 2 327.10 263.24 234.84 232.75 17.97 2.85 7.65 1.22 0.87 0.79 3 345.19 268.75 306.13 249.80 26.34 4.67 8.60 1.87 0.89 0.68 4 428.63 322.28 379.25 295.21 30.58 7.42 8.06 2.51 0.92 0.75 5 264.45 234.81 215.14 206.07 32.79 8.74 15.24 4.24 1.37 1.19 6 351.66 312.57 297.71 283.97 34.80 11.85 11.69 4.17 1.20 1.03 7 298.60 286.63 257.36 250.69 24.51 4.08 9.52 1.63 0.83 0.72 8 367.80 313.60 312.88 288.56 27.59 4.23 8.82 1.47 0.91 0.79 9 194.47 160.90 171.68 144.52 3.50 1.63 2.04 1.13 0.97 0.84 10 357.60 298.38 319.83 272.31 12.621 3.20 3.95 1.18 0.99 0.86 11 421.08 321.36 367.71 315.17 13.66 2.34 3.71 0.74 0.90 0.77 12 411.69 359.78 365.84 310.28 15.06 3.12 4.11 1.01 0.65 0.56
下载: 导出CSV
表 5 鲇繁殖力和卵径
Table 5 Fecundity and egg diameter of S. asotus
月份
Month绝对繁殖力
Absolute fecundity/粒相对繁殖力
Relative fecundity/粒卵径范围
Egg diameter range/μm平均卵径
Average egg diameter/μm5 57 223 266 998.5~1 908.4 1 626.5±101.2a 6 67 744 228 974.2~1 861.6 1 588.6±66.2b 7 42 166 164 856.3~1 601.6 1 259.3±133.6c 8 31 131 99 794.4~1 303.7 1 303.6±88.8c
下载: 导出CSV
-
[1] 龙昱, 刘少军. 本地鲇鱼生物学特性及性腺显微结构初步研究[J]. 生命科学研究, 2006, 10(3): 125-129. [2] 温海深, 毛玉泽. 野生鲇资源开发与养殖[J]. 现代渔业信息, 2003, 18(8): 6-9. [3] 吴士平, 李海洋. 鲇鱼的人工繁殖技术[J]. 水产养殖, 2001, 12(4): 37. [4] 张耀光, 谢小军. 南方鲇的繁殖生物学研究: 性腺发育及周年变化[J]. 水生生物学报, 1996, 20(1): 9-14. [5] 李刚, 胡国宏, 顾权, 等. 怀头鲇人工繁殖及苗种培育技术[J]. 内陆水产, 2006(3): 34-35. [6] 兰国柱, 李海滨, 武金星. 黄河陕西段兰州鲇的人工繁殖技术[J]. 当代畜禽养殖业, 2020(4): 63-64. [7] 韩庆, 孙东方, 彭友林, 等. 洞庭湖鲶鱼的人工繁殖与苗种培育试验[J]. 湖北农业科学, 2009, 48(11): 2801-2803. [8] 刘飞, 韦慧, 顾党恩, 等. 流溪河入侵鱼类豹纹脂身鲶繁殖生物学研究[J]. 淡水渔业, 2017, 47(2): 42-48. [9] 薛晨江, 雷春云, 符世伟. 中华鲱鲇繁殖生物学研究[J]. 水生态学杂志, 2017, 38(6): 102-106. [10] 陈得仿, 王腾, 刘永, 等. 大亚湾黑鲷繁殖特征的生物学研究[J]. 南方水产科学, 2019, 15(5): 41-47. [11] 陈晓娟, 杨寅群, 朱迪, 等. 沅江流域生态环境保护体系研究[J]. 人民长江, 2022, 53(2): 11-14. [12] 殷名称. 鱼类生态学[M]. 北京: 中国农业出版社, 1996: 1-4. [13] 茹辉军, 刘学勤, 黄向荣, 等. 大型通江湖泊洞庭的鱼类物种多样性及其时空变化[J]. 湖泊科学, 2008, 20(1): 93-99. [14] TURGEON K, TURPIN C, GREAVESI E I. Dams have varying impacts on fish communities across latitudes: a quantitative synthesis[J]. Ecol Lett, 2019, 22: 1501-1516. doi: 10.1111/ele.13283
[15] ARTHINGTON A H, DULVY N K, GLADSTONE W, et al. Fish conservation in freshwater and marine realms: status, threats and management[J]. Aquat Conserv: Mar Freshw Ecosys, 2016, 26(5): 838-857. doi: 10.1002/aqc.2712
[16] 朱书礼, 李新辉, 李跃飞, 等. 珠江中下游广东鲂种群结构及生长的空间差异[J]. 华中农业大学学报, 2023, 42(1): 75-81. [17] 郜星晨, 姜伟, 张琪, 等. 长江宜昌段鳙的繁殖生物学初步研究[J]. 水生态学杂志, 2019, 40(6): 92-98. [18] 刘艳超, 陈锋, 翁仕洋, 等. 西藏曲松下洛国家湿地公园短尾高原鳅的体长体质量关系和性腺相关特征[J]. 高原科学研究, 2021(4): 44-49. [19] 黄小林, 杨育凯, 李涛, 等. 池塘养殖黄斑篮子鱼初次性成熟性腺发育研究[J]. 南方水产科学, 2020, 16(5): 99-107. [20] 肖智. 鲇繁殖习性的研究[J]. 中山大学学报论丛, 2000, 20(5): 41-44. [21] 曹希全, 李艳慧, 魏杰, 等. 宽口裂腹鱼体长-体质量关系和肥满度[J]. 西北农业学报, 2019, 28(8): 1380-1386. [22] 田娜娜, 杨瑞斌, 谭博真, 等. 西藏哲古措异尾高原鳅年龄、生长及繁殖特性[J]. 中国水产科学, 2022, 29(7): 1013-1021. [23] 朱邦科, 谢从新, 王明学, 等. 保安湖沙塘鳢的食性、繁殖、年龄及生长的研究[J]. 水生生物学报, 1999, 23(4): 316-323. [24] 何福林, 向建国. 沟鲇的生物性状及生化特性研究[J]. 水利渔业, 2006, 26(6): 54-55. [25] 沈勇, 邱其浚, 孙龙生, 等. 饲料精氨酸与赖氨酸配比对全雄黄颡鱼生长性能、体组成、血清生化指标及氨基酸沉积率的影响[J]. 动物营养学报, 2017, 29(7): 2575-2586. [26] 高小平, 吴金明, 孔赤平, 等. 鄱阳湖刀鲚繁殖群体生物学特征[J]. 水产学杂志, 2022, 35(2): 42-46. [27] 王静安, 陈开健, 刘小燕, 等. 大通湖水域河蟹体长、体质量及头胸甲关系研究[J]. 湖南文理学院学报(自然科学版), 2017, 39(1): 35-39. [28] 杨丽亚, 吕红健, 付梅, 等. 三峡库区短颌鲚年龄和生长特性的研究[J]. 水生生物学报, 2022, 46(1): 17-28. [29] 陆敬波, 姚俊杰, 毛国栋, 等. 网箱养殖福瑞鲤鱼种体长与体质量关系分析[J]. 河北渔业, 2017(7): 28-29, 37. [30] LI L, MA B, WANG J L, et al. Length weight relationships of four fish species from the MoTuo reach of the Yarlung Zangbo River, southwest China[J]. J Appl Ichthyol, 2016, 32(6): 1313-1314. doi: 10.1111/jai.13120
[31] XIE J Y, YANG Y H, LIN S Q. Length-weight and length-length relationships of four fish species from rivers in the Tibetan Plateau, China[J]. J Appl Ichthyol, 2018, 34(6): 1345-1347. doi: 10.1111/jai.13787
[32] 栾静, 徐宾铎, 薛莹, 等. 海州湾方氏云鳚体长与体质量分布特征及其关系[J]. 中国水产科学, 2017, 24(6): 1323-1331. [33] 潘文杰, 高雷, 杨浩, 等. 长江中游宜昌—荆州江段鲢种群年龄结构和生长特征[J]. 中国水产科学, 2019, 26(2): 362-370. [34] 朱立新, 侯刚, 梁振林. 基于贝叶斯方法的黄海北部鳀鱼体长与体质量关系式参数估计[J]. 中国水产科学, 2015, 22(4): 757-769. [35] 张国政, 李显森, 朱建成, 等. 黄海中南部小黄鱼当年幼鱼生长特征的研究[J]. 渔业科学进展, 2010, 31(3): 15-22. [36] FROESE R. Cube law, condition factor and weight-length relationships: history, meta-analysis and recommendations[J]. J Appl Ichthyol, 2006, 22(4): 241-253. doi: 10.1111/j.1439-0426.2006.00805.x
[37] 朱玲, 温海深, 毛玉泽. 野生鲇性腺发育及其与肝重指数关系的研究[J]. 水利渔业, 2002, 22(5): 26-27. doi: 10.3969/j.issn.1003-1278.2002.05.012 [38] 蔡子德, 林岗, 倪家延, 等. 光倒刺鲃的繁殖生物学研究[J]. 广西农业科学, 2007, 38(2): 201-203. [39] LAMBERT Y. Why should we closely monitor fecundity in marine fish populations?[J]. J Northwest Atlantic Fish Sci, 2008, 41: 93-106. doi: 10.2960/J.v41.m628
[40] CAO L, SONG B Y, ZHA J M, et al. Age composition, growth, and reproductive biology of yellow catfish (Peltobagrus fulvidraco, Bagridae) in Ce Lake of Hubei Province, Central China[J]. Environ Biol Fish, 2008, 86(1): 75-88.
[41] 李忠利, 冉辉, 杨马, 等. 锦江翘嘴鲌的繁殖生物学特征[J]. 动物学杂志, 2017, 52(2): 263-270. -
期刊类型引用(1)
1. 赵玉薇，黄成，邢佳敏，陈兴汉，张虹. 疣吻沙蚕(Tylorrhynchus heterochaetus)成熟精子及卵子的超微结构观察. 福建农林大学学报(自然科学版). 2024(06): 792-796 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 314
- HTML全文浏览量: 131
- PDF下载量: 83
- 被引次数: 2
粤公网安备 44010502001741号
