北二星蝽Eysarcoris aeneus的DNA分類學(xué)研究(半翅目：蝽科：二星蝽屬)

趙清，李敏，孫溪，張虎芳*

(1.山西農(nóng)業(yè)大學(xué) 農(nóng)學(xué)院，山西 太谷 030801；2.太原師范學(xué)院，山西 太原 030000；3.南開大學(xué) 昆蟲學(xué)研究所，天津 300071)

?

北二星蝽Eysarcoris aeneus的DNA分類學(xué)研究(半翅目：蝽科：二星蝽屬)

趙清1，李敏2，孫溪3，張虎芳1*

(1.山西農(nóng)業(yè)大學(xué) 農(nóng)學(xué)院，山西 太谷 030801；2.太原師范學(xué)院，山西 太原 030000；3.南開大學(xué) 昆蟲學(xué)研究所，天津 300071)

本文基于形態(tài)和分子數(shù)據(jù)證明北二星蝽Eysarcorisaeneus和尖角二星蝽E.parvus為同一個(gè)種。通過(guò)擴(kuò)增中國(guó)二星蝽屬的4個(gè)種，尖角二星蝽E.parvus，北二星蝽E.aeneus，二星蝽E.guttigerus和廣二星蝽E.ventralis的線粒體COI基因部分片段(約560 bp)，共獲得28條序列。計(jì)算其種間和種內(nèi)遺傳距離，并基于鄰接法(NJ)、最大似然法(ML)和貝葉斯推斷法(BI)進(jìn)行系統(tǒng)進(jìn)化分析；同時(shí)，解剖這些種的生殖節(jié)發(fā)現(xiàn)，尖角二星蝽和北二星蝽之間的生殖節(jié)差異遠(yuǎn)小于與本屬內(nèi)其他種間的生殖節(jié)差異，因此，基于遺傳距離、系統(tǒng)發(fā)育分析以及生殖節(jié)形態(tài)差異的結(jié)果均顯示北二星蝽和尖角二星蝽之間的差異遠(yuǎn)小于這兩個(gè)種與屬內(nèi)其他種之間的差異。我們得出(1)北二星蝽和尖角二星蝽為同一種；(2)DNA分類學(xué)可以利用線粒體COI基因解決有爭(zhēng)議的分類學(xué)問題，同時(shí)綜合利用形態(tài)和分子數(shù)據(jù)能夠?yàn)榇祟悊栴}的解決提供新的方法和視角。

半翅目；蝽科；二星蝽屬；北二星蝽；DNA分類學(xué)

二星蝽屬(EysarcorisHahn,1834)隸屬于半翅目Hemiptera異翅亞目Heteroptera蝽科Pentatomidae蝽亞科Pentatominae。二星蝽成蟲、若蟲可通過(guò)吸食寄主植物的嫩枝、幼莖或葉片的汁液，致植株生長(zhǎng)發(fā)育受阻，籽粒不飽滿，部分種類可危害水稻、麥類、棉花、大豆、胡麻、高粱、甘薯、茄子等作物，是重要的農(nóng)業(yè)害蟲[1]。據(jù)世界主要產(chǎn)稻國(guó)家研究報(bào)道，所有為害稻類的昆蟲中，二星蝽屬的廣二星蝽E.ventralis，北二星蝽E.aeneus和E.trimaculatus位居前列[2,3]。本屬昆蟲分布較廣，在蝽科中由于體小，小盾片基角處具黃白色光滑小圓斑而極易鑒定[4]。

其中，北二星蝽E.aeneus和尖角二星蝽E.parvus在形態(tài)上極為相近，唯有前胸背板側(cè)角不同(圖1中1～8)，前者側(cè)角末端圓鈍且伸出較短，后者末端尖銳呈針狀且伸出較長(zhǎng)；生殖器差異較小，僅在陽(yáng)莖的系膜頂葉、系膜側(cè)葉的長(zhǎng)短、骨化程度有區(qū)別(圖2 Ie，IIe)；另外，這2種在有些地方為同域分布(表1)。因此該2種的分類地位一直有爭(zhēng)議，Josifov(1978)將此2種合并為北二星蝽E.aeneus，將尖角二星蝽E.parvus做次異名處理[5]，但是未提供詳細(xì)的證據(jù)。對(duì)于這些形態(tài)特征有變異的分類處理，除了借助外部形態(tài)特征及雄性外生殖器的比較外，還應(yīng)該提供其他證據(jù)，如分子證據(jù)。

在過(guò)去的幾年中，DNA條形碼作為一個(gè)分子生物學(xué)不可或缺的方法吸引了人們的注意力。它是利用標(biāo)準(zhǔn)化的DNA片段進(jìn)行物種鑒定的方法[6]，能夠幫助非分類學(xué)專家進(jìn)行物種鑒定。線粒體細(xì)胞色素C氧化酶基因COI具有引物通用性高、進(jìn)化速率快、變異速率在種內(nèi)相對(duì)小，在種間又相對(duì)大[6,7]等優(yōu)勢(shì)，其5’端648 bp的片段成為動(dòng)物分類與鑒別的理想DNA 條形碼[8～10]。另外，DNA條形碼在提高物種鑒定效率、發(fā)現(xiàn)新種和隱存種以及研究系統(tǒng)發(fā)育和進(jìn)化關(guān)系等方面提供新的思路和研究工具[11]。因此，DNA條形碼被越來(lái)越多地應(yīng)用于分類，但許多學(xué)者建議COI條形碼被用于新種描述或新異名時(shí)應(yīng)與形態(tài)、地理分布等相結(jié)合。

本文基于形態(tài)和分子證據(jù)探究北二星蝽和尖角二星蝽兩者的關(guān)系。形態(tài)學(xué)方面解剖雄性生殖節(jié)并進(jìn)行比較，分子生物學(xué)方面通過(guò)擴(kuò)增線粒體COI基因序列(約560 bp)，計(jì)算其種內(nèi)種間遺傳距離并構(gòu)建系統(tǒng)進(jìn)化樹。預(yù)期基于兩種數(shù)據(jù)的研究結(jié)果推斷(1)北二星蝽和尖角二星蝽是否為同一種；(2)DNA條形碼是否可以有效地解決有爭(zhēng)議的分類學(xué)問題。

圖1 北二星蝽、尖角二星蝽、廣二星蝽和二星蝽的成蟲背腹面觀Fig.1 Habitus of E.aeneus,E.parvus,E.ventralis and E.guttigerus,dorsal and ventral view注：1～2.北二星蝽(1.背面觀; 2.腹面觀); 3～4.尖角二星蝽(3.背面觀; 4.腹面觀); 5～6.廣二星蝽(5.背面觀; 6.腹面觀); 7～8.二星蝽(7.背面觀; 8.腹面觀)Note:1～2.E.aeneus (1.dorsal view; 2.ventral view); 3～4.E.parvus (3.dorsal view; 4.ventral view); 5～6.E.ventralis (5.dorsal view; 6.ventral view); 7～8.E.guttigerus (7.dorsal view; 8.ventral view)

1 材料與方法

1.1 形態(tài)學(xué)

雄性生殖器取自回軟的標(biāo)本，將其放入2%的NaOH溶液中并在燈下烤30 min，之后置于清水中數(shù)分鐘使陽(yáng)莖充分膨脹。觀察時(shí)需讓陽(yáng)莖膨脹充分至完全伸出陽(yáng)莖鞘之外，必要時(shí)可以輔助一些器械。觀察完之后，將生殖器保存于甘油當(dāng)中，和原標(biāo)本放在一起。成蟲的背腹面照片采用Nikon SMZ1000照相。所用標(biāo)本均來(lái)自于南開大學(xué)昆蟲標(biāo)本館。

1.2 分子生物學(xué)

選取尖角二星蝽16頭、北二星蝽8頭、廣二星蝽和二星蝽各2頭，試驗(yàn)中所得數(shù)據(jù)均來(lái)自于酒精浸制標(biāo)本提取(標(biāo)本詳細(xì)信息見表1)。

基因組的提取采用試劑盒(康為世紀(jì)CW0546)。證據(jù)標(biāo)本除胸部肌肉用于提取基因外，其他部分酒精保存，存放于南開大學(xué)昆蟲所。PCR擴(kuò)增采用如下引物:YT1F(5-AAACTATTAACCTTCAAAG-3)/HCO2198(5-GTCGTGGAAGAAGATTAGTTGTTCTAT-3)。

圖2 北二星蝽、尖角二星蝽、廣二星蝽和二星蝽的雄性生殖節(jié)圖Fig.2 Male genitalia of E.aeneus,E.parvus,E.ventralis and E.guttigerus注：I: 北二星蝽; II: 尖角二星蝽; III: 廣二星蝽; IV: 二星蝽; a: 生殖囊背面觀; b: 生殖囊腹面觀; c:陽(yáng)基側(cè)突端面觀; d: 陽(yáng)基側(cè)突側(cè)面觀; e: 陽(yáng)莖Note:I: E.aeneus; II: E.parvus; III: E.ventralis; IV: E.guttigerus; a: dorsal view of pygophore; b: ventral view of pygophore; c: apical view of paramere; d: lateral view of paramere; e: aedeagus

PCR擴(kuò)增的反應(yīng)體積為50 μL，2 μL上游和下游引物，2～5 μL的基因組模板，25 μL的Mix預(yù)混液(康為世紀(jì)CW0556)，加ddH2O補(bǔ)足至50 μL。PCR反應(yīng)為92 ℃變性40 s，45 ℃退火40 s，72 ℃延伸1 min共5個(gè)循環(huán)，之后再92 ℃變性40 s，51 ℃退火40 s，72 ℃延伸1 min，共33個(gè)循環(huán)，反應(yīng)前94 ℃預(yù)變性2 min，反應(yīng)后72 ℃再延伸8 min，最后4 ℃保溫。

經(jīng)PCR擴(kuò)增后，用1.0%的瓊脂糖凝膠電泳檢測(cè)，若條帶整齊明亮，將PCR產(chǎn)物送至金維智測(cè)序公司直接測(cè)序。

獲得數(shù)據(jù)后，使用BioEdit7.0[12]打開峰圖文件，檢查序列的單一性，對(duì)得到的每條序列進(jìn)行校對(duì)、拼接，并結(jié)合峰圖進(jìn)行修正，確定準(zhǔn)確無(wú)誤后將其另存為Fasta格式文件；然后導(dǎo)入Mega 6.0[13]軟件中，用無(wú)脊椎動(dòng)物線粒體遺傳密碼子表對(duì)COI基因進(jìn)行翻譯，適當(dāng)調(diào)整閱讀框檢查有無(wú)終止密碼子存在，以排除假基因的干擾。將所得序列在NCBI上進(jìn)行BLAST搜索，得到的同源序列大多數(shù)為蝽科線粒體基因組的COI片段。在BioEdit中，用ClastalW進(jìn)行比對(duì)，去掉兩端不整齊的序列。用DAMBE4.0.36[14]計(jì)算單倍型個(gè)數(shù)，不同的單倍型用于遺傳距離的計(jì)算和系統(tǒng)進(jìn)化樹的分析。用Mega6.0基于Kimura-2-Parameter[15]計(jì)算種間和種內(nèi)遺傳距離。

系統(tǒng)進(jìn)化樹基于鄰接法(NJ)[16]，最大似然法(ML)和貝葉斯推斷法(BI)構(gòu)建。最大似然法和鄰接法由PAUP* 4.10b10[17]windows版本執(zhí)行。貝葉斯法由MrBayes 3.1[18]運(yùn)算，最適合數(shù)據(jù)集的進(jìn)化模型利用Modeltest 3.7[19]測(cè)得采用GTR+I模型，運(yùn)行1 000 000代，每1 000代抽樣一次，“burn in”25%的樣本后構(gòu)建系統(tǒng)進(jìn)化樹。所構(gòu)建的進(jìn)化樹均采用TreeView1.6.6查看編輯。

表1 實(shí)驗(yàn)所用標(biāo)本的詳細(xì)信息

2 結(jié)果與分析

2.1 形態(tài)學(xué)數(shù)據(jù)

北二星蝽Eysarcorisaeneus(Scopoli,1763)(圖1)

CimexaeneusScopoli,1763: 122.

CimexfucatusRossi,1790: 235.

CimexperlatusFabricius,1794: 125.

Eysarcorisperlatusvar.spinicollisPuton,1881: 76; 1881: 56

Eysarcorisperlatusvar.ventralisHorvath,1882: 219.

EysarcorisparvusUhler,1896: 258.

Eysarcorisaeneusvar.peeziTamanini,1961: 122.

Eysarcorisaeneusf.nigriventrisStichel,1962: 781.

鑒別特征：頭側(cè)葉稍長(zhǎng)于中葉，具銅綠色金屬光澤。胝區(qū)具一黑斑，方形。前胸背板側(cè)角略伸出，或呈尖刺狀伸出。小盾片倒三角形，基角各具一長(zhǎng)橢圓形的斜置黃白斑，末端外緣常具3個(gè)小黑斑。腹部腹面中央均勻漆黑，邊緣呈鋸齒狀，兩側(cè)各具一長(zhǎng)短不等的黑縱紋。

形態(tài)特征：體卵圓形，淡黃褐色，密被黑色粗刻點(diǎn)。

頭：黑色，具銅綠色光澤，端部略下傾，基部中央具淺色短縱紋；側(cè)葉末端圓鈍，稍長(zhǎng)于中葉；復(fù)眼褐色，單眼紅色；觸角黃褐色，第IV節(jié)端部及第V節(jié)黑。頭部腹面黑，具銅綠色金屬光澤；喙黃褐色，末端伸達(dá)后足基節(jié)處。

胸：前胸背板寬大于長(zhǎng)，前半略下傾，后半隆拱；胝區(qū)具一近方形黑斑；側(cè)角末端圓鈍，略伸出，部分個(gè)體側(cè)角呈尖刺狀伸出較長(zhǎng)。小盾片倒三角形，基角具一斜列的長(zhǎng)橢圓形黃白斑；末端圓鈍，不超過(guò)前翅革片末端，邊緣常具3個(gè)小黑點(diǎn)斑。各胸部側(cè)板黃褐色，具黑色粗刻點(diǎn)；臭腺溝緣呈耳狀，超過(guò)后胸側(cè)板一半。足淡黃褐色，腿節(jié)亞端部具一小黑斑。

腹：腹部腹面淡黃褐色，基部中央無(wú)前指的突起；中央具黑色縱帶，約占腹部1/3，邊緣鋸齒狀，其兩側(cè)各具一長(zhǎng)短不等的黑縱紋。氣門黑色。

雄蟲生殖節(jié)(圖2)：生殖囊杯狀，寬大于長(zhǎng)，腹后緣中央呈寬“V”狀內(nèi)凹，兩側(cè)呈脊?fàn)罴雍瘢韵蚋姑娣?圖2 Ia，IIa,Ib，IIb)。陽(yáng)基側(cè)突較細(xì)長(zhǎng)，末端指狀(圖2 Id，IId)。陽(yáng)莖端較短，向端部漸細(xì)；系膜頂葉深二叉狀，內(nèi)側(cè)骨化；系膜側(cè)葉亦分叉，部分骨化(圖2 Ie，IIe)。

此種的前胸背板側(cè)角變異較大：部分個(gè)體側(cè)角末端圓鈍，伸出體外較短；部分個(gè)體側(cè)角末端呈尖刺狀，水平伸出較長(zhǎng)。解剖其雄性生殖節(jié)發(fā)現(xiàn)這兩類個(gè)體的生殖節(jié)差異與其他種比起來(lái)差異小很多，如生殖囊背腹后緣的彎曲程度(圖2：Ia，Ib，IIa，IIb)，陽(yáng)莖端的長(zhǎng)短、系膜頂葉的形狀及的長(zhǎng)短、系膜側(cè)葉的骨化程度及長(zhǎng)短(圖2：Ie，IIe)等在種間差異很大，而在這兩類個(gè)體之間變化相對(duì)小。

2.2 分子數(shù)據(jù)

2.2.1 遺傳距離分析結(jié)果

經(jīng)測(cè)序共獲得28條序列，用DAMBE統(tǒng)計(jì)得出有17個(gè)單倍型，其中尖角二星蝽16條序列10個(gè)單倍型，北二星蝽8條序列7個(gè)單倍型，而北二星蝽和尖角二星蝽的部分序列出現(xiàn)共享單倍型“HAP9”“HAP15”。利用MEGA6.0分析COI序列(558bp)中共有保守位點(diǎn)(Conserved)453個(gè)，變異位點(diǎn)(Variable)105個(gè)，簡(jiǎn)約信息位點(diǎn)(Parsim-info)101個(gè)，單突變子(Singleton)4個(gè)。統(tǒng)計(jì)堿基組成顯示，密碼子第一位T、C、A、G平均含量分別為40%、8.2%、48.3%、3.5%，密碼子第二位T、C、A、G平均含量分別為25%、17.5%、30.7%、26.9%，密碼子第三位T、C、A、G平均含量分別為45%、25.6%、13.4%、16.2%，而G+C平均含量為33.6%，遠(yuǎn)低于A+T的含量，符合線粒體基因的高A-T含量的特性。

利用MEGA6.0對(duì)COI基因基于K-2-P計(jì)算種內(nèi)/種間遺傳距離(表2)。由表2可以看出廣二星蝽與北二星蝽、尖角二星蝽、二星蝽的遺傳距離分別為0.163 7、0.165 8、0.110 2，而尖角二星蝽和北二星蝽的種間遺傳距離為0.0119。根據(jù)Hebertetal.(2004)提出的標(biāo)準(zhǔn)，種間遺傳距離應(yīng)大于種內(nèi)距離的10倍，所以尖角二星蝽和北二星蝽應(yīng)該為同一個(gè)種。

表2 基于COI基因種內(nèi)-種間遺傳距離Table 2 The intra-and inter-species genetic distance(K2P)of COI

2.2.2 系統(tǒng)發(fā)育分析

我們選取本屬的廣二星蝽E.ventralis作為外群，系統(tǒng)進(jìn)化樹的結(jié)果顯示，基于NJ、ML和BI構(gòu)建的系統(tǒng)進(jìn)化樹結(jié)果基本一致，且節(jié)點(diǎn)支持率較高(圖3～圖5)。在各進(jìn)化樹中二星蝽E.guttigerus及早的分出，北二星蝽和尖角二星蝽聚為一大支，且節(jié)點(diǎn)支持率均較高，為100或1.00，雖支系內(nèi)部略有差異，但總體來(lái)說(shuō)無(wú)論是基于哪種算法的系統(tǒng)進(jìn)化結(jié)果均不支持二星蝽和尖角二星蝽的單系性。

圖3 基于COI基因構(gòu)建的Bayesian系統(tǒng)進(jìn)化樹Fig.3 Phylogenetic tree of Bayesianbased on COI gene

圖4 基于COI基因構(gòu)建的NJ系統(tǒng)進(jìn)化樹Fig.4 Phylogenetic tree of NJ based on COI gene

圖5 基于COI基因構(gòu)建的ML系統(tǒng)進(jìn)化樹Fig.5 Phylogenetic tree of ML based on COI gene

3 結(jié)論與討論

Josifov(1978)雖然將北二星蝽和尖角二星蝽合并為一個(gè)種，但是文中未提供詳細(xì)證據(jù)，亦未提供雄性外生殖節(jié)圖。本文通過(guò)解剖北二星蝽與尖角二星蝽的雄性生殖節(jié)，詳細(xì)比較兩種生殖節(jié)各部分結(jié)構(gòu)，并與廣二星蝽及二星蝽的雄性生殖節(jié)結(jié)構(gòu)進(jìn)行比較，結(jié)果顯示尖角二星蝽和北二星蝽的形態(tài)差異與同屬其他種相比差異極小，不足以達(dá)到種的界限。同時(shí)擴(kuò)增了二星蝽屬內(nèi)4個(gè)種的COI基因，計(jì)算其種內(nèi)種間遺傳距離，結(jié)果顯示北二星蝽與尖角二星蝽間的種間遺傳距離為0.0119，與本屬內(nèi)的種內(nèi)遺傳距離相似，且與其他種的種間遺傳距離在0.139 1～0.165 8之間，如此大的差距，給予Josifov(1978)提出的次異名處理提供了強(qiáng)有力的證據(jù)支持。另外，基于三種算法構(gòu)建的系統(tǒng)進(jìn)化樹結(jié)果均不支持兩者的單系性。因此，無(wú)論從形態(tài)數(shù)據(jù)還是分子數(shù)據(jù)我們均能證明尖角二星蝽和北二星蝽為同一種。

同時(shí)，該研究也證明，DNA分類學(xué)能夠?yàn)閭鹘y(tǒng)分類學(xué)服務(wù)。當(dāng)遇到形態(tài)數(shù)據(jù)與分子數(shù)據(jù)有異議時(shí)，我們需要重新審視鑒定結(jié)果，還應(yīng)該結(jié)合生態(tài)學(xué)、地理分布等來(lái)共同分析其物種界限。毫無(wú)疑問，DNA 條形碼可以幫助一些毫無(wú)頭緒的屬種縮小鑒定范圍；同時(shí)，傳統(tǒng)分類學(xué)又能從形態(tài)上對(duì)分子結(jié)果進(jìn)行檢測(cè)。在研究過(guò)程中，我們必須將兩種方法相互補(bǔ)充、完美結(jié)合才能最大限度發(fā)揮DNA分類學(xué)的優(yōu)勢(shì)。

[1]章士美.中國(guó)經(jīng)濟(jì)昆蟲志(第50冊(cè))[M].北京: 科學(xué)出版社.1995:169.

[2]Lee JG,Hong SS,Kim JY,et al.Occurrence of stink bugs and pecky rice damage by stink bugs in paddy fields in Gyeonggi-do,Korea[J].Korean Journal of Applied Entomology，2009,48(1): 37-44.

[3]Nasiruddin M,Roy RC.Rice field insect pests during the rice growing seasons in two areas of Hathazari,Chittagong[J].Bangladesh Journal of Zoology,2012,40(1): 89-100.

[4]蕭采瑜.中國(guó)蝽類昆蟲鑒定手冊(cè),第I冊(cè)(半翅目,異翅亞目)[M].北京: 科學(xué)出版社,1977: 76-89.

[5]Josifov MV,Kerzhner IM,Heteroptera aus Korea.II.Teil (Aradidae,Berytidae,Lygaeidae,Pyrrhocoridae,Rhopalidae,Alydidae,Coreidae,Urostylidae,Acanthosomatidae,Scutelleridae,Pentatomidae,Cydnidae,Plataspidae)[J].Fragmenta Faunistica,1978,23(9): 137-196.

[6]Hebert PDN,Cywinska A,Ball SL,DeWaard JR.Biological identifications through DNA barcodes[J].Proceedings of the Royal Society of London Series B-Biological Sciences,2003a,270: 313-321.

[7]Hebert PDN,Ratnasingham S,deWaard JR.Barcoding animal life: cytochrome oxidase subunit 1 divergences among closely related species[J].Proceedings of the Royal Society of London Series B-Biological Sciences,2003b,270: S96-S99.

[8]Hebert PDN,Stoeckle MY,Zemlak TS,et al.Identification of birds through DNA barcodes[J].Plos Biology,2004,2(7): 1657-1663.

[9]Ashfaq M,Akhtar S,Khan AM,et al.DNA barcode analysis of butterfly species from Pakistan points towards regional endemism[J].Molecular Ecology Resources,2013,13(5): 832-843.

[10]Yang F,Shi ZY,Bai SL,et al.Comparative studies on species identification of Noctuoidea moths in two nature reserve conservation zones (Beijing,China) using DNA barcodes and thin-film biosensor chips[J].Molecular Ecology Resources,2014,14(1): 50-59.

[11]Agnarsson I,Kuntner M.Taxonomy in a changing world: seeking solutions for a science in crisis[J].Systematic Biology,2007,56: 531-539.

[12]Hall TA.Bioedit: a user-friendly biological sequences alignment editor and analysis program for Windows 95/98/NT [J].Nucleic Acids Symposium Series,1999,41: 95-98.

[13]Tamura K,Stecher G,Peterson D,et al.MEGA6: Molecular Evolutionary Genetics Analysis version6.0 [J].Molecular Biology and Evolution,2013,30(12):2725-2729.

[14]Xia X,Xie Z.DAMBE: software package for data analysis in molecular biology and evolution [J].Journal of Heredity,2001,92(4): 371-373.

[15]Song H,Buhay JE,Whiting MF,et al.Many species in one: DNA barcoding overestimates the number of species when nuclear mitochondrial pseudogenes are coamplified [J].Proceedings of the National Academy of Sciences of the United States of America,2008,105: 13486-13491.

[16]Saitou N,Nei M.The neighbor-joining method: a new method for reconstructing phylogenetic trees [J].Molecular Biology and Evolution,1987,4: 406-425.

[17]Swofford DL.PAUP*: Phylogenetic Analysis Using Parsimony (and other methods).Version 4 (beta 10) [M].Sinauer Associates,Sunderland,Massachusetts,2003：94-96.

[18]Huelsenbeck JP,Ronquist F.MrBayes: Bayesian inference of phylogeny [J].Bioinformatics,2001,17: 754-755.

[19]Posada D,Crandall KA.Modeltest: testing the model of DNA substitution [J].Bioinformatics,1998,14: 817-818.

(編輯：武英耀)

Study on DNA Taxonomy ofEysarcorisaeneus(Hemiptera: Pentatomidae:Eysarcoris)

Zhao Qing1,Li Min2,Sun Xi3,Zhang Hufang1*

(1.CollegeofAgronomy,ShanxiAgriculturalUniversity,TaiguShanxi030801,China;2.TaiyuanNormalUniversity,TaiyuanShanxi030000,China;3.CollegeofLifeSciences,NankaiUniversity,Tianjin300071,China)

This paper based on morphology and molecular data to prove thatEysarcorisaeneus andEysarcorisparvusare the same species.We performed on 28 specimens,amplified about 560bp segment of the mitochondrial cytochrome coxidase I gene (COI) of 4 species,E.parvus,E.aeneus,E.guttigerusandE.ventralis.Computed the inter- and intra-species genetic distances and constructed the phylogenetic trees based on NJ,ML and BI methods.At the same time,we dissected the male genitalia of these species,compared their differences.All the results of the genetic distances,phylogenetic analysis and genitalia morphological characters showed that the difference and the genetic distance betweenE.aeneusandE.parvusmuch less than the difference amongE.aeneus,E.ventralisandE.guttigerus.So we concluded (1)E.aeneusandE.parvusare the same species; (2) DNA barcoding can resolve the argued taxonomy problems using the partial segment of mitochondrialCOIgene.The fact that based on molecular and morphological data can provide the new ideas and new perspective to resolve the taxonomy.

Heteroptera;Pentatomidae;Eysarcoris;Eysarcorisaeneus;DNA taxonomy

22015-03-10

2015-04-04

趙清(1985-)，女(漢)，山東臨沂人，講師，博士后，研究方向：DNA分類學(xué)

*通訊作者：張虎芳，教授，碩士生導(dǎo)師；Tel: 0354-6288225；E-mail: zh_hufang@sohu com

國(guó)家自然科學(xué)基金(31440078);山西省教育廳高等學(xué)?？萍紕?chuàng)新基金(2014131);山西農(nóng)業(yè)大學(xué)科技創(chuàng)新基金(2014017)

S433.3

A

1671-8151(2015)03-0241-08