The species delimitation problem applied to the Agabus bipustulatus complex (Coleoptera, Dytiscidae) in north Scandinavia

Genomic, morphological and physiological data support fast ecotypic differentiation and incipient speciation in an alpine diving beetle

An intricate interplay between evolutionary and demographic processes has frequently resulted in complex patterns of genetic and phenotypic diversity in alpine lineages, posing serious challenges to species delimitation and biodiversity conservation planning. Here we integrate genomic data, geometric morphometric analyses and thermal tolerance experiments to explore the role of Pleistocene climatic changes and adaptation to alpine environments on patterns of genomic and phenotypic variation in diving beetles from the taxonomically complex Agabus bipustulatus species group. Genetic structure and phylogenomic analyses revealed the presence of three geographically cohesive lineages, two representing trans-Palearctic and Iberian populations of the elevation-generalist A. bipustulatus and another corresponding to the strictly-alpine A. nevadensis, a narrow-range endemic taxon from the Sierra Nevada mountain range in southeastern Iberia. The best-supported model of lineage divergence, along with the existence of pervasive genetic introgression and admixture in secondary contact zones, is consistent with a scenario of population isolation and connectivity linked to Quaternary climatic oscillations. Our results suggest that A. nevadensis is an alpine ecotype of A. bipustulatus, whose genotypic, morphological and physiological differentiation likely resulted from an interplay between population isolation and local altitudinal adaptation. Remarkably, within the Iberian Peninsula, such ecotypic differentiation is unique to Sierra Nevada populations and has not been replicated in other alpine populations of A. bipustulatus. Collectively, our study supports fast ecotypic differentiation and incipient speciation processes within the study complex and points to Pleistocene glaciations and local adaptation along elevational gradients as key drivers of biodiversity generation in alpine environments.

Species delimitation in northern European water scavenger beetles of the genus Hydrobius (Coleoptera, Hydrophilidae)

The chiefly Holarctic Hydrobius species complex (Coleoptera, Hydrophilidae) currently consists of Hydrobius arcticus Kuwert, 1890, and three morphological variants of Hydrobius fuscipes (Linnaeus, 1758): var. fuscipes, var. rottenbergii and var. subrotundus in northern Europe. Here molecular and morphological data are used to test the species boundaries in this species complex. Three gene segments (COI, H3 and ITS2) were sequenced and analyzed with Bayesian methods to infer phylogenetic relationships. The Generalized Mixed Yule Coalescent (GMYC) model and two versions of the Bayesian species delimitation method BPP, with or without an a priori defined guide tree (v2.2 & v3.0), were used to evaluate species limits. External and male genital characters of primarily Fennoscandian specimens were measured and statistically analyzed to test for significant differences in quantitative morphological characters. The four morphotypes formed separate genetic clusters on gene trees and were delimited as separate species by GMYC and by both versions of BPP, despite specimens of Hydrobius fuscipes var. fuscipes and Hydrobius fuscipes var. subrotundus being sympatric. Hydrobius arcticus and Hydrobius fuscipes var. rottenbergii could only be separated genetically with ITS2, and were delimited statistically with GMYC on ITS2 and with BPP on the combined data. In addition, six or seven potentially cryptic species of the Hydrobius fuscipes complex from regions outside northern Europe were delimited genetically. Although some overlap was found, the mean values of six male genital characters were significantly different between the morphotypes (p < 0.001). Morphological characters previously presumed to be diagnostic were less reliable to separate Hydrobius fuscipes var. fuscipes from Hydrobius fuscipes var. subrotundus, but characters in the literature for Hydrobius arcticus and Hydrobius fuscipes var. rottenbergii were diagnostic. Overall, morphological and molecular evidence strongly suggest that Hydrobius arcticus and the three morphological variants of Hydrobius fuscipes are separate species and Hydrobius rottenbergii Gerhardt, 1872, stat. n. and Hydrobius subrotundus Stephens, 1829, stat. n. are elevated to valid species. An identification key to northern European species of Hydrobius is provided.

Changing Names with Changed Address: Integrated Taxonomy and Species Delimitation in the Holarctic Colymbetes paykulli Group (Coleoptera: Dytiscidae)

Species delimitation of geographically isolated forms is a long-standing problem in less studied insect groups. Often taxonomic decisions are based directly on morphologic variation, and lack a discussion regarding sample size and the efficiency of migration barriers or dispersal/migration capacity of the studied species. These problems are here exemplified in a water beetle complex from the Bering Sea region that separates North America from Eurasia. Only a few sampled specimens occur from this particular area and they are mostly found in museum and private collections. Here we utilize the theory of integrated taxonomy to discuss the speciation of the Holarctic Colymbetes paykulli water beetle complex, which historically has included up to five species of which today only two are recognized. Three delimitation methods are used; landmark based morphometry of body shape, variation in reticulation patterns of the pronotum exo-skeleton and sequence variation of the partial mitochondrial gene Cyt b. Our conclusion is that the Palearctic and Nearctic populations of C. paykulli are given the status of separate species, based on the fact that all methods showed significant separation between populations. As a consequence the name of the Palearctic species is C. paykulli Erichson and the Nearctic species should be known as C. longulus LeConte. There is no clear support for delineation between Palearctic and Nearctic populations of C. dahuricus based on mtDNA. However, significant difference in size and reticulation patterns from the two regions is shown. The combined conclusion is that the C. dahuricus complex needs a more thorough investigation to fully disentangle its taxonomic status. Therefore it is here still regarded as a Holarctic species. This study highlights the importance to study several diagnosable characters that has the potential to discriminate evolutionary lineage during speciation.

Widespread phenotypic and genetic divergence along altitudinal gradients in animals

Altitudinal gradients offer valuable study systems to investigate how adaptive genetic diversity is distributed within and between natural populations and which factors promote or prevent adaptive differentiation. The environmental clines along altitudinal gradients tend to be steep relative to the dispersal distance of many organisms, providing an opportunity to study the joint effects of divergent natural selection and gene flow. Temperature is one variable showing consistent altitudinal changes, and altitudinal gradients can therefore provide spatial surrogates for some of the changes anticipated under climate change. Here, we investigate the extent and patterns of adaptive divergence in animal populations along altitudinal gradients by surveying the literature for (i) studies on phenotypic variation assessed under common garden or reciprocal transplant designs and (ii) studies looking for signatures of divergent selection at the molecular level. Phenotypic data show that significant between-population differences are common and taxonomically widespread, involving traits such as mass, wing size, tolerance to thermal extremes and melanization. Several lines of evidence suggest that some of the observed differences are adaptively relevant, but rigorous tests of local adaptation or the link between specific phenotypes and fitness are sorely lacking. Evidence for a role of altitudinal adaptation also exists for a number of candidate genes, most prominently haemoglobin, and for anonymous molecular markers. Novel genomic approaches may provide valuable tools for studying adaptive diversity, also in species that are not amenable to experimentation.

Ecotype differentiation in the face of gene flow within the diving beetle Agabus bipustulatus (Linnaeus, 1767) in northern Scandinavia

The repeated occurrence of habitat-specific polyphyletic evolved ecotypes throughout the ranges of widely distributed species implies that multiple, independent and parallel selection events have taken place. Ecological transitions across altitudinal gradients over short geographical distances are often associated with variation in habitat-related fitness, these patterns suggest the action of strong selective forces. Genetic markers will therefore contribute differently to differences between ecotypes in local hybrid zones. Here we have studied the adaptive divergence between ecotypes of the water beetle Agabus bipustulatus along several parallel altitudinal gradients in northern Scandinavia. This water beetle is well known for its remarkable morphological variation associated with mountain regions throughout the western Palaearctic. Two morphological ecotypes are recognised: a montane type with reduced flight muscles and a lowland type with fully developed muscles. Using a multilocus survey of allozyme variation and a morphological analysis with landmark-based morphometrics, across thirty-three populations and seven altitudinal gradients, we studied the local adaptive process of gene flow and selection in detail. Populations were sampled at three different elevations: below, at and above the tree line. The results indicate that the levels of divergence observed between ecotypes in morphology and allele frequencies at α-Glycerophosphate dehydrogenase relative to those shown by neutral molecular markers reflects local diversifying selection in situ. Four main lines of evidence are shown here: (1) A repeated morphological pattern of differentiation is observed across all altitudinal transects, with high reclassification probabilities. (2) Allele and genotype frequencies at the α-Gpdh locus are strongly correlated with altitude, in sharp contrast to the presumable neutral markers. (3) Genetic differentiation is two to three times higher among populations across the tree line than among populations at or below. (4) Genetic differentiation between ecotypes within independent mountain areas is reflected by different sets of allozymes.