欢迎访问《热带作物学报》,
生物技术与组织培养

基于SLAF-seq技术的甘薯种质资源群体遗传进化分析

  • 李慧峰 ,
  • 黄咏梅 ,
  • 李彦青 ,
  • 滑金锋 ,
  • 吴翠荣 ,
  • 范继征 ,
  • 陈天渊 , *
展开
  • 广西农业科学院玉米研究所,广西南宁 530007
陈天渊,E-mail:

李慧峰(1981—),男,博士,副研究员,研究方向:甘薯遗传育种。

收稿日期: 2019-02-20

  要求修回日期: 2019-06-10

  网络出版日期: 2019-12-20

基金资助

广西农业科学院基本科研业务专项(桂农科2017YM12)

国家现代农业产业技术体系广西创新团队建设专项(No.nycytxgxcxtd-11-03)

现代农业产业技术体系建设专项(No.CARS-10-C19)

版权

版权所有,未经授权,不得转载、摘编本刊文章,不得使用本刊的版式设计。

Phylogenetic Analysis of Sweetpotato Germplasm Resources Based on SLAF-seq Technology

  • LI Huifeng ,
  • HUANG Yongmei ,
  • LI Yanqing ,
  • HUA Jinfeng ,
  • WU Cuirong ,
  • FAN Jizheng ,
  • CHEN Tianyuan , *
Expand
  • Maize Research Institute, Guangxi Academy of Agricultural Sciences, Nanning, Guangxi 530007, China
CHEN Tianyuan,E-mail:

Received date: 2019-02-20

  Request revised date: 2019-06-10

  Online published: 2019-12-20

Copyright

Copyright reserved © 2019. Office of Acta Agronomica Sinica All articles published represent the opinions of the authors, and do not reflect the official policy of the Chinese Medical Association or the Editorial Board, unless this is clearly specified.

摘要

基于简化基因组测序技术对122份甘薯种质资源进行分析,开发SNP位点并将其应用于种质资源群体结构和遗传进化分析。结果表明, 从122份甘薯种质资源中共获得563.18 Mb读长,不同材料的读长数量在1 419 809~8 392 785范围之间。测序质量值Q30在90.61%~96.82%之间,平均Q30为92.78%。测序获得的GC含量在36.46%~40.33%之间,平均GC含量为38.17%,对照GC含量为40.72%。根据测序结果共开发出高质量的SLAF标签2 388 759个,平均测序深度为17.45×。其中,多态性的SLAF标签77 761个,占SLAF标签总数的3.26%。根据所获得多态性SLAF标签统计SNP位点信息,共计获得129 063个群体SNP位点。遗传进化分析表明122份种质资源可以分为3个大类,分类结果与它们的地理来源无直接关系,聚类结果可为甘薯育种亲本组配和杂种优势利用提供科学依据。

本文引用格式

李慧峰 , 黄咏梅 , 李彦青 , 滑金锋 , 吴翠荣 , 范继征 , 陈天渊 . 基于SLAF-seq技术的甘薯种质资源群体遗传进化分析[J]. 热带作物学报, 2019 , 40(12) : 2390 -2396 . DOI: 10.3969/j.issn.1000-2561.2019.12.011

Abstract

Based on SLAF-seq technology, 122 sweet potato germplasm resources were used to develop SNP sites, which were applied to the analysis of population structure and genetic evolution on all germplasm resources. Results showed that 563.18 Mb reads length was obtained by sequencing, the reads of the samples varied from 1 419 809 to 8 392 785. The sequencing quality value (Q30) changed from 90.61% to 96.82%, and the average Q30 was 92.78%. The GC content of the samples changed from 36.46%-40.33%, the average value was 38.17%, and that of the control was 40.72%. A total of 2 388 759 SLAF tags were developed, with an average sequencing depth of 17.45×. There were 77 761 polymorphic SLAF tags accounting for 3.26% of the total SLAF tags. Finally, 129 063 SNPs were found. The population structure and genetic evolution analysis showed that 122 germplasm resources could be divided into three groups, which was not directly related to the geographical source, and the results of clustering could provide scientific basis for parents combination and utilization of heterosis in sweetpotato breeding.

甘薯[Ipomoea batatas (L.) Lam.]既是我国保障粮食安全的底线作物,也是重要的保健食品、工业原料和新型能源作物,其种植面积、单产和总产均居世界前列[1,2]。但由于甘薯遗传背景复杂,基因杂合度高,且具有自交不亲和性等特点,给新品种选育工作带来极大困难。同时,甘薯的品质、产量和抗性等相关性状多表现为数量遗传,采用常规育种手段难以较快实现符合市场需求的育种目标,而分子标记辅助选择技术为解决上述难点提供了新的路径。目前甘薯中RAPD[3]、AFLP[4]、ISSR[5]、SSR[6,7]等标记的开发和利用已取得较大进展,但仍难满足育种需求。特异性位点扩增片段测序技术(specific length amplification fragment sequencing,SLAF-seq)是基于第二代高通量测序技术发展而来,在开发大量特异分子标记方面具有通量高、准确性高、成本低、周期短的优势。该技术已应用于遗传图谱构建、基因组关联分析、QTL定位和种质资源鉴定等研究[8,9]。苏文瑾等[10]首次将SLAF-seq技术应用于甘薯SNP位点的开发,共获得795 794个SNP位点,表明其效率远远高于SSR、AFLP、RAPD等分子标记技术。石璇等[11]以8个甘薯品种或其近缘种为材料开发获得了40 765个SNP位点,利用这些位点构建进化树,发现甘薯栽培种和野生种I. trifida的亲缘关系比较近。为弄清122份甘薯种质资源的亲缘关系,本研究以这些种质为实验材料,通过SLAF-seq技术开发大量特异性SNP位点,并利用开发的位点分析材料之间的遗传关系和群体结构,以期为甘薯新品种选育提供科学依据。

1 材料与方法

1.1 材料

根据前期工作基础[12],从广西农业科学院明阳基地甘薯种质资源圃中选取表型差异较大的122份甘薯种质资源作为研究材料,具体情况见表1
表1 试验中选用的甘薯种质资源

Tab. 1 Sweetpotato germplasms in this study

编号
No.
种质名称
Germplasm name
来源地
Origin
编号
No.
种质名称
Germplasm name
来源地
Origin
编号
No.
种质名称
Germplasm name
来源地
Origin
gs-1 阜菜薯1号 中国安徽 gs-42 桂紫薯3号 中国广西 gs-83 冀薯65 中国河北
gs-2 金薯1号 中国福建 gs-43 桂粉3号 中国广西 gs-84 冀薯99 中国河北
gs-3 金薯2号 中国福建 gs-44 桂薯2号 中国广西 gs-85 冀Y25 中国河北
gs-4 福薯30 中国福建 gs-45 桂薯96-8 中国广西 gs-86 商薯19 中国河南
gs-5 泉薯10号 中国福建 gs-46 桂薯5号 中国广西 gs-87 桂薯11号 中国广西
gs-6 泉薯830 中国福建 gs-47 桂薯6号 中国广西 gs-88 漯17-7 中国河南
gs-7 龙薯21 中国福建 gs-48 桂薯9号 中国广西 gs-89 商79291-1 中国河南
gs-8 龙薯28 中国福建 gs-49 桂薯10号 中国广西 gs-90 安薯07 中国河南
gs-9 宁薯10号 中国福建 gs-50 本地菜 中国广西 gs-91 湘薯541 中国湖南
gs-10 配26 中国福建 gs-51 桂菜薯1号 中国广西 gs-92 苏薯17 中国江苏
gs-11 龙薯9号 中国福建 gs-52 紫叶薯 中国广西 gs-93 苏薯16 中国江苏
gs-12 福薯10号 中国福建 gs-53 东皇薯1号 中国广西 gs-94 宁紫薯1号 中国江苏
gs-13 泉薯9号 中国福建 gs-54 桂经薯4号 中国广西 gs-95 徐薯18 中国江苏
gs-14 榕薯201 中国福建 gs-55 桂经薯8号 中国广西 gs-96 徐薯25 中国江苏
gs-15 金山57 中国福建 gs-56 玉林香薯 中国广西 gs-97 徐薯26 中国江苏
gs-16 新普六 中国福建 gs-57 桂薯90-24-2 中国广西 gs-98 徐薯28 中国江苏
gs-17 福薯89 中国福建 gs-58 来宾丹心薯 中国广西 gs-99 徐紫薯2号 中国江苏
gs-18 金山679 中国福建 gs-59 灌阳红薯 中国广西 gs-100 徐紫薯3号 中国江苏
gs-19 红皮黄心薯 中国福建 gs-60 南洋1号 中国广西 gs-101 徐紫薯4号 中国江苏
gs-20 广紫薯8号 中国广东 gs-61 杨圩薯2 中国广西 gs-102 宁菜04-2 中国江苏
gs-22 广薯214 中国广东 gs-63 桂薯3号 中国广西 gs-104 黄金薯 日本
gs-23 广菜薯3号 中国广东 gs-64 富硒11选 中国广西 gs-105 日本引进 日本
gs-24 湛薯118 中国广东 gs-65 桂薯03-143 中国广西 gs-106 日本紫薯 日本
gs-25 湛薯271 中国广东 gs-66 红爪薯 中国广西 gs-107 九州101 日本
gs-26 普薯24 中国广东 gs-67 马山薯 中国广西 gs-108 烟紫薯1号 中国山东
gs-27 普薯32 中国广东 gs-68 武鸣本地红 中国广西 gs-109 烟薯8 中国山东
编号
No.
种质名称
Germplasm name
来源地
Origin
编号
No.
种质名称
Germplasm name
来源地
Origin
编号
No.
种质名称
Germplasm name
来源地
Origin
gs-28 普薯34 中国广东 gs-69 钦南那丽5 中国广西 gs-110 烟薯11 中国山东
gs-29 广薯87 中国广东 gs-70 黄心白薯 中国广西 gs-111 济薯18 中国山东
gs-30 普薯25号 中国广东 gs-71 钦农1号 中国广西 gs-112 绵薯5号 中国四川
gs-31 广紫薯2号 中国广东 gs-72 良圻薯2 中国广西 gs-113 绵薯6号 中国四川
gs-32 广薯42 中国广东 gs-73 桂薯47 中国广西 gs-114 绵薯7号 中国四川
gs-33 广薯135 中国广东 gs-74 红皮薯 中国广西 gs-115 南薯012 中国四川
gs-34 广薯92-66 中国广东 gs-75 永乐薯2 中国广西 gs-116 南薯007 中国四川
gs-35 黄心薯 中国广东 gs-76 黄心薯 中国广西 gs-117 越南黑薯 越南
gs-36 桔黄心薯 中国广东 gs-77 思布黄皮薯 中国广西 gs-118 心香 中国浙江
gs-37 湛江薯 中国广东 gs-78 桂粉2号 中国广西 gs-119 浙菜薯726 中国浙江
gs-38 广薯95-145 中国广东 gs-79 黔薯5号 中国贵州 gs-120 浙薯255 中国浙江
gs-39 广薯79 中国广东 gs-80 黔薯6号 中国贵州 gs-121 渝紫薯7号 中国重庆
gs-40 桂紫薇薯1号 中国广西 gs-81 黔薯2号 中国贵州 gs-122 渝薯27 中国重庆
gs-41 桂紫薯1号 中国广西 gs-82 冀薯98 中国河北

1.2 方法

1.2.1 DNA提取与质量检测 2017年4月在甘薯种质资源圃中采集每份资源的新展开顶叶,利用康维世纪新型植物基因组DNA提取试剂盒(NuClean Plant Genomic DNA Kit)提取DNA,分别采用琼脂糖凝胶电泳和Nanodrop ND-2000超微量核酸蛋白测定仪检测浓度、纯度和完整性,确保提取的DNA符合建库要求。
1.2.2 酶切建库 本研究开展前,未见甘薯或其近缘物种作为参考基因组的报道,故参照苏文瑾等[10]的方法以马铃薯基因组作为参考基因组进行电子酶切预测,选定RsaⅠ为限制性内切酶,通过酶切片段3′端加A处理,连接双接头后进行PCR扩增,选取目的片段回收建库。
1.2.3 建库及产出数据质量评估 采用Illumina HiSeqTM2500对质检合格的文库进行双端测序,测序结束后通过去接头、去低质量阅读框和去污染处理获得干净序列。同时设置水稻‘日本晴’为对照,进行同流程操作,从比对效率、酶切效率和读长插入片段分布情况来评估实验建库的准确性和有效性。最后通过计算GC含量和Q30评估数据质量。
1.2.4 SLAF标签分析与SNP鉴定 依据Sun等[9]关于SLAF标签和多态性SLAF标签定义,根据序列相似性统计各类标签的数量。同时,选取每个SLAF标签中测序深度最高的序列为参考序列,利用BWA软件[13]将测序读长比对到参考序列上,使用GATK[14]和SAM tools[15] 2种方法开发SNP,选取2种方法共同获得的SNP标记作为最终群体SNP标记。
1.2.5 群体的遗传进化分析 根据获得的最终SNP位点,利用Admixture软件[16]分析种质资源的群体结构,并通过MEGA7软件[17]构建样品的群体进化树。

2 结果与分析

2.1 基因组DNA的提取与检测

DNA琼脂糖凝胶电泳检测结果如图1所示,甘薯种质资源DNA样品主带明显,片段大小一致。
图1 甘薯DNA琼脂糖凝胶电泳检测结果

M:Marker;1~23为样品DNA。

Fig. 1 Agar gel electrophoresis test results of genomic DNA in sweetpotato

M: Marker; 1-23: Different samples genomic DNA.

2.2 建库评估

根据马铃薯参考基因组电子酶切预测结果,选取RsaⅠ为限制性内切酶,长度在264~364 bp的酶切片段为SLAF标签,可预测到225 000个SLAF标签。
实验建库评估主要根据对照水稻‘日本晴’的比对效率、酶切效率和读长插入片段分布情况进行判断。从对照的测序数据分析结果来看,其双端比对效率为95.88%,酶切效率为100%,读长插入片段分布在预期范围之内(图2),表明本实验建库比对效率正常,SLAF建库正常以及测序质量正常。
图2 对照序列插入片段分布

Fig. 2 Distribution of insert fragments of control sequences

2.3 测序数据统计与质量评估

为保证分析质量,本研究中均采用读长100 bp× 2作为后续的数据评估和分析数据,同时用对照水稻‘日本晴’的测序数据评估本实验库的准确性。通过Illumina HiSeqTM2500测序平台测序,从122个甘薯种质资源中共获得563.18 Mb读长数据,不同种质材料的读长个数在1 419 809~8 392 785范围内(图3A)。不同材料的测序质量值Q30在90.61%~96.82%之间,平均Q30为92.78%,对照的Q30为95.92%,表明测序碱基错误率低(图3B)。同时,测序获得的GC含量在36.46%~ 40.33%之间,平均GC含量为38.17%,对照GC含量为40.72%,则表明GC含量均较低符合测序要求(图3B)。
图3 不同甘薯种质资源的测序数据比较

Fig. 3 Comparison of sequencing data in different sweetpotato germplasms

2.4 SLAF标签和SNP标记的开发与鉴定

根据122份甘薯种质资源的测序结果,共开发出SLAF标签238 8759个。不同甘薯种质材料获得的SLAF标签数多少不一致,基本处于151 823~290 909之间(图3C),各个材料的测序深度差异较大,其中不同种质材料的测序总深度在1 071 660~6 833 246范围之间,而每个材料的测序平均深度则在7.06%~24.93%之内,平均测序深度为17.45×。通过对所有SLAF标签进行分型,最终获得多态性的SLAF标签77 761个,占SLAF标签总数的3.26%。
此外,从多态性的SLAF标签中获得SNP标记22 4971个,每个材料获得的SNP标记数目在61 887~163 937范围内,这些SNP的平均完整度在27.50%~72.87%之间,SNP的杂合度在19.49 % ~33.42 %之间(图3D)。依据完整度大于50 %和次要基因型频率(MAF)大于5%的选择标准,对22 4971个群体SNP进行过滤,最终得到129 063个群体SNP位点。

2.5 甘薯种质资源的群体遗传结构与遗传进化分析

利用筛选出的群体SNP位点对122份甘薯种质资源进行群体结构分析,通常根据交叉验证错误率来确定分群数,拥有最低交叉验证错误率的分群数为最优分群数,如图4A图4B所示,当K=3时交叉验证错误率最低,进而说明所选择的122份甘薯种质资源群体可以分为3个类群。
图4 基于SNPs的甘薯种质资源群体遗传结构与遗传进化分析

A:不同K值对应的交叉验证错误率;B:不同分群数(K值)对应的聚类结果;C:122份种质资源的进化树。图中所示为分支长度为8.9216的最优树,聚类树是按比例绘制的,其分支长度与用于推断系统发育树的进化距离的单位相同。

Fig. 4 Analysis of population structure and genetic evolution on sweetpotato germplasms by using SNPs

A: Cross validation error rates corresponding to different K-values; B: Cluster results of different structuring number (K-values);C: Phylogenetic tree of 122 sweetpotato germplasms.The picture shows the optimal tree with a branch length of 8.9216. The cluster tree is drawn in proportion,and its branch length is the same as the unit used to deduce the evolutionary distance provided by the experimental results.

从122份甘薯种质资源的系统进化树(图4C)来看,所有资源根据遗传关系可以分为3大类。其中,第Ⅰ大类包括29份资源,又可分为2个亚类,第Ⅰ-1类包括3份来自福建的资源泉薯830、配26、榕薯201和1份来自广西的资源桂菜薯1号,第Ⅰ-2亚类则包括来自福建的新普六、福薯10号,来自广西的思布黄皮薯、来宾丹心薯、杨圩薯2、钦南那丽5、乌邪爪、本地菜、紫叶薯、东皇薯1号、良圻薯2、红皮薯、黄心白薯、永乐薯2、桂薯11号、桂薯2号、南洋1号、黄心薯,来自四川的绵薯7号,来自广东的广薯135、广薯92-66、湛江薯、广菜薯3号、普薯25、普薯24;第Ⅱ大类包括46份资源,第Ⅲ大类包括47份资源,这2个大类又可分为多个小类,但具体分类结果与甘薯种质资源的地理来源无直接关系。

3 讨论

与传统的分子标记相比,SNP标记是更加有效的遗传标记,因为在大多数基因组中它们是最丰富和稳定的遗传变异形式[18]。对于无参考基因组的物种而言,简化基因组测序技术因为能有效克服基因组复杂、序列与标记信息缺乏等难题,更适宜于SNP大规模的开发和分析。目前,已经报道的简化基因组测序技术主要有,限制性酶切位点相关的DNA(RAD)测序技术[19],基于ⅡB型限制性内切酶位点相关的DNA(Double Digest RAD-seq)测序技术[20],基因分型测序(GBS)[21],特异性长度扩增片段测序(SLAF-seq)技术,这些技术有一个共同点,即通过限制性内切酶来降低基因组DNA的复杂程度。本研究运用SLAF-seq技术,从122份甘薯种质资源中获得563.18 Mb读长数据,共开发出多态性的SLAF标签77 761个,最终得到129 063个SNP位点用于遗传进化分析。其中,SNP的平均完整度在27.50 %~72.87 %之间,SNP的杂合度在19.49 %~33.42 %之间。这与石璇等[11]研究结果有一定的差异,这可能与本研究中平均测序深度较高,获得的SNP位点数量多有一定的关系。
利用SNP技术获得的系统进化图与育成品种的系谱作比较,发现部分育成品种之间有较好的吻合度。如桂粉3号和龙薯21优先聚为一类,再与桂经薯8号、广薯205聚为一类,因为前两者均以金山57为亲本,而桂粉3号与桂经薯8号、广薯205则均以广薯87为母本;绵薯6号、商薯19、徐紫薯3号和徐薯18聚为一类,因为绵薯6号的母本为徐薯18,徐紫薯3号的父本为徐薯18,而商薯19的父本中有一部分徐薯18的血缘;广薯214和广紫薯2号聚为一类,其中广薯214的母本为广紫薯2号;龙薯28与广薯95-145聚为一类,两者均含有同一个父本血缘;广薯79、桂薯9号和普薯32号聚为一类,均含有广薯69的血缘;广薯42、桂薯5号、桂紫薇薯1号和广薯87等聚为一类,这些品种均含有广薯88-70的血缘。因此,利用SNP标记技术分析甘薯种质资源的亲缘关系有一定的可信度,其聚类结果与系谱的亲本来源有较好的一致性。同时也表明SNP标记技术在分析甘薯亲缘关系的可行性。从聚类图中可知,桂薯96-8、九州101、浙薯255、黄金薯优先聚为一类,则可推测这些品种之间有一定的亲缘关系,但与它们的地理来源无直接关系。由此表明,根据地理来源判断甘薯种质的亲缘关系有一定的局限性,研究结果与聂立圆等[7]结论一致。利用SNP标记技术获得甘薯种质资源的亲缘关系,可为甘薯育种在亲本组配、杂种优势利用等方面提供直观的技术支持。
[1]
陆漱韵, 刘庆昌, 李惟基 . 甘薯育种学[M]. 北京: 中国农业出版社, 1998: 74.

[2]
马代夫, 李强, 曹清河 , 等. 中国甘薯产业及产业技术的发展与展望[J]. 江苏农业学报, 2012,28(5):969-973.

[3]
Gichuki S T, Berenyi M, Zhang D P , et al. Genetic diversity in sweetpotato [Ipomoea batatas (L.) Lam.] in relationship to geographic sources as assessed with RAPD markers[J]. Genetic Resources and Crop Evolution, 2003,50(4):429-437.

DOI

[4]
Liu D G, Zhao N, Zhai H , et al. AFLP fingerprinting and genetic diversity of main sweetpotato varieties in China[J]. Journal of Integrative Agriculture, 2012,11(9):1424-1433.

[5]
李强, 刘庆昌, 翟红 , 等. 中国甘薯主要亲本遗传多样性的ISSR分析[J]. 作物学报, 2008,34(6):972-977.

DOI

[6]
Yang X S, Su W J, Wang L J , et al. Molecular diversity and genetic structure of 380 sweetpotato accessions as revealed by SSR markers[J]. Journal of Integrative Agriculture, 2015,14(4):633-641.

DOI

[7]
聂立圆, 李爱贤, 秦桢 , 等. 基于SSR标记的132份甘薯种质指纹图谱的构建及遗传多样性分析[J]. 植物遗传资源学报, 2018,19(5):904-911.

[8]
Zhang Y X, Wang L H, Xin H G , et al. Construction of a high-density genetic map for sesame based on large scale marker development by specific length amplified fragment (SLAF) sequencing[J]. BMC Plant Biology, 2013,13:141.

DOI PMID

[9]
Sun X W, Liu D Y, Zhang X F , et al. SLAF-seq: an efficient method of large-scale de novo snp discovery and genotyping using high-throughput sequencing[J]. PLoS One, 2013,8(3):e58700.

DOI PMID

[10]
苏文瑾, 赵宁, 雷剑 , 等. 基于SLAF-seq技术的甘薯SNP位点开发[J]. 中国农业科学, 2016,49(1):27-47.

DOI

[11]
石璇, 王茹媛, 唐君 , 等. 利用简化基因组技术分析甘薯种间单核苷酸多态性[J]. 作物学报, 2016,42(5):641-647.

DOI

[12]
李慧峰, 陈天渊, 黄咏梅 , 等. 甘薯种质资源形态标记遗传多样性分析[J]. 西南农业学报, 2015,28(6):2401-2407.

[13]
Li H, Durbin R . Fast and accurate short read alignment with Burrows-Wheeler transform[J]. Bioinformatics, 2009,25(14):1754-1760.

DOI PMID

[14]
McKenna A, Hanna M, Banks E , et al. The Genome Analysis Toolkit: a Map Reduce framework for analyzing next-generation DNA sequencing data[J]. Genome Research, 2010,20(9):1297-1303.

DOI PMID

[15]
Li H, Handsaker B, Wysoker A , et al. The sequence alignment/map format and SAMtools[J]. Bioinformatics, 2009,25(16):2078-2079.

DOI PMID

[16]
Alexander D H, Novembre J, Lange K . Fast model-based estimation of ancestry in unrelated individuals[J]. Genome Research, 2009,19(9):1655-1664.

DOI PMID

[17]
Kumar S, Stecher G, Tamura K . MEGA7: Molecular Evolutionary Genetics Analysis Version 7.0 for Bigger Datasets[J]. Molecular Biology and Evolution, 2016,33(7):1870-1874.

DOI PMID

[18]
Liu J, Huang S, Sun M , et al. An improved allele-specific PCR primer design method for SNP marker analysis and its application[J]. Plant Methods, 2012,8(1):34.

DOI PMID

[19]
Miller M R, Dunham J P, Amores A , et al. Rapid and cost effective polymorphism identification and genotyping using restriction site associated DNA (RAD) markers[J]. Genome Research, 2007,17(2):240-248.

DOI PMID

[20]
Peterson B K, Weber J N, Kay E H , et al. Double digest RAD seq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species[J]. PLoS One, 2012,7(5):e37135.

DOI PMID

[21]
Poland J A, Brown P J, Sorrells M E , et al. Development of high-density genetic maps for barley and wheat using a novel two-enzyme genotyping-by-sequencing approach[J]. PLoS One, 2012,7(2):e32253.

DOI PMID

文章导航

/