source: https://doi.org/10.7892/boris.89939 | downloaded: 2.10.2021 o:10.1111/zoj.12380 doi: Ichneumonidae) igemre,pr ftemtcodilcyto- mitochondrial the of part marker, single a Diplazon © [email protected] E-mail: author. *Corresponding human of the consequences of of (Pimm other introduction result activities and the a species, habitats, as invasive natural rates of going unprecedented destruction are species on at because biodiversity needed extinct the sorely increase describe is we an planet which our Such at speed 2007). the Monaghan, in & (e.g., Vogler studies 2002; numerous Butcher in 2004; dis- Blaxter, speed demonstrated species greatly facilitate as and to covery, workflow potential taxonomic the the up have data Molecular insect KEYWORDS: on ADDITIONAL phenomenon this of well-established impact between and introgression © prevalence DNA the mitochondrial of evolution. discuss and examples and systematics by few introgression, very facilitated of endosymbiont-mediated species, one is ( insect This gene related. closely identical nuclear share even fast-evolving species exist. examples a of empirical Wolbachia effect few genotyping the the only which explain examine but morphometrics, approach could phenomenon, we endosymbionts this barcoding study, bacterial discussed DNA present have by ( standard the reviews patterns 1 Numerous the In inheritance species. oxidase barcoding. new of mitochondrial of cytochrome for of performance failure rates gene, distortion limited description mitochondrial rate, and a single discovery substitution shown a both up on have speed relies studies to potential several the However, has delimitation species Molecular 2015 November 17 publication for accepted 2015; October 26 revised 2015; July 31 Received 2 1 KLOPFSTEIN SERAINA genus wasp parasitoid Wolbachia Society, Linnean the of Journal Zoological atesrse7 02Br,Switzerland Bern, 3012 7, Baltzerstrasse Switzerland Bern, CH-3005 iiino omnt clg,Isiueo clg n vlto,Uiest fBern, of University Evolution, and Ecology of Institute Ecology, 15, Community Bernastrasse of Bern, Division Burgergemeinde der Museum Naturhistorisches Invertebrates, of Department 06TeLnenSceyo London, of Society Linnean The 2016 06TeLnenSceyo London, of Society Linnean The 2016 lhuhitgaietxnm eonzs1 species, 16 recognizes taxonomy integrative Although . eso httefiueo N acdn onie ihtepeec fteedsmin.Two endosymbiont. the of presence the with coincides barcoding DNA of failure the that show we , INTRODUCTION tal. et noybot itr N acdn nthe in barcoding DNA distort endosymbionts 04.Fraias h s of use the animals, For 2014). , tal. et COI 02 .. Tautz e.g., 2012; , 1,2, Wolbachia yohoeoiae1–hbiiain–itorsin–seisdelimitation. species – introgression – hybridization – 1 oxidase cytochrome altpsand haplotypes ,CRSINKROPF CHRISTIAN *, olgclJunlo h ina Society Linnean the of Journal Zoological 2016, olgclJunlo h ina Society Linnean the of Journal Zoological erve iia eot,poieals fciei oidentify to criteria of list a provide reports, similar review We . 177 Diplazon 4–5.Wt figures 6 With 541–557. , tal. et Wolbachia Wolbachia , 0 Fn mad 03 Meier 2003; and Omland, 10% & approximately different between (Funk threshold-based in range 30% estimated species of to were of studies failures methods prevalence delimitation and the hypotheses accumulate; non-monophyly species to established started and barcoding DNA Schmidt esfcett eii h atmjrt fspecies of Derycke Consortium majority vast 2007; (Hajibabaei the Life delimit to examined reported to was sufficient of locus be ‘barcoding’ this studies, Barcode & earlier Hebert Ratnasingham 2003b; the deWaard, (Hebert, by (www.barcodeoflife.org) cated chrome tan,ee though even strains, atrao acdn nteprsti apgenus wasp parasitoid the in barcoding on bacteria on oee,rprso nossece between inconsistencies on reports however, Soon, 1,2 COI ITS2 COI (Hymenoptera: n ANSBAUR HANNES and .Bsdsicmlt ieg otn ralow a or sorting lineage incomplete Besides ). c n cenn the screening and ) tal. et nyrcvr pt e.Adn multivariate Adding ten. to up recovers only xds uui CI,hsbe advo- been has (COI), 1 subunit oxidase 2016, , 2015). , tal. et 177 2016, , ITS2 541–557 , 08 Smith 2008; , tal. et tal. et ugssta hyaenot are they that suggests 177 06 G 2006; , 03) seilyin Especially 2003a). , 1,2 Diplazon 541 : – 557 tal. et tal. et omez pce for species 2008; , 2006; , tal. et 541 , 542 S. KLOPFSTEIN ET AL.

Bergsten et al., 2012). The most commonly invoked number of species (Heath et al., 1999; Hurst & Jig- explanations for such discordances are inadequate gins, 2005; Raychoudhury et al., 2009); the success of taxonomy, rapid speciation rates relative to the sub- mtDNA-based species delimitation could thus be stitution rates of the marker, incomplete lineage dependent on the infection history of the species sorting, and hybrid introgression. The first issue, under consideration. inadequate taxonomy, either in the form of misiden- The impact of intracellular endosymbionts can be tifications or poorly supported species hypotheses, is manifold. In the best-case scenario with respect to especially important when unverified data are DNA barcoding, the endosymbionts remain within retrieved from databases, or in poorly studied and one host species (or are transferred between species highly diverse taxa (Collins & Cruickshank, 2012). only infectiously, i.e., not via hybridization; see This issue can only be resolved by careful taxonomic below). By manipulating the reproductive biology of work. The second issue, insufficient variability, can their hosts to maximize their own transmission, as be ameliorated by informed marker choice; variabil- demonstrated, for example, in the a-proteobacterium ity has to be sufficiently high to allow distinguishing Wolbachia, they cause a decrease in the intraspecific between very closely-related species, and some sta- mtDNA variability, in the extreme case causing the tistical approaches for species delimitation even spread of a single haplotype within a whole popula- require some variation at the intraspecific level tion or even species (Turelli, Hoffmann & McKech- (Pons et al., 2006; Zhang et al., 2013). Mitochondrial nie, 1992; Jiggins, 2003; Charlat et al., 2009). They DNA (mtDNA) is advantageous in this respect might thus contribute to the ‘barcoding gap’ (i.e., the because of comparatively high evolutionary rates difference between intraspecific and interspecific and a corresponding high per-base-pair information variation) and thus even improve the performance of content at or around the species boundary (Brown DNA barcoding. Endosymbionts only cause problems et al., 1982; Mindell & Thacker, 1996; Lin & Dan- if an infection is passed on between species through forth, 2004; Mueller, 2006; Simon et al., 2006; Zink a (potentially very rare) hybridization event; under & Barrowclough, 2008); on the other hand, nuclear such a scenario, the bacteria are transferred along- markers (especially introns) often still need to be side a foreign mitochondrium and, if the bacterium established for non-model taxa. The third issue, manages to successfully spread through the new host incomplete lineage sorting, leads to ancestral poly- species through vertical transmission, it will drag morphisms being retained in different, closely- the mtDNA with it. related species, resulting in discordances between The numerous reviews that discuss this mecha- gene trees and species trees (Knowles & Carstens, nism of endosymbiont-mediated mtDNA introgres- 2007; Rosenberg & Tao, 2008; Edwards, 2009; Pru- sion (Johnstone & Hurst, 1996; Jiggins, 2003; fer et al., 2012). Once more, mtDNA is at an advan- Ballard & Rand, 2005; Hurst & Jiggins, 2005; Gal- tage here compared to nuclear DNA because it will tier et al., 2009) draw on very few convincing empiri- attain complete lineage sorting and thus reciprocal cal examples. For a study to provide plausible species monophyly more quickly as a result of the evidence for the role of an endosymbiont in facilitat- four-fold smaller effective population size of the ing mtDNA introgression, it needs to include both mitochondrial compared to the nuclear genome. donor and recipient species and demonstrate a strict By contrast, the last issue (i.e., a higher suscepti- association of both their mtDNA and endosymbiont bility of mtDNA to hybrid introgression) has been strains. To our knowledge, there are currently only discussed as a major drawback of mtDNA (Ballard & six studies that fulfil these requirements (Ballard, Rand, 2005; Hurst & Jiggins, 2005; Bachtrog et al., 2000; Jiggins, 2003; Narita et al., 2006; Whitworth 2006; Rubinoff, Cameron & Will, 2006). After a et al., 2007; Gompert et al., 2008; Raychoudhury hybridization event, the recombining nuclear genes et al., 2009). are continuously eliminated with successive back- In the present study, we provide evidence for crossing, but recombination does not or only very endosymbiont-mediated mtDNA introgression from a rarely occur in mtDNA. Hybridization thus some- group of parasitoid wasps. The genus Diplazon has times leads to the complete introgression of a foreign recently been revised on a mainly morphological mitochondrial genome into another species (Paquin basis, with the discovery of four new species in Eu- & Hedin, 2004; Berthier, Excoffier & Ruedi, 2006; rope (Klopfstein, 2014). Twenty Western Palaearctic Edwards & Bensch, 2009; Galtier et al., 2009; Petit species are currently recognized. We analyzed sev- & Excoffier, 2009; Nicholls et al., 2012). Besides the eral specimens from 15 species from Switzerland and lack of recombination, selective sweeps caused by Sweden and one North American species with DNA hitchhiking with selfish genetic elements such as barcoding and found a very poor recovery rate of the endosymbiotic bacteria have been shown to greatly morphologically defined species. Parasitoid wasps influence the population biology of mtDNA in a might be even more prone to Wolbachia infections

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 543 than other insects because of the additional infection specimen did not entirely match the morphological pathway via their arthropod hosts (Cook & Butcher, species concepts in the revision, identifications were 1999; Heath et al., 1999). We thus screened the marked with a ‘cf.’. This was also the case for the wasps for Wolbachia and found that infections specimens of Diplazon tibiatorius with dark hind usually coincide with a failure of DNA barcoding and coxae, which are considered to belong to this species vice versa. To examine the potential role of despite this deviation from the common morphotype Wolbachia in distorting mtDNA inheritance through of D. tibiatorius with orange hind coxae (Klopfstein, hybrid introgression, we compiled a morphometrics 2014). Scientific names and morphological terminol- dataset to support the morphological species ogy are employed sensu Klopfstein (2014). concepts, studied three nuclear and one additional mitochondrial marker for the wasps, and performed MORPHOMETRIC ANALYSIS multilocus strain typing for the endosymbiotic bacteria. The morphological differentiation of Diplazon species is largely based on colour and sculpture of the integument (Klopfstein, 2014). To some extent, such MATERIAL AND METHODS features are subjective and depend on taxonomic expertise. We thus conducted a shape principle com- TAXON SAMPLING AND IDENTIFICATION ponent analysis (PCA) using an independent set of We identified 140 specimens of 15 European species measurements to establish how well the qualitative of Diplazon parasitoid wasps from Switzerland and morphological differentiation is supported by the Sweden (Table 1; see also Supporting information, measurements. We measured 14 characters covering Table S1) using the keys provided in a recent mor- most body parts (see Supporting information, Tables phological revision of the group (Klopfstein, 2014). A S2, S3). The selection corresponds to the measure- North American species, Diplazon bradleyi, was ments used for calculating some of the typically used added together with five outgroup species of the body ratios in Ichneumonidae (Townes, 1969; Klopf- same genus group as Diplazon (Klopfstein, Kropf & stein, 2014). Quicke, 2010; Klopfstein et al., 2011). Because of the We applied the multivariate ratio analysis (MRA) importance of the sculpture of the tergites, mesoscu- of Baur & Leuenberger (2011) to our data. MRA is a tum and mesopleuron for delimiting species in this relatively new approach that is especially suited for genus, we examined all specimens under the same analyzing body measurements in a taxonomy context lighting, a 23-W energy saving lamp. Whenever a because it offers several tools for the analysis of size

Table 1. Summary of taxon and gene sampling

Diplazon species Specimens* Countries† COI ITS2 28S EF1a NADH1

D. albotibialis 2f CH, SE 2 2 2 2 2 D. annulatus 4f, 1 m CH, SF 5 5 D. bradleyi 1f US 1 1 1 1 1 D. deletus 8f CH, SE 8 8 2 1 1 D. flixi 11f, 3 m CH 14 14 2 1 1 D. hyperboreus 2f SE 2 2 1 1 1 D. laetatorius 8f CH, ES, SE, US, ZM 8 8 4 2 2 D. orientalis 3f, 1 m TH 4 4 1 1 1 D. pallicoxa 3f CH, SE 3 3 1 1 1 D. parvus 9f CH, SE 9 6 D. pectoratorius 7f CH, SE, SF 7 7 2 1 1 D. scutatorius 11f CH, SE 11 11 2 2 2 D. tetragonus 19f, 3 m CH, SE 22 21 3 1 1 D. tibiatorius 5f CH, SE 5 5 4 2 2 D. varicoxa 21f, 6 m CH, SE 27 25 3 3 3 D. zetteli 9f, 4 m CH 13 12 1 1 1

*Number of female (f) and male (m) specimens included in the analysis. †Country of origin of the specimens. CH = Switzerland; ES = Spain; SE = Sweden; SF = Finland; TH = Thailand; US = United States of America; ZM = Zambia.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 544 S. KLOPFSTEIN ET AL. and shape in the multivariate geometrical frame- ENDOSYMBIONT SCREENING AND STRAIN TYPING work (Baur et al., 2014; Schweizer, Hertwig & See- hausen, 2014). Here, we computed a shape PCA and To test whether diplazontine wasps are infected with plotted the first shape PC against isometric size, endosymbionts, we performed a preliminary PCR defined as the geometric mean of all variables (see screening of 32 female specimens from a variety of Supporting information, Fig. S1). Graphic visualiza- species, including eight Diplazon species (Diplazon tion of the correlation of size with shape allowed us deletus, D. flixi, D. laetatorius, D. parvus, D. pectora- to estimate the amount of allometry in the data torius, D. tetragonus, D. varicoxa, and D. zetteli), for (Klingenberg, 1998) (i.e., to judge the impact of an Cardinium hertigii, Rickettsia sp., Spiroplasma indirect size effect on the separation of some taxa). ixodetis, S. poulsonii, and Wolbachia pipientis, using Morphometric analyses were calculated with R, the primers listed in Duron et al. (2008). We version 3.1.2 (R Core Team, 2014), using slightly obtained positive results only in twelve cases for modified R-scripts provided by Baur et al. (2014). Wolbachia, and none for any of the other endosym- Scatterplots were generated with the package ‘gg- bionts. We thus focussed on Wolbachia for further plot20 (Wickham, 2009). analyses within the genus Diplazon. We performed PCR screening on one to ten female specimens per species using primers that amplify a MOLECULAR METHODS fragment of the Wolbachia surface protein gene Genomic DNA was extracted from whole specimens (wsp)(wsp81F and wsp691R; (Braig et al., 1998). preserved in 80% ethanol using the Promega Wizard These primers have been used successfully in the kit for blood and tissue extractions. Vouchers and past to test for Wolbachia infections in Hymenotpera DNA samples are kept at the Natural History (Beukeboom & Pijnacker, 2000) and showed a good Museum in Bern and the Swedish Museum of Natu- performance, especially for supergroup A and B Wol- ral History in Stockholm (Table 1; see also Support- bachia, the two groups previously found in ing information, Table S1). Approximately 700 bp hymenopterans (Simoes~ et al., 2011). As a control for from the 50 end of the mitochondrial COI were ampli- the quality of the extracted DNA, we amplified the fied using the primers LCO and HCO designed by ribosomal 28S gene, which is specific to eukaryotes, Folmer et al. (1994). To obtain approximately 800 bp alongside the Wolbachia screening, using the for- of the nuclear ribosomal RNA (rRNA) internal tran- ward primer designed by Belshaw & Quicke (1997) scribed spacer 2 (ITS2), we used the primers and the reverse primer from Mardulyn & Whitfield designed by Quicke et al. (2006). Three additional (1999) and the PCR conditions described above. Spec- markers, part of the nuclear 28S rRNA (28S), the F2 imens tested positive for Wolbachia in the first round copy of elongation factor 1-a (EF1a) (Klopfstein & were used as positive controls in a second round of Ronquist, 2013), and the mitochondrial NADH 1 screening. gene (ND1), were taken from previous studies (see At least one infected individual per species was Supporting information, Table S1) (Klopfstein et al., then typed using a multilocus sequence typing 2010, 2011). approach relying on five house-keeping genes (Baldo Polymerase chain reactions (PCR) were carried out et al., 2006), complemented by the faster-evolving in 20-ll final volumes using Promega GoTaq Flexi wsp gene; the general primers were used as listed in DNA Polymerase kits. Final volumes contained Baldo et al. (2006) (http://pubmlst.org/wolbachia).

30 pmol of MgCl2, 16 pmol of both primers, 4 pmol of Because multiple infections were detected in several each dNTP, 0.3 U Taq polymerase and 2 ll of geno- species, we used molecular cloning in Escherichia mic DNA. PCR conditions were: 94 °C for 2 min, 35 coli with the Topo Ta Cloning kit (Life Technologies) cycles of 30 s at 94 °C, 30 s at the respective anneal- to separate the different products of the wsp gene ing temperature (51 °C for COI and 49 °C for ITS2), in some specimens (Table 2). Four clones were and 1 min at 72 °C, followed by a final extension at sequenced per specimen. Because multilocus 72 °C for 10 min. PCR products were either purified sequence typing is not possible for multiple infections with the GFXTM DNA and Gel Purification kit without making assumptions about the prevalence of (Amersham Biosciences) or by the purification ser- different allele combinations at the five loci (Arthofer vice of Macrogen Korea. The PCR products were et al., 2011), we only resolved the full sequence typ- sequenced on an ABI 377 automated sequencer using ing for the single infections. Nevertheless, we stud- Big Dye Terminator technology (Applied Biosystems). ied some of the double-peaked sequences of the All new sequences (112 new sequences of COI and multiple infections and found that they were in part 134 of ITS2) have been deposited in GenBank under in accordance with multilocus sequence typing accession numbers KR230498 to KR230743 (see (MLST) strain types already found as single infec- Supporting information, Table S1). tions in other specimens (Table 2); however, as a

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 545 result of the highly speculative nature of these asso- (Stamatakis, 2014), with the GTR model and parti- ciations, we only submitted the strain information tioning as above. Clade support was assessed by from single infections to the Wolbachia MLST data- 1000 bootstrap replications. base and to Genbank (accession numbers KR230444 Single-gene trees were also obtained for three to KR230477 for the MLST loci and KR230478 to additional markers (28S, EF1a, ND1) but for a KR230497 for the wsp gene). reduced set of specimens (see Supporting informa- tion, Table S1). For the single-gene analyses, we ran MrBayes as specified before but with two indepen- ALIGNMENT AND CALCULATION OF PAIRWISE dent runs only. To obtain a good estimate of the spe- DISTANCES cies tree needed to test for a correlation between The sequences of the protein-coding COI were aligned Wolbachia infections and a failure of barcoding in a with MUSCLE (Edgar, 2004) after translation into phylogenetic contest, a concatenated analysis of all amino acids using MEGA, version 6.06 (Tamura five genes was performed using a single sequence et al., 2013). Alignment was straightforward because sampled randomly for each species. Settings were no indels were detected. In ITS2, the alignment chosen as for the COI and the ITS2 analyses, with posed more problems because of numerous indels of the protein-coding genes partitioned into combined varying length. We thus only aligned the ingroup first and second versus third positions. All analyses taxa and rooted the tree according to previous multi- were run on the University of Bern Linux Cluster gene phylogenetic analyses (Klopfstein et al., 2010, UBELIX. 2011). Alignments can be downloaded from Tree- BASE (http://purl.org/phylo/treebase/phylows/study/ CORRELATION ANALYSIS TB2:S17676). Uncorrected pairwise distances (p distances) were To determine whether the Wolbachia infection status calculated in MEGA with pairwise deletion. Plots showed any correlation with the performance of the were produced in R, version 3.1.2 (R Core Team, mitochondrial COI in species delimitation, we scored 2014). To measure the performance of barcoding, we each species as being infected or not and as being used the threshold approach because it does not recovered in the most sensitive analysis with the require strict reciprocal monophyly of all the species. threshold set at 0.5% uncorrected p distance. Because We tried thresholds of 2%, 1%, and 0.5% uncorrected the Wolbachia infections did not appear to be indepen- p distance. dent from the phylogeny, we had to use a comparative statistical method that corrects for the phylogenetic relationships. We used BayesDiscrete from the Bayes- PHYLOGENETIC ANALYSIS Traits, version 2.0 (Pagel & Meade, 2006) under an Phylogenetic analyses of the COI and ITS2 genes independent and a dependent model of evolution. We were performed under a Bayesian and a maximum run both a ML and a Bayesian approach. For ML, we likelihood (ML) approach. For the Bayesian analyses, used the best-scoring tree found during the ML search we used MrBayes version 3.2.2 (Ronquist et al., on the five-gene dataset in RaxML as a test phy- 2012) with a ‘mixed’ substitution model (integrating logeny. Likelihoods obtained under the two models over the space of possible submodels of the general with 1000 ML attempts per tree were compared by a time-reversible model; Huelsenbeck, Larget & Alfaro, likelihood-ratio test. For the Bayesian approach, we 2004) and gamma-distributed among-site rate varia- evenly sampled 1000 trees from the post-burn-in tion, including a proportion of invariant sites. COI phase of the four independent runs in MrBayes to was run with the combined first and second codon represent the posterior distribution of topologies. The position in one and the third codon position in a sec- dependent and independent models were compared ond partition, ITS2 was unpartitioned. The four using a Bayes factor test, running both in standard independent Markov chain Monte Carlo (MCMC) mode (i.e., assuming a fixed number of four rates runs with one cold and three heated chains were run under the independent and eight under the dependent for 10 000 000 generations and sampled every model) and in the reversible-jump mode (which tries 1000th generation. As a conservative burn-in, we all combinations of equal and zero rates). used half of the samples. Convergence was good as judged from the average standard deviation of split frequencies (ASDSF) for the topology parameter and RESULTS the potential scale reduction factor (PSRF) for the COI BARCODING scalar parameters (ASDSF = 0.0064 for COI and 0.0043 for ITS2, PSRF < 1.002 for all parameters). The sequences of the barcoding portion of COI For the ML analyses, we used RAXML version 8.1 (50 portion of the gene) of 140 specimens of Diplazon

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 Table 2. Wolbachia infection status, barcoding success, and Wolbachia MLST and wsp profiles in Diplazon and related species 546

Wolbachia Individuals Success

† ‡ § ¶ KLOPFSTEIN S. Species infection inf/test* barcoding Ind # ST CoxA FbpA FtsZ GatB HpcA wsp HVR1 HVR2 HVR3 HVR4

Diplazon None 0/2 0 albotibialis Diplazon Multi 1/2 0 annulatus Diplazon Multi 1/1 1 ? bradleyi © Diplazon Single 1/4 0 706 92 59 17 3 54 68 442 51 55 15 25 TAL. ET 06TeLnenSceyo London, of Society Linnean The 2016 deletus Diplazon Single 1/3 0 832 92 59 17 3 54 68 688** 51 269** 15 25 flixi Diplazon None 0/2 1 hyperboreus Diplazon None 0/4 1 laetatorius Diplazon None 0/3 1 orientalis Diplazon None 0/3 1 pallicoxa Diplazon Triple or 3/5 0 625 23/338/? 1/1/17 12/12/9 21/21/? 19/269/18 parvus more Diplazon None 0/3 0 pectoratorius Diplazon None 0/3 1

olgclJunlo h ina Society Linnean the of Journal Zoological scutatorius Diplazon Single 1/3 0 843 432** 1 408** 3 88 257** 23 1 12 21 19 tetragonus Diplazon Double or 2/2 0 1B2, 1B4 432** 1 408** 3 88 257** 23/? 1/1 12/12 21/? 19/19 tibiatorius more Diplazon Triple or 10/10 0 825 (1/7) (120/249?) (3/6) (88/8?) (185?) 23/? 1/1 12/56 21/15 19/272 varicoxa more Diplazon None 0/4 1 zetteli Campocraspedon None 0/1 N/A caudatus Syrphophilus Single 1/1 N/A 514 433** 84 120 200** 234** 257** 689** 11 9 267** 302** asperatus Tymmophorus Single 2/2 N/A 818 52 201** 54 41 75 11 9 15 25 obscuripes Tymmophorus None 0/1 N/A suspiciosus Xestopelta None 0/1 N/A gracilima 2016, ,

*Number of infected/tested individuals. † 177 Success of DNA barcoding: species delimitation unsuccessful (0) or successfull (1) at the most sensitive threshold of 0.5%. ‡

541–557 , Specimen identifier. § Wolbachia strain type from multilocus sequence typing (MLST), based on the five loci CoxA, FbpA, FtsZ, GatB, and HpcA. ¶ Wolbachia surface protein gene (wsp) allele number, and typing of its four hypervariable regions (HVR). **New alleles found in the present study. ?No sequence obtained. WOLBACHIA DISTORTS BARCODING IN DIPLAZON 547 showed surprisingly little divergence between species sequenced a nuclear gene that is known to evolve at of which some were described more than 100 years a relatively fast rate, the rRNA spacer ITS2. The ago and whose status has subsequently remained resulting gene tree (Fig. 5) corresponds well to COI unchallenged in the taxonomic literature. Of the 15 concerning the deeper nodes but shows a very differ- species included with more than one specimen, nine ent pattern for some of the species. Diplazon deletus were not recovered as monophyletic on the majority- and D. flixi do not appear as closely related species rule consensus tree obtained from the Bayesian anal- in ITS2 but, instead, the former clusters with a ysis (Fig. 1). Even more notable, several species North American species, D. bradleyi. There is very shared identical haplotypes: D. deletus with D. flixi, high support for this grouping and, in addition to the D. annulatus with D. tetragonus, and D. parvus with information in the nucleotide sequence, there are D. varicoxa. three indels in the ITS2 sequence that support this Because gene-tree monophyly is not a necessary relationship (of length 1, 2, and 4 bp, respectively). condition for certain species-delimitation methods to These indels are not present in D. flixi or any other work, we used the threshold method as promoted by species grouped with D. deletus in the COI analyses. many barcoding proponents to assess the extent of Diplazon flixi now clusters with specimens of the two discordance between morphologically defined species morphologically very similar species: Diplazon hyper- and those recovered by DNA barcoding. The species boreus and D. zetteli. Neither D. parvus, nor D. vari- recovered when using the threshold method are coxa are recovered by ITS2 as monophyletic, shown in Figure 2. At 2% uncorrected p distance, a although they appear as clearly separated on the value often used for species delimitation in insects current gene tree, with D. parvus now sharing iden- (Hebert et al., 2004), only six species are recovered, tical sequences with D. tibiatorius. This grouping eight at 1%, and ten at the very low value of 0.5%. again is in better accordance with the morphology of Figure 3 shows that the intra- and interspecific the species (Klopfstein, 2014), and suggests that distances overlap broadly in several species pairs. these are good species after all. The species annula- The use of a relative threshold as sometimes pro- tus and tetragonus, however, are not recovered by posed (e.g., ten times the intraspecific variability) ITS2 either but, instead, show a pattern similar to (Hebert et al., 2004) would thus not have improved the COI tree. Finally, several species that are mono- the situation. We here use the distance to the closest phyletic on the COI tree are not supported here (e.g., other species (i.e., the minimum interspecific dis- the aforementioned D. tibiatorius and the species tance) because reporting the average instead of the pair D. scutatorius and D. orientalis, which share smallest interspecific distances exaggerates the bar- identical ITS2 sequences). The two specimens of coding gap (Meier, Zhang & Ali, 2008). D. hyperboreus do not cluster together; one of them showed some deviations from the morphological diag- nosis of the species (see below). MORPHOMETRY The gene trees obtained from the single-gene The morphometric analyses supported most of the analyses of two additional nuclear markers (28S which species hypotheses derived from discrete morphologi- is in close proximity to ITS2 in the genome, and EF1a; cal characters. For all analyses, only the first shape see Supporting information, Fig. S2) show that these PC was informative, which was then plotted markers both evolve too slowly to contain much infor- against isometric size. A scatterplot confined to the mation about the species in question; although the three species pairs of special interest in the present D. albotibialis and D. pectoratorius pair and the group study (D. deletus–D. flixi, D. annulatus–D. tetragonus, around D. laetatorius are recovered, there is little or D. parvus–D. varicoxa) (Fig. 4; see also Supporting no resolution among the questionable species. Not sur- information, Fig. S1) revealed that some of them can- prisingly, the second mitochondrial gene sequenced, not be separated based on quantitative morphology. ND1, largely confirms the picture recovered in the For example, D. annulatus is nested within D. te- COI tree but with less resolution. Diplazon deletus tragonus, and D. parvus within D. varicoxa. How- and D. flixi once more have identical sequences, as ever, D. deletus is clearly distinct from D. flixi by the have D. bradleyi and D. varicoxa (D. parvus and first shape PC (Fig. 4). The two species are also D. annulatus were not sequenced for this locus). entirely overlapping in size; hence, the shape differ- ence cannot be attributed to an allometric size effect. WOLBACHIA INFECTIONS AND CORRELATION ANALYSIS Twenty out of 54 Diplazon individuals or eight of the NUCLEAR ITS2 16 species tested positive for Wolbachia infections, As morphometry did not provide sufficient support (Table 2). To test whether the infections coincided for all the morphological species hypotheses, we with a failure in DNA barcoding, we used a method

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 548 S. KLOPFSTEIN ET AL.

Xestopelta gracillima 3C3 Tymmophorus obscuripes 1H2 Campocraspedon caudatus 1C2 Syrphophilus asperatus 514 Tymmophorus suspiciosus 2B11 D. albotibialis 901 0.53/61 D. albotibialis 502 0.84/59 D. pectoratorius 1C7 1.0/99 D. pectoratorius 520 1.0/79 D. pectoratorius 1C6 0.63/38 D. pectoratorius 3H2 D. pectoratorius 4F3 0.86/- D. pectoratorius 815 D. pectoratorius 907 D. laetatorius 1C4 D. laetatorius 1C5 D. laetatorius 3A2 1.0/85 D. laetatorius 3B2 0.99/98 D. orientalis 707 D. laetatorius 4B10 D. orientalis 701 1.0/95 D. laetatorius 4G1 D. orientalis 711 D. laetatorius 712 D. orientalis 704 D. laetatorius 703 D. scutatorius 1A10 D. scutatorius 1A11 1.0/100 D. scutatorius 2A3 D. scutatorius 2E12 D. scutatorius 2E3 D. scutatorius 4B12 1.0/99 1.0/95 D. scutatorius 4D4 D. scutatorius 4E3 0.01 substitutions per site D. scutatorius 4F2 D. pallicoxa 503 D. scutatorius 709 D. pallicoxa 601 D. cf. scutatorius 845 1.0/98 D. pallicoxa 618 D. annulatus 1C1 D. tetragonus 2A9 D. annulatus 0.76/- D. annulatus 4G6 D. annulatus 2D6 D. tetragonus 4F6 D. annulatus 2D4 1.0/- D. annulatus 2D7 D. deletus 1A7 D. tibiatorius 1B2 D. deletus 1A8 1.0/98 D. tibiatorius 1B3 D. flixi 4H1 D. cf. tibiatorius 1B4 D. deletus 2E9 D. cf. tibiatorius 3E1 D. deletus 4H10 D. cf. tibiatorius 3F1 D. flixi 1A1 D. zetteli 1A3 D. deletus 714 D. zetteli 1A4 D. flixi 4D12 D. zetteli 2D2 D. deletus 706 D. tetragonus D. zetteli 2D3 D. flixi 1A2 D. zetteli 2D5 D. flixi 2E5 1.0/98 D. zetteli 2E1 D. flixi 2E6 D. zetteli 2E2 D. deletus 2E10 0.97/63 D. zetteli 3A1 D. flixi 4G11 D. zetteli 3B1 0.65/- D. flixi 4G12 D. zetteli 3B4 D. flixi 4G7 D. zetteli 3C1 D. flixi 4G9 D. zetteli 621 D. deletus 4H12 D. zetteli 622 D. flixi 4H3 D. hyperboreus 903 D. flixi 4H5 1.0/68 D. cf. hypberboreus 619 D. flixi 4H8 D. flixi 832 D. bradleyi 2A5 D. tetragonus 1A12 D. parvus 623 D. tetragonus 1B1 D. varicoxa 620 D. tetragonus 2A10 D. parvus 625 D. tetragonus 2A4 D. deletus D. varicoxa 909 D. tetragonus 2A7 D. varicoxa 848 0.58/49 D. tetragonus 2A8 D. parvus 629 D. tetragonus 3B3 D. varicoxa 616 D. tetragonus 3C2 D. cf. parvus 905 0.62/- D. tetragonus 3E2 D. cf. parvus 913 D. parvus D. tetragonus 3G2 D. varicoxa 1A6 D. tetragonus 4B7 D. varicoxa 2E11 D. tetragonus 4E11 D. varicoxa 2E7 D. tetragonus 4G3 D. varicoxa 2E8 D. tetragonus 809 D. varicoxa 3D1 D. tetragonus 843 D. varicoxa 4D3 D. tetragonus 911 D. varicoxa 521 D. tetragonus 3A3 D. parvus 624 D. tetragonus 3D2 1.0/90 D. varicoxa 630 D. tetragonus 3E3 D. varicoxa 631 D. tetragonus 3F2 D. varicoxa 632 D. varicoxa 825 D. flixi D. varicoxa 846 D. parvus 627 D. varicoxa D. varicoxa 921 D. varicoxa 1A5 D. varicoxa 1A9 D. varicoxa 613 D. varicoxa 614 D. varicoxa 615 D. parvus 844 D. varicoxa 633 D. parvus 626 D. cf. varicoxa 902 D. cf. varicoxa 923 D. varicoxa 2F1 D. cf. varicoxa 617

Figure 1. Bayesian majority-rule consensus tree as retrieved from the barcoding fragment of COI mtDNA. Support val- ues close to the nodes represent Bayesian posterior probabilites and bootstrap support based on 1000 replicates. Inlaid photographs show specimens of some of the unresolved species. Part of the tree was cut at the triangle and moved to the left to fit on a single page.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 549 us orius bore eyi tus vus entalis ragonus

Threshold le i utat e r lixi hyper bradl par zetteli tibiatorius d f pallicoxa albotibialis pectoratorius o laetatorius sc annulatus tet varicoxa 2% 1 2 3 4 5 6

1% 1 2 3 4 5 6 7 8

0.5% 1 2 3 4 5 6 7 8 9 10

Figure 2. Numbers of species recovered and identity of lumped species as obtained by the threshold method for three different threshold values. Distances are uncorrected pairwise distances in the CO1 barcoding locus.

Intra-specific distances inter-specific distances (to closest species) orientalis albotibialis pectoratorius albotibialis tetragonus pectoratorius tibiatorius orientalis scutatorius deletus annulatus tetragonus annulatus hyperboreus tetragonus flixi deletus deletus varicoxa parvus zetteli varicoxa orientalis parvus laetatorius flixi scutatorius pallicoxa 0.00 0.01 0.02 0.03 0.04 0.05 0.06 flixi zetteli parvus deletus varicoxa pallicoxa orientalis tibiatorius annulatus laetatorius albotibialis tetragonus scutatorius pectoratorius

Figure 3. Intra- and interspecific uncorrected p distances for each species. The minimum interspecific distances are shown (i.e., those to the closest species included in our dataset; species identities are indicated above the respective box- plot). Grey bars indicate species that do not show any barcoding gap (i.e., for which intra- and interspecific distances overlap). correcting for phylogenetic relationships by compar- dard MCMC = 5.14, from reversible-jump MCMC = ing an evolutionary model assuming independent 2.84; Bayes factors are considered significant from a with one based on dependent evolution. Barcoding value of 2 and highly signficiant from a value of 5; was considered successful if the species were recov- Kass & Raftery, 1995]. ered at the most sensitive threshold of 0.5%. Both To examine the extent to which this method was the ML approach and the Bayesian approach sensitive to taxon sampling, we also analyzed the significantly preferred the dependent over the inde- dataset under the assumption that the aberrant pendent model of evolution in both ML model testing specimen of D. hyperboreus, which was retrieved (likelihood ratio test statistic = 10.29, P < 0.0358) apart from the specimen with the typical morphology and Bayesian testing [Bayes factor as 2 9 (difference in the ITS2 tree, actually represents a different in the logarithm of the harmonic mean) from stan- species. This change had a drastic impact on the

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 550 S. KLOPFSTEIN ET AL.

Species the USA (Wolbachia MLST database, http:// deletus pubmlst.org/wolbachia; accessed 20 May 2015). Dipla- flixi zon tetragonus has two previously unknown alleles for the genes FbpA and HpcA and thus harbours a new strain. Interestingly, the wsp allele found in this spe- cies has previously been found in species as diverse as 0.2 the parasitic wasp Nasonia longicornis (Darling) (Hymenoptera, Pteromalidae), the fruit fly Rhagoletis cerasi (Linnaeus) (Diptera, Tephritidae), and the vine- 0.1 gar fly Leucophenga maculosa (Coquillett) (Diptera, Drosophilidae) (http://pubmlst.org/wolbachia). Dipla- zon tibiatorius carries multiple though closely related 0.0 wsp alleles, although MLST sequence typing only

PC1 (50.9 %) recovered a single strain, which is identical to the one found in D. tetragonus. Finally, the infections of –0.1 D. varicoxa and D. parvus could not be typed by the MLST approach because of multiple infections (triple or more), although one of the cloned wsp sequences

–0.2 corresponds to the same allele (#23) as the infections in D. tetragonus and D. tibiatorius (Table 2). The phy- –0.3 –0.2 –0.1 0.0 0.1 0.2 logeny of the wsp sequences of the single-infected spe- Isometric size cies and the successfully cloned multiple infections (Fig. 6) confirms the strain typing results. Figure 4. Scatterplot of isometric size versus the first shape principle component of the species pair Diplazon deletus–Diplazon flixi. PC, principal component. DISCUSSION

FAILURE OF DNA BARCODING IN DIPLAZON result, which thus has to be considered with caution. With this new dataset, the significance disappeared, We found a rather poor performance of the standard with the likelihood ratio test statistic dropping to DNA barcoding approach in delimiting species in the 5.47 and the P value increasing to 0.24. parasitoid wasp genus Diplazon. Using a 2% sequence divergence threshold as often advocated (Hebert et al., 2004), only six of the 16 species could WOLBACHIA DIVERSITY be recovered, making this approach clearly insuffi- The Wolbachia surface protein gene (wsp)sequences cient in this genus. This result is somewhat in con- indicated single infections in D. deletus, D. flixi, trast to other studies that have reported good D. tetragonus, and the outgroup species Syrphophilus success of barcoding in other groups of parasitoids, asperatus and Tymmophorus obscuripes. Multiple where COI typically recovered many more species infections, as evident from polymorphic peaks and than morphology (Smith et al., 2008; Stigenberg & often also length variation in the wsp sequences, were Ronquist, 2011; Butcher et al., 2012). However, mor- detected for D. annulatus, D. bradleyi, D. parvus, phological examinations were often not very detailed D. tibiatorius,andD. varicoxa. For the latter three, in the past and most studies covered only limited we performed molecular cloning of the wsp gene to geographical regions. Further studies will show assess how many infections were present and whether whether the poor performance in Diplazon is just an they were similar to single infections already detected exception for the group. A simple DNA-based method in the present study (Table 2). for species delimitation would be urgently needed in Multilocus sequence typing of the singly-infected parasitoid wasps, given their enormous and highly species recovered one known and three unknown Wol- understudied diversity (Quicke, 2012). On the other bachia strains (Table 1). Diplazon deletus and D. flixi hand, parasitoids might have extreme population are both infected with the same strain #92 (http:// dynamics because of their high trophic level, show pubmlst.org/wolbachia). The sequences of the fast- unusually biased sex-ratios and even high levels of evolving wsp gene differs only by a single mutation in inbreeding in some groups, and potentially have very the second hypervariable region of the gene, which is fast speciation rates through host switching and further evidence for the very close relationship ecological speciation (Feder & Forbes, 2010; Konig€ between those Wolbachia infections. The same strain et al., 2015). Furthermore, parasitoids might be more has also been found in two lepidopteran species from prone to Wolbachia infections as a result of their

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 551

D. albotibialis 502 0.75/51 D. albotibialis 901 D. pectoratorius 1C7 D. laetatorius 3B2 1.0/100 D. pectoratorius 520 D. laetatorius 1C4 D. pectoratorius 4F3 D. laetatorius 1C5 D. pectoratorius 1C6 0.95/76 D. laetatorius 3A2 D. pectoratorius 3H2 D. laetatorius 4B10 0.51/- D. pectoratorius 815 1.0/100 D. laetatorius 4G1 D. pectoratorius 907 D. laetatorius 712 D. laetatorius 703 D. orientalis 711 D. scutatorius 2E3 1.0/100 D. cf. scutatorius 845 D. scutatorius 709 D. scutatorius 1A10 D. scutatorius 1A11 D. scutatorius 2A3 1.0/100 D. scutatorius 2E12 1.0/96 D. orientalis 701 D. scutatorius 4B12 D. scutatorius 4D4 D. scutatorius 4E3 D. scutatorius 4F2 D. orientalis 707 D. orientalis 704 D. pallicoxa 503 1.0/100 D. pallicoxa 601 D. bradleyi 2A5 D. pallicoxa 618 D. deletus 1A7 1.0/100 D. deletus 1A8 1.0/100 D. deletus 2E10 D. deletus 2E9 1.0/89 D. deletus 4H10 D. deletus 4H12 D. deletus 714 D. deletus 706 D. tetragonus 2A8 0.99/67 D. annulatus 1C1 D. annulatus 2D4 D. annulatus 2D6 D. annulatus 2D7 1.0/86 D. annulatus 4G6 D. tetragonus 2A9 D. tetragonus 4E11 1.0/100 D. tetragonus 4F6 0.005 substitutions per site D. tetragonus 1A12 D. tetragonus 1B1 D. tetragonus 2A10 D. tetragonus 2A4 D. tetragonus 2A7 D. tetragonus 3A3 D. varicoxa 1A5 0.76/53 D. tetragonus 3B3 D. cf. varicoxa 617 D. tetragonus 3C2 D. tetragonus 3D2 D. varicoxa 2E11 D. tetragonus 3E2 D. cf. hyperboreus 619 D. tetragonus 3E3 D. varicoxa 2F1 D. tetragonus 3F2 D. tetragonus 3G2 D. tibiatorius 1B2 D. tetragonus 4B7 D. parvus 624 D. tetragonus 4G3 D. parvus 627 D. tetragonus 8 09 D. cf. tibiatorius 1B4 1.0/87 D. tetragonus 8 43 D. parvus 625 D. flixi 1A1 1.0/86 D. parvus 629 D. flixi 1A2 D. tibiatorius 1B3 D. flixi 2E5 1.0/96 D. parvus 623 D. flixi 2E6 D. parvus 626 D. flixi 4D12 D. cf. tibiatorius 3E1 D. hyperboreus 903 cf. D. tibiatorius 3F1 0.55/- D. flixi 4G12 D. varicoxa 1A6 D. flixi 4G7 D. varicoxa 2E7 D. flixi 4G9 D. varicoxa 4D3 D. flixi 4H1 D. varicoxa 521 D. flixi 4H3 D. varicoxa 631 D. flixi 4H5 D. varicoxa 620 D. varicoxa 630 D. flixi 4H8 D. varicoxa 632 D. flixi 832 D. varicoxa 633 D. flixi 4G11 D. zetteli 2E1 1.0/77 D. varicoxa 846 D. varicoxa 848 1.0/92 D. zetteli 622 D. varicoxa 825 D. zetteli 621 D. varicoxa 921 D. zetteli 3C1 D. varicoxa 1A9 D. zetteli 3B1 D. varicoxa 2E8 D. zetteli 3A1 D. varicoxa 3D1 D. zetteli 2E2 D. varicoxa 613 D. zetteli 2D5 D. varicoxa 614 D. varicoxa 615 0.94/61 D. zetteli 2D3 D. varicoxa 616 D. zetteli 2D2 D. varicoxa 909 D. zetteli 1A4 D. zetteli 1A3

Figure 5. Bayesian majority-rule consensus tree as retrieved from the internal transcribed spacer 2 rRNA. The Dipla- zon species are shown in different colours. Support values close to the nodes represent Bayesian posterior probabilites and the bootstrap support based on 1000 replicates. The branch leading to Diplazon albotibialis and Diplazon pectorato- rius has been shortened to fit on a single page.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 552 S. KLOPFSTEIN ET AL.

Syrphophilus asperatus 514 1.0 Tymmophorus obscuripes 1H2 0.75 Diplazon parvus 624 Clone 1 706 1.0 Diplazon deletus Diplazon flixi 1A2 0.98 Diplazon flixi 832 1.0 Diplazon varicoxa 825 Clone 1 Diplazon parvus 624 Clone 2 1.0 Diplazon parvus 624 Clone 3 Diplazon tetragonus 843 1.0 Diplazon_tibiatorius 1B2 Clone 1 Diplazon tibiatorius 1B2 Clone 2 1.0 Diplazon cf. tibiatorius 1B4 Clone 1 Diplazon cf. tibiatorius 1B4 Clone 2 Diplazon varicoxa 825 Clone 2 0.01 substitutions per site Diplazon varicoxa 825 Clone 3

Figure 6. Bayesian majority-rule consensus tree of the Wolbachia surface protein (wsp) sequences of isolates from nine species of Diplazontinae parasitic wasps. Strains separated by molecular cloning in Escherichia coli were given arbitrary numbers. Values next to nodes represent Bayesian posterior probabilities. intimate relationships with their hosts that might barcoding. However, this result should be treated act as a potent transmission pathway (Cook & with caution because it was very sensitive to the spe- Butcher, 1999; Vavre et al., 1999). All of these fac- cies hypotheses. This was exemplified by our exercise tors might complicate the population biology of their of assuming species status for the aberrant individ- mtDNA and thus impede DNA barcoding. ual of D. hyperboreus, which was sufficient to annihi- In Diplazon, only a combination of the two mark- late our result. The test is probably not very stable ers COI and ITS2 recovered most of the morpholog- because of the very small number of taxa sampled ically defined species; in some cases, for example, (eight infected versus eight non-infected putative the recently described D. parvus, species status is species) and very short branch-lengths in the crown only supported by the combined information because group of the tree. Furthermore, the infection status this species was polyphyletic in both markers but of a species might have changed recently or not be with respect to two different species (Figs 1,4). detected correctly in our few-specimen assay. Given the limited use of COI and ITS2 as markers A recent critique of phylogenetic comparative for species delimitation, the establishment of addi- methods in general highlighted a basic shortcoming tional markers is necessary. A recent bioinformatics in that they often even retrieve a significant correla- approach using comparative genomics (Hartig et al., tion if the character histories involve only a single 2012) has already established a plethora of candi- origin (Maddison & FitzJohn, 2015). Currently, there date loci for the order Hymenoptera, some of which is no way to resolve this issue except for the careful contain introns that might provide sufficient vari- interpretation of such a result. In our case, a single ability to resolve questions at the species level, and origin of both characters can be excluded from the could be analyzed in multispecies coalescent phylogenetic distribution, although the species with approaches that make use of the information inher- both a failure of barcoding and Wolbachia infections ent in independently segregating markers to iden- are certainly concentrated among the crown group of tify reproductively isolated units (Yang & Rannala, the tree. The failure of barcoding might be tightly 2010). In any case, the failure of a single-marker linked to lowered evolutionary rates or rapid specia- identification system in this genus suggests that tion rates, and a phylogenetic component is likely for caution is necessary when using DNA taxonomy in both. By contrast, the observed Wolbachia infections parasitic wasps. certainly do not go back to a single infection in an ancestral species because some strains are very diver- gent and several species harbour multiple strains. CORRELATION BETWEEN WOLBACHIA INFECTIONS AND Furthermore, the retention of Wolbachia over such A FAILURE OF BARCODING time scales and that numbers of species boundaries is We found a significant correlation between the highly unlikely for this endosymbiont, which nor- Wolbachia infection status and the failure of mally shows much faster infection dynamics (Werren,

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 553

Baldo & Clark, 2008; Raychoudhury et al., 2009). cordance between nuclear and mitochondrial gene Even when accepting the correlation as true, there trees in the vinegar fly genus Drosophila (Diptera, might still be other causal explanations than Drosophilidae) with an introgression event from endosymbiont-mediated hybrid introgression (e.g., a Drosphila simulans Sturtevant to Drosphila mauri- role of Wolbachia in increasing speciation rates, tiana Tsacas and David, and such a transfer could which would at the same time decrease the success of even be repeated experimentally (Aubert & Solignac, DNA barcoding) (Werren et al., 2008; Raychoudhury 1990). Jiggins (2003) found a shared mtDNA haplo- et al., 2009). We thus need additional evidence to type in those individuals of two species of Acraea support the introgression scenario. butterflies (Lepidoptera, Nymphalidae) that were infected with the same Wolbachia strain but not in the uninfected individuals, and such a pattern is best A STRONG CASE FOR WOLBACHIA-MEDIATED explained by endosymbiont-mediated hybrid intro- MTDNA TRANSFER gression. Narita et al. (2006) recovered the same Additional evidence for Wolbachia-mediated mtDNA pattern in two recently discovered sibling species of introgression comes from several sources. Under a the butterfly genus Eurema (Lepidoptera, Pieridae). scenario of hybrid introgression, we expect the follow- Whitworth et al. (2007) examined twelve species of ing patterns: (1) very low mtDNA diversity indicative the blowfly genus Protocalliphora (Diptera, Cal- of a recent selective sweep; (2) mtDNA haplotypes liphoridae) and found that four species shared COI that are much more similar than likely given the spe- haplotypes and Wolbachia strains as judged from the cies relationships; (3) identical or at least very similar wsp gene, whereas AFLP markers suggested that Wolbachia infections; and (4) a likely opportunity for these species were not closely related. Similar results hybridization (Hurst & Jiggins, 2005). were obtained by Gompert et al. (2008) in Lycaeides For the species pair D. deletus–D. flixi, without any butterflies (Lepidoptera, Lycaenidae). Finally, Ray- doubt two good biological species that were also choudhury et al. (2009) examined the Wolbachia clearly distinct in our morphometric analysis (Fig. 4), infections in the parasitoid wasp genus Nasonia all four of the above points are fulfilled. Both COI and (Hymenoptera, Pteromalidae) and found the likely ND1 show identical haplotypes between the species co-transmission of Wolbachia and mtDNA from and the Swiss and Swedish populations of D. deletus. Nasonia giraulti Darling to Nasonia oneida There is strong evidence in the nuclear ITS2 marker Raychoudhury & Desjardins. Unusually large that D. deletus is more closely related to the North intraspecific mtDNA variation that coincides with American D. bradleyi; the sharing of an mtDNA hap- infections by different Wolbachia strains has been lotype in D. deletus and D. flixi through incomplete found in several studies (Ballard, Chernoff & James, lineage sorting is thus highly unlikely. The single 2002; Marshall, 2004; Riegler et al., 2005; Charlat Wolbachia infections in both species are of the same et al., 2009; Atyame et al., 2011; Xiao et al., 2012; strain type and only differ by a single mutation in the Ritter et al., 2013) but, without a clear hybridization highly variable wsp gene, an observation in accor- scenario including the species of origin, it might be dance with transmission through hybridization. that the bacterial strains became associated with Finally, D. deletus and D. flixi have been collected in divergent mtDNA haplotypes long after a transfer the same Malaise traps in the Swiss Alps which sug- event (or incomplete lineage sorting). gests geographical and phenological opportunity for hybridization. The case is similar for the species CONCLUSIONS D. parvus and D. varicoxa that share identical COI sequences both in Switzerland and Sweden, whereas Reports of endosymbiont-mediated mtDNA introgres- D. parvus clusters with D. tibiatorius in the ITS2 sion are very rare, with to our knowledge only six tree. Their multiple Wolbachia infections appear to convincing cases currently found in the literature. It partly overlap as well (because both species carry is difficult to estimate the prevalence and thus multiple infections, the strains could not be fully importance of this phenomenon. A recent study by typed). Wolbachia-mediated mtDNA transfer is thus the Consortium for the Barcode of Life (Smith et al., the likely cause for the failure of barcoding in these 2012) addressed Wolbachia and barcoding but cases, but insufficient variability or incomplete lin- mainly focussed on whether the CoxA gene (the bac- eage sorting probably also played a role (e.g., in the terial counterpart of COI) of the endosymbionts had sister species D. albotibialis and D. pectoratorius, been amplified instead of the sequences of the host which do not harbour Wolbachia infections). insects in recent barcoding initiatives; no conclusive There are very few other convincing examples for data was provided about hybrid introgression. endosymbiont-mediated mtDNA introgression in the Recent estimates of the prevalence of Wolbachia literature. Ballard (2000) could best explain the dis- infections on the one hand (20–65% of all insect spe-

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 554 S. KLOPFSTEIN ET AL. cies; Hilgenboeker et al., 2008; Werren & Windsor, Tettelin H, Werren JH. 2006. Multilocus sequence typing 2000) and of hybridization on the other (approxi- system for the endosymbiont Wolbachia pipientis. Applied mately 10% of animal species; Mallet, 2005) suggest and Environmental Microbiology 72: 7098–7110. that endosymbiont-mediated mtDNA introgression Ballard JWO. 2000. When one is not enough: introgression might not be as rare as previously assumed. The of mitochondrial DNA in Drosophila. Molecular Biology and scarcity of empirical studies could simply be the Evolution 17: 1126–1130. result of an inherent difficulty to confirm this sce- Ballard JWO, Rand DM. 2005. The population biology of nario; most studies that include Wolbachia bacteria mitochondrial DNA and its phylogenetic implications. are conducted within species or they are part of a Annual Review of Ecology, Evolution, and Systematics 36: 621–642. broad screening study that does not assess mtDNA Ballard JWO, Chernoff B, James AC. 2002. Divergence of patterns in the hosts. Financial considerations and mitochondrial DNA is not corroborated by nuclear DNA, favourable reports on the success of DNA barcoding morphology, or behavior in Drosophila simulans. Evolution have led to many biodiversity studies including only 56: 527–545. a single mtDNA marker and consequently not Baur H, Leuenberger C. 2011. Analysis of ratiois in multi- detecting any mito-nuclear discordances. Decreasing variate morphometry. Systematic Biology 60: 813–825. sequencing costs, the establishment of additional Baur H, Kranz-Baltensperger Y, Cruaud A, Rasplus J- nuclear markers with sufficient variability, and the Y, Timokhov AV, Gokhman VE. 2014. Morphometric broad realization of the limitation of a single-marker analysis and taxonomic revision of Anisopteromalus species delimitation system will improve this situa- Ruschka (Hymenoptera: Chalcidoidea: Pteromalidae) - an tion and lead to a more precise assessment of the integrative approach. Systematic Entomology 39: 691–709. prevalence of this phenomenon. Belshaw R, Quicke DLJ. 1997. A molecular phylogeny of the Aphidiinae (Hymenoptera: Braconidae). Molecular Phy- logenetics and Evolution 7: 281–293. ACKNOWLEDGEMENTS Bergsten J, Bilton DT, Fujisawa T, Elliott M, Mon- We would like to thank the Swedish Malaise Trap aghan MT, Balke M, Hendrich L, Geijer J, Herrmann Project (Oland,€ Sweden), Victoria Spinas (Alp Flix, J, Foster GN, Ribera I, Nilsson AN, Barraclough TG, Switzerland), and Fridli Marti (Obersand, Switzer- Vogler AP. 2012. The effect of geographical scale of sam- land) for providing specimens. Madelaine Geiger and pling on DNA barcoding. Systematic Biology 61: 851–869. Salome Steiner performed the morphometric mea- Berthier P, Excoffier L, Ruedi M. 2006. Recurrent surements. This research received support from the replacement of mtDNA and cryptic hybridization between Swiss National Science Foundation (PA00P3_145377) two sibling bat species Myotis myotis and Myotis blythii. Proceedings of the Royal Society of London Series B-Biologi- and the SYNTHESYS Project (http://www.synthesys. cal Sciences 273: 3101–3109. info) (grant SE-TAF-663), which is financed by Euro- Beukeboom LW, Pijnacker LP. 2000. Automictic partheno- pean Community Research Infrastructure Action genesis in the parasitoid Venturia canescens (Hymenoptera: under the FP6 ‘Structuring the European Research Ichneumonidae) revisited. Genome 43: 939–944. Area’ Programme. Blaxter ML. 2004. The promise of a DNA taxonomy. Philo- sophical Transactions of the Royal Society of London Series REFERENCES B-Biological Sciences 359: 669–679. Braig HR, Zhou W, Dobson SL, O’Neill SL. 1998. Cloning Arthofer W, Riegler M, Schuler H, Schneider D, Moder and characterization of a gene encoding the major surface K, Miller WJ, Stauffer C. 2011. Allele intersection analy- protein of the bacterial endosymbiont Wolbachia pipientis. sis: a novel tool for multi locus sequence assignment in Journal of Bacteriology 180: 2373–2378. multiply infected hosts. PLoS ONE 6: e22198. Brown WM, Pragera EM, Wang A, Wilson AC. 1982. Atyame CM, Delsuc F, Pasteur N, Weill M, Duron O. Mitochondrial DNA sequences in primates: tempo and mode 2011. Diversification of Wolbachia endosymbiont in the of evolution. Journal of Molecular Evolution 18: 225–239. Culex pipiens mosquito. Molecular Biology and Evolution Butcher BA, Smith MA, Sharkey MJ, Quicke DLJ. 2012. 28: 2761–2772. A turbo-taxonomic study of Thai Aleiodes (Aleiodes) and Aubert J, Solignac M. 1990. Experimental evidence for Aleiodes (Arcaleiodes) (Hymenoptera: Braconidae: Rogadi- mitochondrial DNA introgression between Drosophila spe- nae) based largely on COI barcoded specimens, with rapid cies. Evolution 44: 1272–1282. descriptions of 179 new species. Zootaxa 3457: 1–232. Bachtrog D, Thornton K, Clark A, Andolfatto P. 2006. Charlat S, Duplouy A, Hornett EA, Dyson EA, Davies Extensive introgression of mitochondrial DNA relative to N, Roderick GK, Wedell N, Hurst GDD. 2009. The joint nuclear genes in the Drosophila yakuba species group. Evo- evolutionary histories of Wolbachia and mitochondria in lution 60: 292–302. Hypolimnas bolina. BMC Evolutionary Biology 9: 64. Baldo L, Hotopp JCD, Jolley KA, Bordenstein SR, Collins RA, Cruickshank RH. 2012. The seven deadly sins Biber SA, Choudhury RR, Hayashi C, Maiden MC, of DNA barcoding. Molecular Ecology Resources 13: 969–975.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 555

Cook JM, Butcher RDJ. 1999. The transmission and ings of the Royal Society of London Series B-Biological effects of Wolbachia bacteria in parasitoids. Researches on Sciences 270: 313–321. Population Ecology 41: 15–28. Hebert PDN, Ratnasingham S, deWaard JR. 2003b. Derycke S, Fonseca G, Vierstraete A, Vanfleteren J, Barcoding animal life: cytochrome c oxidase subunit 1 Vincx M, Moens T. 2008. Disentangling taxonomy within divergences among closely related species. Proceedings of the Rhabditis (Pellioditis) marina (Nematoda, Rhabditidae) the Royal Society of London Series B-Biological Sciences species complex using molecular and morhological tools. 270: S96–S99. Zoological Journal of the Linnean Society 152: 1–15. Hebert PDN, Stoeckle MY, Zemlak TS, Francis CM. Duron O, Bouchon D, Boutin S, Bellamy L, Zhou L, 2004. Identification of birds through DNA barcodes. PLoS Engelstadter€ J, Hurst GDD. 2008. The diversity of repro- Biology 2: 1657–1663. ductive parasites among arthropods: Wolbachia do not walk Hilgenboeker K, Hammerstein P, Schlattmann P, Tel- alone. BMC Biology 6: 27. schow A, Werren JH. 2008. How many species are Edgar RC. 2004. MUSCLE: multiple sequence alignment infected with Wolbachia? - a statistical analysis of current with high accuracy and high throughput. Nucleic Acids data. FEMS Microbiology Letters 281: 215–220. Research 32: 1792–1797. Huelsenbeck JP, Larget B, Alfaro ME. 2004. Bayesian Edwards SV. 2009. Is a new and general theory of molecu- phylogenetic model selection using reversible-jump Markov lar systematics emerging? Evolution 63: 1–19. chain Monte Carlo. Molecular Biology and Evolution 21: Edwards SV, Bensch S. 2009. Looking forwards or looking 1123–1133. backwards in avian phylogeography? A comment on Zink Hurst GDD, Jiggins FM. 2005. Problems with mitochon- and Barrowclough 2008. Molecular Ecology 18: 2930–2933. drial DNA as a marker in population, phylogeographic and Feder JL, Forbes AA. 2010. Sequential speciation and the phylogenetic studies: the effects of inherited symbionts. diversity of parasitic insects. Ecological Entomology 35: 67– Proceedings of the Royal Society of London Series B-Biologi- 76. cal Sciences 272: 1525–1534. Folmer O, Black M, Hoeh W, Lutz R, Vrijenhoek R. Jiggins FM. 2003. Male-killing Wolbachia and mitochon- 1994. DNA primers for amplification of mitochondrial drial DNA: selective sweeps, hybrid introgression and para- cytochrome c oxidase subunit I from diverse metazoan site population dynamics. Genetics 164: 5–12. invertebrates. Molecular marine biology and biotechnology Johnstone RA, Hurst GDD. 1996. Maternally inherited 3: 294–299. male-killing microorganisms may confound interpretation Funk DJ, Omland KE. 2003. Species-level paraphyly and of mitochondrial DNA variability. Biological Journal of the polyphyly: frequency, causes, and consequences, with Linnean Society 58: 453–470. insights from animal mitochondrial DNA. Annual Review of Kass RE, Raftery AE. 1995. Bayes factors. Journal of Ecology, Evolution, and Systematics 34: 397–423. American Statistical Association 90: 773–795. Galtier N, Nabholz B, Glemin S, Hurst GDD. 2009. Klingenberg CP. 1998. Heterochrony and allometry: the Mitochondrial DNA as a marker of molecular diversity: a analysis of evolutionary change in ontogeny. Biological reappraisal. Molecular Ecology 18: 4541–4550. Reviews 73: 79–123. Gomez A, Wright PJ, Lunt DH, Cancino JM, Carvalho Klopfstein S. 2014. Revision of the Western Palaearctic GR, Hughes RN. 2007. Mating trials validate the use of Diplazontinae (Hymenoptera, Ichneumonidae). Zootaxa DNA barcoding to reveal cryptic speciation of a marine bry- 3801: 1–143. ozoan taxon. Proceedings of the Royal Society of London Klopfstein S, Ronquist F. 2013. Convergent intron gains in Series B-Biological Sciences 274: 199–207. hymenopteran elongation factor-1a. Molecular Phylogenetics Gompert Z, Forister ML, Fordyce JA, Nice CC. 2008. and Evolution 67: 266–276. Widespread mito-nuclear discordance with evidence for Klopfstein S, Kropf C, Quicke DLJ. 2010. An evaluation introgressive hybridization and selective sweeps in Lycaei- of phylogenetic informativeness profiles and the molecular des. Molecular Ecology 17: 5231–5244. phylogeny of Diplazontinae (Hymenoptera, Ichneumonidae). Hajibabaei M, Janzen DH, Burns JM, Hallwachs W, Systematic Biology 59: 226–241. Hebert PDN. 2006. DNA barcodes distinguish species of Klopfstein S, Quicke DLJ, Kropf C, Frick H. 2011. tropical Lepidoptera. Proceedings of the National Academy Molecular and morphological phylogeny of Diplazontinae of Sciences of the United States of America 103: 968–971. (Hymenoptera, Ichneumonidae). Zoologica Scripta 40: 379– Hartig G, Peters RS, Borner J, Etzbauer C, Misof B, 402. Niehuis O. 2012. Oligonucleotide primers for targeted Knowles LL, Carstens BC. 2007. Delimiting species with- amplification of single-copy nuclear genes in Apocritan out monophyletic gene trees. Systematic Biology 56: 887– Hymenoptera. PLoS ONE 7: e39826. 895. Heath BD, Butcher RDJ, Whitfield WGF, Hubbard SF. Konig€ K, Krimmer E, Brose S, Ganter C, Buschluter€ I, 1999. Horizontal transfer of Wolbachia between phylogenet- Konig€ C, Klopfstein S, Wendt I, Baur H, Krogmann L, ically distant insect species by a naturally occurring mecha- Steidle JLM. 2015. Does early learning drive ecological nism. Current Biology 9: 313–316. divergence during speciation processes in parasitoid wasps? Hebert PDN, Cywinska A, Ball SL, deWaard JR. 2003a. Proceedings of the Royal Society of London Series B-Biologi- Biological identifications through DNA barcodes. Proceed- cal Sciences 282: 20141850.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 556 S. KLOPFSTEIN ET AL.

Lin C-P, Danforth BN. 2004. How do insect nuclear and taxonomy of undescribed insects. Systematic Biology 55: mitochondrial gene substitution patterns differ? Insights 595–609. from Bayesian analyses of combined datasets. Molecular Prufer K, Munch K, Hellmann I, Akagi K, Miller JR, Phylogenetics and Evolution 30: 686–702. Walenz B, Koren S, Sutton G, Kodira C, Winer R, Maddison WP, FitzJohn RG. 2015. The unsolved challenge Knight JR, Mullikin JC, Meader SJ, Ponting CP, Lun- to phylogenetic correlation tests for categorical characters. ter G, Higashino S, Hobolth A, Dutheil J, Karakoc E, Systematic Biology 64: 127–136. Alkan C, Sajjadian S, Catacchio CR, Ventura M, Mar- Mallet J. 2005. Hybridization as an invation of the genome. ques-Bonet T, Eichler EE, Andre C, Atencia R, Trends in Ecology & Evolution 20: 229–237. Mugisha L, Junhold J, Patterson N, Siebauer M, Mardulyn P, Whitfield JB. 1999. Phylogenetic signal in Good JM, Fischer A, Ptak SE, Lachmann M, Symer the COI, 16S, and 28S genes for inferring relationships DE, Mailund T, Schierup MH, Andres AM, Kelso J, among genera of Microgastrinae (Hymenoptera; Bra- Paabo S. 2012. The bonobo genome compared with the conidae): evidence of a high diversification rate in this chimpanzee and human genomes. Nature 486: 527–531. group of parasitoids. Molecular Phylogenetics and Evolution Quicke DLJ. 2012. We know too little about parasitoid wasp 12: 282–294. distributions to draw any conclusions about latitudinal Marshall JL. 2004. The Allonemobius-Wolbachia host-endo- trends in species richness, body size and biology. PLoS symbiont system: evidence for rapid speciation and against ONE 7: e32101. reproductive isolation driven by cytoplasmic incompatibil- Quicke DLJ, Mori M, Zaldivar-Riveron A, Laurenne ity. Evolution 58: 2409–2425. NM, Shaw MR. 2006. Suspended mummies in Aleiodes Meier R, Shiyang K, Vaidya G, Ng PKL. 2006. DNA bar- species (Hymenoptera: Braconidae: Rogadinae) with coding and taxonomy in Diptera: a tale of high intraspecific descriptions of six new species from western Uganda based variability and low identification success. Systematic Biol- largely on DNA sequence data. Journal of Natural History ogy 55: 715–728. 40: 2663–2680. Meier R, Zhang G, Ali F. 2008. The use of mean instead of R Core Team. 2014. R: A language and environment for sta- smallest interspecific distance exaggerates the size of the tistical computing. R Foundation for Statistical Computing, “barcoding gap” and leads to misidentification. Systematic Vienna, Austria. Biology 57: 809–813. Raychoudhury R, Baldo L, Oliveira DCSG, Werren JH. Mindell DP, Thacker CE. 1996. Rates of molecular evolu- 2009. Modes of acquisition of Wolbachia: horizontal trans- tion: phylogenetic issues and applications. Annual Review fer, hybrid introgression, and codivergence in the Nasonia of Ecology and Systematics 27: 279–303. species complex. Evolution 63: 165–183. Mueller RL. 2006. Evolutionary rates, divergence dates, and Riegler M, Sidhu M, Miller WJ, O’Neill SL. 2005. Evi- the performance of mitochondrial genes in Bayesian phylo- dence for a global Wolbachia replacement in Drosophila genetic analysis. Systematic Biology 55: 289–300. melanogaster. Current Biology 15: 1428–1433. Narita S, Nomura M, Kato Y, Fukatsu T. 2006. Genetic Ritter S, Michalski SG, Settele J, Wiemers M, Fric ZF, structure of sibling butterfly species affected by Wolbachia Sielezniew M, Sasic M, Rozier Y, Durka W. 2013. Wol- infection sweep: evolutionary and biogeographical implica- bachia infections mimic cryptic speciation in two parasitic tions. Molecular Ecology 15: 1095–1108. butterfly species, Phengaris teleius and P. nausithous (Lepi- Nicholls JA, Challis RJ, Mutun S, Stone GN. 2012. Mito- dopters: Lycaenidae). PLoS ONE 8: e78107. chondrial barcodes are diagnostic of shared refugia but not Ronquist F, Teslenko M, Van der Mark P, Ayres DL, species in hybridiying oak gallwasps. Molecular Ecology 21: Darling A, Hohna€ S, Larget B, Liu L, Suchard MA, 4051–4062. Huelsenbeck JP. 2012. MrBayes 3.2: efficient Bayesian Pagel M, Meade A. 2006. Bayesian analysis of correlated phylogenetic inference and model choice across a large evolution of discrete characters by reversible-jump Markov model space. Systematic Biology 61: 539–542. chain Monte Carlo. American Naturalist 167: 808–825. Rosenberg NA, Tao R. 2008. Discordance of species trees Paquin P, Hedin M. 2004. The power and perils of ‘molecu- with their most likely gene trees: the case of five taxa. Sys- lar taxonomy’: a case study of eyeless and endangered tematic Biology 57: 131–140. Cicurina (Araneae: Dictynidae) from Texas caves. Molecu- Rubinoff D, Cameron S, Will K. 2006. A genomic perspec- lar Ecology 13: 3239–3255. tive on the shortcomings of mitochondrial DNA for ‘‘barcod- Petit RJ, Excoffier L. 2009. Gene flow and species delimita- ing’’ identification. Journal of Heredity 97: 581–594. tion. Trends in Ecology & Evolution 24: 386–393. Schmidt S, Schmid-Egger C, Moriniere J, Haszprunar Pimm SL, Jenkins CN, Abell R, Brooks TM, Gittleman G, Hebert PDN. 2015. DNA barcoding largely supports JL, Joppa LN, Raven PH, Roberts CM, Sexton JO. 250 years of classical taxonomy: identications for Central 2014. The biodiversity of species and their rates of European bees (Hymenoptera, Apoidea partim). Molecular extinction, distribution, and protection. Science 344: Ecology Resources 15: 985–1000. 1246752. Schweizer M, Hertwig S, Seehausen O. 2014. Diversity Pons J, Barraclough TG, Gomez-Zurita J, Cardoso A, versus disparity and the role of ecological opportunity in a Duran DP, Hazell S, Kamoun S, Sumlin WD, Vogler continental bird radiation. Journal of Biogeography 41: AP. 2006. Sequence-based species delimitation for the DNA 1301–1312.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 WOLBACHIA DISTORTS BARCODING IN DIPLAZON 557

Simoes~ PM, Mialdea G, Reiss D, Sagot M-F, Charlat S. Turelli M, Hoffmann AA, McKechnie SW. 1992. Dynam- 2011. Wolbachia detection: an assessment of standard PCR ics of cytoplasmic incompatibility and mtDNA variation in protocols. Molecular Ecology Resources 11: 567–572. natural Drosophila simulans populations. Genetics 132: Simon C, Buckley TR, Frati F, Stewart JB, Beckenbach 713–723. AT. 2006. Incorporating molecular evolution into phyloge- Vavre F, Fleury F, Lepetit D, Fouillet P, Bouletreau M. netic analysis, and a new compilation of conserved poly- 1999. Phylogenetic evidence for horizontal transmission of merase chain reaction primers for animal mitochondrial Wolbachia in host-parasitoid associations. Molecular Biol- DNA. Annual Review of Ecology, Evolution, and Systemat- ogy and Evolution 16: 1711–1723. ics 37: 545–579. Vogler AP, Monaghan MT. 2007. Recent advances in DNA Smith MA, Rodriguez JJ, Whitfield JB, Deans AR, Jan- taxonomy. Journal of Zoological Systematics and Evolution- zen DH, Hallwachs W, Hebert PDN. 2008. Extreme ary Research 45: 1–10. diversity of tropical parasitoid wasps exposed by iterative Werren JH, Windsor DM. 2000. Wolbachia infection fre- integration of natural history, DNA barcoding, morphology, quencies in insects: evidence of a global equilibrium? Pro- and collections. Proceedings of the National Academy of ceedings of the Royal Society of London Series B-Biological Sciences of the United States of America 105: 12359–12364. Sciences 267: 1277–1285. Smith MA, Bertrand C, Crosby K, Eveleigh ES, Fernan- Werren JH, Baldo L, Clark ME. 2008. Wolbachia: master dez-Triana J, Fisher BL, Gibbs J, Hajibabaei M, Hall- manipulators of invertebrate biology. Nature Reviews - wachs W, Hind K, Hrcek J, Huang D-W, Janda M, Microbiology 6: 741–751. Janzen DH, Li Y, Miller SE, Packer L, Quicke DLJ, Whitworth TL, Dawson RD, Magalon H, Baudry E. Ratnasingham S, Rodriguez JJ, Rougerie R, Shaw 2007. DNA barcoding cannot reliably identify species of the MR, Sheffield C, Stahlhut JK, Steinke D, Whitfield blowfly genus Protocalliphora (Diptera: Calliphoridae). Pro- JB, Wood M, Zhou X. 2012. Wolbachia and DNA barcod- ceedings of the Royal Society of London Series B-Biological ing insects: patterns, potential, and problems. PLoS ONE 7: Sciences 274: 1731–1739. e36514. Wickham H. 2009. ggplot2 - elegant graphics for data analy- Stamatakis A. 2014. RAxML version 8: a tool for phyloge- sis. Dordrecht: Springer. netic analysis and post-analysis of large phylogenies. Bioin- Xiao J-H, Wang N-X, Murphy RW, Cook JM, Jia L-Y, formatics 30: 1312–1313. Huang D-W. 2012. Wolbachia infection and dramatic Stigenberg J, Ronquist F. 2011. Revision of the Western intraspecific mitochondrial DNA divergence in a fig wasp. Palearctic Meteorini (Hymenoptera, Braconidae), with a Evolution 66: 1907–1916. molecular characterization of hidden Fennoscandian species Yang Z, Rannala B. 2010. Bayesian species delimitation diversity. Zootaxa 3084: 1–95. using multilocus sequence data. Proceedings of the National Tamura K, Stecher G, Peterson D, Filipsky A, Kumar S. Academy of Sciences of the United States of America 107: 2013. MEGA6: molecular evolutionary genetics analysis 9264–9269. version 6.0. Molecular Biology and Evolution 30: 2725– Zhang J, Kapli P, Pavlidis P, Stamatakis A. 2013. A gen- 2729. eral species delimitation method with applications to phylo- Tautz D, Arctander P, Minelli A, Thomas RH, Vogler genetic placements. Bioinformatics 29: 2869–2876. AP. 2002. DNA points the way ahead in taxonomy. Nature Zink RM, Barrowclough GF. 2008. Mitochondrial DNA 418: 479. under siege in avian phylogeography. Molecular Ecology 17: Townes HK. 1969. The genera of Ichneumonidae, Part 1. 2107–2121. Memoirs of the American Entomological Institute 11: 1–300.

SUPPORTING INFORMATION Additional Supporting Information may be found in the online version of this article at the publisher’s web- site: Figure S1. Scatterplot of the first shape principle component versus isometric size in selected species pairs of the genus Diplazon. Figure S2. Bayesian consensus trees of the concatenate dataset of five genes, and single-gene trees for the nuclear 28S rDNA and EF1a genes and the mitochondrial NADH1. Table S1. List of specimens and Genbank accession numbers. Table S2. Description of characters used in morphometrics of Diplazon. Table S3. Measurement data from morphometric analysis.

© 2016 The Linnean Society of London, Zoological Journal of the Linnean Society, 2016, 177, 541–557 0.4

0.2

species deletus ixi annulatus 0.0 tetragonus

PC1 (39.1 %) parvus varicoxa −0.2

−0.4 −0.2 0.0 0.2 0.4 Isometr ic siz e Diplazon_albotibialis Diplazon-pectoratorius_1C6 80 Diplazon-albotibialis_5-2 100 Diplazon_pectoratorius Diplazon-pectoratorius_3H2 Diplazon_laetatorius Diplazon-pallicoxa_5-3 100 Diplazon_orientalis Diplazon-tetragonus_3E2 100 Diplazon-tibiatorius_1B2 Diplazon_scutatorius Diplazon-annulatus_1C1 100 Diplazon_pallicoxa 100 Diplazon-tetragonus_3F2 Diplazon_bradleyi Diplazon-annulatus-cf_2D1 100 Diplazon-varicoxa-cf_1A9 60 Diplazon_deletus Diplazon-hyperboreus_9_03 Diplazon_annulatus 76 Diplazon-varicoxa_1A6 100 Diplazon_tetragonus Diplazon-tibiatorius-cf_3F1 Diplazon_flixi Diplazon-varicoxa-cf_6_16 99 Diplazon-tetragonus_1B1 0.02 100 Diplazon_hyperboreus 92 Diplazon-varicoxa-cf_6_17 Diplazon_zetteli Diplazon-varicoxa-cf_3D1 92 Diplazon-deletus_1A7 Diplazon_tibiatorius 88 0.003 Diplazon-deletus_1A8 99 Diplazon_angustus Diplazon-bradleyi_2A5 84 Diplazon_parvus Diplazon-tibiatorius_1B3 100 Diplazon_varicoxa Diplazon-zetteli_1A3 Diplazon-flixi_1A2 88 Diplazon-flixi_1A1 Diplazon-flixi_8_32 all genes, concatenated Diplazon-hyperboreus-cf_6_19 Diplazon-zettali_3C1 Diplazon-tibiatorius-cf_1B4 Diplazon-laetatorius_1C4 98 Diplazon-laetatorius_1C5 71 Diplazon-laetatorius_3A2 Diplazon-laetatorius_3B2 Diplazon-scutatorius_2E3 Diplazon-orientalis_7_1 73 28S rDNA, D2-D3 Diplazon-scutatorius_1A11

Diplazon_albotibialis_5_2 100 Diplazon_pectoratorius_1C6 Diplazon_albotibialis_5_2 Diplazon_tetragonus_4E11 100 Diplazon_pectoratorius_1C6 Diplazon_varicoxa_3D1 Diplazon_laetatorius_1C4 Diplazon_varicoxa_cf_6_17 100 100 Diplazon_laetatorius_1C5 Diplazon_tetragonus_2A10 100 100 Diplazon_orientalis_7_1 Diplazon_varicoxa_cf_1A9 100 Diplazon_scutatorius_1A11 Diplazon_tetragonus_1B1 89 91 Diplazon_scutatorius_2E3 100 Diplazon_deletus_1A8 97 Diplazon_pallicoxa_5_3 Diplazon_deletus_4H12 Diplazon_bradleyi_2A5 Diplazon_zetteli_1A3 86 Diplazon_varicoxa_cf_3D1 Diplazon_varicoxa_6_16 93 Diplazon_wymanni_6_17 79 Diplazon_tibiatorius-cf_1B4 Diplazon_varicoxa_1A6 Diplazon_tibiatorius_1B2 Diplazon_tibiatorius_1B3 Diplazon_bradleyi_2A5 69 Diplazon_varicoxa_cf_6_16 Diplazon_flixi_1A2 Diplazon_hyperboreus_6_19 0.008 Diplazon_flixi_2E5 Diplazon_zetteli_1A3 Diplazon_hyperboreus-cf_6_19 Diplazon_zetteli_1A Diplazon_pallicoxa_5_3 Diplazon_deletus_1A8 0.05 Diplazon_laetatorius_1C4 99 100 Diplazon_flixi_1A2 Diplazon_laetatorius_1C5 99 Diplazon_deletus_1A7 100 Diplazon_orientalis_7_1 Diplazon_tibiatorius_1B2 Diplazon_scutatorius_4F2 76 54 Diplazon_tetragonus_1B1 Diplazon_scutatorius_1A11 83 Diplazon_tibiatoriusCf_1B4 EF1a, F2 copy Diplazon_scutatorius_2E3 NADH1, mtDNA

Supplementary Figure S2. Bayesian consensus trees of the concatenate dataset of five genes, and single-gene trees for the nuclear 28S rDNA and EF1a genes and the mitochondrial NADH1. Support values next to node represent posterior probabilities in percent; values below 50% are ommited. Supplementary table 1. List of specimens and Genbank accession numbers

Species Sex Code Locality Coll1 CO1 ITS2 28S EF1α NADH1 Diplazon f 5_02 Switzerland, GL, Linthal, Obersand, Melchplatz, 2051m. NMBE GU597781 KR230498 GU597541 GU597608 GU597698 albotibialis N46°50.387/E8°55.821. 19.-28.VI.2008, leg. F.Marti Diplazon f 9_01 Sweden, Sm, Nybro kommun, Alsterbro/Alsteran. NRM GU597782 KR230499 GU597542 GU597610 GU597700 albotibialis N56°56'11"/E15°55'12". 25.V.-30.V.2005, leg. SMTP Diplazon f 1-C1 Finland, Muhos. 19.VIII.-3.IX.2005 NMBE KR230632 KR230500 annulatus Diplazon m 2-D4 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230633 KR230501 annulatus N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon f 2-D6 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230634 KR230502 annulatus N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon f 2-D7 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230635 KR230503 annulatus N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon f 4-G6 Switzerland, VS, Champéry, Col de Bretolet 1900m. NMBE KR230636 KR230504 annulatus N46°8.541/E6°47.727. 1.VIII.2007, leg. H.Baur Diplazon f 2-A5 USA, Alaska, Fairbanks, North Star Borough. 20.- NMBE FJ556426 KR230505 FJ556493 GU597602 GU597692 bradleyi 24.VI.2006, leg. D.Fieldling, N.Schiff Diplazon deletus f 1-A7 Switzerland, GR, Sur, Clavenia 1987m. NMBE FJ556427 KR230506 FJ556494 N46°32.005/E8°50.135. 23.-30.VI.2003, leg. H.Baur Diplazon deletus f 1-A8 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE FJ556428 KR230507 FJ556495 GU597603 GU597693 N46°31.482/E9°38.755. 23.-30.VI.2003, leg. H.Baur Diplazon deletus f 2-E10 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230637 KR230508 N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon deletus f 2-E9 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230638 KR230509 N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon deletus f 4-H10 Switzerland, VS, Champéry, Col de Bretolet 1922m. NMBE KR230639 KR230510 N46°8.568/E6°47.742. 2.-9.VIII.2007, leg. S.Klopfstein Diplazon deletus f 4-H12 Switzerland, VS, Champéry, Col de Bretolet 1922m. NMBE KR230640 KR230511 N46°8.568/E6°47.742. 2.-9.VIII.2007, leg. S.Klopfstein Diplazon deletus f 7_06 Sweden, Vr, Munkfors kommun, Ransäter, Rudstorp. NRM KR230642 KR230512 N59°46'22"/E13°28'25". 15.VII.-23.VII.2005, leg. SMTP Diplazon deletus f 7_14 Sweden, Vb, Vindelns kommun, Kulbäcken meadow. NRM KR230641 KR230513 N64°11.413'/E19°36.342'. 01.IX.-22.IX.2003, leg. SMTP Diplazon flixi f 1-A1 Switzerland, GR, Sur, NE Sur, 1770m. N46°31‘/E9°38‘. NMBE FJ556424 KR230514 FJ556491 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon flixi m 1-A2 Switzerland, GR, Sur, NE Sur, 1770m. N46°31‘/E9°38‘. NMBE FJ556425 KR230515 FJ556492 GU597601 GU597691 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon flixi f 2-E5 Switzerland, GR, Sur, NE Sur, 1770m. N46°31‘/E9°38‘. NMBE KR230643 KR230516 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon flixi m 2-E6 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230644 KR230517 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon flixi m 4-D12 Switzerland, VS, Champéry, Col de Bretolet 1900m. NMBE KR230645 KR230518 N46°8.541/E6°47.727. 1.VIII.2007, leg. H.Baur Diplazon flixi f 4-G11 Switzerland, VS, Champéry, Col de Bretolet 1922m. NMBE KR230646 KR230519 N46°8.568/E6°47.742. 23.-30.VIII.2007, leg. S.Klopfstein Diplazon flixi f 4-G12 Switzerland, VS, Champéry, Col de Bretolet 1900m. NMBE KR230647 KR230520 N46°8.541/E6°47.727. 1.VIII.2007, leg. H.Baur Diplazon flixi f 4-G7 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230648 KR230521 N46°31.482/E9°38.755. 23.-30.VI.2003, leg. H.Baur Diplazon flixi f 4-G9 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230649 KR230522 N46°31.482/E9°38.755. 23.-30.VI.2003, leg. H.Baur Diplazon flixi f 4-H1 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230650 KR230523 N46°31.482/E9°38.755. 16.-23.VI.2003, leg. H.Baur Diplazon flixi f 4-H3 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230651 KR230524 N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon flixi f 4-H5 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230652 KR230525 N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Diplazon flixi f 4-H8 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230653 KR230526 N46°31.482/E9°38.755. 9.-16.VI.2003, leg. H.Baur Diplazon flixi f 8_32 Switzerland, GL, Linthal, Obersand, Melchplatz, 2051m. NMBE KR230654 KR230527 N46°50.387/E8°55.821. 10.-19.VII.2008, leg. F.Marti Diplazon f 9_03 Sweden, BD, Kiruna kommun, Abisko NP, 900m. NRM KR230655 KR230528 hyperboreus N68°21.648'/E18°43.245'. 01.VII.-13.VII.2005, leg. SMTP Diplazon cf. f 6_19 Sweden, Up, Älvkarleby kommun, Bätfors. NRM GU597780 KR230529 GU597540 GU597604 GU597694 hyperboreus N60°27.639'/E17°19.069'. 27.VI.-01.VII.2004, leg. SMTP Diplazon f 1-C4 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. 1.- NMBE FJ556429 KR230530 FJ556496 GU597605 GU597695 laetatorius 25.IX.2006, leg. H.Baur & S.Klopfstein Diplazon f 1-C5 USA, MD. VII.2006 NMBE FJ556430 KR230531 FJ556497 GU597606 GU597696 laetatorius Diplazon f 3-A2 Spain, Canary Islands, La Gomera, Ermita Virgen de las NMBE FJ556431 KR230532 FJ556498 laetatorius Nieves, 1147m. N31.10418/E28.28359 Diplazon f 3-B2 Zambia, Choma. 13.-21.V.2006 NMBE FJ556432 KR230533 FJ556499 laetatorius Diplazon f 4-B10 Sweden, Ög, Ödeshögs kommun, Omberg, Stocklycke äng. NRM KR230656 KR230534 laetatorius N58°18.452'/E14°37.859'. 26.VII.-2.VIII.2005, leg. SMTP Diplazon f 4-G1 Sweden, Sm, Nybro kommun, Alsterbro/Alsteran. NRM KR230657 KR230535 laetatorius N56°56'11"/E15°55'12". 22.VI.-28.VI.2006, leg. SMTP Diplazon f 7_03 Switzerland, Bern, Bremgartenwald, clearing, 550m. NMBE KR230658 KR230537 laetatorius N46°57.577/E7°24.939. 23.-30.VII.2008, leg. S.Klopfstein Diplazon f 7_12 Sweden, Öl, Mörbylänga kommun, Frösslunda alvar. NRM KR230659 KR230536 laetatorius N56°32'/E16°34'. 5.VII.-2.VIII.2005, leg. SMTP Diplazon f 7_01 Thailand, Buri Prov., Kanchara, Thongpapoom, Chemical NMBE GU597782 KR230539 GU597542 GU597610 GU597700 orientalis Farm. 13.VI.2009, leg. D.Quicke Diplazon f 7_04 Thailand, Buri Prov., Kanchara, Thongpapoom, Chemical NMBE KR230660 KR230540 orientalis Farm. 13.VI.2009, leg. D.Quicke Diplazon f 7_07 Thailand, Buri Prov., Kanchara, Thongpapoom, Chemical NMBE KR230661 KR230541 orientalis Farm. 13.VI.2009, leg. D.Quicke Diplazon m 7_11 Thailand, Buri Prov., Kanchara, Thongpapoom, Chemical NMBE KR230662 KR230538 orientalis Farm. 13.VI.2009, leg. D.Quicke Diplazon f 5_03 Switzerland, Bern, Bremgartenwald, Nägelisbode, 540m. NMBE GU597783 KR230542 GU597543 GU597611 GU597701 pallicoxa N46°57‘/E7°25.057‘.VI.-8.VII.2008, leg. S.Klopfstein Diplazon f 6_01 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230663 KR230543 pallicoxa N56°55'17"/E16°6'4". 27.VI.-2.VII.2005, leg. SMTP Diplazon f 6_18 Sweden, Sö, Haninge kommun, Tyresta, Urskogsslingan. NRM KR230664 KR230544 pallicoxa N59°10'/E18°14'. 04.VIII.-26.VIII.2003, leg. SMTP Diplazon parvus f 6_23 Switzerland, Bern, Bremgartenwald, Nägelisbode, 540m. NMBE KR230665 KR230545 N46°57.955/E7°25.057. 11.-20.VI.2008, leg. S.Klopfstein Diplazon parvus f 6_24 Switzerland, Bern, Bremgartenwald, Nägelisbode, 540m. NMBE KR230666 KR230546 N46°57‘/E7°25‘. 27.VI.-8.VII.2008, leg. S.Klopfstein Diplazon parvus f 6_25 Switzerland, Bern, Bremgartenwald, clearing, 550m. NMBE KR230667 KR230547 N46°57.577/E7°24.939. 11.-20.VI.2008, leg. S.Klopfstein Diplazon parvus f 6_26 Switzerland, Bern, Bremgartenwald, clearing, 550m. NMBE KR230668 KR230548 N46°57.577/E7°24.939. 11.-20.VI.2008, leg. S.Klopfstein Diplazon parvus f 6_27 Sweden, Sö, Tyresö kommun, Ava, Spirudden. NRM KR230669 KR230549 N59°10.313'/E18°22.197'. 02.VII.-11.VII.2003, leg. SMTP Diplazon parvus f 6_29 Sweden, Up, Älvkarleby kommun, Bätfors. NRM KR230670 KR230550 N60°27.639'/E17°19.069'. 03.VII.-29.VII.2003, leg. SMTP Diplazon parvus f 8_44 Switzerland, Bern, Bremgartenwald, clearing, 550m. NMBE KR230671 N46°57.577/E7°24.939. 11.-20.VI.2008, leg. S.Klopfstein Diplazon cf. f 9_05 Sweden, Ha, Halmstad kommun, Gardshult, Buskastycket. NRM KR230672 parvus N56°41.693'/E13°09.030'. 1.IV.-25.IV.2004, leg. SMTP Diplazon cf. f 9_13 Sweden, Up, Alvkarleby kommun, Marma skjutfält. NRM KR230673 parvus N60°31'/E17°27'. 12.VIII.-26.VIII.2003, leg. SMTP Diplazon f 1-C6 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE FJ556433 KR230551 GU597544 GU597612 GU597702 pectoratorius N46°31.482/E9°38.755. 14.-21.VII.2003, leg. H.Baur Diplazon f 1-C7 Finland, Sipoonkorpi. 10.-16.VI.2006, leg. N.Laurenne NMBE KR230674 KR230552 pectoratorius Diplazon f 3-H2 Finland, Sipoonkorpi. 16.-21.VI.2006, leg. N.Laurenne NMBE FJ556434 KR230553 FJ556501 pectoratorius Diplazon f 4-F3 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230675 KR230554 pectoratorius N56°55'17"/E16°6'4". 24.VIII.-12.IX.2005, leg. SMTP Diplazon f 5_20 Switzerland, VS, Champéry, Col de Bretolet 1922m. NMBE KR230676 KR230555 pectoratorius N46°8.568/E6°47.742. 23.-30.VIII.2007, leg. S.Klopfstein Diplazon f 8_15 Sweden, Dr, Säterdalen, Näsakerspussen. N60°22'/E15°43'. NRM KR230677 KR230556 pectoratorius 23.VI-08.VII.2003, leg. SMTP Diplazon f 9_07 Sweden, Vr, Munkfors, Ransäter, Ransbergs herrgard. NRM KR230678 KR230557 pectoratorius N59°47'25"/E13°24'54". 27.VI.-10.VII.2005, leg. SMTP Diplazon f 1-A10 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230679 KR230558 scutatorius 15.-27.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 1-A11 Finland, Sipoonkorpi. 16.-21.VI.2006, leg. N.Laurenne NMBE FJ556435 KR230559 FJ556502 GU597613 GU597703 scutatorius Diplazon f 2-A3 Finland, Sipoonkorpi. 21.VI.-2.VII.2006, leg. N.Laurenne NMBE KR230680 KR230560 scutatorius Diplazon f 2-E12 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230681 KR230561 scutatorius 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 2-E3 Switzerland, Bern, Bremgartenwald, clearing, 550m. NMBE FJ556436 KR230562 FJ556503 GU597614 GU597704 scutatorius N46°57.577/E7°24.939. 15.IX.-4.X.2006, leg. S.Klopfstein Diplazon f 4-B12 Sweden, An, Ömsköldsviks kommun, Skuleskogen, NRM KR230682 KR230563 scutatorius Langra. N63°05'/E18°29'. 25.VII.-9.VIII.2004, leg. SMTP Diplazon f 4-D4 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230683 KR230564 scutatorius N56°55'17.96"/E16°6'4.42". 2.VII.-12.VII.2005, leg. SMTP Diplazon f 4-E3 Sweden, Up, Knivsta kommun, Rickebasta alsumpskog. NRM KR230684 KR230565 scutatorius N59°44.061'/E17°43.225. 24.VI.-16.VII.2005, leg. SMTP Diplazon f 4-F2 Sweden, Sm, Nybro kommun, Alsterbro/Alsteran. NRM KR230685 KR230566 scutatorius N56°56'11"/E15°55'12". 16.VII.-22.VII.2006, leg. SMTP Diplazon cf. f 7_09 Sweden, Öl, Mörbylanga, Gamla Skogsby (Kalkstad). NRM KR230686 KR230567 scutatorius N56°37'/E16°30'. 7.VIII.-18.IX.2003, leg. SMTP Diplazon cf. f 8_45 Sweden, Dr, Säterdalen, Näsakerspussen. N60°22'/E15°43'. NRM KR230687 KR230568 scutatorius 8.VII.-21.VII.2003, leg. SMTP Diplazon f 1-A12 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230688 KR230569 tetragonus N46°32.005/E8°50.135. 9.-16.VI.2003, leg. H.Baur Diplazon f 1-B1 Finland, Sipoonkorpi. 2.-12.VII.2006, leg. N.Laurenne NMBE FJ556437 KR230570 FJ556504 GU597706 tetragonus Diplazon f 2-A10 Finland, Sipoon, Hindsby. 10.-18.VI.2005, leg. N.Laurenne NMBE KR230689 KR230571 GU597616 tetragonus Diplazon m 2-A4 United Kingdom. 7.-13.VI.2005, leg. NMBE KR230690 KR230572 tetragonus Diplazon f 2-A7 Finland, Sipoonkorpi. 21.VI.-2.VII.2006, leg. N.Laurenne NMBE KR230691 KR230573 tetragonus Diplazon f 2-A8 USA, Alaska, Fairbanks, North Star Borough. 17.- NMBE KR230692 KR230574 tetragonus 20.VI.06, leg. D.Fieldling, N.Schiff Diplazon f 3-A3 Finland, Sipoonkorpi. 16.-21.VI.2006, leg. N.Laurenne NMBE KR230693 KR230575 tetragonus Diplazon f 3-B3 Finland, Sipoonkorpi. 16.-21.VI.2006, leg. N.Laurenne NMBE KR230694 KR230576 tetragonus Diplazon f 3-C2 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230695 KR230577 tetragonus N46°32.005/E8°50.135. 9.-16.VI.2003, leg. H.Baur Diplazon f 3-D2 Russia, Pskov Prov. NP Sebezhsky. 11.-20.VII.2006, leg. NMBE KR230696 KR230578 tetragonus A. Reshikov Diplazon m 3-E2 Russia, Pskov Prov. NP Sebezhsky. 11.-20.VII.2006, leg. NMBE FJ556438 KR230579 FJ556505 tetragonus A. Reshikov Diplazon f 3-E3 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230697 KR230580 tetragonus 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 3-F2 Finland, Sipoon, Hindsby. 13.-29.VI.2006, leg. N.Laurenne NMBE FJ556439 KR230581 FJ556506 tetragonus Diplazon f 3-G2 Germany, Bayern, Nationalpark Bayerischer Wald. VI.- NMBE KR230698 KR230582 tetragonus VII.2001, leg. G.Görsen Diplazon f 4-B7 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230699 KR230583 tetragonus N56°55'53''/E16°5'7". 13.VIII.-24.VIII.2005, leg. SMTP Diplazon f 4-G3 Sweden, Sm, Nybro kommun, Alsterbro/Alsteran. NRM KR230700 KR230584 tetragonus N56°56'11"/E15°55'12". 22.VI.-28.VI.2006, leg. SMTP Diplazon f 8_09 Sweden, Hs, Hudiksvalls kommun, Stensjön-Lomtjärn, NRM KR230701 KR230585 tetragonus Stensjön. N62°08'/E16°17'. 11.-25.VIII.2004, leg. SMTP Diplazon f 8_43 Switzerland, Bern, Bremgartenwald, Nägelisbode, 540m. NMBE KR230702 KR230586 tetragonus N46°57.955/E7°25.057. 20.-27.VI.2008, leg. S.Klopfstein Diplazon f 9_11 Sweden, Up, Alvkarleby kommun, Marma skjutfält. NRM KR230703 tetragonus N60°31'/E17°27'. 12.VIII.-26.VIII.2003, leg. SMTP Diplazon m 2-A9 Finland, Kangasilampi. 7.VII.1998, leg. N.Laurenne NMBE KR230704 KR230587 tetragonus Diplazon f 4-E11 Sweden, Dr, Säterdalen, Näsakerspussen. N60°22'/E15°43'. NRM KR230705 KR230588 tetragonus 23.VI-08.VII.2003, leg. SMTP Diplazon f 4-F6 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230706 KR230589 tetragonus N56°55'17"/E16°6'4". 24.VIII.-12.IX.2005, leg. SMTP Diplazon f 1-B2 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE FJ556440 KR230590 FJ556507 GU597617 GU597707 tibiatorius 19.-27.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 1-B3 Russia, Primorsk flood plain. NMBE FJ556441 KR230591 FJ556508 tibiatorius Diplazon cf. f 1-B4 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. 1.- NMBE FJ556442 KR230592 FJ556509 GU597618 GU597708 tibiatorius 25.IX.2006, leg. H.Baur & S.Klopfstein Diplazon cf. f 3-E1 Switzerland, GR, Sur, Lai Neir, 1920m. N46°31‘/E9°38‘. NMBE KR230707 KR230593 tibiatorius 7.-15.VII.2006, leg. H.Baur & S.Klopfstein Diplazon cf. f 3-F1 Switzerland, GR, Sur, Lai Neir, 1920m. N46°31’/E9°38’. NMBE FJ556443 KR230594 FJ556510 tibiatorius 7.-15.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 1-A5 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230708 KR230595 varicoxa N46°32.005/E8°50.135. 14.-21.VII.2003, leg. H.Baur Diplazon f 1-A6 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230709 KR230596 varicoxa 14.VII.-1.VIII.2006, leg. H.Baur & S.Klopfstein Diplazon m 1-A9 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230710 KR230597 varicoxa 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 2-E11 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230713 KR230600 varicoxa N46°31.482/E9°38.755. 18.-25.VIII.2003, leg. H.Baur Diplazon m 2-E7 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230711 KR230599 varicoxa 14.VII.-1.VIII.2006, leg. H.Baur & S.Klopfstein Diplazon m 2-E8 Switzerland, GR, Sur, NE Sur, 1770m. N46°31’/E9°38’. NMBE KR230712 KR230598 varicoxa 15.-20.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 2-F1 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230714 KR230601 varicoxa N46°31.482/E9°38.755. 4.-11.VIII.2003, leg. H.Baur Diplazon m 3-D1 Switzerland, GR, Sur, Lai Neir, 1920m. N46°31’/E9°38’. NMBE GU597786 KR230602 GU597545 GU597619 GU597709 varicoxa 7.-15.VII.2006, leg. H.Baur & S.Klopfstein Diplazon f 4-D3 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230715 KR230603 varicoxa N56°55’53’’/E16°5’7“. 02.VII.-12.VII.2005, leg. SMTP Diplazon f 5_21 Switzerland, Bern, Bremgartenwald, Nägelisbode, 540m. NMBE KR230716 KR230604 varicoxa N46°57.955/E7°25.057. 23.-30.VII.2008, leg. S.Klopfstein Diplazon m 6_13 Sweden, Sm, Almhults kommun, Stenbrohult, Djäknabygds NRM KR230717 KR230605 varicoxa bokbacke. N56°36'/E14°11'. 26.-15.VII.2003, leg. SMTP Diplazon m 6_14 Sweden, Sm, Almhults kommun, Stenbrohult, Djäknabygds NRM KR230718 KR230606 varicoxa bokbacke. N56°36'/E14°11'. 15.-1.VIII.2003, leg. SMTP Diplazon f 6_15 Sweden, Sm, Almhults kommun, Stenbrohult, Djäknabygds NRM KR230719 KR230607 varicoxa bokbacke. N56°36'/E14°11'. 15.-1.VIII.2003, leg. SMTP Diplazon f 6_16 Sweden, Sk, Klippans kommun, Skäralid, N Liema. NRM GU597785 KR230608 GU597546 GU597607 GU597697 varicoxa N56°01'/E13°13'. 14.VII.-06.VIII.2004, leg. SMTP Diplazon f 6_20 Switzerland, GL, Linthal, Obersand, Melchplatz, 2051m. NMBE KR230720 KR230609 varicoxa N46°50.387/E8°55.821. 28.VI.-10.VII.2008, leg. F.Marti Diplazon f 6_30 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230721 KR230610 varicoxa N56°55’53’’/E16°5’7“. 13.VIII.-24.VIII.2005, leg. SMTP Diplazon f 6_31 Sweden, Up, Älvkarleby kommun, Bätfors. NRM KR230722 KR230611 varicoxa N60°27.639'/E17°19.069'. 27.VI.-01.VII.2004, leg. SMTP Diplazon f 6_32 Sweden, Sm, Almhults kommun, Stenbrohult, Djäknabygds NRM KR230723 KR230612 varicoxa bokbacke. N56°36'/E14°11'. 1.-18.VIII.2003, leg. SMTP Diplazon f 6_33 Sweden, Sk, Klippans kommun, Skäralid, N Liema. NRM KR230724 KR230613 varicoxa N56°01'/E13°13'. 14.VII.-06.VIII.2004, leg. SMTP Diplazon f 8_25 Sweden, Sm, Almhults kommun, Stenbrohult, Djäknabygds NRM KR230725 KR230614 varicoxa bokbacke. N56°36'/E14°11'. 19.V.-27.VI.2004, leg. SMTP Diplazon f 8_46 Switzerland, Bern, Bremgartenwald, clearing, 550m. NMBE KR230726 KR230615 varicoxa N46°57.577/E7°24.939. 11.-20.VI.2008, leg. S.Klopfstein Diplazon f 8_48 Sweden, Sm, Nybro kommun, Bäckebo, Grytsjöns NR. NRM KR230727 KR230616 varicoxa N56°55’53’’/E16°5’7“. 02.VII.-12.VII.2005, leg. SMTP Diplazon f 9_09 Sweden, Sm, Nybro kommun, Alsterbro/Alsteran. NRM KR230729 KR230617 varicoxa N56°56'11"/E15°55'12". 25.V.-30.V.2005, leg. SMTP Diplazon f 9_21 Sweden, Vr, Munkfors kommun, Ransäter, Rudstorp. NRM KR230728 KR230618 varicoxa N59°46'22"/E13°28'25". 7.VII.-15.VII.2005, leg. SMTP Diplazon cf. f 6_17 Sweden, BD, Gällivare kommun, Ätnarova försökspark, NRM GU597784 KR230619 GU597547 GU597609 GU597699 varicoxa Pelttovaara. N67°03'/E20°23'. 29.-13.VIII.2004, leg. SMTP Diplazon cf. f 9_02 Sweden, Up, Alvkarleby kommun, Marma skjutfält. NRM KR230730 varicoxa N60°31'/E17°27'. 29.VII.-12.VIII.2003, leg. SMTP Diplazon cf. f 9_23 Sweden, Up, Alvkarleby kommun, Marma skjutfält. NRM KR230731 varicoxa N60°31.456'. 12.VIII.-26.VIII.2003, leg. SMTP Diplazon zetteli f 1-A3 Switzerland, GR, Sur, Clavenia 1987m. NMBE GU597787 KR230620 GU597548 GU597615 GU597705 N46°32.005/E8°50.135. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli f 1-A4 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230732 KR230621 N46°32.005/E8°50.135. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli f 2-D2 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230733 KR230622 N46°32.005/E8°50.135. 9.-16.VI.2003, leg. H.Baur Diplazon zetteli m 2-D3 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230734 KR230623 N46°32.005/E8°50.135. 9.-16.VI.2003, leg. H.Baur Diplazon zetteli f 2-D5 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230735 KR230624 N46°31.482/E9°38.755. 9.-16.VI.2003, leg. H.Baur Diplazon zetteli f 2-E1 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230736 KR230625 N46°32.005/E8°50.135. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli m 2-E2 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230737 KR230626 N46°32.005/E8°50.135. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli m 3-A1 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230738 KR230627 N46°32.005/E8°50.135. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli f 3-B1 Switzerland, GR, Sur, Clavenia 1987m. NMBE KR230739 KR230628 N46°32.005/E8°50.135. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli m 3-B4 Switzerland, GR, Sur, Piz d'Umblei, 1928m, NMBE KR230740 N46°31.655/E9°38.674. 3.VI.2003, leg. H.Baur Diplazon zetteli f 3-C1 Switzerland, GR, Sur, SE Vauastg Dafora, 1920m. NMBE KR230741 KR230629 N46°31.482/E9°38.755. 2.-9.VI.2003, leg. H.Baur Diplazon zetteli f 6_21 Switzerland, GL, Linthal, Obersand, Melchplatz, 2051m. NMBE KR230742 KR230630 N46°50.387/E8°55.821. 19.-28.VI.2008, leg. F.Marti Diplazon zetteli f 6_22 Switzerland, GL, Linthal, Obersand, Melchplatz, 2051m. NMBE KR230743 KR230631 N46°50.387/E8°55.821. 28.VI.-10.VII.2008, leg. F.Marti Campocraspedon f 1-C2 Switzerland, GR, Sur, Clavenia 1987m. NMBE GU597779 GU597539 GU597600 GU597690 caudatus N46°32.005/E8°50.135. 16.-23.VI.2003, leg. H.Baur Syrphophilus f 5_14 Sweden, Vr, Munkfors kommun, Ransäter, Ransb herrgard. NRM GU597819 GU597580 GU597663 GU597753 asperatus N59°47'25"/E13°24'54". 18.-27.VI.2005, leg. SMTP Tymmophorus f 1-H2 United Kingdom, leg. G.Broad NMBE FJ556476 FJ556543 GU597670 GU597758 obscuripes Tymmophorus f 2-B11 Finland, Muhos. 5.-19.VIII.2005, leg. N.Laurenne NMBE FJ556478 FJ556545 GU597671 GU597761 suspiciosus Xestopelta f 3-C3 United Kingdom, Hilbre Island, on blackthorn. NMBE GU597829 GU597590 (GU597680) GU597770 gracilima 21.VII.2001, leg. G.Broad

1 Collections: NMBE = Natural History Museum Bern, Switzerland. NRM = Naturhistoriska Riksmuseet, Stockholm, Sweden.

Description of characters used in morphometrics of Diplazon

Character Description Magnif Foto

Mesoscutum From front margin to the 32 length prescutellar groove

Mesoscutum Measured in front of the tegulae, 32 width including the margin

Tergite 1 length Front end of tergite, where it 32 inserts into propodeum, to the end where the lateral carinae meet

Tergite 2 length Meeting point of upper and lower 32 lateral carinae at the base of the tergite, to its apical corner, including the folded membrane

Tergite 3 length As for tergite 2 32

Tergite 1 width Width is measured at the point 32 where the tergite is broadest in the apical quarter

Tergite 2 width Width is measured at the end of 32 the tergite, including the folded membranes

Hind tibia Upper and outer end where the 32 length tibia meets the femur, to the round expansion on the outer side of the tibia (longest length)

Flagellomere 1 Length measured on the outer 80 length side of flagellomere

Flagellomere 2 Length measured on the outer 80 length side of flagellomere

Face height Height measured from the 80 tentorial pits to the lower margin of the antennal sockets

Face width Width of face measured just 80 below the antennal sockets

Interocellar Shortest distance between the 80 distance hind ocelli

Ocellocular Shortest distance between the 80 distance hind ocellus and the compound eye

Hindtarsus 2 Length measured on outer side 80 length (longest length; tarsus 4 and 5 shown)

Hindtarsus 4 Length measured on outer side 80 length (longest length; tarsus 4 and 5 shown)

Individual# species mes.l mes.w t1.l t2.l t3.l 1 annulatus 1049.94 1144.24 825.04 703.77 641.22 20 annulatus 925.76 990.06 727.27 580.21 566.55 562 annulatus 976.91 1044.21 701.03 663.68 592.36 10 deletus 820.67 857.45 678.1 599.77 534.22 11 deletus 832.09 881.32 614.69 567.66 494.35 40 deletus 845.64 919.33 665.45 589.23 497.53 46 deletus 930.1 986.36 728.18 649.12 575.73 615 deletus 713.46 749.3 540.67 442.19 413.12 620 deletus 1093.88 1107.38 866.54 727.56 683.91 623 deletus 782.73 826.2 618.62 521.41 457.66 625 deletus 956.77 1005.03 786.63 666.11 563.27 628 deletus 773.41 833.93 617.57 532.39 430.48 631 deletus 905.61 952.37 706.71 619.63 570.14 643 deletus 961.82 998.51 702.37 680.38 563.78 659 deletus 591.39 668.94 485.13 434.17 409.81 552 flixi 969.77 979.65 801.93 567.79 521.47 553 flixi 1034.61 1085.11 777.77 642.85 614.25 554 flixi 866.31 936.14 622.72 559.91 460.46 555 flixi 1102.96 1115.77 844.82 700.66 636.22 556 flixi 869.61 944.44 636.92 563.25 511.09 557 flixi 954.99 1009.99 744.35 616.95 554.52 558 flixi 1111.62 1156.44 856.18 739.83 656.93 559 flixi 1152 1182.26 888.8 685.07 630.27 690 flixi 951.97 1004.98 694.57 567.96 534.91 657 parvus 1002.34 1157.2 772.13 699.19 565.99 658 parvus 1022.81 1124.57 738.18 693.55 608.71 660 parvus 1096.6 1208.78 815.1 702.39 633.68 676 parvus 868.6 1001.91 683.55 558.38 492.99 662 parvus 1126.24 1237.84 905.52 782.55 671.14 663 parvus 944.27 1007.71 771.82 634.56 577.46 16 tetragonus 683.37 778.59 526.05 448.71 399.27 76 tetragonus 793.74 847.96 551.97 529.83 481.97 648 tetragonus 873.07 956.53 641.76 597.41 518.33 655 tetragonus 853.59 893.46 599.47 554.79 494.65 8 tetragonus 1076.57 1147.83 814.55 693.63 587.45 13 tetragonus 858.61 914.87 684.09 593.15 545.08 19 tetragonus 1012.38 1066.69 733.02 681.51 614.17 26 tetragonus 1009.39 1085.16 739.46 668.48 608.74 38 tetragonus 1088.03 1157.97 844.26 732.95 632.19 41 tetragonus 1085.33 1144.37 796.86 718.16 666.75 42 tetragonus 918.09 977.2 684.76 592.83 471.34 55 tetragonus 1025.84 1074.94 814.74 710.67 554.06 64 tetragonus 1086.42 1102.37 814.3 691.58 574.37 66 tetragonus 1114.54 1131.37 840.78 716.52 633.62 73 tetragonus 962.53 1038.95 682.24 612.85 566.66 634 tetragonus 974.18 1007.44 666.12 605.47 497.75 639 tetragonus 1046.13 1072.58 772.81 664.03 602.82 646 tetragonus 1014.29 1120.86 771.93 694.32 555.47 650 tetragonus 958 1048.6 693.6 602.01 535.94 696 tetragonus 896.55 998.34 690.52 650.87 578.73 44 tibiatorius 1401.04 1492.92 1069.02 897.23 776.66 656 tibiatorius 1570.41 1608.33 1159.65 904.45 943.12 9 tibiatorius 1246.26 1252.45 864.22 726.42 685.87 12 tibiatorius 1265.65 1365.96 988.62 869.06 791.75 83 tibiatorius 1365.75 1456.58 1097.2 931.39 863.49 178 tibiatorius 1593.14 1624.79 1224.58 1041.51 930.38 621 tibiatorius 1396.05 1449.12 1023.14 856.19 730.36 626 tibiatorius 1462.14 1512.29 1045.67 882.77 816.25 629 tibiatorius 1512.03 1589.87 1165.44 987.04 845.85 636 tibiatorius 1105.8 1332.91 1021.65 837.25 721.2 638 tibiatorius 1490.2 1488.18 1129.29 934.78 906.81 640 tibiatorius 1357.89 1358.58 987.86 848.05 771.4 641 tibiatorius 1485.17 1656.99 1176.87 994.91 832.53 647 tibiatorius 1444.76 1436.22 1160.05 991.04 914.49 685 tibiatorius 1464.3 1456.91 1184.58 947.84 917.55 39 tibiatorius 1400.39 1439.37 1166.23 826.4 817.81 81 tibiatorius 1316.2 1309.9 979 875.47 748.7 693 tibiatorius 1433.39 1494.77 1141.21 962.92 817.85 694 tibiatorius 1342.66 1431 1097.29 877.26 829.59 28 tibiatoriusCf 1414.72 1459.34 1138.77 915.68 778.74 43 tibiatoriusCf 1466.5 1433.9 1075.95 806.28 883.11 632 tibiatoriusCf 1276.5 1471.53 1042.68 866.37 836.96 56 tibiatoriusCf 1517.65 1496 1206.85 972.34 901.39 58 tibiatoriusCf 1350.22 1407.94 1003.24 767.92 788.99 77 tibiatoriusCf 1303.72 1366.95 980.04 848.89 808.75 27 varicoxa 1461.44 1452.15 1017.93 844.57 706.23 644 varicoxa 1015.61 1157.36 826.99 633.82 535.28 651 varicoxa 935.15 1100.02 739.03 643.16 519.96 79 varicoxa 1358.66 1490.76 1181.54 970.79 793.93 589 varicoxa 1283.22 1490.11 944.63 830.36 674.55 622 varicoxa 1332.85 1490.81 1075.34 852.2 715.09 3 varicoxa 1341.96 1409.31 984.23 828.05 709.73 59 varicoxa 1437.49 1441.13 1034.56 841.65 660.43 75 varicoxa 1242.02 1221.67 885.03 765.9 669.36 617 varicoxa 1489.8 1541.78 1169.93 1004.44 805.32 635 varicoxa 1325.04 1451.1 998.98 948.52 791.73 645 varicoxa 938.98 1066.51 727.43 574.29 537.34 654 varicoxa 1214.82 1309.11 955.34 802.86 631.06 661 varicoxa 1140.09 1247.56 859.89 717.58 539.26 668 varicoxa 1129.74 1219.71 817.33 686.78 567.81 670 varicoxa 1392.58 1414.45 1024.83 877.81 721.86 t1.w t2.w htib.wb htib.l fl1.l fl2.l HOc.d 823.93 1125.77 962.69 1251.54 409.3 290.02 223.16 638.74 956.5 994.88 1649.38 386.14 269.67 197.17 648.62 1040.31 962.69 1658.19 405.69 273.95 196.71 496.64 863.35 876.52 1526.38 365.09 260.09 178.11 493.58 779.24 832.21 1500.79 366.46 272.81 175.5 487.72 798.57 901.79 1538.31 398.56 282.75 185.8 547.41 871.94 1058.81 1649.04 413 300.63 192.52 464.47 740.89 696.14 1269.33 338.1 247.74 166.56 673.08 1090.31 1104.03 1856.38 495.91 347.58 239.53 475.7 765.63 866.03 1375.53 371.4 272.54 164.01 619.97 954.13 966.87 1669.09 431.13 312.51 222.03 492.28 786.16 755.02 1384.94 374.35 258.94 163.4 503.76 895.65 899.06 1498.62 424.87 285.3 188.07 643.58 990.55 1051.86 1704.29 442.08 290.46 219.89 384.27 701.27 677.92 1146.32 309.02 233.48 162.6 692.03 1020.71 789.64 1527.15 372.21 264.75 224.84 728.48 1054 977.84 1660.25 419.64 291.71 227 607.95 896.14 707.83 1366.28 329.98 243.77 227.04 729.62 1086.98 985.52 1774.06 437.98 295.27 259.76 607.12 873.81 827.07 1444.42 363.99 252.19 196.11 668.19 1005.77 878.1 1617.45 398.86 278.22 209.69 779.17 1137.54 1057.41 1742.79 462.52 314.11 247.9 818.61 1207.64 960.06 1751.84 439.47 306.28 248.97 649.12 917.55 809.34 1540.49 388.12 274.18 211.97 721.06 1138.73 816.82 1605.59 409.98 259.75 251.38 717.76 1137.56 826.13 1658.13 401.94 279.21 253.22 817.61 1149.61 905.4 1674.98 429.01 288.03 250.02 608.78 954.43 747.77 1407.07 352.43 245.84 190.49 942.13 1373.55 944.09 1896.5 428.85 285.1 276.39 724.47 1113.58 826.7 1600.14 357.88 259.6 250.64 481.69 690.83 717.45 1236.61 319.13 233.41 164.13 502.2 777.73 976.7 1440.33 352.12 256.56 181.99 613.65 889.48 844.64 1446.58 390.13 255.65 186.6 511.85 783.4 714.25 1326.65 337.13 253.07 194.2 741.61 1110.01 1119.56 1783.91 438.61 313.8 188.19 599.61 940.32 932.06 1527.37 378.68 290.16 201.02 679.39 1046.15 1090.1 1760.93 424.21 298.41 227.45 707.98 1091.97 1084.47 1685.59 372.22 298.93 203.12 745.59 1181.47 1040.12 1795.57 454.55 313.14 201.67 731.18 1137.14 1164.65 1822.17 456.04 322.51 235.62 652.67 996.97 956.81 1525.97 396.01 295.03 203.64 687.66 1129.85 1016.63 1706.55 439.46 314.75 224.58 754.73 1117.71 1090.38 1769.63 420.98 322.91 217.23 761.61 1110.47 1165.47 1805.25 460.81 321.68 229.04 673.09 974.89 1207.04 1631.13 402.02 272.9 198.35 659.41 1028.58 959.59 1542.66 424.89 280.82 217.84 696.58 1170.25 1036.43 1812.12 482.63 324.56 230.36 700.24 1162.26 1155.63 1791.28 474.26 288.49 201.96 650.39 974.76 912.46 1523.73 407.44 302.9 199.38 686 1096.77 817.97 1445.85 367.71 259.3 228.13 1150.99 1566.35 1368.4 2318.32 614.01 382.92 291.16 1133.57 1556.89 1370.19 2488.16 590.68 375.35 289.95 951.49 1250.83 960.35 1874.82 484.22 330.36 305.33 1024.62 1467.55 912.08 1930.54 462.83 328.16 238.36 1087.16 1525.27 1270.31 2367.19 554.64 363.73 323.5 1212.55 1806.1 1306.16 2624.07 563.94 424.58 349.22 1196.56 1720.66 1188.72 2366.62 554.44 398.21 286.68 1175.9 1644.36 1248.66 2366.01 587.87 376 334.74 1209.61 1798.5 1361.44 2554.41 619.42 377.44 335.24 1081.32 1544.83 1153.16 1912.96 501.28 325.93 291.45 1144.9 1634.05 1447.32 2599.12 555.42 370.21 330.76 934.49 1381.48 1243.63 2158.26 517.25 363.6 305.64 1234.91 1705.47 1324.31 2340.92 479.4 380.36 354.57 1149.13 1622.81 1332.13 2375.03 525.57 353.68 323.03 1117.66 1600 1249.71 2332.49 568.73 346.58 274.94 1124.84 1545.1 1297.69 2352.49 553.29 357.57 315.37 985.48 1382.4 1211.9 2123.15 541.25 353.17 313.89 1242.13 1724.63 1366.64 2427.19 567.22 394.41 342.48 1116.35 1593.69 1250.26 2237.95 537.13 337.42 309.76 1080.21 1587.93 1118.06 2238.74 556.79 346.91 312.12 1209.53 1719.33 1308.54 2154.12 612.36 384.63 286.07 1030.7 1574.72 1178.96 2110.07 553.76 358.11 297.98 1104.85 1600.81 1504.96 2547.84 544.72 360.46 304.6 1025.05 1386.12 1168.62 2184.89 549.7 342.46 289.49 1020.95 1474.54 985.84 2021.36 527.19 322.96 281.14 1028.61 1591.74 1097.55 2174.3 497.59 324.52 271.96 762.46 1050.93 1063.79 1654.67 399.73 279.69 242.08 717.42 1031.59 814.47 1514.11 364.58 266.1 236.68 1098.4 1646.79 1178.84 2268.9 533.89 326.17 285.23 1003.62 1413 886.36 1902.78 457.96 309.92 271.96 1065.99 1545.01 1141.82 2005.72 549.58 342.4 294.64 924.68 1321.61 792.74 1965.52 540.29 398.81 268.74 997.32 1427.59 922.27 2097.16 526.22 390.1 249.8 785.21 1119.66 1014.78 1953.04 534.77 349.19 232.64 1074.83 1575.32 1105.12 2313.89 544.02 357.19 286.34 828.35 1320.34 965.86 2064.23 561.73 366.02 227.51 671.73 1075.41 841.15 1533.95 381.22 265.42 206.42 1016.7 1380 918.47 1808.58 430.14 276.63 244.86 921.79 1241.13 959.44 1692.1 403.87 261.01 246.11 708.17 1056.6 918 1783.82 462.6 339.82 212.44 926.42 1441.13 1045.99 2167.11 538.77 400.89 281.36 HOc.eye htar2.l htar4.l 152.83 461.46 186.9 135.44 435.84 199.6 156.63 386.26 165.57 129.9 394.58 158.87 116.39 390.58 151.9 141.97 425.82 176.46 143.02 443.5 186.62 142.51 332.42 129.85 165.54 492.69 183.09 129.05 384.68 153.65 137.04 434.03 171.68 123.15 366.94 148.48 148.69 448.25 167.3 137.81 450.37 179 116.74 300.63 131.97 175.02 383.59 163.14 150.83 450.53 177.38 163.14 368.07 152.98 184.8 475.49 201.2 165.12 380.39 167.59 179.15 423.24 176.54 175.31 473.46 196.72 185.88 495.87 211.08 171.52 415.22 178 162 417.64 163 154.24 394.63 161.67 148.27 417.11 181.5 150.93 334.24 149.34 154.93 394.68 151.4 151.97 394.68 151.4 139.53 296.95 123.22 128.29 346.36 144.49 141.51 376 153.83 127.33 371.13 157.07 142.04 415.7 171.18 128.38 368.76 146.72 137.2 414.21 158.31 142.69 427.54 178.39 148.36 408.18 150.8 155.83 484.51 201.24 133.74 380.83 146.63 129.69 392.71 220.83 143.02 455.37 164.89 139.47 442.19 169.65 150.62 399.78 153.28 210.46 374.67 145.78 147.62 441.33 173.59 137.37 443.2 167.45 127.26 399.96 152.71 127.32 345.74 139.92 253.01 630.47 235.7 209.61 632.1 252.75 169.01 520.93 196.09 186.26 494.19 203.48 191.34 634.44 232.37 222.46 615.23 249.21 204.65 624.75 241.34 222.39 658.2 262.87 255.88 685.34 257.99 195.08 509.13 210.51 194.87 626.44 244.13 210.94 570.07 222.74 204.72 608.03 238.15 227.54 572.05 245.77 212.07 581.01 225.72 221.59 559.39 213.7 208.04 555.52 224.9 224.12 629.05 249.89 199.81 562.09 231.17 221.27 614.27 237.05 220.39 583.92 238.92 221.05 514.73 215.79 238.8 642.86 264.54 231.43 566.89 226.56 218.17 512.62 219.8 202.93 493.38 225.02 164.08 389.85 164.59 152.37 360.38 166.97 214.35 523.78 217.19 176.78 459.01 186.52 179.58 508.6 212.8 181.78 518.6 219.67 173.07 501.04 210.57 166.39 451.07 186.32 287.99 578.74 236.12 186.64 495.47 218.83 144.62 381.54 174.14 179.37 434.93 187.4 172.62 412.19 183.55 165.43 438.22 203.56 194.81 534.87 253.11 Figure S1. Scatterplot of the first shape principle component versus isometric size in selected species pairs of the genus Diplazon.

Table S3. Measurement data from morphometric analysis.