鐘?？? 杜素潔, 潘立婷, 王玉生, 王福蓮, 劉萬學(xué),*

(1. 長(zhǎng)江大學(xué)農(nóng)學(xué)院, 長(zhǎng)江大學(xué)昆蟲研究所, 湖北荊州 434025; 2. 中國(guó)農(nóng)業(yè)科學(xué)院植物保護(hù)研究所, 植物病蟲害生物學(xué)國(guó)家重點(diǎn)實(shí)驗(yàn)室, 北京 100193)

蔥斑潛蠅Liriomyzachinensis隸屬雙翅目(Diptera)潛蠅科(Agromyzidae)斑潛蠅屬Liriomyza，是一種非常重要的潛葉蠅類害蟲，主要為害大蔥Alliumfistulosum、韭菜A.tuberosum、洋蔥A.cepa等百合科(Liliaceae)蔥屬Allium植物，以大蔥受害最重(Spencer, 1973; 康樂, 1996)。其幼蟲、成蟲均可為害，主要以幼蟲潛食葉肉為主，葉肉被取食后，僅剩兩層表皮，呈現(xiàn)白色線狀條斑，形成潛道，嚴(yán)重影響植物正常的光合作用；并且雌成蟲的產(chǎn)卵器可刺破葉片，由此造成的機(jī)械損傷為病原物侵入提供了重要途徑(潘秀美和姜官恒, 2001; 魏國(guó)先等, 2001)。蔥斑潛蠅廣泛分布于東南亞和太平洋地區(qū)，在中國(guó)、日本、韓國(guó)、印度尼西亞、新加坡、朝鮮、越南等國(guó)家危害嚴(yán)重(Tran and Takagi, 2005; Tranetal., 2007)，近年來該蟲也在歐洲部分地區(qū)如法國(guó)、德國(guó)等地暴發(fā)成災(zāi)(Martynovetal., 2016)。該蟲在我國(guó)大部分地區(qū)均具有較高的適應(yīng)性(梁超峰等, 2018)，主要分布于黑龍江、內(nèi)蒙古、遼寧、沈陽、山東、河北、河南、山西、福建、臺(tái)灣等省區(qū)(王永衛(wèi)等, 1992; 王進(jìn)忠等, 2000; Shiao, 2004)。蔥斑潛蠅的溫度適應(yīng)范圍廣，單雌產(chǎn)卵量大，在我國(guó)各地一年發(fā)生4～7代不等，如東北4～5代，山東6～7代，且世代重疊發(fā)生(王立霞等, 2000; 潘秀美和姜官恒, 2001)，有春季和秋季兩個(gè)發(fā)生高峰期，高峰期易在大蔥種植區(qū)暴發(fā)成災(zāi)(潘秀美和姜官恒, 2000; 王立霞等, 2000; Tranetal., 2007)。

大蔥作為蔥斑潛蠅的主要寄主植物，其大規(guī)?；娣e種植一定程度上有利于蔥斑潛蠅的嚴(yán)重發(fā)生。據(jù)2014年農(nóng)業(yè)農(nóng)村部種植業(yè)管理司的統(tǒng)計(jì)結(jié)果，我國(guó)大蔥種植面積已逾57.3萬hm2，主產(chǎn)區(qū)位于秦嶺淮河以北的廣大地區(qū)，涵蓋華北平原、東北平原和黃土高原三大主產(chǎn)區(qū)，其中華北平原的山東、河南和河北三省為我國(guó)大蔥的第一大產(chǎn)區(qū)，總栽培面積占全國(guó)的36%，總產(chǎn)量占全國(guó)的51%(彭帥等, 2017)。山東省作為我國(guó)蔥蒜類蔬菜種植面積最大的省份，2001年蔥蒜類蔬菜種植面積約26.7萬hm2，其80%的面積遭受蔥斑潛蠅的危害，被害株率常達(dá)40%以上(潘秀美和姜官恒, 2001)。傅蘇友和蘇保樂(2011)在調(diào)查山東濰坊市的蔥斑潛蠅發(fā)生情況時(shí)，發(fā)現(xiàn)其在大蔥上為害嚴(yán)重，一株蔥上幼蟲平均為7頭，被害株高達(dá)70%～80%，嚴(yán)重的達(dá)到100%。

目前有關(guān)蔥斑潛蠅的研究多局限于生物學(xué)特性、發(fā)生規(guī)律等基礎(chǔ)性研究，而有關(guān)蔥斑潛蠅的種群遺傳多樣性、種群擴(kuò)張等方面研究還暫未見報(bào)道。蔥斑潛蠅作為一種本地害蟲，相對(duì)于斑潛蠅屬的其他入侵種，如南美斑潛蠅L.huidobrensis、美洲斑潛蠅L.sativae、三葉草斑潛蠅L.trifolii等，其長(zhǎng)期累積的遺傳變異是否與其嚴(yán)重危害存在聯(lián)系，還尚不清楚。分子標(biāo)記技術(shù)的快速發(fā)展，很大程度加快了國(guó)內(nèi)外學(xué)者對(duì)斑潛蠅屬害蟲遺傳分化的研究。Morgan等(2000)首先對(duì)美國(guó)加利福尼亞中部和南部的南美斑潛蠅種群間進(jìn)行RAPD-PCR檢測(cè)，發(fā)現(xiàn)南美斑潛蠅種群間的遺傳變異主要來源于地理分布。王莉萍等(2007)利用rDNA-ITS1標(biāo)記對(duì)我國(guó)8個(gè)不同地理的美洲斑潛蠅種群進(jìn)行分析，得出各地理種群之間的遺傳分化很小，且遺傳分化很有可能與地理分布有關(guān)。本研究選用線粒體COI基因作為分子標(biāo)記，以我國(guó)8省12個(gè)不同地區(qū)的蔥斑潛蠅種群為研究對(duì)象，對(duì)不同地理種群的遺傳多樣性、遺傳分化程度、基因交流等進(jìn)行了分析，以揭示其種群系統(tǒng)發(fā)育地理格局與演變機(jī)制。研究結(jié)果對(duì)蔥斑潛蠅的種群遺傳動(dòng)態(tài)分析、種群擴(kuò)張機(jī)制探討具有一定的指導(dǎo)意義，同時(shí)也為制定蔥斑潛蠅的防治策略提供理論依據(jù)。

1 材料與方法

1.1 供試蟲源

蔥斑潛蠅標(biāo)本于2016-2018年采集自中國(guó)8省12個(gè)地點(diǎn)(表1)，采用五點(diǎn)取樣法，采集露天種植的優(yōu)勢(shì)寄主植物大蔥A.fistulosum和小蔥A.ascalonicum，采集蟲態(tài)為幼蟲或蛹，帶回實(shí)驗(yàn)室飼養(yǎng)。 待蟲源羽化后，浸于無水乙醇中，-40℃保存?zhèn)溆谩?/p>

表1 蔥斑潛蠅不同地理種群的采集信息Table 1 Collecting data of different geographical populations of Liriomyza chinensis

1.2 基因組DNA提取

參照De Barro和Driver(1997)的快速研磨法，并稍作改動(dòng)。具體步驟如下：從標(biāo)本管中取單頭蔥斑潛蠅雌成蟲，置于干凈濾紙上，待無水乙醇揮發(fā)后，吸取10 μL裂解液(0.05 mol/L Tris-HCl, pH 8.4, 50 mmol/L KCl, 0.45% Tween-20, 0.45% NP 40, 0.2% Gelatin, 80 μg/mL蛋白酶K)滴于封口膜上，用滅菌過的0.2 mL PCR管管底充分研磨蟲體，將提取液吸入0.2 mL PCR管中。吸取20 μL裂解液清洗PCR管底和封口膜2次，最后將清洗液吸入0.2 mL PCR管中。將提取液進(jìn)行金屬浴，程序: 65℃下30 min, 25℃ 下2 min, 96℃ 下10 min，滅活蛋白酶K，最終獲得50 μL的蔥斑潛蠅DNA，保存在-20℃冰箱中備用。

1.3 mtCOI基因序列的擴(kuò)增

采用引物組合C1-J-2183 (5′-CAACATTTATTTT GATTTTTTGG-3′)和TL2-N-3014 (5′-TCCATTGCAC TAATCTGCCATATTA-3′)(Simonetal., 1994)擴(kuò)增蔥斑潛蠅mtCOI基因的目的片段。引物由上海生工生物工程股份有限公司合成。

PCR反應(yīng)體系(25 μL): 10×buffer(含Mg2+)2.5 μL, dNTPs(2.5 mmol/L)0.4 μL, 正反引物(10 pmol/L)各0.4 μL,Taq酶(2.5 U/μL) 0.2 μL, 模板DNA 1 μL, ddH2O 20.1 μL。PCR反應(yīng)程序: 95℃預(yù)變性4 min; 95℃變性1 min, 52.5℃退火50 s, 72℃延伸90 s, 35個(gè)循環(huán)；繼續(xù)72℃延伸10 min。取20 μL擴(kuò)增產(chǎn)物在含有GoldViewⅡ的1%瓊脂糖凝膠上以110 V電泳分離30 min，然后以凝膠成像系統(tǒng)分析電泳結(jié)果并拍照保存。

1.4 PCR產(chǎn)物的純化及克隆測(cè)序

由于部分蟲源直接測(cè)序時(shí)出現(xiàn)重疊峰，為了確保測(cè)序的準(zhǔn)確性，對(duì)存在重疊峰的樣品進(jìn)行了克隆測(cè)序。在紫外燈下切下含有目的片段的明亮條帶，使用AxyPrep DNA凝膠回收試劑盒將PCR產(chǎn)物回收純化。將純化產(chǎn)物與克隆質(zhì)粒pEASY-T3 Cloning Vector連接，轉(zhuǎn)化到Trans1-T1感受態(tài)細(xì)胞內(nèi)，篩選陽性單克隆進(jìn)行菌落PCR檢測(cè)，之后送到上海生工生物工程股份有限公司進(jìn)行雙向測(cè)序。

1.5 序列整理及數(shù)據(jù)分析

首先將正反向引物測(cè)得的序列導(dǎo)入MEGA 7.0軟件(Kumaretal., 2016)中，結(jié)合峰圖識(shí)別軟件Chromas 2.6校正序列，再進(jìn)行多序列比對(duì)。利用NCBI數(shù)據(jù)庫(kù)對(duì)獲得的序列進(jìn)行BLAST比對(duì)分析，以確定擴(kuò)增的基因序列來自線粒體COI基因。

利用MEGA 7.0軟件(Kumaretal., 2016)對(duì)蔥斑潛蠅mtCOI基因序列進(jìn)行變異分析，計(jì)算變異位點(diǎn)數(shù)、簡(jiǎn)約信息位點(diǎn)數(shù)、堿基組成等信息。同時(shí)，基于Kimura-2-Parameter(K2P)模型，以美洲斑潛蠅L.sativae、南美斑潛蠅L.huidobrensis、三葉草斑潛蠅L.trifolii、番茄斑潛蠅L.bryoniae和豌豆彩潛蠅Chromatomyiahorticola的相應(yīng)序列為外群，采用非加權(quán)組平均(UPGMA)法構(gòu)建單倍型系統(tǒng)發(fā)育樹，系統(tǒng)樹各分支的置信度(bootstrap)均進(jìn)行1 000次重復(fù)檢驗(yàn)。利用DnaSP 6.1軟件(Rozasetal., 2017)對(duì)不同地理種群的蔥斑潛蠅進(jìn)行單倍型多樣性(Hd)、核苷酸多樣性(Pi)和Tajima’sD中性檢驗(yàn)等分析。應(yīng)用Network 4.6基于Median-joining算法構(gòu)建單倍型間的網(wǎng)絡(luò)關(guān)系圖。

應(yīng)用Arlequin 3.5軟件(Excoffier and Lischer, 2010)計(jì)算種群間的遺傳分化系數(shù)(FST)(Hudsonetal., 1992)、遺傳距離[FST/(1-FST)]、基因流(Nm)等參數(shù)；并進(jìn)行分子方差分析(AMOVA)，計(jì)算遺傳變異在種群間和種群內(nèi)的方差貢獻(xiàn)率。

根據(jù)采集地經(jīng)緯度計(jì)算出成對(duì)種群間的地理距離，同時(shí)結(jié)合遺傳分化系數(shù)(FST)，用TFPGA軟件(Miller, 1997)中的Mantel算法檢測(cè)不同種群間的遺傳距離和地理距離的相關(guān)性，進(jìn)行10 000次置換，并做顯著性檢驗(yàn)。

2 結(jié)果

2.1 蔥斑潛蠅線粒體COI基因的堿基組成及單倍型分布

將中國(guó)蔥斑潛蠅12個(gè)地理種群共253條mtCOI基因序列比對(duì)、剪切后得到759 bp的片段，分析表明，無堿基的插入或缺失，具明顯的A/T堿基偏倚性，A+T平均含量為70.3%；序列共有20個(gè)變異位點(diǎn)，占?jí)A基總數(shù)的2.64%，包括13個(gè)單變異多態(tài)性位點(diǎn)和7個(gè)簡(jiǎn)約信息位點(diǎn)。

單倍型的數(shù)量及其分布和頻次的結(jié)果顯示(表2)，供試種群中共發(fā)現(xiàn)13種單倍型Hap1-Hap13 (GenBank登錄號(hào): MN700917-MN700929)，包括7個(gè)共享單倍型(Hap1, Hap3, Hap6, Hap8-Hap11)，其余6個(gè)單倍型為相應(yīng)地理種群獨(dú)享單倍型，其中Hap2屬于北京海淀種群(BJHD), Hap4和Hap5屬于甘肅白銀種群(GSBY), Hap7屬于海南陵水種群(HNLS), Hap12和Hap13屬于吉林通化種群(JLTH)。Hap1出現(xiàn)頻率最高，占總樣本的81.82%，且分布最廣，為所有地理種群共有，其余12種單倍型出現(xiàn)頻率均較低(≤3.95%)。各地理種群中，吉林通化種群(JLTH)的單倍型種類最多(5種)，其次為河北石家莊種群(HBSJZ)、河南許昌種群(HNXC)和山東濟(jì)南種群(SDJN)有4種，而河北保定種群(HBBD)與內(nèi)蒙烏海種群(NMWH)僅有1種單倍型。

2.2 蔥斑潛蠅種群的遺傳多樣性及中性檢驗(yàn)

蔥斑潛蠅各種群mtCOI基因單倍型多樣性(Hd)、核苷酸多樣性(Pi)、核苷酸平均差異數(shù)(K)、Tajima’sD中性檢測(cè)分析結(jié)果顯示(表2)，總?cè)后w單倍型多樣性(Hd)為0.327，核苷酸多樣性(Pi)為0.00159，核苷酸平均差異數(shù)(K)為1.21011，中性檢驗(yàn)結(jié)果不顯著(Tajima’sD=-1.64799,P>0.05)，說明供試蔥斑潛蠅總?cè)后w在過去較近的一段歷史時(shí)期內(nèi)未經(jīng)歷明顯的群體擴(kuò)張和持續(xù)增長(zhǎng)模式，群體大小保持相對(duì)穩(wěn)定狀態(tài)。各地理種群?jiǎn)伪缎投鄻有?Hd)在0～0.763間，核苷酸多樣性(Pi)在0～0.00716間，核苷酸平均差異數(shù)在0.09524～5.43684間，吉林通化種群(JLTH)的單倍型多樣性最高(Hd=0.763)，其次是河南許昌種群(HNXC,Hd=0.562)和河北石家莊種群(HBSJZ,Hd=0.489)，而河北保定種群(HBBD,Hd=0)和內(nèi)蒙烏海種群(NMWH,Hd=0)遺傳多樣性最低。各種群間的Tajima’sD中性檢測(cè)分析結(jié)果表明，甘肅白銀種群(GSBY)、河北石家莊種群(HBSJZ)和河南許昌種群(HNXC)表現(xiàn)顯著，山東濟(jì)南種群(SDJN)、山東臨沂種群(SDLY)和山東煙臺(tái)種群(SDYT)表現(xiàn)極顯著，說明這6個(gè)地區(qū)的蔥斑潛蠅種群在較近的歷史時(shí)期經(jīng)歷了種群擴(kuò)張，而其他地區(qū)的群體表現(xiàn)均不顯著。

表2 蔥斑潛蠅12個(gè)不同地理種群基于mtCOI基因序列的遺傳多樣性分析及中性檢驗(yàn)Table 2 Genetic diversity and neutrality test of Liriomyza chinensis among 12 geographical populations based on mtCOI gene sequence

2.3 蔥斑潛蠅mtCOI基因單倍型遺傳距離及系統(tǒng)進(jìn)化分析

遺傳分化的結(jié)果顯示(表3)，各單倍型間的K2P遺傳距離范圍為0.001～0.020，單倍型Hap11與其他單倍型遺傳距離最大，遺傳距離達(dá)到了0.012～0.020；其次是單倍型Hap2, Hap3, Hap8 和Hap13，這4個(gè)單倍型與其他單倍型之間的遺傳距離在0.003～0.019之間；其余7個(gè)單倍型之間遺傳距離均很接近(≤0.003)。

表3 基于蔥斑潛蠅地理種群mtCOI基因序列不同單倍型的遺傳距離Table 3 Genetic distances of different haplotypes of Liriomyza chinensis geographical populations based on mtCOI gene sequence

從蔥斑潛蠅mtCOI基因單倍型系統(tǒng)發(fā)生樹可見(圖1)，蔥斑潛蠅13種單倍型中，Hap11單獨(dú)形成一支系，其他12種單倍型聚為另一支系。這12種單倍型聚成的大支又可分為兩小支系，第一小支包含Hap2, Hap3, Hap8和Hap13共4個(gè)單倍型，第二小支系包含Hap1, Hap4-Hap7, Hap9-10及Hap12共8種單倍型。

圖1 UPGMA法構(gòu)建的基于mtCOI基因序列的蔥斑潛蠅地理種群 及近緣種的單倍型系統(tǒng)發(fā)育樹Fig. 1 Phylogenetic tree of haplotypes of Liriomyza chinensis geographical populations and related species constructed by UPGMA method based on mtCOI gene sequence Hap1-Hap13: 單倍型Haplotypes.

為了進(jìn)一步直觀體現(xiàn)出蔥斑潛蠅mtCOI基因單倍型之間的遺傳分化及在種群中的分布頻率，構(gòu)建13種單倍型的Network網(wǎng)絡(luò)關(guān)系圖(圖2)。結(jié)果顯示，以單倍型Hap1為中心，其他單倍型均圍繞這一單倍型呈星形輻射分布狀，據(jù)此推測(cè)分布廣泛且分布頻次最高的Hap1為祖先單倍型。

圖2 蔥斑潛蠅地理種群mtCOI基因單倍型中介網(wǎng)絡(luò)關(guān)系圖Fig. 2 Median-joining haplotype network of Liriomyza chinensis geographical populations based on mtCOI gene sequence 每個(gè)圓形表示一種單倍型，圓形面積反映出其單倍型在總?cè)后w中的出現(xiàn)頻率，H_1-H_13分別與單倍型Hap1-Hap13相對(duì)應(yīng)，mv1和mv2表示可能的突變位點(diǎn)。Each circle indicates one haplotype, and the circle size represents the frequency of haplotypes appearing in all populations tested. H_1-H_13 correspond to haplotype Hap1-Hap13, respectively, and mv1 and mv2 indicate possible mutation sites.

2.4 不同地理種群間的遺傳分化、基因交流和分子方差分析

蔥斑潛蠅總?cè)后w的基因流(Nm)為3.33629， 遺傳分化系數(shù)(FST)為0.06971, 符合1≤Nm<4, 0.05

表4 蔥斑潛蠅成對(duì)種群間遺傳分化系數(shù)(FST, 下三角)與種群間的基因流(Nm, 上三角)Table 4 Pairwise genetic differentiation coefficient (FST, below the diagonal) and gene flow (Nm, above the diagonal) between geographical populations of Liriomyza chinensis

對(duì)供試蔥斑潛蠅種群進(jìn)行的分子方差分析(AMOVA)結(jié)果(表5)表明，蔥斑潛蠅種群內(nèi)遺傳變異顯著(P<0.05)，占總變異的93.23%，說明遺傳變異主要來自于種群內(nèi)部；種群間遺傳變異水平低，僅占總變異的6.97%；不同寄主間(FCT=-0.00485<0.05)表現(xiàn)出較低的遺傳分化水平，各寄主植物間的遺傳變異不明顯(P>0.05)；同一寄主不同地理種群間(FSC=0.07216)表現(xiàn)出中等程度的遺傳分化，遺傳變異顯著(P<0.05)。

表5 基于mtCOI基因的蔥斑潛蠅種群間分子方差分析Table 5 Analysis of molecular variance (AMOVA) within and among populations of Liriomyza chinensis based on mtCOI gene

對(duì)所有種群進(jìn)行Mantel檢驗(yàn)，發(fā)現(xiàn)遺傳距離與地理距離之間不存在顯著的相關(guān)性(r=0.3820,P>0.05)(圖3)。

圖3 蔥斑潛蠅12個(gè)地理種群間遺傳距離與地理距離的相關(guān)性Fig. 3 Correlation between the genetic distance and the geographical distance for 12 geographical populations of Liriomyza chinensis FST/(1-FST)表示遺傳距離。The genetic distance is expressed as FST/(1-FST).

Va,Vb,Vc: 分別表示方差組分的數(shù)量Number of variance components;FCT: 寄主種群間變異Variation among host populations;FSC: 寄主種群內(nèi)變異Variation within host populations;FST: 種群內(nèi)變異Variation within a population.

3 討論

線粒體DNA(mtDNA)是常用的分子遺傳標(biāo)記，其中的線粒體細(xì)胞色素氧化酶亞基I(COI)基因在昆蟲種類鑒定、隱存種的發(fā)現(xiàn)、種群遺傳多樣性、物種的系統(tǒng)進(jìn)化等研究方面均得到廣泛應(yīng)用(DeSalleetal., 1987; 成新躍等, 2000; Hebertetal., 2003; 何恒果, 2008; Norgateetal., 2009; Yeetal., 2018)?；趍tCOI基因，本研究對(duì)我國(guó)8個(gè)省/直轄市的12個(gè)不同地區(qū)的蔥斑潛蠅種群的遺傳多樣性、遺傳分化程度、基因交流等進(jìn)行了深入分析。

本研究結(jié)果表明，蔥斑潛蠅mtDNA COI基因序列片段的堿基組成表現(xiàn)較強(qiáng)的A/T 偏倚性，與典型的昆蟲線粒體堿基組成特點(diǎn)(Simonetal., 1994; Wangetal., 2014; Chenetal., 2018)相符。在12個(gè)地理種群中共定義了13個(gè)單倍型，包括7個(gè)共享單倍型和6個(gè)獨(dú)享單倍型，說明不同地理區(qū)域的蔥斑潛蠅在保持著穩(wěn)定遺傳變異的同時(shí)，也產(chǎn)生了適應(yīng)特定地區(qū)環(huán)境的遺傳變異。在7個(gè)共享單倍型中，單倍型Hap1在所有地理種群中均有分布且出現(xiàn)頻率最高，同時(shí)結(jié)合單倍型網(wǎng)絡(luò)結(jié)構(gòu)圖來看，可認(rèn)為單倍型Hap1為蔥斑潛蠅群體中適應(yīng)性較強(qiáng)并廣泛分布的原始單倍型。

本研究中蔥斑潛蠅總?cè)后w的單倍型多樣性較低(表2)。有研究報(bào)道，地理距離、寄主植物種類以及生境等是影響昆蟲的種群遺傳多樣性和遺傳結(jié)構(gòu)的重要因素(Morganetal., 2000; Scheffer and Lewis, 2006; Tudaetal., 2014)。而本研究Mantel檢驗(yàn)結(jié)果顯示，地理距離與蔥斑潛蠅的遺傳距離間沒有相關(guān)性，即地理距離并不影響蔥斑潛蠅的遺傳分化(圖3)。從寄主植物的角度來講，蔥斑潛蠅是寡食性昆蟲，僅取食百合科蔥屬，且主要寄主植物為大蔥、小蔥，在其他植物上少見發(fā)生，種類單一的寄主植物可能在一定程度上限制了蔥斑潛蠅種群的遺傳分化。在12個(gè)地理種群中，僅吉林通化種群的單倍型多樣性表現(xiàn)較高(Hd=0.763)(表2)。結(jié)合采樣地特征分析其原因，吉林通化的蔥斑潛蠅為農(nóng)家自留地采集，與商業(yè)化大規(guī)模種植的蔥相比，化學(xué)農(nóng)藥的使用量較少，較低的藥劑選擇壓力可能增加了蔥斑潛蠅的遺傳變異?？偟膩碚f，作者認(rèn)為蔥斑潛蠅在吉林通化可能受到人工選擇壓力較小以及其他原因，造成了其較高的遺傳多樣性。

遺傳分化系數(shù)(FST)和基因流(Nm)是研究種群遺傳分化的重要參數(shù)，可以解釋影響種群發(fā)生遺傳分化的因素(Rousset, 1997)。Tajima’sD檢驗(yàn)可揭示物種在一定歷史時(shí)期內(nèi)的種群動(dòng)態(tài)等情況，通過檢驗(yàn)值來判斷種群是否存在擴(kuò)張及瓶頸效應(yīng)(Tajima, 1989)。在本研究中，蔥斑潛蠅總?cè)后w的遺傳分化系數(shù)(FST)為0.06971(0.05

本研究?jī)H以mtDNA COI基因長(zhǎng)為759 bp的單個(gè)片段作為分子遺傳標(biāo)記，為研究蔥斑潛蠅的系統(tǒng)進(jìn)化等提供了基礎(chǔ)性遺傳資料。本研究發(fā)現(xiàn)蔥斑潛蠅的遺傳多樣性水平低，后續(xù)還需要增加不同生境類型、寄主植物種類、蟲源采集量等數(shù)據(jù)，才能更深層次地系統(tǒng)闡析不同地理區(qū)域蔥斑潛蠅的遺傳分化。