中國青鳳蝶3個亞種的遺傳變異分析

楊斯琦　王亞楠　王方方　何婭　靳亞麗

摘要[目的]探討中國青鳳蝶（Graphium sarpedon Linnaeus）3個亞種的遺傳變異。[方法]對青鳳蝶3個亞種的COⅠ基因和EF-1α基因部分序列進行測定，并對序列的堿基組成、轉換顛換數和遺傳距離等進行分析。[結果]在測得的COⅠ基因（709 bp）和EF-1α（1 064 bp）基因中，有44個變異位點，6個簡約信息位點，COⅠ基因A+T平均含量為69.3%，存在較強的含量偏向性，亞種間的遺傳距離相差甚小。[結論]青鳳蝶3個亞種之間沒有明顯序列差異，仍然適宜采用傳統的基于形態學的青鳳蝶亞種分類體系。

關鍵詞青鳳蝶；COⅠ基因；EF-1α基因；遺傳變異

中圖分類號Q969.42文獻標識碼A文章編號0517-6611（2017）31-0156-04

Abstract [Objective]To investigate the genetic variation in three subspecies of Graphium sarpedon Linnaeus in China.[Method] The COⅠ gene and the EF1α gene were partially sequenced.The nucleotide composition， the number of transition and transversion， genetic distance of the segment had been analyzed.[Result]There were 44 variation sites and 6 parsimonyinformative sites in combined sequence of the COⅠ gene（709 bp）and EF1α（1 064 bp） gene in the examined species. The average A+T content of COⅠ gene was 69.3%，which showed strong A+T bias. But genetic distance among subspecies were very small.[Conclusion]It has no significant sequence differences between the three subspecies， which is still suitable to use the traditional morphological classification system for Graphium sarpedon subspecies.

Key wordsGraphium sarpedon；COⅠ gene； EF1α gene；Genetic variation

青鳳蝶（Graphium sarpedon Linnaeus），又名樟青鳳蝶或青帶鳳蝶，屬于鱗翅目（Lepidoptera），鳳蝶科（Papilionidae），青鳳蝶屬（Graphium），分布于我國陜西、湖北、湖南、四川、云南、貴州、西藏、江西、浙江、福建、廣西、廣東、海南、臺灣、香港等地區[1]。青鳳蝶在我國有3個亞種，即指名亞種[G.sarpedon sarpedon（Linnaeus）]、藍斑亞種[G. sarpedon connectens（Fruhstorfer）]、斑帶亞種[Graphium sarpedon semifasciatum（Honrath）]。目前對青鳳蝶斑帶亞種的分類尚存有爭議，魏忠民等[2]通過對中國青鳳蝶標本的觀察，歸納出6種中國青鳳蝶的變異型，描述其主要形態特征，認為青鳳蝶的斑帶亞種是一種變型，而不是亞種。對于青鳳蝶亞種的爭議問題，以DNA序列為基礎的分辨依據研究較少。

目前國內外學者結合線粒體基因和核基因等分子標記

對昆蟲類群的系統學研究已有大量的文獻報道。線粒體DNA具有結構簡單、進化速度快、高拷貝數、易于分離純化等特點，是研究系統發育、種群遺傳變異和分化、區分近緣種的重要工具[3]。由于線粒體基因是母系遺傳，其中所含的進化信息并不能完全反映雙親進化的歷史，核基因也含有豐富的生物學信息，結合適當的核基因將更好地分析蝶類的系統發育[4]。

筆者對我國分布的青鳳蝶3個亞種的COⅠ基因和EF-1α基因部分序列進行測定，研究它們在DNA序列上的差異，以期為我國青鳳蝶3個亞種分類地位提供分子證據。

1材料與方法

1.1材料

青鳳蝶3個亞種的標本來自我國廣西、上海，包括作為外群的木蘭青鳳蝶，共 15只成蟲標本（表1）。

1.2方法

1.2.1提取基因組DNA及PCR擴增。

選取部分腹部或足部組織作為試驗材料，采用試劑盒insect kit（Omega）提取基因組DNA，樣品保存在-20 ℃備用。試驗采用的引物序列如下。

COⅠ上游引物：5′-GGTCAACAAATCATAAAGATATTG-3′；

下游引物：5′-TAAACTTCAGGGTGACCAAAAAAT-3′。

EF-1α上游引物：5′-TCGATATCGCTTTGTGGAAGTT-3′；

下游引物：5′-ACGCACGGCAAAACGACCGAGAG-3′。

PCR 反應體積為50 μL，反應體系包括1.25 U的 Taq DNA聚合酶，0.2 mmol/L dNTP，上下游引物分別為0.4 μmol/L，10×PCR Buffer 為5 μL，模板 DNA 溶液 2 μL（50～100 ng DNA）。擴增條件為95 ℃預變性3 min；95 ℃變性30 s，55 ℃復性30 s，72 ℃延伸1 min，30個循環；最后在72 ℃延伸7 min，4 ℃保溫。用 Eppendorf 梯度PCR儀進行擴增。

1.2.2PCR產物純化、測序。

用1%的瓊脂糖凝膠電泳檢測 PCR 產物，用DNA 凝膠回收試劑盒，對目的條帶進行回收和純化后，連接到T載體，經過轉化挑克隆，提取質粒再委托上海賽默飛（Thermo Fisher）公司進行雙向測序。

1.2.3序列分析及差異比較。

測定15只成蟲標本的COⅠ和EF-1α基因部分序列，正反鏈序列拼接后，再Blast進行序列同源性比較，基于Kimura-2參數，采用MEGA6.0軟件進行序列組成分析，如序列堿基組成（nucleotide composition）、保守位點（conservedsjtes）、變異位點（variable sjtes）、簡約信息位點（parsimony information sites）、自裔位點（singletonsites）等。選擇Tamurar Nei模型，計算各分類單元之間的遺傳距離、轉換數與顛換數比值（Ts/Tv）。

2結果與分析

2.1基因序列的確定

根據所選引物的位置，COⅠ基因和EF-1α基因預測擴增的片段分別是709 bp和1 064 bp，電泳檢測的結果與預期相符（圖1）。運用Blast對序列進行分

析，顯示其與青鳳蝶的COⅠ基因和EF-1α基因具有很高

的同源性。

2.2序列組成分析

除外群外，COⅠ基因部分序列共709 bp，變異位點6 bp，無簡約信息位點，自裔位點6 bp，堿基組成上 T、C、A、G的含量分別是39.1%、16.5%、30.2%、144%，A+T平均含量為69.3%，其中密碼子第2位點A+T含量最高，達90.5%，第2位點G含量最低，平均為0.4%，T含量最高，達49.0%，這表明密碼子的堿基使用頻率存在明顯的偏向性（表2）。

核EF-1α基因部分序列共1 064 bp，變異位點38 bp，簡約信息位點6 bp，自裔位點32 bp，T、C、A、G的含量分別是22.1%、28.6%、24.5%、24.8%，A+T平均含量為46.6%，其中第1位點含量最高，達到59.2%，第1位點G的含量最低，平均為15.2%，T的含量最高，達到28%，表明密碼子的堿基使用頻率也存在一定的偏向性（表3）。

通過Paup4.0b10軟件對COⅠ基因和EF-1α基因進行PHT檢驗（同質性檢驗partition homogeneity test），結果顯示基因進化水平不具有顯著差異（P=1.0>0.01），因而可以將2組數據整合在一起進行分析。利用MEGA6.0軟件，基于Kimura-2參數對COⅠ和 EF-1α基因部分序列的整合數據進行分析，統計堿基轉換數與顛換數比值（R）。

由表4可知，核苷酸替換發生率比較低，第3位點略高，轉換多于顛換，第1位點最保守。轉換數與顛換數比值（R）平均為1.5<2，表現出較為明顯的替換飽和。

2.3遺傳距離分析

基于Kimura-2參數，采用MEGA 6.0計算個體間及亞種間的遺傳距離（表5）。不含外群的個體兩兩之間的遺傳距離在0.1%～0.5%，平均遺傳距離是03%，將3個亞種進行分組后計算遺傳距離，并與外群對照（表6），我國青鳳蝶3個亞種間遺傳距離是0.30%～034%，兩兩之間的遺傳距離相差很小，其中指名亞種[G.sarpedon sarpedon（Linnaeus）]與斑帶亞種[G.sarpedon semifasciatum（Honrath）]差異相對更小，與藍斑亞種[G.sarpedon connectens（Fruhstorfer）]的差異相對更大。

2.4系統發育信號檢測

采用堿基替換飽和性分析，對數據進行系統發育信號強弱的檢測，以遺傳距離為橫坐標，以堿基轉換數和顛換數為縱坐標作散點圖（圖2）。

從圖2可以看出，隨著遺傳距離的增大，Ts增加的速度大于Tv增加的速度，且與遺傳距離之間有良好的線性關系，Tv則逐漸趨于穩定，達到飽和。

通過Paup4.0b10軟件對數據做PTPtest序列信息檢驗，結果顯示P為0.002，表明該數據具有較顯著的系統發育信息，而非隨機數據，可用來進行系統發育的推斷[5-6]。

3討論

我國青鳳蝶分布有3個亞種，據《中國蝶類志》中描述，青鳳蝶的形態特征為翅窄長，底色黑，無尾狀突起，前后翅中央貫穿1列略呈方形的藍綠色斑（后翅前面的1個為白色），后翅外緣有1列綠藍色的新月斑。后翅反面近翅基有1條紅色短線，翅中部至后緣處有數條紅色斑紋。

3個亞種之間形態特征存在一些差異，青鳳蝶藍斑亞種的綠色斑偏藍，前翅頂角有1個綠斑特別小；斑帶亞種的后翅色帶不全。有學者提出青鳳蝶的斑帶亞種是一種變型，而不是亞種。為此，該研究聯合2個基因序列的整合數據對青鳳蝶3亞種間遺傳變異進行了分析。之所以選擇COⅠ和EF-1α基因聯合分析，一方面是增加系統發生分析的置信度[7-8]，另一方面也與2個基因各自的特點有關。線粒體COⅠ是編碼線粒體細胞色素氧化酶亞基Ⅰ的基因，是常用的分子標記之一，它既相對保守又存在高變區，高變區內遺傳進化速率更快，種間的遺傳差異也更明顯，適合于種間和種內分類鑒定及系統發生研究[9]。EF-1α基因進化速度比較快，它也用于分析低階元的系統發育，已有學者做了相關的研究，如Reed等[10]用線粒體COⅠ、COⅡ基因的全長和核基因的EF-1α基因聯合分析了鳳蝶科鳳蝶屬23個種和亞種的個體。

在鱗翅目同種個體之間，COⅠ基因的差異為025%[11]，從試驗結果來看，我國青鳳蝶3個亞種間COⅠ基因相似

度為 99.15%，核苷酸差異很小，表明COⅠ基因在個體間變異率較低。3個亞種間EF-1α基因相似度96.43%，表明存在一定的變異率。3亞種青鳳蝶在形態上的差異可能

是由于地理隔離形成，其自身的遺傳變異性非常低[12]。

從系統發育樹的拓撲結構看，其未顯示3個亞種間有很好的單系性，而出現了并系性，且各支自舉檢驗值較低（<50%），因而沒有列出發育樹，置信度低下可能是因為分子數據所提供的信息位點較少[13]，或者是COⅠ基因和EF-1α基因序列突變達到飽和。結合對3個亞種基因序列兩兩距離進行P-distance計算，并按照亞種分組后計算遺傳距離，可以看出指名亞種和斑帶亞種的親緣關系較近，但不足以認定斑帶亞種是屬于指名亞種的一種變型。

綜上所述，結合目前已有的2個基因部分序列和各亞種之間形態特征的差異，仍然適宜采用傳統的基于形態學的青鳳蝶亞種分類體系。

參考文獻

[1] 武春生.中國動物志：昆蟲綱 第25卷 鱗翅目 鳳蝶科[M].北京：科學出版社，2001.

[2] 魏忠民，武春生.中國青鳳蝶變異型初步觀察[J].昆蟲知識，2006，43（3）：431-432.

[3] 壽建新.國內外蝴蝶分類認識總結[J].西安文理學院學報（自然科學版），2014，17（4）：21-27.

[4] 袁鋒，袁向群.蝶類分子系統學研究進展[J].西北農業學報，2013，22（12）：1-14.

[5] KIM K J，JANSEN R K.A chloroplast DNA phylogeny of lilacs（Syringa，Oleaceae）：Plastome groups show a strong correlation with crossing groups[J].American J Bot，1998，85（9）：1338-1351.

[6] AUSTIN C C，ZUG G R.Molecular and morphological evolution in the southcentral Pacific skink Emoia tongana（Reptilia：Squamata）：Uniformity and humanmediated dispersal[J].Australian J Zool，1999，47（5）：425-437.

[7] MONTEIRO A，PIERCE N E.Phylogeny of Bicyclus（Lepidoptera：Nymphalidae）inferred from COI，COII and EF1α gene sequences[J].Mol Phylogenet Evol，2001，18（2）：264-281.

[8] NORMARK B B.Molecular systematics and evolution of the aphid family Lachnidae[J].Mol Phylogenet Evol，2000，14（1）：131-140.

[9] 吳冬霞，朱國萍，陳娜，等.線蛺蝶亞科蝶類部分類群線粒體COI 基因的分子系統發生分析[J].生命科學研究，2007，11（1）：64-71.

[10] REED R D，SPERLING F A.Interaction of process partitions in phylogenetic analysis：An example from the swallowtail butterfly genus Papilio[J].Mol Biol Evol，1999，16（2）：286-297.

[11] HEBERT P D N，STOECKLE M Y，ZEMLAK T S，et al.Identification of birds through DNA barcodes[J].PLoS Biol，2004，2（10）：1657-1663.

[12] 葉維萍，葉海燕，盧慧甍，等.基于線粒體12S rRNA和ND5基因序列的中國飛蝗屬3亞種系統發育關系研究[J].昆蟲分類學報，2005，27（1）：5-12.

[13] 劉曉燕，吳孝兵，諸立新.中國黃粉蝶亞科六屬間基于COⅡ和EF-1α基因部分序列的系統發育關系（鱗翅目：粉蝶科）[J].昆蟲學報，2007，50（6）：604-609.