Characterization of microsatellites and polymorphic marker development in ragworm (Tylorrhynchus heterochaetus) based on genome survey data
-
摘要:
为了解疣吻沙蚕 (Tylorrhynchus heterochaetus) 基因组信息并高效地开发微卫星标记，指导其种质资源保护与新品种的遗传改良研究，采用低深度高通量测序开展全基因组survey，k-mer分析估计疣吻沙蚕基因组大小为759.53 Mb，杂合率1.41%，重复序列比例45.92%；初步组装获得2 181 621条scaffold，全长为840 375 821 bp。在基因组序列中检测到130 216个微卫星位点，丰度为154.9 个·Mb−1。微卫星重复次数集中在4~18拷贝；单碱基重复比例最高 (35.00%)，二碱基 (32.48%)、三碱基 (14.42%) 次之；二碱基、三碱基优势基序分别是AT/AT、AAT/ATT，表现出A/T碱基优势。从随机挑选的50对引物中筛选到15对多态标记，在30尾样本中共检测到87个等位基因，等位基因数 (Na) 为2.000~12.000 (平均5.800)，有效等位基因数 (Ne) 为1.164~6.713 (平均3.328)，期望杂合度 (He) 为0.141~0.789 (平均0.561)，多态信息含量 (PIC) 为0.136~0.776 (平均0.511)；其中13个为高度或中度多态性位点，在遗传分析中有较高的实用价值。结果表明，疣吻沙蚕基因组为复杂基因组，其微卫星位点类型丰富且具备良好的多态性潜能，可为种质资源评价、群体遗传学及分子育种研究提供有效的标记资源。
-
关键词:
- 疣吻沙蚕 /
- 全基因组survey /
- 简单重复序列 /
- 多态性位点
Abstract:In order to understand the genomic information of Tylorrhynchus heterochaetus and efficiently develop microsatellite markers, so as to guide the conservation of its germplasm resources and genetic improvement of new varieties, we conducted a whole-genome survey by using low depth high-throughput sequencing. A total of 57.48 Gb of clean data were generated after the quality control of raw data. K-mer analysis estimates that the genome size of T. heterochaetus was 759.53 Mb; the heterozygosity rate was 1.41%; the proportion of repetitive sequences was 45.92%. Preliminary assembly obtained 2 181 621 scaffolds with a total length of 840 375 821 bp. A total of 130 216 microsatellite loci were detected with a density of 154.9 loci per Mb. The repeated number of microsatellite units largely ranged from 4 to 18. The ratio of mononucleotide loci was the highest (35.00%), followed by those of dinucleotide (32.48%) and trinucleotide (14.42%) loci. AT/AT and AAT/ATT motifs were dominant in dinucleotide and trinucleotide loci, respectively, indicating an A/T dominance. Fifteen polymorphic loci were identified from 50 randomly selected primers, and 87 alleles were amplified in a T. heterochaetus population containing 30 individuals. The number of alleles per locus ranged from 2.000 to 12.000, with an mean of 5.800. The effective allele number (Ne) and expected heterozygosity (He) ranged from 1.164 to 6.713 and from 0.141 to 0.789, with means of 3.328 and 0.561, respectively. The polymorphic information content (PIC) ranged from 0.136 to 0.776, with a mean of 0.511. Thirteen loci were found to be highly or moderately polymorphic, having high practical value in genetic analysis. In conclusion, T. heterochaetus genome is a complex genome, and its microsatellites have a rich variety and high polymorphic potential. The results can provide effective marker resources for germplasm resource evaluation, population genetics and molecular breeding research.
-
疣吻沙蚕 (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 疣吻沙蚕基因组 k-mer 种类频率分布
Figure 1. Frequency distribution of k-mer species in genome of T. heterochaetus
![]()
图 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
![]()
图 4 疣吻沙蚕基因组微卫星重复基序类型分布特征
Figure 4. Distribution pattern of microsatellite motif types in genome of T. heterochaetus
![]()
图 5 疣吻沙蚕基因组微卫星长度分布特征
注：a. 不同长度区间微卫星数量及比例；b. 不同类型微卫星长度分布特征。
Figure 5. Distribution pattern of length of microsatellite loci genome of T. heterochaetus
Note: a. Number and percentage of microsatellite loci at different length intervals; b. Length distribution of the six motif types of microsatellite loci.
![]()
图 6 部分多态微卫星标记的毛细管电泳分型结果
Figure 6. A set of polymorphic microsatellite loci visualized by high-resolution capillary electrophoresis
表 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 
下载: 导出CSV
表 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.
下载: 导出CSV
表 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.
下载: 导出CSV
-
[1] YANG Z Q, SUNIL C, JAYACHANDRAN M, et al. Anti-fatigue effect of aqueous extract of Hechong (Tylorrhynchus heterochaetus) via AMPK linked pathway[J]. Food Chem Toxicol, 2020, 135: 111043. doi: 10.1016/j.fct.2019.111043
[2] 苏跃朋, 黄啟, 崔阔鹏. 珠江河口区禾虫产业技术现状及增养殖效益分析[J]. 海洋与渔业, 2016(10): 64-67. [3] ZHANG W X, WANG Z X, GANESAN K, et al. Antioxidant activities of aqueous extracts and protein hydrolysates from marine worm Hechong (Tylorrhynchus heterochaeta)[J]. Foods, 2022, 11(13): 1837. doi: 10.3390/foods11131837
[4] 杨尉, 陈兴汉. 疣吻沙蚕-水稻生态复合种养技术要点及效益分析[J]. 南方农业, 2022, 16(20): 17-20, 24. [5] CHEN X H, YANG S, YANG W, et al. First genetic assessment of brackish water polychaete Tylorrhynchus heterochaetus: mitochondrial COI sequences reveal strong genetic differentiation and population expansion in samples collected from southeast China and north Vietnam[J]. Zool Res, 2020, 41(1): 61-69. doi: 10.24272/j.issn.2095-8137.2020.006
[6] CHEN X H, LI M M, LIU H P, et al. Mitochondrial genome of the polychaete Tylorrhynchus heterochaetus (Phyllodocida, Nereididae)[J]. Mitochondrial DNA A, 2016, 27(5): 3372-3373. doi: 10.3109/19401736.2015.1018226
[7] CHEN H, LI X, WANG Y, et al. De novo transcriptomic characterization enables novel microsatellite identification and marker development in Betta splendens[J]. Life, 2021, 11(8): 803. doi: 10.3390/life11080803
[8] 孙效文, 张晓锋, 赵莹莹, 等. 水产生物微卫星标记技术研究进展及其应用[J]. 中国水产科学, 2008, 15(4): 689-703. [9] 张永德, 文露婷, 罗洪林, 等. 卵形鲳鲹基因组调研及其SSR分子标记的开发应用[J]. 南方农业学报, 2020, 51(5): 983-994. [10] 上官清, 陈昆慈, 刘海洋, 等. 斑鳢基因组中微卫星分布特征及野生种群遗传结构分析[J]. 南方水产科学, 2020, 16(3): 47-60. [11] LIU B H, SHI Y J, YUAN J Y, et al. Estimation of genomic characteristics by analyzing k-mer frequency in de novo genome projects[J]. Quant Biol, 2013, 35(s1-3): 62-67.
[12] LUO R B, LIU B H, XIE Y L, et al. SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler[J]. GigaScience, 2012, 1(1): 18. doi: 10.1186/2047-217X-1-18
[13] LALITHA S. Primer premier 5[J]. Biotech Softw Internet Rep, 2000, 1(6): 270-272. doi: 10.1089/152791600459894
[14] 刘玉萍, 王棋, 黄新芯, 等. 基于高通量测序的带鱼肌肉组织转录组微卫星信息分析[J]. 南方农业学报, 2022, 53(3): 725-734. [15] PEAKALL R, SMOUSE P E. GenAlEx 6.5: genetic analysis in Excel. Population genetic software for teaching and research: an update[J]. Bioinformatics, 2012, 28(19): 2537-2539. doi: 10.1093/bioinformatics/bts460
[16] TEMNYKH S, DECLERCK G, LUKASHOVA A, et al. Computational and experimental analysis of microsatellites in rice (Oryza sativa L.): frequency, length variation, transposon associations, and genetic marker potential[J]. Genome Res, 2001, 11(8): 1441-1452. doi: 10.1101/gr.184001
[17] SIMAKOV O, MARLETAZ F, CHO S J, et al. Insights into bilaterian evolution from three spiralian genomes[J]. Nature, 2013, 493(7433): 526-531. doi: 10.1038/nature11696
[18] TONG L, DAI S X, KONG D J, et al. The genome of medicinal leech (Whitmania pigra) and comparative genomic study for exploration of bioactive ingredients[J]. BMC Genom, 2022, 23(1): 76. doi: 10.1186/s12864-022-08290-5
[19] MARTÍN-ZAMORA F M, LIANG Y, GUYNES K, et al. Annelid functional genomics reveal the origins of bilaterian life cycles[J]. Nature, 2023, 615(7950): 105-110. doi: 10.1038/s41586-022-05636-7
[20] de OLIVEIRA A L, MITCHELL J, GIRGUIS P, et al. Novel insights on obligate symbiont lifestyle and adaptation to chemosynthetic environment as revealed by the giant tubeworm genome[J]. Mol Biol Evol, 2022, 39(1): msab347. doi: 10.1093/molbev/msab347
[21] LI Y, TASSIA M G, WAITS D S, et al. Genomic adaptations to chemosymbiosis in the deep-sea seep-dwelling tubeworm Lamellibrachia luymesi[J]. BMC Biol, 2019, 17(1): 91. doi: 10.1186/s12915-019-0713-x
[22] JIN F, ZHOU Z L, GUO Q, et al. High-quality genome assembly of Metaphire vulgaris[J]. PeerJ, 2020, 8: e10313. doi: 10.7717/peerj.10313
[23] ZAKAS C, HARRY N D, SCHOLL E H, et al. The genome of the poecilogonous Annelid Streblospio benedicti[J]. Genome Biol Evol, 2022, 14(2): evac008. doi: 10.1093/gbe/evac008
[24] SHAO Y, WANG X B, ZHANG J J, et al. Genome and single-cell RNA-sequencing of the earthworm Eisenia andrei identifies cellular mechanisms underlying regeneration[J]. Nat Commun, 2020, 11(1): 2656. doi: 10.1038/s41467-020-16454-8
[25] ZWARYCZ A S, NOSSA C W, PUTNAM N H, et al. Timing and scope of genomic expansion within Annelida: evidence from homeoboxes in the genome of the earthworm Eisenia fetida[J]. Genome Biol Evol, 2016, 8(1): 271-281. doi: 10.1093/gbe/evv243
[26] KENNY N J, NAMIGAI E K O, MARLÉTAZ F, et al. Draft genome assemblies and predicted microRNA complements of the intertidal lophotrochozoans Patella vulgata (Mollusca, Patellogastropoda) and Spirobranchus (Pomatoceros) lamarcki (Annelida, Serpulida)[J]. Mar Genom, 2015, 24(2): 139-146.
[27] SUN Y N, SUN J, YANG Y, et al. Genomic signatures supporting the symbiosis and formation of chitinous tube in the deep-sea tubeworm Paraescarpia echinospica[J]. Mol Biol Evol, 2021, 38(10): 4116-4134. doi: 10.1093/molbev/msab203
[28] 高胜寒, 禹海英, 吴双阳, 等. 复杂基因组测序技术研究进展[J]. 遗传, 2018, 40(11): 944-963. [29] 徐杰杰, 毕宜慧, 程景颢, 等. 中华绒螯蟹 (Eriocheir sinensis) 全基因组微卫星分布特征研究[J]. 基因组学与应用生物学, 2021, 40(Z2): 2422-2429. [30] 梁霞, 王慧琪, 马宇璇, 等. 鲤鱼(Cyprinus carpio)全基因组微卫星分布特征研究[J]. 南京师大学报 (自然科学版), 2021, 44(3): 103-111. [31] ZHANG Q, ZHANG C S, YU Y, et al. Characteristic analysis of simple sequence repeats in the ridgetail white prawn Exopalaemon carinicauda genome and its application in parentage assignment[J]. J World Aquacult Soc, 2020, 51(3): 690-701. doi: 10.1111/jwas.12650
[32] SRIVASTAVA S, KUSHWAHA B, PRAKASH J, et al. Development and characterization of genic SSR markers from low depth genome sequence of Clarias batrachus (Magur)[J]. J Genet, 2016, 95(3): 603-609. doi: 10.1007/s12041-016-0672-8
[33] 彭冶, 李杰, 王涛, 等. 瓦氏黄颡鱼全基因组微卫星的分布特征及其定位的初步研究[J]. 南方水产科学, 2022, 18(1): 90-98. [34] XU S Y, SONG N, XIAO S J, et al. Whole genome survey analysis and microsatellite motif identification of Sebastiscus marmoratus[J]. Biosci Rep, 2020, 40(2): BSR20192252. doi: 10.1042/BSR20192252
[35] 王九龙, 李洪莉, 尹硕, 等. 绿鳍马面鲀全基因组微卫星分布特征[J]. 烟台大学学报 (自然科学与工程版), 2022, 35(3): 285-293. [36] 王佳佳, 王琼, 秦桢, 等. 凡纳滨对虾全基因组SSR标记开发及不同养殖群体的遗传多样性分析[J]. 水产学报, 2023, 47(6): 64-74. [37] SUN J X, PENG G H, XIONG L J, et al. Genome-wide SSR marker development and application in genetic diversity analysis of the red swamp crayfish, Procambarus clarkii (Girard, 1852) in China[J]. Crustaceana, 2021, 94(2): 189-205. doi: 10.1163/15685403-bja10076
[38] 倪守胜, 杨钰, 柳淑芳, 等. 基于高通量测序的虾夷扇贝基因组微卫星特征分析[J]. 渔业科学进展, 2018, 39(1): 107-113. [39] 熊良伟, 王帅兵, 岳丽佳, 等. 宽体金线蛭基因组SSR序列特征分析及其分子标记开发[J]. 南方农业学报, 2018, 49(11): 2298-2303. [40] LIU H Y, ZHANG Y F, WANG G B, et al. Development and characterization of microsatellite markers in the earthworm Drawida gisti Michaelsen, 1931 and cross-amplification in two other congeners[J]. Mol Biol Rep, 2020, 47(10): 8265-8269. doi: 10.1007/s11033-020-05799-4
[41] 王斌, 孙静, 刘凌云, 等. 蛭类转录组中EST-SSR分析及抗凝血相关分子标记的挖掘[J]. 中草药, 2017, 48(1): 172-178. [42] MADUNA S N, VIVIAN-SMITH A, JÓNSDÓTTIR Ó D B, et al. Genome- and transcriptome-derived microsatellite loci in lumpfish Cyclopterus lumpus: molecular tools for aquaculture, conservation and fisheries management[J]. Sci Rep, 2020, 10(1): 559. doi: 10.1038/s41598-019-57071-w
[43] 李强勇, 李旻, 曾地刚, 等. 凡纳滨对虾微卫星分子标记的开发及不同养殖家系遗传多态性分析[J]. 南方农业学报, 2020, 51(2): 429-436. [44] WIERDL M, DOMINSKA M, PETES T D. Microsatellite instability in yeast: dependence on the length of the microsatellite[J]. Genetics, 1997, 146(3): 769-779. doi: 10.1093/genetics/146.3.769
[45] JO E, LEE S J, CHOI E, et al. Whole genome survey and microsatellite motif identification of Artemia franciscana[J]. Biosci Rep, 2021, 41(3): BSR20203868. doi: 10.1042/BSR20203868
[46] SCHLÖTTERER C, TAUTZ D. Slippage synthesis of simple sequence DNA[J]. Nucleic Acids Res, 1992, 20(2): 211-215. doi: 10.1093/nar/20.2.211
[47] 马军, 刘嘉鑫, 江智景, 等. 基于RNA-seq数据的密斑刺鲀SSR分子标记开发及鉴定[J]. 南方水产科学, 2020, 16(1): 127-136. [48] 朱维岳, 周桃英, 钟明, 等. 基于遗传多样性和空间遗传结构的野生大豆居群采样策略[J]. 复旦学报 (自然科学版), 2006, 45(3): 321-327. -
期刊类型引用(1)
1. 赵玉薇，黄成，邢佳敏，陈兴汉，张虹. 疣吻沙蚕(Tylorrhynchus heterochaetus)成熟精子及卵子的超微结构观察. 福建农林大学学报(自然科学版). 2024(06): 792-796 . 
百度学术
其他类型引用(1)
计量
- 文章访问数: 281
- HTML全文浏览量: 115
- PDF下载量: 59
- 被引次数: 2
粤公网安备 44010502001741号
