基于CO I基因全長(zhǎng)序列的假眼小綠葉蟬地理種群遺傳分化研究

周寧寧, 王夢(mèng)馨, 崔 林, 潘 鋮, 張新亭, 韓寶瑜

(中國(guó)計(jì)量學(xué)院浙江省生物計(jì)量及檢驗(yàn)檢疫技術(shù)重點(diǎn)實(shí)驗(yàn)室， 杭州 310018)

假眼小綠葉蟬Empoascavitis(G?the) 隸屬于半翅目(Hemiptera) 、葉蟬科(Cicadellidae)、小綠葉蟬屬(Empoasca)，廣泛分布于我國(guó)大陸茶區(qū)，繁殖快，世代重疊嚴(yán)重，抗藥性強(qiáng)，防治難度大，是我國(guó)茶樹(shù)上為害最重的害蟲(chóng)。常年假眼小綠葉蟬可致長(zhǎng)江中下游茶區(qū)夏、秋茶減產(chǎn)10%—15%，嚴(yán)重發(fā)生年份達(dá)50%以上[1]。研究者們對(duì)于假眼小綠葉蟬生物學(xué)習(xí)性、地理分布、危害損失估計(jì)和防治等進(jìn)行了深入的研究[2- 3]，但關(guān)于其地理種群遺傳和分化情況的報(bào)道甚少，僅應(yīng)用 RAPD-PCR 技術(shù)鑒定了茶園小綠葉蟬類的優(yōu)勢(shì)種[4]。我國(guó)幅員遼闊，茶樹(shù)廣泛種植于秦嶺淮河以南的祖國(guó)南半壁，各茶區(qū)的自然地理?xiàng)l件差別大，茶樹(shù)種植制度不同，可以造就不同地理區(qū)域的種群或生態(tài)型；另一方面，假眼小綠葉蟬成蟲(chóng)善爬善跳還具有一定的飛行能力，若蟲(chóng)也善于爬行和跳躍，成、若蟲(chóng)能適應(yīng)各種茶園生境。如果這種非遠(yuǎn)距離的擴(kuò)散連續(xù)地或漸次地發(fā)生，則有可能使得不同地理區(qū)域的種群發(fā)生基因交流[5]。

隨著分子系統(tǒng)地理學(xué)的發(fā)展，分子標(biāo)記技術(shù)越來(lái)越多地應(yīng)用到分子系統(tǒng)學(xué)研究中。mtDNA具有嚴(yán)格遵循母系遺傳，不發(fā)生重組、倒位、易位等突變現(xiàn)象以及進(jìn)化速率快等特點(diǎn)，而常被用于分析昆蟲(chóng)物種間及物種內(nèi)種群的分化程度[6- 7]。COI基因是線粒體蛋白編碼基因之一，密碼子第三位堿基可相對(duì)自由變異[8]。且COI的演化速率是其他線粒體蛋白編碼基因的3倍，兩端序列相對(duì)保守，是較好的分子標(biāo)記[9- 10]，已被用于研究玉米黃翅葉蟬Dalbulusmaidis、黃條脊冠葉蟬Aphrodesleafhoppers等葉蟬的遺傳多樣性[11- 12]。

目前主要以COI基因的部分序列進(jìn)行遺傳多樣性研究。本實(shí)驗(yàn)首次擴(kuò)增得到我國(guó)13個(gè)主要茶葉種植省份的茶園假眼小綠葉蟬線粒體COI基因的全長(zhǎng)序列，通過(guò)滑動(dòng)窗口分析確定其突變位點(diǎn)的分布，并發(fā)現(xiàn)了茶園假眼小綠葉蟬線粒體COI基因存在不完整的終止密碼子。再利用此序列分析了茶園假眼小綠葉蟬種群間的遺傳多樣性、遺傳分化程度、基因流水平及分子變異情況等，以探討長(zhǎng)期的地理阻隔對(duì)其種群遺傳特性的影響及其生態(tài)適生性、以及該葉蟬自身擴(kuò)散和人工助遷等因素造成的種群變異，為針對(duì)不同地理區(qū)域假眼小綠葉蟬種群的控制策略提供參考。

1 材料和方法

1.1 供試蟲(chóng)源

本實(shí)驗(yàn)所用茶園假眼小綠葉蟬，在2011年10月—2013年10月期間采集于山東青島、河南西九華山、江蘇蘇州、安徽敬亭山、湖北英山、四川樂(lè)山、重慶永川、江西婺源、浙江蘭溪、湖南長(zhǎng)沙、福建武夷山、廣西桂林、和云南勐海等13個(gè)產(chǎn)茶省的重點(diǎn)產(chǎn)茶縣市，詳情見(jiàn)表1。采集的成蟲(chóng)和若蟲(chóng)浸泡于無(wú)水乙醇中，放置于-20℃長(zhǎng)期保存?zhèn)溆谩?/p>

1.2 基因組DNA的提取

取單頭假眼小綠葉蟬，放在加有適量液氮的1.5 mL離心管中，用無(wú)菌的研磨棒磨成粉末，采用動(dòng)物組織快速提取試劑盒(VWI)，按照試劑盒說(shuō)明書(shū)步驟進(jìn)行基因組DNA提取。從提取的DNA樣品中取1 μL用Nanodrop 2000c 儀器(Thermo Scientific公司)和1%瓊脂糖凝膠電泳檢測(cè)其濃度和純度，基因組DNA保存于-20℃?zhèn)溆谩?/p>

1.3 CO I基因全長(zhǎng)擴(kuò)增和序列測(cè)定

根據(jù)葉蟬科線粒體COI基因兩端的保守區(qū)域，設(shè)計(jì)特異性引物用于擴(kuò)增COI基因全長(zhǎng)序列(擴(kuò)增產(chǎn)物包括COI基因兩端的側(cè)翼序列，5′包括基因NADH 脫氫酶亞基2的部分序列，以及tRNATrp、tRNACys、tRNATyr3個(gè)相連的tRNA序列；3′端包括tRNALeu(UUR)序列以及COII基因部分序列)。引物序列：COI-F(正向引物)：5′- AACAACAGC TTTCTAGAAATTTGC- 3′,COI-R(反向引物)：5′-TTCCATAGTTGGTGAAGAAGC- 3′。擴(kuò)增片段的長(zhǎng)度約為1773 bp。

表1 假眼小綠葉蟬地理種群采集信息及供試個(gè)體數(shù)量

1.4 數(shù)據(jù)整理與分析

測(cè)序結(jié)果通過(guò)Chromas軟件讀取，用DNAMAN7.0對(duì)序列進(jìn)行編輯，去除兩端的側(cè)翼序列。用Clustal X 1.8軟件進(jìn)行多序列同源比對(duì)，并輔以人工校對(duì)，生成用于MEGA5.1軟件分析的數(shù)據(jù)格式[13]。將整理后的序列在NCBI網(wǎng)站上進(jìn)行BLAST比對(duì)，確定目標(biāo)序列。使用MEGA5.1軟件分析堿基組成、變異位點(diǎn)及堿基轉(zhuǎn)換顛換比；基于Kimura- 2-Parameter模型計(jì)算種群間的遺傳距離；以琉璃葉蟬HomalodiscavitripennisCOI基因全長(zhǎng)序列作為外群，采用NJ鄰接法構(gòu)建單倍型系統(tǒng)發(fā)育樹(shù)，系統(tǒng)發(fā)育樹(shù)各分支的置信度(bootstrap)均進(jìn)行1000次重復(fù)檢驗(yàn)[14]。采用NETWORK4.6.1.2軟件繪制單倍型間的進(jìn)化網(wǎng)絡(luò)圖[15]。利用DnaSP5.0軟件分析13個(gè)種群內(nèi)的單倍型多樣性Hd、核苷酸多樣性Pi、核苷酸平均差異數(shù)k，種群間核苷酸平均差異數(shù)Kxy、核苷酸歧義度Dxy、遺傳分化系數(shù)Gst與固定系數(shù)Fst等分子遺傳學(xué)系數(shù)，同時(shí)生成滑動(dòng)窗口，分析突變位點(diǎn)在全長(zhǎng)序列中的分布，并進(jìn)行Tajima′sD中性檢驗(yàn)[16- 17]。應(yīng)用Arlequin 3軟件進(jìn)行AMOVA分子變異分析，判斷遺傳分化的主要來(lái)源是種群內(nèi)或種群間[18]。通過(guò)GenALEx 6.41軟件檢驗(yàn)種群間的遺傳距離與地理距離的自然對(duì)數(shù)矩陣之間的相關(guān)性[19]。

2 結(jié)果與分析

2.1 茶園假眼小綠葉蟬線粒體CO I基因全長(zhǎng)序列及滑動(dòng)窗口分析

通過(guò)整理去除測(cè)序結(jié)果兩端側(cè)翼序列，得到茶園假眼小綠葉蟬線粒體COI基因全長(zhǎng)序列，為1534 bp。該基因沒(méi)有完整的終止密碼子，最后一個(gè)堿基為T(mén)。從已報(bào)道的研究中可知，許多半翅目昆蟲(chóng)，如琉璃葉蟬Homalodiscavitripennis(GenBank登錄號(hào)：NC_006899.1)、青蛾蠟蟬Geishadistinctissima(GenBank登錄號(hào)：NC_012617.1)[20- 21]以及紅頭鳳沫蟬Paphnutiusruficeps(GenBank登錄號(hào)：NC_021100.1)[22]，其線粒體COI基因以不完全終止密碼子(T)作為停止信號(hào)，與本研究所得結(jié)果一致。在13個(gè)種群176個(gè)個(gè)體中共發(fā)現(xiàn)突變位點(diǎn)113個(gè)，約占全長(zhǎng)的7.37%，其中有7個(gè)位點(diǎn)的突變引起了氨基酸的改變(359，798，967，1175，1316，1330，1422)。單一多態(tài)位點(diǎn)40個(gè)，簡(jiǎn)約信息位點(diǎn)73個(gè)，未發(fā)現(xiàn)堿基缺失或插入現(xiàn)象。由表2可知，在所有核苷酸替換中，轉(zhuǎn)換的發(fā)生頻率占87.61%，而顛換僅占12.39%。其中腺嘌呤(A)與鳥(niǎo)嘌呤(G)之間的轉(zhuǎn)換占總轉(zhuǎn)換數(shù)的30.8%，胞嘧啶(C)與胸腺嘧啶(T)之間的轉(zhuǎn)換占56.81%，值得注意的是，T轉(zhuǎn)換為C的比例高達(dá)43.68%。總體轉(zhuǎn)換/顛換偏倚率R值為7.01。據(jù)研究發(fā)現(xiàn)，通常親緣關(guān)系較近的分類階元之間核苷酸替換主要為轉(zhuǎn)換，而在親緣關(guān)系較遠(yuǎn)的分類階元之間，核苷酸替換率則以顛換為主[23- 24]。本研究發(fā)現(xiàn)茶園假眼小綠葉蟬種內(nèi)線粒體COI基因核苷酸突變位點(diǎn)的變異類型以轉(zhuǎn)換為主，驗(yàn)證了上述規(guī)律。

表2 假眼小綠葉蟬CO I基因堿基替換情況

通過(guò)軟件DnaSP5.0滑動(dòng)窗口分析COI全長(zhǎng)序列，結(jié)果見(jiàn)圖1。從圖中可以看到，序列間的突變位點(diǎn)在絕大部分區(qū)段內(nèi)均有分布。其中，在400 bp左右、700 bp到800 bp之間，1200 bp和1300 bp左右有較多的突變位點(diǎn)；前200 bp，900 bp左右，1100 bp左右的突變位點(diǎn)較少。由此可見(jiàn)，與使用基因的部分序列進(jìn)行遺傳分化分析相比，使用序列的全長(zhǎng)能夠包含更豐富的信息。

圖1 序列間突變位點(diǎn)的滑動(dòng)窗口分析

2.2 茶園假眼小綠葉蟬CO I基因單倍型系統(tǒng)進(jìn)化分析

從已獲得的176條COI基因全長(zhǎng)序列中定義了105種單倍型：Hap1—Hap105。其中Hap5為共享單倍型，在所有種群中出現(xiàn)24次，占所有種群的13.64%。其余頻率較高的單倍型有：Hap16出現(xiàn)13次(7.39%)，Hap3和Hap26出現(xiàn)9次(5.11%)。在105種單倍型中，有84種單倍型是由單個(gè)種群所特有的獨(dú)享單倍型。由單倍型系統(tǒng)發(fā)育樹(shù)(圖2，圖中僅顯示大于30%的置信度值bootstrap)可以看出，在幾個(gè)大分支中置信度較低，說(shuō)明這些單倍型之間的差異較小，不能形成可靠的分支。而各種群的單倍型分散在不同的分布群中，未形成明顯的系統(tǒng)地理結(jié)構(gòu)。圖3是茶園假眼小綠葉蟬mtDNACOI基因105 種單倍型間的進(jìn)化網(wǎng)絡(luò)圖，表現(xiàn)了各地理種群?jiǎn)伪缎偷姆植记闆r，以及單倍型之間的演化關(guān)系。每種顏色代表一個(gè)種群的單倍型組成。沒(méi)有單倍型存在的那些節(jié)點(diǎn)是單倍型進(jìn)化過(guò)程中的中間節(jié)點(diǎn)。單倍型Hap71在進(jìn)化網(wǎng)絡(luò)的較為中心位置，而其他幾種出現(xiàn)頻率較高的單倍型則處于邊緣。由圖3中可以看到，各地理種群的單倍型散布于整個(gè)網(wǎng)絡(luò)圖中，每個(gè)種群由多種較為原始的單倍型演化產(chǎn)生的單倍型共同組成。說(shuō)明種群間存在基因交流，分化程度較小。

2.3 茶園假眼小綠葉蟬不同地理種群的遺傳多樣性、種群間的基因流及遺傳分化分析

各種群COI基因單倍型、核苷酸多樣性分析及Tajima′sD中性檢驗(yàn)結(jié)果如表2所示。由表可知，13個(gè)種群的總體單倍型多樣性指數(shù)Hd為0.9720，種群間核苷酸多樣性指數(shù)Pi為0.00607。種群內(nèi)核苷酸多樣性在0.00413—0.00651范圍內(nèi)，平均為0.00546。其中單倍型多樣性最高的種群是福建武夷山(FJWYS)，而云南勐海(YNMH)最低。核苷酸多樣性最高和最低的種群分別為湖南長(zhǎng)沙(HNCS)和安徽敬亭山(AHJTS)。

圖2 利用鄰接法構(gòu)建的假眼小綠葉蟬不同地理種群CO I單倍型間的系統(tǒng)發(fā)育樹(shù)

對(duì)各個(gè)地理種群的茶園假眼小綠葉蟬COI基因全長(zhǎng)序列基于Tajima′sD值進(jìn)行中性檢驗(yàn)(表3)，檢測(cè)結(jié)果均不顯著，全部種群均為P>0.10，說(shuō)明茶園假眼小綠葉蟬種群的進(jìn)化過(guò)程遵循中性進(jìn)化模型，各地理種群在進(jìn)化過(guò)程中沒(méi)有出現(xiàn)群體擴(kuò)張，群體大小保持相對(duì)穩(wěn)定。

表3 假眼小綠葉蟬CO I基因單倍型、核苷酸多樣性分析及Tajima′s D中性檢驗(yàn)

圖3 105 種 mtDNA CO I基因 單倍型的進(jìn)化網(wǎng)絡(luò)關(guān)系

應(yīng)用DnaSP5.0軟件計(jì)算反映茶園假眼小綠葉蟬種群間基因差異程度的遺傳學(xué)參數(shù)。各種群間核苷酸差異數(shù)Kxy在6.47778—13.37500之間，均值為8.98383；核苷酸歧義度Dxy在0.00422—0.00872之間，均值為0.00586(表4)。各地理種群總遺傳分化系數(shù)Gst為0.03652，總固定系數(shù)Fst為0.10876，總基因流Nm為4.097，固定系數(shù)Fst在-0.04168—0.43061之間。其中云南勐海(YNMH)種群與其他種群的遺傳分化程度最大(表5)。說(shuō)明種群間分化程度低，基因交流頻繁。

通過(guò)AMOVA分子變異分析結(jié)果可以看出，不同地理種群之間的變異方差組分為0.59490(占方差比率的12.66%)，而種群內(nèi)的方差組分為4.10510(占方差比率的87.33%)且差異顯著(P<0.05)(表6)。說(shuō)明這13個(gè)茶園假眼小綠葉蟬地理種群的遺傳分化主要來(lái)自于種群內(nèi)部，種群間的差異較小。

從表7可知，各個(gè)地理種群之間的遺傳距離在0.0042—0.0088之間。應(yīng)用Mantel相關(guān)性檢驗(yàn)計(jì)算兩個(gè)矩陣之間的相關(guān)系數(shù)r為-0.0697(P=0.5373>0.05，1000次隨機(jī)抽樣)，表明這兩個(gè)矩陣之間沒(méi)有顯著的相關(guān)性，即不同種群之間的遺傳分化程度與地理隔離沒(méi)有顯著關(guān)系。

3 討論

同行專家們也發(fā)表了少量關(guān)于茶園葉蟬類種群分化和基因流的研究報(bào)道。在印度茶小綠葉蟬EmpoascaflavescensF. 是茶樹(shù)上的優(yōu)勢(shì)種，印度孟加拉邦的茶小綠葉蟬因茶園管理模式的不同引起了種群的遺傳分化[26]。在普通茶園中主要依靠噴施化學(xué)農(nóng)藥進(jìn)行防治，而有機(jī)茶園主要依靠茶樹(shù)品種的抗蟲(chóng)性與施用微生物農(nóng)藥防治，在這兩種管理模式下，茶小綠葉蟬面臨不同的選擇壓力，造成了不同管理模式的種群間差異大于相同管理模式下的種群間差異。在我國(guó)也有關(guān)于茶園假眼小綠葉蟬種群分化的初步研究，采用RAPD分子標(biāo)記技術(shù)，探討了我國(guó)7個(gè)省份的茶園假眼小綠葉蟬的遺傳差異及其親緣關(guān)系[4]。

基于前人的相關(guān)成果，本研究采用線粒體COI基因全長(zhǎng)序列，對(duì)我國(guó)主要產(chǎn)茶省份的茶園假眼小綠葉蟬不同地理種群的遺傳多樣性、遺傳分化程度、基因流水平及分子變異進(jìn)行了解析，建立了單倍型系統(tǒng)發(fā)育樹(shù)和進(jìn)化網(wǎng)絡(luò)圖。依據(jù)幾十年來(lái)的茶園葉蟬類發(fā)生和為害觀察、標(biāo)本收集及鑒定，茶園中約有10余種葉蟬，但優(yōu)勢(shì)種是假眼小綠葉蟬，其個(gè)體數(shù)占茶園葉蟬個(gè)體數(shù)的99.99%以上[27]，這個(gè)百分率很穩(wěn)定，不隨時(shí)間、地點(diǎn)的變化而改變，本研究的供試蟲(chóng)態(tài)都是成蟲(chóng)，供試標(biāo)本采集時(shí)間從2011年10月至2013年10月，采到的標(biāo)本立即浸入無(wú)水乙醇中，數(shù)小時(shí)之后就保存于-20℃?zhèn)溆?，可盡量消除采集時(shí)間因素對(duì)于試驗(yàn)結(jié)果的影響。在13個(gè)種群176個(gè)個(gè)體中共發(fā)現(xiàn)突變位點(diǎn)113個(gè)，核苷酸替換以轉(zhuǎn)換為主。突變位點(diǎn)在整個(gè)區(qū)段內(nèi)均有分布，但不同位置的進(jìn)化速率有所不同。共產(chǎn)生105種單倍型，具有較高的單倍型多態(tài)性。其中單倍型Hap5為共享單倍型，這種單倍型可認(rèn)為是能夠適應(yīng)環(huán)境變化并在種群中穩(wěn)定存在的優(yōu)勢(shì)單倍型，該單倍型很有可能就是祖先單倍型[28- 29]；在所有的單倍型中，有84種是獨(dú)享單倍型，是由單個(gè)種群所特有的。各地理種群內(nèi)具有較為豐富的單倍型多樣性，表明茶園假眼小綠葉蟬對(duì)于不同環(huán)境具有較強(qiáng)的生存能力，這與假眼小綠葉蟬廣泛存在于各種茶園環(huán)境中的情況一致。從單倍型系統(tǒng)發(fā)育樹(shù)來(lái)看，單倍型的分布與其所在的地理單元之間不存在明顯的對(duì)應(yīng)關(guān)系。Tajima′sD檢驗(yàn)是基于種內(nèi)多態(tài)性的一種中性檢驗(yàn)方法，可反映出物種種群變化動(dòng)態(tài)的歷史。統(tǒng)計(jì)值為正值時(shí)說(shuō)明序列進(jìn)化方式為平衡選擇，且有一些單倍型分化；負(fù)值時(shí)為負(fù)向選擇，種群的擴(kuò)張存在瓶頸效應(yīng)。本實(shí)驗(yàn)中總?cè)后w和各個(gè)地理種群的檢測(cè)結(jié)果均不顯著，說(shuō)明茶園假眼小綠葉蟬種群的進(jìn)化過(guò)程遵循中性進(jìn)化模型，過(guò)去沒(méi)有出現(xiàn)群體擴(kuò)張和持續(xù)性增長(zhǎng)，群體大小保持相對(duì)穩(wěn)定[17]。固定系數(shù)Fst是一種種群間遺傳距離的測(cè)度參數(shù)[30]，從本實(shí)驗(yàn)的Fst值上看不同地理種群之間存在一定程度的遺傳分化?；蛄骷慈后w間基因交流，是種群遺傳結(jié)構(gòu)均質(zhì)化的主要因素之一，具有高水平基因流的種群間往往比具有有限基因流的種群間的遺傳分化程度小[31]。本實(shí)驗(yàn)中總基因流Nm為4.097>4，當(dāng)Nm>4時(shí)，種群間的基因交流比較充分，均質(zhì)化作用足以抵制遺傳漂變的作用，一定程度上抵消種群間的遺傳分化[32]。通過(guò)AMOVA分子變異分析可以看出，假眼小綠葉蟬的遺傳分化主要來(lái)自于種群內(nèi)部，種群間的差異較小。通過(guò)Mnatel檢驗(yàn)可知，不同種群之間的遺傳分化程度與地理距離沒(méi)有顯著關(guān)系。

表4 假眼小綠葉蟬成對(duì)種群間的核苷酸平均差異數(shù) (Kxy) (上三角) 與核苷酸歧義度(Dxy)(下三角)

表5 假眼小綠葉蟬成對(duì)種群間的遺傳分化系數(shù) (Gst) (上三角) 與固定系數(shù)(Fst)(下三角)

表6 假眼小綠葉蟬13個(gè)地理種群CO I基因全長(zhǎng)序列的分子變異分析

表7 假眼小綠葉蟬13個(gè)地理種群的地理距離自然對(duì)數(shù)值(km) (上三角)與遺傳距離(下三角)

綜上所述，認(rèn)為假眼小綠葉蟬種群的進(jìn)化過(guò)程遵循中性進(jìn)化模型，群體大小保持相對(duì)穩(wěn)定，這與茶園中假眼小綠葉蟬每年保持高發(fā)生量的情況相吻合，因?yàn)樵S多茶區(qū)每年假眼小綠葉蟬發(fā)生時(shí)間、發(fā)生量和施藥次數(shù)基本不變，每年都是10多次。假眼小綠葉蟬種群間存在一定程度的遺傳分化，這是因?yàn)槲覈?guó)幅員廣大，不同茶園所在的自然地理?xiàng)l件存在較大差異，不同地理種群之間的假眼小綠葉蟬可能存在著生物學(xué)和生態(tài)學(xué)特征的差異而有一定遺傳分化。而較為充分的基因交流帶來(lái)的均質(zhì)化作用，在一定程度上抵消了種群間的遺傳分化。種群間假眼小綠葉蟬充分的基因交流主要源于假眼小綠葉蟬自然遷移和人為傳播兩方面的原因。假眼小綠葉蟬成蟲(chóng)本身具有一定的飛行擴(kuò)散能力，能適應(yīng)各種地理環(huán)境，這種近距離的擴(kuò)散行為一直在進(jìn)行中，能夠引起不同種群間一定程度的基因交流。相比于茶園假眼小綠葉蟬的自然遷移，人為助遷的作用十分突出。在許多茶人們的印象中，20世紀(jì)50年代省際茶區(qū)之間假眼小綠葉蟬的體色等特征是有差異的，可能意味著大尺度的自然地理景觀條件下該葉蟬種群存在生態(tài)型的差異。然而，20世紀(jì)90年代以后，省內(nèi)茶區(qū)之間、省際之間茶樹(shù)鮮葉的異地加工十分活躍，幫助了茶園假眼小綠葉蟬的大規(guī)模遷移，比如，近些年來(lái)從早春至深秋，黔、川、渝、皖南、豫南茶區(qū)鮮葉源源不斷地調(diào)往浙、蘇茶區(qū)加工以獲得高價(jià)出售，這些鮮葉夾帶了大量的假眼小綠葉蟬成蟲(chóng)、若蟲(chóng)和卵，它們隨著茶樹(shù)鮮葉在1—3日內(nèi)遷移幾百—幾千km；閩南、閩北茶區(qū)茶苗不斷調(diào)往浙、蘇南、皖南、魯南和豫南等茶區(qū)，浙西茶區(qū)的安吉白茶茶苗持續(xù)調(diào)往周邊各省茶區(qū)。到達(dá)新區(qū)的假眼小綠葉蟬很快適應(yīng)了當(dāng)?shù)夭鑸@環(huán)境，經(jīng)過(guò)幾百代的世代更替，使得假眼小綠葉蟬不同地理種群生態(tài)型的差異逐漸消失。這也解釋了為何假眼小綠葉蟬的遺傳分化主要來(lái)自于種群內(nèi)部、種群間的差異較小、以及種群之間的遺傳分化程度與地理距離沒(méi)有顯著關(guān)系。

需要一提的是：來(lái)自西南邊陲的云南勐海種群(YNMH)在所有種群中單倍型多樣性Hd最小，與其他種群的遺傳分化程度最大，與其他種群之間的遺傳距離也最大，而除云南以外其他種群之間的差異比較一致。這說(shuō)明了云南的茶園假眼小綠葉蟬相比較于其他種群來(lái)說(shuō)處于一個(gè)相對(duì)封閉的環(huán)境中，與其他種群基因交流較小。這是因?yàn)樵颇系牟铇?shù)多為喬木型、半喬木型，茶園地理環(huán)境、氣候和植被、以及茶葉內(nèi)含物等顯著區(qū)別于其他假眼小綠葉蟬種群所在茶區(qū)的相應(yīng)因素。而且采集YNMH種群的云南的茶園處于原始森林，交通不便、人跡罕至，該假眼小綠葉蟬種群較為獨(dú)立，與其他種群的交流較少。

本研究首次擴(kuò)增得到茶園假眼小綠葉蟬COI基因的全長(zhǎng)序列，并將其應(yīng)用到種群間的遺傳分化研究中。與使用COI基因的部分序列進(jìn)行遺傳分化分析相比，使用序列的全長(zhǎng)包含了更豐富的信息，獲得了充分的變異位點(diǎn)及較為全面的遺傳學(xué)數(shù)據(jù)。在下一步的研究中，還需要配合其他分子標(biāo)記，比如線粒體的控制區(qū)以及核基因上的分子標(biāo)記，增加試驗(yàn)采樣點(diǎn)及樣本數(shù)量，進(jìn)一步加深對(duì)中國(guó)茶園假眼小綠葉蟬種群間關(guān)系及遺傳機(jī)制的研究。

:

[1] Wang Q S, Wang D F, Wu G Y. Research advances onEmpoascavitis(G?the) in tea trees in China. Fujian Journal of Agricultural Sciences, 2013, 28(6): 615- 623.

[2] Wang Y G, Li L J, Ran L X, Yu X S, Wu W W, Han B Y, Li Z M. Investigation on population dynamics ofEmpoascavitis(G?the) of different planting patterns. Southwest China Journal of Agricultural Science, 2010, 23(2): 413- 415.

[3] Liu L F, Xu D L, Mu D, Han B Y. Analysis of the feeding behaviours of tea green leafhoppers on resistant to sensible cultivars of tea plants by electrical penetration graphy techniques. Journal of Anhui Agricultural University, 2011, 38(2): 281- 286.

[4] Fu J Y, Han B Y. Studies on genetic relationships among populations ofEmpoascavitis(G?the) from tea gardens in seven provinces based on RAPD analysis. Acta Agriculturae Zhejiangensis, 2007, 19(1): 11- 14.

[5] Johnson K B, Radcliffe E B, Teng P S. Effects of interacting populations of Alternaria solani, Verticillium dahliae, and the potato leafhopper (Empoascafabae) on potato yield. Phytopathology, 1986, 76: 1046- 1052.

[6] Wang D M, Lin Y. Application of mitochondrial DNA gene order in molecular systematics of insects. Guangdong Journal of Agricultural Sciences, 2010, 37(6): 188- 190.

[7] Wang B X, Yang L F. Phylogenetic utilites of mitochondrial DNA sequences in the study of insect systematic. Entomological Knowledge, 2002, 39(2): 88- 92.

[8] Lv H J, Huang Y. Phylogenetic relationship among some groups of orthopteran based on complete sequences of the mitochondrial COI gene. Zoological Research, 33(3): 319- 328.

[9] Pavan-Kumar A, Gireesh-Babu P, Babu P S, Jaiswar A K, Krishna V H, Prasasd K P, Chaudhari A, Raje S G, Chakraborty S K, Krishna Gopal, Lakra W S. Molecular phylogeny of elasmobranchs inferred from mitochondrial and nuclear markers. Molecular Biology Reports, 2014, 41(1): 447- 457.

[10] Zanol J, Halanych K M, Struck T H, Fauchald K. Phylogeny of the bristle worm family Eunicidae (Eunicida, Annelida) and the phylogenetic utility of noncongruent 16S,COIand 18S in combined analyses. Molecular Phylogenetics and Evolution, 2010, 55(2): 660- 676.

[11] Palomera V, Bertin S, Rodríguez A, Bosco D, Virla E, Moya-Raygoza G. Is there any genetic variation among native mexican and argentinian populations ofDalbulusmaidis(Hemiptera: Cicadellidae)? Florida Entomologist, 2012, 95(1): 150- 155.

[12] luemel J K, Andrew K R, Virant-Doberlet M, Symondson W O C. Primers for identification of type and other archived specimens ofAphrodesleafhoppers(Hemiptera, Cicadellidae). Molecular Ecology Resources, 2011, 11(5): 770- 774.

[13] Thompson J D, Gibson T J, Plewniak F, Jeanmougin F, Higgins D G. The CLUSTAL_X windows interface: flexible strategies for multiple sequence alignment aided by quality analysis tools. Nucleic Acids Research, 1997, 25(24): 4876- 4882.

[14] Tamura K, Dudley J, Nei M, Kumar S. MEGA4: molecular evolutionary genetics analysis (MEGA) software version 4. 0. Molecular Biology and Evolution, 2007, 24(8): 1596- 1599.

[15] Polzin T, Daneshmand S V. On Steiner trees and minimum spanning trees in hypergraphs. Operations Research Letters, 2003, 31(1): 12- 20.

[16] Librado P, Rozas J. DnaSP v5: a software for comprehensive analysis of DNA polymorphism data. Bioinformatics, 2009, 25(11): 1451- 1452.

[17] Tajima F. Statistical method for testing the neutral mutation hypothesis by DNA polymorphism. Genetics, 1989, 123(3): 585- 595.

[18] Excoffier L, Laval G, Schneider S. Arlequin (version 3. 0): an integrated software package for population genetics data analysis. Evolutionary Bioinformatics Online, 2007, 1: 47- 50.

[19] Peakall R, Smouse P E. GENALEX 6: genetic analysis in Excel. Population genetic software for teaching and research. Molecular Ecology Notes, 2006, 6(1): 288- 295.

[20] Song N, Liang A P. The complete mitochondrial genome sequence ofGeishadistinctissima(Hemiptera: Flatidae) and comparison with other hemipteran insects. Acta Biochimica et Biophysica Sinica, 2009, 41(3): 206- 216.

[21] Bengtsson-Palme J, Ryberg M, Hartmann M, Branco S, Wang Z, Godhe A, De Wit P, Sánchez-García M, Ebersberger I, de Sousa F, Amend A, Jumpponen A, Unterseher M, Kristiansson E, Abarenkov K, Bertrand Y J K, Sanli K, Eriksson K M, Vik U, Veldre V, Nilsson R H. Improved software detection and extraction of ITS1 and ITS2 from ribosomal ITS sequences of fungi and other eukaryotes for analysis of environmental sequencing data. Methods in Ecology and Evolution, 2013, 4(10): 914- 919.

[22] Liu J, Liang A P. The complete mitochondrial genome of spittlebugPaphnutiusruficeps(Insecta: Hemiptera: Cercopidae) with a fairly short putative control region. Acta Biochimica et Biophysica Sinica, 2013, 45(4): 309- 319.

[23] Simon C, Frati F, Beckenbach A, Crespi B, Liu H, Flook P. Evolution, weighting, and phylogenetic utility of mitochondrial gene sequences and a compilation of conserved polymerase chain reaction primers. Annals of the Entomological Society of America, 1994, 87(6): 651- 701.

[24] Frati F, Simon C, Sullivan J, Swofford D L. Evolution of the mitochondrial cytochrome oxidase II gene in Collembola. Journal of Molecular Evolution, 1997, 44(2): 145- 158.

[25] Slatkin M. A measure of population subdivision based on microsatellite allele frequencies. Genetics, 1995, 139(1): 457- 462.

[26] Saha D, Mukhopadhyay A, Bahadur M. Genetic diversity ofEmpoascaflavescensFabricius (Homoptera: Cicadellidae), an emerging pest of tea from sub-himalayan plantations of West Bengal, India. Proceedings of the Zoological Society, 2012, 65(2): 126- 131.

[27] Zhao D X, Chen Z M, Cheng J A. Belongingness of tea leafhopper dominant species. Journal of Tea Science, 2000, 20(2): 101- 104.

[28] Wen L Y, Zhang L X, Liu N F. Phylogenetic relationship ofPerdixdauuricaeinferred from mitochondrial Cytochrome b gene. Zoological Research, 2005, 26(1): 69- 75.

[29] Li J B, Ren Z M. Genetic diversity among schlechtendalia chinensis individuals revealed by Cyt b sequences. Journal of Fudan University: Natural Science, 2009, 48(5): 680- 686.

[30] Rousset F. Genetic differentiation and estimation of gene flow from F-statistics under isolation by distance. Genetics, 1997, 145(4): 1219- 1228.

[31] Millar C I, Libby W J. Strategies for conserving clinal, ecotypic, and disjunct population diversity in widespread species // Genetics and Conservation of Rare Plants. New York: Oxford University Press, 1991: 149- 170.

[32] Whitlock M C, McCauley D E. Indirect measures of gene flow and migration: FST≠1/(4Nm+1). Heredity, 1999, 82(2): 117- 125.

參考文獻(xiàn):

[1] 王慶森, 王定鋒, 吳光遠(yuǎn). 我國(guó)茶樹(shù)假眼小綠葉蟬研究進(jìn)展. 福建農(nóng)業(yè)學(xué)報(bào), 2013, 28(6): 615- 623.

[2] 汪云剛, 李良靜, 冉隆珣, 玉香甩, 吳文偉, 韓寶瑜, 李忠美. 不同種植模式下茶假眼小綠葉蟬種群動(dòng)態(tài)的調(diào)查. 西南農(nóng)業(yè)學(xué)報(bào), 2010, 23(2): 413- 415.

[3] 劉麗芳, 徐德良, 穆丹, 韓寶瑜. EPG 技術(shù)分析不同品種茶樹(shù)抗假眼小綠葉蟬取食行為的差異. 安徽農(nóng)業(yè)大學(xué)學(xué)報(bào), 2011, 38(2): 281- 286.

[4] 付建玉, 韓寶瑜. 七省茶園假眼小綠葉蟬的 RAPD 分析及其親緣關(guān)系探討. 浙江農(nóng)業(yè)學(xué)報(bào), 2007, 19(1): 11- 14.

[6] 王德明, 林楊. 線粒體DNA基因序列在昆蟲(chóng)分子系統(tǒng)學(xué)研究中的應(yīng)用. 廣東農(nóng)業(yè)科學(xué), 2010, 37(6): 188- 190.

[7] 王備新, 楊蓮芳. 線粒體DNA序列特點(diǎn)與昆蟲(chóng)系統(tǒng)學(xué)研究. 昆蟲(chóng)知識(shí), 2002, 39(2): 88- 92.

[27] 趙冬香, 陳宗懋, 程家安. 茶小綠葉蟬優(yōu)勢(shì)種的歸屬. 茶葉科學(xué), 2000, 20(2): 101- 104.

[28] 文隴英, 張立勛, 劉迺發(fā). 以mtDNA細(xì)胞色素b基因探討斑翅山鶉的分類地位. 動(dòng)物學(xué)研究, 2005, 26(1): 69- 75.

[29] 李繼變, 任竹梅. 角倍蚜mtDNA Cyt b 基因遺傳多樣性分析. 復(fù)旦學(xué)報(bào): 自然科學(xué)版, 2009, 48(5): 680- 686.