Physiological response of Sargassum hemiphyllum to cadmium stress
-
摘要:
实验选取广东省沿海常见物种半叶马尾藻 (Sargassum hemiphyllum)为研究对象，采用室内半静态的实验方法，研究了半叶马尾藻对不同浓度重金属镉离子 (Cd2+) 胁迫的生理响应和耐受性，为广东省沿海海区重金属污染修复提供依据。结果显示，当ρ(Cd2+)为0.1 mg·L−1时，半叶马尾藻藻体中叶绿素a (Chl a)、类胡萝卜素 (Car)、可溶性蛋白 (SP)、可溶性糖 (SS)、抗氧化酶[超氧化物歧化酶 (SOD) 和过氧化氢酶 (CAT)]活性与对照组之间均无显著差异 (P>0.05)；当ρ(Cd2+)为0.5 ~12.5 mg·L−1时，藻体中Chl a、Car、SP、SS、SOD和CAT活性均显著增加 (P<0.05)。当ρ(Cd2+)为0.1~0.5 mg·L−1时，藻体的丙二醛 (MDA) 摩尔质量浓度无显著差异 (P>0.05)；当ρ(Cd2+)为2.5 ~12.5 mg·L−1时，胁迫组藻体的MDA摩尔质量浓度显著升高 (P<0.05)。因此，半叶马尾藻在ρ(Cd2+)≤0.5 mg·L−1下具有较强的抗逆性，而过高浓度Cd2+胁迫 (>0.5 mg·L−1)会对藻体的光合作用和抗氧化功能产生负面影响。
Abstract:Sargassum hemiphyllum, a common species in coastal waters of Guangdong Province, was selected to study its physiological response to different concentrations of heavy metal (Cd2+) stress and its tolerance to Cd2+ stress in laboratory semi-static experiments, so as to provide a basis for remediation of heavy metal pollution in that area. The results show that chlorophyll a (Chl a), carotenoid (Car), soluble protein (SP), soluble sugar (SS) and antioxidant enzyme activity (SOD and CAT) of S. hemiphyllum at Cd2+ concentration of 0.1 mg∙L−1
had no significant difference with the control (P>0.05). The contents of Chl a, Car, SP, SS, SOD and CAT at ρ(Cd2+) of 0.5−12.5 mg∙L−1, were significantly higher than those of the control (P<0.05); there was no significant change in MDA contents under Cd2+ concentration of 0.1−0.5 mg∙L−1 (P>0.05), and significant increase at Cd2+ concentration of 2.5−12.5 mg∙L−1 (P<0.05). Therefore, S. hemiphyllum had stronger stress resistance at Cd2+ concentration of ≤0.5 mg∙L−1, while excessive Cd2+ stress of >0.5 mg∙L−1 will be harmful for its photosynthesis and antioxidant function.
-
疣吻沙蚕 (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. 结论
疣吻沙蚕基因组属于高杂合、高重复的复杂基因组，其微卫星位点的类型丰富且具备较好的多态性潜能，可作为有效资源用于微卫星标记的大规模开发，对种质资源评价与保护利用、种群遗传学以及分子育种研究具有实际价值。
-
![]()
图 1 重金属镉离子对半叶马尾藻叶绿素a和类胡萝卜素的影响
不同字母表示同一胁迫时间内不同处理组间有显著性差异 (P<0.05)，后图同此
Figure 1. Effects of Cd2+ on Chl a and Car of S. hemiphyllum
Different letters represent significant difference among different treatments at the same time (P<0.05). The same case in the following figures.
![]()
图 2 重金属镉离子对半叶马尾藻可溶性糖和可溶性蛋白的影响
Figure 2. Effects of Cd2+ on soluble sugar and soluble protein of S. hemiphyllum
-
[1] PAN K, WANG W X. Trace metal contamination in estuarine and coastal environments in China[J]. Sci Total Environ, 2012: 3-16.
[2] KAPPEL C V, HALPERN B S, SELKOE K A, et al. A global map of human impact on marine ecosystems[J]. Science, 2008, 319(5865): 948-952. doi: 10.1126/science.1149345
[3] GU Y G, HUANG H H, LIU Y, et al. Non-metric multidimensional scaling and human risks of heavy metal concentrations in wild marine organisms from the Maowei Sea, the Beibu Gulf, South China Sea[J]. Environ Toxicol Pharmacol, 2018, 59: 119-124. doi: 10.1016/j.etap.2018.03.002
[4] WANG S L, XU X R, SUN Y X, et al. Heavy metal pollution in coastal areas of South China: a review[J]. Mar Pollut Bull, 2013, 76(1): 7-15.
[5] GU Y G, LI Q S, FANG J H, et al. Identification of heavy metal sources in the reclaimed farmland soils of the pearl river estuary in China using a multivariate geostatistical approach[J]. Ecotox Environ Safe, 2014, 105: 7-12. doi: 10.1016/j.ecoenv.2014.04.003
[6] 李慧卿, 孙泽臣, 常丹, 等. 毛建对重金属Ni2+和Cu2+吸附的研究[J]. 食品研究与开发, 2016, 37(10): 42-45, 46. doi: 10.3969/j.issn.1005-6521.2016.10.011 [7] 王荣霞, 黄敏, 陈傅晓, 等. 棒叶蕨藻变种对重金属Cu2+、Pb2+和Cd2+胁迫的生理响应[J]. 渔业现代化, 2017, 44(2): 25-29, 37. doi: 10.3969/j.issn.1007-9580.2017.02.005 [8] 黄鹤忠. 江蓠对N、P和重金属Cd2+、Cr6+、Cu2+、Ni2+污染的去除效应及其生理生化响应[D]. 扬州: 扬州大学, 2013: 92-104. [9] 罗正明, 贾雷坡, 刘秀丽, 等. 水环境镉对水生动物毒性的研究进展[J]. 食品工业科技, 2015, 36(15): 376-381. [10] LU G Y, KE C H, ZHU A, et al. Oyster-based national mapping of trace metals pollution in the Chinese coastal waters[J]. Environ Pollut, 2017, 224: 658-669. doi: 10.1016/j.envpol.2017.02.049
[11] 巩慧敏, 刘永, 肖雅元, 等. 近海海水和表层沉积物重金属污染与生态风险评价——以海南新村港为例[J]. 农业现代化研究, 2018, 39(4): 700-708. [12] 雷富, 陈宪云, 许铭本, 等. 广西茅尾海海水和表层沉积物中重金属污染的调查及评价[J]. 广西科学院学报, 2013, 29(3): 176-180, 185. doi: 10.3969/j.issn.1002-7378.2013.03.007 [13] 刘慧杰, 刘文君, 刘继平, 等. 南中国海表层海水重金属含量及其潜在生态风险分析[J]. 中国环境科学, 2017, 37(10): 3891-3898. doi: 10.3969/j.issn.1000-6923.2017.10.034 [14] 吕芳, 詹冬梅, 郭文, 等. 重金属镉对铜藻生长及其生物活性成分积累的影响[J]. 大连海洋大学学报, 2017(6): 664-670. [15] 尹文珂, 程金凤, 肖婉露, 等. 四尾栅藻对重金属镉胁迫的响应[J]. 农业环境科学学报, 2015, 34(4): 633-638. [16] 韩婷婷, 付贵权, 齐占会, 等. 充气培养对半叶马尾藻生长、营养盐吸收和生化组成的影响[J]. 中国水产科学, 2015, 22(2): 311-318. [17] 韩婷婷, 齐占会, 黄洪辉, 等. 3种氮源加富对半叶马尾藻幼苗生长和生化组成的影响[J]. 渔业科学进展, 2018, 39(4): 159-166. [18] 付贵权, 韩婷婷, 巩秀玉, 等. 重金属Cu2+胁迫对半叶马尾藻生长及生理生化特性的影响[J]. 南方水产科学, 2015, 11(4): 34-39. doi: 10.3969/j.issn.2095-0780.2015.04.005 [19] 韩婷婷, 方建光, 张继红, 等. 静止和充气培养条件下光强对龙须菜光合生理效应的影响[J]. 水产学报, 2013, 37(12): 1801-1806. [20] WELLBUM A R. The spectral determination of chlorophylls a and b, as well as total carotenoids, using various solvents with spectrophotometers of different resolution[J]. J Plant Physiol, 1994, 144(3): 307-313. doi: 10.1016/S0176-1617(11)81192-2
[21] KOEHERT G. Carbohydrate determination by phenol-sulfuric acid method [M]. London: Cambridge University Press, 1978: 95-97.
[22] 李恒, 李美真, 曹婧, 等. 温度对几种大型海藻硝氮吸收及其生长的影响[J]. 渔业科学进展, 2013, 34(1): 159-165. doi: 10.3969/j.issn.1000-7075.2013.01.024 [23] 黄中坚, 钟志海, 宋志民, 等. 不同营养盐水平对芋根江蓠的生长及生化组分的影响[J]. 南方水产科学, 2014, 10(5): 30-38. doi: 10.3969/j.issn.2095-0780.2014.05.005 [24] 王增焕, 林钦, 李刘冬, 等. 大型海藻对重金属镉、铜的富集动力学研究[J]. 中国环境科学, 2013, 33(1): 154-160. doi: 10.3969/j.issn.1000-6923.2013.01.022 [25] 吕利云. 几种大型海藻对海水中重金属污染的生物修复能力研究[D]. 青岛: 中国海洋大学, 2013: 24-51. [26] 屈肖肖. 重金属 Pb、Cd 胁迫下龙须菜的富集效应及差异 [D]. 汕头: 汕头大学, 2012: 38-33. [27] 余江, 杨宇峰, 聂湘平. 大型海藻龙须菜对重金属镉胁迫的响应[J]. 四川大学学报(工程科学版), 2007, 39(3): 83-90. doi: 10.3969/j.issn.1009-3087.2007.03.016 [28] 朱喜锋. 重金属汞、铜和镉对三种大型经济海藻毒性效应的研究[D]. 汕头: 汕头大学, 2010: 4-6. [29] 周文彬, 邱保胜. 藻类对重金属的耐性与解毒机理[J]. 湖泊科学, 2004, 16(3): 265-272. doi: 10.3321/j.issn:1003-5427.2004.03.012 [30] 戴灵豪, 史冬玲, 王林燕, 等. 镉胁迫下植物膜蛋白的研究进展[J]. 农业工程, 2017(6): 153-156. doi: 10.3969/j.issn.2095-1795.2017.06.045 [31] TUDOREANU L, PHILLIPS C J C. Modeling cadmium uptake and accumulation in plants[J]. Adv Agron, 2004, 84(4): 121-157.
[32] COLLEN J, PINTO E, PEDERSEN M. Induction of oxidative stress in the red macroalga gracilaria tenuistipitata, by pollutant metals[J]. Arch Environ Contam Toxicol, 2003, 45(3): 337-342.
-
期刊类型引用(1)
1. 赵玉薇，黄成，邢佳敏，陈兴汉，张虹. 疣吻沙蚕(Tylorrhynchus heterochaetus)成熟精子及卵子的超微结构观察. 福建农林大学学报(自然科学版). 2024(06): 792-796 . 
百度学术
其他类型引用(1)
下载:
计量
- 文章访问数: 4604
- HTML全文浏览量: 1874
- PDF下载量: 35
- 被引次数: 2
粤公网安备 44010502001741号
