房守敏

(西華師范大學(xué)生命科學(xué)學(xué)院,四川南充 637002)

隨著分子生物學(xué)的發(fā)展,DNA序列用于昆蟲(chóng)和其他動(dòng)物系統(tǒng)發(fā)育和分子進(jìn)化的研究越來(lái)越受到人們的重視。除了一些快速進(jìn)化的核基因如無(wú)翅(wingless,Wg)[1]和延伸因子1α(elongation factor 1-alpha,EF-1α)[2]等能用于分子系統(tǒng)學(xué)研究外,線粒體DNA(mitochondrial DNA,mtDNA)由于其進(jìn)化速率快而備受關(guān)注。自Nass等[3]于1962年發(fā)現(xiàn)mtDNA以來(lái),近年來(lái)人們對(duì)其基因組結(jié)構(gòu)和進(jìn)化進(jìn)行了深入細(xì)致的研究,發(fā)現(xiàn)線粒體基因組分子量小、拷貝數(shù)多、復(fù)合母系遺傳、不發(fā)生重組、突變率高等特點(diǎn)。動(dòng)物和昆蟲(chóng)的線粒體基因往往比核基因的進(jìn)化速率分別快5～10倍[3]和2～9倍[4]。不同的線粒體基因有不同的進(jìn)化速率,通常進(jìn)化速率最大的是控制區(qū);蛋白編碼基因中的細(xì)胞色素氧化酶I(cytochrome oxidase I,COI)、COII和NADH脫氫酶亞基5(NADH dehydrogenase subunit 5,ND5)次之;核糖體RNA(ribosomal RNA,rRNA)的大小亞基(12S和16S)較保守;進(jìn)化速率最慢的是轉(zhuǎn)運(yùn)RNA(transfer RNA,tRNA)。由于線粒體進(jìn)化速率快和基因組相對(duì)較小等優(yōu)點(diǎn),mtDNA已被廣泛用于系統(tǒng)發(fā)育、分子進(jìn)化、昆蟲(chóng)分類(lèi)學(xué)和群體遺傳學(xué)等研究[5]。

鱗翅目(Lepidoptera)有40個(gè)以上的總科、約120個(gè)以上的科及18萬(wàn)種,是僅次于鞘翅目的第二大目,包括蝶類(lèi)和蛾類(lèi)。鱗翅目昆蟲(chóng)是完全變態(tài)生物,發(fā)育過(guò)程分為卵、幼蟲(chóng)、蛹及成蟲(chóng)四個(gè)時(shí)期,是研究昆蟲(chóng)擬態(tài)、與寄主植物關(guān)系、翅膀模式發(fā)育生物學(xué)和昆蟲(chóng)生理等研究的重要模式生物。目前,GenBank(http://www.ncbi.nlm.nih.gov/)中鱗翅目昆蟲(chóng)約有68條線粒體全基因組序列,其中47條來(lái)自蠶蛾總科(Bombycoidea)[6-9],這可能是由于家蠶是重要的經(jīng)濟(jì)昆蟲(chóng)并且也是鱗翅目昆蟲(chóng)的模式生物的緣故。另外,有成千上萬(wàn)條的鱗翅目線粒體基因序列被測(cè)定,這些序列已被廣泛應(yīng)用于群體遺傳結(jié)果和系統(tǒng)發(fā)生分析等研究。本文主要對(duì)線粒體DNA在鱗翅目昆蟲(chóng)分子系統(tǒng)學(xué)研究的應(yīng)用,以及不同線粒體DNA區(qū)段在鱗翅目昆蟲(chóng)不同分類(lèi)階元的適用性進(jìn)行了綜述,為今后鱗翅目昆蟲(chóng)系統(tǒng)發(fā)生分析提供參考。

1 鱗翅目昆蟲(chóng)mtDNA結(jié)構(gòu)特征

鱗翅目昆蟲(chóng)mtDNA與其他動(dòng)物一樣,是閉合雙連環(huán)狀DNA分子。目前GenBank中鱗翅目昆蟲(chóng)mt全基因組的大小為14 535 bp(玉米螟Ostrinia nubilalis,Accession no.AF442957)～16 094 bp(瑪愛(ài)杰鳳蝶Papiliomaraho,Accession no.NC_003395)。線粒體全基因組包括13個(gè)蛋白編碼基因(COI、COII和COIII;細(xì)胞色素氧化酶亞基b基因Cytb;ATPase復(fù)合物亞基ATP6和ATP8基因;ND1～6和ND4L)、12S rRNA、16S rRNA、22個(gè)tRNA及1個(gè)A+T富集區(qū)。A+T富集區(qū)有控制mtDNA的復(fù)制和轉(zhuǎn)錄的功能,又叫控制區(qū)(control region,CR)。線粒體基因組各基因間排列緊密,僅存在較短的基因間隔序列。例如中國(guó)野桑蠶mt基因組含有17個(gè)基因,間區(qū)僅391 bp[10]。而且線粒體基因不含內(nèi)含子,部分基因還存在重疊現(xiàn)象[10]。因此,線粒體堿基的使用節(jié)約、高效。

2 mtDNA應(yīng)用于鱗翅目昆蟲(chóng)分子系統(tǒng)學(xué)的研究

2.1 蛋白編碼基因

蛋白編碼基因是線粒體的重要組織部分,對(duì)于完成線粒體全基因組測(cè)序的鱗翅目昆蟲(chóng)而言,常用13個(gè)基因的核酸或氨基酸連接序列進(jìn)行分子進(jìn)化的研究[10-11]。但更多的是選取線粒體的部分代表性基因進(jìn)行相關(guān)的分子系統(tǒng)學(xué)的研究。

2.1.1 COI與COII基因

我們對(duì)PubMed(http://www.ncbi.nlm.nih.gov/)中截至2009年的文獻(xiàn)進(jìn)行分析與整理,發(fā)現(xiàn)鱗翅目昆蟲(chóng)分子系統(tǒng)學(xué)的研究中COI與COII基因的應(yīng)用最廣,至少有51項(xiàng)研究涉及利用了COI和(或)COII基因,而且COI基因的應(yīng)用更為廣泛(表1)。這些研究中還常常結(jié)合線粒體的其他基因或區(qū)段(16S rDNA、ND5和控制區(qū)等)及核基因Wg(wingless)和EF-1α(elongationfactor 1-alpha)共同揭示鱗翅目昆蟲(chóng)的進(jìn)化關(guān)系。表1還顯示,利用COI/COII基因分析鱗翅目蛺蝶科Nymphalidae、鳳蝶科Papilionidae和灰蝶科Lycaenidae等昆蟲(chóng)系統(tǒng)發(fā)生關(guān)系的研究較多。

表1 線粒體COI/COI I基因在鱗翅目昆蟲(chóng)分子系統(tǒng)學(xué)的研究中的應(yīng)用

續(xù)表1 線粒體COI/COII基因在鱗翅目昆蟲(chóng)分子系統(tǒng)學(xué)的研究中的應(yīng)用

COI基因已被較多地應(yīng)用于鱗翅目昆蟲(chóng)種內(nèi)、種間及屬間分子系統(tǒng)學(xué)研究。DeChaine和Martini[12]對(duì)落基山脈20個(gè)不同區(qū)域帕納塞斯蝶Parnassius smintheus的385個(gè)個(gè)體COI基因進(jìn)行了測(cè)序,進(jìn)化分析發(fā)現(xiàn)在落基山脈的北部和南部至少存在兩個(gè)多樣化的中心,并且利用巢式進(jìn)化支分析(nested clade analysis)揭示出該物種經(jīng)歷了多次的群體擴(kuò)增和分裂事件。Rand等[13]以COI基因?qū)业艭hrysoritis、Poecilmitis和Oxychaeta三個(gè)屬的昆蟲(chóng)進(jìn)行分子系統(tǒng)發(fā)生分析時(shí)發(fā)現(xiàn)3個(gè)屬應(yīng)歸屬于同一個(gè)Chrysoritis屬,這與基于形態(tài)特征的系統(tǒng)假設(shè)一致。Hwang等[14]分析發(fā)現(xiàn)家蠶COI基因與野桑蠶、日本柞蠶Antheraea yamamai和中國(guó)柞蠶A.pernyi的序列相似性分別為97%、85%和87%,基于COI構(gòu)建的系統(tǒng)發(fā)生樹(shù)表明大蠶蛾科和蠶蛾科呈單源系起。Behere等[15]以COI基因部分序列對(duì)夜娥科Noctuidae的昆蟲(chóng)作了進(jìn)化分析,表明澳洲棉鈴蟲(chóng)Helicoverpa punctigera在進(jìn)化上先于煙草夜娥H.assulta,而后相繼分化形成棉鈴蟲(chóng)H.armigera和玉米夜蛾H.zea,并且北美和南美大陸的H.zea從H.armigera祖先分歧年代非常近＜1.5百萬(wàn)年前(million years ago,Mya)。

COII基因在鱗翅目昆蟲(chóng)亞種及種間的分子進(jìn)化研究中有一定的應(yīng)用。Long等[16]對(duì)大二尾蛺蝶Polyura eudamippus 5個(gè)亞種的50個(gè)樣本COII基因的單倍型進(jìn)行了分析,結(jié)果表明在405bp的COII部分序列中共發(fā)現(xiàn)了9種單倍型,5個(gè)亞種的地理分布與單倍型空間分布基本一致,并且推測(cè)在冰河期以后大二尾蛺蝶采取了兩條線路的擴(kuò)張,一條是向北至四川、重慶和湖北,另一條則是向東至中國(guó)大陸的東南海和臺(tái)灣。Ohno等[17]以COII基因?qū)⒄让鳲strinia latipennis群體的兩個(gè)種O.latipennis(Warren)與O.ovalipennis Ohno進(jìn)行分子進(jìn)化研究,發(fā)現(xiàn)其親緣關(guān)系較近,約在0.3百萬(wàn)年前發(fā)生的分歧。王戎疆等[18]利用線粒體COII基因?qū)ξ餐惖麑?個(gè)種的12個(gè)樣品分析,在405bp長(zhǎng)的COII片段有11.4%的位點(diǎn)為多態(tài)性位點(diǎn),在進(jìn)化上尾蛺蝶屬蝴蝶分為兩大分支,一支包括大二尾蛺蝶、二尾蛺蝶和忘憂尾蛺蝶,另一分支包括窄斑鳳尾蛺蝶和黑鳳尾蛺蝶。這些分子系統(tǒng)學(xué)的結(jié)果均與形態(tài)學(xué)的結(jié)果相一致。

相對(duì)于COI或COII單基因研究而言,COI和COII基因的連接序列有著更廣泛的應(yīng)用,它可用于屬、種及亞種分類(lèi)階元的進(jìn)化分析,而且還常用于分歧年代的估算。1994年,Brown等[19]對(duì)絲蘭蛾科(Prodoxidae)Greya屬16個(gè)種的COI+COII基因765 bp序列作了進(jìn)化分析,其結(jié)果與形態(tài)分類(lèi)大致相似。1999年,Caterino和Sperling[20]初次對(duì)鳳蝶22個(gè)種和1個(gè)亞種的COI和COII基因測(cè)序,研究者分別以COI、COII和COI+COII的核酸序列作進(jìn)化分析,結(jié)果發(fā)現(xiàn)僅以COI基因不能支持garamas與scamander世系單系起源,COII則不能解決Papilionini科單系起源問(wèn)題,而COI+COII則能較好地消除這兩個(gè)問(wèn)題,因此非常適宜于鳳蝶亞屬及種間的分子進(jìn)化研究。2004年,Zakharov等[21]利用COI+COII對(duì)鳳蝶的51個(gè)種進(jìn)行了系統(tǒng)發(fā)育分析,參考考古學(xué)證據(jù),估算出了系統(tǒng)發(fā)生樹(shù)每個(gè)節(jié)點(diǎn)的分化年代,同時(shí)計(jì)算出了鳳蝶COI+COII的每年每位點(diǎn)的替換率為7.8×10-9～10.2×10-9,由此認(rèn)為鳳蝶COI+COII的進(jìn)化速率比昆蟲(chóng)的“標(biāo)準(zhǔn)”速率[22]慢2-4倍。2008年,Pan等[10]利用鳳蝶COI+COII的進(jìn)化速率作為分子鐘對(duì)蠶的進(jìn)化年代進(jìn)行了估算,結(jié)果表明家蠶與中國(guó)野桑蠶的分歧年代為1.08～1.41百萬(wàn)年。

COI/COII與其他一些線粒體基因和核基因結(jié)合對(duì)鱗翅目昆蟲(chóng)系統(tǒng)發(fā)育和分子進(jìn)化的研究也是非常有用的工具之一。其他的線粒體基因或區(qū)域主要包括12S rRNA、16S rRNA、ND5及控制區(qū);核基因則主要是EF-1α和Wg基因(表1)。2001年,Monteiro和Pierce[23]將COI、COII和EF-1α基因應(yīng)用于Bicyclus屬54個(gè)種的系統(tǒng)發(fā)育分析,他們分別建立了單基因和COI+COII+EF-1α基因連接序列的系統(tǒng)發(fā)育樹(shù),結(jié)果發(fā)現(xiàn)COI+COII+EF-1α基因連接序列的系統(tǒng)發(fā)育樹(shù)與形態(tài)分類(lèi)的結(jié)果較一致。在此后幾年里,有更多的研究人員將COI/COII和EF-1α/Wg等結(jié)合起來(lái)研究鱗翅目昆蟲(chóng)一些種屬間的關(guān)系(表1)。因此,mt基因與核基因結(jié)合也是研究昆蟲(chóng)系統(tǒng)發(fā)育和分子進(jìn)化一個(gè)切實(shí)可行的方法。

2.1.2 ND5基因

ND5基因是線粒體中快速進(jìn)化的基因之一,其基因長(zhǎng)度也比較適合分子進(jìn)化[24]。1999年,Yagi等[25]以ND5基因?qū)θ毡绝P蝶不同科的代表性物種進(jìn)化關(guān)系進(jìn)行研究時(shí),發(fā)現(xiàn)鳳蝶科Papilionidae與其他科形成不同的進(jìn)化類(lèi)群,而且在Papilionidae中分為兩個(gè)亞族,這一結(jié)論與傳統(tǒng)分類(lèi)一致。2001年,Yagi等[26]又以ND5基因?qū)θ毡救簫u和亞洲大陸東部的鳳蝶P.glacialis和P.stubbendorf ii的進(jìn)化關(guān)系作了詳細(xì)分析,發(fā)現(xiàn)這兩個(gè)種起源于共同的祖先,其大致分歧年代在1.7～2.0 Mya。Omoto等[27]利用ND5基因主要對(duì)絹蝶Parnassius的進(jìn)化關(guān)系作了探討,結(jié)果表明該屬為單系起源,其內(nèi)部可以分為多個(gè)類(lèi)群,即顯示在進(jìn)化過(guò)程中該物種發(fā)生了相對(duì)快速的擴(kuò)散。2002年,Yukuhiro等[7]以甲蟲(chóng)ND5基因的進(jìn)化速率(遺傳距離為0.01,分化時(shí)間等于3.6百萬(wàn)年)作為分子鐘,估算出了家蠶和日本野桑蠶分歧時(shí)間約為7.1百萬(wàn)年前。

除單獨(dú)用于進(jìn)化研究外,ND5基因也可結(jié)合其他mt基因同時(shí)用于鱗翅目昆蟲(chóng)分子系統(tǒng)學(xué)研究。2008年,Albre等[28]以ND5和COII基因?qū)rebia tyndarus種群的種和亞種的親緣關(guān)系進(jìn)行了研究,結(jié)果發(fā)現(xiàn)以ND5或COII的進(jìn)化樹(shù)拓?fù)浣Y(jié)構(gòu)基本一致,ND5+COII連接序列構(gòu)建的進(jìn)化樹(shù)表明Erebia tyndarus種群可以分為兩個(gè)不同的進(jìn)化枝,且鑒別出了11個(gè)新種。因此ND5基因已被用于鱗翅目昆蟲(chóng)的亞種、種、屬和科等分類(lèi)階元的分子系統(tǒng)學(xué)研究。

2.1.3 Cytb基因

2005年,Li等[29]利用Cytb基因探討了蠶蛾屬中的家蠶和野桑蠶的起源與進(jìn)化關(guān)系,結(jié)果發(fā)現(xiàn)14個(gè)家蠶品種與日本野桑蠶Cytb基因核苷酸水平的差異為5.4%～5.8%,而與中國(guó)鎮(zhèn)江野桑蠶的差異僅為0.8%～1.9%,進(jìn)化樹(shù)也顯示家蠶與中國(guó)野桑蠶的親緣關(guān)系更近,這與公認(rèn)的家蠶起源于中國(guó)野桑蠶一致。陳永久等[30]測(cè)定了4種云南白馬雪山和1種新疆天山的珍稀絹蝶Cytb基因部分序列,研究表明西猴絹蝶Parnassius simo分化較早,其次是珍珠絹蝶Parnassiusorleans和阿波羅絹蝶Parnassius apollo,巴裔絹蝶Parnassiusbaileyi和愛(ài)坷絹蝶Parnassius acco的遺傳關(guān)系較近,且分化較晚,與形態(tài)學(xué)結(jié)果相吻合。

2.2 核糖體RNA基因

線粒體核糖體RNA分為大小兩個(gè)亞基,其中16S rRNA應(yīng)用于分子系統(tǒng)學(xué)的研究相對(duì)較為廣泛。1999年,Hwang等[31]用12S和16S rRNA基因?qū)πQ蛾科與大蠶蛾科的親緣關(guān)系進(jìn)行研究,分別以12S和16S rRNA構(gòu)建的系統(tǒng)發(fā)生樹(shù)的拓?fù)浣Y(jié)構(gòu)相同,進(jìn)化樹(shù)顯示蠶蛾屬與大蠶蛾屬為單系起源。2008年,劉彥群等[32]以12S rRNA基因序列對(duì)鱗翅目柞蠶屬、樗蠶屬和桑蠶屬等9種絹絲昆蟲(chóng)之間的系統(tǒng)發(fā)生關(guān)系進(jìn)行了研究,從分子水平證實(shí)了所發(fā)現(xiàn)的野生柞蠶仍屬于中國(guó)柞蠶種,進(jìn)化分析支持柞蠶屬、樗蠶屬和桑蠶屬的單系起源。Mahendran等[33]利用16S rRNA對(duì)大蠶蛾科Satumiidae和蠶蛾科14個(gè)不同物種的系統(tǒng)發(fā)生分析也佐證了大蠶蛾科和蠶蛾科呈單系起源的觀點(diǎn)。2007年,Sobti等[34]以16S rRNA對(duì)6個(gè)北印度蝴蝶種的遺傳關(guān)系研究發(fā)現(xiàn),6個(gè)種按屬聚為2個(gè)不同的類(lèi)群,并指出結(jié)合rRNA與蛋白編碼基因?qū)⒛軓牟煌姆诸?lèi)階元上更好地理解鱗翅目昆蟲(chóng)的系統(tǒng)發(fā)生關(guān)系。另外,國(guó)內(nèi)也有利用16S rRNA對(duì)蛺蝶科、鳳蝶科內(nèi)亞科主要類(lèi)群進(jìn)化關(guān)系的研究報(bào)道[35-36]。

目前,已有部分研究將rRNA與線粒體蛋白編碼基因結(jié)合起來(lái)研究鱗翅目昆蟲(chóng)的進(jìn)化關(guān)系。Aubert等[37]研究了鳳蝶屬(Papilio)種間的系統(tǒng)發(fā)生關(guān)系,當(dāng)分別利用16S rRNA、ND1(僅密碼子第一二位堿基)以及16S rRNA+ND1(僅密碼子第一二位堿基)連接序列構(gòu)建了系統(tǒng)發(fā)生樹(shù),發(fā)現(xiàn)16S rRN A+ND1的系統(tǒng)發(fā)生樹(shù)最切合實(shí)際。另外,還有少量應(yīng)用16S rRNA+ND1連接序列在蝴蝶種屬分子進(jìn)化的研究報(bào)道[38-39]。Wahlberg和Zimmermann[40]測(cè)定了蛺蝶科網(wǎng)蛺蝶族Melitaeini 13屬77個(gè)物種的16S rRNA和COI基因的部分序列,系統(tǒng)發(fā)生樹(shù)顯示大多數(shù)前人所定義的屬都是并系的,但發(fā)現(xiàn)Euphydryas、Phyciodes、Chlosyne和Melitaea為單系群,進(jìn)化樹(shù)并暗示melitaeines起源于新北區(qū),然后曾三次擴(kuò)散到新熱帶區(qū)和兩次擴(kuò)散到古北區(qū)。2004年,Lange等[41]利用16S rRNA和COII的部分序列對(duì)夜蛾科Noctuidae和螟蛾總科Pyraloidea 10屬和6個(gè)族的26個(gè)種進(jìn)行分子進(jìn)化研究,結(jié)果發(fā)現(xiàn)夜蛾科為單系起源,而螟蛾總科呈并系發(fā)生;并且夜蛾科內(nèi)Sesamia和Bathytricha也是單系發(fā)生,而且Busseola是Bathytricha發(fā)生的基礎(chǔ)。

2.3 控制區(qū)

在過(guò)去的十年間,線粒體控制區(qū)已被廣泛應(yīng)用于動(dòng)物群體遺傳和分子進(jìn)化等研究。但是昆蟲(chóng)控制區(qū)是否能用于分子系統(tǒng)學(xué)等研究一直爭(zhēng)論不休[42]。2004年,Vila等[43]對(duì)控制區(qū)是否能用于鱗翅目昆蟲(chóng)分子進(jìn)化研究進(jìn)行了新的評(píng)價(jià)。他們對(duì)蛺蝶科的15個(gè)種進(jìn)行系統(tǒng)發(fā)生分析時(shí)發(fā)現(xiàn),控制區(qū)構(gòu)建的系統(tǒng)發(fā)生樹(shù)比用COI基因更能反映蛺蝶科內(nèi)真實(shí)的起源進(jìn)化關(guān)系。并進(jìn)一步指出,在解決快速進(jìn)化的種分類(lèi)水平的系統(tǒng)發(fā)生關(guān)系時(shí),CR+COI基因可能是一個(gè)非常有用的工具。2005年,Vila等[44]就以CR+COI基因序列對(duì)蛺蝶科Erebia triaria蝴蝶世系的初始分化與雜交進(jìn)行了研究,其結(jié)論支持在第四紀(jì)期間伊比利亞半島不是唯一的避難所,而是山脈與不同冰河時(shí)期的復(fù)雜分布促進(jìn)了群體的隔離。

3 展望

mtDNA已被廣泛應(yīng)用于鱗翅目昆蟲(chóng)的近緣種、屬、科甚至是目等分類(lèi)階元的系統(tǒng)發(fā)生分析和分子進(jìn)化等研究。但是不同的mt基因進(jìn)化速率不同,而且編碼基因密碼子的第三位堿基和控制區(qū)容易發(fā)生替換而達(dá)到飽和,加之外群選擇不當(dāng),常會(huì)引起系統(tǒng)發(fā)生樹(shù)長(zhǎng)枝吸收效應(yīng)直接影響分析的結(jié)果[45]。因此mt基因用于系統(tǒng)發(fā)育和分子進(jìn)化等研究也有其適合的領(lǐng)域,在不同的物種中可能選擇mt基因不同,在不同分類(lèi)階元也有各自適合的基因。通常親緣關(guān)系較遠(yuǎn)的物種,一般選取比較保守的基因;親緣關(guān)系較近的物種,則選取選擇進(jìn)化速率較快的序列?？v觀鱗翅目的研究,COI和COII基因應(yīng)用最廣,并且結(jié)合線粒體其他基因或區(qū)段(ND5、12S rRNA、16S rRNA及CR)或核基因(EF-1α、Wg)等進(jìn)行鱗翅目昆蟲(chóng)的分子系統(tǒng)學(xué)研究也有著不可替代的優(yōu)勢(shì)。相信mtDNA序列將在昆蟲(chóng)分子系統(tǒng)學(xué)研究中有著更為廣闊的空間,并隨著各物種進(jìn)化研究的相繼開(kāi)展必將為瀕危物種的保護(hù),多樣性的保持和種質(zhì)資源的保存帶來(lái)新的希望。

[1]Brower AV,DeSalle R.Patterns of mitochondrial versus nuclear DNA sequence divergence among nymphalid butterflies:the utility of wingless as a source of characters for phylogenetic inference[J].Insect Mol.Biol.,1998,7:73-82.

[2]Friedlander TP,Regier JC,Mitter C.Nuclear gene sequences for higher level phylogenetic analysis:fourteen promising candidates[J].Syst.Biol.,1992,41:483-490.

[3]Brown WM,Prager EM,Wang A,et al.Mitochondrial DNA sequences of primates:tempo and mode of evolution[J].J.Mol.Evol.,1982,18:225-239.

[4]DeSalle R,Freedman T,Prager EM,et al.Tempo and mode of sequence evolution in mitochondrial DNA of Hawaiian Drosophila[J].J.M ol.Evol.,1987,26:157-164.

[5]Pakendorf B,Stoneking M.Mitochondrial DNA and human evolution[J].Annu Rev Genomics Hum.Genet.,2005,6:165-183.

[6]魯成,廖順堯,劉運(yùn)強(qiáng),等.家蠶線粒體基因組全序列測(cè)定與分析[J].農(nóng)業(yè)生物技術(shù)學(xué)報(bào),2002,10:163-170.

[7]Yukuhiro K,Sezutsu H,Itoh M,et al.Significant levels of sequence divergence and gene rearrangements have occurred between the mitochondrial genomes of the wild mulberry silkmoth,Bombyx mandarina,and its close relative,the domesticated silkmoth,Bombyx mori[J].Mol.Biol.Evol.,2002,19:1385-1389.

[8]Hu XL,Cao GL,Xue RY,et al.The complete mitogenome and phylogenetic analysis of Bombyx mandarina strain Qingzhou[J].Mol.Biol.Rep.,2009(In press)

[9]Li D,Guo Y,Shao H,et al.Genetic diversity,molecular phylogeny and selection evidence of the silkworm mitochondria implicated by complete resequencing of 41 genomes[J].BMC Evol.Biol.,2010,10:81.

[10]Pan MH,Yu QY,Xia YL,et al.Characterization of mitochondrial genome of Chinese wild mulberry silkworm,Bomyx mandarina(Lepidoptera:Bombycidae)[J].Sci.China Ser.C-Life Sci.,2008,51:693-701.

[11]Yang L,Wei ZJ,Hong GY,et al.The complete nucleotide sequence of the mitochondrial genome of Phthonandria atrilineata(Lepidoptera:Geometridae)[J].Mol.Biol.Rep.,2008,36(6):1441-9

[12]DeChaine EG,Martini AP.Historic cycles of fragmentation and expansion in Parnassius smintheus(papilionidae)inferred using mitochondrial DNA[J].Evolution Int.J.Org.Evolution,2004,58:113-127.

[13]Rand DB,Heath A,Suderman T,et al.Phylogeny and life history evolution of the genus Chrysoritis within the Aphnaeini(Lepidoptera:Lycaenidae),inferred from mitochondrial cytochrome oxidase I sequences[J].Mol.Phylogenet.Evol.,2000,17:85-96.

[14]Hwang J,Lee JS,Goo TW,et al.The comparative molecular study between Bombycidae and Saturniidae based on mtDNA RFLP and cytochrome oxidase I gene sequences:implication for molecular evolution[J].Z.Naturforsch C,,1999,54:587-594.

[15]Behere GT,Tay WT,Russell DA,et al.Mitochondrial DNA analysis of field populations of Helicoverpa armigera(Lepidoptera:Noctuidae)and of its relationship to H.zea[J].BMC Evol.Biol.,2007,7:117.

[16]Long Y,Wan H,Yan F,et al.Glacial effects on sequence divergence of mitochondrial COII of Polyura eudamippus(Lepidoptera:Nymphalidae)in China[J].Biochem Genet,2006,44:361-77.

[17]Ohno S,Ishikawa Y,Tatsuki S,et al.Variation in mitochondrial COII gene sequences among two species of Japanese knotweed-boring moths,Ostrinia latipennis and O.ovalipennis(Lepidoptera:Crambidae)[J].Bull.Entomol.Res.,2006,96:243-9.

[18]王戎疆,萬(wàn)宏,龍玉,等.利用線粒體CO II基因序列對(duì)中國(guó)尾蛺蝶屬系統(tǒng)分化的研究(鱗翅目:蛺蝶科)[J].昆蟲(chóng)學(xué)報(bào),2004,47:243-247.

[19]Brown JM,Pellmyr O,Thompson JN,et al.Phylogeny of Greya(Lepidoptera:Prodoxidae),based on nucleotide sequence variation in mitochondrial cytochrome oxidase I and II:congruence with morphological data[J].Mol.Biol.Evol.,1994,11:128-141.

[20]Caterino MS,Sperling FA.Papilio phylogeny based on mitochondrial cytochrome oxidase I and II genes[J].Mol.Phylogenet.Evol.,1999,11:122-137.

[21]Zakharov EV,Caterino MS,Sperling FA.M olecular phylogeny,historical biogeography,and divergence time estimates for swallowtail butterflies of the genus Papilio(Lepidoptera:Papilionidae)[J].Syst.Biol.,2004,53:193-215.

[22]Beckenbach AT,Wei YW,Liu H.Relationships in the Drosophila obscura species group inferred from mitochondrial cytochrome oxidase II sequences[J].Mol.Biol.Evol,1993,10:619-634.

[23]M onteiro A,Pierce NE.Phylogeny of Bicyclus(Lepidoptera:Nymphalidae)inferred from COI,COII,and EF-1alpha gene sequences[J].Mol.Phylogenet.Evol.,2001,18:264-281.

[24]Clary DO,Wolstenholme DR The mitochondrial DNA molecule of Drosophila yakuba:nucleotide sequence,gene organization,and genetic code J M ol Evol.,1985,22:252-271.

[25]Yagi T,Sasaki G,Takebe H.Phylogeny of Japanese papilionid butterflies inferred from nucleotide sequences of the mitochondrial ND5 gene[J].J.Mol.Evol.,1999,48:42-48.

[26]Yagi T,Katoh T,Chichvarkhin A,et al.M olecular phylogeny of butterflies Parnassius glacialis and P.stubbendorfii at various localities in East Asia[J].Genes Genet.Syst.,2001,76:229-234.

[27]Omoto K,Katoh T,Chichvarkhin A,et al.M olecular systematics and evolution of the"Apollo"butterflies of the genus Parnassius(Lepidoptera:Papilionidae)based on mitochondrial DNA sequence data[J].Gene,2004,326:141-147.

[28]Albre J,Gers C,Legal L.Molecular phylogeny of the Erebia tyndarus(Lepidoptera,Rhopalocera,Nymphalidae,Satyrinae)species group combining CoxII and ND5 mitochondrial genes:a case study of a recent radiation[J].Mol.Phylogenet.Evol.,2008,47:196-210.

[29]Li A,Zhao Q,Tang S,et al.Molecular phylogeny of the domesticated silkworm,Bombyx mori,based on the sequences of mitochondrial cytochrome b genes[J].J.Genet.,2005,84:137-142.

[30]陳永久,張亞平,沈發(fā)榮,等.中國(guó)5種珍稀絹蝶非損傷性取樣mtDNA序列及系統(tǒng)進(jìn)化[J].遺傳學(xué)報(bào),1999,26:203-207.

[31]Hwang JS,Lee JS,Goo TW,et al.Molecular genetic relationships between Bombycidae and Saturniidae based on the mitochondria DNA encoding of large and small rRNA[J].Genet.Anal.,1999,15:223-8.

[32]劉彥群,靳向東,秦利,等.九種絹絲昆蟲(chóng)線粒體12S rRNA基因的序列特征和系統(tǒng)發(fā)育分析[J].昆蟲(chóng)學(xué)報(bào),2008,51:307-314.

[33]Mahendran B,Ghosh SK,Kundu SC.Molecular phylogeny of silk-producing insects based on 16S ribosomal RNA and cytochrome oxidase subunit I genes[J].J.Genet.,2006,85:31-38.

[34]Sobti RC,Sharma VL,Kumari M,et al.Genetic relatedness of six North-Indian butterfly species(Lepidoptera:Pieridae)based on 16S rRNA sequence analysis[J].Mol.Cell.Biochem.,2007,295:145-51.

[35]陳娜,朱國(guó)萍,郝家勝,等.基于線粒體rDNA序列探討蛺蝶科(鱗翅目,蝶亞目)主要分類(lèi)群的系統(tǒng)發(fā)生關(guān)系[J].動(dòng)物學(xué)報(bào),2007,53:106-115.

[36]蘇成勇,朱國(guó)萍,郝家勝,等.鳳蝶亞科(鳳蝶科,鱗翅目)16S rRNA基因的分子系統(tǒng)發(fā)生分析[J].動(dòng)物分類(lèi)學(xué)報(bào),2007,32:335-342.

[37]Aubert J,Legal L,Descimon H,et al.Molecular phylogeny of swallowtail butterflies of the tribe Papilionini(Papilionidae,Lepidoptera)[J].Mol.Phylogenet.Evol.,1999,12:156-167.

[38]Katoh T,Chichvarkhin A,Yagi T,et al.Phylogeny and evolution of butterflies of the genus Parnassius:inferences from mitochondrial 16S and ND1 sequences[J].Zoolog.Sci.,2005,22:343-51.

[39]M artin J,Gilles A,Descimon H.Molecular phylogeny and evolutionary patterns of the european satyrids(Lepidoptera:Satyridae)as revealed by mitochondrial gene sequences[J].Mol.Phylogenet.Evol.,2000,15:70-82.

[40]Wahlberg N,Zimmermann M.Pattern of Phylogenetic Relationships among Members of the Tribe M elitaeini(Lepidoptera:Nymphalidae)Inferred from Mitochondrial DNA Sequences[J].Cladistics,2000,16:347-363.

[41]Lange CL,Scott KD,Graham GC,et al.Sugarcane moth borers(Lepidoptera:Noctuidae and Pyraloidea):phylogenetics constructedusing COII and 16S mitochondrial partial gene sequences[J].Bull.Entomol.Res.,2004,94:457-64.

[42]Zhang DX,Hewitt GM.Insect mitochondrial control region:A review of its structure,evolution,and usefulness in evolutionary studies[J].Biochem.Syst.Ecol.,1997,25:99-120.

[43]Vila M,Bjorklund M.The utility of the neglected mitochondrial control region for evolutionary studies in lepidoptera(insecta)[J].J.Mol.Evol.,2004,58:280-290.

[44]Vila M,Vidal-Romani JR,Bjorklund M.The importance of time scale and multiple refugia:incipient speciation and admixture of lineages in the butterfly Erebia triaria(Nymphalidae)[J].Mol.Phylogenet.Evol.,2005,36:249-60.

[45]Wheeler WC.Nucleic acid sequence phylogeny and random outgroups[J].Cladistics,1990,6:363-367.