免费文献传递   相关文献

Amplifying stability and genetic diversity of microsatellite loci in different Rhopalosiphum padi (Linnaeus) populations in China

我国禾谷缢管蚜微卫星位点扩增稳定性及遗传多样性



全 文 :植物保护学报 Journal of Plant Protection, 2015, 42(3): 297 - 303 DOI: 10􀆰 13802 / j. cnki. zwbhxb. 2015􀆰 03􀆰 002
基金项目: 国家自然科学基金(31272036,31471766),教育部高等学校博士学科点专项科研基金(20110204110001)
∗通讯作者(Author for correspondence), E⁃mail: maohua. chen@ nwsuaf. edu. cn
收稿日期: 2014 - 06 - 30
我国禾谷缢管蚜微卫星位点扩增稳定性及遗传多样性
段辛乐1   乔宪凤2   彭  雄1   左亚运1   刘高鸣1   陈茂华1∗
(1.西北农林科技大学植物保护学院, 农业部西北黄土高原作物有害生物综合治理重点实验室,
陕西 杨凌 712100; 2.西北农林科技大学动物医学院, 陕西 杨凌 712100)
摘要: 为筛选用于我国禾谷缢管蚜种群遗传学研究的微卫星位点,从 8 个省(市)共 9 个地区采集
282 头试虫,采用微卫星 PCR产物荧光标记与自动扫描分型方法,研究了 8 个欧洲禾谷缢管蚜微卫
星位点在试虫个体中的扩增稳定性和遗传多样性。 结果显示:位点 R1􀆰 35 在 9 个种群中均不能扩
增;位点 R5􀆰 29b只在 7 个种群的少数样本中能扩增;位点 R2􀆰 73、R3􀆰 171、R5􀆰 10、R5􀆰 138、R5􀆰 50 和
R6􀆰 3 在各种群中均能稳定扩增,这 6 个位点的无效等位基因频率为 0􀆰 0044 ~ 0􀆰 2663,平均等位基
因数为 2􀆰 9 ~ 9􀆰 3 个,观测杂合度为 0􀆰 047 ~ 0􀆰 912,其中位点 R6􀆰 3 具有较低的观测杂合度(0􀆰 047)
和等位基因数(2􀆰 9),不适合用于种群遗传多样性研究,而位点 R2􀆰 73、R3􀆰 171、R5􀆰 10、R5􀆰 138 和
R5􀆰 50 均具有较高的杂合度和等位基因数,可用于我国禾谷缢管蚜的种群遗传学研究。
关键词: 禾谷缢管蚜;微卫星; 遗传多样性;扩增稳定性;多态性
Amplifying stability and genetic diversity of microsatellite loci in different
Rhopalosiphum padi (Linnaeus) populations in China
Duan Xinle1   Qiao Xianfeng2   Peng Xiong1   Zuo Yayun1   Liu Gaoming1   Chen Maohua1∗
(1. Key Laboratory of Crop Pest Integrated Pest Management on the Loess Plateau of Ministry of Agriculture,
College of Plant Protection, Northwest A & F University, Yangling 712100, Shaanxi Province, China;
2. College of Veterinary Medicine, Northwest A & F University, Yangling 712100, Shaanxi Province, China)
Abstract: In order to screen polymorphic microsatellite loci which can be used for the population genetics
research of Rhopalosiphum padi populations in China, the amplifying stability and genetic diversity of
eight previously reported microsatellite loci were analyzed in 282 individuals sampled in nine regions of
eight provinces in China. The PCR fragments of the microsatellite loci were labeled fluorescently and
scanned automatically. The results showed that the locus R1􀆰 35 could not be amplified in any samples of
the nine populations, and the locus R5􀆰 29b could only be amplified in a few samples of seven
populations. Loci R2􀆰 73, R3􀆰 171, R5􀆰 10, R5􀆰 138, R5􀆰 50 and R6􀆰 3 could be amplified stably in
samples of all populations. The null allele frequency (Na) of the six stably amplified loci was between
0􀆰 0044 and 0􀆰 2663; the mean number of allele (NA ) ranged from 2􀆰 9 to 9􀆰 3, and the observed
heterozygosity ( Ho ) was from 0􀆰 047 to 0􀆰 912. Locus R6􀆰 3 showed lower observed heterozygosity
(0􀆰 047) and mean number of allele (2􀆰 9) than other loci, indicating that it could not be used for
further population genetics analysis in China. Loci R2􀆰 73, R3􀆰 171, R5􀆰 10, R5􀆰 138 and R5􀆰 50 showed
high observed heterozygosity and mean number of allele, suggesting that they are usable for the genetic
diversity and genetic structure research of R. padi populations in China.
Key words: Rhopalosiphum padi; microsatellite locus; genetic diversity; amplifying stability; polymorphism
    禾谷缢管蚜 Rhopalosiphum padi (Linnaeus)是
麦类作物重要的世界性害虫,在我国各小麦种植区
均有分布 (张广学和钟铁森,1983; Simon et al. ,
1991),该虫以刺吸式口器吮吸麦株茎、叶和嫩穗汁
液,影响小麦正常生长发育,麦苗受害轻度时叶色发
黄、生长停滞、分蘖减少,严重时则导致小麦生长停
滞,最后枯黄甚至死亡。 该虫还能够传播大麦黄矮
病毒(Barley yellow dwarf virus,BYDV),每年给世界
各地的小麦生产造成重大经济损失 (张向才等,
1985;费晓东等,2011)。 禾谷缢管蚜主要在我国南
方冬麦区发生为害,由于近年来农田生态条件的急
剧变化,禾谷缢管蚜种群数量骤增,且成灾区呈现向
北方扩散的趋势 (曹雅忠等,2006; Chen & Han,
2006;鲁艳辉等,2009)。
微卫星标记具有多态性丰富、数量多且分布均
匀、共显性遗传等特点 ( Balloux & Lugon⁃Moulin,
2002;Selkoe & Toonen,2006),近年来已被广泛应用
于生物种群遗传结构、遗传分化、基因交流以及种群
微进化等相关研究(Chen & Dorn,2010;冯彦霞等,
2013;Zheng et al. ,2013)。 微卫星标记已被应用于
不同蚜虫遗传结构、生殖模式、不同生殖模式种群间
的遗传差异和蚜虫迁飞规律等研究。 研究发现,欧
洲地区麦长管蚜 Sitobion avenae 种群具有不同的遗
传结构和生殖方式,且与气候条件相关(Papura et
al. ,2003),而英国的麦长管蚜 S. avenae 具有较高
的迁飞性(Llewellyn et al. ,2003)。 通过对欧洲地区
葡萄根瘤蚜 Phylloxera vitifoliae 种群的微卫星标记
研究显示,欧洲地区的葡萄根瘤蚜种群间缺乏基因
交流(Vorwerk & Forneck,2006),而英国的桃蚜 My⁃
zus persicae 种群之间基因交流频繁(Fenton et al. ,
2010)。 Delmotte et al. (2002)利用 7 对微卫星引物
对法国地区不同生殖方式的禾谷缢管蚜遗传结构进
行了研究;Gilabert et al. (2009)运用微卫星标记研
究了环境因素对法国禾谷缢管蚜生殖方式的影响。
目前,关于我国禾谷缢管蚜种群遗传学的微卫星标
记相关研究未见报道。
虽然微卫星引物是根据特定微卫星位点的侧翼
序列设计的,具有较高的保守性,但侧翼序列位点的
突变或者碱基缺失,都会导致微卫星位点的扩增失
败。 门秋雷等 (2012 ) 发现欧洲苹果蠹蛾 Cydia
pomonella微卫星位点中部分位点在我国苹果蠹蛾
种群中不能稳定扩增;Ravel et al. (2002)研究埃及
伊蚊 Aedes aegypti时也出现微卫星位点扩增失败的
现象,因此开展相关研究必须进行大量的预备试验,
筛选出能够稳定扩增且多态性较好的微卫星位点。
Simon et al. (2001)从欧洲禾谷缢管蚜种群中筛选
得到 8 个多态性较高的微卫星位点,并被运用于欧
洲地区禾谷缢管蚜生殖方式和遗传结构的研究
(Delmotte et al. ,2002;Gilabert et al. ,2009),但是这
些微卫星位点是否适用于我国禾谷缢管蚜种群尚不
清楚。
关于我国禾谷缢管蚜种群的遗传多样性和遗传
结构的研究未见报道,而微卫星标记是蚜虫种群遗
传学研究中应用最广泛的分子标记之一。 本研究从
我国 8 个省(市)的 9 个不同地区采集禾谷缢管蚜,
分析已报道的欧洲种群中 8 个禾谷缢管蚜微卫星位
点在这 9 个种群试虫上的扩增稳定性及遗传多样
性,进而筛选出多态性较高且能够稳定扩增的微卫
星位点,以期为进一步研究我国不同地区禾谷缢管
蚜种群的遗传多样性、遗传分化、遗传结构及基因流
动奠定基础。
1 材料与方法
1􀆰 1 材料
供试昆虫:于 2013 年禾谷缢管蚜发生期,分别
在重庆北碚、陕西汉中、贵州贵阳、湖北武汉、安徽滁
州、吉林白城、青海西宁、山东菏泽和淄博等地小麦
田采集禾谷缢管蚜试虫。 每个地点采集试虫 50 头
以上,每株麦苗采集 1 头试虫,且同一地点各单头蚜
虫采集距离相隔 30 m以上,同一地点采集的蚜虫视
为一个种群,各种群代码(供试样本数)分别为 CBB
(33)、SXH(31)、GGY(34)、HWH(34)、ACZ(32)、
JBC(28)、QXN(28)、SHZ(31)和 SZB(31)。 将采集
的试虫分别放入装有小麦叶片的指型管内,每个指
型管放 1 头试虫带回实验室,置于温度 23 ± 1 ℃、相
对湿度(60 ± 5)% 、光周期 L ∶ D = 16 h ∶ 8 h 条件隔
离饲养于小偃 22 麦苗上备用。
试剂:基因组 DNA 提取试剂盒,北京全式金生
物技术有限公司;Taq MasterMix,北京康为世纪生物
科技有限公司。
仪器:伯乐 C1000 PCR 仪,美国伯乐(Bio⁃Rad)
生命医学产品(上海)有限公司;超微量核酸定量仪
Aosheng Nano⁃200, 杭 州 奥 盛 仪 器 有 限 公 司;
ABI3730XL型 DNA 自动分型仪,美国 Applied Bio⁃
892 植  物  保  护  学  报 42 卷
systems公司; DYY⁃12 型电泳仪,北京六一仪器厂。
1􀆰 2 方法
1􀆰 2􀆰 1 DNA的提取
将单头禾谷缢管蚜试虫置于离心管中并加入液
氮充分碾磨。 参照 DNA 提取试剂盒说明书分别提
取各头试虫的基因组 DNA,依次向离心管中加入试
剂盒各组分:裂解液( lysis buffer,LB)180 μL、蛋白
酶 K 20 μL、RNase A 20 μL、连接液(binding buffer
2,BB2) 200 μL、清洗液 ( clean buffer 2, CB2) 400
μL、漂洗液(wash buffer 2,WB2)400 μL 等,分别进
行离心,最后用去离子水定容至 30 μL。 取 3 ~ 4 μL
溶液提取基因组 DNA 并用 1􀆰 0%琼脂糖凝胶电泳
检测 DNA提取的质量,并用超微量核酸定量仪检测
基因组 DNA OD260值以及 OD260 / 280值,经检测合格的
样本基因组 DNA于 - 20 ℃下保存备用。
1􀆰 2􀆰 2 引物合成
根据目前已有的 8 对用于欧洲禾谷缢管蚜的微
卫星引物进行引物合成(Simon et al. ,2001;Wilson
et al. ,2004),各个微卫星位点的正向引物(F)和反
向引物(R)序列见表 1。 参照微卫星 PCR产物荧光
标记方法(Schuelke,2000),分别采用 3 条引物对各
样本中每个微卫星位点进行 PCR扩增,第 1 条引物
为各微卫星的正向引物 5′端加上M13 正向引物(5′⁃
TGTAAAACGACGGCCAGT⁃3′)而形成的上游引物;
第 2 条引物为不加任何修饰的各微卫星的反向引
物;第 3 条引物是 5′端带有羧基荧光素(6⁃carboxy⁃
fluorescine,FAM)标记的 M13 正向引物。 所有引物
均由生工生物工程(上海)股份有限公司合成。
表 1 禾谷缢管蚜 8 个微卫星位点的特性
Table 1 Characterizations of eight microsatellite loci in Rhopalosiphum padi
位点(GenBank登录号)
Locus (GenBank accession number)
引物序列(5′⁃3′)
Sequence of primers (5′⁃3′)
重复结构
Repeat motif
退火温度 (℃)
Anneal temperature
R1􀆰 35 (AF277461) F: CGTAGACCGCCGCGGG
R: GTCGTTTCTGGTCAGCGGCC
(TC) 13 60
R2􀆰 73 (AF277466) F: CGTAGACCGCCGCGGG
R: GTCGTTTCTGGTCAGCGGCC
(GT) 25GG(ATT) 4 60
R3􀆰 171 (AF277463) F: ACATCGTAAGACGTAAAACGAC
R: CAAAGCAATACCGCATAACG
(AT) 2T(AT) 9 …(AC) 11 51
R5􀆰 10 (AF277462) F: CCGACTAAGCTTAATATTGTTTG
R: CGGTTCGGAGAACATAAGAG
(CA) 3AC(AG) 6
(ATT) 5 …(GA) 15
52
R5􀆰 138 (AF349565) F: TATACACGCTCGCGCTTACG
R: CCGAGCACGAATTGTTCC
(CA) 35 51
R5􀆰 29b (AF277464) F: CATGAGTGTGTCCCTTTTAAC
R: GATGGACGAGGGGACAC
(AC) 13 58
R5􀆰 50 (AF349550) F: TGTTACGCGGAGTGTGTAGG
R: CCACAGAGCGTTGTCATC
(AC) 22 56
R6􀆰 3 (AF277465) F: CGAAATGTACCCACTATAAAC
R:CAAATTTAAATGTATAATCAATG
(AT) 7A(AT) 6(TA) 3 46
    F: 各微卫星位点的正向引物; R:各微卫星位点的反向引物。 F: The forward primer for microsatellite locus; R: the reverse primer for
microsatellite locus.
1􀆰 2􀆰 3 微卫星等位基因扩增及基因型检测
微卫星等位基因扩增的 PCR 反应体系为 25
μL,其中包括 Taq MasterMix 12􀆰 5 μL、10 μmol / L含
有 M13 序列的上游引物 0􀆰 5 μL、10 μmol / L 反向引
物和 M13 引物各 2􀆰 0 μL、20 ~ 50 ng / μL 模板 DNA
1􀆰 5 μL,ddH2O补足 25 μL。 PCR 扩增反应由两步
PCR构成,反应条件为:94 ℃ 预变性 2 min,94 ℃
变性 20 s,各引物特异性退火温度(表 1)20 s,72 ℃
延伸 20 s,共 30 个循环;然后 95 ℃ 变性 30 s,53 ℃
(M13 引物退火温度)退火 45 s,72 ℃延伸 45 s,共 8
个循环;最后 72 ℃补偿延伸 10 min。 PCR 产物经
1􀆰 0%琼脂糖凝胶电泳检测后,送生工生物工程(上
海)股份有限公司进行基因型自动扫描检测。
1􀆰 2􀆰 4 等位基因分析
利用 GeneMapper 4􀆰 0 软件(Applied Biosystems,
CA,USA)将荧光检测获取的峰型电泳图谱转换成
扩增产物的片段长度,并进行等位基因分析。 利用
FREENA软件计算各位点在各种群中的无效等位基
因频率 ( null allele frequency,Na ) ( Chapuis & Es⁃
toup,2007)。 利用 Excel Microsatellite Toolkit软件计
9923 期 段辛乐等: 我国禾谷缢管蚜微卫星位点扩增稳定性及遗传多样性
算各位点的多态信息含量(polymorphism information
content,PIC)(Park,2001)。 利用 Arlequin 3􀆰 5 软件
计算在各位点上的平均等位基因数(mean number of
allele,NA )、观察杂合度 ( observed heterozygosity,
Ho)、期望杂合度( expected heterozygosity,He )及近
交系数( inbreeding index,F is) (Excoffier & Lischer,
2010)。 利用 Genepop 4􀆰 0􀆰 1 软件对各种群在各位
点上的哈迪 -温伯格平衡(Hardy⁃Weinberg equilib⁃
rium)P值进行无偏估测(Rousset,2008)。
2 结果与分析
2􀆰 1 微卫星位点扩增稳定性
对所选取禾谷缢管蚜的 8 个微卫星位点进行
PCR扩增,其中 R2􀆰 73、R3􀆰 171、R5􀆰 10、R5􀆰 138 和
R5􀆰 50 这 5 个微卫星位点在各种群中均能稳定扩
增,R6􀆰 3 在青海西宁、吉林白城种群和贵州贵阳
的少数样本中无法扩增,但对数据分析的有效性
无影响;R1􀆰 35 在各种群中均不能扩增;R5􀆰 29b 仅
能在重庆北碚、陕西汉中、贵州贵阳、湖北武汉、安
徽滁州和山东菏泽、淄博种群的少数样本中
扩增。  
2􀆰 2 微卫星位点多态信息含量
位点 R2􀆰 73、R5􀆰 138 和 R5􀆰 50 在 9 个种群中的
多态信息含量 ( PIC) 分别为 0􀆰 5048 ~ 0􀆰 7967、
0􀆰 6220 ~ 0􀆰 8941 和 0􀆰 5497 ~ 0􀆰 8198,这 3 个微卫星
位点的 PIC均大于 0􀆰 5;位点 R3􀆰 171 在青海西宁种
群的 PIC 值为 0􀆰 2620,而在其它 8 个种群中的 PIC
值为 0􀆰 5572 ~ 0􀆰 8327;位点 R5􀆰 10 在湖北武汉、山
东淄博和山东菏泽的 PIC 值分别为 0􀆰 1197、0􀆰 3000
和 0􀆰 4674,其它种群的 PIC 值为 0􀆰 5231 ~ 0􀆰 7142;
位点 R6􀆰 3 在各种群中 PIC 值均较低,为 0􀆰 2106 ~
0􀆰 4338(表 2)。
表 2 六个稳定扩增的禾谷缢管蚜微卫星位点的多态信息含量
Table 2 The polymorphism information content (PIC) of six stably amplified Rhopalosiphum padi microsatellite loci
位点
Locus
地理种群 Geographic population
JBC QXN SXH CBB GGY HWH SZB SHZ ACZ
R2􀆰 73 0􀆰 7967 0􀆰 7562 0􀆰 7653 0􀆰 7026 0􀆰 5048 0􀆰 5759 0􀆰 5421 0􀆰 5936 0􀆰 6040
R3􀆰 171 0􀆰 7982 0􀆰 2620 0􀆰 8327 0􀆰 8222 0􀆰 5655 0􀆰 6431 0􀆰 7121 0􀆰 5572 0􀆰 6307
R5􀆰 10 0􀆰 7142 0􀆰 5812 0􀆰 6145 0􀆰 6383 0􀆰 6504 0􀆰 1197 0􀆰 3000 0􀆰 4674 0􀆰 5231
R5􀆰 138 0􀆰 8733 0􀆰 8295 0􀆰 8905 0􀆰 8941 0􀆰 7531 0􀆰 6565 0􀆰 6834 0􀆰 6897 0􀆰 6220
R5􀆰 50 0􀆰 7386 0􀆰 6781 0􀆰 7523 0􀆰 8198 0􀆰 5914 0􀆰 7336 0􀆰 7862 0􀆰 6911 0􀆰 5497
R6􀆰 3 0􀆰 2106 0􀆰 3711 0􀆰 4338 0􀆰 4103 0􀆰 3744 0􀆰 3769 0􀆰 2554 0􀆰 3645 0􀆰 2608
    JBC: 吉林白城种群; QXN: 青海西宁种群; SXH: 陕西汉中种群; CBB:重庆北碚种群; GGY: 贵州贵阳种群; HWH: 湖北武汉种
群; SZB:山东淄博种群; SHZ:山东菏泽种群; ACZ: 安徽滁州种群。 JLB: Population sampled in Baicheng of Jilin Province; QXN: popu⁃
lation sampled in Xining of Qinghai Province; SXH: population sampled in Hanzhong of Shaanxi Province; CBB: population sampled in Beibei
of Chongqing; GGY: population sampled in Guiyang of Guizhou Province; HWH: population sampled in Wuhan of Hubei Province; SZB: popu⁃
lation sampled in Zibo of Shandong Province; SHZ: population sampled in Heze of Shandong Province; ACZ: population sampled in Chuzhou of
Anhui Province.
2􀆰 3 微卫星位点的遗传多样性
6 个稳定扩增的微卫星位点无效等位基因频率
为 0􀆰 0044 ~ 0􀆰 2663,其中位点 R6􀆰 3 的无效等位基
因频率最高,为 0􀆰 2663,位点 R2􀆰 73 的无效等位基
因频率最低, 为 0􀆰 0044, 位点 R3􀆰 171、 R5􀆰 10、
R5􀆰 138 和 R5􀆰 50 的无效等位基因频率分别为
0􀆰 0622、0􀆰 0922、0􀆰 0667 和 0􀆰 0100。 各微卫星位点
的平均等位基因数为 2􀆰 9 ~ 9􀆰 3 个,其中位点
R5􀆰 138 的平均等位基因数最多,为 9􀆰 3 个,位点
R6􀆰 3 的平均等位基因数最少,为 2􀆰 9 个。 各位点观
测杂合度为 0􀆰 047 ~ 0􀆰 912,期望杂合度为 0􀆰 421 ~
0􀆰 815,其中位点 R5􀆰 50 观测杂合度最高,为 0􀆰 912,
位点 R6􀆰 3 最低,为 0􀆰 047;位点 R5􀆰 138 期望杂合度
最高,为 0􀆰 815,位点 R6􀆰 3 最低,为 0􀆰 421。 在稳定
扩增的 6 个位点中,位点 R6􀆰 3 平均等位基因数较少
(2􀆰 9),且观察杂合度很低(0􀆰 047),表现出明显的
杂合子不足(表 3)。
2􀆰 4 微卫星位点哈迪 -温伯格平衡检测
本研究中能稳定扩增的 6 个微卫星位点在陕西
汉中、贵州贵阳、湖北武汉、山东淄博、山东菏泽和安
徽滁州种群均偏离哈迪 -温伯格平衡(P < 0􀆰 05),
在吉林白城种群中只有位点 R5􀆰 10 偏离哈迪 -温
伯格平衡,在青海西宁种群只有位点 R6􀆰 3 偏离哈
迪 -温伯格平衡,在重庆北碚种群除了位点 R5􀆰 10
以外均偏离哈迪 -温伯格平衡(表 4)。
003 植  物  保  护  学  报 42 卷
表 3 禾谷缢管蚜 6 个微卫星位点的遗传多样性
Table 3 The genetic diversity of six microsatellite loci in Rhopalosiphum padi
位点
Locus
平均等位基因数
NA
观察杂合度
Ho
期望杂合度
He
近交系数
F is
无效等位基因频率
Na
R2􀆰 73 6􀆰 1 0􀆰 899 0􀆰 718 - 0􀆰 253 0􀆰 0044
R3􀆰 171 7􀆰 2 0􀆰 635 0􀆰 699 0􀆰 092 0􀆰 0622
R5􀆰 10 5􀆰 3 0􀆰 466 0􀆰 564 0􀆰 175 0􀆰 0922
R5􀆰 138 9􀆰 3 0􀆰 803 0􀆰 815 0􀆰 015 0􀆰 0667
R5􀆰 50 6􀆰 4 0􀆰 912 0􀆰 762 - 0􀆰 197 0􀆰 0100
R6􀆰 3 2􀆰 9 0􀆰 047 0􀆰 421 0􀆰 889 0􀆰 2663
    NA: Mean number of allele; Ho: observed heterozygosity; He: expected heterozygosity; F is: inbreeding index; Na: null allele frequency.
表 4 不同种群在不同位点上的哈迪 -温伯格平衡 P值
Table 4 P values of Hardy⁃Weinberg equilibrium in each Rhopalosiphum padi population at each microsatellite locus
位点
Locus
地理种群 Geographic population
JBC QXN SXH CBB GGY HWH SZB SHZ ACZ
R2􀆰 73 0􀆰 8408 0􀆰 1395 0􀆰 0022∗∗ 0􀆰 0013∗∗ 0􀆰 0000∗∗ 0􀆰 0003∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗
R3􀆰 171 0􀆰 1046 0􀆰 3059 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0001∗∗ 0􀆰 0024∗∗ 0􀆰 0000∗∗
R5􀆰 10 0􀆰 0017∗∗ 0􀆰 0919 0􀆰 0147∗ 0􀆰 2878 0􀆰 0000∗∗ 0􀆰 0222∗ 0􀆰 0227∗ 0􀆰 0051∗∗ 0􀆰 0000∗∗
R5􀆰 138 0􀆰 5185 0􀆰 9569 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0317∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗
R5􀆰 50 0􀆰 3997 0􀆰 1383 0􀆰 0113∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0490∗ 0􀆰 0071∗∗ 0􀆰 0010∗∗ 0􀆰 0000∗∗
R6􀆰 3 0􀆰 2727 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0000∗∗ 0􀆰 0003∗∗ 0􀆰 0000∗∗ 0􀆰 0001∗∗
    JBC: 吉林白城种群; QXN: 青海西宁种群; SXH:陕西汉中种群; CBB: 重庆北碚种群; GGY: 贵州贵阳种群; HWH: 湖北武汉种
群; SZB: 山东淄博种群; SHZ: 山东菏泽种群; ACZ: 安徽滁州种群。 ∗和∗∗分别表示在 P <0􀆰 05和 P <0􀆰 01水平上显著偏离哈迪 -温
伯格平衡。 JLB: Population sampled in Baicheng of Jilin Province; QXN: population sampled in Xining of Qinghai Province; SXH: population
sampled in Hanzhong of Shaanxi Province; CBB: population sampled in Beibei of Chongqing; GGY: population sampled in Guiyang of Guizhou
Province; HWH: population sampled in Wuhan of Hubei Province; SZB: population sampled in Zibo of Shandong Province; SHZ: population
sampled in Heze of Shandong Province; ACZ: population sampled in Chuzhou of Anhui Province. ∗ or ∗∗ indicates significant departure from
Hardy⁃Weinberg equilibrium at P < 0􀆰 05 or P < 0􀆰 01 level, respectively.
3 讨论
由于新微卫星位点的开发费用较高,且工作量
大,因此利用已报道相同种或者近缘物种的微卫星
引物进行相关研究是一种非常有效的途径,但研究
中发现,相同物种的微卫星引物在不同种群中的个
体中具有不同的扩增稳定性与遗传多样性(Ravel et
al. ,2002;门秋雷等,2012),这可能是因为相同物种
受到不同环境的影响使微卫星位点侧翼序列发生突
变,或者微卫星位点产生缺失。 本研究中的 8 个禾
谷缢管蚜微卫星位点在法国禾谷缢管蚜种群中均能
够稳定扩增,并且都表现出较高的多态性(Delmotte
et al. ,2002;Gilabert et al. ,2009),但在我国种群中
只有 R2􀆰 73、R3􀆰 171、R5􀆰 10、R5􀆰 138、R5􀆰 50 和 R6􀆰 3
这 6 个位点能够稳定扩增,而位点 R1􀆰 35 在我国各
种群均不能扩增,位点 R5􀆰 29b仅能在部分种群的少
数样本中扩增,但反应条件不稳定,这 2 个位点在我
国种群中不能扩增的原因有待进一步研究。
多态信息含量(PIC)是衡量微卫星多态性的重
要指标,一般认为,当 PIC > 0􀆰 5 时,该微卫星位点为
高度多态性位点,当 0􀆰 25 < PIC < 0􀆰 5 时为中度多态
性位点,PIC < 0􀆰 25 时为低度多态性位点(Botstein
et al. ,1980)。 微卫星位点的多态信息含量与该位
点的可用性及使用效率相关,多态信息含量越高,种
群中该微卫星位点的杂合子比例越高,提供的遗传
信息就越多。 本研究中除 R6􀆰 3 位点外, R2􀆰 73、
R3􀆰 171、R5􀆰 10、R5􀆰 138 和 R5􀆰 50 位点整体上都表
现出高度多态性(PIC > 0􀆰 5),说明这 5 个微卫星标
记在禾谷缢管蚜种群中的变异程度较高,能为禾谷
缢管蚜遗传多样性研究提供充分的信息,适于我国
不同地区禾谷缢管蚜种群的遗传变异研究。
微卫星标记中无效等位基因的存在对群体遗传
结果有显著影响(Paetkau & Strobeck,1995)。 由于
引物侧翼序列变异、大片段等位基因缺失等因素导
致无效等位基因的产生,使得该位点的目的条带缺
失,而被认为种群中纯合子数量过剩,从而降低观测
杂合度和群体总体遗传多样性,进而影响对种群遗
传结构分析的正确性(Brookfield,1996)。 本研究中
1033 期 段辛乐等: 我国禾谷缢管蚜微卫星位点扩增稳定性及遗传多样性
R2􀆰 73、R3􀆰 171、R5􀆰 10、R5􀆰 138 和 R5􀆰 50 的无效等
位基因频率在 0􀆰 0044 ~ 0􀆰 0922 之间,这与 Gliabert
et al. (2009)在欧洲禾谷缢管蚜种群中报道的无效
等位基因频率相似。 而 R6􀆰 3 的无效等位基因频率
为 0􀆰 2663,无效等位基因可能就是导致 R6􀆰 3 位点
观测杂合度低于其它位点的原因之一。
等位基因数和杂合度是微卫星位点遗传多样性
的重要指标,二者与位点的遗传多样性呈正相关
(Maudet et al. ,2002;闫路娜和张德兴,2004)。 本
研究中 R2􀆰 73、R3􀆰 171、R5􀆰 10、R5􀆰 138 和 R5􀆰 50 的
平均等位基因数为 6􀆰 86,表明这 5 对微卫星引物具
有较丰富的遗传信息以及均匀的基因频率分布,用
其研究我国禾谷缢管蚜的遗传多样性具有较高的可
靠性和有效性。 期望杂合度则能够反映出不同遗传
结构之间的变异程度,一般将杂合度高于 0􀆰 1 的位
点称为多态位点 ( Ott, 1999 ),本研究中 R2􀆰 73、
R3􀆰 171、R5􀆰 10、R5􀆰 138 和 R5􀆰 50 位点的平均期望
杂合度为 0􀆰 741,均具有较高的多态性和丰富的遗
传多样性。 而位点 R6􀆰 3 则表现出较低的观察杂合
度,说明该位点杂合子缺乏,同时该位点等位基因数
和多态信息含量较低,无效等位基因频率较高,说明
该位点在我国禾谷缢管蚜种群中具有较低的遗传多
样性,不适于我国禾谷缢管蚜种群遗传多样性研究。
哈迪 -温伯格平衡是指随机交配的足够大的种
群在理想状态下基因频率保持不变。 等位基因的杂
合子不足或杂合子过剩都能导致种群偏离哈迪 -温
伯格平衡。 在本研究中采自陕西汉中、重庆北碚、湖
北武汉、贵州贵阳、安徽滁州和山东菏泽、淄博等地
的种群显著偏离哈迪 -温伯格平衡,而吉林白城和
青海西宁种群基本处于哈迪 -温伯格平衡状态。 除
位点 R6􀆰 3 以外,本研究中其余 5 个稳定扩增的微卫
星位点的无效等位基因频率都比较低,因此无效等
位基因的存在不是导致种群偏离哈迪 -温伯格平衡
的主要原因。 已有研究表明,孤雌生殖的蚜虫因为
缺少基因交流,基因突变容易在种群中积累,因此具
有较高的等位基因杂合度,而同种蚜虫的有性生殖
种群等位基因杂合度相对较低(Ott,1999;Simon et
al. ,2002;Delmotte et al. ,2002)。 禾谷缢管蚜在我
国不同地区主要是孤雌生殖(未发表),因此分析认
为孤雌生殖导致的杂合子过剩可能是本研究中采自
重庆北碚、湖北武汉、贵州贵阳、安徽滁州、山东菏泽
和淄博等地理种群偏离哈迪 -温伯格平衡的主要原
因。 基于以上的遗传多样性指标分析,本试验表明
微卫星位点 R2􀆰 73、 R3􀆰 171、 R5􀆰 138、 R5􀆰 10 和
R5􀆰 50 能应用于我国禾谷缢管蚜种群遗传多样性和
遗传结构等相关研究。
参 考 文 献 (References)
Balloux F, Lugon⁃Moulin N. 2002. The estimation of population dif⁃
ferentiation with microsatellite markers. Molecular Ecology, 11
(2): 155 - 165
Botstein D, White RL, Skolnick M, Davis RW. 1980. Construction
of a genetic linkage map in man using restriction fragment length
polymorphisms. American Journal of Human Genetics, 32(3):
314 - 331
Brookfield JFY. 1996. A simple new method for estimating null al⁃
lele frequency from heterozygote deficiency. Molecular Ecology,
5(3): 453 - 455
Cao YZ, Yin J, Li KB, Zhang KC, Li XQ. 2006. The reasons of
resurgence and control strategies of wheat aphids. Plant Protec⁃
tion, 32(5): 72 - 75 (in Chinese) [曹雅忠, 尹娇, 李克斌,
张克诚, 李贤庆. 2006. 小麦蚜虫不断猖獗原因及控制对策
的探讨. 植物保护, 32(5): 72 - 75]
Chapuis MP, Estoup A. 2007. Microsatellite null alleles and estima⁃
tion of population differentiation. Molecular Biology and Evolu⁃
tion, 24(3): 621 - 631
Chen MH, Dorn S. 2010. Microsatellites reveal genetic differentia⁃
tion among populations in an insect species with high genetic
variability in dispersal, the codling moth, Cydia pomonella
(L. ) ( Lepidoptera: Tortricidae). Bulletin of Entomological
Research, 100(1): 75 - 85
Chen MH, Han ZJ. 2006. Cloning and sequence analysis of 2 differ⁃
ent acetylcholinesterase genes in Rhopalosiphum padi and Sito⁃
bion avenae. Genome, 49(3): 239 - 243
Delmotte F, Leterme N, Gauthier JP, Rispe C, Simon JC. 2002.
Genetic architecture of sexual and asexual populations of the a⁃
phid Rhopalosiphum padi based on allozyme and microsatellite
markers. Molecular Ecology, 11(4): 711 - 723
Excoffier L, Lischer HEL. 2010. Arlequin suite ver 3􀆰 5: a new se⁃
ries of programs to perform population genetics analyses under
Linux and Windows. Molecular Ecology Resources, 10 (3):
564 - 567
Fei XD, Li C, Zhang QW, Zhao ZW. 2011. The effects of wheat
planted adjacent to rape on wheat aphids population dynamic.
Journal of Plant Protection, 38 (4): 339 - 343 ( in Chinese)
[费晓东, 李川, 张青文, 赵章武. 2011. 油菜 -小麦邻作
模式对麦蚜种群动态的影响. 植物保护学报, 38 (4):
339 - 343]
Feng YY, Li W, Sun HY, Deng YY, Yu HS, Chen HG. 2013. Ge⁃
netic diversity of Gaeumannomyces graminis var. tritici popula⁃
tions in Huang⁃Huai winter wheat production region of China.
Journal of Plant Protection, 40(6): 495 - 501 ( in Chinese)
[冯彦霞, 李伟, 孙海燕, 邓渊钰, 于汉寿, 陈怀谷. 2013.
203 植  物  保  护  学  报 42 卷
黄淮麦区小麦全蚀病菌群体的遗传多样性分析. 植物保护
学报, 40(6): 495 - 501]
Fenton B, Margaritopoulos J, Malloch GL, Foster SP. 2010. Micro⁃
evolutionary change in relation to insecticide resistance in the
peach⁃potato aphid, Myzus persicae. Ecological Entomology, 35
(S1): 131 - 146
Gilabert A, Simon JC, Mieuzet L, Halkett F, Stoeckel S, Plante⁃
genest M, Dedryver CA. 2009. Climate and agricultural context
shape reproductive mode variation in an aphid crop pest. Mo⁃
lecular Ecology, 18(14): 3050 - 3061
Llewellyn KS, Loxdale HD, Harrington R, Brookes CP, Clark SJ,
Sunnucks P. 2003. Migration and genetic structure of the grain
aphid (Sitobion avenae) in Britain related to climate and clonal
fluctuation as revealed using microsatellites. Molecular Ecolo⁃
gy, 12(1): 21 - 34
Lu YH, Yang T, Gao XW. 2009. Establishment of baseline suscep⁃
tibility data to various insecticides for aphids Rhopalosiphum pa⁃
di (Linnaeus) and Sitobion avenae (Fabricius). Acta Entomo⁃
logica Sinica, 52(1): 52 - 58 (in Chinese) [鲁艳辉, 杨婷,
高希武. 2009. 禾谷缢管蚜和麦长管蚜玻璃管药膜法敏感
毒力基线的建立. 昆虫学报, 52(1): 52 - 58]
Maudet C, Miller C, Bassano B, Breitenmoser⁃Würsten, Gauthier
D, Obexer⁃Ruff G, Michallet J, Taberlet P, Luikart G. 2002.
Microsatellite DNA and recent statistical methods in wildlife
conservation management: applications in Alpine ibex [Capra
ibex ( ibex)]. Molecular Ecology, 11(3): 421 - 436
Men QL, Chen MH, Zhang YL, Feng JN. 2012. Amplifying stabili⁃
ty and genetic diversity of microsatellite loci in codling moth
Cydia pomonella (Lepidoptera: Tortricidae) populations in Chi⁃
na. Journal of Plant Protection, 39 (4): 341 - 346 ( in Chi⁃
nese) [门秋雷, 陈茂华, 张雅林, 冯纪年. 2012. 中国疫区
内苹果蠹蛾微卫星位点的扩增稳定性及遗传多样性. 植物
保护学报, 39(4): 341 - 346]
Ott J. 1999. Analysis of human genetic linkage. Baltimore: Johns
Hopkins University Press, pp. 24 - 36
Paetkau D, Strobeck C. 1995. The molecular basis and evolutionary
history of a microsatellite null allele in bears. Molecular Ecolo⁃
gy, 4(4): 519 - 520
Papura D, Simon JC, Halkett F, Delmotte F, Le Gallic JF, De⁃
dryver CA. 2003. Predominance of sexual reproduction in Ro⁃
manian populations of the aphid Sitobion avenae inferred from
phenotypic and genetic structure. Heredity, 90(5): 397 - 404
Park SDE. 2001. Trypanotolerance in West African cattle and the
population genetic effects of selection. Ph. D Thesis. Dublin:
University of Dublin
Ravel S, Hervé JP, Diarrassouba S, Kone A, Cuny G. 2002. Mic⁃
rosatellite markers for population genetic studies in Aedes aegypti
(Diptera: Culicidae) from Côte d′Ivoire: evidence for a micro⁃
geographic genetic differentiation of mosquitoes from Bouaké.
Acta Tropica, 82(1): 39 - 49
Rousset F. 2008. GENEPOP’007: a complete reimplementation of
the GENEPOP software for Windows and Linux. Molecular E⁃
cology Resources, 8(1): 103 - 106
Schuelke M. 2000. An economic method for the fluorescent labeling
of PCR fragments. Nature Biotechnology, 18(2): 233 - 234
Selkoe KA, Toonen RJ. 2006. Microsatellites for ecologists: a prac⁃
tical guide to using and evaluating microsatellite markers. Ecol⁃
ogy Letters, 9(5): 615 - 629
Simon JC, Blackman RL, Le Gallic JF. 1991. Local variability in
the life cycle of the bird cherry⁃oat aphid, Rhopalosiphum padi
(Homoptera: Aphididae) in western France. Bulletin of Ento⁃
mological Research, 81(3): 315 - 322
Simon JC, Leterme N, Delmotte F, Martin O, Estoup A. 2001. Iso⁃
lation and characterization of microsatellite loci in the aphid
species, Rhopalosiphum padi. Molecular Ecology Notes, 1(1 /
2): 4 - 5
Simon JC, Rispe C, Sunnucks P. 2002. Ecology and evolution of
sex in aphids. Trends in Ecology & Evolution, 17(1): 34 - 39
Vorwerk S, Forneck A. 2006. Reproductive mode of grape phylloxe⁃
ra ( Daktulosphaira vitifoliae, Homoptera: Phylloxeridae ) in
Europe: molecular evidence for predominantly asexual popula⁃
tions and a lack of gene flow between them. Genome, 49(6):
678 - 688
Wilson ACC, Massonnet B, Simon JC, Prunier⁃Leterme N, Dolatti
L, Llewellyn KS, Figueroa CC, Ramirez CC, Blackman RL,
Estoup A, et al. 2004. Cross⁃species amplification of microsat⁃
ellite loci in aphids: assessment and application. Molecular E⁃
cology Notes, 4(1): 104 - 109
Yan LN, Zhang DX. 2004. Effects of sample size on various genetic
diversity measures in population genetic study with microsatellite
DNA markers. Acta Zoologica Sinica, 50(2): 279 - 290 ( in
Chinese) [闫路娜, 张德兴. 2004. 种群微卫星 DNA分析中
样本量对各种遗传多样性度量指标的影响. 动物学报, 50
(2): 279 - 290]
Zhang GX, Zhong TS. 1983. Economic insect fauna of China fasc.
25 Homoptera, Aphidinea. Beijing: Science Press, pp. 356 -
362 ( in Chinese) [张广学, 钟铁森. 1983. 中国经济昆虫
志,第二十五册,同翅目蚜虫类(一). 北京: 科学出版社,
pp. 356 - 362]
Zhang XC, Zhou GH, Shi M, Fang JZ, Zhao ZP, Li SH, Dong QZ,
Wei K. 1985. The rules of long distance migration and virus
transmission of wheat aphids. Journal of Plant Protection, 12
(1): 9 - 16 (in Chinese) [张向才, 周广和, 史明, 方建中,
赵争平, 李淑华, 董庆周, 魏凯. 1985. 麦蚜远距离迁飞和
传毒规律的研究. 植物保护学报, 12(1): 9 - 16]
Zheng Y, Peng X, Liu GM, Pan HY, Dorn S, Chen MH. 2013.
High genetic diversity and structured populations of the oriental
fruit moth in its range of origin. PLoS ONE, 8(11): e78476
(责任编辑:高  峰)
3033 期 段辛乐等: 我国禾谷缢管蚜微卫星位点扩增稳定性及遗传多样性