Development and evaluation of di-/tri-nucleotide-repeated microsatellites by RAD-seq in Decapterus macrosoma
-
摘要:
通过对长体圆鲹(Decapterus macrosoma)基因组进行RAD-Seq高通量测序，共获得58 180条微卫星序列，选取112条二、三核苷酸重复的微卫星序列设计引物，经筛选后，共获得27个具有多态性的微卫星标记。利用一个长体圆鲹群体对通过筛选的微卫星标记的种群遗传学特征进行评价。结果显示，27对引物扩增的序列中18个位点为二核苷酸重复，重复次数为9~14次，9个位点为三核苷酸重复，重复次数为6~10次，等位基因数(Na)为5~17 (平均10.6)，表观杂合度(Ho)为0.342 9~0.857 1 (平均0.631 7)，期望杂合度(He)为0.538 3~0.911 8 (平均0.796 8)，多肽信息含量(PIC)为0.497~0.886 (平均0.780 9)，除1个位点外，其他位点PIC值均大于0.500，表明开发的微卫星位点具有较高的多态性。“哈迪-温伯格”平衡(HWE)检测结果显示，19个标记等位基因频率符合HWE。连锁不平衡检测表明各位点间无连锁不平衡现象。该研究开发的27个微卫星标记可为长体圆鲹种群遗传学研究提供基础。
Abstract:A total of 58 180 microsatellites were developed by RAD-Seq high-throughput sequencing technology in Decapterus macrosoma. One hundred and twelve microsatellite loci were randomly selected from di-/tri-nucleotide repeat microsatellite loci. Finally, twenty-seven highly polymorphic markers were developed. The population genetic analysis finds 18 dinucleotide-repeated microsatellite loci (9−14 repeated) and 9 trinucleotide-repeated microsatellite loci (6−10 repeated). The number of alleles (Na) ranged from 5 to 17 (mean 10.6). The observed and expected heterozygosities ranged from 0.342 9 to 0.857 1 (mean 0.631 7) and 0.538 3 to 0.911 8 (mean 0.796 8), respectively. The polymorphism information content (PIC) ranged from 0.497 to 0.886 (mean 0.780 9), and the PIC at all loci except for one was geater than 0.500, indicating that the microsatellites were highly polymorphic. The Hardy-Weinberg equilibrium (HWE) test shows that 19 microsatellite markers conformed to HWE. There is no linkage imbalance among the 27 loci. The 27 microsatellites developed by this study will be useful for further population genetic analysis.
-
长体圆鲹(Decapterus macrosoma)，又名长身圆鲹，隶属于鲈形目、鲹科、圆鲹属，主要分布于中国南海、印度尼西亚、澳洲和日本南部沿海等地[1]。长体圆鲹在中国南海分布较广，是南海灯光围网渔业的主要捕捞对象，具有较高的经济价值[2-5]。目前，国内外学者关于长体圆鲹的研究主要集中在生长繁殖[6]和资源评估[7-8]方面，与种群遗传和分子标记相关的研究报道较少。微卫星标记仅见翟云等[9]开发蓝圆鲹微卫星标记中获得5个跨物种通用标记可于用长体圆鲹，并无专门针对长体圆鲹开发的微卫星标记。种群遗传信息的匮乏，将大大影响对其资源的评估和长期有效的管理。
微卫星分子标记因是共显性标记，具有多态性高、变异性强、数据易统计等突出优点[10]，广泛应用于海洋生物遗传结构及遗传多样性分析[11-12]。但由于微卫星标记通用性较差，常常具有极强的种属特异性。鱼类微卫星标记开发中多以二核苷酸重复为主[13-15]，普遍认为它们具有较高的遗传变异[16]，但是也有部分学者研究认为三、四核苷酸重复位点较二核苷酸重复具有更高的筛选效率和多态性[17-19]。
本研究通过RAD-Seq高通量测序方法开发长体圆鲹二、三核苷酸微卫星分子标记，并对测试群体进行多样性分析，同时比较二、三核苷酸的筛选效率和多态性差异，旨在为长体圆鲹种群遗传结构及遗传多样性分析提供技术基础，并为该资源的评估和管理提供帮助。
1. 材料与方法
1.1 样品采集与基因组DNA提取
长体圆鲹样品采集于中国南海中沙群岛东部海域，共35尾。剪取部分肌肉样品加入无水乙醇保存。每个样品剪取少量肌肉组织，使用海洋动物组织基因组DNA提取试剂盒(天根，北京) 提取基因组DNA，0.8%的琼脂糖凝胶电泳检测DNA提取质量，−20 ℃保存备用。
1.2 高通量测序与引物合成
使用HiSeq2000高通量测序仪(Illumina，USA) 对长体圆鲹基因组DNA进行RAD-seq (测序服务由广州基迪奥生物科技有限公司提供)，经生物信息学搜索出微卫星位点[20]。使用Premier 5.0软件在重复单元侧翼序列上选择性设计出112条引物，主要参数为：G-C含量为40%~60%，引物长度为18~25 bp，退火温度为45~60 ℃，预期产物长度180~320 bp。送上海英潍捷基贸易有限公司合成引物。
1.3 引物筛选与分型检测
选取3个样本混合成的基因组DNA为模板，优化PCR反应条件，对引物进行首轮筛选，琼脂糖电泳检测是否能扩增出稳定且均一的目的片段。之后选取8尾个体的基因组DNA作为模板，使用三引物法[21]，利用M13荧光接头引物进行PCR扩增，扩增产物送华大基因公司经毛细管电泳进行等位基因分型，检测引物是否具有多态性。PCR反应体系为15 μL，其中包括10×PCR Buffer 1.5 μL，2.5 mmol·L–1 MgCL2 1.2 μL，2 mmol·L–1 dNTPs 2 μL，M13正向引物(10 μmol·L–1) 0.2 μL，M13反向引物(10 μmol·L–1) 0.6 μL，M13通用荧光引物(10 μmol·L–1) 0.5 μL，Taq酶(5 U·μL–1) 0.15 μL，DNA模版1 μL，加双蒸水至15 μL。PCR扩增程序为：94 ℃预变性5 min；94 ℃变性45 s，55~60 ℃退火45 s，72 ℃延伸45 s，35个循环；94 ℃变性45 s，53 ℃退火45 s，72 ℃延伸45 s，8个循环；72 ℃延伸30 min。
1.4 长体圆鲹群体遗传学评价
使用35尾长体圆鲹个体的基因组DNA为模板，对通过筛选的微卫星标记的种群遗传学特征进行评价。PCR反应体系和条件、等位基因分型方法如上。使用软件Genepop 4.0[22]对每个标记的种群遗传学特征值进行计算，包括等位基因数(Na)、表观杂合度(Ho)和期望杂合度(He)，进行“哈迪-温伯格”平衡(HWE)检验和连锁不平衡检测，并对P值进行Bonferroni校正。使用Cervus 3.0.7[23]软件计算多态信息含量(PIC)。
2. 结果
2.1 高通量测序结果与微卫星位点分析
RAD-seq高通量测序共获得长体圆鲹基因组原始数据2.06 G，GC含量为41.43%，Q30达93.05%。说明测序结果质量较好，可用于后续分析。搜索后共获得微卫星序列58 180条，一至六核苷酸重复微卫星位点70 508个，其中二核苷酸重复微卫星位点最多(37 646个)，占总数的53.39% (表1)，说明二核苷酸重复为主要的微卫星类型。二核苷酸重复微卫星位点共有4种重复类型，4种类型重复微卫星数量相差较大，AC/GT类有29 754个，占二核苷酸重复的68.4%；AG/CT类有6 487个，占17.2%；AT/TA有1 340个，占3.6%；GC/CG仅有65个，占0.17%。
表 1 长体圆鲹基因组中不同类型SSR统计Table 1. Different types of SSR statistics in D. macrosoma genome重复单元
repeat unit微卫星数量/个
number of microsatellite占比/%
ratio一核苷酸 mono-nucleotide 8 184 11.61 二核苷酸 di-nucleotide 37 646 53.39 三核苷酸 tri-nucleotide 13 960 19.80 四核苷酸 tetra-nucleotide 7 741 10.98 五核苷酸 penta-nucleotide 2 255 3.20 六核苷酸 hexa-nucleotide 722 1.02 合计 total 70 508 100.00 2.2 PCR引物设计和筛选
选取112条二、三核苷酸重复序列设计引物，其中二核苷酸重复为81对，三核苷酸重复为31对。经过筛选后，共有27对引物通过筛选(表2)，27对引物扩增的序列中18个位点为二核苷酸重复，重复次数为9~14次；8个位点为三核苷酸重复，重复次数为6~10次。二核苷酸重复位点检出效率为22.2%，三核苷酸重复位点检出效率为25.8%。
表 2 27对长体圆鲹微卫星引物信息Table 2. Information of 27 pairs of primers in D. macrosoma位点
locus引物序列 (5'−3')
primer sequence重复单元
repeat motif退火温度/℃
annealing temperature期望长度/bp
allele sizeDma03 F:CCACGCCTATTGAGTTACAGA (CA)9 60 186 R:GAGCCAGTGGATGAACAGAGT Dma07 F:GCCCCTGTGGGTGTGTGA (CA)9 60 225 R:GGGTGGTGGGTTCGGTTT Dma12 F:GAACCAGTGCCTACAATAGA (AC)9 60 243 R:CTGCTCACGGTAAGTCCA Dma15 F:ACAGGAAGGAACAGGACAG (TG)10 55 254 R:TATTGAAGTGAAAAAGCCG Dma22 F:CGCTGTTGAAATGAAGAAGA (GT)10 60 317 R:AGTGATGTCGCCTCATAAAT Dma23 F:AAACTGAGGGCGAGATAGAGG (AC)10 55 190 R:CCACAGGCTGAGTAAACCAAC Dma26 F:ATCCCATTCACCGACATAG (TG)10 58 258 R:CTGTGGTATCGTTCCCTGT Dma28 F:TGATTGGCTTCTACTCTGC (AC)10 55 281 R:AGTGGCTTGTTTGACTCTTAT Dma36 F:GGATGTAGTGAAGAGGGGAG (GT)11 55 239 R:CACAATCAGTGTTATGGCAG Dma38 F:GCCAATAAAGGCAAACAGT (CA)11 60 227 R:ATCCGAGACAAAGACATACAA Dma39 F:AGTGTGCTGACTTTTCTCTG (CA)11 55 241 R:TTATTGTTTGTTGTCTGGGT Dma45 F:CTCCTTTTTCTTCTTCCTCT (CA)11 60 281 R:CTACCTGCTCTTCAACTCAT Dma51 F:TGACAGCCTCCACTACTCC (GA)12 55 225 R:GCTAACCAGACACGCAAA Dma54 F:AAAGCCCATCTGTCTCGT (GT)12 60 202 R:TGTTTCAGTCCGTTCCTG Dma58 F:TCAAGAGGGAGTGGGAGC (AC)12 58 279 R:TCAAATGGGTGTTTAGCG Dma64 F:GCTCAGACTGCGTGGACA (TG)13 55 314 R:GCTGGTGAACAACAGGACA Dma72 F:TTCCGCAGGCATAAAAAC (CT)13 58 301 R:CCAAGGTCCGCTACACTA Dma76 F:TTCTCGCTGACCTGCTTG (TG)14 55 253 R:GCGTCCTCGTCGTCTTTC Dma81 F:GAGACACGGTCAGAAAACA (TGC)6 60 216 R:GGAAGTAGGACTCTAGGGG Dma82 F:CTGTCACTCCATTCCTATTCC (GTT)6 58 244 R:CCTACATTTGTGCTTTTGTTC Dma83 F:CTCTAAAGCCGACCTAACC (CTT)6 58 239 R:TGTCTCAACACAGCGAAAC Dma84 F:AAACTAACTCATCACCAG (TGT)6 55 283 R:AAACGACAGGAACTCAAT Dma85 F:CTCACTTTGACCCAACCAG (AGG)6 55 256 R:CCTTTCACCGAGACACCAG Dma131 F:TGCGGATGGGTGGTAGTGT (GGT)8 55 208 R:ATTGCTGGTAGTCGGTGGC Dma132 F:CCCAGTGAGACCAGAACCA (GCT)8 55 268 R:GACCCGTAGACAGGAGAGT Dma135 F:GTTGTTGTTTTTTTCCTT (GCA)9 55 301 R:CATCAGTCTGGCTTTATA Dma145 F:ACGATACAGCAGCCGAAG (TCA)10 60 197 R:AGTGATGTCGCCTCATAAAT 2.3 微卫星标记的种群遗传学评价
使用1个采集自南海东南部海域的长体圆鲹群体对筛选合格的微卫星标记进行种群遗传学评价。所有27个标记在测试群体中共检测到285个等位基因，等位基因数为5~17，Ho为0.342 9~0.857 1，平均为0.631 7；He为0.538 3~0.911 8，平均为0.7968。PIC为0.497~0.886，平均为0.780 9 (表3)，表明开发的微卫星位点具有较高的多态性。共有19个标记等位基因频率符合“哈迪-温伯格”平衡。连锁不平衡检测表明各位点间无连锁不平衡现象。
表 3 长体圆鲹微卫星标记的种群遗传学特征Table 3. Characteristics of microsatellite loci in D. macrosoma位点
locusN Na Ho He PHWE PIC Dma03 35 8 0.857 1 0.790 9 0.042 0 0.746 0 Dma07 35 9 0.771 4 0.855 9 0.553 4 0.825 0 Dma12 34 11 0.685 7 0.816 6 0.106 4 0.809 0 Dma15* 35 14 0.542 9 0.911 8 0.000 0 0.890 0 Dma22* 32 14 0.485 7 0.813 2 0.000 0 0.864 0 Dma23 34 17 0.714 3 0.864 8 0.003 3 0.868 0 Dma26* 34 14 0.428 6 0.869 1 0.000 0 0.870 0 Dma28 35 11 0.742 9 0.837 7 0.035 3 0.803 0 Dma36 35 9 0.685 7 0.786 3 0.054 5 0.743 0 Dma38* 34 12 0.628 6 0.851 2 0.000 0 0.849 0 Dma39* 33 12 0.514 3 0.782 0 0.000 0 0.793 0 Dma45 34 12 0.771 4 0.847 8 0.139 2 0.844 0 Dma51* 31 12 0.485 7 0.784 5 0.000 0 0.859 0 Dma54 35 15 0.771 4 0.900 6 0.139 0 0.878 0 Dma58 34 12 0.771 4 0.845 6 0.449 7 0.843 0 Dma64* 31 13 0.428 6 0.788 3 0.000 0 0.865 0 Dma72 35 8 0.628 6 0.713 9 0.155 7 0.659 0 Dma76 35 12 0.685 7 0.908 1 0.003 6 0.886 0 Dma81 35 8 0.628 6 0.795 0 0.045 7 0.752 0 Dma82 35 5 0.485 7 0.538 3 0.023 6 0.497 0 Dma83* 34 8 0.342 9 0.705 8 0.000 0 0.674 0 Dma84 35 7 0.542 9 0.704 8 0.005 9 0.649 0 Dma85 35 8 0.828 6 0.747 8 0.876 9 0.704 0 Dma131 35 9 0.628 6 0.717 6 0.646 6 0.661 0 Dma132 35 8 0.657 1 0.746 6 0.365 5 0.692 0 Dma135 34 10 0.742 9 0.840 5 0.241 8 0.835 0 Dma145 34 7 0.600 0 0.749 3 0.077 4 0.725 0 注：N. 有效样品数；Na. 等位基因数；Ho. 表观杂合度；He. 期望杂合度；PHWE. “哈迪-温伯格”平衡显著性检验P值；PIC. 多态信息含量；*. 经Bonferroni校正后显著背离“哈迪-温伯格”平衡 (校正P<0.001 85) Note: N. effective number of samples; Na. number of alleles; Ho. observed heterozygosity; He. expected heterozygosity; PHWE. Hardy–Weinberg probability test; PIC. polymorphism information content; *. significant deviation from HWE after Bonferroni's correction (adjusted P-value<0.001 85) 3. 讨论
3.1 高通量测序发掘微卫星序列的技术优势
传统微卫星标记开发方法耗时长、花费高、技术难度大。以磁珠富集法为例，标记开发过程中基因组DNA浓度、接头连接效率、富集过程中的杂交温度以及洗脱条件的控制等因素都会影响微卫星筛选的效率[24-25]，且最终获得的有效微卫星序列仅几百条[26-27]。相比较而言，高通量测序技术开发微卫星标记，省略了建库、克隆、筛选等，只需提取基因组DNA测序，利用生物信息学手段可直接获取微卫星序列，通常是传统方法获得微卫星序列数目的几百倍[15,28-29]，具有高效、便捷、准确的特点，能够满足短时间内大批量微卫星位点的开发需求，比如连锁图谱构建、QTL定位等[30-31]。
本次RAD-seq高通量测序共获得长体圆鲹基因组原始数据2.06 G，GC含量为41.43%，测序质量Q30达93.05%；共获得微卫星序列58 180条，一至六核苷酸重复微卫星位点70 508个。说明测序质量稳定高效，并获得了数量庞大、类型丰富的长体圆鲹微卫星序列，可用于后续长体圆鲹微卫星标记的大规模开发和相关遗传学研究。
3.2 不同核苷酸重复微卫星位点比较
本次高通量测序结果表明在长体圆鲹微卫星位点中二核苷酸重复为主要重复类型，AC/GT类重复数量最为丰富，GC/CG重复较为少见。此结果与大量水产动物微卫星位点研究结果相一致[32-34]，差异仅在于比例多少，以及除二核苷酸重复占主要地位外其他核苷酸重复的含量差异。熊良伟等[33]对中华鳑鲏(Rhodeus sinensis)微卫星的分析中，二核苷酸占总微卫星位点的53.59%，其次为单核苷酸重复，二核苷酸重复中AC/GT类占60.63%，GC/CG仅占0.32%。在裸体异鳔鳅鮀 (Xenophysogobio nudicorpa)中[32]，二核苷酸重复占总微卫星位点比例高达83.15%，AC/GT类重复占49.36%，GC/CG重复仅有4个。
多数鱼类开发的微卫星标记以二核苷酸重复为主，但研究表明，在人类基因组中三核苷酸重复序列与遗传疾病的发生有关，并且具有较高的多态性和遗传稳定性[35]。部分学者对三、四核苷酸重复微卫星标记的研究结果存在差异。房祖业等[28]对大刺鳅 (Mastacembelus armatus) 二、三、四核苷酸重复微卫星标记的筛选发现二核苷酸重复较三、四核苷酸重复具有更高的筛选效率和多态性；鲁翠云等[17]、谭照君等[18]、李文升等[19]的研究认为三、四核苷酸具有更高的多态性和分型效果。长体圆鲹二、三核苷酸的筛选效率分别为22.2%和29.0%，PIC分别为0.827 4和0.687 7 (表4)。就筛选效率而言，三核苷酸重复略高于二核苷酸重复，但二者相差不大。PIC为衡量种群遗传变异程度的重要指标[36]，二核苷酸重复多态性明显高于三核苷酸重复。本文中长体圆鲹二核苷酸重复筛选效率低于三核苷酸重复，但多态性二核苷酸重复明显高于三核苷酸重复。因此，筛选效率和多态性的差异可能由种属差异或其他多种因素导致。
表 4 长体圆鲹二、三核苷酸重复微卫星标记的比较Table 4. Comparison on di- and trinucleotide-repeated microsatellite loci in D. macrosoma序列
sequence引物数
primer number重复次数
repeat number筛选效率
efficiencyPIC 二核苷酸重复
di-nucleotide-repeated18 9~14 22.2% 0.827 4 三核苷酸重复
tri-nucleotide-repeated9 6~10 29.0% 0.687 7 通过筛选的27对引物中18个位点为二核苷酸重复，重复次数为9~14次不等；9个位点为三核苷酸重复，重复次数为6~10次不等，符合Ellegren[37]提出的真核生物微卫星位点重复大部分在30次重复以下。但Ellegren[37]认为二核苷酸重复以15~19次为主，本文中高通量测序获得的二核苷酸重复主要在6~15次。基于Weber[38]的研究结果，重复次数高的微卫星在种群中表现出的多态性较高，龚小玲等[39]对澳洲鳗鲡 (Anguilla australis) 进行标记开发时发现，微卫星重复序列的重复次数过高会影响PCR效果，应选择居中的重复次数为宜。长体圆鲹二核苷酸重复PIC为0.827 4，具有较高多态性，表明选择6~15次的二核苷酸重复是合适的。
3.3 微卫星标记的种群遗传学特征
群体杂合度的高低反映了群体在多个基因座上的遗传变异及群体遗传多样性丰富度[19]。本研究中长体圆鲹中沙群体的平均Ho为0.631 7，平均He为0.796 8，说明长体圆鲹该群体的遗传多样性较高。平均Ho和He存在差异，说明存在杂合子缺失或者纯合子过剩的情况。PIC也是衡量群体遗传多样性的重要指数，Botstein等[36]认为基因标记PIC>0.5为高度多态位点，0.25<PIC<0.5为中度多态位点，PIC<0.25为低度多态性位点，通常不作为遗传多样性分析。本文中长体圆鲹位点除1个为中度多态外，其他位点均为高度多态位点。表明开发所得的长体圆鲹微卫星标记在中沙群体中具有较好的遗传稳定性和丰富的遗传多样性。
在所有27个位点中有8个位点偏离了“哈迪-温伯格”平衡，这些位点不适合进一步的遗传分析。近亲杂交、无效等位基因、种群退化和自然选择等因素皆可能导致微卫星位点偏离HWE[15]。
-
表 1 长体圆鲹基因组中不同类型SSR统计
Table 1 Different types of SSR statistics in D. macrosoma genome
重复单元
repeat unit微卫星数量/个
number of microsatellite占比/%
ratio一核苷酸 mono-nucleotide 8 184 11.61 二核苷酸 di-nucleotide 37 646 53.39 三核苷酸 tri-nucleotide 13 960 19.80 四核苷酸 tetra-nucleotide 7 741 10.98 五核苷酸 penta-nucleotide 2 255 3.20 六核苷酸 hexa-nucleotide 722 1.02 合计 total 70 508 100.00 
下载: 导出CSV
表 2 27对长体圆鲹微卫星引物信息
Table 2 Information of 27 pairs of primers in D. macrosoma
位点
locus引物序列 (5'−3')
primer sequence重复单元
repeat motif退火温度/℃
annealing temperature期望长度/bp
allele sizeDma03 F:CCACGCCTATTGAGTTACAGA (CA)9 60 186 R:GAGCCAGTGGATGAACAGAGT Dma07 F:GCCCCTGTGGGTGTGTGA (CA)9 60 225 R:GGGTGGTGGGTTCGGTTT Dma12 F:GAACCAGTGCCTACAATAGA (AC)9 60 243 R:CTGCTCACGGTAAGTCCA Dma15 F:ACAGGAAGGAACAGGACAG (TG)10 55 254 R:TATTGAAGTGAAAAAGCCG Dma22 F:CGCTGTTGAAATGAAGAAGA (GT)10 60 317 R:AGTGATGTCGCCTCATAAAT Dma23 F:AAACTGAGGGCGAGATAGAGG (AC)10 55 190 R:CCACAGGCTGAGTAAACCAAC Dma26 F:ATCCCATTCACCGACATAG (TG)10 58 258 R:CTGTGGTATCGTTCCCTGT Dma28 F:TGATTGGCTTCTACTCTGC (AC)10 55 281 R:AGTGGCTTGTTTGACTCTTAT Dma36 F:GGATGTAGTGAAGAGGGGAG (GT)11 55 239 R:CACAATCAGTGTTATGGCAG Dma38 F:GCCAATAAAGGCAAACAGT (CA)11 60 227 R:ATCCGAGACAAAGACATACAA Dma39 F:AGTGTGCTGACTTTTCTCTG (CA)11 55 241 R:TTATTGTTTGTTGTCTGGGT Dma45 F:CTCCTTTTTCTTCTTCCTCT (CA)11 60 281 R:CTACCTGCTCTTCAACTCAT Dma51 F:TGACAGCCTCCACTACTCC (GA)12 55 225 R:GCTAACCAGACACGCAAA Dma54 F:AAAGCCCATCTGTCTCGT (GT)12 60 202 R:TGTTTCAGTCCGTTCCTG Dma58 F:TCAAGAGGGAGTGGGAGC (AC)12 58 279 R:TCAAATGGGTGTTTAGCG Dma64 F:GCTCAGACTGCGTGGACA (TG)13 55 314 R:GCTGGTGAACAACAGGACA Dma72 F:TTCCGCAGGCATAAAAAC (CT)13 58 301 R:CCAAGGTCCGCTACACTA Dma76 F:TTCTCGCTGACCTGCTTG (TG)14 55 253 R:GCGTCCTCGTCGTCTTTC Dma81 F:GAGACACGGTCAGAAAACA (TGC)6 60 216 R:GGAAGTAGGACTCTAGGGG Dma82 F:CTGTCACTCCATTCCTATTCC (GTT)6 58 244 R:CCTACATTTGTGCTTTTGTTC Dma83 F:CTCTAAAGCCGACCTAACC (CTT)6 58 239 R:TGTCTCAACACAGCGAAAC Dma84 F:AAACTAACTCATCACCAG (TGT)6 55 283 R:AAACGACAGGAACTCAAT Dma85 F:CTCACTTTGACCCAACCAG (AGG)6 55 256 R:CCTTTCACCGAGACACCAG Dma131 F:TGCGGATGGGTGGTAGTGT (GGT)8 55 208 R:ATTGCTGGTAGTCGGTGGC Dma132 F:CCCAGTGAGACCAGAACCA (GCT)8 55 268 R:GACCCGTAGACAGGAGAGT Dma135 F:GTTGTTGTTTTTTTCCTT (GCA)9 55 301 R:CATCAGTCTGGCTTTATA Dma145 F:ACGATACAGCAGCCGAAG (TCA)10 60 197 R:AGTGATGTCGCCTCATAAAT
下载: 导出CSV
表 3 长体圆鲹微卫星标记的种群遗传学特征
Table 3 Characteristics of microsatellite loci in D. macrosoma
位点
locusN Na Ho He PHWE PIC Dma03 35 8 0.857 1 0.790 9 0.042 0 0.746 0 Dma07 35 9 0.771 4 0.855 9 0.553 4 0.825 0 Dma12 34 11 0.685 7 0.816 6 0.106 4 0.809 0 Dma15* 35 14 0.542 9 0.911 8 0.000 0 0.890 0 Dma22* 32 14 0.485 7 0.813 2 0.000 0 0.864 0 Dma23 34 17 0.714 3 0.864 8 0.003 3 0.868 0 Dma26* 34 14 0.428 6 0.869 1 0.000 0 0.870 0 Dma28 35 11 0.742 9 0.837 7 0.035 3 0.803 0 Dma36 35 9 0.685 7 0.786 3 0.054 5 0.743 0 Dma38* 34 12 0.628 6 0.851 2 0.000 0 0.849 0 Dma39* 33 12 0.514 3 0.782 0 0.000 0 0.793 0 Dma45 34 12 0.771 4 0.847 8 0.139 2 0.844 0 Dma51* 31 12 0.485 7 0.784 5 0.000 0 0.859 0 Dma54 35 15 0.771 4 0.900 6 0.139 0 0.878 0 Dma58 34 12 0.771 4 0.845 6 0.449 7 0.843 0 Dma64* 31 13 0.428 6 0.788 3 0.000 0 0.865 0 Dma72 35 8 0.628 6 0.713 9 0.155 7 0.659 0 Dma76 35 12 0.685 7 0.908 1 0.003 6 0.886 0 Dma81 35 8 0.628 6 0.795 0 0.045 7 0.752 0 Dma82 35 5 0.485 7 0.538 3 0.023 6 0.497 0 Dma83* 34 8 0.342 9 0.705 8 0.000 0 0.674 0 Dma84 35 7 0.542 9 0.704 8 0.005 9 0.649 0 Dma85 35 8 0.828 6 0.747 8 0.876 9 0.704 0 Dma131 35 9 0.628 6 0.717 6 0.646 6 0.661 0 Dma132 35 8 0.657 1 0.746 6 0.365 5 0.692 0 Dma135 34 10 0.742 9 0.840 5 0.241 8 0.835 0 Dma145 34 7 0.600 0 0.749 3 0.077 4 0.725 0 注：N. 有效样品数；Na. 等位基因数；Ho. 表观杂合度；He. 期望杂合度；PHWE. “哈迪-温伯格”平衡显著性检验P值；PIC. 多态信息含量；*. 经Bonferroni校正后显著背离“哈迪-温伯格”平衡 (校正P<0.001 85) Note: N. effective number of samples; Na. number of alleles; Ho. observed heterozygosity; He. expected heterozygosity; PHWE. Hardy–Weinberg probability test; PIC. polymorphism information content; *. significant deviation from HWE after Bonferroni's correction (adjusted P-value<0.001 85)
下载: 导出CSV
表 4 长体圆鲹二、三核苷酸重复微卫星标记的比较
Table 4 Comparison on di- and trinucleotide-repeated microsatellite loci in D. macrosoma
序列
sequence引物数
primer number重复次数
repeat number筛选效率
efficiencyPIC 二核苷酸重复
di-nucleotide-repeated18 9~14 22.2% 0.827 4 三核苷酸重复
tri-nucleotide-repeated9 6~10 29.0% 0.687 7
下载: 导出CSV
-
[1] 中国科学院动物研究所, 中国科学院海洋研究所, 上海水产学院. 南海鱼类志[M]. 北京: 科学出版社, 1962: 387-389. [2] 杨吝, 张旭丰, 谭永光, 等. 南海北部灯光围网渔获组成分析[J]. 南方水产, 2009, 5(6): 65-70. doi: 10.3969/j.issn.1673-2227.2009.06.012 [3] 张立, 李渊, 林龙山, 等. 南海中南部主要经济种类渔业资源声学评估[J]. 海洋渔业, 2016, 38(6): 577-587. doi: 10.3969/j.issn.1004-2490.2016.06.003 [4] 张俊, 张鹏, 陈作志, 等. 南海外海鲹科鱼类资源量及其分布[J]. 南方水产科学, 2016, 12(4): 38-48. doi: 10.3969/j.issn.2095-0780.2016.04.005 [5] 张鹏, 张俊, 李渊, 等. 秋季南海中南部海域的一次灯光罩网探捕调查[J]. 南方水产科学, 2016, 12(2): 67-74. doi: 10.3969/j.issn.2095-0780.2016.02.010 [6] SHIRAISHI T, TANAKA H, OHSHIMO S, et al. Age, growth and reproduction of two species of scad, Decapterus macrosoma and D. macarellus in the waters off southern Kyushu[J]. Jpn Agr Res Q, 2010, 44(2): 197-206. doi: 10.6090/jarq.44.197
[7] BALASUBRAMANIAN N K, NATARAJAN P. Resource characteristics of the scads, Decapterus russelli and D. macrosoma of the Vizhinjam area, southwest coast of India[J]. Indian J Fish, 1999, 46(2): 111-122.
[8] POTIE M, DRAPEAU L. Modelling and forecasting the catch of the scads (Decapterus macrosoma, Decapterus russellii) in the Javanese purse seine fishery using ARIMA time series models[J]. Asian Fish Sci, 2000, 13: 75-85.
[9] 翟云, 吴仁协, 牛素芳, 等. 基于SLAF-seq技术开发蓝圆鲹微卫星标记及跨物种扩增检测[J]. 应用海洋学学报, 2018, 37(3): 426-434. doi: 10.3969/J.ISSN.2095-4972.2018.03.016 [10] DEWOODY J A, AVISE J C. Microsatellite variation in marine, freshwater and anadromous fishes compared with other animals[J]. J Fish Biol, 2000, 56(3): 461-473. doi: 10.1111/jfb.2000.56.issue-3
[11] EUSTICE M, YU Q Y, LAI C W. Development and application of microsatellite markers for genomic analysis of papaya[J]. Tree Genet Genome, 2008, 4(2): 333-341. doi: 10.1007/s11295-007-0112-2
[12] VARELA A I, RITCHIE P A, SMITH P J. Global genetic population structure in the commercially exploited deep-sea teleost orange roughy (Hoplostethus atlanticus) based on microsatellite DNA analyses[J]. Fish Res, 2013, 140: 83-90. doi: 10.1016/j.fishres.2012.12.011
[13] 海萨, 李家乐, 郭焱, 等. 河鲈微卫星引物筛选[J]. 水生态学杂志, 2008, 1(6): 90-94. [14] 赵哲霞, 蒋珊, 王滨花, 等. 黄颡鱼属SSR分子鉴定及其遗传多样性[J]. 南昌大学学报(理科版), 2014, 38(5): 498-501. doi: 10.3969/j.issn.1006-0464.2014.05.019 [15] 翟云, 吴仁协, 牛素芳, 等. 采用高通量技术开发花鲈二碱基重复微卫星标记[J]. 基因组学与应用生物学, 2018(9): 1-11. [16] CASTOE T A, STREICHER J W, MEIK J M, et al. Thousands of microsatellite loci from the venomous coralsnake Micrurus fulvius and variability of select loci across populations and related species[J]. Mol Ecol Resour, 2012, 12(6): 1105-1113. doi: 10.1111/1755-0998.12000
[17] 鲁翠云, 毛瑞鑫, 李鸥, 等. 鲤鱼三、四核苷酸重复微卫星座位的筛选及特征分析[J]. 农业生物技术学报, 2009, 17(6): 979-987. doi: 10.3969/j.issn.1674-7968.2009.06.007 [18] 谭照君, 张天奇, 鲁翠云, 等. 鲢三、四核苷酸重复微卫星标记的筛选及其特征分析[J]. 上海海洋大学学报, 2011, 20(3): 328-335. [19] 李文升, 刘翠, 鲁翠云, 等. 草鱼三、四核苷酸重复微卫星标记的分离与特征分析[J]. 中国水产科学, 2011, 18(4): 742-750. [20] THIEL T, MICHALEK W, VARSHNEY R K, et al. Exploiting EST databases for the development and characterization of gene-derived SSR-markers in barley (Hordeum vulgare L.)[J]. Theor Appl Genet, 2003, 106(3): 411-422. doi: 10.1007/s00122-002-1031-0
[21] SCHUELKE M. An economic method for the fluorescent labeling of PCR fragments[J]. Nat Biotechnol, 2000, 18(2): 233-234. doi: 10.1038/72708
[22] ROUSSET F. genepop'007: a complete re-implementation of the genepop software for Windows and Linux[J]. Mol Ecol Resour, 2008, 8(1): 103-106. doi: 10.1111/j.1471-8286.2007.01931.x
[23] KALINOWSKI S T, TAPER M L, MARSHALL T C. Revising how the computer program CERVUS accommodates genotyping error increases success in paternity assignment[J]. Mol Ecol, 2007, 16(5): 1099-1106. doi: 10.1111/j.1365-294X.2007.03089.x
[24] 杨兵, 林琳, 李纯厚, 等. 基于高通量测序的二长棘鲷微卫星标记开发与评价[J]. 南方水产科学, 2015, 11(4): 116-120. doi: 10.3969/j.issn.2095-0780.2015.04.017 [25] 秦海峰, 龙宁, 吴建国, 等. 甜叶菊微卫星富集文库的构建与多态性标记的筛选[J]. 作物学报, 2014, 40(3): 447-456. [26] LIN L, LI C H, CHEN Z Z, et al. Development and characterization of twenty-three microsatellite markers for the purpleback flying squid (Symplectoteuthis oualaniensis)[J]. Conserv Genet Resour, 2015, 7(1): 161-163. doi: 10.1007/s12686-014-0318-1
[27] KONG X L, CHEN Z Z, LIN L, et al. Polymorphic micro-satellite loci isolated from the yellowbelly threadfin bream, Nemipterus bathybius[J]. Genet Mol Res, 2014, 13(3): 5254-5257. doi: 10.4238/2014.July.24.3
[28] 房祖业, 陈晓东, 吴咏诗, 等. 大刺鳅 (Mastacembelus armatus) 二、三、四碱基重复微卫星标记的筛选和特征分析[J]. 海洋与湖沼, 2018, 49(1): 174-182. [29] 高峰涛, 邵长伟, 崔忠凯, 等. 基于高通量测序的青石斑鱼基因组微卫星开发及评价[J]. 中国海洋大学学报(自然科学版), 2017, 47(4): 52-57. [30] SONG W T, PANG R Y, NIU Y Z, et al. Construction of high-density genetic linkage maps and mapping of growth-related quantitative trail loci in the Japanese flounder (Paralichthys olivaceus)[J]. PLoS One, 2012, 7(11): 50404. doi: 10.1371/journal.pone.0050404
[31] SONG W, LI Y, ZHAO Y, et al. Construction of high-density microsatellite genetic linkage maps and mapping of sexual and growthrelated traits in half-smooth tongue sole (Cynoglossus semilaevis)[J]. PLoS One, 2012, 7(12): e52097. doi: 10.1371/journal.pone.0052097
[32] 曾晓芸, 杨宗英, 田辉伍, 等. 基于Mi-Seq高通量测序分析裸体异鳔鳅鮀微卫星组成[J]. 淡水渔业, 2015, 45(1): 3-7. doi: 10.3969/j.issn.1000-6907.2015.01.001 [33] 熊良伟, 王帅兵, 封琦, 等. 基于高通量测序的中华鳑鲏基因组微卫星特征分析及标记开发[J]. 江苏农业科学, 2018, 46(18): 164-168. [34] XIONG L W, WANG Q, QIU G F. Large-scale isolation of microsatellites from Chinese mitten crab Eriocheir sinensis via a solexa genomic survey[J]. Int J Mol Sci, 2012, 13(12): 16333-16345. doi: 10.3390/ijms131216333
[35] 高焕, 孔杰. 串联重复序列的物种差异及其生物功能[J]. 动物学研究, 2005, 26(5): 555-564. doi: 10.3321/j.issn:0254-5853.2005.05.017 [36] BOTSTEIN D, WHITE R L, SKOLNICK M, et al. Construction of a genetic linkage map in man using restriction fragment length polymorphisms[J]. Am J Hum Genet, 1980, 32(3): 314-331.
[37] ELLEGREN H. Microsatellite evolution: a battle between replication slippage and point mutation[J]. Trends Genet, 2002, 18(2): 70.
[38] WEBER J L. Informativeness of human (dC-dA) n·(dG-dT) n polymorphisms[J]. Genomics, 1990, 7(4): 524-530. doi: 10.1016/0888-7543(90)90195-Z
[39] 龚小玲, 李思发, 蔡完其, 等. 澳洲鳗鲡微卫星分子标记的筛选与检测[J]. 中国水产科学, 2009, 16(1): 133-138. doi: 10.3321/j.issn:1005-8737.2009.01.019 -
期刊类型引用(6)
1. 唐棣，杨明秋，邹雄，刘洪涛. 基于锈斑蟳简化基因组数据的微卫星DNA标记开发. 广西科学. 2024(04): 763-772 . 
百度学术
2. 刘峰林，王泽松，王太，邓嘉俊，王小东，金燕君，宋昭彬. 基于高通量测序的厚唇裸重唇鱼微卫星分子标记筛选. 动物学杂志. 2023(06): 817-826 . 百度学术
3. 张芹，王延晖，王冰柯，屈长义. 团头鲂基因组微卫星特征及SSR位点开发. 江苏农业科学. 2022(03): 79-85 . 百度学术
4. 上官清，陈昆慈，刘海洋，欧密，罗青，王亚坤，徐晟云，赵建. 斑鳢基因组中微卫星分布特征及野生种群遗传结构分析. 南方水产科学. 2020(03): 47-60 . 本站查看
5. 张永德，文露婷，罗洪林，林勇，杜雪松，余艳玲，韦孜娜，黄姻. 卵形鲳鲹基因组调研及其SSR分子标记的开发应用. 南方农业学报. 2020(05): 983-994 . 百度学术
6. 王锦秀，宋勇，王新月，陈生熬，任道全. 基于de novo高通量测序的叶尔羌高原鳅微卫星位点筛选与多态性分析. 塔里木大学学报. 2020(02): 25-32 . 百度学术
其他类型引用(7)
计量
- 文章访问数: 5533
- HTML全文浏览量: 2028
- PDF下载量: 39
- 被引次数: 13
粤公网安备 44010502001741号
