论文部分内容阅读
二化螟Chilo suppressalis（Walker）属鳞翅目Lepidoptera,草螟蛾科Crambidae,其分布范围广、经济危害重,是水稻种植区的最主要害虫之一。由于其食性较杂,寄主较多,分布范围广泛,二化螟与众多寄主之间以及与复杂多变的环境之间的交互作用可能会带来其较强的遗传变异,从而增加对该种害虫防治工作的难度。因此,了解二化螟不同地理种群的遗传结构,分析各个种群间的亲缘关系,再现演化历史,探讨分化成因及预防机制就显得非常必要。本研究首次采用扩增二化螟4个微卫星位点（microsatellite loci）和4个线粒体DNA（mtDNA）基因片断等两种分子标记相结合的手段,运用系统发育分析方法、谱系生物地理学方法以及遗传变异参数分析等多种生物信息学的研究方法,对采自我国18个不同地区的二化螟种群进行系统而深入的比较研究,获得以下结果与结论:1)扩增了我国18个不同地区的二化螟种群中381个个体的4个微卫星位点（microsatellite loci）和4个线粒体DNA（mtDNA）基因片断。这18种群基本覆盖了二化螟在我国的分布区域。在4个微卫星位点上,共检测到39个等位基因。在4个线粒体基因中,共获得并向GenBank中注册了98个单倍型,其中包括:22个16S基因的单倍型,27个COⅠ基因的单倍型,24个COⅡ基因的单倍型,和25个ND1基因的单倍型。2)微卫星标记（microsatellite marker）和线粒体（mtDNA）基因序列数据均表明18个取样点的二化螟种群间出现了明显的遗传分化,只是基于前者获得的结果较后者的为小。以微卫星标记得出的种群间分化参数FSTS在-0.00011-0.12946之间,以16S基因序列数据求得的种群间分化参数FSTS在-0.00358-0.41001之间,COⅠ数据的FSTS在-0.01918-0.55016之间变化,COⅡ数据的FSTS在-0.02904-0.38614之间波动,ND1数据的FSTS在-0.03202-0.48485间变动。其中COⅠ和COⅡ基因具有很高的遗传变异度。平均核酸多样度分别为0.86%和0.69%,变异位点占序列长度的比例依次为13.9%和9.6%。3)基于微卫星标记数据构建的NJ树和依据线粒体基因序列数据生成的MP树以及ML树,以及4个基因的单倍型网状图分析均可得出如下结论:中国的二化螟C.suppressalis很可能分化为遗传上不同的、与地理位置相一致的3个分支,即华中分支、西南分支和华北东北分支。这3大分支在南北方向上被秦岭、伏牛山、长江、淮河等山脉和河流分隔开来;长江以南区域,则以罗霄山脉为界分为东西两支。而且,二化螟的遗传分化大致与其在3大地理区域内每年发生的代次相一致。在华中地区,二化螟一年发生3到4代,而在西南地区和华北东北地区,二化螟分别发生2到3代和1到2代。4)中性检验结果表明华中、西南和华北东北等3大地理区域间存在明显的遗传分化。当18个种群放在一起考虑时,4个基因的Tajima’s D均为显著性的负值（从-1.85到-1.47,P＜0.05）,Fu’s F值为显著性的负值或接近显著性的负值（-7.59--2.90,P=0.03-0.26）;当把华中、西南以及华北东北等3大地理区域的种群各自单独作为一组分析时,虽然Tajima’s D未达到显著性水平,但仍为负值（从-1.37到-0.42）,而Fu’s F值在不同的区域中出现负值甚至出现正值（-0.49-2.27）,表明3大地理区域间遗传分化明显。5)基于微卫星标记数据和mtDNA基因序列数据的遗传变异分析（AMOVA）结果均显示我国的二化螟种群遗传分化明显。种群分化参数FST为:微卫星标记数据计算结果是0.06004(P＜10-4),16S、COⅠ、COⅡ和ND1等4个mtDNA基因序列数据计算结果分别为0.19485(P＜10-4)、0.27607(P＜10-4)、0.22949(P＜10-4)和0.29285(P＜10-4)。华中、西南和华北东北等3大地理区域间种群分化参数FCT为:微卫星标记数据计算结果是0.03855(P＜10-4),16S、COⅠ、COⅡ和ND1等4个基因序列数据计算的FCT分别为0.10682(P＜10-4)、0.21533(P＜10-4)、0.16862(P＜10-4)和0.23845(P＜10-4)。再次证明了中国的二化螟种群已分化为华中、西南和华北东北等3个分支。各区域内种群分化参数FSC为:微卫星标记数据计算结果是0.02236(P＜10-4),16S、COⅠ、COⅡ和ND1等4个基因序列数据计算的FSC分别为0.09855(P＜10-4)、0.07741(P＜10-4)、0.07322(P＜10-4)和0.07143(P＜10-4)。两种分子标记手段获得了相同的结果,即区域间种群分化参数FCT是同一标记方法的区域内种群分化参数FSC的二倍或更多,表明华中、西南和华北东北等3大区域间基因交流有限。6)二化螟两两种群间的基因交流水平较低（M＜100,81.7%Nem＜1.0）,即使在地理距离最近的种群间也是如此。这可能是由于地理隔离以及该害虫的非迁飞性特性引起的。在本研究涉及的整个取样范围内,没有发现与经度或纬度相一致的单向性基因流动。不过,当分析范围缩小至华中、西南或华北东北这些较小的区域时,基因流动的趋向变得明显起来。在华中地区,基因流倾向于从宁波流向西部,如衢州、南昌等地。在西南地区,基因流可能是从柳州向西流向贵阳、雅安等地。在华北东北地区,基因则可能流自武汉或驻马店,向北到沈阳、长春等地。在3个小范围内基因交流比较明显,而在全国整个取样点的范围内缺乏基因交流进一步确认了地理隔离作用,而且暗示着中国的二化螟可能起源于3个不同的避难所（refugia）,并经历了平行进化从而形成了目前的遗传结构。本研究所得的二化螟有效种群（effective populaion size）较小、距离最近的种群间基因交流有限以及很强的距离隔离关系（isolation-by-distance relationship）等结果表明二化螟处于遗传平衡状态（genetic equilibrium）,3个地区之间的遗传分化很可能是迁徙-漂变平衡的结果（migration-drift equilibrium）。7)二化螟4个mtDNA基因的错配分布（mismatch distribution）均为单峰（unimodal）。按照2.3%/1百万年的分子变异率,依据上述4个基因序列数据计算的合并时间（相当于群体扩张的启动时间）依次大约是65000或128800年,100000或282000年,100000或282000年,和100000或133000年。总的来说,在这4个基因中,每个基因在3大区域中的绝大多数τ值基本相近,各区域间τ值的95%置信值（95%C.I.）相互重叠,表明各区域的种群扩张时间大体相同。在华中地区,COⅠ,COⅡ,和ND1基因数据得出的合并时间早于60000年。在华北东北地区,16S,COⅠ,和ND1基因的分化时间大于120000年。在西南地区,16S,COⅠ,COⅡ,和ND1基因分化早于140000年前。只有华北东北地区的COⅡ基因和华中地区的16S基因例外。两者计算结果为种群分化时间晚于10000年前。考虑到昆虫分子钟的不精确性和缺乏必要的昆虫化石资料,3大分支分隔于3个避难所的时间大约开始于60000年以前,这个时间远早于水稻的人工栽培（domestication）时间。这可能表示二化螟的种群分化很可能是气候因素和/或地质事件（如,最后一次冰盛期）单独或共同作用的结果,二化螟在野生稻的人工栽培过程中经历了寄主转换。伴随着野生稻的人工栽培过程,栽培稻为蛀食于其植株体内的蛀虫提供了新的生境,进而强化了业已形成的遗传分化。