High resolution DNA barcode library for European butterflies reveals continental patterns of mitochondrial genetic diversity

A DNA barcode library of Austrian geometridae (Lepidoptera) reveals high potential for DNA-based species identification

Situated in the Eastern section of the European Alps, Austria encompasses a great diversity of different habitat types, ranging from alpine to lowland Pannonian ecosystems, and a correspondingly high level of species diversity, some of which has been addressed in various DNA barcoding projects. Here, we report a DNA barcode library of all the 476 species of Geometridae (Lepidoptera) that have been recorded in Austria. As far as possible, species were sampled from different Austrian regions in order to capture intraspecific genetic variation. In total, 2500 DNA barcode sequences, representing 438 species, were generated in this study. For complete coverage of Austrian geometrid species in the subsequent analyses, the dataset was supplemented with DNA barcodes from specimens of non-Austrian origin. Species delimitations by ASAP, BIN and bPTP methods yielded 465, 510 and 948 molecular operational taxonomic units, respectively. Congruency of BIN and ASAP partitions with morphospecies assignments was reasonably high (85% of morphospecies in unique partitions), whereas bPTP appeared to overestimate the number of taxonomic units. The study furthermore identified taxonomically relevant cases of morphospecies splitting and sharing in the molecular partitions. We conclude that DNA barcoding and sequence analysis revealed a high potential for accurate DNA-based identification of the Austrian Geometridae species. Additionally, the study provides an updated checklist of the geometrid moths of Austria.

Concentration-dependent effect of plant secondary metabolites on bacterial and fungal microbiomes in caterpillar guts

IMPORTANCE

The caterpillar gut is an excellent model system for studying host-microbiome interactions, as it represents an extreme environment for microbial life that usually has low diversity and considerable variability in community composition. Our study design combines feeding caterpillars on a natural and artificial diet with controlled levels of plant secondary metabolites and uses metabarcoding and quantitative PCR to simultaneously profile bacterial and fungal assemblages, which has never been performed. Moreover, we focus on multiple caterpillar species and consider diet breadth. Contrary to many previous studies, our study suggested the functional importance of certain microbial taxa, especially bacteria, and confirmed the previously proposed lower importance of fungi for caterpillar holobiont. Our study revealed the lack of differences between monophagous and polyphagous species in the responses of microbial assemblages to plant secondary metabolites, suggesting the limited role of the microbiome in the plasticity of the herbivore diet.

Principles for Constructing DNA Barcode Reference Libraries

All DNA barcode methods rely on reference sequences linked to well-curated voucher specimens. Definitions for and locations of DNA barcode reference libraries are not standardized, and vary throughout the literature. Standardizing, and centralizing reference specimens would provide an unambiguous source, analogous to reference genomes, to reproduce identifications and improve a library. This chapter proposes a working definition of a DNA barcode reference library, consistent with DNA barcode data standards, along with principles and methods to consider when producing or using such a library. These methods allow explicit traceback to sequence-sources which elevate the value of voucher specimens, and create a potential for community curation.

Taxonomic Structure and Wing Pattern Evolution in the Parnassius mnemosyne Species Complex (Lepidoptera, Papilionidae)

In our study, using the analysis of DNA barcodes and morphology (wing color, male genitalia, and female sphragis shape), we show that the group of species close to P. mnemosyne comprises the western and eastern phylogenetic lineages. The eastern lineage includes P. stubbendorfii, P. glacialis, and P. hoenei. The western lineage includes three morphologically similar species: P. mnemosyne (Western Eurasia), P. turatii (southwestern Europe), and P. nubilosusstat. nov. (Turkmenistan and NE Iran), as well as the morphologically differentiated P. ariadne (Altai). The latter species differs from the rest of the group in the presence of red spots on the wings. Parnassius mnemosyne s.s. is represented by four differentiated mitochondrial clusters that show clear association with specific geographic regions. We propose to interpret them as subspecies: P. mnemosyne mnemosyne (Central and Eastern Europe, N Caucasus, N Turkey), P. mnemosyne adolphi (the Middle East), P. mnemosyne falsa (Tian Shan), and P. mnemosyne gigantea (Gissar-Alai in Central Asia). We demonstrate that in P. ariadne, the red spots on the wing evolved as a reversion to the ancestral wing pattern. This reversion is observed in Altai, where the distribution areas of the western lineage, represented by P. ariadne, and the eastern lineage, represented by P. stubbendorfii, overlap. These two species hybridize in Altai, and we hypothesize that the color change in P. ariadne is the result of reinforcement of prezygotic isolation in the contact zone. The lectotype of Parnassius mnemosyne var. nubilosus Christoph, 1873, is designated.

Genomics of the relict species Baronia brevicornis sheds light on its demographic history and genome size evolution across swallowtail butterflies

Relict species, like coelacanth, gingko, tuatara, are the remnants of formerly more ecologically and taxonomically diverse lineages. It raises the questions of why they are currently species-poor, have restrained ecology, and are often vulnerable to extinction. Estimating heterozygosity level and demographic history can guide our understanding of the evolutionary history and conservation status of relict species. However, few studies have focused on relict invertebrates compared to vertebrates. We sequenced the genome of Baronia brevicornis (Lepidoptera: Papilionidae), which is an endangered species, the sister species of all swallowtail butterflies, and is the oldest lineage of all extant butterflies. From a dried specimen, we were able to generate both long-read and short-read data and assembled a genome of 406 Mb for Baronia. We found a fairly high level of heterozygosity (0.58%) compared to other swallowtail butterflies, which contrasts with its endangered and relict status. Taking into account the high ratio of recombination over mutation, demographic analyses indicated a sharp decline of the effective population size initiated in the last million years. Moreover, the Baronia genome was used to study genome size variation in Papilionidae. Genome sizes are mostly explained by transposable elements activities, suggesting that large genomes appear to be a derived feature in swallowtail butterflies as transposable elements activity is recent and involves different transposable elements classes among species. This first Baronia genome provides a resource for assisting conservation in a flagship and relict insect species as well as for understanding swallowtail genome evolution.

A multilevel social-ecological network approach for reconciling coastal saltmarsh conservation and development

In a large-scale region, governance for connectivity in an ecological system often conflicts with management boundaries, causing inefficiencies. Collaboration among management organizations in different areas can help overcome this problem. However, few studies quantified the collaborations' practical relationship with connectivity, considering that some potentially connected paths are easy to neglect by managers. In this paper, collaborations among government agencies in project application process were analyzed, and a multilevel social-ecological network analysis (SENA) approach was developed to identify the collaboration's effect on genetically connected coastal areas. The network framework and methods were shown in a case of coastal saltmarsh conservation and development in the Yellow River Delta, China. Collaboration patterns in conservation and development networks were analyzed and compared among local, subregional, and regional government agencies working in genetically connected coastal areas. Project information flow, reflecting communication frequency and decision-making chances among government agencies was quantified and correlated with ecological connectivity to inform governance effects. Results showed areas with the potential to realize social-ecological alignment, where collaborative networks were measured by network density (percentage of connected network edges). The current reveals that development has more significant potential than conservation at most levels to overcome the misalignment of the social-ecological system, also known as scale mismatch. Empirical evidence also showed a correlation between communication capacity in development networks and improved ecological conditions. The multilevel SENA advanced in this paper can be used for natural resource management when connectivity plays a major role.

Cryptic Species in Ecotoxicology

The advent of genetic methods has led to the discovery of an increasing number of species that previously could not be distinguished from each other on the basis of morphological characteristics. Even though there has been an exponential growth of publications on cryptic species, such species are rarely considered in ecotoxicology. Thus, the particular question of ecological differentiation and the sensitivity of closely related cryptic species is rarely addressed. Tackling this question, however, is of key importance for evolutionary ecology, conservation biology, and, in particular, regulatory ecotoxicology. At the same time, the use of species with (known or unknown) cryptic diversity might be a reason for the lack of reproducibility of ecotoxicological experiments and implies a false extrapolation of the findings. Our critical review includes a database and literature search through which we investigated how many of the species most frequently used in ecotoxicological assessments show evidence of cryptic diversity. We found a high proportion of reports indicating overlooked species diversity, especially in invertebrates. In terrestrial and aquatic realms, at least 67% and 54% of commonly used species, respectively, were identified as cryptic species complexes. The issue is less prominent in vertebrates, in which we found evidence for cryptic species complexes in 27% of aquatic and 6.7% of terrestrial vertebrates. We further exemplified why different evolutionary histories may significantly determine cryptic species' ecology and sensitivity to pollutants. This in turn may have a major impact on the results of ecotoxicological tests and, consequently, the outcome of environmental risk assessments. Finally, we provide a brief guideline on how to deal practically with cryptic diversity in ecotoxicological studies in general and its implementation in risk assessment procedures in particular. Environ Toxicol Chem 2023;42:1889-1914. © 2023 The Authors. Environmental Toxicology and Chemistry published by Wiley Periodicals LLC on behalf of SETAC.

Global determinants of insect mitochondrial genetic diversity

Understanding global patterns of genetic diversity is essential for describing, monitoring, and preserving life on Earth. To date, efforts to map macrogenetic patterns have been restricted to vertebrates, which comprise only a small fraction of Earth's biodiversity. Here, we construct a global map of predicted insect mitochondrial genetic diversity from cytochrome c oxidase subunit 1 sequences, derived from open data. We calculate the mitochondrial genetic diversity mean and genetic diversity evenness of insect assemblages across the globe, identify their environmental correlates, and make predictions of mitochondrial genetic diversity levels in unsampled areas based on environmental data. Using a large single-locus genetic dataset of over 2 million globally distributed and georeferenced mtDNA sequences, we find that mitochondrial genetic diversity evenness follows a quadratic latitudinal gradient peaking in the subtropics. Both mitochondrial genetic diversity mean and evenness positively correlate with seasonally hot temperatures, as well as climate stability since the last glacial maximum. Our models explain 27.9% and 24.0% of the observed variation in mitochondrial genetic diversity mean and evenness in insects, respectively, making an important step towards understanding global biodiversity patterns in the most diverse animal taxon.

More hidden diversity in a cryptic species complex: a new subspecies of Leptideasinapis (Lepidoptera, Pieridae) from Northern Iran

A new subspecies of Leptideasinapis from Northern Iran, discovered by means of DNA barcoding, is described as Leptideasinapistabarestanassp. nov. The new subspecies is allopatric with respect to other populations of L.sinapis and is genetically distinct, appearing as a well-supported sister clade to all other populations in COI-based phylogenetic reconstructions. Details on karyotype, genitalia, ecology and behaviour for the new subspecies are given and a biogeographical speciation scenario is proposed.

DNA Barcode Library of Megadiverse Lepidoptera in an Alpine Nature Park (Italy) Reveals Unexpected Species Diversity

Species inventories are a prerequisite for biodiversity monitoring and conservation, particularly in protected areas. However, the possibilities of a standardized survey of species diversity using DNA barcoding have so far hardly been implemented, especially in species-rich groups. A first-time molecular-based and nearly complete inventory of the megadiverse insect order Lepidoptera in a protected area in the Alps (Cottian Alps, Italy) was intended to test the possibilities and reliability of DNA-based identifications. From voucher material collected between 2019 and 2022, we successfully sequenced 1213 morphospecies that grouped into 1204 BINs (barcode index numbers), whereas DNA barcoding failed for another 18 species. A total of 35 species shared a BIN with one or more taxa, but a majority of 19 species could still be discriminated by divergent sequences. A total of 12 morphospecies split into two BINs. These species and a further 22 taxa with unique BINs and barcode divergences >2% to the nearest neighbor require taxonomic re-assessment. Two additional cryptic species from the study area were described recently. Finally, 16 species are newly recorded for Italy. Our study, therefore, demonstrates the importance of DNA barcoding for both faunistics and the discovery of cryptic diversity, even in apparently well-studied protected areas.

Genomic evidence for three distinct species in the Erebia manto complex in Central Europe (Lepidoptera, Nymphalidae)

A problem to implement conservation strategies is that in many cases recognized taxa are in fact complexes of several cryptic species. Failure to properly delineate species may lead to misplaced priorities or to inadequate conservation measures. One such species complex is the yellow-spotted ringlet Erebia manto, which comprises several phenotypically distinct lineages, whose degree of genomic isolation has so far not been assessed. Some of these lineages are geographically restricted and thus possibly represent distinct units with conservation priorities. Using several thousand nuclear genomic markers, we evaluated to which degree the bubastis lineage from the Alps and the vogesiaca lineage from the Vosges, are genetically isolated from the widespread manto lineage. Our results suggest that both lineages are genetically as strongly differentiated from manto as other taxonomically well separated sibling species in this genus from each other, supporting a delineation of bubastis and vogesiaca as independent species. Given the restricted and isolated range of vogesiaca as well as the disjunct distribution of bubastis, our findings have significant implication for future conservation efforts on these formerly cryptic species and highlight the need to investigate the genomic identity within species complexes.

Supplementary Information

The online version contains supplementary material available at 10.1007/s10592-023-01501-w.

A DNA barcode library for katydids, cave crickets, and leaf-rolling crickets (Tettigoniidae, Rhaphidophoridae and Gryllacrididae) from Zhejiang Province, China

Barcode libraries are generally assembled with two main objectives in mind: specimen identification and species discovery/delimitation. In this study, the standard COI barcode region was sequenced from 681 specimens belonging to katydids (Tettigoniidae), cave crickets (Rhaphidophoridae), and leaf-rolling crickets (Gryllacrididae) from Zhejiang Province, China. Of these, four COI-5P sequences were excluded from subsequent analyses because they were likely NUMTs (nuclear mitochondrial pseudogenes). The final dataset consisted of 677 barcode sequences representing 90 putative species-level taxa. Automated cluster delineation using the Barcode of Life Data System (BOLD) revealed 118 BINs (Barcodes Index Numbers). Among these 90 species-level taxa, 68 corresponded with morphospecies, while the remaining 22 were identified based on reverse taxonomy using BIN assignment. Thirteen of these morphospecies were represented by a single barcode (so-called singletons), and each of 19 morphospecies were split into more than one BIN. The consensus delimitation scheme yielded 55 Molecular Operational Taxonomic Units (MOTUs). Only four morphospecies (I _max > DNN) failed to be recovered as monophyletic clades (i.e., Elimaeaterminalis, Phyllomimusklapperichi, Sinochloraszechwanensis and Xizicushowardi), so it is speculated that these may be species complexes. Therefore, the diversity of katydids, cave crickets, and leaf-rolling crickets in Zhejiang Province is probably slightly higher than what current taxonomy would suggest.

Analysis of the Holarctic Dictyoptera aurora Complex (Coleoptera, Lycidae) Reveals Hidden Diversity and Geographic Structure in Müllerian Mimicry Ring

The elateroid family Lycidae is known for limited dispersal propensity and high species-level endemism. The red net-winged beetle, Dictyoptera aurora (Herbst, 1874), differs from all relatives by the range comprising almost the entire Holarctic region. Based on a five-marker phylogeny and 67 barcode entries (cox1-5' mtDNA) from the whole range, we recovered two genetically distinct species within traditionally defined D. aurora and resurrected the name D. coccinata (Say, 1835) as the oldest available synonym for Nearctic populations. Yet, no reliable morphological trait distinguishes these species except for minute differences in the male genitalia. D. coccinata is a monophylum resulting from a single Miocene dispersal event, ~15.8 million years ago, and genetic divergence implies long-term isolation by the Bering Strait. Far East Asian and west European populations are also genetically distinct, although to a lower extent. Two independent colonization events established the Fennoscandian populations after the last glacial maximum. Besides intrinsic factors, the high morphological similarity might result from stabilizing selection for shared aposematic signals. The rapidly accumulating barcode data provide valuable information on the evolutionary history and the origins of regional faunas.

One Species, Hundreds of Subspecies? New Insight into the Intraspecific Classification of the Old World Swallowtail (Papilio machaon Linnaeus, 1758) Based on Two Mitochondrial DNA Markers

The Old World swallowtail Papilio machaon Linnaeus, 1758 is one of the most well-known and most characteristic members of the family Papilionidae. Over the past two centuries, the butterfly has been the subject of many studies. P. machaon is characterised by a tendency to change the wing colour pattern. In turn, due to the great interest of collectors and amateur entomologists, these studies have been converted into the description of over 100 colour forms, aberrations and subspecies. In this study, mitochondrial DNA (mtDNA), 16S rDNA and cytochrome b sequences were used to examine the correlation between the intraspecific classification and genetic structure of P. machaon. The study used 87 specimens from 59 different localities covering the geographic distribution of this species in the Palaearctic. The phylogenetic relationships within and between the Old World swallowtail subspecies showed that the intraspecific classification proposed by various authors does not correlate with the variability in mitochondrial DNA sequences. In addition, populations occurring at the species distribution borders in the Palaearctic Region (i.e., Japan, Kamchatka, Morocco and Sakhalin) are genetically distinct from other species.

Butterfly–parasitoid–hostplant interactions in Western Palaearctic Hesperiidae: a DNA barcoding reference library

The study of ecological interactions between plants, phytophagous insects and their natural enemies is an essential but challenging component for understanding ecosystem dynamics. Molecular methods such as DNA barcoding can help elucidate these interactions. In this study, we employed DNA barcoding to establish hostplant and parasitoid interactions with hesperiid butterflies, using a complete reference library for Hesperiidae of continental Europe and north-western Africa (53 species, 100% of those recorded) based on 2934 sequences from 38 countries. A total of 233 hostplant and parasitoid interactions are presented, some recovered by DNA barcoding larval remains or parasitoid cocoons. Combining DNA barcode results with other lines of evidence allowed 94% species-level identification for Hesperiidae, but success was lower for parasitoids, in part due to unresolved taxonomy. Potential cases of cryptic diversity, both in Hesperiidae and Microgastrinae, are discussed. We briefly analyse the resulting interaction networks. Future DNA barcoding initiatives in this region should focus attention on north-western Africa and on parasitoids, because in these cases barcode reference libraries and taxonomy are less well developed.

A DNA barcode reference library for endemic Ponto-Caspian amphipods

The Ponto-Caspian region is an endemicity hotspot that harbours several crustacean radiations, among which amphipods are the most diverse. These poorly known species are severely threatened in their native range, while at the same time they are invading European inland waters with significant ecological consequences. A proper taxonomic knowledge of this fauna is paramount for its conservation within the native region and monitoring outside of it. Here, we assemble a DNA barcode reference library for nearly 60% of all known Ponto-Caspian amphipod species. We use several methods to define molecular operational taxonomic units (MOTUs), based on two mitochondrial markers (COI and 16S), and assess their congruence with current species-level taxonomy based on morphology. Depending on the method, we find that 54–69% of species had congruent morpho-molecular boundaries. The cases of incongruence resulted from lumping distinct morphospecies into a single MOTU (7–27%), splitting a morphospecies into several MOTUs (4–28%), or both (4–11%). MOTUs defined by distance-based methods without a priori divergence thresholds showed the highest congruence with morphological taxonomy. These results indicate that DNA barcoding is valuable for clarifying the diversity of Ponto-Caspian amphipods, but reveals that extensive work is needed to resolve taxonomic uncertainties. Our study advances the DNA barcode reference library for the European aquatic biota, paving the way towards improved taxonomic knowledge needed to enhance monitoring and conservation efforts.

Underestimated cryptic diversity in the Caryocolumtricolorella species complex (Lepidoptera, Gelechiidae)

The taxonomy of the Caryocolumtricolorella species complex, an informal subsection of the diverse Caryocoluminteralbicella species group, is revised and four species are separated from DNA barcodes of the mitochondrial COI (cytochrome c oxidase subunit 1) gene and adult morphology: C.tricolorella (Haworth, 1812), C.fibigerium Huemer, 1988, C.herwigvanstaai sp. nov., and C.olekarsholti sp. nov. These species show a vicariant distribution pattern, with C.tricolorella widely distributed in Central and Northern Europe, C.fibigerium restricted to the Iberian Peninsula and southern France, C.herwigvanstaai sp. nov. to the Italian Peninsula, and C.olekarsholti sp. nov. to the Balkans. All species are described in detail, and the adults and genitalia of both sexes are illustrated.

Pronounced mito-nuclear discordance and various Wolbachia infections in the water ringlet Erebia pronoe have resulted in a complex phylogeographic structure

Several morphological and mitochondrial lineages of the alpine ringlet butterfly species Erebia pronoe have been described, indicating a complex phylogenetic structure. However, the existing data were insufficient and allow neither a reconstruction of the biogeographic history, nor an assessment of the genetic lineages. Therefore, we analysed mitochondrial (COI, NDI) and nuclear (EF1α, RPS5) gene sequences and compared them with sequences from the sister species Erebia melas. Additionally, we combined this information with morphometric data of the male genitalia and the infection patterns with Wolbachia strains, based on a WSP analysis. We obtained a distinct phylogeographic structure within the E. pronoe-melas complex with eight well-distinguishable geographic groups, but also a remarkable mito-nuclear discordance. The mito-nuclear discordance in E. melas and E. pronoe glottis can be explained by different ages of Wolbachia infections with different Wolbachia strains, associated selective sweeps, and hybridisation inhibition. Additionally, we found indications for incipient speciation of E. pronoe glottis in the Pyrenees and a pronounced range dynamic within and among the other high mountain systems of Europe. Our results emphasize the importance of combined approaches in reconstructing biogeographic patterns and evaluating phylogeographic splits.

Hybridization fuelled diversification in Spialia butterflies

The importance of hybridization and introgression is well documented in the evolution of plants but, in insects, their role is not fully understood. Given the fact that insects are the most diverse group of organisms, assessing the impact of reticulation events on their evolution may be key to comprehend the emergence of such remarkable diversity. Here, we used an insect model, the Spialia butterflies, to gather genomic evidence of hybridization as a promoter of novel diversity. By using double‐digest RADseq (ddRADseq), we explored the phylogenetic relationships between Spialia orbifer, S. rosae and S. sertorius, and documented two independent events of interspecific gene flow. Our data support that the Iberian endemism S. rosae probably received genetic material from S. orbifer in both mitochondrial and nuclear DNA, which could have contributed to a shift in the ecological preferences of S. rosae. We also show that admixture between S. sertorius and S. orbifer probably occurred in Italy. As a result, the admixed Sicilian populations of S. orbifer are differentiated from the rest of populations both genetically and morphologically, and display signatures of reproductive character displacement in the male genitalia. Additionally, our analyses indicated that genetic material from S. orbifer is present in S. sertorius along the Italian Peninsula. Our findings add to the view that hybridization is a pervasive phenomenon in nature and in butterflies in particular, with important consequences for evolution due to the emergence of novel phenotypes.

Local adaptation to climate anomalies relates to species phylogeny

Climatic anomalies are increasing in intensity and frequency due to rapid rates of global change, leading to increased extinction risk for many species. The impacts of anomalies are likely to vary between species due to different degrees of sensitivity and extents of local adaptation. Here, we used long-term butterfly monitoring data of 143 species across six European bioclimatic regions to show how species’ population dynamics have responded to local or globally-calculated climatic anomalies, and how species attributes mediate these responses. Contrary to expectations, degree of apparent local adaptation, estimated from the relative population sensitivity to local versus global anomalies, showed no associations with species mobility or reproductive rate but did contain a strong phylogenetic signal. The existence of phylogenetically-patterned local adaptation to climate has important implications for forecasting species responses to current and future climatic conditions and for developing appropriate conservation practices.

Melero et al. investigate butterfly responses to climatic anomalies from long-term monitoring observations in the field. They found the degree of adaptation to local fluctuations in climate had a strong phylogenetic signal but was not associated with mobility or reproductive rate of a species.

Genomics Reveal Admixture and Unexpected Patterns of Diversity in a Parapatric Pair of Butterflies

We studied the evolutionary relationship of two widely distributed parapatric butterfly species, Melitaea athalia and Melitaea celadussa, using the ddRAD sequencing approach, as well as genital morphology and mtDNA data. M. athalia was retrieved as paraphyletic with respect to M. celadussa. Several cases of mito-nuclear discordance and morpho-genetic mismatch were found in the contact zone. A strongly diverged and marginally sympatric clade of M. athalia from the Balkans was revealed. An in-depth analysis of genomic structure detected high levels of admixture between M. athalia and M. celadussa at the contact zone, though not reaching the Balkan clade. The demographic modelling of populations supported the intermediate genetic make-up of European M. athalia populations with regards to M. celadussa and the Balkan clade. However, the dissimilarity matrix of genotype data (PCoA) suggested the Balkan lineage having a genetic component that is unrelated to the athalia-celadussa group. Although narrowly sympatric, almost no signs of gene flow were found between the main M. athalia group and the Balkan clade. We propose two possible scenarios on the historical evolution of our model taxa and the role of the last glacial maximum in shaping their current distribution. Finally, we discuss the complexities regarding the taxonomic delimitation of parapatric taxa.

Demographic inferences and climatic niche modelling shed light on the evolutionary history of the emblematic cold-adapted Apollo butterfly at regional scale

Cold-adapted species escape climate warming by latitudinal and/or altitudinal range shifts, and currently occur in Southern Europe in isolated mountain ranges within "sky islands". Here, we studied the genetic structure of the Apollo butterfly in five such sky islands (above 1,000 m) in France, and infer its demographic history since the last interglacial, using single nucleotide polymorphisms (ddRADseq SNPs). The Auvergne and Alps populations show strong genetic differentiation but not alpine massifs, although separated by deep valleys. Combining three complementary demographic inference methods and species distribution models (SDMs) we show that the LIG period was highly unfavourable for Apollo that probably survived in small population in the highest summits of Auvergne. The population shifted downslope and expanded eastward between LIG and LGM throughout the large climatically suitable Rhône valley between the glaciated summits of Auvergne and Alps. The Auvergne and Alps populations started diverging before the LGM but remained largely connected till the mid-Holocene. Population decline in Auvergne was more gradual but started before (~7 kya vs. 800 ya), and was much stronger with current population size ten times lower than in the Alps. In the Alps, the low genetic structure and limited evidence for isolation by distance suggest a nonequilibrium metapopulation functioning. The core Apollo population experienced cycles of contraction-expansion with climate fluctuations with largely interconnected populations overtime according to a "metapopulation-pulsar" functioning. This study demonstrates the power of combining demographic inferences and SDMs to determine past and future evolutionary trajectories of an endangered species at a regional scale.

A Meta-Analytical Investigation of the Gap between Measured and Predicted Inter-Population Genetic Diversity in Species of High Conservation Concern-The Case of the Critically Endangered European Mink Mustela lutreola L., 1761

Although properly designed sampling in population genetic studies is of key importance for planning evidence-informed conservation measures, sampling strategies are rarely discussed. This is the case for the European mink Mustela lutreola, a critically endangered species. In order to address this problem, a meta-analysis aiming to examine the completeness of mtDNA haplotype sampling in recent studies of M. lutreola inter-population genetic diversity was conducted. The analysis was performed using the sample-size-based rarefaction and extrapolation sampling curve method for three populations—the Northeastern (Russia, Belarus and Estonia), the Western (France and Spain), and the Southeastern (Romania). The extrapolated values of the Shannon–Wiener index were determined, assuming full sample coverage. The gap between the measured and predicted inter-population genetic diversity was estimated, indicating that the identified level of sample coverage was the lowest for the NE population (87%), followed by the SE population (96%) and the W population (99%). A guide for sampling design and accounting for sampling uncertainty in future population genetic studies on European mink is provided. The relatively low sample coverage for the Russian population clearly indicates an urgent need to take conservation measures for European mink in this country.

A molecular-based identification resource for the arthropods of Finland

To associate specimens identified by molecular characters to other biological knowledge, we need reference sequences annotated by Linnaean taxonomy. In this study, we (1) report the creation of a comprehensive reference library of DNA barcodes for the arthropods of an entire country (Finland), (2) publish this library, and (3) deliver a new identification tool for insects and spiders, as based on this resource. The reference library contains mtDNA COI barcodes for 11,275 (43%) of 26,437 arthropod species known from Finland, including 10,811 (45%) of 23,956 insect species. To quantify the improvement in identification accuracy enabled by the current reference library, we ran 1000 Finnish insect and spider species through the Barcode of Life Data system (BOLD) identification engine. Of these, 91% were correctly assigned to a unique species when compared to the new reference library alone, 85% were correctly identified when compared to BOLD with the new material included, and 75% with the new material excluded. To capitalize on this resource, we used the new reference material to train a probabilistic taxonomic assignment tool, FinPROTAX, scoring high success. For the full-length barcode region, the accuracy of taxonomic assignments at the level of classes, orders, families, subfamilies, tribes, genera, and species reached 99.9%, 99.9%, 99.8%, 99.7%, 99.4%, 96.8%, and 88.5%, respectively. The FinBOL arthropod reference library and FinPROTAX are available through the Finnish Biodiversity Information Facility (www.laji.fi) at https://laji.fi/en/theme/protax. Overall, the FinBOL investment represents a massive capacity-transfer from the taxonomic community of Finland to all sectors of society.

Evidence for multiple colonisations and Wolbachia infections shaping the genetic structure of the widespread butterfly Polyommatus icarus in the British Isles

The paradigm of isolation in southern refugia during glacial periods followed by expansions during interglacials, producing limited genetic differentiation in northern areas, dominates European phylogeography. However, the existence of complex structured populations in formerly glaciated areas, and islands connected to mainland areas during glacial maxima, call for alternative explanations. We reconstructed the mtDNA phylogeography of the widespread Polyommatus Icarus butterfly with an emphasis on the formerly glaciated and connected British Isles. We found distinct geographical structuring of CO1 haplogroups, with an ancient lineage restricted to the marginal European areas, including Northern Scotland and Outer Hebrides. Population genomic analyses, using ddRADSeq genomic markers, also reveal substantial genetic structuring within Britain. However, there is negligble mito-nuclear concordance consistent with independent demographic histories of mitochondrial versus nuclear DNA. While mtDNA-Wolbachia associations in northern Britain could account for the geographic structuring of mtDNA across most of the British Isles, for nuclear DNA markers (derived from ddRADseq data) butterflies from France cluster between northern and southern British populations - an observation consistent with a scenario of multiple recolonisation. Taken together our results suggest that contemporary mtDNA structuring in the British Isles (and potentially elsewhere in Europe) largely results from Wolbachia infections, however, nuclear genomic structuring suggests a history of at least two distinct colonisations. This two-stage colonisation scenario has previously been put forth to explain genetic diversity and structuring in other British flora and fauna. Additionally, we also present preliminary evidence for potential Wolbachia-induced feminization in the Outer Hebrides.

The Pleistocene species pump past its prime: Evidence from European butterfly sister species

The Pleistocene glacial cycles had a profound impact on the ranges and genetic make-up of organisms. While it is clear that the contact zones that have been described for many sister taxa are secondary and have formed in the current interglacial, it is unclear when the taxa involved began to diverge. Previous estimates based on small numbers of loci are unreliable given the stochasticity of genetic drift and the contrasting effects of incomplete lineage sorting and gene flow on gene divergence. Here, we use genome-wide transcriptome data to estimate divergence for 18 sister species pairs of European butterflies showing either sympatric or contact zone distributions. We find that in most cases, species divergence predates the mid-Pleistocene transition or even the entire Pleistocene period. We also show that although post-divergence gene flow is restricted to contact zone pairs, they are not systematically younger than sympatric pairs. This suggests that contact zones are not limited to the initial stages of the speciation process, but can involve notably old taxa. Finally, we show that mitochondrial divergence and nuclear divergence are only weakly correlated and mitochondrial divergence is higher for contact zone pairs.

A DNA barcode library for the butterflies of North America

Although the butterflies of North America have received considerable taxonomic attention, overlooked species and instances of hybridization continue to be revealed. The present study assembles a DNA barcode reference library for this fauna to identify groups whose patterns of sequence variation suggest the need for further taxonomic study. Based on 14,626 records from 814 species, DNA barcodes were obtained for 96% of the fauna. The maximum intraspecific distance averaged 1/4 the minimum distance to the nearest neighbor, producing a barcode gap in 76% of the species. Most species (80%) were monophyletic, the others were para- or polyphyletic. Although 15% of currently recognized species shared barcodes, the incidence of such taxa was far higher in regions exposed to Pleistocene glaciations than in those that were ice-free. Nearly 10% of species displayed high intraspecific variation (>2.5%), suggesting the need for further investigation to assess potential cryptic diversity. Aside from aiding the identification of all life stages of North American butterflies, the reference library has provided new perspectives on the incidence of both cryptic and potentially over-split species, setting the stage for future studies that can further explore the evolutionary dynamics of this group.