全 文 ：作物学报 ACTA AGRONOMICA SINICA 2013, 39(6): 1134−1139 http://www.chinacrops.org/zwxb/
ISSN 0496-3490; CODEN TSHPA9 E-mail: xbzw@chinajournal.net.cn

本研究由国家现代农业产业技术体系建设专项(CARS-19)资助。
* 通讯作者(Corresponding author): 关凤芝, E-mail: kj-gfz@163.com
第一作者联系方式: E-mail: wujianzhong176@163.com
Received(收稿日期): 2012-10-19; Accepted(接受日期): 2013-01-15; Published online(网络出版日期): 2013-03-22.
URL: http://www.cnki.net/kcms/detail/11.1809.S.20130322.1739.016.html
DOI: 10.3724/SP.J.1006.2013.01134
亚麻遗传连锁图谱的构建
吴建忠 黄文功 康庆华 赵东升 袁红梅 于 莹 刘 岩 姜卫东
程莉莉 宋喜霞 赵 茜 吴广文 关凤芝*
黑龙江省农业科学院经济作物研究所, 黑龙江哈尔滨 150086
摘 要: 利用DIANE (纤用亚麻栽培种)和宁亚 17 (油用亚麻栽培种)为杂交亲本, 构建 30个 F2单株作为作图群体, 选
用 71对 SRAP和 24对 SSR共显性标记构建了全长为 546.5 cM, 含 12个连锁群(LGs)的亚麻遗传连锁图谱, 标记均
匀分布于 12个连锁群, 每个连锁群有 4~15个标记, 标记间平均距离为 5.75 cM。结果表明, SRAP标记和 SSR标记
中共显性标记适合于亚麻遗传连锁图谱的构建, 但该图谱覆盖的基因组范围较小, 需继续图谱的完整性工作。本研究
为今后亚麻在分子生物学方面的研究提供了基础信息。
关键词: 亚麻; SRAP; SSR; 遗传连锁图谱
Construction of a Genetic Linkage Map in Flax (Linum usitatissimum L.)
WU Jian-Zhong, HUANG Wen-Gong, KANG Qing-Hua, ZHAO Dong-Sheng, YUAN Hong-Mei, YU Ying,
LIU Yan, JIANG Wei-Dong, CHENG Li-Li, SONG Xi-Xia, ZHAO Qian, WU Guang-Wen, and GUAN Feng-Zhi*
Industrial Crops Institute, Heilongjiang Academy of Agricultural Sciences, Harbin 150086, China
Abstract: Using a population with 30 F2 individuals from a cross between DIANE (cultivar for fibre) and NINGYA 17 (cultivar
for oil) was established flax genetic linkage map with 71 SRAP primers and 24 SSR primers. The map contained 12 linkage
groups which spanned a total map length of 546.5 cM, with four to five marker loci per linkage group. The markers on the map
distributed evenly in the 12 linkage groups with an average distance of 5.75 cM between loci. The map merely covered a narrow
range of the flax genome and needed further improvement. The results showed that the codominant markers in SRAP and SSR
were more suitable for the construction of genetic map in flax. Constructing a integrated genetic linkage map of the whole genome
in flax provides important informations for future molecular biology research in flax.
Keywords: Flax; SRAP; SSR; Genetic Linkage Map
亚麻(Linum usitatissimum L.)是一种世界性的纤维作
物和油料作物, 在经济作物生产中具有非常重要的地位,
纤维亚麻在我国的收获面积占世界第一 , 黑龙江省纤维
亚麻的种植面积和总产量占全国 80%以上[1], 同时, 亚麻
又是我国主要的纺织原料 , 其纺织产品是出口创汇的主
要商品之一。油用亚麻主产于印度、加拿大、中国、美国
等, 我国油用亚麻的种植面积仅次于印度和加拿大, 居世
界第 3位, 主要分布于甘肃、内蒙古、新疆等地区。然而,
相比较其他作物 , 亚麻的基础科学研究还有待进一步深
入。传统的育种手段已经无法满足生产中对亚麻品种的需
要, 培育优质、高产、高纤、高抗和具有显著优点的专用
亚麻品种已经是当务之急 , 而分子标记辅助育种将是现
有育种手段的必要补充。迄今, 在水稻、小麦、大豆、油
菜、玉米等大田作物都有以遗传连锁图谱构建为基础的分
子标记辅助品种改良的大量报道 , 这对亚麻的分子研究
起到了很好的借鉴作用。截止目前, 亚麻遗传连锁图谱的
构建在国内还未见有相关的报道, 国外报道过 3 例 [2-4],
图谱情况见表 1, 但均未整合, 无法应用于分子标记辅助
育种的实践。基于亚麻大多数农艺性状都属于数量性状,
其 QTL的研究报道较少, 进行亚麻相关性状 QTL定位和
分子标记辅助育种的前提就是构建一张高密度标记且适
合的亚麻遗传连锁图谱。本研究旨在构建一张适合密度的
亚麻遗传连锁图谱, 为下一步进行亚麻的基因定位、分子
标记辅助育种奠定基础。
第 6期 吴建忠等: 亚麻遗传连锁图谱的构建 1135


表 1 亚麻遗传连锁图谱构建情况
Table 1 Construction of genetic linkage map in flax
参考文献
Reference
年份
Year
作图群体
Mapping population
图谱长度
Map length (cM)
标记类型
Marker type
连锁群
No. of LG
标记数
No. marker
Spielmeyer et al. [2] 1998 DH 1400 AFLP 18 213
Oh et al. [3] 2000 F2 1000 RFLP+RAPD 15 94
Cloutier et al. [4] 2011 DH 833.8 SSR 24 113

1 材料与方法
1.1 试验材料
以纤用亚麻栽培种 DIANE和油用亚麻栽培种宁亚 17
配制杂交组合, 获得 F2代分离群体, 选取其中 30 株作为
作图群体 , 2012 年种植于黑龙江省农业科学院国家高
新技术产业示范园区(民主乡), 大田管理标准, 苗期取单
株叶片 , 参照本实验室方法提取双亲及群体各单株叶片
总 DNA [5], 琼脂糖凝胶电泳检测 DNA质量, 合格样品于
–20℃冰箱保存待用。
1.2 标记分析
100对 SSR引物全部来自加拿大学者Yong-Bi Fu的学术
交流所赠(学者介绍详见 http://www4.agr.gc.ca/AAFC-AAC/
display-afficher.do?id=1181940775526&lang=eng), 编 号 为
Lu1~ Lu100 (表 2), PCR体系含 10 μmol L–1的正、反向 SSR
引物各 1 μL, 10 μmol L–1 dNTPs 0.25 μL, 25 ng μL–1模板
DNA 1 μL, 10×PCR缓冲液 2 μL, 2.5 U μL–1 Taq DNA聚合
酶 0.3 μL, 以 ddH2O 补至 10 μL。PCR 反应在 Biometra
PCR扩增仪上进行。程序为 94℃预变性 4 min; 94℃变性
1 min, 55℃退火 30 s, 72℃延伸 1 min, 35个循环; 72℃延
伸 7 min, 4℃保存。
根据 Li 和 Quiros [6]的原则, 设计 10条锚定外显子
区域的引物和 17 条锚定非编码区的 SRAP 引物, 共得
到 170 对引物组合(表 3)。PCR 体系已经过优化[7], 总体
积为 15 μL, 含 100 ng模板 DNA、引物各 50 ng、1×反应
缓冲液、200 μmol L–1 dNTPs、1.5 μmol L–1 MgCl2和 1 U
Taq DNA聚合酶(MBI公司)。参照 Li和 Quiros [6]的 PCR
程序。
参照厉朝龙等[8]的方法, 2种引物的 PCR扩增后的产
物都经 6%聚丙烯酰胺凝胶电泳分离检测及银染。

表 2 SSR标记编号及序列
Table 2 Identifier and sequence of SSR markers
编号
Identifier
序列
Sequence (5′–3′)
编号
Identifier
序列
Sequence (5′–3′)
F: TGCAATTCGCGATACTGTTC F: AGAATAATGGCCTTGGCTACA Lu10
R: TCCAGATGCAATCTTCTCCA
Lu54
R: TGAACTTGAACCATGGCAGA
F: CATTTCTCCCCTTTTCGTCA F: AATGAGCACAACAACAGCAAG Lu12
R: GTTCGTCCGGAGAAACCTC
Lu55
R: AGCAGCTCTGGACTTGAGGA
F: CCATGAAAAGACGATGACGA F: CAGAAACGAAAGTCTGGAGCA Lu13
R: CGGTTGGTCCACAAAAGAAC
Lu59
R: GCGCTCTCTTCTGCTAATGC
F: GAATTCCCGGGATCTCTCAT F: GGTTTCACTTCATCCGCTTT Lu16
R: AAGTAGACGTACGGGCCTTG
Lu5
R: ACGACCCAGATTTCACTTGG
F: AACCTGAACCAGACGAGCAT F: TATGGGCTCTGCTGCTTTG Lu17
R: CATGGTTTTGAAGGGGTGAT
Lu67
R: GCTAAAGCGCCTTTCTTTCC
F: GCACGTAGGTGTTCTTATGCAA F: CGAGGATCAGAAGTCGAAGG Lu18
R: AGCTTCTTTCAGACGCCATC
Lu68
R: TTCCTTGTCTGAGTGCGAAA
F: AACCTGAACCAGACGAGCAT F: CCTCCTAAAACCGGTCAACA Lu19
AGGTGGATCCAGCAAGCTAA
Lu73
R: AAACCCTTCTGACGCTCATC
F: GCTTTGTCCTCAGGCTTCAC F: CAAGTAATATTCCTGACCACCCTAA Lu21
R: AAGCAAAAGATGCAGAAGCA
Lu79
R: TGGGTAGCATAAGCCTGCTG
F: TCCATCCTCTGCATTTGTGA F: GGGCACGAATTGATTCTCAT Lu22
R: AAGACGAGTGCCCATTCCTA
Lu86
R: AACGACACATGCACAAGCTC
F: TACTGGACGACACCAGCATC F: TTGCAGGGTTCTTCACATCA Lu25
R: GGGAGCAAATCGATCAAGAG
Lu90
R: AAAACCCAGAAAAAGCAGCA
F: CGATGATCACTGGACGGATA F: GGTAAACACAGTGGCTGCTAGA Lu31
R: CATAGCTTCAAAGGCAGCAC
Lu92
R: TTGCTCTCCAATCCAGGTTC
F: GCAGCTTCTGTCATTAGTGCTG F: TTTGCCAACCAAGCTAGACA Lu32
R: GGAGACAAAAGGCCAAGTAAAA
Lu94
R: GATCATCGTCCAGCAGCAG
1136 作 物 学 报 第 39卷

1.3 标记命名及数据记载
SSR标记以所用编号命名, SRAP标记以引物组合名称
命名, 标记名后面加小写字母后缀表示由同一引物对的多
个多态性标记。分子标记的记载方式参照 Mapmaker/Exp
3.0说明书[9], 按显性遗传处理所有扩增出的位点。
1.4 数据分析
统计 F2 群体各位点的基因型, 计算群体各位点的基
因型频率, 用 χ2 测验检验各标记位点基因型的分离情况
(显性 3∶1, 共显性 1∶2∶1)。采用 SPSS 11.5 软件中的
Correlate中的 Bivariate进行相关分析[10]。
1.5 分子标记连锁图构建
采用 Mapmaker/Exp 3.0软件构建亚麻遗传连锁群[11],
以“group”命令对所有的标记分组, LOD 值设定在 3.0, 最
大遗传距离为 50 cM。分组后分别对每组标记用“order”
命令自动排序。对于无法自动排序的 group, 手工排序, 即
先选择 6 个标记用“compare”命令比较排序, 选定最有可
能的一种序列 , 用“try”命令或“build”命令将该组内剩余
标记逐一添加到该序列中, 最后用“ripple”命令确定最佳
排列顺序。采用 Kosambi函数用“map”命令计算连锁距离。
用 MapDraw V2.1绘制亚麻遗传连锁图[12]。
2 结果与分析
2.1 亲本间多态性标记的筛选
分别用 100 对 SSR 引物和 170 对 SRAP 引物扩增亲
本, 表现共显性差异明显且清晰可辨的引物共 56 对, 其
中 SSR引物有 24对, SRAP引物有 32对, 分别占多态性
条带的 42.9%和 57.1%, 共显性多态性标记的总筛出率达
20.7%, 其中 SSR标记共显性筛出率为 24.0%, SRAP标记
为 18.8%。
2.2 群体标记分析
2.2.1 SSR 分析 利用亲本间筛出的 24对 SSR引物分
析亚麻 F2 群体, 群体中能重复检出多态性且带型清晰的
引物有 24对, 共检测出 24个多态性位点, 即每对引物仅
产生 1个多态性位点。图 1所示为引物 Lu10在群体中的
扩增情况, 对图谱上 SSR标记基因型频率进行 χ2检验, 当
P=0.01时, 有 11个 SSR标记表现偏分离, 占 36.7%。

表 3 本研究应用的 SRAP引物
Table 3 SRAP primers used in this study
上游 Upper primer 下游 Lower primer
编号 Code 序列 Sequence (5–3) 编号 Code 序列 Sequence (5–3)
E1 GACTGCGTACGAATT AAT M1 TGAGTCCAAACCGG ATA
E2 GACTGCGTACGAATT TGC M2 TGAGTCCAAACCGG AGC
E3 GACTGCGTACGAATT GAC M3 TGAGTCCAAACCGG AAT
E4 GACTGCGTACGAATT TGA M4 TGAGTCCAAACCGG ACC
E5 GACTGCGTACGAATT AAC M5 TGAGTCCAAACCGG AAG
E6 GACTGCGTACGAATT GCA M6 TGAGTCCAAACCGG TAG
E7 GACTGCGTACGAATT ATG M7 TGAGTCCAAACCGG TTG
E8 GACTGCGTACGAATT AGC M8 TGAGTCCAAACCGG TGT
E9 GACTGCGTACGAATT ACG M9 TGAGTCCAAACCGG TCA
E10 GACTGCGTACGAATT TAG M10 TGAGTCCAAACCGG TAC
E11 GACTGCGTACGAATT TCG
E12 GACTGCGTACGAATT GTC
E13 GACTGCGTACGAATT GGT
E14 GACTGCGTACGAATT CAG
E15 GACTGCGTACGAATT CTG
E16 GACTGCGTACGAATT CGG
E17 GACTGCGTACGAATT CCA


图 1 引物 Lu10在群体中的扩增
Fig. 1 Amplification results of the primer E2M10 in F2 individuals
P1: 母本; P2: 父本; 1~30: F2株系。P1: female parent; P2: male parent; 1–30: F2 plant lines.
第 6期 吴建忠等: 亚麻遗传连锁图谱的构建 1137


2.2.2 SRAP 分析 32 对 SRAP 引物组合共扩增出 99
个带型清晰可辨且可重复检测的共显性位点 , 每个组合
的多态性条带 1~5个不等, 平均每个引物组合产生 3.1个
多态性条带。图 2所示为引物 E4M4在群体中的扩增情况。
在 P=0.01的条件下, 对 99个标记进行卡方检测, 其中 28
个为偏分离标记, 占 28.3%。
2.3 遗传连锁图谱的构建
利用 Mapmaker 软件对 123 个多态性标记(24 个 SSR
标记和 99个 SRAP标记)进行连锁遗传分析。在 LOD≥3.0
的条件下, 共有 95 个标记(包括 24 个 SSR 标记和 71 个
SRAP标记)进入连锁群, 其中 SSR和 SRAP标记进入连锁
群的比例分别是 96%和 73%。最终得到一张包含 12个连
锁群(LGs)的亚麻分子遗传图谱。该图谱总长 546.5 cM
(Kosambi函数), 标记间平均图距为 5.75 cM。不同连锁群
的遗传距离和连锁群上的标记数目变化都较大。连锁群长
度在 20.6~98.7 cM范围内, 标记数介于 5~15间(表 4和图
3)。
2.4 标记在图谱上的分布
经相关分析 , 本研究所得到的标记数目和图谱长度
的相关系数为 0.88, 达到极显著相关水平, 表明这些标记
在亚麻图谱上的分布是均匀的。
3 讨论
3.1 关于用 F2家系构建图谱的可行性
F2 群体是常用的作图群体, 迄今大多数植物的 DNA
标记连锁图谱都是用 F2 群体构建的。对于亚麻这样严格
的自花授粉作物, 建立 F2 群体是很容易的, 这是使用 F2
群体进行亚麻遗传作图的最大优点。但对于 F2 群体标记
分析时表现显性的标记 , 将无法识别显性纯合基因型和
杂合基因型, 由于这种基因型信息简并现象的存在, 会降
低作图的精度。为了提高精度, 减小误差, 一方面可以使
用较大的群体, 另一方面可以选用表现共显性的标记, 从
而会增加DNA标记分析的费用。本研究选用 SRAP和 SSR
标记中表现共显性的部分标记进行亚麻遗传图谱的构建,
这样大大提高了该图谱的精确度 , 但较小的群体数量可
能对本次构建的图谱造成一定偏差。F2群体的一个缺点是
不易长期保存, 有性繁殖一代后, 群体的遗传结构就会发
生变化, 为了延长亚麻 F2 群体的使用时间, 本研究将使
用 F2单株的衍生系, 构建永久 F2群体, 继续进行图谱的
加密及其他图谱完善性工作。


图 2 引物 E4M4在群体中的扩增
Fig. 2 Amplification results of the primer E4M4 in F2 individuals
P1: 母本; P2: 父本; 1~30: F2株系。P1: female parent; P2: male parent; 1–30: F2 plant lines.

表 4 连锁图谱中各个连锁群的标记分布
Table 4 Genetic distances and marker distribution in linkage groups (LGs)
连锁群
Linkage group
标记数目
Number of marker
长度
Distance (cM)
平均图距
Average distance (cM)
LG1 5 21.6 4.32
LG2 12 56.3 4.69
LG3 15 98.7 6.58
LG4 14 76.7 5.48
LG5 7 26.5 3.79
LG6 7 30.6 4.37
LG7 5 31.1 6.22
LG8 8 71.5 8.94
LG9 4 25.6 6.40
LG10 8 63.7 7.96
LG11 5 23.7 4.74
LG12 5 20.6 4.12
1138 作 物 学 报 第 39卷


图 3 亚麻 SSR和 SRAP标记遗传连锁图谱
Fig. 3 A genetic linkage map of SSR and SRAP morphological markers in flax

3.2 关于引物筛选和标记选择
本研究中对引物的筛选全部选择共显性表现的标记,
其中 SSR 标记的共显性表现率为 24%, 且这些标记都为
单位点标记, 最后均进入图谱, 其所揭示的多态性较丰富,
重复率和可信度高, 已成为遗传标记中的热点, 但在亚麻
中开发引物费时且昂贵, 因此开发的 SSR 标记不多; 对
于 SRAP 标记, 17 条上游引物和 10 条下游引物搭配组成
170对 SRAP引物组合, 只有 32对引物组合中的 99个标
记位点表现出共显性, 最终有 71 个 SRAP 标记成功进入
连锁群, 有 28个 SRAP标记未进入连锁群, 占所有标记的
22.76%。Li等[2]认为 SRAP标记具有操作简单、多态性产
率高、广覆盖度和高重复性, 在遗传图谱的构建上有最佳
的标记分布。SRAP 标记是对基因组中 ORFs 进行扩增,
因而对基因相对较少的着丝粒附近以及端粒的扩增会较
少, 因此本研究结合可扩增这些区域的 SSR标记, 将可获
得覆盖整个基因组的连锁图, 利用 SRAP标记和 SSR标记
结合使用可能更有利于构建标记均匀分布的连锁群。
Lorieux 等[13-14]指出选择作用(配子选择和合子选择)对共
第 6期 吴建忠等: 亚麻遗传连锁图谱的构建 1139


显性标记重组值估计的影响小于显性标记 , 标记偏分离
的最简单的克服方法是作图时将偏分离标记剔除掉 , 但
这样会降低对基因组的覆盖率。一个折衷的方法是先利用
正常分离的标记作图, 然后将偏分离标记添加进去, 根据
对原来图谱的影响大小决定是否保留偏分离标记。但是缺
乏一个判断影响大小的标准 , 对于容易出现遗传偏分离
的物种, 遗传作图应尽可能选用共显性标记类型。本研究
用 SRAP 和 SSR 标记进行群体分析中全部选择共显性表
现的标记, 此类标记可鉴别纯合基因型和杂合基因型, 提
供完整的遗传信息 , 因此选择共显性标记在亚麻遗传连
锁图谱的构建中具有重要作用。
3.3 关于图谱质量
本研究构建了一张总长为 546.5 cM的亚麻遗传连锁
图谱, 95 个标记分布于 12 个连锁群上, 标记间平均距离
为 5.75 cM。迄今为止, Spielmeyer等[2]基于 213个 AFLP
标记构建的图谱总长为 1400 cM, 含 18个连锁群, 标记间
平均距离为 10 cM; Oh等[3]用 RFLP和 RAPD共 94个标
记构建了一张含 15个连锁群, 总长达 1000 cM, 标记间距
约 10.6 cM 的遗传连锁图谱; 最新的一张图谱是 Cloutier
等[4]首次使用 SSR标记构建, 含 113个标记, 总长为 833.8
cM, 共 24个连锁群, 标记间平均距离为 7.3 cM。与已发
布图谱比较 , 本研究构建的图谱在标记密度上有了很大
的提高, 标记分布比较均匀, 但图谱总长度还有待增加标
记数量来实现。
本研究所构建的图谱覆盖基因组范围较小。一个重要
原因可能是图谱构建中大多数(74.74%)用 SRAP 标记, 其
扩增区域主要是基因组的编码区 , 对基因分布相对较少
的非编码区检测机会降低, 使得所构建的连锁图谱缩短,
虽然在图谱构建过程中结合了 SSR 标记以扩增着丝粒及
端粒区域, 但从最后结果来看, 标记分布比较集中, 未能
完全覆盖亚麻全基因组(2n=30), 仅集中于 12 个连锁群,
还是远远不能满足饱和图谱构建的需要。造成基因组覆盖
度偏小的另一个原因可能是本研究所用作图群体偏小 ,
作图群体的大小很大程度上决定了遗传图谱的分辩率和
精确度, 群体越大, 作图精度越高, 但费时费力, 构建分
子标记骨架连锁图 , 则可基于大群体中的一个随机小群
体, 本研究将继续在此基础上扩大群体规模, 构建合适的
群体 , 增加多种类型的分子标记进行遗传图谱的加密工
作。图谱长度较小的另一个可能原因是双亲间的差异位点
较少 , 构建的图谱无法覆盖到整个亚麻基因组的实际长
度。本研究下一步拟通过增加现有标记类型及数量和加大
作图群体的数量来弥补遗传连锁图谱对整个基因组的实
际覆盖长度不足这一弊端。

致谢: 感谢华中农业大学周永明教授及其研究课题组对
本研究的指导和帮助。
References
[1] FAO Databank. http://faostat.fao.org/site/567/default.aspx [2008-
12-01]
[2] Spielmeyer W, Green A G, Bittisnich D, Mendham N, Lagudah E
S. Identification of quantitative trait loci contributing to Fusarium
wilt resistance on an AFLP linkage map of flax (Linum usitatis-
simum L.). Theor Appl Genet, 1998, 97: 633–641
[3] Oh T J, Gorman M, Cullis C A. RFLP and RAPD map-ping in
flax (Linum usitatissimum). Theor Appl Genet, 2000, 101:
590–593
[4] Cloutier S, Ragupathy R, Niu Z, Duguid S. SSR-based linkage
map of flax (Linum usitatissimum L.) and mapping of QTLs un-
derlying fatty acid composition traits. Mol Breed, 2011 28:
437–451
[5] Wu J-Z(吴建忠). Analysis of DNA extraction method in flax ge-
nomic. Heilongjiang Agric Sci (黑龙江农业科学), 2011, (7):
18–19 (in Chinese with English abstract)
[6] Li G, Quiros C F. Sequence-related amplified polymorphism
(SRAP), a new marker system based on a simple PCR reaction:
its application to mapping and gene tagging in Brassica. Theor
Appl Genet, 2001, 103: 455–461
[7] Wu J-Z(吴建忠), Huang W-G(黄文功), Zhao D-S(赵东升), Kang
Q-H(康庆华), Liu Y(刘岩), Zhao Q(赵茜). Optimization of
SRAP reaction system and screening of polymorphic SRAP
markers in flax. Plant Fiber Sci China (中国麻业科学), 2011, (6):
281–284 (in Chinese with English abstract)
[8] Li C-L(历朝龙), Chen S-Q(陈枢青), Liu Z-Y(刘子骀), Shen
Q-G(沈奇桂), Zhao L-H(赵鲁杭). Biochemistry and Molecular
Biology Experiment Technology (生物化学和分子生物学实验
技术). Zhejiang: Zhejiang University Press, 1999 (in Chinese)
[9] Lander E S, Green P, Abrahamson J, Barlow A, Daly M J, Lin-
coln S E, Newburg L. MAPMAKER: An interactive computer
package for constructing primary genetic linkage maps of ex-
perimental and natural populations. Genomics, 1987, 1: 174–181
[10] Zhang W-T(张文彤), Yan J(闫洁). SPSS Statistical Analysis
Foundation Course (SPSS统计分析基础教程). Beijing: Higher
Education Press, 2004. pp 257–278 (in Chinese)
[11] Lincoln S E, Daly M E, Lander E S. Constructing Genetic Maps
with Mapmaker/Exp 3.0. Cambridge, MA, USA: Whitehead In-
stitute for Biomedical Research, 1992
[12] Liu R-H(刘仁虎), Meng J-L(孟金陵). MapDraw: a Microsoft
Excel macro for drawing genetic linkage maps based on given
genetic linkage data. Hereditas (遗传), 2003, 25(3): 317–321 (in
Chinese with English abstract)
[13] Lorieux M, Goffinet B, Perrier X, Gonzalez de Leon, Lanaud C.
Maximum-likelihood models for mapping genetic markers
showing segregation distortion: 1. Backcross populations. Theor
Appl Genet, 1995, 90: 73–80
[14] Lorieux M, Perrier X, Goffinet B, Lanaud C, Gonzalez de Leon.
Maximum-likelihood modles for mapping genetic markers
showing segregation distortion: 2. F2 populations. Theor Appl
Genet, 1995, 90: 81–89