楊美瓊， 張建慶， 蔡怡珊， 王飛鵬， 黃恩炯

(1.福州國(guó)際旅行衛(wèi)生保健中心，福建 福州350001；2.福建農(nóng)林大學(xué)植物保護(hù)學(xué)院，福建 福州350002)

庫(kù)蠓屬(Culicoides)屬節(jié)肢動(dòng)物門(mén)(Arthropoda)昆蟲(chóng)綱(Insecta)雙翅目(Diptera)蠓科(Ceratopogonidae)，是蠓科中最大的一個(gè)屬，也是分布最廣、種類(lèi)最多、與人畜關(guān)系最密切的一類(lèi)吸血蠓，還是多種人獸共患病的傳播媒介[1-2].庫(kù)蠓的種類(lèi)鑒定是動(dòng)植物檢疫和疾病預(yù)防控制的核心環(huán)節(jié).

目前，蠓類(lèi)的鑒定仍主要依賴于傳統(tǒng)形態(tài)學(xué)方法.但庫(kù)蠓不僅體型微小、鑒別特征少，而且具有表型可塑性和遺傳可變性，同時(shí)鑒定常受蟲(chóng)態(tài)、性別和標(biāo)本完整性的影響，所以容易出現(xiàn)誤判[3].為彌補(bǔ)傳統(tǒng)形態(tài)學(xué)鑒定的不足，亟需尋找一種快速準(zhǔn)確的鑒定技術(shù).

近年來(lái)，隨著以PCR 為基礎(chǔ)的DNA 序列測(cè)定方法的建立和廣泛使用，以DNA 條形碼技術(shù)為主的分子分類(lèi)技術(shù)在昆蟲(chóng)的鑒定中發(fā)揮了重要作用，極大地彌補(bǔ)了傳統(tǒng)形態(tài)學(xué)鑒定的不足.線粒體COⅠ、COⅡ和核糖體DNA(rDNA)的內(nèi)轉(zhuǎn)錄間隔區(qū)ITS 等基因也在蠓類(lèi)的分子鑒定中得到應(yīng)用，澄清了蠓類(lèi)一些近緣種和復(fù)合組的分類(lèi)問(wèn)題[4].然而，其主要以單基因檢測(cè)為主，仍難以實(shí)現(xiàn)對(duì)近緣種和復(fù)合組的準(zhǔn)確鑒定.

本研究選擇COⅠ、Cyt b、16S rDNA 和ITS 等多基因?qū)Σ勺晕覈?guó)14 個(gè)省(自治區(qū))的77 種庫(kù)蠓進(jìn)行分子鑒定，以期為庫(kù)蠓的準(zhǔn)確鑒定提供依據(jù).

1 材料與方法

1.1 材料

1.1.1 試驗(yàn)標(biāo)本 庫(kù)蠓采自我國(guó)福建、海南、廣西、黑龍江、河北、安徽、湖南、新疆、青海、甘肅、四川、云南、西藏及臺(tái)灣等14 個(gè)省(自治區(qū))的36 個(gè)采樣點(diǎn)，共77 種284 份標(biāo)本.將每份標(biāo)本拍照并編號(hào)后用于試驗(yàn).所有標(biāo)本均保存于軍事科學(xué)院軍事醫(yī)學(xué)研究院醫(yī)學(xué)昆蟲(chóng)標(biāo)本館.

1.1.2 主要試劑 昆蟲(chóng)基因組DNA 提取試劑盒(Insect gDNA Minipre Kit)購(gòu)自Biomiga 公司，凝膠回收試劑盒(E.Z.N.A.Gel Extraction Kit)購(gòu)自O(shè)MEGA 公司，引物由Invitrogen 公司合成，克隆試劑盒(pMD 18-T Vector Cloning Kit)購(gòu)自TaKaRa 公司.

1.2 方法

1.2.1 標(biāo)本采集 采用燈誘法和牛帳誘法[5]采集庫(kù)蠓.采獲的庫(kù)蠓經(jīng)麻醉后保存于95%乙醇中.

1.2.2 標(biāo)本制作和形態(tài)鑒定 參照文獻(xiàn)[5]進(jìn)行標(biāo)本制作和形態(tài)鑒定.鑒定結(jié)果由軍事科學(xué)院軍事醫(yī)學(xué)研究院虞以新研究員復(fù)核確認(rèn).

1.2.3 基因組DNA 提取 取庫(kù)蠓翅基部約1 mm3的組織，研碎，按照DNA 提取試劑盒說(shuō)明書(shū)提取基因組DNA，-20 ℃保存?zhèn)溆?

1.2.4 PCR 擴(kuò)增 選用COⅠ、Cyt b、16S rDNA 和ITS 基因進(jìn)行PCR 擴(kuò)增，引物序列及相關(guān)信息見(jiàn)表1.除用于切膠回收的PCR 反應(yīng)體系為50 μL 外，其余均為25 μL 擴(kuò)增體系.PCR 擴(kuò)增體系(25 μL):DNA 模板1 μL，上、下游引物各0.5 μL，10×緩沖液2.5 μL，dNTPs 0.5 μL，Taq 酶1.25 μL，ddH2O 18.75 μL.反應(yīng)條件:94 ℃預(yù)變性3 min，94 ℃45 s，40～48 ℃45 s，72 ℃60 s，35 個(gè)循環(huán)，72 ℃延伸10 min.

表1 引物信息Table 1 Information on the PCR primers

1.2.5 PCR 產(chǎn)物純化、克隆及測(cè)序 PCR 產(chǎn)物經(jīng)凝膠電泳檢測(cè)后，用凝膠回收試劑盒進(jìn)行純化，需要克隆的樣品使用克隆試劑盒進(jìn)行DNA 克隆，陽(yáng)性克隆菌液用60%甘油保存.直接測(cè)序的產(chǎn)物，所用引物與PCR擴(kuò)增所用引物相同；克隆的樣品，采用相應(yīng)載體的通用引物作為測(cè)序引物.為確保所獲基因序列的準(zhǔn)確性，所有產(chǎn)物均送至上海生工生物工程有限公司進(jìn)行雙向測(cè)序.

1.2.6 序列分析 利用DNAMAN、BioEdit、Clustal W、MEGA 5.2、MrBayes 等生物信息學(xué)軟件對(duì)庫(kù)蠓基因序列進(jìn)行分析；同時(shí)，應(yīng)用MEGA 5.2 計(jì)算庫(kù)蠓種間與種內(nèi)各基因基于Kimura 雙參數(shù)模型的遺傳距離，并以與庫(kù)蠓屬親緣關(guān)系最近的細(xì)蠓屬中的北域細(xì)蠓(Leptoconops borealis)[9-10]作為外群，通過(guò)貝葉斯推理法(Bayesian inference， BI)構(gòu)建系統(tǒng)發(fā)育樹(shù).

2 結(jié)果與分析

2.1 各基因的遺傳距離

2.1.1 COⅠ 93.48%庫(kù)蠓的種間遺傳距離為0.15～0.30，平均距離為0.22(圖1A).其中，安酋庫(kù)蠓(C.andrewsi)與印度庫(kù)蠓(C.indianus)、美麗庫(kù)蠓(C.bellulus)與長(zhǎng)喙庫(kù)蠓(C.longirostris)、大和庫(kù)蠓(C.japonicus)與泊地庫(kù)蠓(C.toyamaruae)、馬來(lái)庫(kù)蠓(C.malayae)與類(lèi)牛庫(kù)蠓(C.parabubalus)之間的遺傳距離最小，僅為0.002，而琉球庫(kù)蠓(C.actoni)與新竹庫(kù)蠓(C.liui)之間的遺傳距離最大(0.37).65.85%庫(kù)蠓的種內(nèi)遺傳距離小于0.02，平均距離為0.03(圖1B).其中，美麗庫(kù)蠓、長(zhǎng)喙庫(kù)蠓、怒江庫(kù)蠓(C.nujiangensis)、異域庫(kù)蠓(C.peregrinus)、邊緣庫(kù)蠓(C.pictimargo)、色莖庫(kù)蠓(C.pictipennis)等6 個(gè)庫(kù)蠓的種內(nèi)遺傳距離為0，褐肩庫(kù)蠓(C.parahumeralis)的種內(nèi)遺傳距離最大(0.17).

圖1 庫(kù)蠓COⅠ基因序列K2P 遺傳距離頻度分布Fig.1 Frequency distribution of K2P distance of the COⅠgenes from Culicoides species

2.1.2 Cyt b 94.13%庫(kù)蠓的種間遺傳距離為0.15～0.30，平均距離為0.23(圖2A).其中，大和庫(kù)蠓與泊地庫(kù)蠓、馬來(lái)庫(kù)蠓與類(lèi)牛庫(kù)蠓之間的遺傳距離最小(0.002)，類(lèi)牛庫(kù)蠓與褐肩庫(kù)蠓之間的遺傳距離最大(0.35).69.44%庫(kù)蠓的種內(nèi)遺傳距離小于0.02，平均距離為0.04(圖2B).其中，安酋庫(kù)蠓、美麗庫(kù)蠓、沙庫(kù)蠓(C.desertorum)、赫氏庫(kù)蠓(C.hegneri)、屏東庫(kù)蠓(C.hui)、金娜庫(kù)蠓(C.kinari)、邊緣庫(kù)蠓、細(xì)須庫(kù)蠓(C.tenuipalpis)等8 種庫(kù)蠓的種內(nèi)遺傳距離為0，霍飛庫(kù)蠓(C.huffi)的種內(nèi)遺傳距離最大(0.19).

圖2 庫(kù)蠓Cyt b 基因序列K2P 遺傳距離頻度分布Fig.2 Frequency distribution of K2P distance of the Cyt b genes from Culicoides species

2.1.3 16S rDNA 96.15%庫(kù)蠓的種間遺傳距離為0.10～0.25，平均距離為0.16(圖3A).其中，大和庫(kù)蠓與泊地庫(kù)蠓之間的遺傳距離最小(0)，類(lèi)牛庫(kù)蠓與李拭庫(kù)蠓(C.riethi)之間的遺傳距離最大(0.27).72.09%庫(kù)蠓的種內(nèi)遺傳距離小于0.02，平均距離為0.02(圖3B).其中，美麗庫(kù)蠓、沙庫(kù)蠓、屏東庫(kù)蠓、赫氏庫(kù)蠓、大和庫(kù)蠓、金娜庫(kù)蠓、馬來(lái)庫(kù)蠓、迷你庫(kù)蠓(C.minimus)、異域庫(kù)蠓、泊地庫(kù)蠓等10 個(gè)蠓的種內(nèi)遺傳距離為0，細(xì)須庫(kù)蠓的種內(nèi)遺傳距離最大(0.09).

2.1.4 ITS 89.92%庫(kù)蠓的種間遺傳距離為0.15～0.30，平均距離為0.21(圖4A).其中，藍(lán)腹庫(kù)蠓(C.cyancus)與華鎣庫(kù)蠓(C.huayingensis)之間的遺傳距離最小(0.003)，克彭庫(kù)蠓(C.kepongensis)與異域庫(kù)蠓之間的遺傳距離最大(0.39).97.22%庫(kù)蠓的種內(nèi)遺傳距離在0.14 以內(nèi)，平均距離為0.04(圖4B).其中，盔狀庫(kù)蠓(C.cassideus)的種內(nèi)遺傳距離為0.003，短須庫(kù)蠓(C.brevipalpis)的種內(nèi)遺傳距離最大(0.17).

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

為確保系統(tǒng)發(fā)育分析結(jié)果的準(zhǔn)確性，對(duì)各個(gè)數(shù)據(jù)集的系統(tǒng)發(fā)育信號(hào)進(jìn)行評(píng)估.結(jié)果表明，各數(shù)據(jù)集的堿基替換均未飽和，可用于系統(tǒng)發(fā)育樹(shù)的分析.

本研究獲得線粒體COⅠ+16S rDNA+Cyt b 聯(lián)合基因158 條(片段長(zhǎng)度為1 530 bp)和COⅠ+16S rDNA+Cyt b+ITS 完全聯(lián)合基因118 條(片段長(zhǎng)度為2 787 bp)，隸屬于屋室亞屬(subgenusOecactaPoey)、帶紋亞屬(subgenusBeltranmyiaVargas)、庫(kù)蠓亞屬(subgenusCulicoidesLatreille)、二囊亞屬(subgenusAvaritiaFox)、血色亞屬(subgenusHaemophoructusMacfie)、三囊亞屬(subgenusTrithecoidesWirth et Hubert)和單囊亞屬(subgenusMonoculicoidesKhalaf).

圖3 庫(kù)蠓16S rDNA 基因序列K2P 遺傳距離頻度分布Fig.3 Frequency distribution of K2P distance of the 16S rDNA genes from Culicoides species

圖4 庫(kù)蠓ITS 基因序列K2P 遺傳距離頻度分布Fig.4 Frequency distribution of K2P distance of the ITS genes from Culicoides species

根據(jù)AIC(Akaike information criterion)原則選擇最佳核苷酸替換模型進(jìn)行系統(tǒng)發(fā)育分析，各單基因及聯(lián)合基因的最適模型都為GTR+I+G.從線粒體基因聯(lián)合數(shù)據(jù)集的系統(tǒng)發(fā)育樹(shù)(圖5)可知，不同亞屬之間構(gòu)成了13 個(gè)平行支系.其中以二囊亞屬為最大的支系，主要包含琉球庫(kù)蠓、迷你庫(kù)蠓組成的琉球庫(kù)蠓組(actoni group)以及東方庫(kù)蠓(C.orientalis)、牧場(chǎng)庫(kù)蠓(C.pastus)、暗色庫(kù)蠓(C.fulvus)、連斑庫(kù)蠓(C.jacobsoni)、美麗庫(kù)蠓、長(zhǎng)喙庫(kù)蠓、黑斑庫(kù)蠓(C.migromaculatus)、怒江庫(kù)蠓等組成的東方庫(kù)蠓組(orientalis group)；琉球庫(kù)蠓組與東方庫(kù)蠓組構(gòu)成姐妹群，并有較高的支持率(PP 值為0.73).此外，在另外一個(gè)較大的支系上，屋室亞屬與帶紋亞屬構(gòu)成姐妹群(PP 值為0.51).

基于完全聯(lián)合數(shù)據(jù)集的系統(tǒng)發(fā)育樹(shù)(圖6)分為2 個(gè)支系:其中一個(gè)小支系僅包含麥?zhǔn)蠋?kù)蠓(C.mcdonaldi)和刺螫庫(kù)蠓(C.punctatus)兩種；另一個(gè)支系拓?fù)浣Y(jié)構(gòu)基本同于線粒體基因聯(lián)合數(shù)據(jù)集，不同的是牧場(chǎng)庫(kù)蠓未與東方庫(kù)蠓組的其他蠓種聚在一起，而是與尖喙庫(kù)蠓(C.oxystoma)聚在一支.

在不同數(shù)據(jù)集的系統(tǒng)發(fā)育樹(shù)中，僅單囊亞屬是單系群，其余亞屬均不具有單系性.

3 討論與結(jié)論

本研究以近6 年來(lái)在我國(guó)14 個(gè)省(自治區(qū))采獲的庫(kù)蠓標(biāo)本為研究對(duì)象，初步建立了基于線粒體COⅠ、Cyt b、16S rDNA 和核糖體ITS 多基因的庫(kù)蠓分子分類(lèi)體系.該體系與形態(tài)學(xué)鑒定結(jié)果基本吻合，并從中發(fā)現(xiàn)了庫(kù)蠓屬7 個(gè)新種，9 個(gè)我國(guó)新記錄種(待發(fā)表).

種間與種內(nèi)的遺傳距離分析結(jié)果顯示，基于各基因序列K2P 校正的種間與種內(nèi)平均遺傳距離不同，以16S rDNA 的種間(0.16) 和種內(nèi)(0.02) 的平均遺傳距離最小，COⅠ(種間0.22、種內(nèi)0.03)和ITS(種間0.21、種內(nèi)0.04)的平均遺傳距離次之，而Cyt b 的種間(0.23)和種內(nèi)(0.04)的平均遺傳距離最大.在不同的基因標(biāo)記中，分布范圍較廣的蠓種的種內(nèi)遺傳距離常較大，如褐肩庫(kù)蠓基于COⅠ的種內(nèi)遺傳距離達(dá)0.17；而分布范圍較窄的蠓種的種內(nèi)遺傳距離較小，如美麗庫(kù)蠓基于COⅠ、Cyt b 以及16S rDNA 的種內(nèi)遺傳距離都為0.這提示大的種內(nèi)變異可能產(chǎn)生于地理隔離的不同族群[11]或存在隱存種[12].一些形態(tài)近似的庫(kù)蠓的種間遺傳距離較小，如安酋庫(kù)蠓與印度庫(kù)蠓、美麗庫(kù)蠓與長(zhǎng)喙庫(kù)蠓、馬來(lái)庫(kù)蠓與類(lèi)牛庫(kù)蠓、大和庫(kù)蠓與泊地庫(kù)蠓、藍(lán)腹庫(kù)蠓與華鎣庫(kù)蠓，基于不同基因的遺傳距離僅為0～0.003.

圖5 基于COⅠ+Cyt b+16S rDNA 聯(lián)合基因的BI 系統(tǒng)發(fā)育樹(shù)Fig.5 Phylogenetic tree based on COⅠ+Cyt b+16S rDNA by Bayesian inference

系統(tǒng)發(fā)育樹(shù)表明，僅單囊亞屬是單系群，屋室亞屬、帶紋亞屬、庫(kù)蠓亞屬、二囊亞屬、血色亞屬、三囊亞屬都形成并系群，這與早先的研究結(jié)果[12]一致.

圖6 基于COⅠ+Cyt b+16S rDNA+ITS 聯(lián)合基因的BI 系統(tǒng)發(fā)育樹(shù)Fig.6 Phylogenetic tree based on COⅠ+Cyt b+16S rDNA+ITS by Bayesian inference

不同基因數(shù)據(jù)集的聚類(lèi)結(jié)果提示:(1)安酋庫(kù)蠓(同龍溪庫(kù)蠓C.lungchiensis)、印度庫(kù)蠓、異域庫(kù)蠓可能為同種異名.龍溪庫(kù)蠓與安酋庫(kù)蠓沒(méi)有明顯區(qū)別，因后者發(fā)表在先，故統(tǒng)一更名為安酋庫(kù)蠓；安酋庫(kù)蠓與印度庫(kù)蠓也非常近似[13-14]，且兩者間基于COⅠ的遺傳距離僅為0.002，明顯低于學(xué)者基于DNA 條形碼對(duì)種間遺傳距離大于2%的界定[15]；Borkent et al[16]將龍溪庫(kù)蠓(安酋庫(kù)蠓)作為異域庫(kù)蠓的同種異名，這與本研究中安酋庫(kù)蠓與異域庫(kù)蠓形成姐妹群的結(jié)果相一致.(2)馬來(lái)庫(kù)蠓、類(lèi)牛庫(kù)蠓、南山庫(kù)蠓(C.lansangensis)、蘇島庫(kù)蠓(C.sumatrae)可能為一個(gè)種團(tuán).在不同基因的系統(tǒng)發(fā)育樹(shù)上，馬來(lái)庫(kù)蠓與類(lèi)牛庫(kù)蠓都為姐妹群，且兩種間遺傳距離小(COⅠ和16S rDNA 均為0.002)；從形態(tài)上看，這兩種蠓的鑒定特征基本相似，用于辨別的主要特征僅為前者的受精囊一大一小[17]，后者的兩個(gè)受精囊基本等大[18].系統(tǒng)發(fā)育樹(shù)上，南山庫(kù)蠓、蘇島庫(kù)蠓與馬來(lái)庫(kù)蠓、類(lèi)牛庫(kù)蠓聚于同一支且有高的支持率，PP 值均為1.00；形態(tài)上，這4 種庫(kù)蠓的翅斑也較為近似.(3)美麗庫(kù)蠓與長(zhǎng)喙庫(kù)蠓可能是種內(nèi)的形態(tài)變異體或同種多態(tài)型.美麗庫(kù)蠓翅面徑脈端部2個(gè)明斑明顯呈不同程度的傾斜，而長(zhǎng)喙庫(kù)蠓的明斑偏直，且前者的翅面明斑明顯小于后者；美麗庫(kù)蠓雄蟲(chóng)的尾器第9 背板不呈弧形凹陷，而長(zhǎng)喙庫(kù)蠓有寬而深的凹陷；前者雌蟲(chóng)的2 個(gè)受精囊不等大，而后者雌蟲(chóng)的2 個(gè)受精囊基本等大.但兩者在系統(tǒng)發(fā)育樹(shù)上為姐妹群，且基于不同基因數(shù)據(jù)集的遺傳距離均很小(COⅠ為0.002，16S rDNA 為0.001，Cyt b 和ITS 均為0.008).(4)泊地庫(kù)蠓可能是大和庫(kù)蠓的近緣種.系統(tǒng)發(fā)育樹(shù)顯示，泊地庫(kù)蠓與大和庫(kù)蠓是單系群，且PP 值為1.00，這與兩者間較小的遺傳距離(COⅠ和Cyt b均為0.002，16S rDNA 和ITS 均為0)相符.但兩者在形態(tài)上有明顯區(qū)別:大和庫(kù)蠓翅面中肘叉和縱二室有明斑，而泊地庫(kù)蠓通常沒(méi)有，且前者的翅面明斑大于后者；前者雌蟲(chóng)的受精囊呈卵形，而后者呈圓形.(5)注意迷你庫(kù)蠓與琉球庫(kù)蠓的準(zhǔn)確鑒別.兩者都屬于小型蠓種(體長(zhǎng)約為1 mm)，是琉球庫(kù)蠓組里僅有的2 個(gè)種.形態(tài)上除受精囊大小外并無(wú)明顯區(qū)別，琉球庫(kù)蠓的兩個(gè)受精囊近乎等大，而迷你庫(kù)蠓的兩個(gè)受精囊大小不一；基于COⅠ基因的系統(tǒng)發(fā)育樹(shù)顯示，兩者形成姐妹群.因此，單靠形態(tài)學(xué)可能會(huì)造成錯(cuò)誤的鑒定.

利用線粒體COⅠ、16S rDNA、Cyt b 和核糖體ITS 基因不僅能準(zhǔn)確鑒定蠓種，還能澄清庫(kù)蠓屬里的同種異名問(wèn)題.但本研究仍有一些方面待改進(jìn):部分標(biāo)本因形態(tài)不完整而無(wú)法鑒定到具體種，或因保存方式不當(dāng)而未能提取出核酸；同時(shí)，所有試驗(yàn)標(biāo)本均為單只處理，且僅有胸段翅基部組織用于核酸提取，提取的總DNA 濃度低、體積小，導(dǎo)致部分標(biāo)本的某些基因片段擴(kuò)增失敗，造成數(shù)據(jù)不全；再者，吉林亞屬(subgenusJilinocoides)、橋莖亞屬(subgenusPontoculicoides)等個(gè)別亞屬種類(lèi)數(shù)太少，分布地點(diǎn)較為局限，本研究未能采獲這些亞屬的標(biāo)本.因此，今后有待補(bǔ)充標(biāo)本并做進(jìn)一步研究.