Complete DNA barcode reference library for a country's butterfly fauna reveals high performance for temperate Europe

Looking back on a decade of barcoding crustaceans

Species identification represents a pivotal component for large-scale biodiversity studies and conservation planning but represents a challenge for many taxa when using morphological traits only. Consequently, alternative identification methods based on molecular markers have been proposed. In this context, DNA barcoding has become a popular and accepted method for the identification of unknown animals across all life stages by comparison to a reference library. In this review we examine the progress of barcoding studies for the Crustacea using the Web of Science data base from 2003 to 2014. All references were classified in terms of taxonomy covered, subject area (identification/library, genetic variability, species descriptions, phylogenetics, methods, pseudogenes/numts), habitat, geographical area, authors, journals, citations, and the use of the Barcode of Life Data Systems (BOLD). Our analysis revealed a total number of 164 barcoding studies for crustaceans with a preference for malacostracan crustaceans, in particular Decapoda, and for building reference libraries in order to identify organisms. So far, BOLD did not establish itself as a popular informatics platform among carcinologists although it offers many advantages for standardized data storage, analyses and publication.

Correlation between the green-island phenotype and Wolbachia infections during the evolutionary diversification of Gracillariidae leaf-mining moths

Internally feeding herbivorous insects such as leaf miners have developed the ability to manipulate the physiology of their host plants in a way to best meet their metabolic needs and compensate for variation in food nutritional composition. For instance, some leaf miners can induce green-islands on yellow leaves in autumn, which are characterized by photosynthetically active green patches in otherwise senescing leaves. It has been shown that endosymbionts, and most likely bacteria of the genus Wolbachia, play an important role in green-island induction in the apple leaf-mining moth Phyllonorycter blancardella. However, it is currently not known how widespread is this moth-Wolbachia-plant interaction. Here, we studied the co-occurrence between Wolbachia and the green-island phenotype in 133 moth specimens belonging to 74 species of Lepidoptera including 60 Gracillariidae leaf miners. Using a combination of molecular phylogenies and ecological data (occurrence of green-islands), we show that the acquisitions of the green-island phenotype and Wolbachia infections have been associated through the evolutionary diversification of Gracillariidae. We also found intraspecific variability in both green-island formation and Wolbachia infection, with some species being able to form green-islands without being infected by Wolbachia. In addition, Wolbachia variants belonging to both A and B supergroups were found to be associated with green-island phenotype suggesting several independent origins of green-island induction. This study opens new prospects and raises new questions about the ecology and evolution of the tripartite association between Wolbachia, leaf miners, and their host plants.

Macaria mirthae Vargas et al (Lepidoptera: Geometridae): Confirmation of the Use of an Invasive Host Plant in the Northern Atacama Desert of Chile Based on DNA Barcodes

Macaria mirthae Vargas et al (Lepidoptera: Geometridae) is a geometrid moth native to the northern Atacama Desert of Chile. Its oligophagous larvae are associated with native hosts of the plant family Fabaceae, the most important of which is Acacia macracantha. The invasive tree Leucaena leucocephala (Fabaceae) was recently recorded as a host plant for M. mirthae based on morphology. The taxonomic status of larvae collected on A. macracantha and L. leucocephala was assessed using sequences of the DNA barcode fragment of the cytochrome c oxidase subunit I (COI) gene. Genetic divergence between samples from the host plants was found to be 0%-0.8% (Kimura 2-parameter model). Neighbor-joining and maximum likelihood analyses were also performed, including additional barcode sequences of Neotropical geometrid moths from GenBank and BOLD databases. Sequences of the larvae from both host plants clustered in a single clade with high statistical support in both analyses. Based on these results, it is concluded that M. mirthae has effectively expanded its host range and its larvae are currently feeding on the exotic tree L. leucocephala. Additionally, the importance of this new host association in a highly disturbed habitat is briefly discussed in terms of the field biology of this native geometrid moth.

DNA barcode reference library for Iberian butterflies enables a continental-scale preview of potential cryptic diversity

How common are cryptic species--those overlooked because of their morphological similarity? Despite its wide-ranging implications for biology and conservation, the answer remains open to debate. Butterflies constitute the best-studied invertebrates, playing a similar role as birds do in providing models for vertebrate biology. An accurate assessment of cryptic diversity in this emblematic group requires meticulous case-by-case assessments, but a preview to highlight cases of particular interest will help to direct future studies. We present a survey of mitochondrial genetic diversity for the butterfly fauna of the Iberian Peninsula with unprecedented resolution (3502 DNA barcodes for all 228 species), creating a reliable system for DNA-based identification and for the detection of overlooked diversity. After compiling available data for European butterflies (5782 sequences, 299 species), we applied the Generalized Mixed Yule-Coalescent model to explore potential cryptic diversity at a continental scale. The results indicate that 27.7% of these species include from two to four evolutionary significant units (ESUs), suggesting that cryptic biodiversity may be higher than expected for one of the best-studied invertebrate groups and regions. The ESUs represent important units for conservation, models for studies of evolutionary and speciation processes, and sentinels for future research to unveil hidden diversity.

Phylogenetic molecular species delimitations unravel potential new species in the pest genus Spodoptera Guenée, 1852 (Lepidoptera, Noctuidae)

Nowadays molecular species delimitation methods promote the identification of species boundaries within complex taxonomic groups by adopting innovative species concepts and theories (e.g. branching patterns, coalescence). As some of them can efficiently deal with large single-locus datasets, they could speed up the process of species discovery compared to more time consuming molecular methods, and benefit from the existence of large public datasets; these methods can also particularly favour scientific research and actions dealing with threatened or economically important taxa. In this study we aim to investigate and clarify the status of economically important moths species belonging to the genus Spodoptera (Lepidoptera, Noctuidae), a complex group in which previous phylogenetic analyses and integrative approaches already suggested the possible occurrence of cryptic species and taxonomic ambiguities. In this work, the effectiveness of innovative (and faster) species delimitation approaches to infer putative species boundaries has been successfully tested in Spodoptera, by processing the most comprehensive dataset (in terms of number of species and specimens) ever achieved; results are congruent and reliable, irrespective of the set of parameters and phylogenetic models applied. Our analyses confirm the existence of three potential new species clusters (for S. exigua (Hübner, 1808), S. frugiperda (J.E. Smith, 1797) and S. mauritia (Boisduval, 1833)) and support the synonymy of S. marima (Schaus, 1904) with S. ornithogalli (Guenée, 1852). They also highlight the ambiguity of the status of S. cosmiodes (Walker, 1858) and S. descoinsi Lalanne-Cassou & Silvain, 1994. This case study highlights the interest of molecular species delimitation methods as valuable tools for species discovery and to emphasize taxonomic ambiguities.

One species in eight: DNA barcodes from type specimens resolve a taxonomic quagmire

Each holotype specimen provides the only objective link to a particular Linnean binomen. Sequence information from them is increasingly valuable due to the growing usage of DNA barcodes in taxonomy. As type specimens are often old, it may only be possible to recover fragmentary sequence information from them. We tested the efficacy of short sequences from type specimens in the resolution of a challenging taxonomic puzzle: the Elachista dispunctella complex which includes 64 described species with minuscule morphological differences. We applied a multistep procedure to resolve the taxonomy of this species complex. First, we sequenced a large number of newly collected specimens and as many holotypes as possible. Second, we used all >400 bp examine species boundaries. We employed three unsupervised methods (BIN, ABGD, GMYC) with specified criteria on how to handle discordant results and examined diagnostic bases from each delineated putative species (operational taxonomic units, OTUs). Third, we evaluated the morphological characters of each OTU. Finally, we associated short barcodes from types with the delineated OTUs. In this step, we employed various supervised methods, including distance‐based, tree‐based and character‐based. We recovered 658 bp barcode sequences from 194 of 215 fresh specimens and recovered an average of 141 bp from 33 of 42 holotypes. We observed strong congruence among all methods and good correspondence with morphology. We demonstrate potential pitfalls with tree‐, distance‐ and character‐based approaches when associating sequences of varied length. Our results suggest that sequences as short as 56 bp can often provide valuable taxonomic information. The results support significant taxonomic oversplitting of species in the Elachista dispunctella complex.

Testing DNA barcode performance in 1000 species of European lepidoptera: large geographic distances have small genetic impacts

This study examines the performance of DNA barcodes (mt cytochrome c oxidase 1 gene) in the identification of 1004 species of Lepidoptera shared by two localities (Finland, Austria) that are 1600 km apart. Maximum intraspecific distances for the pooled data were less than 2% for 880 species (87.6%), while deeper divergence was detected in 124 species. Despite such variation, the overall DNA barcode library possessed diagnostic COI sequences for 98.8% of the taxa. Because a reference library based on Finnish specimens was highly effective in identifying specimens from Austria, we conclude that barcode libraries based on regional sampling can often be effective for a much larger area. Moreover, dispersal ability (poor, good) and distribution patterns (disjunct, fragmented, continuous, migratory) had little impact on levels of intraspecific geographic divergence. Furthermore, the present study revealed that, despite the intensity of past taxonomic work on European Lepidoptera, nearly 20% of the species shared by Austria and Finland require further work to clarify their status. Particularly discordant BIN (Barcode Index Number) cases should be checked to ascertain possible explanatory factors such as incorrect taxonomy, hybridization, introgression, and Wolbachia infections.

Effects of phylogenetic reconstruction method on the robustness of species delimitation using single-locus data

Coalescent-based species delimitation methods combine population genetic and phylogenetic theory to provide an objective means for delineating evolutionarily significant units of diversity. The generalised mixed Yule coalescent (GMYC) and the Poisson tree process (PTP) are methods that use ultrametric (GMYC or PTP) or non-ultrametric (PTP) gene trees as input, intended for use mostly with single-locus data such as DNA barcodes.

Here, we assess how robust the GMYC and PTP are to different phylogenetic reconstruction and branch smoothing methods. We reconstruct over 400 ultrametric trees using up to 30 different combinations of phylogenetic and smoothing methods and perform over 2000 separate species delimitation analyses across 16 empirical data sets. We then assess how variable diversity estimates are, in terms of richness and identity, with respect to species delimitation, phylogenetic and smoothing methods.

The PTP method generally generates diversity estimates that are more robust to different phylogenetic methods. The GMYC is more sensitive, but provides consistent estimates for BEAST trees. The lower consistency of GMYC estimates is likely a result of differences among gene trees introduced by the smoothing step. Unresolved nodes (real anomalies or methodological artefacts) affect both GMYC and PTP estimates, but have a greater effect on GMYC estimates. Branch smoothing is a difficult step and perhaps an underappreciated source of bias that may be widespread among studies of diversity and diversification.

Nevertheless, careful choice of phylogenetic method does produce equivalent PTP and GMYC diversity estimates. We recommend simultaneous use of the PTP model with any model-based gene tree (e.g. RAxML) and GMYC approaches with BEAST trees for obtaining species hypotheses.

'Direct PCR' optimization yields a rapid, cost-effective, nondestructive and efficient method for obtaining DNA barcodes without DNA extraction

Macroinvertebrates that are collected in large numbers pose major problems in basic and applied biodiversity research: identification to species via morphology is often difficult, slow and/or expensive. DNA barcodes are an attractive alternative or complementary source of information. Unfortunately, obtaining DNA barcodes from specimens requires many steps and thus time and money. Here, we promote a short cut to DNA barcoding, that is, a nondestructive PCR method that skips DNA extraction ('direct PCR') and that can be used for a broad range of invertebrate taxa. We demonstrate how direct PCR can be optimized for the larvae and adults of nonbiting midges (Diptera: Chironomidae), a typical invertebrate group that is abundant, contains important bioindicator species, but is difficult to identify based on morphological features. After optimization, direct PCR yields high PCR success rates (>90%), preserves delicate morphological features (e.g. details of genitalia, and larval head capsules) while allowing for the recovery of genomic DNA. We also document that direct PCR can be successfully optimized for a wide range of other invertebrate taxa that need routine barcoding (flies: Culicidae, Drosophilidae, Dolichopodidae, Sepsidae; sea stars: Oreasteridae). Key for obtaining high PCR success rates is optimizing (i) tissue quantity, (ii) body part, (iii) primer pair and (iv) type of Taq polymerase. Unfortunately, not all invertebrates appear suitable because direct PCR has low success rates for other taxa that were tested (e.g. Coleoptera: Dytiscidae, Copepoda, Hymenoptera: Formicidae and Odonata). It appears that the technique is less successful for heavily sclerotized insects and/or those with many exocrine glands.

The slugs of Britain and Ireland: undetected and undescribed species increase a well-studied, economically important fauna by more than 20%

The slugs of Britain and Ireland form a well-studied fauna of economic importance. They include many widespread European species that are introduced elsewhere (at least half of the 36 currently recorded British species are established in North America, for example). To test the contention that the British and Irish fauna consists of 36 species, and to verify the identity of each, a species delimitation study was conducted based on a geographically wide survey. Comparisons between mitochondrial DNA (COI, 16S), nuclear DNA (ITS-1) and morphology were investigated with reference to interspecific hybridisation. Species delimitation of the fauna produced a primary species hypothesis of 47 putative species. This was refined to a secondary species hypothesis of 44 species by integration with morphological and other data. Thirty six of these correspond to the known fauna (two species in Arion subgenus Carinarion were scarcely distinct and Arion (Mesarion) subfuscus consisted of two near-cryptic species). However, by the same criteria a further eight previously undetected species (22% of the fauna) are established in Britain and/or Ireland. Although overlooked, none are strictly morphologically cryptic, and some appear previously undescribed. Most of the additional species are probably accidentally introduced, and several are already widespread in Britain and Ireland (and thus perhaps elsewhere). At least three may be plant pests. Some evidence was found for interspecific hybridisation among the large Arion species (although not involving A. flagellus) and more unexpectedly in species pairs in Deroceras (Agriolimacidae) and Limacus (Limacidae). In the latter groups, introgression appears to have occurred in one direction only, with recently-invading lineages becoming common at the expense of long-established or native ones. The results show how even a well-studied, macroscopic fauna can be vulnerable to cryptic and undetected invasions and changes.

A transcontinental challenge--a test of DNA barcode performance for 1,541 species of Canadian Noctuoidea (Lepidoptera)

This study provides a first, comprehensive, diagnostic use of DNA barcodes for the Canadian fauna of noctuoids or "owlet" moths (Lepidoptera: Noctuoidea) based on vouchered records for 1,541 species (99.1% species coverage), and more than 30,000 sequences. When viewed from a Canada-wide perspective, DNA barcodes unambiguously discriminate 90% of the noctuoid species recognized through prior taxonomic study, and resolution reaches 95.6% when considered at a provincial scale. Barcode sharing is concentrated in certain lineages with 54% of the cases involving 1.8% of the genera. Deep intraspecific divergence exists in 7.7% of the species, but further studies are required to clarify whether these cases reflect an overlooked species complex or phylogeographic variation in a single species. Non-native species possess higher Nearest-Neighbour (NN) distances than native taxa, whereas generalist feeders have lower NN distances than those with more specialized feeding habits. We found high concordance between taxonomic names and sequence clusters delineated by the Barcode Index Number (BIN) system with 1,082 species (70%) assigned to a unique BIN. The cases of discordance involve both BIN mergers and BIN splits with 38 species falling into both categories, most likely reflecting bidirectional introgression. One fifth of the species are involved in a BIN merger reflecting the presence of 158 species sharing their barcode sequence with at least one other taxon, and 189 species with low, but diagnostic COI divergence. A very few cases (13) involved species whose members fell into both categories. Most of the remaining 140 species show a split into two or three BINs per species, while Virbia ferruginosa was divided into 16. The overall results confirm that DNA barcodes are effective for the identification of Canadian noctuoids. This study also affirms that BINs are a strong proxy for species, providing a pathway for a rapid, accurate estimation of animal diversity.

Phylogeny, systematics and biogeography of the genus panolis (lepidoptera: noctuidae) based on morphological and molecular evidence

The genus Panolis is a small group of noctuid moths with six recognized species distributed from Europe to East Asia, and best known for containing the widespread Palearctic pest species P. flammea, the pine beauty moth. However, a reliable classification and robust phylogenetic framework for this group of potentially economic importance are currently lacking. Here, we use morphological and molecular data (mitochondrial genes cytochrome c oxidase subunit I and 16S ribosomal RNA, nuclear gene elongation factor-1 alpha) to reconstruct the phylogeny of this genus, with a comprehensive systematic revision of all recognized species and a new one, P. ningshan sp. nov. The analysis results of maximum parsimony, maximum likelihood and Bayesian inferring methods for the combined morphological and molecular data sets are highly congruent, resulting in a robust phylogeny and identification of two clear species groups, i.e., the P. flammea species group and the P. exquisita species group. We also estimate the divergence times of Panolis moths using two conventional mutation rates for the arthropod mitochondrial COI gene with a comparison of two molecular clock models, as well as reconstruct their ancestral areas. Our results suggest that 1) Panolis is a young clade, originating from the Oriental region in China in the Late Miocene (6-10Mya), with an ancestral species in the P. flammea group extending northward to the Palearctic region some 3-6 Mya; 2) there is a clear possibility for a representative of the Palearctic clade to become established as an invasive species in the Nearctic taiga.

Unraveling the sequence information in COI barcode to achieve higher taxon assignment based on Indian freshwater fishes

Efficacy of cytochrome c oxidase subunit I (COI) DNA barcode in higher taxon assignment is still under debate in spite of several attempts, using the conventional DNA barcoding methods, to assign higher taxa. Here we try to understand whether nucleotide and amino acid sequence in COI gene carry sufficient information to assign species to their higher taxonomic rank, using 160 species of Indian freshwater fishes. Our results reveal that with increase in the taxonomic rank, sequence conservation decreases for both nucleotides and amino acids. Order level exhibits lowest conservation with 50% of the nucleotides and amino acids being conserved. Among the variable sites, 30-50% were found to carry high information content within an order, while it was 70-80% within a family and 80-99% within a genus. High information content shows sites with almost conserved sequence but varying at one or two locations, which can be due to variations at species or population level. Thus, the potential of COI gene in higher taxon assignment is revealed with validation of ample inherent signals latent in the gene.

DNA barcodes identify Central Asian Colias butterflies (Lepidoptera, Pieridae)

A majority of the known Colias species (Lepidoptera: Pieridae, Coliadinae) occur in the mountainous regions of Central-Asia, vast areas that are hard to access, rendering the knowledge of many species limited due to the lack of extensive sampling. Two gene regions, the mitochondrial COI ‘barcode’ region and the nuclear ribosomal protein RpS2 gene region were used for exploring the utility of these DNA markers for species identification. A comprehensive sampling of COI barcodes for Central Asian Colias butterflies showed that the barcodes facilitated identification of most of the included species. Phylogenetic reconstruction based on parsimony and Neighbour-Joining recovered most species as monophyletic entities. For the RpS2 gene region species-specific sequences were registered for some of the included Colias spp. Nevertheless, this gene region was not deemed useful as additional molecular ‘barcode’. A parsimony analysis of the combined COI and RpS2 data did not support the current subgeneric classification based on morphological characteristics.

Genetic patterns in European geometrid moths revealed by the Barcode Index Number (BIN) system

Background

The geometrid moths of Europe are one of the best investigated insect groups in traditional taxonomy making them an ideal model group to test the accuracy of the Barcode Index Number (BIN) system of BOLD (Barcode of Life Datasystems), a method that supports automated, rapid species delineation and identification.

Methodology/Principal Findings

This study provides a DNA barcode library for 219 of the 249 European geometrid moth species (88%) in five selected subfamilies. The data set includes COI sequences for 2130 specimens. Most species (93%) were found to possess diagnostic barcode sequences at the European level while only three species pairs (3%) were genetically indistinguishable in areas of sympatry. As a consequence, 97% of the European species we examined were unequivocally discriminated by barcodes within their natural areas of distribution. We found a 1:1 correspondence between BINs and traditionally recognized species for 67% of these species. Another 17% of the species (15 pairs, three triads) shared BINs, while specimens from the remaining species (18%) were divided among two or more BINs. Five of these species are mixtures, both sharing and splitting BINs. For 82% of the species with two or more BINs, the genetic splits involved allopatric populations, many of which have previously been hypothesized to represent distinct species or subspecies.

Conclusions/Significance

This study confirms the effectiveness of DNA barcoding as a tool for species identification and illustrates the potential of the BIN system to characterize formal genetic units independently of an existing classification. This suggests the system can be used to efficiently assess the biodiversity of large, poorly known assemblages of organisms. For the moths examined in this study, cases of discordance between traditionally recognized species and BINs arose from several causes including overlooked species, synonymy, and cases where DNA barcodes revealed regional variation of uncertain taxonomic significance.

Building a DNA barcode reference library for the true butterflies (Lepidoptera) of Peninsula Malaysia: what about the subspecies?

The objective of this study was to build a DNA barcode reference library for the true butterflies of Peninsula Malaysia and assess the value of attaching subspecies names to DNA barcode records. A new DNA barcode library was constructed with butterflies from the Museum of Zoology, University of Malaya collection. The library was analysed in conjunction with publicly available DNA barcodes from other Asia-Pacific localities to test the ability of the DNA barcodes to discriminate species and subspecies. Analyses confirmed the capacity of the new DNA barcode reference library to distinguish the vast majority of species (92%) and revealed that most subspecies possessed unique DNA barcodes (84%). In some cases conspecific subspecies exhibited genetic distances between their DNA barcodes that are typically seen between species, and these were often taxa that have previously been regarded as full species. Subspecies designations as shorthand for geographically and morphologically differentiated groups provide a useful heuristic for assessing how such groups correlate with clustering patterns of DNA barcodes, especially as the number of DNA barcodes per species in reference libraries increases. Our study demonstrates the value in attaching subspecies names to DNA barcode records as they can reveal a history of taxonomic concepts and expose important units of biodiversity.

A retrospective approach to testing the DNA barcoding method

A decade ago, DNA barcoding was proposed as a standardised method for identifying existing species and speeding the discovery of new species. Yet, despite its numerous successes across a range of taxa, its frequent failures have brought into question its accuracy as a short-cut taxonomic method. We use a retrospective approach, applying the method to the classification of New Zealand skinks as it stood in 1977 (primarily based upon morphological characters), and compare it to the current taxonomy reached using both morphological and molecular approaches. For the 1977 dataset, DNA barcoding had moderate-high success in identifying specimens (78-98%), and correctly flagging specimens that have since been confirmed as distinct taxa (77-100%). But most matching methods failed to detect the species complexes that were present in 1977. For the current dataset, there was moderate-high success in identifying specimens (53-99%). For both datasets, the capacity to discover new species was dependent on the methodological approach used. Species delimitation in New Zealand skinks was hindered by the absence of either a local or global barcoding gap, a result of recent speciation events and hybridisation. Whilst DNA barcoding is potentially useful for specimen identification and species discovery in New Zealand skinks, its error rate could hinder the progress of documenting biodiversity in this group. We suggest that integrated taxonomic approaches are more effective at discovering and describing biodiversity.

Wide-ranging barcoding aids discovery of one-third increase of species richness in presumably well-investigated moths

Rapid development of broad regional and international DNA barcode libraries have brought new insights into the species diversity of many areas and groups. Many new species, even within well-investigated species groups, have been discovered based initially on differences in DNA barcodes. We barcoded 437 collection specimens belonging to 40 pre-identified Palearctic species of the Elachista bifasciella group of moths (Lepidoptera, Elachistidae). Although the study group has been a subject of several careful morphological taxonomic examinations, an unexpectedly high number of previously undetected putative species is revealed, resulting in a 34% rise in species number in the study area. The validity of putative new species was subsequently supported with diagnostic morphological traits. We show that DNA barcodes provide a powerful method of detecting potential new species even in taxonomic groups and geographic areas that have previously been under considerable morphological taxonomic scrutiny.

In the shadow of phylogenetic uncertainty: the recent diversification of Lysandra butterflies through chromosomal change

The phylogeny of the butterfly genus Lysandra (Lycaenidae, Polyommatinae) has been intractable using both molecular and morphological characters, which could be a result of speciation due to karyotype instability. Here we reconstruct the phylogeny of the group using multi-locus coalescent-based methods on seven independent genetic markers. While the genus is ca. 4.9 Mya old, the diversification of the extant lineages was extremely recent (ca. 1.5 Mya) and involved multiple chromosomal rearrangements. We find that relationships are uncertain due to both incomplete lineage sorting and hybridization. Minimizing the impact of reticulation in inferring the species tree by testing for mitochondrial introgression events yields a partially resolved tree with three main supported clades: L. punctifera+L. bellargus, the corydonius taxa, and L. coridon+the Iberian taxa, plus three independent lineages without apparently close relatives (L. ossmar, L. syriaca and L. dezina). Based on these results and new karyotypic data, we propose a rearrangement recognizing ten species within the genus. Finally, we hypothesize that chromosomal instability may have played a crucial role in the Lysandra recent diversification. New chromosome rearrangements might be fixed in populations after severe bottlenecks, which at the same time might promote rapid sorting of neutral molecular markers. We argue that population bottlenecks might be a prerequisite for chromosomal speciation in this group.

Separation possibilities and genital measurement variations in two cryptic species of European pierid butterflies, Leptidea juvernica Williams, 1946 and L. sinapis (Linnaeus, 1758)

Of the three recently separated cryptic butterfly species of the Leptidea sinapis complex, the two species L. sinapis and L. juvernica occur sympatrically and syntopically in central Europe. As the separation of these species requires genital or genetic characters, their correct identification is a fundamental problem to be solved prior to any biological and distributional studies. In the present study a morphometric approach was applied to test for separation possibilities and to examine genital measurement variations based on large population samples of the two species (347 females and 636 males). Butterflies were collected at 456 localities distributed across Poland. Specimens of both sexes could be separated using either the shape or the length of the antrum bursae (females), phallus and saccus (males). Intraspecific seasonal differences in these measurements may affect the discrimination of males only. Genital characters were significantly larger in butterflies of the spring brood than in those of the second brood. In females, antrum bursae length ranges were separated by an interval, allowing for the exact identification of each specimen. In males, no such intervals were found between phallus and saccus length ranges, which slightly overlapped. Discriminant analysis resulted in 100% of males correctly classified. The present study confirmed the validity of interspecific differences in the shape of phallus and saccus for species identification without measurements. The length of these two characters discriminates the two species when separation is limited to specimens from the same generations. As an alternative for discriminant analysis, either phallus length and vinculum width or phallus length and the ratio of the ventral edge of genital capsule length and saccus length (VEL/SL ratio) is recommended for the practical separation of the two species. The differences in the shape and size of the copulatory organs of the two species seem to indicate reproductive isolation due to mechanical incompatibility of their genitals.

When the rule becomes the exception. no evidence of gene flow between two Zerynthia cryptic butterflies suggests the emergence of a new model group

There is increasing evidence that most parapatric cryptic/sister taxa are reproductively compatible across their areas of contact. Consequently, the biological species concept, which assumes absence of interbreeding, is becoming a not so effective criterion in evolutionary ecology. Nevertheless, the few parapatric sister taxa showing complete reproductive barriers represent interesting models to study speciation processes and the evolution of reproductive isolation. In this study, we examined contact populations in northwestern Italy of two butterfly species, Zerynthia polyxena and Z. cassandra, characterized by different genitalic morphotypes. We studied levels of divergence among 21 populations distributed from Sicily to France using three genetic markers (the mitochondrial COI and ND1 genes and the nuclear wingless gene) and genitalic geometric morphometrics. Moreover, we performed species distribution modelling to estimate different climatic requirements of Z. polyxena and Z. cassandra. We projected climatic data into glacial maximum scenarios in order to verify if and to which extent glacial cycles could have contributed to speciation processes. Genetic and morphometric analyses identified two main groups. All specimens showed a concordant pattern of diversification, including those individuals sampled in the contact area. Haplotype distribution and climatic models showed that during glacial maxima both species experienced a strong range contraction and presumably remained separated into different microrefugia in southern France, in the Italian Peninsula and on the islands of Elba and Sicily. Long term separation was probably favoured by reduced dispersal ability and high phylopatry, while genitalic diversification probably favoured interbreeding avoidance. Conversely, the aposematic wing pattern remained almost identical. We compared our results with those obtained in other species and concluded that Z. polyxena and Z. cassandra represent a valuable model in the study of speciation.

Quantifying species diversity with a DNA barcoding-based method: Tibetan moth species (Noctuidae) on the Qinghai-Tibetan Plateau

With the ongoing loss of biodiversity, there is a great need for fast and effective ways to assess species richness and diversity: DNA barcoding provides a powerful new tool for this. We investigated this approach by focusing on the Tibetan plateau, which is one of the world's top biodiversity hotspots. There have been few studies of its invertebrates, although they constitute the vast majority of the region's diversity. Here we investigated species diversity of the lepidopteran family Noctuidae, across different environmental gradients, using measurements based on traditional morphology as well as on DNA barcoding. The COI barcode showed an average interspecific K2P distance of 9.45±2.08%, which is about four times larger than the mean intraspecific distance (1.85±3.20%). Using six diversity indices, we did not detect any significant differences in estimated species diversity between measurements based on traditional morphology and on DNA barcoding. Furthermore, we found strong positive correlations between them, indicating that barcode-based measures of species diversity can serve as a good surrogate for morphology-based measures in most situations tested. Eastern communities were found to have significantly higher diversity than Western ones. Among 22 environmental factors tested, we found that three (precipitation of driest month, precipitation of driest quarter, and precipitation of coldest quarter) were significantly correlated with species diversity. Our results indicate that these factors could be the key ecological factors influencing the species diversity of the lepidopteran family Noctuidae on the Tibetan plateau.

A simple 2D non-parametric resampling statistical approach to assess confidence in species identification in DNA barcoding--an alternative to likelihood and bayesian approaches

In the recent worldwide campaign for the global biodiversity inventory via DNA barcoding, a simple and easily used measure of confidence for assigning sequences to species in DNA barcoding has not been established so far, although the likelihood ratio test and the bayesian approach had been proposed to address this issue from a statistical point of view. The TDR (Two Dimensional non-parametric Resampling) measure newly proposed in this study offers users a simple and easy approach to evaluate the confidence of species membership in DNA barcoding projects. We assessed the validity and robustness of the TDR approach using datasets simulated under coalescent models, and an empirical dataset, and found that TDR measure is very robust in assessing species membership of DNA barcoding. In contrast to the likelihood ratio test and bayesian approach, the TDR method stands out due to simplicity in both concepts and calculations, with little in the way of restrictive population genetic assumptions. To implement this approach we have developed a computer program package (TDR1.0beta) freely available from ftp://202.204.209.200/education/video/TDR1.0beta.rar.

Biogeography and systematics of Aricia butterflies (Lepidoptera, Lycaenidae)

Butterflies of the Aricia species group represent a paradigm of unresolved taxonomy, both at the genus and species levels. We studied phylogenetic relationships, biogeography, and systematics based on genetic--nuclear and mitochondrial--and morphometric--external (wings) and internal (genitalia) - data. We show that Aricia is a monophyletic genus comprising the taxa Pseudoaricia, Ultraaricia and Umpria, which are here considered junior synonyms of Aricia. The taxa allous, inhonora, issekutzi, mandzhuriana, myrmecias and transalaica, which have often been raised to species rank, are shown to probably represent subspecies or synonyms. We show that montensis is likely a good species that is sister to all A. artaxerxes populations across the Palearctic region. The species A. anteros and A. morronensis are shown to display deep intraspecific divergences and they may harbor cryptic species. We also discovered that A. cramera and A. agestis exhibit a pattern of mutual exclusion on islands, and a parapatric distribution in mainland with a narrow contact zone where potential hybrids were detected. The lack of a prezygotic barrier that prevents their coexistence could explain this phenomenon. This study will hopefully contribute to the stability of the systematics of Aricia, a group with potential for the study of the link between speciation and biogeography.

Allopatry as a gordian knot for taxonomists: patterns of DNA barcode divergence in arctic-alpine lepidoptera

Many cold adapted species occur in both montane settings and in the subarctic. Their disjunct distributions create taxonomic complexity because there is no standardized method to establish whether their allopatric populations represent single or different species. This study employs DNA barcoding to gain new perspectives on the levels and patterns of sequence divergence among populations of 122 arctic-alpine species of Lepidoptera from the Alps, Fennoscandia and North America. It reveals intraspecific variability in the barcode region ranging from 0.00-10.08%. Eleven supposedly different species pairs or groups show close genetic similarity, suggesting possible synonymy in many cases. However, a total of 33 species show evidence of cryptic diversity as evidenced by the presence of lineages with over 2% maximum barcode divergence in Europe, in North America or between the two continents. Our study also reveals cases where taxonomic names have been used inconsistently between regions and exposes misidentifications. Overall, DNA barcodes have great potential to both increase taxonomic resolution and to make decisions concerning the taxonomic status of allopatric populations more objective.

A trans-Amazonian screening of mtDNA reveals deep intraspecific divergence in forest birds and suggests a vast underestimation of species diversity

The Amazonian avifauna remains severely understudied relative to that of the temperate zone, and its species richness is thought to be underestimated by current taxonomy. Recent molecular systematic studies using mtDNA sequence reveal that traditionally accepted species-level taxa often conceal genetically divergent subspecific lineages found to represent new species upon close taxonomic scrutiny, suggesting that intraspecific mtDNA variation could be useful in species discovery. Surveys of mtDNA variation in Holarctic species have revealed patterns of variation that are largely congruent with species boundaries. However, little information exists on intraspecific divergence in most Amazonian species. Here we screen intraspecific mtDNA genetic variation in 41 Amazonian forest understory species belonging to 36 genera and 17 families in 6 orders, using 758 individual samples from Ecuador and French Guiana. For 13 of these species, we also analyzed trans-Andean populations from the Ecuadorian Chocó. A consistent pattern of deep intraspecific divergence among trans-Amazonian haplogroups was found for 33 of the 41 taxa, and genetic differentiation and genetic diversity among them was highly variable, suggesting a complex range of evolutionary histories. Mean sequence divergence within families was the same as that found in North American birds (13%), yet mean intraspecific divergence in Neotropical species was an order of magnitude larger (2.13% vs. 0.23%), with mean distance between intraspecific lineages reaching 3.56%. We found no clear relationship between genetic distances and differentiation in plumage color. Our results identify numerous genetically and phenotypically divergent lineages which may result in new species-level designations upon closer taxonomic scrutiny and thorough sampling, although lineages in the tropical region could be older than those in the temperate zone without necessarily representing separate species. In-depth phylogeographic surveys are urgently needed to avoid underestimating tropical diversity, and the use of mtDNA markers can be instrumental in identifying and prioritizing taxa for species discovery.

Haplotype and minimum-chimerism consensus determination using short sequence data

Background

Assembling haplotypes given sequence data derived from a single individual is a well studied problem, but only recently has haplotype assembly been considered for population-sampled data. We discuss a software tool called Hapler, which is designed specifically for low-diversity, low-coverage data such as ecological samples derived from natural populations. Because such data may contain error as well as ambiguous haplotype information, we developed methods that increase confidence in these assemblies. Hapler also reconstructs full consensus sequences while minimizing and identifying possible chimeric points.

Results

Experiments on simulated data indicate that Hapler is effective at assembling haplotypes from gene-sized alignments of short reads. Further, in our tests Hapler-generated consensus sequences are less chimeric than the alternative consensus approaches of majority vote and viral quasispecies estimation regardless of error rate, read length, or population haplotype bias.

Conclusions

The analysis of genetically diverse sequence data is increasingly common, particularly in the field of ecoinformatics where transcriptome sequencing of natural populations is a cost effective alternative to genome sequencing. For such studies, it is important to consider and identify haplotype diversity. Hapler provides robust haplotype information and identifies possible phasing errors in consensus sequences, providing valuable information for population studies and downstream usage of resulting assemblies.

Barcoding and border biosecurity: identifying cyprinid fishes in the aquarium trade

BACKGROUND

Poorly regulated international trade in ornamental fishes poses risks to both biodiversity and economic activity via invasive alien species and exotic pathogens. Border security officials need robust tools to confirm identifications, often requiring hard-to-obtain taxonomic literature and expertise. DNA barcoding offers a potentially attractive tool for quarantine inspection, but has yet to be scrutinised for aquarium fishes. Here, we present a barcoding approach for ornamental cyprinid fishes by: (1) expanding current barcode reference libraries; (2) assessing barcode congruence with morphological identifications under numerous scenarios (e.g. inclusion of GenBank data, presence of singleton species, choice of analytical method); and (3) providing supplementary information to identify difficult species.

METHODOLOGY/PRINCIPAL FINDINGS

We sampled 172 ornamental cyprinid fish species from the international trade, and provide data for 91 species currently unrepresented in reference libraries (GenBank/Bold). DNA barcodes were found to be highly congruent with our morphological assignments, achieving success rates of 90-99%, depending on the method used (neighbour-joining monophyly, bootstrap, nearest neighbour, GMYC, percent threshold). Inclusion of data from GenBank (additional 157 spp.) resulted in a more comprehensive library, but at a cost to success rate due to the increased number of singleton species. In addition to DNA barcodes, our study also provides supporting data in the form of specimen images, morphological characters, taxonomic bibliography, preserved vouchers, and nuclear rhodopsin sequences. Using this nuclear rhodopsin data we also uncovered evidence of interspecific hybridisation, and highlighted unrecognised diversity within popular aquarium species, including the endangered Indian barb Puntius denisonii.

CONCLUSIONS/SIGNIFICANCE

We demonstrate that DNA barcoding provides a highly effective biosecurity tool for rapidly identifying ornamental fishes. In cases where DNA barcodes are unable to offer an identification, we improve on previous studies by consolidating supplementary information from multiple data sources, and empower biosecurity agencies to confidently identify high-risk fishes in the aquarium trade.

DNA barcoding of recently diverged species: relative performance of matching methods

Recently diverged species are challenging for identification, yet they are frequently of special interest scientifically as well as from a regulatory perspective. DNA barcoding has proven instrumental in species identification, especially in insects and vertebrates, but for the identification of recently diverged species it has been reported to be problematic in some cases. Problems are mostly due to incomplete lineage sorting or simply lack of a 'barcode gap' and probably related to large effective population size and/or low mutation rate. Our objective was to compare six methods in their ability to correctly identify recently diverged species with DNA barcodes: neighbor joining and parsimony (both tree-based), nearest neighbor and BLAST (similarity-based), and the diagnostic methods DNA-BAR, and BLOG. We analyzed simulated data assuming three different effective population sizes as well as three selected empirical data sets from published studies. Results show, as expected, that success rates are significantly lower for recently diverged species (∼75%) than for older species (∼97%) (P<0.00001). Similarity-based and diagnostic methods significantly outperform tree-based methods, when applied to simulated DNA barcode data (P<0.00001). The diagnostic method BLOG had highest correct query identification rate based on simulated (86.2%) as well as empirical data (93.1%), indicating that it is a consistently better method overall. Another advantage of BLOG is that it offers species-level information that can be used outside the realm of DNA barcoding, for instance in species description or molecular detection assays. Even though we can confirm that identification success based on DNA barcoding is generally high in our data, recently diverged species remain difficult to identify. Nevertheless, our results contribute to improved solutions for their accurate identification.

DNA barcodes for insects

DNA barcoding refers to the technique of sequencing a short fragment of the mitochondrial cytochrome c oxidase subunit I (COI) gene, the "DNA barcode," from a taxonomically unknown specimen and performing comparisons with a reference library of barcodes of known species origin to establish a species-level identification. The library barcodes gain their value due to an intimate association-through the vouchered specimens from where they came-with other data; particularly Linnaean names, collection localities, and morphology in the form of digital images. Consequently, this chapter details means of efficiently obtaining barcodes along two general streams: rapid barcode assembly to populate the library and retrieval of barcodes from highly prized specimens, but also emphasizes organization and collection of the barcode collaterals.

Future directions

It is a risky task to attempt to predict the direction that DNA barcoding and its applications may take in the future. In a very short time, the endeavor of DNA barcoding has gone from being a tool to facilitate taxonomy in difficult to identify species, to an ambitious, global initiative that seeks to tackle such pertinent and challenging issues as quantifying global biodiversity, revolutionizing the forensic identifications of species, advancing the study of interactions among species, and promoting the reconstruction of evolutionary relationships within communities. The core element of DNA barcoding will always remain the same: the generation of a set of well-identified samples collected and genotyped at one or more genetic barcode markers and assembled into a properly curated database. But the application of this body of data will depend on the creativity and need of the research community in using a "gold standard" of annotated DNA sequence data at the species level. We foresee several areas where the application of DNA barcode data is likely to yield important evolutionary, ecological, and societal insights, and while far from exclusive, provide examples of how DNA barcode data will continue to empower scientists to address hypothesis-driven research. Three areas of immediate and obvious concern are (1) biodiversity inventories, (2) phylogenetic applications, and (3) species interactions.

Host tracking or cryptic adaptation? Phylogeography of Pediobius saulius (Hymenoptera, Eulophidae), a parasitoid of the highly invasive horse-chestnut leafminer

Classical biological control is often advocated as a tool for managing invasive species. However, accurate evaluations of parasitoid species complexes and assessment of host specificity are impeded by the lack of morphological variation. Here, we study the possibility of host races/species within the eulophid wasp Pediobius saulius, a pupal generalist parasitoid that parasitize the highly invasive horse-chestnut leaf-mining moth Cameraria ohridella. We analysed the population genetic structure, host associations and phylogeographic patterns of P. saulius in Europe using the COI mitochondrial gene. This marker strongly supports a division into at least five highly differentiated parasitoid complexes, within two of which clades with differing degrees of host specialization were found: a Balkan clade that mainly (but not only) attacks C. ohridella and a more generalist European group that attacks many hosts, including C. ohridella. The divergence in COI (up to 7.6%) suggests the existence of cryptic species, although this is neither confirmed by nuclear divergence nor morphology. We do not find evidence of host tracking. The higher parasitism rates observed in the Balkans and the scarcity of the Balkan–Cameraria haplotypes out of the Balkans open the possibility of using these Balkan haplotypes as biological control agents of C. ohridella elsewhere in Europe.

A fuzzy-set-theory-based approach to analyse species membership in DNA barcoding

Reliable assignment of an unknown query sequence to its correct species remains a methodological problem for the growing field of DNA barcoding. While great advances have been achieved recently, species identification from barcodes can still be unreliable if the relevant biodiversity has been insufficiently sampled. We here propose a new notion of species membership for DNA barcoding-fuzzy membership, based on fuzzy set theory-and illustrate its successful application to four real data sets (bats, fishes, butterflies and flies) with more than 5000 random simulations. Two of the data sets comprise especially dense species/population-level samples. In comparison with current DNA barcoding methods, the newly proposed minimum distance (MD) plus fuzzy set approach, and another computationally simple method, 'best close match', outperform two computationally sophisticated Bayesian and BootstrapNJ methods. The new method proposed here has great power in reducing false-positive species identification compared with other methods when conspecifics of the query are absent from the reference database.

A combined genetic-morphometric analysis unravels the complex biogeographical history of Polyommatus icarus and Polyommatus celina common blue butterflies

Widespread species have the potential to reveal large-scale biogeographical patterns, as well as responses to environmental changes possibly unique to habitat generalists. This study presents a continental-scale phylogeographical analysis of Polyommatus icarus, one of the most common Palaearctic butterflies, and the morphologically and ecologically similar Polyommatus celina, a recently discovered cryptic species. By combining data from mitochondrial [cytochrome c oxidase subunit I (COI)] and nuclear [internal transcribed spacer (ITS2)] molecular markers with geometric morphometrics, we document a complex phylogeographical history for the two species. Despite morphological similarities, the genetic divergence between these two species is high (more than 5% at COI) and they are not sister species. For the first time, we show that P. celina occurs not only in North Africa but also in Europe, where it inhabits several west Mediterranean islands, as well as large parts of Iberia, where it occurs in parapatry with P. icarus. The two species appear to completely exclude each other on islands, but we provide morphological and molecular evidence that introgression occurred in the Iberian Peninsula. We discovered strongly diverged lineages that seem to represent relict populations produced by past range expansions and contractions: Crete and Iberian isolates for P. icarus, Balearics-Sardinia and Sicily-Lipari for P. celina. This study shows that a combined genetic-morphometric approach can shed light on cryptic diversity while providing the necessary resolution to reconstruct a fine-scale phylogeographical history of species at both spatial and temporal levels.

Reading the complex skipper butterfly fauna of one tropical place

Background

An intense, 30-year, ongoing biodiversity inventory of Lepidoptera, together with their food plants and parasitoids, is centered on the rearing of wild-caught caterpillars in the 120,000 terrestrial hectares of dry, rain, and cloud forest of Area de Conservacion Guanacaste (ACG) in northwestern Costa Rica. Since 2003, DNA barcoding of all species has aided their identification and discovery. We summarize the process and results for a large set of the species of two speciose subfamilies of ACG skipper butterflies (Hesperiidae) and emphasize the effectiveness of barcoding these species (which are often difficult and time-consuming to identify).

Methodology/Principal Findings

Adults are DNA barcoded by the Biodiversity Institute of Ontario, Guelph, Canada; and they are identified by correlating the resulting COI barcode information with more traditional information such as food plant, facies, genitalia, microlocation within ACG, caterpillar traits, etc. This process has found about 303 morphologically defined species of eudamine and pyrgine Hesperiidae breeding in ACG (about 25% of the ACG butterfly fauna) and another 44 units indicated by distinct barcodes (n = 9,094), which may be additional species and therefore may represent as much as a 13% increase. All but the members of one complex can be identified by their DNA barcodes.

Conclusions/Significance

Addition of DNA barcoding to the methodology greatly improved the inventory, both through faster (hence cheaper) accurate identification of the species that are distinguishable without barcoding, as well as those that require it, and through the revelation of species “hidden” within what have long been viewed as single species. Barcoding increased the recognition of species-level specialization. It would be no more appropriate to ignore barcode data in a species inventory than it would be to ignore adult genitalia variation or caterpillar ecology.

DNA barcoding of lichenized fungi demonstrates high identification success in a floristic context

• Efforts are currently underway to establish a standard DNA barcode region for fungi; we tested the utility of the internal transcribed spacer (ITS) of nuclear ribosomal DNA for DNA barcoding in lichen-forming fungi by sampling diverse species across eight orders. • Amplification of the ITS region (ITS1-5.8S-ITS2) was conducted for 351 samples, encompassing 107, 55 and 28 species, genera and families, respectively, of lichenized fungi. We assessed the ability of the entire ITS vs the ITS2 alone to discriminate between species in a taxonomic dataset (members of the genus Usnea) and a floristic dataset. • In the floristic dataset, 96.3% of sequenced samples could be assigned to the correct species using ITS or ITS2; a barcode gap for ITS is present in 92.1% of species. Although fewer species have a barcode gap in the taxonomic dataset (73.3% with ITS and 68.8% with ITS2), up to 94.1% of samples were assigned to the correct species using BLAST. • While discrimination between the most closely related species will remain challenging, our results demonstrate the potential to identify a high percentage of specimens to the correct species, and the remainder to the correct genus, when using DNA barcoding in a floristic context.

A phylogenetic revision of the Glaucopsyche section (Lepidoptera: Lycaenidae), with special focus on the Phengaris-Maculinea clade

Despite much research on the socially parasitic large blue butterflies (genus Maculinea) in the past 40 years, their relationship to their closest relatives, Phengaris, is controversial and the relationships among the remaining genera in the Glaucopsyche section are largely unresolved. The evolutionary history of this butterfly section is particularly important to understand the evolution of life history diversity connected to food-plant and host-ant associations in the larval stage. In the present study, we use a combination of four nuclear and two mitochondrial genes to reconstruct the phylogeny of the Glaucopsyche section, and in particular, to study the relationships among and within the Phengaris-Maculinea species. We find a clear pattern between the clades recovered in the Glaucopsyche section phylogeny and their food-plant associations, with only the Phengaris-Maculinea clade utilising more than one plant family. Maculinea is, for the first time, recovered with strong support as a monophyletic group nested within Phengaris, with the closest relative being the rare genus Caerulea. The genus Glaucopsyche is polyphyletic, including the genera Sinia and Iolana. Interestingly, we find evidence for additional potential cryptic species within the highly endangered Maculinea, which has long been suspected from morphological, ecological and molecular studies.

Unprecedented within-species chromosome number cline in the Wood White butterfly Leptidea sinapis and its significance for karyotype evolution and speciation

Background

Species generally have a fixed number of chromosomes in the cell nuclei while between-species differences are common and often pronounced. These differences could have evolved through multiple speciation events, each involving the fixation of a single chromosomal rearrangement. Alternatively, marked changes in the karyotype may be the consequence of within-species accumulation of multiple chromosomal fissions/fusions, resulting in highly polymorphic systems with the subsequent extinction of intermediate karyomorphs. Although this mechanism of chromosome number evolution is possible in theory, it has not been well documented.

Results

We present the discovery of exceptional intraspecific variability in the karyotype of the widespread Eurasian butterfly Leptidea sinapis. We show that within this species the diploid chromosome number gradually decreases from 2n = 106 in Spain to 2n = 56 in eastern Kazakhstan, resulting in a 6000 km-wide cline that originated recently (8,500 to 31,000 years ago). Remarkably, intrapopulational chromosome number polymorphism exists, the chromosome number range overlaps between some populations separated by hundreds of kilometers, and chromosomal heterozygotes are abundant. We demonstrate that this karyotypic variability is intraspecific because in L. sinapis a broad geographical distribution is coupled with a homogenous morphological and genetic structure.

Conclusions

The discovered system represents the first clearly documented case of explosive chromosome number evolution through intraspecific and intrapopulation accumulation of multiple chromosomal changes. Leptidea sinapis may be used as a model system for studying speciation by means of chromosomally-based suppressed recombination mechanisms, as well as clinal speciation, a process that is theoretically possible but difficult to document. The discovered cline seems to represent a narrow time-window of the very first steps of species formation linked to multiple chromosomal changes that have occurred explosively. This case offers a rare opportunity to study this process before drift, dispersal, selection, extinction and speciation erase the traces of microevolutionary events and just leave the final picture of a pronounced interspecific chromosomal difference.