中國縞蠅科(雙翅目：縞蠅總科)高級階元的幾何形態學分析

戚凌， 何向怡， 李文亮

(河南科技大學林學院，河南洛陽471023)

雙翅目Diptera的分類較為復雜，迄今未形成統一的分類體系。目前一般按照Delfinado和Hardy(1977)、Richards和Daries(1977)的分類標準，將縞蠅科Lauxaniidae劃屬于雙翅目環裂亞目Cyclorrhapha無瓣類Acalyptrata有縫組Schizophora縞蠅總科Lauxanioidea(高彩霞，2003)；目前世界已知 3亞科172屬約2 150種，我國已知32屬343種(李文亮，2014)。縞蠅總科是無瓣蠅類各總科系統進化過程中較為原始的類群之一，對縞蠅科昆蟲的研究是揭示無瓣蠅類各總科間系統發育關系的重要環節(李文亮，2014)。在縞蠅科分類中，往往使用傳統的形態學方法鑒定，即觀察記述標本的外部分類特征，然后進一步解剖、鑒定、繪圖等，需要花費大量的時間和精力。

近年來幾何形態分析方法越來越成熟，并在分類學、系統發育學、生物進化學等方向有所應用。20世紀早期，生物學從定性描述向定量描述過渡；20世紀中期，傳統的形態測量學產生；20世紀80年代，形態結構數值化以及數據分析技術得到重要突破，生物體形態的幾何信息能夠被獲得，并用于綜合分析(Rohlfetal.，1998；Perrardaetal.，2016)。這種方法利用了變量之間的角度轉換、物體表面的紋理甚至實物標本的物理性能等數據信息，能夠將生物體形態信息中的大小和形狀分開(白明等，2014)。目前在幾何形態學分析中，最常見的方法是標點法，該方法能夠找到研究樣本之間的形態“同源信號”和“簡約信息位點”，排除樣本之間由于位置排列產生的差異(Yanetal.，2011)。

翅是昆蟲獲得飛翔能力的重要形態結構，翅的形狀與昆蟲的生存方式、擴散能力有著密切的關系，翅形的變化有可能對物種的分化形成起到關鍵的作用(白明，楊星科，2007)。李玲等(2017)曾基于幾何形態學對3種實蠅翅脈形態進行了分析，并證明幾何形態學方法可用于實蠅種間的分類。段琛等(2018)也通過幾何形態學方法對3種庫蠓的翅形變化進行了分析，并成功找出了翅的主要變異區域。

目前尚無用幾何形態學方法對縞蠅科翅脈進行研究的報道，而本試驗運用幾何形態學分析中的標點法分析了翅脈在縞蠅科分類標準上的科學依據，并利用幾何形態學的原理對縞蠅科高級階元進行初步分析，探討翅位點的變異度，為縞蠅科系統分類開拓新的途徑，并為縞蠅科幾何形態學系統發育研究打下初步基礎。

1 材料與方法

1.1 標本材料

縞蠅科昆蟲標本來源于河南科技大學昆蟲標本館和中國農業大學昆蟲博物館，共2亞科20屬252種(表1)。

1.2 試驗方法

1.2.1 獲取圖像用EOS450D/5D相機拍攝照片，再用Adobe Photoshop CS將圖片設置為灰度模式，并調整翅脈的角度。

1.2.2 數據提取及標準化處理用tps Util將翅脈圖片格式文件轉化為tps格式文件，然后分別對樣本進行同源標點，每個樣本均按照同樣的位置、順序標注14個標點(圖1)：(1)前緣脈(C)與肩脈(h)交點，(2)C與亞前緣脈(Sc)交點，(3)Sc與第一徑脈(R1)交點，(4)第二、三合徑脈(R2+3)與翅緣交點，(5)第四、五合徑脈(R4+5)與翅緣交點，(6)第一中脈(M1)與翅緣交點，(7)第一肘臀合脈(CuA1)與翅緣交點，(8)第二臀脈(A2)末端，(9)第二肘臀合脈(A1+CuA2)末端，(10)徑中橫脈(r-m)與R4+5交點，(11)r-m與M1交點，(12)中肘橫脈(dm-cu)與M1交點，(13)dm-cu與CuA1交點，(14)翅外緣黑色短鬃的末端。

表1 縞蠅科標本種類統計Table 1 Statistics of specimens in Lauxaniidae family

1.2.3 數據分析利用MorphoJ基于主成分分析(principal component analysis，PCA)和典型變量分析(canonical variate analysis，CVA)對翅脈的幾何形態學數據進行分析。先進行原始數據的主成分分析，再從不同方向分組進行不同階級的主成分分析和典型變量分析，并通過PCA和CVA的薄板樣條圖，詳細解析每個標點所代表的翅脈的差異性。

2 結果

2.1 縞蠅科總體分析

全部樣本的PCA結果(圖2)顯示：在提取的24個主成分中，第一主成分(PC1)占總變異量的30.010%，第二主成分(PC2)占總變異量的18.698%，第三主成分(PC3)占總變異量的15.377%，前三個主成分的特征值占總變異量的64.085%，能夠說明各樣本之間的主要差異。

薄板樣條分析結果(圖3)顯示：在PC1中，扭曲度最大的是標點4和14；在PC2中，標點10和11的扭曲度相對較大。說明在縞蠅科的進化與變異中，R2+3、黑色短鬃和r-m是變異度較高的翅脈，在進化中相對較不穩定。

圖1 縞蠅科昆蟲翅的標點Fig. 1 14 landmarks on the wing of Lauxaniidae

C. 前緣脈， h. 肩脈， Sc. 亞前緣脈， R1. 第一徑脈， R2+3. 第二、三合徑脈， R4+5. 第四、五合徑脈， M1. 第一中脈， CuA1. 肘脈(Cu)與第一臀脈(A1)的部分段合并而成的第一肘臀合脈， A2. 第二臀脈， r-m. 徑中橫脈， dm-cu. 中肘橫脈， A1+CuA2. Cu與第一、二臀脈(A1、A2)合并而成的第二肘臀合脈

C. costa， h. humeral vein， Sc. subcosta， R1. first radial vein， R2+3. combined vein of second and third radial vein， R4+5. combined vein of fourth and fifth radial vein， M1. first medial vein， CuA1. combined vein of cubital and first anal vein， A2. second anal， r-m. radial-median crossvein， dm-cu. median-cubital crossvein， A1+CuA2. combined vein of cubital， first and second anal vein

圖2 縞蠅科昆蟲翅脈主成分分析Fig. 2 Principal component analysis for the wing veins of Lauxaniidae

圖3 縞蠅科昆蟲翅脈PC1與PC2的薄板樣條Fig. 3 Thin plate spline of PC1 and PC2 for the wing veins of Lauxaniidae

1～14見圖1

1-14 see Fig.1

2.2 同脈縞蠅亞科與縞蠅亞科的分析

在總體分析的基礎上，已知前3個主成分的特征值占總變異量的64.085%。縞蠅科2亞科的PCA結果(圖4)顯示：同脈縞蠅亞科主要聚集在PC1軸的0點以左(-0.10～0)，而縞蠅亞科主要聚集在PC1軸的0點以右(0～0.01)；兩者在PC2方向上并沒有明顯差異。縞蠅科2亞科的CVA結果(圖5)顯示：同脈縞蠅亞科與縞蠅亞科僅在橫軸的0點附近有小部分重合，而絕大部分區間是非重合的。同脈縞蠅亞科與縞蠅亞科在PC1和CV1方向上具有明顯的差異，說明這2個亞科能夠通過幾何形態學方法分析翅脈特征來進行分類鑒定。

圖4 同脈縞蠅亞科(紅色)與縞蠅亞科(藍色)翅脈的主成分分析散點
Fig. 4 Principal component analysis for the wing veins of Homeoneuriae (red) and Lauxaniinae (blue)

圖5 同脈縞蠅亞科(紅色)與縞蠅亞科(藍色)翅脈的典型變量分析柱形Fig. 5 Canonical variate analysis for the wing veins of Homeoneuriae (red) and Lauxaniinae (blue)

CV1與PC1的薄板樣條圖(圖6)顯示：扭曲度最大的點均為標點4、14，扭曲度較小的均為標點1、5和6。說明R2+3和黑色短鬃的變異度最大，而C、R4+5和M1的變異度較大。

2.3 縞蠅科4屬間的分析

為了找出縞蠅科在屬級特征上的翅脈變化，本試驗選擇了4個樣本數量相對較多的屬(表2)，即同脈縞蠅屬Homoneura、黑縞蠅屬Minettia、四帶縞蠅屬Tetroxyrhina和斑縞蠅屬Trypetisoma進行CVA。在CVA中，前3個典型變量分別占總變異量的74.767%、18.956%和6.277%，共計100%。

CVA結果(圖7)顯示：不同屬各自聚集在不同區域，說明這4屬在翅脈上能作出明顯區分。

薄板樣條圖(圖8)顯示：在CV1薄板樣條中，標點4和14相對扭曲度較大，且它們的變異趨勢相反，說明R2+3與翅緣的交點向上變異，黑色短鬃則越來越靠近標點5，R2+3與短鬃間的距離是逐漸增大的。在CV2薄板樣條中，扭曲度較大的是標點10和11，說明r-m有相對右移且增長的變異趨勢；標點1、2、6和7的變異度極小，說明這些點所對應的C、Sc、M1和CuA1是相對穩定、不易變異的。

圖6 同脈縞蠅亞科與縞蠅亞科翅脈CV1和PC1的薄板樣條
Fig. 6 Thin plate spline of CV1 and PC1 for the wing veins of Homeoneuriae and Lauxaniinae

表2 縞蠅科4屬間分析的樣本數量Table 2 Number of samples for the 4 genera of Lauxaniidae

3 結論與討論

3.1 結論

總體分析、兩大亞科分析和4個常見屬分析結果顯示，R2+3與翅緣交點、黑色短鬃端點和r-m端點是變異度較大的位點，R4+5與翅緣交點、M1與翅緣交點、CuA1與翅緣交點為穩定的位點；翅緣黑色短鬃端點可作為亞科的重要鑒別特征，這與傳統的縞蠅科分類標準相吻合，證明了本次試驗的可行性與科學性；4屬分析中屬間離散度較高，說明標點法能較好地進行屬級階元的鑒別。

3.2 討論

本試驗得到的幾何形態學分析結果與傳統形態學鑒定結果相一致，即二者的相互印證度較好，證明了幾何形態學方法在縞蠅分類中應用的可行性。但本試驗采用的部分縞蠅屬的樣本數量較少，且每種僅選取了1個樣本，可能會對試驗結果有一定影響。在屬間分析中，4個常見屬不足以代表整個縞蠅科，下一步可以增加各個種的樣本數量，擴充樣本后進一步在低級階元如屬間或種間進行分析研究。此外，本試驗選取的標點數量較少，且未結合翅輪廓的變化進行比較，研究結果只能確定一些差異性較大的翅脈。而翅的輪廓特征也是相對較重要的翅特征之一，曾有用半標點法結合翅輪廓對蟻蛉科Myrmeleontidae昆蟲進行幾何形態學分析的研究(詹慶斌，2014)，所以可以更進一步結合縞蠅翅輪廓的研究來補充證明翅在分類依據上的重要性。

圖7 縞蠅科4屬翅脈的典型變量分析散點
Fig. 7 Canonical variate analysis for the wing veins of the 4 genera of Lauxaniidae

H. 同脈縞蠅屬Homoneura， M. 黑縞蠅屬Minettia， TET. 四帶縞蠅屬Tetroxyrhina， TRY. 斑縞蠅屬Trypetisoma

圖8 縞蠅科4屬翅脈的CV1與CV2薄板樣條Fig. 8 Thin plate spline of CV1 and CV2 for the wing veins of 4 genera of Lauxaniidae

1～14見圖1

1-14 see Fig.1

除了將翅作為分析對象，幾何形態學方法還能夠通過昆蟲其他部位的特征來進行分類鑒定等研究。在關于昆蟲分類地位的研究中，李榮榮等(2019)通過幾何形態學方法對二星蝽屬Eysarcoris及其近緣屬昆蟲的小盾片形狀進行了分析，并根據試驗結果將二星蝽屬劃歸于二星蝽族。在有關昆蟲進化發育的研究中，曾有利用幾何形態學方法對鍬甲的前胸背板與鞘翅進行分析從而探究其演變進化的試驗，該研究使用半標點法標注了鍬甲前胸背板與鞘翅的外輪廓，并結合分子數據建立了系統發育樹，對鍬甲的外部形態演變作出了推斷(Zhangetal.，2019)。在昆蟲形態與功能關系的研究中，佟一杰等(2016)對蜣螂后胸叉骨的形態多樣性進行了幾何形態學定量評估，證明滾糞球或直接掘洞的行為會對蜣螂后胸叉骨的形態有顯著性影響。基于上述成功的研究案例，本試驗下一步考慮將除了翅以外的其他部位特征作為幾何形態學分析對象。在縞蠅科昆蟲中，背中鬃著生位點是一個顯著的形態分類特征，所以本試驗將嘗試用幾何形態學方法分析中鬃著生位點，以期能夠與翅的分析結果相結合，對縞蠅科昆蟲進行更準確的分類。

近年來，昆蟲學依據其獨有的特性及后發優勢，將幾何形態學的應用推向一個新的高度，未來將不斷開拓新的研究模式與方向。昆蟲占領了眾多生態位，并產生了極為多樣的形態適應，通過幾何學與發育生物學、支序系統學、古生物學、生態學等方面證據的深度結合與相互印證，可從獨特的角度對相關科學問題進行解讀(Baylacetal.，2003)。