,,, ,
庫(kù)蠓(CulicoidesLatreille)俗稱“墨蚊”、“小咬”,隸屬于昆蟲(chóng)綱(insecta)、雙翅目(diptera)、蠓科(ceratopogonidae),全世界現(xiàn)有1 368種、中國(guó)已知348種[1-2],是引起野生和家養(yǎng)反芻動(dòng)物流行病的重要媒介昆蟲(chóng),其通過(guò)刺吸動(dòng)物血液來(lái)傳播藍(lán)舌病(BT)、家畜流行病出血熱(EHD)、非洲馬瘟病(AHS)等諸多動(dòng)物傳染病[3]。上述疾病的病原通過(guò)媒介庫(kù)蠓和易感宿主之間的一系列復(fù)制循環(huán)而持續(xù)存活于自然界中[4],因而其易于引發(fā)醫(yī)療衛(wèi)生和獸醫(yī)等領(lǐng)域的家畜流行病,從而造成牛、羊等養(yǎng)殖業(yè)的巨大經(jīng)濟(jì)損失,歐洲曾發(fā)生過(guò)藍(lán)舌病大流行[5-6],國(guó)內(nèi)云南、廣西、新疆和內(nèi)蒙古等29個(gè)省也曾先后在牛羊血液中檢測(cè)出藍(lán)舌病病毒(BTV)[7]。目前,BT和AHS已被世界動(dòng)物衛(wèi)生組織(OIE)列為A類(lèi)傳染病[8],中國(guó)也將BT列為一類(lèi)動(dòng)物疫病[9]。
庫(kù)蠓體型微小、鑒別特征稀少,傳統(tǒng)分類(lèi)鑒定尚存在諸多問(wèn)題。首先是非成蟲(chóng)階段的蟲(chóng)態(tài)和蟲(chóng)體殘缺的成蟲(chóng)均難以鑒定,其次受庫(kù)蠓性別限制,雄蟲(chóng)借助外生殖器等形態(tài)特征易于鑒別,而雌蟲(chóng)缺乏良好的尾器特征而不易鑒別,最后庫(kù)蠓表型可塑性和遺傳可變性等也容易導(dǎo)致誤判[10-12]。因此,庫(kù)蠓分類(lèi)中亟需快速、準(zhǔn)確鑒定種類(lèi)的技術(shù)方法,對(duì)蠓傳疾病的監(jiān)控、防治具有重要意義。
隨著分子生物學(xué)和生物信息學(xué)的快速發(fā)展,分子生物學(xué)技術(shù)已應(yīng)用于昆蟲(chóng)分類(lèi)鑒定和系統(tǒng)演化等領(lǐng)域,這在很大程度上彌補(bǔ)了傳統(tǒng)形態(tài)分類(lèi)的不足,尤其是近似種或隱種的分類(lèi)[13-15]。在昆蟲(chóng)核基因中,串聯(lián)重復(fù)是核糖體DNA(rDNA)的存在形式,同時(shí)其一個(gè)轉(zhuǎn)錄單元包括18S、5.8S、28S三個(gè)部分,而ITS1和ITS2兩個(gè)可變區(qū)則分別位于18S和5.8S以及5.8S 和28S之間[16]。目前,昆蟲(chóng)分子系統(tǒng)學(xué)研究也主要集中在rDNA的ITS、28S rRNA和18S rRNA等區(qū)域。高度保守的序列通常應(yīng)用于高級(jí)分類(lèi)階元的系統(tǒng)發(fā)育研究,如28S rRNA和18S rRNA[17];而變異較大、選擇壓力小的序列則適用于種屬水平上的研究,如ITS1和ITS2序列[18],且已廣泛應(yīng)用于媒介昆蟲(chóng)的種類(lèi)鑒定[19-20]。然而,蠓科昆蟲(chóng)的分子系統(tǒng)學(xué)研究,自1992年Tabachnick等[21]報(bào)道變翅庫(kù)蠓(Culicoidesvariipennis)種群的基因差異開(kāi)始,之后國(guó)外關(guān)于蠓科分子系統(tǒng)學(xué)的研究與日劇增[22-26],但國(guó)內(nèi)相關(guān)研究嚴(yán)重不足。
因此,本研究以埃蠓屬的環(huán)紋埃蠓Alloheleaannulata為外群,對(duì)庫(kù)蠓屬中的凹緣庫(kù)蠓C.holcus、連斑庫(kù)蠓C.jacobsoni、印度庫(kù)蠓C.indianus、荒川庫(kù)蠓C.arakawai、霍飛庫(kù)蠓C.huffi和肩宏庫(kù)蠓C.humeralis等6種庫(kù)蠓進(jìn)行分類(lèi)研究(凹緣庫(kù)蠓C.holcus、和連斑庫(kù)蠓C.jacobsoni為二囊亞屬近似種),通過(guò)聯(lián)合形態(tài)學(xué)特征和核糖體DNA-ITS1序列數(shù)據(jù),彌補(bǔ)傳統(tǒng)形態(tài)鑒定的不足,進(jìn)一步證實(shí)ITS1序列在庫(kù)蠓及其近似種的分類(lèi)鑒別方面的適用性,為解決蠓科昆蟲(chóng)分類(lèi)中形態(tài)學(xué)數(shù)據(jù)難以解決的問(wèn)題提供更多的數(shù)據(jù)支撐。
1.1.1庫(kù)蠓標(biāo)本 本研究共使用7種23只蠓蟲(chóng)標(biāo)本的數(shù)據(jù),其中18只為本研究組2015年7月采自廣西大瑤山自然保護(hù)區(qū),其余5只的數(shù)據(jù)來(lái)自GenBank數(shù)據(jù)庫(kù),詳見(jiàn)表1。
1.1.2主要試劑及儀器 基因組提取試劑盒(QIAGEN, 德國(guó))、pEASY-T1克隆試劑盒(北京全式金公司, 中國(guó))、PCR MasterMix(TIANGEN, 中國(guó))、引物(上海捷瑞生物工程有限公司)、瓊脂糖凝膠回收試劑盒(TIANGEN, 中國(guó))、瓊脂糖(Biowest Agarose, 西班牙)、PCR擴(kuò)增儀(美國(guó)BOI-RAD公司)、UVP凝膠成像系統(tǒng)(美國(guó)BOI-RAD公司)等。
1.2.1庫(kù)蠓采集 采用網(wǎng)捕法和燈誘法,在廣西大瑤山自然保護(hù)區(qū)不同生境設(shè)置采樣點(diǎn)進(jìn)行蠓蟲(chóng)收集,并將其保存于無(wú)水乙醇中。
1.2.2標(biāo)本鑒定 按照常規(guī)方法進(jìn)行標(biāo)本的分選、解剖、制片、鑒定以及拍照,且胸部標(biāo)本單獨(dú)編號(hào)保存在無(wú)水乙醇中。
1.2.3基因組DNA 將庫(kù)蠓標(biāo)本研磨至無(wú)明顯組織塊,按照QIAGEN DNeasy Blood& Tissue Kit說(shuō)明書(shū)進(jìn)行其基因組DNA提取,并保存至-20 ℃?zhèn)溆谩?/p>
1.2.4PCR擴(kuò)增 PCR擴(kuò)增的上游引物為PanCulF (5′-GTAGGTGAACCTGCGGAAGG-3′),下游引物為PanCulR (5′-TGCGGTCTTCATCGACCCAT-3′),由上海捷瑞生物工程有限公司合成。PCR擴(kuò)增的體系如表2所示,其反應(yīng)條件為94 ℃預(yù)變性3 min;94 ℃ 30 s,40 ℃ 60 s,68 ℃ 60 s,循環(huán)5次;94 ℃ 60 s,44 ℃ 60 s,68 ℃ 60 s,循環(huán)35次;最后68 ℃延伸10 min。
表1 供試庫(kù)蠓的詳細(xì)信息Tab.1 Detailed information of the biting midges in this study
表2 PCR擴(kuò)增體系Tab.2 PCR amplification system
1.2.5PCR擴(kuò)增產(chǎn)物純化、克隆及測(cè)序 將PCR擴(kuò)增產(chǎn)物進(jìn)行瓊脂糖凝膠電泳,同時(shí)利用DNA回收試劑盒進(jìn)行DNA產(chǎn)物的純化,之后使用pEASY-T1克隆試劑盒進(jìn)行DNA克隆,最后用60%甘油保存陽(yáng)性克隆菌液并送至上海生物工程股份有限公司進(jìn)行雙向測(cè)序。
1.2.6rDNA-ITS1序列分析 利用Clustal X、DNA MAN、MEGA 6.06等軟件對(duì)本研究中rDNA-ITS1序列進(jìn)行比對(duì)、編輯、比較等分析,統(tǒng)計(jì)序列中各堿基(A、T、C、G)的組成比例以及保守位點(diǎn)(Conserved sites, C)、變異位點(diǎn)(Variable sites, V)、簡(jiǎn)約信息位點(diǎn)(Parsim-Imformative sites, Pi)和自裔位點(diǎn)(Singleton sites, S)的數(shù)量。同時(shí),以Kimura雙參數(shù)模型計(jì)算遺傳距離,并以環(huán)紋埃蠓A.annulata為外群,通過(guò)鄰接法(NJ)、最大似然法(ML)構(gòu)建分子系統(tǒng)樹(shù),并通過(guò)Bootstrap 1 000次自舉檢驗(yàn)系統(tǒng)樹(shù)中各結(jié)點(diǎn)的置信值。
1.2.7統(tǒng)計(jì)學(xué)方法 采用SPSS 21.0軟件進(jìn)行統(tǒng)計(jì)分析,數(shù)據(jù)以均數(shù)±標(biāo)準(zhǔn)差表示。兩兩比較采用兩獨(dú)立樣本t檢驗(yàn),P<0.05表示差異有統(tǒng)計(jì)學(xué)意義。
2.1形態(tài)學(xué)特征對(duì)比 本研究?jī)?nèi)群的6種庫(kù)蠓形態(tài)上比較接近,尤其是凹緣庫(kù)蠓C.holcus和連斑庫(kù)蠓C.jacobsoni最為近似,以凹緣庫(kù)蠓C.holcus、連斑庫(kù)蠓C.jacobsoni和荒川庫(kù)蠓C.arakawai為例(表3、圖1)進(jìn)行形態(tài)特征比較。上述3種庫(kù)蠓中雄蟲(chóng)可根據(jù)其尾器的4個(gè)特征進(jìn)行有效鑒別,但雌蟲(chóng)中僅荒川庫(kù)蠓C.arakawai可直接根據(jù)頭部、翅部和腹部的明顯差異進(jìn)行準(zhǔn)確區(qū)分,另外二囊亞屬的兩近似種除翅部特征的細(xì)微差別外,余未見(jiàn)差異,故鑒定非常困難。
表3 本研究三種庫(kù)蠓形態(tài)特征對(duì)比Tab.3 Comparison of morphological features of three Culicoides
1:凹緣庫(kù)蠓C.holcus;2:連斑庫(kù)蠓C.jacobsoni;3:印度庫(kù)蠓C.indianus;4:荒川庫(kù)蠓C.arakawai(雌);5:荒川庫(kù)蠓C.arakawai(雄);6:環(huán)斑埃蠓A.annulata;7:霍飛庫(kù)蠓C.huffi;8:肩宏庫(kù)蠓C.humeralis圖1 庫(kù)蠓近似種一側(cè)翅圖Fig.1 Wing of one side of sibling species
2.2DNA序列分析 本研究獲得的DNA序列經(jīng)Clustal W比對(duì)和人工校對(duì)后,得到rDNA-ITS1的序列長(zhǎng)度為352 bp,該序列中存在較多的插入和缺失位點(diǎn)。通過(guò)MEGA 6.06軟件分析ITS1序列中T、C、A、G四種堿基的平均含量分別為36.8%、13.0%、30.0%和20.1%,其中A+T平均含量約為66.8%,C+G平均含量約為33.1%;而在3個(gè)密碼子位點(diǎn)中,第1位點(diǎn)的A+T含量相對(duì)最高(67.7%),第3位點(diǎn)A+T平均含量相對(duì)最低(64.6%)。此外,ITS1序列中保守位點(diǎn)(C)、變異位點(diǎn)(V)、簡(jiǎn)約信息位點(diǎn)(Pi)和自裔位點(diǎn)(S)的個(gè)數(shù)分別為163、189、182和7,其中變異位點(diǎn)(V)占比約為53.7%。
2.3序列同源性比對(duì) 將本研究所獲全部庫(kù)蠓ITS1序列與GenBank數(shù)據(jù)庫(kù)中的序列進(jìn)行同源比對(duì),搜索到部分同源序列,即荒川庫(kù)蠓C.arakawai(登錄號(hào):AB462265、AJ489503,序列相似度:99%)、連斑庫(kù)蠓C.jacobsoni(登錄號(hào):AB462262,序列相似度:97%)、印度庫(kù)蠓C.indianus(登錄號(hào)為AB462268,序列相似度:93%)。
2.4遺傳距離 利用MEGA 6.06軟件計(jì)算本研究中蠓蟲(chóng)ITS1序列的遺傳距離,種內(nèi)和種間的遺傳距離范圍分別為0.000~0.020和0.113~0.621,通過(guò)統(tǒng)計(jì)軟件分析得出種內(nèi)組和種間組的遺傳距離(均數(shù)±標(biāo)準(zhǔn)差)分別為0.009±0.008和0.391±0.169,兩者差異具有統(tǒng)計(jì)學(xué)意義(t=32.430,P<0.05)。
2.5系統(tǒng)發(fā)育分析 以環(huán)紋埃蠓A.annulata為外群,采用鄰接法(NJ)和最大似然法(ML)分別構(gòu)建6種庫(kù)蠓的系統(tǒng)發(fā)育樹(shù),得到基本一致的拓?fù)浣Y(jié)構(gòu)(圖2、3)。系統(tǒng)發(fā)育樹(shù)聚類(lèi)分析表明:1)各物種間分別構(gòu)成單系(群),同種不同地理種群聚為一支,且種間無(wú)交叉。2)物種間親緣關(guān)系明顯,二囊亞屬的凹緣庫(kù)蠓C.holcus和連斑庫(kù)蠓C.jacobsoni兩種首先聚在一起,親緣關(guān)系較近并構(gòu)成一個(gè)姊妹群,再與庫(kù)蠓亞屬的印度庫(kù)蠓C.indianus并為一大支,同時(shí)帶紋亞屬的荒川庫(kù)蠓C.arakawai與三囊亞屬的肩宏庫(kù)蠓C.humeralis親緣關(guān)系較近先聚為一小支,再與屋室亞屬的霍飛庫(kù)蠓C.huffi結(jié)合構(gòu)成一大支,之后兩大支聚合形成內(nèi)群,而外群埃蠓屬的環(huán)紋埃蠓A.annulata在樹(shù)的根部,外群與內(nèi)群形成完整的有根樹(shù)。
各分支上的數(shù)字為自舉檢驗(yàn)置信值1 000次圖2 基于rDNA-ITS1序列構(gòu)建的NJ樹(shù)Fig.2 Neighbor-Joining tree constructed based on rDNA-ITS1 sequence
各分支上的數(shù)字為自舉檢驗(yàn)置信值1 000次圖3 基于rDNA-ITS1序列構(gòu)建的ML樹(shù)Fig.3 Maximum Likelihood tree constructed based on rDNA-ITS1 sequence
核糖體DNA內(nèi)轉(zhuǎn)錄間隔區(qū)1 (nuclear ribosomal DNA internal transcribed spacer 1, rDNA-ITS1)基因位于核糖體DNA的18S和5.8S之間,其特點(diǎn)為:①屬于中度重復(fù)序列;②不承擔(dān)編碼蛋白質(zhì)的功能,也不參與核糖體的形成,故承受的選擇壓力較小,進(jìn)化速度相對(duì)較快,從而具有高變性;③在ITS1序列兩側(cè)的18S和5.8S序列具有保守性,便于通用引物的設(shè)計(jì)、PCR擴(kuò)增和測(cè)序;④序列的組成、長(zhǎng)度以及折疊方式存在多態(tài)性,可用于屬、種和種群水平的研究[27-28]。目前,應(yīng)用ITS1解決親緣關(guān)系較近的庫(kù)蠓分類(lèi)以及系統(tǒng)關(guān)系的研究較多,如Mathieu等[29]通過(guò)設(shè)計(jì)ITS1-5.8S-ITS2 的TaqMan 探針對(duì)不顯庫(kù)蠓(Culicoidesobsoletus)和蘇格蘭庫(kù)蠓(Culicoidesscoticus)進(jìn)行雙重實(shí)時(shí)熒光定量PCR實(shí)驗(yàn),Deblauwe等[30]通過(guò)設(shè)計(jì)ITS1的DNA微陣列形式的探針雜交對(duì)不顯庫(kù)蠓種群(Obsoletusgroup)的庫(kù)蠓進(jìn)行分子鑒定,說(shuō)明ITS1是庫(kù)蠓種間和種內(nèi)遺傳分化程度及系統(tǒng)進(jìn)化研究的良好指標(biāo)。因此,基于庫(kù)蠓形態(tài)特征和rDNA-ITS1序列數(shù)據(jù)的基礎(chǔ)上對(duì)6種吸血蠓進(jìn)行鑒定和分析,其研究結(jié)果表明形態(tài)學(xué)特征與分子數(shù)據(jù)具有一致性,共同驗(yàn)證了rDNA-ITS1基因序列在吸血庫(kù)蠓(包括近似、近緣種)鑒別和系統(tǒng)發(fā)育研究中的適用性。
本研究rDNA-ITS1序列的堿基組成中A+T的平均含量(76.8%)遠(yuǎn)高于C+G的平均含量(33.2%),具有AT偏向顯著。同時(shí),該序列的變異位點(diǎn)(189個(gè))較多,占比約為53.7%,且相對(duì)于線粒體基因(如COI和Cytb)其還存在較多的插入和缺失位點(diǎn)。Cêtre-Sossah等[31]利用ITS1序列對(duì)盛行于非洲和歐洲的藍(lán)舌病及非洲馬瘟病的主要傳播媒介庫(kù)蠓——?dú)埑釒?kù)蠓C.imicola進(jìn)行分子鑒定,通過(guò)設(shè)計(jì)特異引物可得到長(zhǎng)度各不相同的ITS1擴(kuò)增產(chǎn)物,根據(jù)電泳條帶的位置即可進(jìn)行種類(lèi)鑒別;同時(shí),Gomulski等[32]利用rDNA-ITS2序列探究庫(kù)蠓亞屬及近緣種的系統(tǒng)發(fā)生,研究發(fā)現(xiàn)其存在明顯的多態(tài)性,表明絕大多數(shù)庫(kù)蠓的ITS序列存在多態(tài)性,而并非是由于PCR擴(kuò)增過(guò)程中產(chǎn)生的錯(cuò)誤所致。
進(jìn)一步對(duì)庫(kù)蠓遺傳距離研究發(fā)現(xiàn),各種間遺傳距離(0.113~0.621)明顯大于各種內(nèi)遺傳距離(0.000~0.020),其中最大種內(nèi)遺傳距離(印度庫(kù)蠓C.indianus:0.020)亦小于最小種間遺傳距離(凹緣庫(kù)蠓C.holcus與連斑庫(kù)蠓C.jacobsoni:0.113),進(jìn)一步說(shuō)明庫(kù)蠓種間和種內(nèi)遺傳距離不存在重疊區(qū),足以進(jìn)行準(zhǔn)確區(qū)分。同時(shí),種間平均遺傳距離(0.411)為種內(nèi)平均遺傳距離(0.010)的41倍,也完全符合Hebert 等[33-34]提出的物種鑒定規(guī)則。因此,rDNA-ITS1基因可應(yīng)用于庫(kù)蠓的分子分類(lèi)。然而,本研究中印度庫(kù)蠓C.indianus的種內(nèi)遺傳距離相對(duì)其他種類(lèi)略大(0.020),Matsumoto等[35]在基于ITS1和ITS2對(duì)25種庫(kù)蠓進(jìn)行研究的過(guò)程中發(fā)現(xiàn)印度庫(kù)蠓C.indianus存在ITS1長(zhǎng)和短兩種序列類(lèi)型,這可能就是造成本研究中其種內(nèi)遺傳距離偏大的原因。另外,關(guān)于昆蟲(chóng)近緣種雜交的問(wèn)題,Pesson等{36]在對(duì)白蛉亞科昆蟲(chóng)進(jìn)行研究時(shí)發(fā)現(xiàn)其確實(shí)存在近緣物種雜交,但本研究中采自統(tǒng)一地點(diǎn)的凹緣庫(kù)蠓C.holcus和連斑庫(kù)蠓C.jacobsoni并未發(fā)現(xiàn)任何基因滲透現(xiàn)象,這與Ritchie等[37]和Augot等[15]關(guān)于庫(kù)蠓的研究結(jié)果相一致。
此外,在對(duì)6種庫(kù)蠓的DNA序列聚類(lèi)分析研究中發(fā)現(xiàn),其與形態(tài)學(xué)鑒定和遺傳距離計(jì)算的結(jié)果基本相同,如凹緣庫(kù)蠓C.holcus和連斑庫(kù)蠓C.jacobsoni兩近似種同屬二囊亞屬,其形態(tài)學(xué)特征相似且種間遺傳距離(0.113)最小,在系統(tǒng)發(fā)育樹(shù)中優(yōu)先聚在一起構(gòu)成姊妹群。然而,來(lái)自日本的連斑庫(kù)蠓C.jacobsoni(AB462262.1)與我國(guó)廣西的種類(lèi)首先分開(kāi),從拓?fù)浣Y(jié)構(gòu)上來(lái)看我國(guó)廣西的凹緣庫(kù)蠓C.holcus和連斑庫(kù)蠓C.jacobsoni聚為一支,但是從時(shí)間節(jié)點(diǎn)和分支長(zhǎng)短來(lái)看同屬連斑庫(kù)蠓C.jacobsoni的種類(lèi)具有更近的親緣關(guān)系。此外,構(gòu)建的NJ樹(shù)和ML樹(shù)的拓?fù)浣Y(jié)構(gòu)也基本一致,在種間和種內(nèi)的分支上均具較高的置信值,說(shuō)明應(yīng)用rDNA-ITS1基因構(gòu)建的分子進(jìn)化樹(shù)能明顯區(qū)分庫(kù)蠓及其近似種,可靠程度較高。然而,荒川庫(kù)蠓C.arakawai和肩宏庫(kù)蠓C.humeralis的種內(nèi)存在置信值高低不等的差異,這與Gomulski 等[38]在基于同為核糖體DNA-ITS區(qū)的 ITS2序列對(duì)包括不顯庫(kù)蠓Culicoidesobsoletus在內(nèi)的二囊亞屬進(jìn)行系統(tǒng)發(fā)育關(guān)系的研究過(guò)程中的情況類(lèi)似,究其原因可能主要是兩種建樹(shù)方法之間存在建樹(shù)理論基礎(chǔ)和建樹(shù)偏向性方面的差異。與此同時(shí),在種內(nèi)個(gè)別分支上亦存在相對(duì)較低的置信值,如凹緣庫(kù)蠓C.holcus和連斑庫(kù)蠓C.jacobsoni之間的置信值僅為51,這與李瑋瑋[39]等對(duì)小蔗螟的研究以及秦芳[40]等對(duì)利川馬遺傳多態(tài)性的研究得到的結(jié)果相似,其可能是由于ITS1序列在堿基組成及變異速度方面存在差異以及ITS1序列片段短(352 bp,且存在插入和缺失位點(diǎn))和樣本數(shù)量少(n=23)等導(dǎo)致提供的有效分子遺傳學(xué)信息有限有關(guān)。
在廣泛使用ITS1等核糖體內(nèi)轉(zhuǎn)錄間隔區(qū)的基因序列進(jìn)行分析時(shí),發(fā)現(xiàn)它存在基因的橫向轉(zhuǎn)移、重復(fù)以及外顯子拼接、域拼接等非線性的現(xiàn)象,說(shuō)明其并非對(duì)所有物種都具有適用性[28]。同時(shí),在核基因組研究中發(fā)現(xiàn)其存在線粒體假基因(NUMTs)的插入現(xiàn)象,這毋庸置疑是影響核基因的分子鑒定過(guò)程中準(zhǔn)確性的又一不利因素[41]。目前,蠓科學(xué)者在進(jìn)行庫(kù)蠓的系統(tǒng)發(fā)育分析時(shí),不僅僅是從傳統(tǒng)形態(tài)學(xué)角度進(jìn)行單因素的分析,同時(shí)還聯(lián)合線粒體基因、核基因,甚至是三維形態(tài)重建等多方面進(jìn)行綜合比對(duì)分析,如Augot等[15]聯(lián)合Cytb、CO I、ITS1和形態(tài)學(xué)特征可進(jìn)行翅痣庫(kù)蠓(Culicoidesstigma)和帕羅庫(kù)蠓(Culicoidesparroti)等近似種的準(zhǔn)確鑒別。同時(shí),分類(lèi)學(xué)家們也在不斷發(fā)掘新的分子標(biāo)記和分類(lèi)技術(shù),如CAD基因和幾何形態(tài)學(xué)等在庫(kù)蠓的分子系統(tǒng)學(xué)研究中的應(yīng)用[42-43]。因此,在后續(xù)的昆蟲(chóng)分子系統(tǒng)學(xué)研究領(lǐng)域中的趨勢(shì)將是結(jié)合生態(tài)學(xué)、形態(tài)學(xué)特征和多個(gè)分子標(biāo)記數(shù)據(jù)等多重分類(lèi)變量以及開(kāi)發(fā)新分子標(biāo)記和分類(lèi)技術(shù)來(lái)進(jìn)行更為全面、系統(tǒng)的分析。