与小麦抗白粉病基因 Pm48紧密连锁分子标记的开发
付必胜1,**, 刘颖1,2,**, 张巧凤1, 吴小有1, 高海东3, 蔡士宾1, 戴廷波2,*, 吴纪中1,*
1江苏省农业科学院粮食作物研究所 / 江苏省农业种质资源保护与利用平台, 江苏南京 210014
2南京农业大学农学院, 江苏南京 210095
3南京集思慧远生物科技有限公司, 江苏南京 210014
*通讯作者(Corresponding authors): 吴纪中, E-mail: wujz@jaas.ac.cn, Tel: 025-84391667; 戴廷波, E-mail: tingbod@njau.edu.cn, Tel: 025-84395033

第一作者联系方式: E-mail: fbs1006@126.com **同等贡献(Contributed equality to this work)

摘要

Pm48为本实验室鉴定的一个抗白粉病新基因。为精细定位该基因, 利用混池ddRAD测序鉴定了81个与该基因关联的序列, 开发了STS标记 Xmp931, 转化了CAPS标记 Xmp928 Xmp930 Xmp936; 同时, 利用粗山羊草基因组序列开发了71个基因组SSR标记, 定位了其中的 Xmp1089 Xmp1112。在115个宁糯麦1号×Tabasco衍生的 F2:3家系中, Xmp928与目的基因共分离, Xmp1112位于近着丝粒方向处距抗病基因3.1 cM。在671个纯合感病家系中, 标记 Xmp928仍与目的基因共分离。利用3个中国春5DS缺失系, 最终将 Pm48定位在小麦5DS上0.63-0.67的臂区段中。

关键词: 小麦; 抗白粉病基因; 分子标记; 混池ddRAD测序
Development of Markers Closely Linked with Wheat Powdery Mildew Resistance Gene Pm48
FU Bi-Sheng1,**, LIU Ying1,2,**, ZHANG Qiao-Feng1, WU Xiao-You1, GAO Hai-Dong3, CAI Shi-Bin1, DAI Ting-Bo2,*, WU Ji-Zhong1,*
1 Institute of Food Crops, Jiangsu Academy of Agricultural Sciences / Jiangsu Provincial Platform for Conservation and Utilization of Agricultural Germplasm, Nanjing 210014, China
2 College of Agriculture, Nanjing Agricultural University, Nanjing 210095, China
3 Genepioneer Biotechnologies Co. Ltd., Nanjing 210014, China
Abstract

Pm48 is a novel powdery mildew resistance gene identified previously in our laboratory. This study aimed at developing close molecular markers for fine mapping of the gene. The ddRAD-sequencing assay revealed 81 SNPs associated with the target gene, in which one converted into the STS marker Xmp931 and three converted into the CAPS markers Xmp928, Xmp930, and Xmp936. We also developed 71 genomic SSR markers according to the genome sequence of Aegilops tauschii. And mapped two of them, Xmp1089 and Xmp1112. Using the 115 F2:3 families derived from the cross of Ningnuomai 1 × Tabasco, the target gene was found to be co-segregated with Xmp928 and distal to Xmp1112 with the genetic distance of 3.1 cM towards centromere. In the 671 homozygous susceptible families, Xmp928 also showed co-segregated with the target gene. We also physically mapped Pm48 to the bin of 5DS 0.63-0.67 by using three Chinese Spring 5DS deletion lines.

Keyword: Wheat; Powdery mildew resistance gene; Molecular markers; Bulked ddRAD-seq

由禾布氏白粉病菌(Blumeria graminis f. sp. tritici)引起的白粉病是小麦的重要病害, 在严重流行年份可导致小麦减产34%以上[1]。抗病品种在小麦白粉病防治中发挥了重要作用, 但白粉病菌的变异会导致抗病品种丧失抗性, 因此需要不断挖掘新的抗病基因来培育小麦抗白粉病新品种。

随着分子标记技术的发展, 分子标记的类型逐渐丰富, 如初期的RFLP (restriction fragment length polymorphism)、AFLP (amplified fragment length polymorphism), 以及根据序列开发的SSR (simple sequence repeats)、STS (sequence tagged sites)、CAPS (cleaved amplified polymorphic sequence)、dCAPS (derived cleaved amplified polymorphic sequence)、SNP (single nucleotide-acid polymorphism)等。分子标记类型的丰富以及越来越成熟的分子标记技术极大地推动了小麦抗白粉病新基因的发掘, 越来越多的小麦抗白粉病基因通过分子标记技术得到鉴定, 如Pm47[2]Pm50[3]Pm51[4]等。迄今为止, 已经在50个位点报道了80多个抗白粉病基因[5], 分布在小麦除4D和5A以外的所有染色体上。

由于小麦还没有完整的基因组序列, 在特定染色体区域开发标记一般需要借助物理定位的EST、共线性关系、基因组序列、SNP芯片序列信息开发标记。在小麦中, 不同倍性小麦的遗传和物理图谱信息、网上公布的水稻、短柄草、大量小麦EST序列信息和小麦基因组序列信息等是用于饱和目标基因区域遗传图谱的重要标记来源。这些信息为开发新的标记奠定了良好的基础。

近年蓬勃发展的高通量测序也是基因型分析的一种新策略。为了降低复杂基因组的测序费用, 研究者已经提出了基于第二代测序技术的简化代表文库测序策略[6], 其中最简单的方法是分离纯化特定片段大小范围的限制性酶切片段。如, 利用RAD (restriction site-associated DNA, 限制性酶切位点关联DNA)测序方法已成功构建了若干高密度遗传图谱[7, 8]。而ddRAD (double-digest RAD, 双酶切RAD)测序是RAD衍生的一种测序策略[9], 与RAD的单末端测序相比, ddRAD的双末端测序能够增加标记的特异性和准确性。近年, 在花生[10]、草莓[11]、油菜[12]等物种中已经利用该测序方法构建了若干高密度遗传图谱。

目前, 已经克隆的小麦抗病基因有Lr10[13]Vrn1[14]Yr36[15]Lr34/Yr18/Pm38[16]等, 而抗白粉病基因只克隆了Pm3[17]Pm21[18]。本实验室在德国小麦品种Tabasco中鉴定到抗白粉病基因Pm46, 位于5DS染色体上[19], 后由于基因重名被重新命名为Pm48[20]。本实验室鉴定了与Pm48紧密连锁的分子标记, 但只包括1个显性STS标记和5个SSR标记[19]。显性标记在分离群体中的带型信息相对较少, 不利于该基因的有效利用; 另外, SSR标记由于是基于重复序列的标记类型, 没有足够的序列信息, 不利于抗病基因的比较基因组学分析及进一步精细定位。因此, 本实验室深入研究并开发了与Pm48紧密连锁的新类型分子标记, 为该基因的分离和有效利用奠定基础。

1 材料与方法
1.1 试验材料

利用抗小麦白粉病的德国品种Tabasco与本实验室培育的糯性感病品种糯麦1号, 经杂交和自交获得4129个宁糯麦1号× Tabasco F2单株, 用于抗性遗传分析; 从436个F2:3家系[19]中随机选择115个家系组成分离群体用于新开发分子标记的多态性分析和定位。

利用“ 中国春” 及其缺失系del5DS-1、del5DS-2和del5DS-5进行目标分子标记的染色体臂定位, 其中del5DS-1保留有5DS末端63%长度, 缺失了37%的染色体长度, del5DS-5保留有染色体5DS末端67%长度, 缺失了33%的染色体长度, del5DS-2保留有5DS末端78%长度, 缺失了22%的染色体长度[21]

1.2 抗病性评价方法

参考高海东等[19]报道的方法评价白粉病抗性。将小麦种子播在72孔的穴盘, 亲本品种及家系各播种15~20粒, 每穴盘随机种植9粒感病对照苏麦3号。于幼苗一叶期人工接种, 将充分发病的足量苏麦3号叶片悬于待接叶片上方, 轻轻抖落病菌孢子, 菌种为南京地区流行的Bgt18。在人工智能气候室(14 h光照/10 h黑暗, 18~22℃, 相对湿度80%~90%)中培养接种后的麦苗。接种后7 d当感病对照表现出明显病症时调查病情, 按盛宝钦[22]的0~4级分类方法记录, 即0级为免疫, 植株无病斑, 0; 级为坏死反应, 叶片有枯死斑, 1级为高抗, 病斑直径小于l mm, 菌丝层稀薄可见绿色叶面, 偶见较大病斑, 但仍透绿, 产孢量极少, 2级为中抗, 叶片病斑直径小于l mm, 但菌丝层较厚, 不透绿, 能产生一定量孢子, 3级为中感, 叶片病斑多, 且直径大于l mm, 菌丝层厚, 产孢量大, 但病斑不连片, 4级为高感, 叶片病斑直径大于l mm, 菌丝层厚, 产孢量多, 病斑连片。

1.3 分子标记开发

1.3.1 基于ddRAD测序获得的序列 根据前期F2:3家系的分子鉴定结果[19], 从中随机选取50个纯合抗病家系和50个纯合感病家系, 每家系选取8株, 将其DNA等量混合构建抗池(BR)和感池(BS), 由南京集思慧远生物科技有限公司进行ddRAD测序。根据抗感池序列SNP是否产生酶切位点差异, 决定转化成CAPS或dCAPS标记。如果抗、感池间的SNP存在酶切位点的差异, 则直接将该SNP标记转化成CAPS标记, 即在SNP位点两侧设计引物, 并利用相应的限制性内切酶进行酶切; 如果抗、感池间的SNP不存在酶切位点的差异, 则利用软件dCAPS Finder 2.0 (http://helix.wustl.edu/dcaps/dcaps.html)和Primer 5.0将SNP转化成dCAPS标记, 即在引物中引入错配与SNP位点形成新的酶切位点。新开发的多态标记信息见表1

表1 基于ddRAD测序和粗山羊草基因组序列新开发的标记信息 Table 1 Information of primers newly developed based on ddRAD-seq and Aegilops tauschii genome sequence

1.3.2 基于粗山羊草基因组序列 高海东等[19]鉴定出Pm48的紧密连锁标记Xmp510。本研究利用Xmp510序列及ddRAD测序获得的多态标记序列, 在线搜索比对粗山羊草标记序列(http://probes.pw.usda.gov/WheatDMarker/ phpblast/blast.php)。先对相应的粗山羊草基因组序列进行SSR筛查, 然后用MACVECTOR V8.0 (Accelrys, UK)设计引物, 并进行亲本及抗感池间多态性筛选。新开发的多态标记信息见表1

1.4 分子标记分析

参照Ma等[23]报道的SDS法, 从小麦嫩叶中提取总DNA, -4℃保存。在SENSOQUEST LABCYCLER PCR仪(德国)上扩增, 反应体系为10 µ L, 包括1× buffer 10 µ L (含MgCl2 15 nmol), DNA模板1 ng, 上、下游引物各2 pmol, 2 nmol的dNTPs, 1 U Taq酶。PCR程序为94℃预变性5 min, 然后按94℃变性30 s、50~60℃退火45 s、72℃延伸50 s进行36个循环扩增, 最后72℃延伸7 min。扩增产物经8%非变性聚丙烯酰胺凝胶电泳, 银染法显色。

CAPS标记选择特定的限制性内切酶(TaKaRa), 如Alu I、Acc I、Mae II、Hpy 188I和Taq I, 按生产商说明书进行酶切反应, 反应体系8 µ L, 含PCR产物约50 ng, 于37° C或65° C (依酶特性而定)过夜, 用8%非变性聚丙烯酰胺凝胶电泳分离检测酶切产物。

1.5 遗传图谱构建

采用MapMaker 3.0b/EXP[24]和MapDraw V2.1[25]构建连锁图谱, 选用Kosambi函数[26]计算遗传距离, 连锁判断LOD阈值设为3.0。

2 结果与分析
2.1 与Pm48连锁分子标记的开发及定位

对抗池、感池进行ddRAD测序, 获得81个可能与Pm48关联的SNP序列, 将其中14个关联性强的SNP标记转化为STS或者CAPS标记。其中, Xmp931在抗、感亲本间有多态性, 且与抗、感池表现一致; Xmp928(图1-A)、Xmp930(图1-B)和Xmp936经相应的限制性内切酶酶切后, 在抗、感亲本以及池间表现一致多态性。利用这4个新开发的标记检测115个F2:3家系, 结果4个标记均被定位于Pm48的遗传图谱上, 其中Xmp928与目标基因共分离, Xmp931Xmp936Xmp930位于远着丝粒的一侧, 与目标基因分别相距7.8、9.7和15.7 cM (图2-A)。

图1 Xmp928(A)、Xmp930(B)和Xmp1089(C)在两亲本及抗、感池的扩增带型
M: 50 bp DNA ladder; 1: Tabasco; 2: 宁糯麦1号; 3: 抗病池; 4: 感病池。
Fig. 1 Polymorphism patterns between the parents and the bulks detected with markers Xmp928(A)、Xmp930(B), and Xmp1089(C)
M: 50 bp DNA ladder; 1: Tabasco; 2: Ningnuomai 1; 3: resistant pool; 4: susceptible pool.

图2 抗白粉病基因Pm48的遗传图谱(A)和物理图谱(B)
A图中左侧数字为遗传距离(cM); B图中右侧括号内数字是着丝粒到染色体断点处长度与染色体短臂总长度的比值, 下方椭圆代表着丝粒。
Fig. 2 Genetic (A) and physical (B) maps of the powdery mildew resistance gene Pm48
In linkage map A, numbers on the left indicate genetic distances (cM). In physical map B, numbers in the parentheses indicate the ratio of length from centromere to chromosome breakpoint/total length of the short arm of chromosome and the ellipse at the bottom shows the position of centromere.

利用已定位Pm48连锁标记序列信息, 比对粗山羊草基因组序列, 开发了71个gSSR标记, 其中Xmp1089(图1-C)和Xmp1112 (图3)在抗、感基因型间存在多态性, 并且被整合到Pm48连锁图谱中, 距Pm48分别为3.1 cM和10.8 cM (图2-A)。

图3 Xmp1112在杂交亲本和中国春5DS缺体中的扩增带型
M: 50 bp DNA ladder; 1: Tabasco; 2: 宁糯麦1号; 3~5: 中国春缺体系del5DS-1、del5DS-2和del5DS-5; 6: 中国春。箭头指示目标带缺失。
Fig. 3 Profiles of Xmp1112 amplified in cross parents and Chinese Spring 5DS deletion lines
M: 50 bp DNA ladder; 1: Tabasco; 2: Ningnuomai 1; 3-5: Chinese Spring deletion lines del5DS-1, del5DS-2 and del5DS-5; 6: Chinese Spring. The arrow shows vacancy of target band.

2.2 Pm48的染色体臂定位

利用3个中国春5DS缺失系对Pm48进行染色体臂定位, 结果显示Xmp1112位于小麦5DS 0.63~0.67的染色体区段(图3), 在该染色体区段中已定位标记Xcfd81[19], 同时连锁图谱显示Pm48位于这2个标记之间(图2-A), 因此认为Pm48位于小麦5DS 0.63~0.67的染色体区段(图2-B)。

2.3 重组体的筛选

利用4129个F2单株, 通过接种Bgt18鉴定抗病性, 结果有3103个抗病单株和1026个感病单株, 符合单个显性基因的3︰1分离(χ 2=0.04; χ 20.05, 1=3.84)。这些感病单株移栽后完成整个生育期并收获种子的有671个家系, 对每个家系衍生的15~20个植株进行表型鉴定后, 推测均为纯合感病家系。利用位于Pm48两侧的分子标记对这671个纯合感病家系进行基因型检测, 其中共显性标记Xmp1112Xcfd81分别筛选到12个和8个重组家系; 标记Xmp928Xmp510对这19个重组体进行基因型的分析结果表明, Xmp510检测到5个重组体, Xmp928没有检测到重组体, 与目标基因共分离(表2)。

表2 19个纯合感病重组体基因型分析和分类 Table 2 Genotype classification of the 19 homozygous susceptible recombinants
3 讨论

为进一步精细定位小麦抗白粉病基因Pm48, 本研究利用简化测序和粗山羊草基因组序列信息开发了6个新标记。与实验室之前发表的遗传图谱[19]相比, 本研究在Xgwm205与基因之间新增了2个gSSR标记(Xmp1089Xmp1112), 将Pm48界定在Xmp510Xmp1112之间的4.6 cM之间, 为基因的进一步分离奠定了基础。此外, 还开发了一个与基因共分离的标记Xmp928, 该标记无论在小群体中还是在扩大的纯合感病家系中均与抗病基因共分离, 将作为抗病基因Pm48的重要选择标记。Xmp928是重复序列连接(repeat junction)类型的标记, 这类标记在重复序列超过80%的小麦基因组中非常丰富[27], 将是小麦基因定位的重要标记类型, 然而该类型标记无法与基因序列关联, 很难为基因分离提供帮助。

简化测序可以降低测序成本, 根据性状混池后测序更进一步降低成本。本研究利用ddRAD测序成功开发了3个CAPS标记和1个STS标记, 转化了与基因关联性强的14个SNP标记, 其中4个(28.6%)为多态性标记, 而其余标记没有多态性。对于没有转化成功的10个SNP标记, 推测是由于ddRAD测序后获得的重复序列导致非特异性以及扩增引入的错配。本研究还利用粗山羊草基因组序列开发了71个gSSR标记, 但只有2个被定位, 多态率为2.8%, 远低于其他研究报道的gSSR的多态率[28, 29, 30]。推测可能的原因, 一是D基因组的多态性偏低[28, 29, 30], 二是Tabasco、宁糯麦1号两亲本的分子标记多态率低。对比上述两种方法, 利用混池的ddRAD测序效率较高。

本研究利用5DS缺失系成功将抗病基因Pm48定位在小麦5DS 0.63-0.67的染色体区段中, 这一区段的物理距离大约为55 Mb [31]。标记Xcfd81Xmp1112相距8.0 cM, 假定两标记位于边界, 该区域的重组率估算为0.145 cM/Mb, 低于Erayman等[32]在该区域估算的重组率(0.341 cM/Mb), 并非重组热点区域。尽管如此, 该区段较小的物理距离仍然将有助于抗病基因的分离。本实验室将在以下两方面进一步研究, 一是解决标记多态性低的问题, 准备利用不同的感病亲本与Tabasco构建多个分离群体, 构建Pm48的饱和遗传图谱; 二是利用混池进行转录组测序, 发掘Pm48候选基因。

The authors have declared that no competing interests exist.

作者已声明无竞争性利益关系。The authors have declared that no competing interests exist.

参考文献
[1] Johnson J, Baenziger P, Yamazaki W, Smith R. Effects of powdery mildew on yield and quality of isogenic lines of ‘Chancellor’ wheat. Crop Sci, 1979, 19: 349-352 [本文引用:1]
[2] Xiao M, Song F, Jiao J, Wang X, Xu H, Li H. Identification of the gene Pm47 on chromosome 7BS conferring resistance to powdery mildew in the Chinese wheat land race Hongyanglazi. Theor Appl Genet, 2013, 126: 1397-1403 [本文引用:1]
[3] Mohler V, Bauer C, Schweizer G, Kempf H, Hartl L. Pm50: a new powdery mildew resistance gene in common wheat derived from cultivated emmer. J Appl Genet, 2013, 54: 259-263 [本文引用:1]
[4] Zhan H, Li G, Zhang X, Li X, Guo H, Gong W, Jia J, Qiao L, Ren Y, Yang Z, Chang Z. Chromosomal location and comparative genomics analysis of powdery mildew resistance gene Pm51 in a putative wheat-Thinopyrum ponticum introgression line. PLoS One, 2014, 9: e113455 [本文引用:1]
[5] Xu H, Yi Y, Ma P, Qie Y, Fu X, Xu Y, Zhang X, An D. Molecular tagging of a new broad-spectrum powdery mildew resistance allele Pm2c in Chinese wheat land race Niaomai. Theor Appl Genet, 2015, 128: 2077-2084 [本文引用:1]
[6] Hyten D, Cannon S, Song Q, Weeks N, Fickus E, Shoemaker R, Specht J, Farmer A, May G, Cregan P. Highthroughput SNP discovery through deep resequencing of a reduced representation library to anchor and orient scaffolds in the soybean whole genome sequence. BMC Genomics, 2010, 11: 38 [本文引用:1]
[7] Pfender W, Saha M, Johnson E, Slabaugh M. Mapping with RAD (restriction-site associated DNA) markers to rapidly identify QTL for stem rust resistance in Lolium perenne. Theor Appl Genet, 2011, 122: 1467-1480 [本文引用:1]
[8] Wang N, Fang L, Xin H, Wang L, Li S. Construction of a high-density genetic map for grape using next generation restriction-site associated DNA sequencing. BMC Plant Biol, 2012, 12: 148 [本文引用:1]
[9] Peterson B K, Weber J N, Kay E H, Fisher H S, Hoekstra H E. Double digest RADseq: an inexpensive method for de novo SNP discovery and genotyping in model and non-model species. PLoS One, 2012, 7: e37135 [本文引用:1]
[10] Davik J, Sargent D J, Brurberg M B, Lien S, Kent M, Alsheikh M. A ddRAD based linkage map of the cultivated strawberry, Fragaria xananassa. PLoS One, 2015, 10: e0137746 [本文引用:1]
[11] Zhou X, Xia Y, Ren X, Chen Y, Huang L, Huang S, Liao B, Lei Y, Yan L, Jiang H. Construction of a SNP-based genetic linkage map in cultivated peanut based on large scale marker development using next-generation double-digest restriction-site-associated DNA sequencing (ddRADseq). BMC Genomics, 2014, 15: 351 [本文引用:1]
[12] Wu Z, Wang B, Chen X, Wu J, King G, Xiao Y, Liu K. Evaluation of linkage disequilibrium pattern and association study on seed oil content in Brassica napus using ddRAD sequencing. PLoS One, 2016, 11: e0146383 [本文引用:1]
[13] Feuillet C, Travella S, Stein N, Albar L, Nublat A, Keller B. Map-based isolation of the leaf rust disease resistance gene Lr10 from the hexaploid wheat (Triticum aestivum L. ) genome. Proc Natl Acad Sci USA, 2003, 100: 15253-15258 [本文引用:1]
[14] Yan L, Loukoianov A, Tranquilli G, Helguera M, Fahima T, Dubcovsky J. Positional cloning of the wheat vernalization gene VRN1. Proc Natl Acad Sci USA, 2003, 100: 6263-6268 [本文引用:1]
[15] Fu D, Uauy C, Distelfeld A, Blechl A, Epstein L, Chen X, Sela H, Fahima T, Dubcovsky J. A kinase-START gene confers temperature-dependent resistance to wheat stripe rust. Science, 2009, 323: 1357-1360 [本文引用:1]
[16] Krattinger S, Lagudah E, Spielmeyer W, Singh R, Huerta-Espino J, McFadden H, Bossolini E, Selter L, Keller B. A putative ABC transporter confers durable resistance to multiple fungal pathogens in wheat. Science, 2009, 323: 1360-1363 [本文引用:1]
[17] Yahiaoui N, Srichumpa P, Dudler R, Keller B. Genome analysis at different ploidy levels allows cloning of the powdery mildew resistance gene Pm3b from hexaploid wheat. Plant J, 2004, 37: 528-538 [本文引用:1]
[18] Cao A, Xing L, Wang X, Yang X, Wang W, Sun Y, Qian C, Ni J, Chen Y, Liu D, Wang X, Chen P. Serine/threonine kinase gene Stpk-V, a key member of powdery mildew resistance gene Pm21, confers powdery mildew resistance in wheat. Proc Natl Acad Sci USA, 2011, 108: 7727-7732 [本文引用:1]
[19] Gao H, Zhu F, Jiang Y, Wu J, Yan W, Zhang Q, Jacobi A, Cai S. Genetic analysis and molecular mapping of a new powdery mildew resistant gene Pm46 in common wheat. Theor Appl Genet, 2012, 125: 967-973 [本文引用:8]
[20] McIntosh R, Dubcovsky J, Rogers W, Morris W, Appels R, Xia X. Catalogue of gene symbols for wheat: 2013-2014 supplement, http://www.wheat.pw.usda.gov/GG2/pubs.shtml [本文引用:1]
[21] Endo T R, Gill B S. The deletion stocks of common wheat. J Hered, 1996, 87: 295-307 [本文引用:1]
[22] 盛宝钦. 用反应型记载小麦苗期白粉病. 植物保护, 1988, (1): 14
Sheng B Q. Infection reaction types against wheat powdery mildew at seedling stage. Plant Prot, 1988, (1): 14 (in Chinese) [本文引用:1]
[23] Ma Z Q, Sorrells M E, Tanksley S D. RFLP markers linked to powdery mildew resistance genes Pm1, Pm2, Pm3, and Pm4 in wheat. Genome, 1994, 37: 871-875 [本文引用:1]
[24] Land er E, Green P, Abrahamson J, Barlow A, Daley M, Lincoln S, Newburg L. Mapmaker: an interactive computer package for constructing primary genetic linkage maps of experimental and natural populations. Genomics, 1987, 1: 174-181 [本文引用:1]
[25] 刘仁虎, 孟金陵. MapDraw, 在Excel中绘制遗传连锁图的宏. 遗传, 2003, 25: 317-321
Liu R H, Meng J L. MapDraw: a Microsoft Excel macro for drawing genetic linkage maps based on given genetic linkage data. Hereditas (Beijing), 2003, 25: 317-321 (in Chinese with English abstract) [本文引用:1]
[26] Kosambi D. The estimation of map distances from recombination values. Ann Eugenics, 1943, 12: 172-175 [本文引用:1]
[27] You F M, Wanjugi H, Huo N, Lazo G R, Luo M C, Anderson O D, Dvorak J, Gu Y Q. RJPrimers: unique transposable element insertion junction identification and primer design for marker development. Nucl Acid Res, 2010, 38: 313-320 [本文引用:1]
[28] Xue S, Zhang Z, Lin F, Kong Z, Cao Y, Li C, Yi H, Mei M, Zhu H, Wu J, Xu H, Zhao D, Tian D, Zhang C, Ma Z. A high-density intervarietal map of the wheat genome enriched with markers derived from expressed sequence tags. Theor Appl Genet, 2008, 117: 181-189 [本文引用:2]
[29] Suenaga K, Khairallah M, William H, Hoisington D. A new intervarietal linkage map and its application for quantitative trait locus analysis of “gigas” features in bread wheat. Genome, 2005, 48: 65-75 [本文引用:2]
[30] Paillard S, Schnurbusch T, Winzeler M, Messmer M, Sourdille P, Abderhalden O, Keller B, Schachermayr G. An integrative genetic linkage map of winter wheat (Triticum aestivum L. ). Theor Appl Genet, 2003, 107: 1235-1242 [本文引用:2]
[31] Akpinar B, Magni F, Yuce M, Lucas S, Šimková H, Šafář J, Vautrin S, Bergès H, Cattonaro F, Doležel J, Budak H. The physical map of wheat chromosome 5DS revealed gene duplications and small rearrangements. BMC Genomics, 2015, 16: 453 [本文引用:1]
[32] Erayman M, Sand hu D, Sidhu D, Dilbirligi M, Baenziger P S, Gill K S. Demarcating the gene-rich regions of the wheat genome. Nucl Acids Res, 2004, 32: 3546-3565 [本文引用:1]