2urn:lsid:arphahub.com:pub:45048D35-BB1D-5CE8-9668-537E44BD4C7Eurn:lsid:zoobank.org:pub:91BD42D4-90F1-4B45-9350-EEF175B1727AZooKeysZK1313-29891313-2970Pensoft Publishers10.3897/zookeys.502.91915222Research ArticleCantharidaeTheory & MethodologyNeogeneChinaSignificance of hind wing morphology in distinguishing genera and species of cantharid beetles with a geometric morphometric analysisSuJunyan1GuanKaile1WangJiaxu1YangYuxiayuxia0305@126.comhttps://orcid.org/0000-0002-3118-66591Hebei UniversityBaodingChina
Corresponding author: Yuxia Yang (yxyang@hbu.edu.cn)
Academic editor: L. Penev
2015040520155021125FFC5FFE3-FF8A-FF97-2B73-6F0C7B099778DE5DFB51-580E-4CA6-A3A5-DAB51EA7ECD25788373012201423042015Junyan Su, Kaile Guan, Jiaxu Wang, Yuxia YangThis is an open access article distributed under the terms of the Creative Commons Attribution License (CC BY 4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.http://zoobank.org/DE5DFB51-580E-4CA6-A3A5-DAB51EA7ECD2
There remain some difficulties in delimitation of related genera or sibling species for cantharid beetles, because the traditionally taxonomic method and morphological characters have not been updated or introduced. In the present study, we firstly use the landmark-based geometric morphometrics to analyze and compare the hind wings of nine species belonging to three genera of Cantharinae to ascertain whether this approach may be used as a reliable method in the study of the taxonomy of this group. The results show that the shape differences of the hind wings among genera seem more variable than that within each genus, and the variations for each species are different from one another, as shown in the principal component analyses. And the canonical variates analyses show that there are significant differences among the genera and the species of each genus, which demonstrates that the hind wing shape can be diagnostic for both generic and specific identification of the cantharid beetles. This study sheds new light into clarifying the taxonomic uncertainties of Cantharidae, and lays a foundation for further studies on the evolution of the cantharid hind wing shape.
Guan K, Su J, Wang J, Yang Y (2015) Significance of hind wing morphology in distinguishing genera and species of cantharid beetles with a geometric morphometric analysis. ZooKeys 502: 11–25. doi: 10.3897/zookeys.502.9191
Introduction
The Cantharinae represents a subfamily of beetles belonging to the family Cantharidae (Bouchard et al. 2011). To date, it has approximately 2000 species belonging to 43 genera (Yang 2010, Švihla 2011), which are widely distributed in the Holoarctic and Oriental regions (Brancucci 1980). Traditionally, the taxonomy of this group is mainly based on the structure of male genitalia and tarsal claws. However, it is impossible to accurately identify all species by only using these characters, especially for the morphologically similar sibling species, such as Falsopodabrushimalaicus species complex (Yang et al. in press). Moreover, it is not easy to clarify the status of some species among the related genera, such as Habronychus (Monohabronychus) multilimbatus (Pic, 1910), which was transferred several times (Okushima 2003, Švihla 2004, Brancucci 2007) in the Stenothemus genera complex (Švihla 2004). These difficulties underline the need for further studies to clarify the taxonomy of cantharid beetles either by searching for new morphological characters of high diagnostic value or applying alternative effective methods.
It is well-known that wing shape of insects exhibits a high heritability in nature (Bitner-Mathé and Klaczko 1999, Moraes et al. 2004), wing morphology is of a primary importance to entomologists interested in systematics. It was Comstock (1893) who first popularized the use of insect wing venation for traditional classification (Kunkel 2004). Since the 1970’s, several authors have begun to use the insect wings especially 2D morphometrical studies in systematics and phylogeny (Plowright and Stephen 1973, Rohlf 1993, Klingenberg 2003, Gumiel et al. 2003). Geometric morphometrics utilizes powerful and comprehensive statistical procedures to analyze shape differences of a morphological feature, using either homologous landmarks or outlines of the structure (Rohlf and Marcus 1993, Marcus and Corti 1996, Adam et al. 2004), and it is considered to be the most rigorous morphometric method (Gilchrist et al. 2000, Debat et al. 2003). Wings are excellent structure for studying morphological variation because they are basically 2-dimensional and the venation provides many well-defined morphological landmarks (Gumiel et al. 2003), the interactions of the veins, which are easy for identification and able to capture the general shape of the wing (Bookstein 1991). Among insects, the use of geometric morphometric analysis to study wing venation has been useful in identification at the individual level (Baylac et al. 2003, Dujardin et al. 2003, Sadeghi et al. 2009), in distinguishing sibling species (Matias et al. 2001, De la Riva et al. 2001, Villegas et al. 2002, Klingenberg and Savriama 2002, Roggero and Dentrèves 2005, Aytekin et al. 2007, Francuski et al. 2009, Tüzün 2009) and in delimitation among the genera (Baracchi et al. 2011). However, this modern effective methodology has not been applied in the studies of cantharid beetles until now.
In Cantharidae, the venation of hind wings was suggested to be of diagnostic value in the subfamily level based on the comparative morphology by Brancucci (1980). But within the subfamily, the variables of the veins are shown to be quantitative in metric properties, which can not be studied well by the traditional morphometrics, so it remains unknown whether the hind wing morphology contributes to the delimitation of genera or species or not. Thus in the present study, we apply the landmark-based geometric morphometric method to quantify and analyze wing morphological features in nine species belonging to three genera of Cantharinae, including Lycocerus Gorham, 1889 (sensu Okushima 2005, more than 300 species in the world), Prothemus Champion, 1926 (60 species in total), and Themus Motschulsky, 1838 (approximately 250 species worldwide), which are all mostly distributed in the Oriental region. The central aim of the study is to evaluate wing shape variation and test the possible use of wing shape patterns for generic or specific taxonomy of Cantharinae.
Material and methodSample collections
Hind wings of the following Cantharinae species (Table 1) are used in this study. Prior to geometric morphometric analysis, identification of specimens was performed using other morphological characters of adults (Yang 2010). The materials of the representative species are deposited in the Museum of Hebei University, Baoding, China (MHBU) and Institute of Zoology, Chinese Academy of Sciences, Beijing, China (IZAS) respectively. The left hind wing of each specimen (215 wings in total) was removed from the body and mounted in neutral balsam between a microscope slide and a cover slip. For each species, the chosen male and female specimens are subequal in number.
The number of specimens of each species used in the GM analysis.
The images of hind wings were captured using a stereomicroscope Nikon SMZ1500 and attached video camera Canon 450D connected to a HP computer. They were annotated using the TpsUtil software (Rohlf 2010a). The coordinates of the landmarks (13 landmarks in total, Table 2) were digitized by the TpsDig2.16 software (Rohlf 2010b) as shown in Fig. 1.
Hind wing of Lycocerusasperipennis showing digitizing landmarks.
https://binary.pensoft.net/fig/43028
Landmarks of hind wing (according to veins nomenclature system by Kukalová-Peck and Lawrence (1993).
No.
Junctions of veins
No.
Junctions of veins
1
ScP (Subcosta Posterior) and RA
8
MP_{1+2} and MP_{3+4}
2
RA (Radius Anterior) and RA_{3+4}
9
MP_{3+4} and CuA_{1} (Cubitus Anterior)
3
RA_{1+2} and RA_{3+4}
10
MP_{4} and MP_{3}
4
RA_{3+4} and r3 (radial crossvein)
11
CuA_{1} and CuA_{2}
5
RA_{3+4} and r4
12
CuA and CuA_{1+2}
6
r4 and RP (Radius Posterior)
13
AA (Anal Anterior) and CuA_{3+4}
7
RP and MP_{1+2} (Media Posterior)
Geometric morphometric analyses
To examine the wing shape variation, the digitized landmark data is analyzed using MorphoJ software (Klingenberg 2011). The variability in the shape space is assessed using a Principal Component Analysis (PCA). To better visualize the shape variation, we present the average configuration of landmarks for each genus or species. Deformation grids are used to portray the resulting shape variations.
The relative similarity and discrimination of the genera or species is analyzed using Canonical Variates Analysis (CVA). CVA finds shape values that maximize group means relative to variation within groups, by assuming that covariate matrices are identical (Klingenberg 2010). This is an effective method for detecting differences among taxa. The statistical significance of pairwise differences in mean shapes is determined using permutation tests (10 000 replications) with Procrustes and Mahalanobis distances. Both tests are used to assess significance because p-values can differ due to the anisotropy (direction dependency) of shape variation (Klingenberg and Monteiro 2005).
To evaluate the role of wing size in discrimination among different genera or species, the centorid size (CS) was compared. In the absence of allometry, the CS is the only size measure uncorrelated with all the shape variables (Bookstein 1991). The CS values are compared for genera and species respectively, because as a measurement of overall size variation of wings, they are far more sensitive than conventional measurements (Klingenberg et al. 1998). One-way analysis of variance (ANOVA) and Tukey HSD pairwise comparisons are employed to determine significant differences among genera or species. For visualizing size differences among groups, a 95% confidence intervals of the mean is computed using SPSS 13.0 and plotted in EXCEL.
Results
The shape variations of the hind wings in the genera Lycocerus, Prothemus and Themus is shown by the first two principal components of PCA (Fig. 2A). The thin plate spline visualizations show that the medial area (around by junctions Nos 9‒13) contributes most to the shape differences among the genera, especially the situation of the junction of MP_{4} and MP_{3} (No. 10) is most variable in Themus, while least in Lycocerus, and similar for the junction of ScP and RA (No. 1). Also, the junctions of r4 and RP (No. 6) and RP and MP_{1+2} (No. 7) appear more variable in Themus than in Lycocerus or Prothemus. Besides, the hind wing shape is more elongate in Themus than the other two genera. The centroid size (Fig. 6A, Table 7) is significantly different among the three genera (all p<0.05). The CVA scatterplot of shape differences for these genera (Fig. 2B) shows that each genus occupies different area. Mahalanobis distances among the three genera are significantly different in all pairwise comparisons (p<0.05), and Procrustes distances (p<0.05) are similar (Tables 3).
Shape variables of the hind wings in the genera of Lycocerus, Prothemus and Themus. A principal component analysis (PCA) of hind wing configuration. Plot of PC1 (74.39% of total variation) and PC2 (8.52% variation) showing 90% confidence ellipses of population means B canonical variate analysis (CVA) of same matrix, also showing 90% confidence ellipses of population means. The averaged shape of each genus is depicted as deformations using thin plate splines.
https://binary.pensoft.net/fig/43029
Difference in the hind wing shapes among the genera Lycocerus, Pothemus and Themus. Mahalanobis distances (left) & Procrustes distances (right): p-values (above); distances between populations (below).
Lycocerus
Pothemus
Themus
Lycocerus
Pothemus
Themus
Lycocerus
—
<.0001
<.0001
—
<.0001
<.0001
Pothemus
4.6396
—
<.0001
0.0456
—
<.0001
Themus
10.8932
10.446
—
0.1323
0.1088
—
In Lycocerus (Fig. 3A), the thin plate spline visualizations show that the junction of MP_{4} and MP_{3} (No. 10) is less variable in L.orientalis than in L.asperipennis or L.metallescens, and MP_{3+4} and CuA_{1} (No. 9) is more variable in L.asperipennis than the other two. In Prothemus (Fig. 4A), the junction of MP_{4} and MP_{3} (No. 10) is most variable in P.kiukiangensis, while least in P.purpuripennis, and AA and CuA_{3+4} (No. 13) is less variable in P.chinensis than the other two. In Themus (Fig. 5A), the junction of ScP and RA (No. 1) is most variable in T.licenti, while least in T.impressipennis. The centroid size (Fig. 6B, Table 7) is significantly different between L.asperipennis and L.orientalis (p=0.001) or L.metallescens (p=0.001), P.chinensis and P.kiukiangensis (p=0.005) or P.purpuripennis (p=0.002), but others are not (p>0.05). The CVA scatterplots of shape differences for each genus (Fig. 3B, 4B, 5B) all show that each species occupies different area. Mahalanobis distances among the three species of each genus are significantly different in all pairwise comparisons (p<0.05), and Procrustes distances are similar (p<0.05) (Tables 4, 5, 6).
Shape variables of the hind wings in the Lycocerus species. A principal component analysis (PCA) of hind wing configuration. Plot of PC1 (49.02% of total variation) and PC2 (14.92% variation) showing 90% confidence ellipses of population means B canonical variate analysis (CVA) of same matrix, also showing 90% confidence ellipses of population means. The averaged shape of each species is depicted as deformations using thin plate splines.
Shape variables of the hind wings in the Prothemus species. A principal component analysis (PCA) of hind wing configuration. Plot of PC1 (38.40% of total variation) and PC2 (15.88% variation) showing 90% confidence ellipses of population means B canonical variate analysis (CVA) of same matrix, also showing 90% confidence ellipses of population means. The averaged shape of each species is depicted as deformations using thin plate splines.
Shape variables of the hind wings in the Themus species. A principal component analysis (PCA) of hind wing configuration. Plot of PC1 (32.87% of total variation) and PC2 (16.48% variation) showing 90% confidence ellipses of population means B canonical variate analysis (CVA) of same matrix, also showing 90% confidence ellipses of population means. The averaged shape of each species is depicted as deformations using thin plate splines.
Comparisons of centroid size variables among different groups: ALycocerus, Prothemus and ThemusBLycocerusasperipennis, L.metallescens and L.orientalis; Prothemuschinensis, P.kiukiangensis and P.purpuripennis; Themuslicenti, T.coelestis and T.impressipennis.
https://binary.pensoft.net/fig/43033
Difference in the hind wing shapes among the species of Lycocerus. Mahalanobis distances (left) & Procrustes distances (right): p-values (above); distances between populations (below).
L.metallescens
L.asperipennis
L.orientalis
L.metallescens
L.asperipennis
L.orientalis
L.metallescens
—
<.0001
<.0001
—
<.0001
0.0466
L.asperipennis
5.6866
—
<.0001
0.0413
—
0.0003
L.orientalis
4.2970
4.4457
—
0.0182
0.0321
—
Difference in the hind wing shapes among the species of Prothemus. Mahalanobis distances (left) & Procrustes distances (right): p-values (above); distances between populations (below).
P.chinensis
P.kiukiangensis
P.purpuripennis
P.chinensis
P.kiukiangensis
P.purpuripennis
P.chinensis
—
<.0001
<.0001
—
<.0001
0.0002
P.kiukiangensis
5.7352
—
<.0001
0.0376
—
<.0001
P.purpuripennis
4.8174
5.5146
—
0.0247
0.0381
—
Difference in the hind wing shapes among the species of Themus. Mahalanobis distances (left) & Procrustes distances (right): p-values (above); distances between populations (below).
T.licenti
T.coelestis
T.impressipennis
T.licenti
T.coelestis
T.impressipennis
T.licenti
—
<.0001
<.0001
—
<.0001
<.0001
T.coelestis
6.7942
—
<.0001
0.0363
—
0.0001
T.impressipennis
6.8548
3.9959
—
0.0311
0.016
—
Discussion
The result of PCA shows that the shape differences of the hind wings among the genera Lycocerus, Prothemus and Themus (Fig. 2A) are mostly associated with the junctions of MP_{4} and MP_{3} (No. 10), ScP and RA (No. 1), r4 and RP (No. 6) and RP and MP_{1+2} (No. 7), and the shape of Themus is much more different from that of Lycocerus than Prothemus. And those variations within each genus (Figs 3A, 4A, 5A) appear in one or two junctions, which are either same to that of the genera or not, such as MP_{3+4} and CuA_{1} (No. 9) in Lycocerus and AA and CuA_{3+4} (No. 13) in Pothemus. This demonstrates that the shape differences among genera are much more variable than that within genus, and the variations among the species of each genus are different from one another.
The CVA results (Figs 2B, 3B, 4B, 5B) show that the three genera and the species of each genus are all successfully discriminated, since that Mahalanobis and Procrustes distances (Tables 3–6) for each group are significantly different (p<0.05). It suggests that the hind wing shape is useful for discrimination of both genus and species in Cantharinae by the geometric morphometrics. Also, the hind wing size is considered to be valuable in delineating the genera, but its role is uncertain for the species because of the inconsistent results in the three genera (Table 7).
Tukey HSD for the CS among different groups: p-values (above); mean differences (below). Asterisk (*) indicates the mean difference is significant at the 0.05 level.
CS among different genera
Lycocerus
Prothemus
Themus
Lycocerus
—
0.006
0
Prothemus
218.52316401(*)
—
0.001
Themus
483.54109456(*)
-265.01793055(*)
—
CS among the species of Lycocerus
L.asperipennis
L.metallescens
L.orientails
L.asperipennis
—
0.001
0.001
L.metallescens
474.67493257(*)
—
1
L.orientails
489.29359311(*)
14.61866054
—
CS among the species of Prothemus
P.chinensis
P.kiukiangensis
P.purpuripennis
P.chinensis
—
0.005
0.002
P.kiukiangensis
-456.74308033(*)
—
1
P.purpuripennis
-460.37428735(*)
-3.63E+00
—
CS among the species of Themus
T.coelestis
T.impressipennis
T.licenti
T.coelestis
—
0.711
0.998
T.impressipennis
-183.8607895
—
0.992
T.licenti
-79.25669086
104.6040987
—
Herein it can be concluded that the hind wing shape is useful for the discriminations of genera and species of Cantharinae. The geometric morphometrics represents a reliable tool not only in the taxonomic research but also in further study on the evolution of the hind wing shape of cantharid beetles.
Acknowledgment
We are grateful to Prof. Xingke Yang (IZAS) for his great support in providing the studied material. Thanks are due to the anonymous referee for their valuable comments on our manuscript and Mr. John MacDermott (USA) for correcting the English.
The present study was supported by the National Natural Science Foundation of China (Nos. 31172135, 41401064), the Knowledge Innovation of Chinese Academy of Sciences (Nos. KSCX2-EW-G-4, KSCX2-EW-Z-8), the Natural Science Foundation of Hebei Province (No. C2013201261) and the Foundation of the Key Laboratory of Zoological Systematics and Application of Hebei Province (No. 14967611D).
ReferencesAdamDCRohlfFJSliceDE (2004) Geometric Morphometrics: Ten Years of Progress following the “Revolution”.AytekinMATerzoMRasmontPçağatayN (2007) Landmark based geometric morphometric analysis of wing shape in Sibiricobombus Vogt (Hymenoptera: Apidae: Bombus Latreille).BaracchiDDapportoLTurillazziS (2011) Relevance of wing morphology in distinguishing and classifying genera and species of Stenogasterinae wasps.BaylacMVillemantCSimbolottiG (2003) Combing geometric morphometrics with pattern recognition for the investigation of species complex.Bitner-MathéBCKlaczkoLB (1999) Heritability, phenotypic and genetic correlations of size and shape of Drosophilamediopuncatata wings.BooksteinFL (1991) BouchardPBousquetYDaviesAEAlonso-ZarazagaMALawrenceJFLyalCHNewtonAFReidCASchmittMSlipińskiSASmithAB (2011) Family-group names in Coleoptera (Insecta).BrancucciM (1980) Morphologie comparee, evolution et systematique des Cantharidae (Insecta: Coleoptera).ComstockJH (1893) DebatVBéginMLegoutHDavidJR (2003) Allometric and nonallometric components of Drosophila wing shape respond differently to developmental temperature.De la RivaJLe PontFAliVMatiasAMollinedoSDujardinJP (2001) Wing geometry as a tool for studying the Lutzumyialongipalpis (Diptera: Pychodidae) complex.DujardinJPLe PontFBaylacM (2003) Geometric versus interspecific differentiation of sand flies: a landmark data analysis.FairmaireL (1886) Descriptions de coléoptères de l’intérieur de la Chine (2. partie).FairmaireL (1891) Coléoptères de l’intérieur de la Chine. (7 partie).FairmaireL (1900) Description de coléoptères nouveaux recueillis en Chine par M. de Latouche.FrancuskiLjVujićAKovačevićALudoškiJMilankovV (2009) Identification of the species of the Cheilosiavariabilis group (Diptera, Syrphidae) from the Balkan Peninsula using wing geometric morphometrics, with the revision of the status of C.melanoparedi Vujić, 1996.GilchristASAzevedoRBRPartridgeLO’HigginsP (2000) Adaption and constraint in the evolution of Drosophilamelanogaster wing shape.GorhamHS (1889) Descriptions of new speicies and a new genus of Coleoptera of the family Telephoridae.GumielMCataláSNoireauFde AriasARGarciaADujardinJP (2003) Wing geometry in Triatomainfestans (Klug) and T.melanosoma Martinez, Olmedo & Carcavallo (Hemiptera: Reduviidae).KlingenbergCP (2003) Developmental instability as a research tool: using patterns of fluctuating asymmetry to infer the developmental origins of morphological integration. In: PolakM (Ed.) KlingenbergCP (2010) Evolution and development of shape: integrating quantitative approaches.KlingenbergCP (2011) MorphoJ: an integrated software package for geometric morphometrics.KlingenbergCPMonteiroLR (2005) Distances and directions in multidimensional shape spaces: implications for morphometric applications.KlingenbergCPSavriamaY (2002) Geometric morphometrics of complex symmetric structures: Shape analysis of symmetry and asymmetry with Procrustes methods.KlingenbergCPMcintyreGSZaklanSD (1998) Left-right asymmetry of fly wings and the evolution of body axes.Kukalová-PeckJLawrenceJF (1993) Evolution of the hind wings in Coleoptera.KunkelJG (2004) MarcusLFCortiM (1996) Overview of the new, or geometric morphometrics. In: MarcusLFCortiMLoyANaylorGJPSliceDE (Eds) MatiasADe la RivaJTorrezMDujardinJP (2001) Rhodniusrobustus in Bolivia identified by its wings.MoraesEMSpressolaVLPradoPRRCostaLFSeneFM (2004) Divergence in wing morphology among sibling species of the Drosophilabuzzatii cluster.PicM (1938) Malacodermes exotiques.PlowrightRCStephenWP (1973) Evolutionary relationships in northern European Bombus and Psithyrus species (Apidae: Hymenoptera).RoggeroADentrèvesPP (2005) Geometric morphometric analysis of wing variation between two populations of the Scythrisobscurella species-group: geographic or interspecific differences? (Lepidoptera: Scythrididae).RohlfFJ (1993) Relative warp analysis and an example of its application to mosquito wings. In: MarcusLFBelloEGarcia-ValdecasasA (Eds) RohlfFJ (2010a) RohlfFJ (2010b) RohlfFJMarcusLF (1993) A revolution in morphometrics.SadeghiSAdriaensDDumontHJ (2009) Geometric morphometric analysis of wing shape variation in ten European populations of Calopteryxsplendens (Harris, 1782).ŠvihlaV (2004) New taxa of the subfamily Cantharinae (Coleoptera, Cantharidae) from southeastern Asia with notes on other species.ŠvihlaV (2011) New taxa of the subfamily Cantharinae (Coleoptera: Cantharidae )from south-eastern Asia, with notes on other species III.TüzünA (2009) Significance of wing morphometry in distinguishing some of the Hymenoptera species.VillegasJFeliciangeliMDDujardinJP (2002) Wing shape divergence between Rhodniusprolixus from Cojedes (Venezuela) and Rhodniusrobustus from Merida (Venezuela).YangYX (2010) YangYXSuJYYangXKBaiM (in press) New insight to solve the confusing Taxonomy problem using geometric morphometrics: an Example from a new species of Falsopodabrus Pic (Coleoptera, Cantharidae).