High molecular and phenotypic diversity in theMerodon aviduscomplex (Diptera, Syrphidae): cryptic speciation in a diverse insect taxon

Re-description, systematics and complete mitochondrial genome of Philheliuscoreanus (Shiraki, 1930) (Diptera, Syrphidae) in the Republic of Korea

Background

The hoverfly Philheliuscoreanus (Shiraki, 1930) was first described, based on only Korean male specimens and subsequent descriptions of the female from Russia did not include discussions of phenotypic variation. Furthermore, full-length mitochondrial genome sequences for the genus are lacking.

New information

To address these gaps, we here provide a diagnosis, re-description and mitochondrial genome of Philheliuscoreanus (Shiraki, 1930). We evaluated genitalic characters of both males and females with colour photographs and they showed intraspecific variation. There was significant variation in the yellow spots on the pleuron, particularly in females. After obtaining the complete mitochondrial genome of P.coreanus, we performed a phylogenetic analysis using Maximum Likelihood, based on 13 protein-coding genes, with a focus on relationships within the tribe Syrphini. Our results supported the monophyly of Syrphini, showing a sister-group relationship between Philhelius and Doros Meigen, 1822. Furthermore, the Philhelius + Doros clade was closely related to the Chrysotoxum + Dideopsis clade, with relatively high support. The newly-obtained mitochondrial genome of P.coreanus and high-resolution phylogenetic analysis provide essential resources for further analyses of the genus and relationships within Syrphini.

Molecular differentiation analysis of ten putative species of Fannia (Diptera: Fanniidae) collected in carrion-baited traps from Colombia

The genus Fannia is the most representative of the Fannidae family of true flies with worldwide distribution. Some species are attracted to decomposing materials and live vertebrate animals, which makes them important in forensics, medical and veterinary fields. However, identifying Fannia species can be difficult due to the high similarity in the external morphology of females and limited descriptions and morphological keys. Herein, molecular markers could provide a complementary tool for species identification. However, molecular identification has still limited application since databases contain few data for neotropical species of Fannia. This study assessed the potential of two molecular markers, the COI-3' region and internal transcribed spacer 2 (ITS2), to differentiate 10 putative species of the genus Fannia from Colombia using distance-based and tree-based approaches. The partial ITS2 and/or COI-3' regions allowed molecular diagnosis of six species, while pairs of species Fannia colazorrensis + F. dodgei and F. laclara + F. aburrae are conflicting. Although these results might suggest that conflicting pair species are conspecific, consistent morphological differences between males do not support this hypothesis. The lack of differentiation at the nuclear and mitochondrial molecular markers for the conflicting species may be due to incomplete evolutionary lineage separation, hybridization, or introgression events. In addition, sexual selection on male morphological traits before species-specific differences in molecular markers emerge may partially explain the results. Our study provides a valuable dataset to identify and confirm some Fannia species molecularly. Further, they could be used to associate females and immature stages with their conspecifics as a baseline to deep into their biology, ecology, distribution and potential applications in forensic and medico-veterinary entomology.

Unveiling the Wing Shape Variation in Northern Altiplano Ecosystems: The Example of the Butterfly Phulia nymphula Using Geometric Morphometrics

The Andean Altiplano, characterized by its extreme climatic conditions and high levels of biodiversity, provides a unique environment for studying ecological and evolutionary adaptations in insect morphology. Butterflies, due their large wing surface compared to body surface, and wide distribution among a geographical area given the flight capabilities provided by their wings, constitute a good biological model to study morphological adaptations following extreme weathers. This study focuses on Phulia nymphula, a butterfly species widely distributed in the Andes, to evaluate wing shape variation across six localities in the Northern Chilean Altiplano. The geometric morphometrics analysis of 77 specimens from six locations from the Chilean Altiplano (Caquena, Sorapata Lake, Chungará, Casiri Macho Lake, Surire Salt Flat, and Visviri) revealed significant differences in wing shape among populations. According to the presented results, variations are likely influenced by local environmental conditions and selective pressures, suggesting specific adaptations to the microhabitats of the Altiplano. The first three principal components represented 60.92% of the total wing shape variation. The detected morphological differences indicate adaptive divergence among populations, reflecting evolutionary responses to the extreme and fragmented conditions of the Altiplano. This study gives insights into the understanding of how high-altitude species can diversify and adapt through morphological variation, providing evidence of ecological and evolutionary processes shaping biodiversity in extreme environments.

The Extraordinary Diversity of Merodon avidus Complex (Diptera: Syrphidae)-Adding New Areas, New Species and a New Molecular Marker

In this paper, the Merodon avidus (Diptera, Syrphidae) species complex was revised, whereupon we discovered and described four new species for science: Merodon atroavidus Vujić, Radenković et Likov sp. nov., M. magnus Vujić, Kočiš Tubić et Ačanski sp. nov., M. nigroscutum Vujić, Radenković et Likov sp. nov. and M. pseudomoenium Vujić, Kočiš Tubić et Ačanski sp. nov. An integrative taxonomy approach was used to delimit species boundaries. Two molecular markers (the mitochondrial COI gene and nuclear 28S rRNA gene-newly analysed marker for the complex) and geometric morphometry of the wing shape, together with morphological data and distribution, successfully separated all species from the complex. The morphological variability of the analysed species is described and discussed and an illustrated diagnostic key for typical morpho-forms of species from the M. avidus complex is presented. A distribution map of all investigated species from the complex is provided. The level of endemicity of the M. avidus complex was discussed.

Sexual size and shape dimorphism, and allometric scaling in the pupal and adult traits of Eristalis tenax

The patterns and amount of variation in size, shape, and/or life history traits between females and males are fundamentally important to gain the comprehensive understanding of the evolution of phenotypic diversity. In addition, the covariation of phenotypic traits can significantly contribute to morphological diversification and sexual dimorphism (SD). Using linear and geometric morphometrics, 237 Eristalis tenax specimens sampled from five populations were, therefore, comparatively assessed for the variation in sexual size dimorphism (SSD), sexual shape dimorphism (SShD), and life history traits, as well as for trait covariation (ontogenetic and static allometry). Pupal body, adult wing, and body mass traits were analyzed. Female-biased SSD was observed for pupal length, width, and centroid size, adult wing centroid size, mass, wing loading, and wing area. Conversely, pupal length/width ratio, developmental time, and mass were not found to be sexually dimorphic. Next, wing SShD, but not pupal body SShD was revealed, while allometry was found to be an important "determinant of SD" at the adult stage, with only a minor impact at the pupal stage. By comparing the patterns of covariance (based on allometric slope and intercept) between respective body mass and morphometric traits of pupae and adults, greater variation in allometric slopes was found in adult traits, while static allometries of the two stages significantly differed, as well. Finally, the results indicate that changes in the allometric intercept could be an important source of intraspecific variation and SD in drone fly adults.

Taxonomic complexity in the genus Merodon Meigen, 1803 (Diptera, Syrphidae)

The genus Merodon Meigen, 1803 is distributed across the Palaearctic and Afrotropical Regions. The present work summarizes the knowledge from recent taxonomic and systematic revisions and includes an identification key for the five monophyletic lineages (namely albifrons, aureus, avidus-nigritarsis, desuturinus and natans), 24 species groups, two species subgroups and 10 unplaced species, along with diagnosis and illustrations. A list of 234 taxa, including 194 described and 40 undescribed species, is appended. Most of the species are distributed in the Palaearctic (209 taxa, 181 described), while 27 species (14 described) are known from the Afrotropical Region. Three lineages (aureus, desuturinus and natans) are present in the Afrotropical Region, as well as in the Palaearctic. The Afrotropical melanocerus species group of the desuturinus lineage and the bombiformis species group of the aureus lineage are endemic to the Afrotropical Region, and all other species groups belong to the Palaearctic fauna. The albifrons lineage contains six species groups (albifrons, constans, equestris, geniculatus, ruficornis and rufus) and two unplaced taxa. The aureus lineage includes five species groups (aureus, bombiformis, funestus, nanus and spinitarsis). The avidus-nigritarsis lineage is divided into 10 species groups (aberrans, aurifer, avidus, clavipes, fulcratus, italicus, nigritarsis, pruni, serrulatus and tarsatus) and eight unplaced taxa. The desuturinus lineage contains two species groups: the Afrotropical melanocerus group, with the melanocerus and planifacies subgroups plus the species M. cuthbertsoni Curran, 1939, and the Palaearctic murorum species group. The natans lineage consists of the natans species group plus the species M. segetum Fabricius, 1794.

Comparative Phylogeography and Integrative Taxonomy of Ochlerotatus caspius (Dipera: Culicidae) and Ochlerotatus dorsalis

Given that accurately identifying pathogen vectors is vital for designing efficient mosquito control programs based on the proper surveillance of the epidemiologically important species, it has been suggested the complementary use of independently evolving genes and morphometric traits as a reliable approach for the characterization and delimitation of related species. Hence, we examined the spatial distribution of COI mtDNA and ITS2 rDNA variation from the historical perspective of Ochlerotatus caspius (Pallas, 1771) and O. dorsalis (Meigen, 1830), while simultaneously testing the utility of the two markers in integrative species delimitation when combined with phenotypic character analyses of larvae and adults. Despite the striking difference in haplotype diversity (high in COI mtDNA, low in ITS2 rDNA), no evident phylogeographic structure was apparent in the Palearctic O. caspius. The Holarctic O. dorsalis species was subdivided into two highly distinctive COI mtDNA phylogroups which corresponded to the Nearctic and Palearctic regions. Strong support for the independence of the two allopatric evolutionary lineages suggested that geographical barrier and climatic changes during Pleistocene caused vicariance of the ancestral range. COI mtDNA reliably distinguished O. caspius and O. dorsalis, while ITS2 rDNA yet again lacked the proper resolution for solving this problem. An integrative approach based on the larval and adult morphological traits have varying taxonomic applications due to their differential diagnostic values. Thus, by the implementation of an integrative taxonomic approach, we successfully detected species borders between the two epidemiologically relevant species and uncovered the presence of cryptic diversity within O. dorsalis.

Toward the conservation of the endemic monotypic fish genus Aulopyge from the Balkan Dinaric karst: Integrative assessment of introduced and natural population

The complex biogeographical history of the Balkan Peninsula caused remarkable freshwater fish diversity and endemism, among which Cyprinidae fish dominate. The Dinaric karst was a Pleistocene refugium and it harbors ancient and endemic cyprinids, including Aulopyge huegelii, a sole representative of its genus. Being highly distributionally restricted, it faces various threats that promote a critical decline in population abundance and even population extinction. Phenotypic and molecular diversity of the introduced (Šator Lake, Šator Mountain) and natural (Studena River, Duvanjsko Polje) populations of Dalmatian barbelgudgeon from Bosnia and Herzegovina was studied by using two mitochondrial genes and morphometric traits (linear and geometric morphometrics). Nonparametric ANOVA showed that two analyzed populations significantly differed in six linear measurements, except snout length and postorbital head length. Contrary to centroid size, two populations were found to be significantly different in body shape. Deformation grids indicated that individuals from Studena River are characterized by wider and slightly shorter body comparing to individuals from Šator Lake. Incongruence in cytochrome c oxidase subunit I (COI) and cytochrome b (cyt b) mitochondrial DNA (mtDNA) variation was observed since a common COI haplotype was observed, while four and three cyt b haplotypes were registered in Šator Lake and Studena River, respectively. Since it was demonstrated that cyt b mtDNA was a faster evolving gene, we encourage its use in intraspecies studies, especially for evaluating the connectivity of fragmented populations and for studying the evolutionary footprint of the processes incorporated into the distinctive evolution of Aulopyge. Finally, findings herewith provide a firm basis for designing a long-term sustainable conservation strategy for endemic species in Dinaric karst.

An Integrative Approach Using Phylogenomics and High-Resolution X-Ray Computed Tomography for Species Delimitation in Cryptic Taxa

Morphologically cryptic taxa have proved to be a long-standing challenge for taxonomists. Lineages that show strong genomic structuring across the landscape but are phenotypically similar pose a conundrum, with traditional morphological analyses of these cryptic lineages struggling to keep up with species delimitation advances. Micro X-ray computed tomography (CT) combined with geometric morphometric analyses provides a promising avenue for identification of morphologically cryptic taxa, given its ability to detect subtle differences in anatomical structures. However, this approach has yet to be used in combination with genomic data in a comparative analytical framework to distinguish cryptic taxa. We present an integrative approach incorporating genomic and geometric morphometric evidence to assess the species delimitation of grassland earless dragons (Tympanocryptis spp.) in north-eastern Australia. Using mitochondrial and nuclear genes (ND2 and RAG1, respectively), along with $>$8500 SNPs (nuclear single nucleotide polymorphisms), we assess the evolutionary independence of target lineages and several closely related species. We then integrate phylogenomic data with osteological cranial variation between lineages using landmark-based analyses of three-dimensional CT models. High levels of genomic differentiation between the three target lineages were uncovered, also supported by significant osteological differences. By incorporating multiple lines of evidence, we provide strong support for three undescribed cryptic lineages of Tympanocryptis in north-eastern Australia that warrant taxonomic review. Our approach demonstrates the successful application of CT with integrative taxonomic approaches for cryptic species delimitation, which is broadly applicable across vertebrates containing morphologically similar yet genetically distinct lineages. Additionally, we provide a review of recent integrative taxonomic approaches for cryptic species delimitation and an assessment of how our approach can value-add to taxonomic research.

Checklist of hover flies (Diptera, Syrphidae) of the Republic of Georgia

A checklist of the Syrphidae species of the Republic of Georgia is presented. New hover fly (Diptera: Syrphidae) records from Georgia are provided as a result of field work conducted in 2018. At the same time, published syrphid records for the country are here reviewed and updated. A total of 357 species of hoverflies are now documented from Georgia, 40 of which are reported for the first time. Moreover, DNA barcodes were sequenced for 238 specimens, representing 74 species from this country.

Conflict and congruence between morphological and molecular data: revision of the Merodon constans group (Diptera : Syrphidae)

We revise the Merodon constans species group of the genus Merodon Meigen, 1803 (Diptera: Syrphidae), provide morphological diagnosesand descriptions, as well as an illustrated key and a discussion on the different taxonomic characters used. In total, 15 species were studied, their geographic distributions are presented on maps, and nine new species are described. Two species are redefined and neotypes are designated, lectotypes are designated for five species, and onespeciesis reinstated as valid. Following a detailed study of type material in different entomological collections, the status of several species is revised and three new synonymies are proposed. The M. constans species group was resolved as being monophyletic within the M. albifrons lineage based on molecular analyses using COI and 28S rRNA gene sequences. Three species morphologically similar to M. constans (Rossi, 1794) but occurring outside its distributional rangewere supported as being valid and distinct species on the basis of molecular data, but they were not distinguishable based on morphological characters. By contrast, continental populations of M. analis Meigen, 1822 could not be separated from Mediterranean M. constans based on differences in COI or 28S rRNA genes. The same molecular markers could not discriminate between two other species pairs. We conclude that these molecular markers only partially resolve species within the M. constans group. Geometric morphometry of wing shape successfully separated M. analis and M. constans, as well as M. spineus Vujić, Šašić Zorić & Likov, sp. nov. in both species and population analyses.

The diagnostic utility of sequence-based assays for the molecular delimitation of the epidemiologically relevant Culex pipiens pipiens taxa (Diptera: Culicidae)

The northern house mosquito (Culex pipiens pipiens L.) is a vector of several important pathogens and comprises two epidemiologically distinct ecotypes (molestus Forskål and pipiens). The delimitation of its ecotypes is a crucial, yet controversial step in vector surveillance due to varying diagnostic values of different characters. Therefore, we reviewed the success of a diagnostic assay based on the mitochondrial cytochrome c oxidase subunit I locus (COI) by analyzing previously published sequences of molestus and pipiens sampled in different geographical areas. Next, by genotyping individuals from Northern Serbia at this locus, we additionally assessed whether genetic structure of urban and rural Cx. p. pipiens ecotypes corresponded to the admixture pattern. Finally, to account for the different susceptibility of genetic markers to introgression, we also analyzed genetic structuring based on the ribosomal internal transcribed spacer 2 (ITS2). No latitude-dependent differentiation of Cx. p. pipiens ecotypes was found at a global level, with the COI assay further failing to accurately identify molestus and pipiens ecotypes. Likewise, both individual- (BAPS) and population-based (analysis of molecular variance and FST estimates) methods showed no significant urban/rural genetic differentiation in Serbia, indicating unhindered gene flow between different Cx. p. pipiens habitat types. The findings challenge the previous instances of Cx. p. pipiens ecotype identification, while also spotlighting the vectorial capacity of their hybrid offspring.

The First Molecular and Phenotypic Characterization of the Invasive Population of Aedes albopictus (Diptera: Culicidae) from the Central Balkans

Aedes (Stegomya) albopictus (Skuse 1984), the Asian tiger mosquito, represents the most invasive and one of the medically most important mosquito vectors. Although native to South East Asia, the species has recently spread globally, and was registered in the city of Novi Sad (Serbia, Central Balkans) in August 2018. We characterized the invasive population using phenotypic (wing size and shape) and molecular (nuclear, internal transcribed spacer 2- ITS2, and mitochondrial, cytochrome c oxidase subunit I- COI) markers. The results of phenotypic analyses indicated that the Serbian population could be differentiated from the native (Thailand) and invasive (Hawaii and Florida) populations due to restricted gene flow, founder effect, and supposed different strain origin. The Serbian population showed genetic homogeneity, indicative of a small founder number (bottleneck invasion model). Despite the incorporation of ITS2 GenBank sequences into the data set, neither spatial (Geneland) nor nonspatial (BAPS) genetic structuring analyses helped infer the Serbian population origin. However, the comparison of the retrieved COI haplotype with previously characterized mitogenomes indicated a temperate strain origin, capable of overwintering. Such findings suggest that the newly registered Ae. albopictus population could be able to establish itself since previous studies outlined Novi Sad as a suitable area.

Coexistence of genetically different Varroa destructor in Apis mellifera colonies

The aim of this study was to investigate the genetic diversity of Varroa destructor parasitizing Apis mellifera colonies and to test for possible host-parasite association at the mitochondrial DNA (mtDNA) level. Six A. mellifera haplotypes (including a novel C2aa) and five haplotypes of V. destructor were detected in 29 analyzed colonies from eight sampling sites in Serbia. We revealed the presence of the K and S1 haplotypes as well as KS1 and KP1 heteroplasmic mite individuals in all localities, while the P1 haplotype was only found in four sampling sites. Significant differences in V. destructor genetic diversity were found at both apiary and colony levels, with mite haplotypes coexisting in almost all tested colonies. In addition, a significant correlation between the number of analyzed mites per colony and the number of identified V. destructor haplotypes was observed. However, no significant host-parasite relationship was found, suggesting that mites bearing different haplotypes as well as those heteroplasmic individuals are well adapted to the host, A. mellifera, independently of the identified haplotype present in each colony. Our results will contribute to future population and biogeographic studies concerning V. destructor infesting A. mellifera, as well as to better understanding their host-parasite relationship.

Instant taxonomy: choosing adequate characters for species delimitation and description through congruence between molecular data and quantitative shape analysis

The lack of university funding is one of the major impediments to taxonomy, partly because traditional taxonomic training takes longer than a PhD course. Understanding ranges of phenotypic variability for different morphological structures, and their use as characters for delimitation and description of taxa, is a tedious task. We argue that the advent of molecular barcoding and quantitative shape analysis makes it unnecessary. As an example, we tackle a problematic species-complex of marine copepods from Korea and Japan, approaching it as a starting taxonomist might. Samples were collected from 14 locations and the mitochondrial COI gene was sequenced from 42 specimens. Our phylogenetic analyses reveal four distinct clades in Korea and Japan, and an additional nine belonging to a closely related complex from other parts of the Northern Pacific. Twenty different morphological structures were analysed for one Japanese and two Korean clades using landmark-based two-dimensional geometric morphometrics. Although there is no single morphological character that can distinguish with absolute certainty all three cryptic species, most show statistically significant interspecific differences in shape and size. We use five characters to describe two new species from Korea and to re-describe Tigriopus japonicus Mori, 1938 from near its type locality.

Distribution and predictors of wing shape and size variability in three sister species of solitary bees

Morphological traits can be highly variable over time in a particular geographical area. Different selective pressures shape those traits, which is crucial in evolutionary biology. Among these traits, insect wing morphometry has already been widely used to describe phenotypic variability at the inter-specific level. On the contrary, fewer studies have focused on intra-specific wing morphometric variability. Yet, such investigations are relevant to study potential convergences of variation that could highlight micro-evolutionary processes. The recent sampling and sequencing of three solitary bees of the genus Melitta across their entire species range provides an excellent opportunity to jointly analyse genetic and morphometric variability. In the present study, we first aim to analyse the spatial distribution of the wing shape and centroid size (used as a proxy for body size) variability. Secondly, we aim to test different potential predictors of this variability at both the intra- and inter-population levels, which includes genetic variability, but also geographic locations and distances, elevation, annual mean temperature and precipitation. The comparison of spatial distribution of intra-population morphometric diversity does not reveal any convergent pattern between species, thus undermining the assumption of a potential local and selective adaptation at the population level. Regarding intra-specific wing shape differentiation, our results reveal that some tested predictors, such as geographic and genetic distances, are associated with a significant correlation for some species. However, none of these predictors are systematically identified for the three species as an important factor that could explain the intra-specific morphometric variability. As a conclusion, for the three solitary bee species and at the scale of this study, our results clearly tend to discard the assumption of the existence of a common pattern of intra-specific signal/structure within the intra-specific wing shape and body size variability.

Molecular species delimitation in the genus Eumerus (Diptera: Syrphidae)

Eumerus is one of the most diverse genera of hoverfly worldwide. Species delimitation within genus is considered to be difficult due to: (a) lack of an efficient key; (b) non-defined taxonomical status of a large number of species; and (c) blurred nomenclature. Here, we present the first molecular study to delimit species of the genus by using a fragment of the mitochondrial cytochrome-c oxidase subunit I gene (COI) gene. We assessed 75 specimens assigned to 28 taxa originating from two biogeographic zones: 22 from the western Palaearctic and six from the Afrotropical region. Two datasets were generated based on different sequence lengths to explore the significance of availability of more polymorphic sites for species delimitation; dataset A with a total length of 647 bp and dataset B with 746 bp. Various tree inference approaches and Poisson tree processes models were applied to evaluate the putative 'taxonomical' vs. 'molecular' taxa clusters. All analyses resulted in high taxonomic resolution and clear species delimitation for both the dataset lengths. Furthermore, we revealed a high number of mitochondrial haplotypes and high intraspecific variability. We report two major monophyletic clades, and seven 'molecular' groups of taxa formed, which are congruent with morphology-based taxonomy. Our results support the use of the mitochondrial COI gene in species diagnosis of Eumerus.

Formal nomenclature and description of cryptic species of the Encyrtus sasakii complex (Hymenoptera: Encyrtidae)

With the recent development of molecular approaches to species delimitation, a growing number of cryptic species have been discovered in what had previously been thought to be single morpho-species. Molecular methods, such as DNA barcoding, have greatly enhanced our knowledge of taxonomy, but taxonomy remains incomplete and needs a formal species nomenclature and description to facilitate its use in other scientific fields. A previous study using DNA barcoding, geometric morphometrics and mating tests revealed at least two cryptic species in the Encyrtus sasakii complex. (Hymenoptera: Encyrtidae). To describe these two new species formally (Encyrtus eulecaniumiae sp. nov. and Encyrtus rhodococcusiae sp. nov.), a detailed morphometric study of Encyrtus spp. was performed in addition to the molecular analysis and evaluation of biological data. Morphometric analyses, a multivariate ratio analysis (MRA) and a geometric morphometric analysis (GMA) revealed a great number of differences between the species, but reliable characteristics were not observed for diagnosing the cryptic species. We thus diagnosed these three Encyrtus species on the basis of the characteristics that resulted from genetic markers (mitochondrial cytochrome c oxidase subunit I and nuclear 28S rRNA) and biological data. A formal nomenclature and description of cryptic species was provided on the basis of an integrated taxonomy.

Phylogeographic patterns of Merodon hoverflies in the Eastern Mediterranean region: revealing connections and barriers

We investigated the phylogeographic patterns of Merodon species (Diptera, Syrphidae) in the Eastern Mediterranean. Ten species were sampled on five different islands and mainland sites as a minimum. All samples were screened for their mtDNA COI barcode haplotype diversity, and for some samples, we additionally generated genomic fingerprints. The recently established zoogeographic distribution categories classify these species as having (1) Balkan distribution; (2) Anatolian distribution; (3) continental areas and large islands distribution; and (4) with wide distribution. The ancestral haplotypes and their geographical localities were estimated with statistical parsimony (TCS). TCS networks identified as the ancestral haplotype samples that originated from localities situated within the distributional category of the species in question. Strong geographical haplotype structuring was detected for many Merodon species. We were particularly interested to test the relative importance of current (Aegean Sea) and past Mid-Aegean Trench) barriers to dispersal for Merodon flies in the Aegean. We employed phylogenetic β-diversity (Pβ total) and its partition in replacement (Pβ repl) and richness difference (Pβ rich) to test the importance of each explanatory variable (interisland distance, MAT, and island area) in interisland differences using partial Mantel tests and hierarchical partitioning of variation. β-Analyses confirmed the importance of both current and past barriers to dispersal on the evolution of group. Current interisland distance was particularly important to explain the replacement of haplotypes, while the MAT was driving differences in richness of haplotypes, revealing the MAT as a strong past barrier whose effects are still visible today in the phylogenetic history of the clade in the Aegean. These results support the hypothesis of a highly restricted dispersal and gene flow among Merodon populations between islands since late Pleistocene. Additionally, patterns of phylogeographic structure deduced from haplotype connections and ISSR genome fingerprinting data revealed a few putative cases of human-mediated transfers of Merodon spp.

Constrained body shape among highly genetically divergent allopatric lineages of the supralittoral isopod Ligia occidentalis (Oniscidea)

Multiple highly divergent lineages have been identified within Ligia occidentalis sensu lato, a rocky supralittoral isopod distributed along a ~3000 km latitudinal gradient that encompasses several proposed marine biogeographic provinces and ecoregions in the eastern Pacific. Highly divergent lineages have nonoverlapping geographic distributions, with distributional limits that generally correspond with sharp environmental changes. Crossbreeding experiments suggest postmating reproductive barriers exist among some of them, and surveys of mitochondrial and nuclear gene markers do not show evidence of hybridization. Populations are highly isolated, some of which appear to be very small; thus, the effects of drift are expected to reduce the efficiency of selection. Large genetic divergences among lineages, marked environmental differences in their ranges, reproductive isolation, and/or high isolation of populations may have resulted in morphological differences in L. occidentalis, not detected yet by traditional taxonomy. We used landmark‐based geometric morphometric analyses to test for differences in body shape among highly divergent lineages of L. occidentalis, and among populations within these lineages. We analyzed a total of 492 individuals from 53 coastal localities from the southern California Bight to Central Mexico, including the Gulf of California. We conducted discriminant function analyses (DFAs) on body shape morphometrics to assess morphological variation among genetically differentiated lineages and their populations. We also tested for associations between phylogeny and morphological variation, and whether genetic divergence is correlated to multivariate morphological divergence. We detected significant differences in body shape among highly divergent lineages, and among populations within these lineages. Nonetheless, neither lineages nor populations can be discriminated on the basis of body shape, because correct classification rates of cross‐validated DFAs were low. Genetic distance and phylogeny had weak to no effect on body shape variation. The supralittoral environment appears to exert strong stabilizing selection and/or strong functional constraints on body shape in L. occidentalis, thereby leading to morphological stasis in this isopod.

Cryptic Species or Inadequate Taxonomy? Implementation of 2D Geometric Morphometrics Based on Integumental Organs as Landmarks for Delimitation and Description of Copepod Taxa

Discovery of cryptic species using molecular tools has become common in many animal groups but it is rarely accompanied by morphological revision, creating ongoing problems in taxonomy and conservation. In copepods, cryptic species have been discovered in most groups where fast-evolving molecular markers were employed. In this study at Yeelirrie in Western Australia we investigate a subterranean species complex belonging to the harpacticoid genus Schizopera Sars, 1905, using both the barcoding mitochondrial COI gene and landmark-based two-dimensional geometric morphometrics. Integumental organs (sensilla and pores) are used as landmarks for the first time in any crustacean group. Complete congruence between DNA-based species delimitation and relative position of integumental organs in two independent morphological structures suggests the existence of three distinct evolutionary units. We describe two of them as new species, employing a condensed taxonomic format appropriate for cryptic species. We argue that many supposedly cryptic species might not be cryptic if researchers focus on analyzing morphological structures with multivariate tools that explicitly take into account geometry of the phenotype. A perceived supremacy of molecular methods in detecting cryptic species is in our view a consequence of disparity of investment and unexploited recent advancements in morphometrics among taxonomists. Our study shows that morphometric data alone could be used to find diagnostic morphological traits and gives hope to anyone studying small animals with a hard integument or shell, especially opening the door to assessing fossil diversity and rich museum collections. We expect that simultaneous use of molecular tools with geometry-oriented morphometrics may yield faster formal description of species. Decrypted species in this study are a good example for urgency of formal descriptions, as they display short-range endemism in small groundwater calcrete aquifers in a paleochannel, where their conservation may be threatened by proposed mining.

Species delimitation using morphology, morphometrics, and molecules: definition of the Ophionscutellaris Thomson species group, with descriptions of six new species (Hymenoptera, Ichneumonidae)

The diverse genus Ophion is almost entirely undescribed in the Nearctic region. In this paper we define the Ophionscutellaris species group. This species group is well-supported by analysis of DNA (ITS2, COI, and 28S D2-D3) and morphology. It includes the Palearctic species Ophionscutellaris and the Nearctic species Ophionidoneus. An integrative analysis of DNA, geometric wing morphometrics, classical morphometrics and qualitative morphology indicates that this species group contains a minimum of seven species in North America, although the full diversity of the group has likely not been sampled. Ophionclave Schwarzfeld, sp. n., Ophionaureus Schwarzfeld, sp. n., Ophionbrevipunctatus Schwarzfeld, sp. n., Ophiondombroskii Schwarzfeld, sp. n., Ophionkeala Schwarzfeld, sp. n. and Ophionimportunus Schwarzfeld, sp. n. are described, and a key to the known Nearctic species of the Ophionscutellaris group is provided.

A complex evolutionary history in a remote archipelago: phylogeography and morphometrics of the Hawaiian endemic Ligia isopods

Compared to the striking diversification and levels of endemism observed in many terrestrial groups within the Hawaiian Archipelago, marine invertebrates exhibit remarkably lower rates of endemism and diversification. Supralittoral invertebrates restricted to specific coastal patchy habitats, however, have the potential for high levels of allopatric diversification. This is the case of Ligia isopods endemic to the Hawaiian Archipelago, which most likely arose from a rocky supralittoral ancestor that colonized the archipelago via rafting, and diversified into rocky supralittoral and inland lineages. A previous study on populations of this isopod from O'ahu and Kaua'i revealed high levels of allopatric differentiation, and suggested inter-island historical dispersal events have been rare. To gain a better understanding on the diversity and evolution of this group, we expanded prior phylogeographic work by incorporating populations from unsampled main Hawaiian Islands (Maui, Moloka'i, Lana'i, and Hawai'i), increasing the number of gene markers (four mitochondrial and two nuclear genes), and conducting Maximum likelihood and Bayesian phylogenetic analyses. Our study revealed new lineages and expanded the distribution range of several lineages. The phylogeographic patterns of Ligia in the study area are complex, with Hawai'i, O'ahu, and the Maui-Nui islands sharing major lineages, implying multiple inter-island historical dispersal events. In contrast, the oldest and most geographically distant of the major islands (Kaua'i) shares no lineages with the other islands. Our results did not support the monophyly of all the supralittoral lineages (currently grouped into L. hawaiensis), or the monophyly of the terrestrial lineages (currently grouped into L. perkinsi), implying more than one evolutionary transition between coastal and inland forms. Geometric-morphometric analyses of three supralittoral clades revealed significant body shape differences among them. A taxonomic revision of Hawaiian Ligia is warranted. Our results are relevant for the protection of biodiversity found in an environment subject to high pressure from disturbances.

Wing geometric morphometric inferences on species delimitation and intraspecific divergent units in the Merodon ruficornis group (Diptera, Syrphidae) from the Balkan Peninsula

A study of the Merodon taxa on the Balkan Peninsula, a region with a number of Pleistocene refugia, provides a useful framework for examining evolutionary processes and detecting hidden biodiversity. The phenotypic diversity of 22 samples of the Merodon ruficornis group on the Balkan Peninsula was examined using landmark-based geometric morphometrics. The boundaries of the species M. ruficornis, M. trebevicensis, M. auripes, M. armipes, and M. Ioewi were well defined based on wing shape and size. Canonical variate analysis showed that wing shape possessed sufficient differences to discriminate the species with a successful classification rate of 75-92% for males and 82-100% for females. The observed interspecific differentiation is generally in agreement with a previous study of the M. ruficornis group using a traditional morphological approach and molecular markers (allozyme loci, COI mtDNA). The spatial variability between conspecific populations and interpopulation variation were assessed based on both wing shape and size for male specimens. Phenotypically divergent units were delineated within previously defined species of the M. ruficornis group, indicating the possible presence of evolutionary independent units within the taxa analysed.