全 文 :植物病理学报
ACTA PHYTOPATHOLOGICA SINICA 46(1): 91 ̄102(2016)
收稿日期: 2015 ̄01 ̄28ꎻ 修回日期: 2015 ̄06 ̄26
基金项目: 国家高技术研究发展计划(863计划)项目(2011AA10A101)
通讯作者: 邓其明ꎬ副研究员ꎬ主要从事作物遗传育种研究ꎻE ̄mail:dengqmsc@163.com
李 平ꎬ教授ꎬ主要从事作物遗传育种研究ꎻE ̄mail:liping6575@163.com
共同第一作者: 颜学海ꎬ四川眉山人ꎬ硕士研究生ꎬ主要从事分子植物病理研究ꎻE ̄mail: yanxuehai2012@163.comꎮ
doi:10.13926 / j.cnki.apps.2016.01.011
基于毒性相关基因序列的稻瘟病菌
群体遗传多样性分析
邓其明∗ꎬ颜学海#ꎬ何建美ꎬ邓元宝ꎬ韩 冬ꎬ杨 阳ꎬ李 彬ꎬ朱 军ꎬ 李 平∗
(四川农业大学水稻研究所ꎬ温江 611130)
摘要:选用 16对毒性相关基因特异性引物对四川和重庆 9个县(市)分离到的 200个稻瘟病菌单孢菌株进行 PCR扩增ꎬ并
采用最长距离法进行聚类分析ꎬ结果显示各引物均能扩增出其目的条带ꎬ多态位点百分率(P)高达 93.75%ꎬ扩增频率差异
较大ꎻ200个菌株可归为 70个不同的单元型ꎬ其中单元型 SCH13 为优势单元型ꎻ在 0.86 遗传相似水平上ꎬ200 个菌株可划
分为 27个遗传宗谱ꎬ包括 1个优势宗谱ꎬ3个亚优势宗谱ꎬ14个次要宗谱ꎬ9个小宗谱ꎬ层次丰富ꎻ在群体平均水平上ꎬ病菌
群体具有丰富的遗传多样性(H= 0.324 4ꎬI= 0.484 2)ꎬ且群体间差异较大ꎻ9 个种群在遗传距离为 0.05 水平上可分为 4 个
类群ꎬ种群遗传谱系与地理区域分布呈一定相关性ꎮ 同时ꎬ该地区的群体存在一定的遗传分化(HT = 0.320 0)ꎬ群体内多样
性大于群体间多样性(Hs= 0.179 6ꎬDst= 0.140 4)ꎬ总遗传变异的 56.13%存在于群体内(Gst = 0.438 7)ꎬ群体间基因流动性
较小(Nm= 0.639 6)ꎮ 本研究揭示了四川和重庆部分区域稻瘟病菌群体遗传结构、遗传多样性及其与地理分布之间的关
系ꎬ为抗病育种和品种布局奠定了基础ꎮ
关键词:稻瘟病菌ꎻ致病基因ꎻ无毒基因ꎻDNA指纹ꎻ遗传多样性
Analysis of genetic diversity among populations of rice blast fungus based on the
fingerprint of virulence gene DENG Qi ̄mingꎬ YAN Xue ̄haiꎬ HE Jian ̄meiꎬ DENG Yuan ̄baoꎬ
HAN Dongꎬ YANG Yangꎬ LI Binꎬ ZHU JunꎬLI Ping (Rice Research Institute of Sichuan Agricultural Universityꎬ
Wenjiang 611130ꎬ China)
Abstract: Two hundred single spore isolates of Magnaporthe oryzea from nine counties / cities in Sichuan and
Chongqing were amplified with 16 pairs of virulence gene ̄specific primers. The results showed that the target
bands could be clearly amplified with all the primersꎬ and the percentage of polymorphic loci (P) reached
93.75%. Those isolates were classified into 70 different haplotypes and 27 genetic lineages at the genetic similarity
coefficient level of 0.86 by cluster analysis. Among the 70 haplotypesꎬ SCH13 was the dominant haplotype.
The twenty ̄seven genetic lineages included 1 advantage lineageꎬ 3 subdominant lineagesꎬ 14 secondary lineages and
9 minor lineages. At population levelꎬ it was showed that there existed rich genetic diversity among 200
isolatesꎬ and the diversity existed among fungal populations from different geographic regions (H=0.324 4ꎬI=
0.484 2) . The 9 populations were clustered into 4 groups at 0. 05 genetic distance levelꎬ and the genetic
lineages of the populations from different areas are related to their geographic distributions. Meanwhileꎬ certain
genetic differentiation existed in the same fungal population. The genetic diversity of within ̄population was
higher than that of among ̄population (Hs = 0.179 6ꎬDst = 0.140 4)ꎬ and accounted for 56.13% of the total
variation existed in within ̄population (Gst = 0.438 7) . The gene flow was less in within ̄population (Nm =
植物病理学报 46卷
0.639 6) . The study reveals the genetic structureꎬ genetic diversity of M. oryzea populations and the
relationship of the geographic distributions and the fungal populations from Sichuan and Chongqingꎬ which
could provide the basis for rice blast resistance breeding and rational arrangement of rice varieties in these regions.
Key words: Magnaporthe griseaꎻ pathogenicity geneꎻ avirulence geneꎻ DNA fingerprintingꎻ genetic diversity
中图分类号: S435.111.41 文献标识码: A 文章编号: 0412 ̄0914(2016)01 ̄0091 ̄12
由稻瘟病菌 Magnaporthe oryzae (无性态:
Pyricularia oryzae)引起的稻瘟病是世界水稻生产
中三大毁灭性病害之一ꎬ对水稻的产量和品质均造
成严重破坏[1~3]ꎬ已成为水稻高产、稳产的主要障
碍ꎮ 在西南地区ꎬ稻瘟病是一季稻主要流行性病害
之一ꎬ四川、重庆等地稻瘟病发生率极高ꎮ 实践证
明ꎬ利用品种自身的抗性控制该病的发生和危害是
目前最为经济有效的防治方法ꎮ 但是ꎬ稻瘟菌毒性
相关基因的频繁变异使得病原菌小种极不稳定ꎬ常
导致新培育出的抗稻瘟病品种在推广 3~ 5 年后抗
性逐渐减弱ꎬ甚至丧失[4ꎬ5]ꎮ 由于水稻寄主与稻瘟
病菌之间的互作符合 Flor 的“基因对基因”学说ꎬ
它们之间存在典型的特异性识别和互作的规
律[6~8]ꎮ 因此ꎬ开展稻瘟病菌群体遗传多样性的研
究、掌握病原菌群体变化规律及其毒性信息ꎬ对于
品种选育和合理布局以及延长抗病品种的使用寿
命都具有重要的意义[9~12]ꎮ 研究稻瘟病菌群体结
构的传统方法主要是通过致病型鉴定来实现的ꎮ
该方法在抗病育种上具有一定的应用价值ꎬ但易受
环境或其他因素的影响ꎬ因此ꎬ其结果难以准确地
反映菌株内在的变异性[13ꎬ14]ꎮ 目前ꎬ国内外有许
多学者应用分子标记技术研究稻瘟病菌群体遗传
多样性ꎬ但结果大多显示稻瘟菌群体遗传多样性与
病原菌致病型之间无明显相关性[15ꎬ16]ꎮ 因此ꎬ有
必要寻找出一种新的方法ꎬ能探究病原菌群体的遗
传多样性与毒性相关性ꎮ 有学者提出ꎬ通过分析无
毒基因在基因组中的多态性研究稻瘟病菌群体的
遗传结构变异是一种新方法ꎬ其结果可能会揭示稻
瘟菌群体遗传结构与致病型的相关性[17]ꎮ
经过世界各国科学家近 10 年的努力ꎬ目前已
克隆和分析了近 40个与稻瘟病菌生活史中各个阶
段相关的功能基因[18]ꎮ 其中ꎬ致病基因 MPG1ꎬ
MPS1ꎬMagBꎬ CPKAꎬMNH1ꎬMgYCA1ꎬMgS11 和
MgATG5在稻瘟病菌的发病过程中是非常重要的ꎮ
如 CPKA 编码蛋白激酶 A 的催化亚基 PKA ̄cꎬ是
稻瘟病菌致病性所必需的基因[19]ꎻMPG1 是稻瘟
病菌重要的疏水蛋白基因ꎬ主要功能是参与稻瘟病
菌对外界信号的识别[20]ꎻMPS1 不但影响着附着
胞的侵入过程ꎬ还与稻瘟病菌的营养生长、产孢和
子囊形成有关[21]ꎮ 目前ꎬ已有近 85 个水稻稻瘟病
抗性基因被鉴定ꎬ而对稻瘟病菌无毒基因研究较
少ꎬ已成功克隆的无毒基因仅有几个[22]ꎬ其中ꎬ
Avr ̄piaꎬ Avr ̄pitꎬ Avr ̄pikꎬ Avr ̄pitaꎬ Avr ̄piztꎬ Avr ̄
co39ꎬPRE1和 ACE1 在参与寄主与病原菌互作过
程中是非常重要的ꎮ 如无毒基因 Avr ̄pita编码 223
个氨基酸的金属蛋白酶ꎬ能与抗性基因 Pi ̄ta 直接
互作ꎬ激活特异性免疫反应[23]ꎻ再如含有 AVR1 ̄
CO39的稻瘟病菌可以特异的激发有抗病基因 Pi ̄
CO39( t) 水稻的防卫反应[24]ꎮ
本研究通过选取上述 16个已克隆或鉴定的具
有重要功能的稻瘟病菌毒性相关基因的特异性引
物ꎬ对 200个稻瘟病菌菌株进行 PCR检测ꎬ分析病
菌群体的遗传结构和遗传多样性、以及群体间的亲
缘关系和遗传分化关系ꎬ以期在毒性相关基因遗传
信息的基础上掌握稻瘟病菌群体遗传变化规律ꎬ为
抗病品种的合理布局和使用有效期的延长提供一
定的参考ꎮ
1 材料与方法
1.1 供试菌株
供试菌株来源于四川雅安、温江、南充、自贡、
宜宾、绵阳以及重庆永川、垫江、梁平水稻主产区ꎬ
由四川农业大学水稻研究所近年来对这些稻瘟病
易发区进行菌株的收集、分离及稻草保存ꎬ合计
200个单孢菌株(表 1)ꎮ
1.2 菌丝体培养
将供试菌株活化后再纯化一次ꎬ挑取适量菌丝
在 PDA液体培养基中于转速 120 rmin ̄1摇床上
28℃培养 3~6 dꎬ待菌丝体生长至一定量后用灭菌
纱布过滤收集、烘干ꎬ ̄20℃保存备用ꎮ
29
1期 邓其明ꎬ等:基于毒性相关基因序列的稻瘟病菌群体遗传多样性分析
Table 1 Magnaporthe oryzae isolates tested in this study
Origin Isolate No. of isolates
Yaan 8、12、14、19、25、29、45、21 ̄1 ̄3、44 ̄2 ̄3、95 ̄82 ̄1、62 ̄1 ̄2、62 ̄2 ̄2、91 ̄9 ̄2、91 ̄13 ̄2、97 ̄28 ̄1、
97 ̄41 ̄1、97 ̄51 ̄1、 97 ̄53 ̄1、 97 ̄63 ̄2、 95 ̄83 ̄1、H8 ̄1、H8 ̄2、D56 ̄1、D56 ̄2、 E57 ̄1、 E57 ̄2、
LJ84 ̄121、lJ86 ̄12、LJ94 ̄5、LJ95 ̄1、LJ95 ̄2、LJ96 ̄2、LJ96 ̄3、LJ98 ̄2、LJ98 ̄4、LJ99 ̄2、LJ99 ̄
7、LJ100 ̄1、LJ100 ̄5、LJ101 ̄5、LJ104 ̄8、LJ104 ̄9、LJ105 ̄5、LJ106 ̄6、LJ106 ̄12、LJ107 ̄8、
LJ107 ̄11、LJ109 ̄3、115 ̄6、115 ̄10、116 ̄1、117 ̄12
52
Wenjiang C29 ̄1、C29 ̄2、WJ ̄1、WJ ̄2、WJ ̄3、WJ ̄4、WJ ̄5、WJ ̄6、WJ ̄7、WJ ̄8、WJ ̄9、WJ ̄10、WJ ̄11、
WJ ̄13、WJ ̄14、WJ ̄15、WJ ̄16、WJ ̄17、WJ ̄19、WJ ̄20、WJ ̄21、WJ ̄22、WJ ̄23、WJ ̄24、WJ ̄
25、WJ ̄26、WJ ̄27、WJ ̄28、WJ ̄29、WJ ̄30、WJ ̄31、WJ ̄32、WJ ̄33、WJ ̄34、WJ ̄36、WJ ̄38、
WJ ̄39
37
Nanchong NC ̄1、NC ̄2、NC ̄3、NC ̄4、NC ̄5、NC ̄6、NC ̄7、NC ̄9、NC ̄10、NC ̄11、NC ̄12、NC ̄13、NC ̄
14、NC ̄15、NC ̄16、NC ̄18、NC ̄22、NC ̄24、NC ̄25、NC ̄26、NC ̄27、NC ̄28、NC ̄29、NC ̄30、
NC ̄31、NC ̄32、NC ̄33、NC ̄34、NC ̄36、G33 ̄1、G33 ̄2、H34 ̄1、H34 ̄2
33
Zigong B2 ̄1、B2 ̄2、C3 ̄1、C3 ̄2、V22 ̄1、V22 ̄2 6
Yibin E5、Z52 ̄1、Z52 ̄2、Z52 ̄3 4
Mianyang P16、A27、X50 ̄1、X50 ̄2 4
Yongchuan N14、YC ̄2、YC ̄3、YC ̄4、YC ̄5、YC ̄6、YC ̄7、YC ̄8、YC ̄9、YC ̄10、YC ̄22、YC ̄23、YC ̄25、
YC ̄26、YC ̄27、YC ̄28、YC ̄29、YC ̄30、YC ̄31、YC ̄32、YC ̄33、YC ̄34、YC ̄35、YC ̄36、YC ̄39
25
Dianjiang DJ ̄1、DJ ̄2、DJ ̄3、DJ ̄4、DJ ̄5、DJ ̄6、DJ ̄7、DJ ̄8、DJ ̄9、DJ ̄10、DJ ̄12、DJ ̄13、DJ ̄14、DJ ̄15、
DJ ̄16、DJ ̄17、DJ ̄18、DJ ̄19、DJ ̄20、DJ ̄22、DJ ̄23、DJ ̄24、DJ ̄26、DJ ̄27
24
Liangping LP ̄1、LP ̄2、LP ̄3、LP ̄4、LP ̄5、LP ̄6、LP ̄7、LP ̄8、LP ̄9、LP ̄11、LP ̄12、LP ̄14、LP ̄16、LP ̄
18、LP ̄19
15
1.3 病菌基因组 DNA提取
将收集的菌丝体分装于 2 mL 的 Eppendorf 管
中ꎬ加 2颗 5 mm灭菌碳钢珠后迅速加入 65℃预热
的 CTAB提取液ꎬ应用上海净信科技生产的全自
动样品快速研磨仪进行磨样ꎬ使用 CTAB[25]法进
行基因组 DNA 提取ꎬ最后用 RNAase 进行处理后
抽提一次ꎬ ̄20℃保存备用ꎮ
1.4 引物的选择及 PCR检测
参照 Zhao[26]所选用的稻瘟病菌致病相关基
因 MPG1、MPS1、MagB、CPKA、MNH1、MgYCA1、
MgS11、MgATG5和 Li等[27]]所选用的稻瘟病菌无
毒基因 Avr ̄pia、Avr ̄pit、Avr ̄pik、Avr ̄pita、Avr ̄pizt、
Avr ̄co39、ACE1特异性引物以及 Zhang 等[28]所选
用的与无毒基因 PRE1 紧密连锁的 SSR 标记引物
共 16对ꎮ 以上引物均由英潍捷基(上海)贸易有
限公司合成ꎬ通过 PCR 扩增检测特异性引物对不
同来源菌株的扩增情况ꎮ
1.5 数据统计与分析
根据 PCR 扩增反应的结果ꎬ建立基于稻瘟病
菌毒性相关基因指纹的“0 ~ 1”数据库ꎮ 根据每个
样品电泳条带的有无ꎬ有赋值为“1” ꎬ无赋值为
“0”ꎮ 统计各引物的扩增频率ꎬ并应用 DPS数据处
理系统的 Nei & Li(1979)方法[29]计算遗传距离ꎬ
用最长距离法对 16 对引物的扩增情况进行聚类ꎬ
分析稻瘟病菌群体遗传结构ꎮ 利用 POPGENE 32
软件对群体进行遗传多样性进行分析ꎬ计算有效等
位基因数 (Ne)、 Nei s 基因多样性指数 (H)、
Shannon信息指数( I)、多态性位点数(NP)、多态
位点百分率(P)、群体遗传多样性均值(HT)、群体
内遗传多样性均值(Hs)、群体间遗传多样性均值
(Dst)、遗传分化系数均值(Gst)、基因流(Nm)等
值ꎮ 根据 Neis遗传距离ꎬ利用 MEGA 5.1 软件的
UPGMA法分析群体间的遗传分化关系ꎮ
39
植物病理学报 46卷
2 结果与分析
2.1 供试菌株的 PCR扩增结果
2.1.1 致病基因特异性引物扩增结果 8 对致病
基因的特异性引物均能扩增出符合目的片段大小
的特异性条带ꎬ部分代表性结果如图 1所示ꎮ 各特
异性引物扩增频率差异较大ꎮ 其中ꎬ致病基因
MPS1特异性引物扩增频率最高(92.00%)ꎬMNH1
特异性引物扩增频率最低(10.50%)ꎬ其余 6 个扩
增频率在 65%~90%之间(图 2)ꎮ
Fig. 1 PCR amplification of some isolates by the four pairs of specific primers
M:Markerꎻ CK: Negative contro1 (ddH2O)ꎻ 1~23: DNA samples of blast fungus isolates.
Fig. 2 The PCR amplification frequency of different pathogenic gene ̄specific primers
49
1期 邓其明ꎬ等:基于毒性相关基因序列的稻瘟病菌群体遗传多样性分析
2.1.2 无毒基因特异性引物扩增结果 8 对无毒
基因的特异性引物扩增结果显示均能扩增出符合
目的片段大小的指纹条带ꎬ部分代表性结果如图 3
所示ꎮ 各特异性引物扩增频率差异较大ꎮ 其中ꎬ无
毒基 因 Avr ̄pia 特 异 性 引 物 扩 增 频 率 最 高
(98.50%)ꎻAvr ̄pik、Avr ̄pizt、PRE1、ACE1等特异性
引物扩增频率均在 95.00%以上ꎻAvr ̄pit 特异性引
物扩增频率最低(58.50%)(图 4)ꎮ
Fig. 3 The PCR results of some isolates by four pairs of specific primers
M: Markerꎻ CK: Negative contro1(ddH2O)ꎻ 1~23: DNA samples of blast fungus isolates.
Fig. 4 The PCR amplification frequency of different avirulence gene ̄specific primers
59
植物病理学报 46卷
2.2 稻瘟病菌群体遗传结构
2.2.1 稻瘟病菌群体单元型分析 根据 PCR扩增
结果建立基于毒性相关基因序列指纹条带的 0 ~ 1
数据库ꎬ16 个毒性相关基因 MPG1、MPS1、MagB、
CPKA、MNH、 MgYcA1、 MgS11、 MgATG5、 Avr ̄pia、
Avr ̄pit、Avr ̄pik、Avr ̄pita、Avr ̄pizt、Avr ̄co39、PRE1、
ACE1特异性指纹片段分别用 a、b、c、d、e、f、g、h、i、
j、k、l、m、n、o、p表示ꎬ将扩增出的 DNA带型完全相
同的菌株归为一个单元型ꎬ结果显示 200 个供试菌
株可以归为 70个不同的单元型(表 2)ꎮ 其中单元
型 SCH13包含 55个菌株ꎬ占供试菌株的 27.50%ꎬ此
类型 DNA谱带出现的频率最高ꎬ为优势单元型ꎻ其
次是单元型 SCH15 和 SCH36ꎬ分别包含 12 和 8 个
菌株ꎬ所占供试菌株比例分别为 6.00%和 4.00%ꎻ其
余单元型菌株数较少ꎬ所占比例均在 4.00%以下ꎬ说
明供试菌株在群体遗传结构上层次丰富ꎬ较为复杂ꎮ
Table 2 Haplotypes of Magnaporthe oryzae isolates
No. of
haplotypes
Haplotype
No. of
isolates
Percentage
/ %
No. of
haplotypes
Haplotype
No. of
isolates
Percentage
/ %
SCH01 aaeijklp 3 1.5 SCH36 bfghijklmnop 8 4
SCH02 aeijklmp 2 1 SCH37 bfghijkmnop 1 0.5
SCH03 aeijklmop 1 0.5 SCH38 bfhijklmnop 1 0.5
SCH04 ijklmp 1 0.5 SCH39 bfghikmnop 2 1
SCH05 eijklmp 1 0.5 SCH40 bfghiklmnop 2 1
SCH06 bcdeghijklmnop 1 0.5 SCH41 abfghijkmnop 1 0.5
SCH07 bcdefghijklmnop 1 0.5 SCH42 abcdfghijkmnop 1 0.5
SCH08 abcdefghijklmop 1 0.5 SCH43 abefhijklmop 1 0.5
SCH09 bcdghijkmop 1 0.5 SCH44 bfghiklno 1 0.5
SCH10 bcdefgiklmop 1 0.5 SCH45 abfghiklmnop 3 1.5
SCH11 bcdfghikmop 1 0.5 SCH46 abdfghiklmnop 1 0.5
SCH12 bcdeghijklmnop 4 2 SCH47 abghiklmnop 1 0.5
SCH13 abcdfghijklmnop 55 27.5 SCH48 abcfghiklmnop 1 0.5
SCH14 abcdefghiklmnop 3 1.5 SCH49 abcefghiklmnop 1 0.5
SCH15 abcdfghiklmnop 12 6 SCH50 abfghikmno 1 0.5
SCH16 abcdfgijklmnop 3 1.5 SCH51 abcfghikmno 1 0.5
SCH17 abcdfgiklmnop 1 0.5 SCH52 abcfghikmnop 2 1
SCH18 bcfghijklmnop 6 3 SCH53 hiklmnop 1 0.5
SCH19 bcdfghijklmnop 1 0.5 SCH54 hikmop 1 0.5
SCH20 abfghijklmnop 7 3.5 SCH55 cdhikmnop 1 0.5
SCH21 abcfghijklmnop 5 2.5 SCH56 bcdfghiklmnop 3 1.5
SCH22 abcfghijklmno 1 0.5 SCH57 bcdfghiklmnop 2 1
SCH23 abcdfghijklmno 1 0.5 SCH58 bcdefgiklmop 1 0.5
SCH24 adfghijklmnop 1 0.5 SCH59 bcdfgikmop 5 2.5
SCH25 acdfghijklmnop 1 0.5 SCH60 bcdfghikmop 6 3
SCH26 bcdfghikmnop 7 3.5 SCH61 abcdfghikmop 2 1
SCH27 abcdfghikmnop 6 3 SCH62 abcdfghiklmop 1 0.5
SCH28 bcdfghijkmnop 1 0.5 SCH63 abcdfgiklmop 1 0.5
SCH29 abcdfghijkmnop 4 2 SCH64 abcdfgikmop 1 0.5
SCH30 abcfghikmnop 1 0.5 SCH65 bdfhikmop 1 0.5
SCH31 bcdghikmnop 1 0.5 SCH66 bcdfhikmop 2 1
SCH32 abcdfghijkmnop 1 0.5 SCH67 bdfgikmnop 1 0.5
SCH33 abcdfhikmnop 1 0.5 SCH68 abcdimop 1 0.5
SCH34 adghijkmop 1 0.5 SCH69 bcdimop 1 0.5
SCH35 adegijklmnop 1 0.5 SCH70 cdfimop 1 0.5
69
1期 邓其明ꎬ等:基于毒性相关基因序列的稻瘟病菌群体遗传多样性分析
2.2.2 稻瘟病菌群体遗传宗谱分析 对所有单元
型进行聚类分析ꎬ结果显示 200 个菌株在 0.17 遗
传相似水平上可划分为 1 个遗传宗谱ꎻ在 0.50 遗
传相似水平上可划分为 5个遗传宗谱ꎬ包括 1个绝
对优势宗谱、2 个亚优势宗谱和 2 个次要宗谱ꎻ在
0.76遗传相似水平上可划分为 14 个遗传宗谱ꎬ包
括 1个绝对优势宗谱、3个亚优势宗谱、8个次要宗
谱和 2 个小宗谱ꎻ在 0.86 遗传相似水平下可划分
为 27个宗谱ꎬ包括 1个绝对优势宗谱、3 个亚优势
宗谱、14个次要宗谱和 9 个小宗谱ꎮ 可以看出ꎬ遗
传相似水平越高ꎬ供试菌株划分得越细ꎬ宗谱层次
越丰富(表 3)ꎮ
在本研究中以 0.86遗传相似水平下分析稻瘟
病菌群体结构特征层次比较丰富(表 4)ꎮ 在 0.86
遗传相似水平下:宗谱 SCL5 包含 78 个菌株ꎬ占供
试菌株的 39.00%ꎬ为优势宗谱ꎻ宗谱 SCL6、8、22
分别包含 21、20、17 个菌株ꎬ属亚优势宗谱ꎻ宗谱
SCL1、2、3、4、7、9、12、13、14、17、18、23、24、26 分别
包含 2 ~ 10 个菌株ꎬ为次要宗谱ꎻ其余 9 个宗谱均
只包含 1个菌株ꎬ为小宗谱ꎮ
Table 3 Analysis of population genetic lineages of Magnaporthe oryzae at
different genetic similar levels
Genetic
similar level
No. of
lineages
Absolute
advantage lineage
Subdominant
advantage lineage
Secondary lineage Minor lineage
Lineage
No. of
isolates
Lineage
No. of
isolates
Lineage
No. of
isolates
Lineage
No. of
isolates
0.17 1 SCL1 200 — — — — — —
0.50 5 SCL2 131 SCL3、5 29 SCL1、4 3~8 — —
0.76 14 SCL4 101 SCL5、7、12 16~22
SCL1、 2、 3、 6、
10、11、13、14
2~11 SCL8、9 1
0.86 27 SCL5 78 SCL6、8、22 17~21
SCL1、2、3、4、7、
9、12、13、14、17、
18、23、24、26
2~10
SCL10、 11、
15、 16、 19、
20、 21、 22、
25
1
Table 4 Genetic lineages of Magnaporthe oryzae isolates at the
genetic similarity coefficient level of 0.86
Lineage No. of isolates Percentage / % Lineage No. of isolates Percentage / %
SCL1 6 3.00 SCL15 1 0.50
SCL2 2 1.00 SCL16 1 0.50
SCL3 3 1.50 SCL17 7 3.50
SCL4 3 1.50 SCL18 4 2.00
SCL5 78 39.00 SCL19 1 0.50
SCL6 21 10.50 SCL20 1 0.50
SCL7 2 1.00 SCL21 1 0.50
SCL8 20 10.00 SCL22 17 8.50
SCL9 2 1.00 SCL23 5 2.50
SCL10 1 0.50 SCL24 3 1.50
SCL11 1 0.50 SCL25 1 0.50
SCL12 10 5.00 SCL26 2 1.00
SCL13 4 2.00 SCL27 1 0.50
SCL14 2 1.00
79
植物病理学报 46卷
2.2.3 稻瘟病菌遗传结构的地理分布特点 根据
稻瘟病菌毒性相关基因指纹特征ꎬ对不同来源菌株
遗传结构的地理分布特点进行统计分析(表 5)ꎮ
从表 5可以看出ꎬ来源于四川 136个稻瘟病菌株可
划归为 49 个单元型ꎬ分属 20 个遗传宗谱ꎻ而来源
于重庆的 64个稻瘟病菌株可划归为 26个单元型ꎬ
分属 14个遗传宗谱ꎮ 可以看出ꎬ9 个区域稻瘟病
菌遗传宗谱非常丰富ꎬ除永川未检测到宗谱 SCL5
外ꎬ其余 8个区域均检测到相同的遗传宗谱 SCL5ꎬ
同时各区域又具有自己的特异性宗谱ꎮ
2.3 西南地区稻瘟病菌群体遗传多样性分析
2.3. 1 稻瘟病菌群体遗传多样性水平 利用
POPGEME32软件对 9 个区域的种群进行遗传参
数分析ꎬ结果显示 9个区域的群体遗传多样性比较
丰富(表 6)ꎮ 在群体平均水平上ꎬ有效等位基因数
目(Ne)为 1.564 2 ꎬNeis 基因多样性指数(H)为
0.324 4ꎬShannon 信息指数( I)为 0.484 2ꎬ多态性
位点数 (NP)为 15 个ꎬ多态位点百分率 ( P)为
93.75%ꎮ 不同区域之间ꎬ各群体的遗传多样性水
Table 5 Geographic distribution of the genetic diversity of Magnaporthe oryzae population
Origin No. of isolates Haplotype Genetic lineage
Sichuan Yaan 52 2、8、9、10、13、14、15、16、19、21、24、25、26、
28、34、43、56、57、58
1、3、4、5、6、7、8、10、15、22
Wenjiang 37 7、12、18、20、36、37、38、39、40、41、44、45、
50、51
3、5、6、12、13、14、16、17、18、19
Nanchong 33 1、12、13、15、21、22、29、30、32、35、42、46、
47、48、49、52
1、5、6、8、9、11、14、17、18
Zigong 6 1、2、3、4、6、12 1、2、3、5
Yibin 4 13 5
Mianyang 4 13、16 5
Chongqing Yongchuan 25 5、11、26、27、54、59、60、61、63、64、65、66、
68、69、70
2、4、8、9、20、22、23、24、26、27
Dianjiang 24 14、15、17、26、27、31、33、55、62、66、67 5、8、21、22、23、24、25
Liangping 15 13、15、23、27、29 5、6、8
Table 6 Genetic diversity within 9 populations of Magnaporthe oryzae in Southwest region
Population aNa bNe cH dI eNP fP / %
Zigong 1.687 5 1.392 9 0.239 1 0.362 1 11 68.75
Mianyang 1.062 5 1.044 2 0.025 9 0.037 8 1 6.25
Yibin 1.000 0 1.000 0 0.000 0 0.000 0 0 0.00
Wenjiang 1.687 5 1.324 7 0.202 0 0.312 3 11 68.75
Yaan 1.750 0 1.467 9 0.275 9 0.410 4 12 75.00
Nanchong 1.812 5 1.582 0 0.329 5 0.480 5 13 81.25
Liangping 1.187 5 1.168 0 0.085 7 0.121 6 3 18.75
Yongchuan 1.812 5 1.438 8 0.254 5 0.382 8 13 81.25
Dianjiang 1.562 5 1.341 5 0.204 2 0.304 3 9 56.25
Population average level 1.937 5 1.564 2 0.324 4 0.484 2 15 93.75
ANa: Observed number of allelesꎻ BNe: Effective number of allelesꎻ CH: Neis gene diversity indexꎻ DI: Shannons informa ̄
tion indexꎻ ENP: Number of polymorphic lociꎻ FP: Percentage of polymorphic loci.
89
1期 邓其明ꎬ等:基于毒性相关基因序列的稻瘟病菌群体遗传多样性分析
平具有一定差异ꎮ 来源于南充的群体具有较高的
遗传多样性(H = 0.329 5ꎬI = 0.480 5)ꎬ其次是雅
安、永川、自贡相对较高ꎬ而温江和垫江相对较低ꎬ
绵阳、宜宾和梁平最低ꎮ
2.3.2 稻瘟病菌群体遗传分化 通过遗传参数的
分析ꎬ结果显示四川和重庆稻瘟病菌群体存在一定
的遗传分化ꎮ 来自 9 个区域的稻瘟病菌群体遗传
多样性(HT)均值为 0.320 0ꎬ群体内遗传多样性
(Hs)均值为 0.179 6ꎬ群体间遗传多样性(Dst)均
值为 0.140 4ꎬ可以看出ꎬ这些区域稻瘟病菌群体内
遗传多样性大于群体间遗传多样性ꎮ 遗传分化系
数(Gst)均值为 0.438 7ꎬ表明这些区域稻瘟病菌总
遗传变异的 43.87%存在于群体间ꎬ56.13%存在于
群体内ꎮ 基因流(Nm)值为 0.639 6ꎬ说明四川和重
庆稻瘟病菌群体间基因流动性较小ꎮ
2.4 西南地区稻瘟病菌群体间的聚类分析
根据 Neis遗传距离ꎬ应用 MEGA5.10 软件构
建 8个稻瘟病菌种群的 UPGMA 系统树(图 5)ꎮ
在遗传距离为 0.05水平上ꎬ9个种群可分为 4个类
群ꎬ其中自贡、温江单独构成一个类群ꎬ湄潭、宜宾、
梁平、雅安、南充 5个地区的种群构成一个类群ꎬ而
南充和垫江地区的种群构成另一个类群ꎻ在遗传距
离为 0.10水平上ꎬ9 个种群可以分为 3 个类群ꎻ而
遗传距离在 0.15水平上只能分为 2 个类群ꎮ 可以
看出ꎬ来源于自贡的菌株与其余 8个区域的菌株在
遗传结构上差异较大ꎬ但总体上供试稻瘟病菌种群
大都聚为一类ꎬ亲缘关系较近ꎬ同时也反映出种群
遗传谱系与地理区域分布呈现一定相关性ꎮ
3 结论与讨论
稻瘟病菌在长期的进化过程中ꎬ其群体遗传结
构随品种、地理环境、栽培方式等因素变化而发生
变异ꎬ病原菌毒性相关基因及其致病型的组成和变
异直接关系到稻瘟病的发生和流行[30]ꎮ 因此ꎬ在
传统的植物病理方法研究病原菌群体结构的基础
上ꎬ致力于群体遗传结构多样性的研究对于合理布
局、延长抗病品种使用寿命等方面具有重要意义ꎮ
以菌株 DNA 带的位置相似水平为标准进行
聚类ꎬ将菌株划分成不同的密切相关的组ꎬ每一个
组称为一个遗传宗谱ꎬ以此对稻瘟病菌群体进行描
述[31~33]ꎮ 本研究结果表明 16 对毒性相关基因引
物均能扩增出目的条带ꎬ多态位点百分率高达
93.75%ꎬ说明选用的引物具有较高的多态性ꎬ可以
用于稻瘟病菌群体遗传分析ꎮ 聚类分析表明 200
个菌株可以归为 70 个不同的单元型和 27 个遗传
宗谱ꎬ供试菌株在群体遗传结构上出现了许多次要
小宗谱ꎬ层次十分丰富ꎬ且区域内的菌株在遗传结
构上变异度非常大ꎮ Li 等[34]采用 Pot2 ̄rep ̄PCR
指纹识别技术将宁夏 406 个稻瘟病菌菌株鉴定为
65个单元型ꎻTong 等[29]利用 13 对 SSR 引物将来
源于湖南 19 个县(市)的 169 个稻瘟病单孢菌株
在 0.80相似水平上划分为 8个遗传宗谱ꎻWu[35]等
Fig. 5 UPGMA dendrogram among populations of Magnaporthe oryzae from Sichuan and Chongqing
99
植物病理学报 46卷
通过基于 SRAP标记的分子指纹分析法对来源于
广东省和云南省的菌株进行遗传结构分析ꎬ在相似
系数取 0.83 时将广东的 40 个菌株划分为 9 个宗
谱ꎬ云南的 40 个菌株划分为 16 个宗谱ꎮ 可以看
出ꎬ四川和重庆等西南地区稻瘟病菌的群体遗传结
构较其他区域更具明显的优势宗谱ꎬ同时还存在着
许多特异性的小宗谱ꎬ层次比较丰富ꎬ这可能与这
些区域适温、高湿、寡照的气象条件、复杂的生态地
理环境及寄主多样化等因素有关ꎮ
生物群体的遗传多样性是评价生物资源状况
的重要依据之一ꎬNei s 基因多样性指数(H)和
Shannon信息指数( I)是衡量生物遗传多样性的重
要参数[36]ꎮ 本研究中ꎬ西南地区遗传多样性具有
较高水平(H=0.324 4ꎬI = 0.484 2)ꎬ且群体间的遗
传多样性水平存在差异ꎮ 来源于南充的菌株在群
体水平上基因多样性指数(H)和 Shannon 信息指
数( I)最高ꎬ说明南充的病菌群体发生变异的潜力
最大ꎬ其次是雅安、永川、自贡相对较高ꎬ应加强对
这些区域病菌群体的监测ꎬ掌握其动态变化规律ꎬ
为科学的品种布局、有效防止稻瘟病提供更充分的
依据ꎮ
在群体遗传分化水平上ꎬ四川和重庆稻瘟病菌
群体内多样性(Hs = 0.179 6)大于群体间多样性
(Dst=0.140 4)ꎮ Wright[37]提出群体间遗传分化
系数在 0.05~ 0.15 之间表示遗传分化程度居于中
等水平ꎻ当群体间基因流大于 1 时ꎬ群体间就存在
一定的基因流动ꎬ能发挥匀质化作用ꎮ 本研究 9 个
稻瘟病菌群体间的遗传分化系数均值为 0.438 7ꎬ
说明四川和重庆稻瘟病菌不同群体间存在较大的
遗传分化ꎬ总遗传变异的 56.13% 存在于群体内ꎬ
进一步表明群体内变异和群体间变异均是这些区
域稻瘟病菌遗传变异的来源ꎮ 不同群体间的基因
流为 0.639 6ꎬ说明区域病菌群体间基因交流较少ꎬ
可能与各区域病菌群体的特殊地理位置有关ꎮ
本研究应用 16对毒性相关基因特异性引物对
200个稻瘟病菌株进行 PCR 检测ꎬ分析其群体遗
传多样性ꎬ在病菌毒性信息的基础上揭示了 9个区
域稻瘟病菌群体遗传规律及其与地理分布之间的
关系ꎬ为抗病育种和品种布局奠定了基础ꎮ
参考文献
[1] Couch B Cꎬ Hohn L M. A multilocus gene genealogy
concordant with host preference indicates segregation of
a new speciesꎬ Magnaporthe oryzaeꎬ from M. grisea
[J] . Mycologiaꎬ 2002ꎬ 94: 683 ̄693.
[2] Jeong J Sꎬ Mitchell T Kꎬ Dean R A. The Magnaporthe
grisea snodprot1 homologꎬ MSP1ꎬ is required for viru ̄
lence [ J] . FEMS Microbiology Lettersꎬ 2007ꎬ 273:
157 ̄165.
[3] Ashikawa Iꎬ Hayashi Nꎬ Yamane Hꎬ et al. Two adja ̄
cent nucleotide ̄binding site ̄leucine ̄rich repeat class
genes are required to confer Pikm ̄specific rice blast re ̄
sistance [J] . Geneticsꎬ 2008ꎬ 180(4): 2267 ̄2276.
[4] Li Yꎬ Liu E Mꎬ Dai L Yꎬ et al. Genetic diversity a ̄
mong populations as related to pathotypes for Magna ̄
porthe grisea in Hunan province ( in Chinese) [ J] .
Chinese Journal of Rice Science (中国水稻科学)ꎬ
2007ꎬ 21(3): 304 ̄308.
[5] Zhou J Hꎬ Wang J Lꎬ Jiang W R ꎬ et al. Virulence
genes diversity and geographic distribution of
Pyricularia grisea ( in Chinese) [J] . Acta Agronomica
Sinica (作物学报)ꎬ 2003ꎬ 29(5): 646 ̄651.
[6] Wang S Wꎬ Zheng W Jꎬ Zhao J Mꎬ et al.
Identification and analysis of Magnaporthe oryzae
avirulence genes in Liaoning province ( in Chinese)
[J] . Scientia Agricultura Sinica (中国农业科学)ꎬ
2014ꎬ 47(3): 462 ̄472.
[7] Liu J Lꎬ Wang X Jꎬ Mitchell Tꎬ et al. Recent progress
and understanding of the molecular mechanisms of the
rice Magnaporthe oryzae interaction [ J] . Molecular
Plant Pathologyꎬ 2010ꎬ 11(3): 419 ̄427.
[8] Flor H H. Current status of the gene ̄for ̄gene concept
[ J ] . Annual Review of Physiologyꎬ 1971ꎬ 9:
275 ̄296.
[9] Lei C Lꎬ Wang J Lꎬ Jiang W Rꎬ et al. Study on
pathologic races and virulence of blast fungus and their
movement in japonica rice ̄growing region of northern
China ( in Chinese) [J] . Acta Agronomica Sinica (作
物学报)ꎬ 2000ꎬ 26(6): 769 ̄776.
[10] Dong L Yꎬ Wang Qꎬ Liu S Fꎬ et al. Pathogenicity a ̄
nalysis of Magnaporthe oryzae populations of Yunnan
on monogenic lines for resistance to rice blast ( in Chi ̄
nese) [ J] . Southwest China Journal of Agricultural
Sciences (西南农业学报)ꎬ 2012ꎬ 25(2): 467 ̄473.
001
1期 邓其明ꎬ等:基于毒性相关基因序列的稻瘟病菌群体遗传多样性分析
[11] Mekwatanakarn Pꎬ Kositratana Wꎬ Levy Mꎬ et al.
Pathotype and virulence gene diversity of Pyricularia
grisea in Thailand as determined by rice lines near iso ̄
genic for major resistance genes [ J] . Plant Diseaseꎬ
2000ꎬ 84: 60 ̄70.
[12] Xiao D Fꎬ Zhang P Sꎬ Wang Lꎬ et al. Research pro ̄
gress on populations and physiological rice distribution
of rice blast pathogen (Magnaporthe grisea) in China
( in Chinese) [J] . Chinese Journal Rice Sciences (中
国水稻科学)ꎬ 2013ꎬ 27(3): 312 ̄320.
[13] Chen D Hꎬ Zeigler R Sꎬ Leung Hꎬ et al. Population
structure of Pyricularia grisea at two screening sites in
the Philippines [ J ] . Phytopathologyꎬ 1995ꎬ 85:
1011 ̄1020.
[14] Han S Sꎬ Ra D Sꎬ Cheoi S Hꎬ et al. Population dy ̄
namics of Pyricularia grisea during leaf and panicle
blast stages in the same field [ J] . Korean Journal of
Plant Pathologyꎬ 1993ꎬ 13: 408 ̄415.
[15] Javan ̄Nikkhah Mꎬ McDonald B Aꎬ Banke Sꎬ et al.
Genetic structure of Iranian Pyricularia grisea popula ̄
tions based on rep ̄PCR fingerprinting [ J] . European
Journal of Plant Pathologyꎬ 2004ꎬ 110(9): 909 ̄919.
[16] Park S Yꎬ Milgroom M Gꎬ Han S Sꎬ et al. Diversity
of pathotypes and DNA fingerprint haplotypes in popu ̄
lations of Magnaporthe grisea in Korea over two deca ̄
des[J] . Phytopathologyꎬ 2003ꎬ 93(11): 1378 ̄1385.
[17] Mu H Mꎬ Jiang Hꎬ Wang Y Lꎬ et al. Methodology of
genetic diversity research on rice blast pathogen Mag ̄
naporthe grisea ( in Chinese) [ J] . Chinese Journal
Rice Sciences (中国水稻科学)ꎬ 2013ꎬ 27 ( 5):
545 ̄552.
[18] Clergeot P Hꎬ Gourgues Mꎬ Cots Jꎬ et al. PLS1ꎬ a
gene encoding a tetraspanin ̄like proteinꎬ is required for
penetration of rice leaf by the fungal pathogen Magna ̄
porthe grisea [ J ] . Proceedings of the National
Academy of Sciences. 2001ꎬ 98: 6963 ̄6986.
[19] Zhou X Yꎬ Zhao X Hꎬ Xue C Yꎬ et al. Bypassing
both surface attachment and surface recognition require ̄
ments for appressorium formation by overactive Ras
signaling in Magnaporthe oryzae [J] . Molecular Plant ̄
microbe Interactionsꎬ 2014ꎬ 27(9):996 ̄1004.
[20] Rey A Aꎬ Hocher Aꎬ Kwan A Hꎬ et al. Backbone and
sidechain(1)Hꎬ (13)C and (15)N chemical shift as ̄
signments of the hydrophobin MPG1 from the rice blast
fungus Magnaporthe oryzae [ J] . Biomolecular NMR
Assignmentsꎬ 2013ꎬ 7(1): 109 ̄112.
[21] Kong L Aꎬ Li G Tꎬ Liu Yꎬ et al. Differences between
appressoria formed by germ tubes and appressorium ̄
like structures developed by hyphal tips in
Magnaporthe oryzae [ J ] . Fungal Genetics and
Biologyꎬ 2013ꎬ 56: 33 ̄41.
[22] Huang Jꎬ Si W Nꎬ Deng Q Mꎬ et al. Rapid evolution
of avirulence genes in rice blast fungus Magnaporthe
oryzae [J] . BMC Geneticsꎬ 2014ꎬ 15: 45.
[23] Chuma Iꎬ Isobe Cꎬ Hotta Yꎬ et al. Multiple transloca ̄
tion of the AVR ̄Pita effector gene among chromosomes
of the rice blast fungus Magnaporthe oryzae and related
species [J] . PLoS Pathogensꎬ 2011ꎬ 7(7): e1002147.
[24] Ribot Cꎬ Césari Sꎬ Abidi Iꎬ et al. The Magnaporthe
oryzae effector AVR1 ̄CO39 is translocated into rice
cells independently of a fungal ̄derived machinery [J] .
Plant Journalꎬ 2013ꎬ 74(1): 1 ̄12.
[25] Guo Mꎬ Chen Yꎬ Du Yꎬ et al. The bZIP transcription
factor MoAp1 mediates the oxidative stress response
and is critical for pathogenicity of the rice blast fungus
Magnaporthe oryzae [ J ] . PLoS Pathogensꎬ 7
(2): e1001302.
[26] Zhao W W. Detection of pathogenicity genes of Mag ̄
napoethe grisea and distribution in Heilongjiang pro ̄
vince ( in Chinese) [D]ꎬ Heilongjiang: Heilongjiang
Bayi Agricultural Universityꎬ 2010.
[27] Li X Xꎬ Wang Q Luo S Xꎬ et al. Analyzing avirulence
genes of Magnaporthe oryzae from Heilongjiang pro ̄
vince and screening rice germplasm with resistance to
blast fungus ( in Chinese ) [ J ] . Acta Agronomica
Sinica (作物学报)ꎬ 2012ꎬ 38(12): 2192 ̄2197.
[28] Zhang Q Fꎬ Jin X Hꎬ C X Xꎬ et al. Identification of
Magnaporthe Oryzae avirulence gene in Heilongjiang
province with SSR marker ( in Chinese) [ J] . Chinese
Agricultural Science Bulletin(中国农学通报)ꎬ 2010ꎬ
26(10): 250 ̄254.
[29] Tong J Xꎬ Ren Z Hꎬ Liu Y Xꎬ et al. Analysis of ge ̄
netic diversity of Magnaporthe oryzae in Hunan ( in
Chinese) [ J] . Hybrid Rice (杂交水稻)ꎬ 2012ꎬ 27
101
植物病理学报 46卷
(3): 66 ̄70.
[30] Liu Z Hꎬ Li X Fꎬ Wang Y Yꎬ et al. Investigation on
population genetic structure of Magnaporthe grisea in
Yunnan ( in Chinese) [J] . Journal of Yunnan Agricul ̄
tural University (云南农业大学学报)ꎬ 2013ꎬ 28(1):
9 ̄15.
[31] Lu D Hꎬ Ye H Lꎬ Liao H Mꎬ et al. Population genetic
structure and its variation of rice blast fungus in long ̄
grained nonglutinous rice growing regions in Sichuan
( in Chinese) [ J] . Acta Phytophylacica Sinica (植物
保护学报)ꎬ 2005ꎬ 32(1): 23 ̄28.
[32] Hamer J E ꎬ Farrall Lꎬ Orbach M Jꎬ et al. Host
species ̄specific conservation of a family of repeated
DNA sequences in the genome of a fungal plant patho ̄
gen [J] . Proceedings of the National Academy of Sci ̄
encesꎬ 1989ꎬ 86: 9981 ̄9985.
[33] Chen D Hꎬ Zeigeler R Sꎬ Leung Hꎬ et al. Population
structure of Pyricularia grisea at two screening sites in
the Philippines [ J] . Phytopathologyꎬ 1995ꎬ 85(9):
1011 ̄1020.
[34] Li W Qꎬ Wang Y Cꎬ Zheng X B. Genetic structure of
Magnaporthe grisea population from rice in Ningxia
Hui Autonomous Regionꎬ China ( in Chinese) [ J] .
Journal of Nangjing Agricultural University (南京农业
大学学报)ꎬ 2008ꎬ 31(4): 66 ̄72.
[35] Wu W Hꎬ Wang Lꎬ Cheng G Zꎬ et al. Studies on mo ̄
lecular genetics of rice blast fungus population ̄compari ̄
son of genetic and pathotypic structures of two rice
blast fungus populations derived from Guangdong and
Yunnan provinces of China ( in Chinese) [J] . Scientia
Agricultura Sinica (中国农业科学)ꎬ 2004ꎬ 37(5):
675 ̄680.
[36] Guo B Yꎬ Zhou Cꎬ Lü Z Mꎬ et al. Genetic diversity of
different geographical populations in Octopus variabilis
revealed by ISSR analysis ( in Chinese) [ J] . Oceano ̄
logia et Limnologia Sinica (海洋与湖沼)ꎬ 2011ꎬ 42
(6): 868 ̄873.
[37] Wright S. Evolution in mendelian populations [J] . Ge ̄
neticsꎬ 1931ꎬ 16(2): 97 ̄159.
责任编辑:李晖
201