2017, Vol. 48
中国海洋湖沼学会主办。
文章信息
- 吕振明, 侯龙, 龚理, 刘立芹, 陈永久, 郭宝英, 董迎辉, 吴常文. 2017.
- LÜ Zhen-Ming, HOU Long, GONG Li, LIU Li-Qin, CHEN Yong-Jiu, GUO Bao-Ying, DONG Ying-Hui, WU Chang-Wen. 2017.
- 基于de novo高通量测序的曼氏无针乌贼(Sepiella japonica)ESTs中微卫星位点筛选与特征分析
- ISOLATION AND ANALYSIS ON EST MICROSATELLITES OF SEPIELLA JAPONICA BY DE NOVO HIGH-THROUGHPUT TRANSCRIPTOME SEQUENCING
- 海洋与湖沼, 48(4): 877-883
- Oceanologia et Limnologia Sinica, 48(4): 877-883.
- http://dx.doi.org/10.11693/hyhz20170300079
-
文章历史
- 收稿日期：2017-03-29
- 收修改稿日期：2017-04-23
2. 浙江万里学院 浙江省水产种质资源高效利用技术研究重点实验室 宁波 315100
2. Zhejiang Key Laboratory of Aquatic Germplasm Resources, Zhejiang Wanli University, Ningbo 315100, China
曼氏无针乌贼(Sepiella japonica)属软体动物门、头足纲、十腕目、乌贼科、无针乌贼属种类, 是我国重要的海洋渔业资源之一, 20世纪70年代, 我国曼氏无针乌贼年产量曾达6万多吨, 是我国著名的“四大海产”之一(董正之, 1988)。但自20世纪80年代以来, 由于过渡捕捞或海域环境的破坏, 其资源量严重衰退, 几近枯竭的边缘(Lü et al, 2016)。为了恢复我国曼氏无针乌贼的渔业资源和产量, 近年来我国陆续攻克了其规模化苗种繁育和养殖技术, 人工养殖已在国内沿海省份悄然兴起, 而其人工增殖放流更是在国内得到迅猛开展, 并取得显著成效(Jiang et al, 2014; Yan et al, 2016)。然而, 为更好地开展曼氏无针乌贼的繁养及增殖养护工作, 迫切需要对其现有种群的种质状况和遗传组成进行精确的了解。
微卫星标记技术是研究物种种质状况和遗传组成的重要方法, 由于其高度多态性和共显性等优势而被广泛运用于物种的群体遗传结构(An et al, 2014)、种质状况评估(Lü et al, 2013)及遗传改良研究中(Fuji et al, 2007)。近年来, 有关曼氏无针乌贼微卫星标记的开发已有少量相关报道(Wu et al, 2010; Guo et al, 2013), 但仍需要开发出更多的微卫星标记用于其种质和遗传组成的评估, 同时以往的研究均以DNA建库和探针杂交富集法进行微卫星的筛选(Zane et al, 2002), 但该方法不仅费时、费力, 而且需要预先对物种的基因组信息有较好的了解(Wang et al, 2012)。近年来, 基于二代高通量测序的出现为快速、高效的进行微卫星的开发提供了可能(Csencsics et al, 2010), 该技术可通过一次高通量的测序产生海量的微卫星序列, 因而被广泛运用于许多非模式动植物的微卫星开发中(Yu et al, 2011; Zhu et al, 2012)。特别是转录组高通量测序技术, 由于成本低, 微卫星序列检出率高, 因而被广泛运用于包括鱼类(Teacher et al, 2012)、虾蟹类(Ma et al, 2014)、贝类(Patnaik et al, 2016)等水生生物的微卫星标记开发中, 目前采用高通量测序进行乌贼等头足类动物的微卫星标记开发还未见相关报道, 本研究拟采用转录组测序技术从海量的EST序列中筛选微卫星序列, 分析其序列组成和特征, 开发多态性较强的微卫星位点用于其种群遗传变异分析, 相关研究成果将为今后我国曼氏无针乌贼资源科学合理保护和开发奠定基础。
1 材料与方法 1.1 材料用于转录组de novo高通量测序分析的曼氏无针乌贼(Sepiella japonica)鲜活样品取自浙江省海洋水产研究所西闪养殖基地。活体解剖后取肌肉、肝脏、鳃、脑、性腺等组织, 迅速放入液氮中冷冻保存备用。用于多态微卫星位点分离与鉴定及群体遗传多态度检测的曼氏无针乌贼样本取自福建漳州海域的自然群体, 共计35只, 取肌肉组织用70%的酒精保存备用。
1.2 方法 1.2.1 转录组EST微卫星序列获取曼氏无针乌贼组织RNA的提取、cDNA文库的构建参照Teacher等(2012)的方法进行。构建好的cDNA文库在CBot上进行簇生成反应(TruSeq PE Cluster Kit V3-cBot-HS, Illumina)后, 采用Illumina Hiseq 2000高通量测序平台(TruSeq SBS KIT-HS V3, Illumina)进行双末端测序, 由深圳华大基因科技有限公司完成。获得的clean reads进行De novo序列拼接和组装, 采用MISA软件对拼接得到的序列进行微卫星的搜索, SSR位点确认设置条件为单碱基类型重复至少12次, 2碱基类型重复至少6次, 3、4碱基类型重复至少5次, 6碱基以上类型重复至少4次。对转录组EST微卫星的类型, 各微卫星类型比例及序列特征等进行分析和统计。
1.2.2 多态微卫星位点分离与鉴定选取二、三碱基重复单元为主的微卫星序列, 采用primer3设计微卫星引物, 引物合成由上海英俊生物技术有限公司完成。微卫星位点分离时DNA的提取采用传统的酚、氯仿抽提法; PCR扩增采用25μL的反应体系, 内含50ng基因组DNA, 1×buffer (Promega, USA), 2mmol/L MgCl2, 0.2mmol/L dNTPs, 0.2μmol/L各引物和1U Taq酶(Promega, USA)。反应程序条件为94℃预变性5min; 94℃变性1min, 50—65℃退火1min, 72℃延伸1min, 35个循环; 72℃延伸10min, 各引物退火条件在50—65℃范围内进行优化。利用10个个体进行多态性标记的检测, PCR产物采用毛细管电泳检测, 用Genemapper 3.5 (Applied Biosystem)分析产物长度。
1.2.3 群体遗传多态度检测与统计采用筛选出的多态性微卫星引物对采自福建漳州天然海域的35只曼氏无针乌贼进行SSR扩增, 产物经电泳、分型后, 采用ARLEQUIN3.5.1.3 (Excoffier et al, 2010)软件统计位点等位基因数(NA)、观测杂合度(Ho)、预期杂合度(He)、多态信息含量(PIC)等遗传变异参数。采用CERVUS3.0.3 (Marshall et al, 1998)和GENEPOP 4.0.10 (Raymond et al, 1995)检测各位点等位基因分布是否符合Hardy-Weinberg平衡(HWE)及各位点是否存在连锁不平衡。
2 结果 2.1 曼氏无针乌贼转录组EST中微卫星的序列特征采用Illumina测序和De novo拼接后共获得了85674个Unigene, 总长66760294nt, 共搜索到47604个微卫星序列, 平均每1402nt的EST序列中就含有1个至少12nt的微卫星序列, 微卫星分布相当丰富。对微卫星的组成类型和数量进行分析, 结果表明:单碱基重复是曼氏无针乌贼转录组EST最主要的微卫星类型(38.69%); 其次为三碱基(31.14%)和二碱基(26.35%)重复类型, 上述三种短序列重复类型占微卫星总量的96.18%。而长序列类型微卫星相对占比很低, 四碱基(3.29%)、五碱基(0.38%)、六碱基(0.14%)的重复类型总和只占微卫星总量的3.82%, 如图 1所示。各碱基类型的微卫星中, 不同序列组成的微卫星数量迥异。如图 2所示, 二碱基类型中, AC (54.51%)和AG (31.22%)是最常见的二碱基重复形式, 两者在二碱基类型中的占比达85.73%;而CG (0.10%)的重复形式最少, 较低频率的GC类型可能与胞嘧啶易甲基化, 使其因脱氨基作用而变成胸腺嘧啶有关(Schorderet et al, 1992)。而三碱基类型中各碱基组成形式差异相对较小, AGC (16.37%)和AAC (14.06%)是最常见的三碱基重复形式, 其次为ATC (13.72%)、ACC (13.07%)、AAT (10.67%)、ACT (10.27%)、AGG (9.62%)、AAG (9.52%), 而ACG (1.99%)和CCG (0.72%)类型则最少。
 |
| 图 1 曼氏无针乌贼转录组EST中1-6碱基类型微卫星的组成 Fig. 1 The percentages of mono-, di-, tri-, tetra-, quad-, penta-, and hexa-nucleotide repeats in SSR motif sequences in S. japonica transcriptome EST |
 |
| 图 2 曼氏无针乌贼转录组EST中二、三碱基类型微卫星的序列组成情况 Fig. 2 Observed counts of identified microsatellite loci for different repeat sequence motifs of di-and trinucleotide repeats in S. japonica transcriptome EST |
各碱基类型的微卫星在重复次数变异上也存在巨大差别。如表 1所示, 二碱基类型中重复次数变异6—11次不等, 而三碱基类型的重复次数变异范围为5—8次; 四碱基和五碱基类型的重复次数仅5—6次和4—5次; 六碱基类型中未发现有重复次数的变异, 均为重复4次, 表明各碱基类型随着碱基数的增加存在的变异也逐渐减少。同时, 转录组中微卫星均以最低重复数的序列最多, 如在二碱基类型中, 最低重复数6的微卫星数量最多, 占所有二碱基类型的35.04%。三碱基类型中, 最低重复数5的微卫星数量最多, 占所有三碱基类型的53.85%。随着重复数的增加, 相应微卫星的数量逐级减少。
| 重复数 | 二碱基 | 三碱基 | 四碱基 | 五碱基 | 六碱基 |
| 4 | 0 | 0 | 0 | 85 | 33 |
| 5 | 0 | 3948 | 713 | 4 | 0 |
| 6 | 2174 | 1934 | 61 | 0 | 0 |
| 7 | 1292 | 1281 | 0 | 0 | 0 |
| 8 | 921 | 169 | 0 | 0 | 0 |
| 9 | 866 | 0 | 0 | 0 | 0 |
| 10 | 665 | 0 | 0 | 0 | 0 |
| 11 | 278 | 0 | 0 | 0 | 0 |
| 12 | 8 | 0 | 0 | 0 | 0 |
| 总数 | 6204 | 7332 | 774 | 89 | 33 |
共筛选了121对二、三碱基类型为主的微卫星引物, 仅65对引物能稳定扩增出清晰、可辨的微卫星条带, 其中具多态位点的微卫星数达24对。采用24对微卫星多态位点对采自福建漳州海域的35只曼氏无针乌贼个体进行遗传变异和多态性检测, 结果如表 2所示: 24对引物共扩增得到142个等位基因, 每个位点的等位基因数(Na)3—10不等, 平均5.92个等位基因。杂合度观测值(Ho)在0.030—0.857之间, 均值为0.518, 杂合度预期值(He)在0.420—0.826之间, 均值为0.681。多态信息含量(PIC)在0.3863—0.8728之间。除位点CL2235-2 (PIC=0.3863)、CL899(0.4009) 外, 所有位点的PIC值都大于0.5, 显示了较高的多态性。二、三碱基类型微卫星相比, 二碱基微卫星标记的多态性无论是位点等位基因数(Na)、杂合度(He)、多态信息含量(PIC)等参数都明显比三碱基高, 并且在He和PIC上的差异达到显著(P < 0.05)。不同重复次数的微卫星相比, 重复次数越高的微卫星其多态性也越高, 但仅在位点等位基因数(Na)差异上达到显著水平(P < 0.05), 如表 3所示。经检测, 有4个位点明显偏离Hardy-Weinberg平衡(P < 0.05);经连锁不平衡检测, 所有位点均未检测到明点的连锁不平衡现象(P > 0.05)。
| 位点 | 重复序列 | 引物序列(5′—3′) | Ta(℃) | Na (范围) | Ho | He | PHWE | PIC | GenBank登录号 |
| CL168 | (AAC)6 | F: ACAATCAACGGCTGTAAAGTCA R: GACTATGGTTTGGATTTGGCAT | 55 | 3(172-178) | 0.647 | 0.599 | 0.633 | 0.5219 | KU306816 |
| CL218 | (TA)6 | F: CCTGGTTCACAACAACATTTTC R: CAAAGACAGGACACGACCAATA | 55 | 5(194-208) | 0.030 | 0.757 | 0.000* | 0.7128 | KU306820 |
| CL2235-2 | (CAC)5 | F: CGAGATTCTGTACAAGGCTC R: GGAATTGTGGCAGAAAGACT | 58 | 5(201-216) | 0.382 | 0.620 | 0.908 | 0.3863 | KU306831 |
| CL227 | (TTC)5 | F: TGTCATTTCTGATTCGTCATCC R: GATCCAGCAACAATGGACTTG | 55 | 5(211-229) | 0.686 | 0.620 | 0.765 | 0.5656 | KU306817 |
| CL259 | (CT)7 | F: CTGTGTTTGAGCACGAATTAGG R: CAACAGGTAGGTCTTGTGCAAC | 55 | 7(193-207) | 0.771 | 0.826 | 0.8025 | 0.8189 | KU306827 |
| CL2635 | (TG)8 | F: TAGACTTGTATGGGTGGGTG R: TAAATCCACTGCCTGCTTTC | 61 | 6(202-220) | 0.143 | 0.591 | 0.000* | 0.5217 | KU306829 |
| CL334 | (AC)6 | F: CTCACTGGCACTTCTTGTCTTA R: ATGTACTTAGAGGAGGAAACGTATG | 55 | 4(198-208) | 0.857 | 0.693 | 0.069 | 0.6317 | KU306821 |
| CL3354 | (CTG)5… (TGC)5 | F: CCTCGGCTTCTGATGAAAAT R: AGCCTTACTTCTGCAACATG | 55 | 6(213-228) | 0.485 | 0.699 | 0.071 | 0.6550 | KU306828 |
| CL402 | (CTG)6 | F: CATGGGTGACAGAGGACATAAA R: GGCACGAACAGAACATGAAAT | 55 | 5(218-230) | 0.571 | 0.602 | 0.813 | 0.5450 | KU306818 |
| CL416 | (GT)11 | F: CTCACGTACCAAGACTCCATACC R: TTGCACCAGAATCAAATGAGAG | 55 | 8(162-182) | 0.486 | 0.767 | 0.000* | 0.7347 | KU306822 |
| CL5062-1 | (TGT)7 | F: GGTGGTGAGCTAACATGTTT R: CAACGGAATATTGCCGAAGA | 55 | 6(174-189) | 0.457 | 0.579 | 0.169 | 0.5448 | KU306832 |
| CL5062-2 | (ATG)6… (TGC)5 | F: GCATCCCTTTGTCTGTCTTT R: GTCAATTCGGTGGCCTATAC | 52 | 5(203-215) | 0.500 | 0.626 | 0.050 | 0.5848 | KU306832 |
| CL529 | (GA)6 | F: ATAGAACGTGGTGGGAGAGAAA R: TCCTTTTGTGTGTTTGTGGTTC | 56 | 7(192-214) | 0.429 | 0.660 | 0.643 | 0.6004 | KU306819 |
| CL662-2 | (GTA)6 | F: GCTGTTTAGGAGTAGCACTG R: CGGCTACTCATCAAAGACTG | 55 | 7(194-212) | 0.543 | 0.722 | 0.242 | 0.6782 | KU306830 |
| CL885 | (GA)8 | F: CCTTGTACAACTGCCACATT R: CAGCTGCCCAAAATATCTGT | 55 | 8(196-214) | 0.514 | 0.579 | 0.246 | 0.5574 | KU306825 |
| CL899 | (AAT)6 | F: ATCCAAGTCTTCTGCCTGAT R: CCATGTAACGATTCTCCCAC | 55 | 6(188-203) | 0.382 | 0.420 | 0.112 | 0.4009 | KU306824 |
| CL904 | (AT)8 | F: TCTAGGCCTGTGGTTAATGT R: TGATCGTTACTTGATGGCAG | 55 | 8(222-236) | 0.486 | 0.752 | 0.058 | 0.7136 | KU306823 |
| CL933 | (CCT)6 | F: CTGGGACGGGCAAATAAATT R: GTCTGAGGTTAGAGTCTGCT | 60 | 8(182-206) | 0.600 | 0.702 | 0.051 | 0.6658 | KU306826 |
| CL38 | ACC*5 | F: CAAATTCCTCTGCATTGGGA R: GAGTAAAGAAGGCAGGGTTG | 55 | 4(154-163) | 0.571 | 0.689 | 0.089 | 0.6371 | KX839253 |
| CL1360 | CT*8 | F: TCTCTCCCACACCAATATCC R: GTTACCGTTCTGATGTCTGG | 55 | 10(190-208) | 0.314 | 0.884 | 0.000* | 0.8728 | KX839254 |
| CL327 | TA*6 | F: ACAGCATCTTCTGGTAAGCCAT R: TAGTCCTGTCACCACAGTTATGC | 58 | 5(156-164) | 0.571 | 0.749 | 0.111 | 0.7050 | KX839255 |
| CL710 | CTC*5 | F: GGTCGTTATCCATTTGGTCATT R:AAAGAAGAAGGAGAGGAGGAGG | 55 | 5(184-198) | 0.857 | 0.738 | 0.085 | 0.6945 | KX839256 |
| CL7421 | CA*7…CA*7 | F:TACACGTCCTAGAGAGCCTA R: TAATACCGCCTCACTCTCTC | 55 | 4(180-186) | 0.543 | 0.729 | 0.061 | 0.6793 | KX839257 |
| CL357 | GGT*5 | F: CGTTCCTCATGTCTTCCTCTTC R: AAGCAACTGGAATTGGCTCTTA | 55 | 5(186-208) | 0.600 | 0.743 | 0.058 | 0.6971 | KX839258 |
| 注: Na为位点等位基因数; Ho为杂合度观测值; He为杂合度预期值；PHWE为Hardy-Weinberg平衡显著性P值; PIC为多态信息含量; Ta为PCR退火温度; *代表显著偏离Hardy-Weinberg平衡的位点 | |||||||||
| 碱基数 | Na | Ho | He | PIC |
| 二碱基 | 6.55 | 0.47 | 0.73a | 0.69a |
| 三碱基 | 5.38 | 0.56 | 0.64b | 0.58b |
| 重复数 | Na | Ho | He | PIC |
| 5次 | 5.00a | 0.60 | 0.68 | 0.61 |
| 6次 | 5.50b | 0.51 | 0.65 | 0.60 |
| 7次以上 | 7.13c | 0.46 | 0.71 | 0.68 |
| 注: Na为位点等位基因数; Ho为杂合度观测值; He为杂合度预期值; PIC为多态信息含量; 表中同列数值右上角字母不同的代表有显著性差异(P < 0.05) | ||||
从本研究的结果可以看出, 基于转录组的高通量测序是非常快速、高效的海洋头足类微卫星分离方法。本研究之前, 仅有约25对曼氏无针乌贼的微卫星标记被成功开发(Wu et al, 2010; Guo et al, 2013), 而利用de novo转录组高通量测序技术, 我们共获得了85674个Ungene, 基因序列总长66760294nt, 共搜索到47604个微卫星序列, 平均每1402nt的EST就有一个不少于12nt的微卫星序列。大部分序列都含有完整的适于微卫星引物设计的两翼序列, 并可用于今后的微卫星标记大规模开发中, 这些微卫星序列几乎是之前采用DNA建库和探针杂交富集法获得的微卫星序列的上千倍。这也许就是为什么近年来高通量测序技术越来越多地被用于微卫星标记开发之中的原因。
除快速、高效之外, 高通量测序技术还提供了传统方法无法提供的微卫星基因组序列特征, 而这些微卫星基因组序列信息可为更有效地筛选微卫星目标序列、设计引物, 更快速地筛选高多态性的微卫星标记提供了良好的基础。如由研究结果可知, 一、二、三碱基是最主要基因组微卫星碱基类型, 四、五、六碱基则相对稀少; 呈现出随着碱基数量的增加, 微卫星数量逐步减少的现象, 这与已知的许多真核生物的基因组微卫星组成情况类似(Chistiakov et al, 2006; Zhu et al, 2012), 这就提示我们在曼氏无针乌贼微卫星标记规模化开发时, 二、三碱基的微卫星类型可能是最主要的筛选对象。同时二、三碱基的微卫星类型相比, 二碱基微卫星的重复次数从6—12共7种类型, 而三碱基微卫星的重复次数仅5—8共4种类型, 二碱基比三碱基微卫星具有更多的变异类型, 暗示二碱基也可能比三碱基微卫星标记具有更高的多态性, 而该结果也得到了后续实验结果的证实, 二碱基微卫星无论在He和PIC上都显著高于三碱基(P < 0.05), 这也提示为了更有效地筛选多态性微卫星标记, 二碱基类型将是首选的筛选对象。另外, 不同重复次数的微卫星相比, 重复次数越高的微卫星其多态性也越高, 如图 2所示。该结果与前人在哺乳动物包括人类中的研究结果相似(Kelkar et al, 2008), 这也告诉我们, 在曼氏无针乌贼的微卫星开发中, 选择重复次数高的微卫星将是提高微卫星多态性的有效方法。
本研究共筛选出24对多态性的微卫星位点, 并用于检测来自福建漳州的曼氏无针乌贼自然群体的遗传多态性。结果表明, 漳州曼氏无针乌贼自然群体仍存在较高的遗传多样性。24个位点的等位基因数(NA)3—10不等, 平均5.92个等位基因。杂合度观测值(Ho)在0.030—0.857之间, 杂合度预期值(He)在0.420—0.826之间, 多态信息含量(PIC)在0.3863— 0.8728之间。该结果与Guo等(2013)以DNA测序为基础的生物素磁珠法开发的微卫星标记多态性大致相当, 而略低于Wu等(2010)以尼龙膜杂交法开发的微卫星标记的多态性。但应该指出的是Guo等(2013)研究的群体是福建宁德的自然群体, 与本研究的福建漳州自然群体所处海区接近, 而Wu等(2010)研究的群体为广东湛江群体。该结果说明, 采用转录组测序技术获得的微卫星多态性也可与DNA测序获得的微卫星多态性相当, 而与过去认为的编码区域的微卫星多态性一般较非编码区域的微卫星多态性低的看法相异(Mochmann et al, 2004; Pearson et al, 2005), 而与Kelkar等(2008)哺乳动物如人与大猩猩中的结果类似。因此该结果再一次说明转录组测序技术是一种高效快捷的微卫星标记开发方法。经检测24个微卫星位点无明显的连锁不平衡现象(P > 0.05), 仅4个显著偏离Hardy-Weinberg平衡(P < 0.05), 这些新开发的多态性微卫星位点可为今后曼氏无针乌贼渔业资源的科学管理和保护提供有力的遗传学评估手段。
4 结论本研究的结果表明, 基于转录组的高通量测序技术是海洋头足类微卫星标记高效开发的理想方法。与传统的DNA建库和探针杂交富集法开发微卫星标记相比, 其不仅快速、高效, 而且开发的微卫星多态性也与传统方法相当。本研究筛选出的24对多态性微卫星位点可用于今后乌贼资源遗传评估和科学有效管理过程。
| 董正之, 1988. 中国动物志:软体动物门, 头足纲. 北京: 科学出版社, |
| An H S, Lee J W, Kim H Y, et al, 2014. Population genetic structure of the sea bass (Lateolabrax japonicus) in Korea based on multiplex PCR assays with 12 polymorphic microsatellite markers. Genes & Genomics, 36(2): 247–259 |
| Chistiakov D A, Hellemans B, Volckaert F A M, et al, 2006. Microsatellites and their genomic distribution, evolution, function and applications:a review with special reference to fish genetics. Aquaculture, 255(1-4): 1–29 DOI:10.1016/j.aquaculture.2005.11.031 |
| Csencsics D, Brodbeck S, Holderegger R, 2010. Cost-effective, species-specific microsatellite development for the endangered dwarf bulrush (Typha minima) using next-generation sequencing technology. Journal of Heredity, 101(6): 789–793 DOI:10.1093/jhered/esq069 |
| Excoffier L, Lischer H E L, 2010. Arlequin suite ver 3.5:a new series of programs to perform population genetics analyses under Linux and Windows. Molecular Ecology Resources, 10(3): 564–567 |
| Fuji K, Hasegawa O, Honda K, et al, 2007. Marker-assisted breeding of a lymphocystis disease-resistant Japanese flounder (Paralichthys olivaceus). Aquaculture, 272(1-4): 291–295 DOI:10.1016/j.aquaculture.2007.07.210 |
| Guo B Y, Qi P Z, Zhu A Y, et al, 2013. Isolation and characterization of new polymorphic microsatellite markers from the cuttlefish Sepiella maindroni (Cephalopoda; Sepiidae). Genetics and Molecular Research, 12(3): 2376–2379 DOI:10.4238/2013.July.11.1 |
| Jiang L H, Zhu A Y, Wu C W, et al, 2014. Tetracycline immersion tagging of cuttlefish, Sepiella japonica, Larvae. Journal of the World Aquaculture Society, 45(3): 342–349 DOI:10.1111/jwas.2014.45.issue-3 |
| Kelkar Y D, Tyekucheva S, Chiaromonte F, et al, 2008. The genome-wide determinants of human and chimpanzee microsatellite evolution. Genome Research, 18(1): 30–38 |
| Lü Z M, Li H M, Liu L Q, et al, 2013. Rapid development of microsatellite markers from the large yellow croaker (Pseudosciaena crocea) using next generation DNA sequencing technology. Biochemical Systematics and Ecology, 51: 314–319 DOI:10.1016/j.bse.2013.09.019 |
| Lü Z M, Liu W, Liu L Q, et al, 2016. Cloning, characterization, and expression profile of estrogen receptor in common Chinese cuttlefish, Sepiella japonica. Journal of Experimental Zoology Part A, 325(3): 181–193 DOI:10.1002/jez.v325.3 |
| Ma H Y, Jiang W, Liu P, et al, 2014. Identification of transcriptome-derived microsatellite markers and their association with the growth performance of the mud crab (Scylla paramamosain). PLoS One, 9(2): e89134 DOI:10.1371/journal.pone.0089134 |
| Marshall T C, Slate J, Kruuk L E B, et al, 1998. Statistical confidence for likelihood-based paternity inference in natural populations. Molecular Ecology, 7(5): 639–655 DOI:10.1046/j.1365-294x.1998.00374.x |
| Mochmann L H, Wells R D, 2004. Transcription influences the types of deletion and expansion products in an orientation-dependent manner from GAC·GTC repeats. Nucleic Acids Research, 32(15): 4469–4479 DOI:10.1093/nar/gkh787 |
| Patnaik B B, Wang T H, Kang S W, et al, 2016. Sequencing, de novo assembly, and annotation of the transcriptome of the endangered freshwater pearl bivalve, Cristaria plicata, provides novel insights into functional genes and marker discovery. PLoS One, 11(2): e0148622 DOI:10.1371/journal.pone.0148622 |
| Pearson C E, Edamura K N, Cleary J D, et al, 2005. Repeat instability:mechanisms of dynamic mutations. Nature Reviews Genetics, 6(10): 729–742 DOI:10.1038/nrg1689 |
| Raymond M, Rousset F, 1995. An exact test for population differentiation. Evolution, 49(6): 1280–1283 DOI:10.1111/evo.1995.49.issue-6 |
| Schorderet D F, Gartler S M, 1992. Analysis of CpG suppression in methylated and nonmethylated species. Proceedings of the National Academy of Sciences of the United States of America, 89(3): 957–961 DOI:10.1073/pnas.89.3.957 |
| Teacher A G F, Khknen K, Meril J, 2012. Development of 61 new transcriptome-derived microsatellites for the Atlantic herring (Clupea harengus). Conservation Genetics Resources, 4(1): 71–74 DOI:10.1007/s12686-011-9477-5 |
| Wang J J, Yu X M, Zhao K, et al, 2012. Microsatellite development for an endangered bream megalobrama pellegrini (Teleostei, Cyprinidae) using 454 sequencing. International Journal of Molecular Sciences, 13(12): 3009–3021 DOI:10.3390/ijms13033009 |
| Wu C W, Chi C F, He G Y, et al, 2010. Isolation via enrichment and characterization of ten polymorphic microsatellite loci in the cuttlefish, Sepiella maindroni de Rochebruns. Acta Oceanologica Sinica, 29(6): 121–124 DOI:10.1007/s13131-010-0083-2 |
| Yan Y J, Wang T M, Liu W, et al, 2016. Identification and expression profile of the gonadotropin-releasing hormone receptor in common Chinese cuttlefish, Sepiella japonica. Journal of Experimental Zoology Part A, 325(7): 453–466 DOI:10.1002/jez.v325.7 |
| Yu J N, Won C, Jun J, et al, 2011. Fast and cost-effective mining of microsatellite markers using NGS technology:an example of a Korean water deer Hydropotes inermis argyropus. PLoS One, 6(11): e26933 DOI:10.1371/journal.pone.0026933 |
| Zane L, Bargelloni L, Patarnello T, 2002. Strategies for microsatellite isolation:a review. Molecular Ecology, 11(1): 1–16 DOI:10.1046/j.0962-1083.2001.01418.x |
| Zhu H, Senalik D, McCown B H, et al, 2012. Mining and validation of pyrosequenced simple sequence repeats (SSRs) from American cranberry (Vaccinium macrocarpon Ait. ). Theoretical and Applied Genetics, 124(1): 87–96 DOI:10.1007/s00122-011-1689-2 |