Reproductive isolation and patterns of genetic differentiation in a cryptic butterfly species complex

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.

Integrative analysis reveals cryptic speciation linked to habitat differentiation within Albanian populations of the anomalous blues (Lepidoptera, Lycaenidae, Polyommatus Latreille, 1804)

The Balkan Peninsula is one of the greatest hotspots for biodiversity in Europe. While the region has been investigated thoroughly, some parts remain understudied and may still harbour undiscovered diversity, even in well-studied organisms such as Lepidoptera. Here we investigated the group of the so-called anomalous blue butterflies, also known as 'brown complex' of the subgenus Agrodiaetus Hübner, 1822 including the taxa of the entire Polyommatusaroaniensis (Brown, 1976) species complex. This species complex is distributed in the southern part of the Balkan Peninsula and known to be represented by three closely related allopatric species, differentiated by their chromosome numbers (n) and mitochondrial (mt) DNA. These are P.aroaniensis sensu stricto (Southern Greece, Peloponnese, n=47-48; mt haplogroup aroa1), P.timfristos Lukhtanov, Vishnevskaya et Shapoval, 2016 (Central Greece, Attika, n=38, aroa2) and P.orphicus Kolev, 2005 (North-Eastern Greece, Southern Bulgaria, n=41-42, orph1). Based on an analysis of chromosomal, molecular and morphological markers, we demonstrate that a fourth taxon of this species complex exists in Albania. This taxon possesses the mt haplogroup aroa3, which is the most differentiated within the entire P.aroaniensis species complex, and the karyotype (n=42-43), which differs by one fixed chromosome fission from P.orphicus. The Albanian taxon seems to be ecologically specialised (habitat on dark-coloured, ophiolitic substrate soils) and differs in colouration (wing reflectance) from the others taxa of the P.aroaniensis species group. Based on the evidence here presented and following the current view of the taxonomy of the group, we propose considering the Albanian taxon as a new species, here described as Polyommatuslurae sp. nov. At the contact zone between the new species and P.orphicus, in addition to typical ones, we detected specimens with haplogroup orph2, karyotype n=43 and intermediate morphology, which seem to represent P.lurae × P.orphicus hybrids.

Phylogeny and Integrative Taxonomy of the Genera Gymnaetoides and Pseudotachycines (Orthoptera: Rhaphidophoridae)

The genera Gymnaetoides and Pseudotachycines are endemic to China and are morphologically homogeneous. The few available diagnostic characters make species identification particularly challenging. Species cannot be classified according to the given generic diagnosis, and phylogenetic analyses have not been reported. Here, we reconstruct the phylogeny using Bayesian inference and maximum likelihood and employ four approaches to delimit species. The results suggest that both Gymnaetoides and Pseudotachycines are paraphyletic. Therefore, we revise their taxonomy based on the combination of morphological characters and molecular data. A new genus Homotachycines Zhu & Shi gen. nov. is erected, and six new combinations are proposed. Species delimitation identifies 15 new species and one new subspecies: Gymnaetoides huangshanensis, G. petalus, G. yangmingensis, G. lushanensis, Pseudotachycines procerus, P. procerus guizhouensis, P. zhengi, P. nephrus, P. sagittus, P. fengyangshanensis, Homotachycines triangulus, H. quadratus, H. baokangensis, H. fusus, H. concavus, and H. qinlingensis sp. nov. Moreover, we find that the shapes of the dorsal lateral lobes and the dorsal median lobe of the male genitalia are also important characters for identifying these genera and that the shapes of the dorsal and lateral sclerites of the male genitalia are suitable for the classifications of species.

Reproductive interference and Satyrisation: mechanisms, outcomes and potential use for insect control

Reproductive Interference occurs when interactions between individuals from different species disrupt reproductive processes, resulting in a fitness cost to one or both parties involved. It is typically observed between individuals of closely related species, often upon secondary contact. In both vertebrates and invertebrates, Reproductive Interference is frequently referred to as 'Satyrisation'. It can manifest in various ways, ranging from blocking or reducing the efficacy of mating signals, through to negative effects of heterospecific copulations and the production of sterile or infertile hybrid offspring. The negative fitness effects of Satyrisation in reciprocal matings between species are often asymmetric and it is this aspect, which is most relevant to, and can offer utility in, pest management. In this review, we focus on Satyrisation and outline the mechanisms through which it can operate. We illustrate this by using test cases, and we consider the underlying reasons why the reproductive interactions that comprise Satyrisation occur. We synthesise the key factors affecting the expression of Satyrisation and explore how they have potential utility in developing new routes for the management and control of harmful insects. We consider how Satyrisation might interact with other control mechanisms, and conclude by outlining a framework for its use in control, highlighting some of the important next steps.

The scent chemistry of butterflies

Butterflies use structurally highly diverse volatile compounds for communication, in addition to visual signals. These compounds originate from plants or a formed de novo especially by male butterflies that possess specific scent organs.

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.

The Effects of GC-Biased Gene Conversion on Patterns of Genetic Diversity among and across Butterfly Genomes

Recombination reshuffles the alleles of a population through crossover and gene conversion. These mechanisms have considerable consequences on the evolution and maintenance of genetic diversity. Crossover, for example, can increase genetic diversity by breaking the linkage between selected and nearby neutral variants. Bias in favor of G or C alleles during gene conversion may instead promote the fixation of one allele over the other, thus decreasing diversity. Mutation bias from G or C to A and T opposes GC-biased gene conversion (gBGC). Less recognized is that these two processes may-when balanced-promote genetic diversity. Here, we investigate how gBGC and mutation bias shape genetic diversity patterns in wood white butterflies (Leptidea sp.). This constitutes the first in-depth investigation of gBGC in butterflies. Using 60 resequenced genomes from six populations of three species, we find substantial variation in the strength of gBGC across lineages. When modeling the balance of gBGC and mutation bias and comparing analytical results with empirical data, we reject gBGC as the main determinant of genetic diversity in these butterfly species. As alternatives, we consider linked selection and GC content. We find evidence that high values of both reduce diversity. We also show that the joint effects of gBGC and mutation bias can give rise to a diversity pattern which resembles the signature of linked selection. Consequently, gBGC should be considered when interpreting the effects of linked selection on levels of genetic diversity.

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.

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

The study of global biodiversity will greatly benefit from access to comprehensive DNA barcode libraries at continental scale, but such datasets are still very rare. Here, we assemble the first high-resolution reference library for European butterflies that provides 97% taxon coverage (459 species) and 22,306 COI sequences. We estimate that we captured 62% of the total haplotype diversity and show that most species possess a few very common haplotypes and many rare ones. Specimens in the dataset have an average 95.3% probability of being correctly identified. Mitochondrial diversity displayed elevated haplotype richness in southern European refugia, establishing the generality of this key biogeographic pattern for an entire taxonomic group. Fifteen percent of the species are involved in barcode sharing, but two thirds of these cases may reflect the need for further taxonomic research. This dataset provides a unique resource for conservation and for studying evolutionary processes, cryptic species, phylogeography, and ecology.

Evaluating DNA Barcoding for Species Identification and Discovery in European Gracillariid Moths

Gracillariidae is the most species-rich leaf-mining moth family with over 2,000 described species worldwide. In Europe, there are 263 valid named species recognized, many of which are difficult to identify using morphology only. Here we explore the use of DNA barcodes as a tool for identification and species discovery in European gracillariids. We present a barcode library including 6,791 COI sequences representing 242 of the 263 (92%) resident species. Our results indicate high congruence between morphology and barcodes with 91.3% (221/242) of European species forming monophyletic clades that can be identified accurately using barcodes alone. The remaining 8.7% represent cases of non-monophyly making their identification uncertain using barcodes. Species discrimination based on the Barcode Index Number system (BIN) was successful for 93% of species with 7% of species sharing BINs. We discovered as many as 21 undescribed candidate species, of which six were confirmed from an integrative approach; the other 15 require additional material and study to confirm preliminary evidence. Most of these new candidate species are found in mountainous regions of Mediterranean countries, the South-Eastern Alps and the Balkans, with nine candidate species found only on islands. In addition, 13 species were classified as deep conspecific lineages, comprising a total of 27 BINs with no intraspecific morphological differences found, and no known ecological differentiation. Double-digest restriction-site associated DNA sequencing (ddRAD) analysis showed strong mitonuclear discrepancy in four out of five species studied. This discordance is not explained by Wolbachia-mediated genetic sweeps. Finally, 26 species were classified as “unassessed species splits” containing 71 BINs and some involving geographical isolation or ecological specialization that will require further study to test whether they represent new cryptic species.

Overlooked cryptic diversity in Muschampia (Lepidoptera: Hesperiidae) adds two species to the European butterfly fauna

Cryptic species represent a challenge for documenting global biodiversity. Even in well-studied groups, such as European butterflies, the application of integrative approaches has allowed the recognition of an unexpected number of cryptic taxa. Here, we combine the analysis of mitochondrial (cytochrome c oxidase I, COI) and nuclear (internal transcribed spacer 2, ITS2) markers with geometric morphometrics of the male genitalia to study diversity within the butterfly Muschampia proto. The nuclear marker reveals three well-supported and deeply diverged lineages, which are also detected based on mitochondrial DNA, although the latter recovers one of them as paraphyletic with poor support. These lineages also present distinct male genital characters, which allow blind assignment of > 97% of specimens when applying a jackknife procedure. We conclude that M. proto comprises three cryptic species that started to differentiate ~2 Mya: M. proto, distributed in northern Africa, the Iberian Peninsula and southern France; Muschampia alta comb. & stat. nov., occurring in southern Italy and the Balkan Peninsula; and Muschampia proteides, present in the easternmost part of Europe, the Near East and Iran. This discovery adds two new species to the European butterfly fauna and highlights the necessity to continue investigating potential cryptic diversity.

Incomplete Sterility of Chromosomal Hybrids: Implications for Karyotype Evolution and Homoploid Hybrid Speciation

Heterozygotes for major chromosomal rearrangements such as fusions and fissions are expected to display a high level of sterility due to problems during meiosis. However, some species, especially plants and animals with holocentric chromosomes, are known to tolerate chromosomal heterozygosity even for multiple rearrangements. Here, we studied male meiotic chromosome behavior in four hybrid generations (F1–F4) between two chromosomal races of the Wood White butterfly Leptidea sinapis differentiated by at least 24 chromosomal fusions/fissions. Previous work showed that these hybrids were fertile, although their fertility was reduced as compared to crosses within chromosomal races. We demonstrate that (i) F1 hybrids are highly heterozygous with nearly all chromosomes participating in the formation of trivalents at the first meiotic division, and (ii) that from F1 to F4 the number of trivalents decreases and the number of bivalents increases. We argue that the observed process of chromosome sorting would, if continued, result in a new homozygous chromosomal race, i.e., in a new karyotype with intermediate chromosome number and, possibly, in a new incipient homoploid hybrid species. We also discuss the segregational model of karyotype evolution and the chromosomal model of homoploid hybrid speciation.

Evolution of multiple sex-chromosomes associated with dynamic genome reshuffling in Leptidea wood-white butterflies

Sex-chromosome systems tend to be highly conserved and knowledge about their evolution typically comes from macroevolutionary inference. Rapidly evolving complex sex-chromosome systems represent a rare opportunity to study the mechanisms of sex-chromosome evolution at unprecedented resolution. Three cryptic species of wood-white butterflies—Leptidea juvernica, L. sinapis and L. reali—have each a unique set of multiple sex-chromosomes with 3–4 W and 3–4 Z chromosomes. Using a transcriptome-based microarray for comparative genomic hybridisation (CGH) and a library of bacterial artificial chromosome (BAC) clones, both developed in L. juvernica, we identified Z-linked Leptidea orthologs of Bombyx mori genes and mapped them by fluorescence in situ hybridisation (FISH) with BAC probes on multiple Z chromosomes. In all three species, we determined synteny blocks of autosomal origin and reconstructed the evolution of multiple sex-chromosomes. In addition, we identified W homologues of Z-linked orthologs and characterised their molecular differentiation. Our results suggest that the multiple sex-chromosome system evolved in a common ancestor as a result of dynamic genome reshuffling through repeated rearrangements between the sex chromosomes and autosomes, including translocations, fusions and fissions. Thus, the initial formation of neo-sex chromosomes could not have played a role in reproductive isolation between these Leptidea species. However, the subsequent species-specific fissions of several neo-sex chromosomes could have contributed to their reproductive isolation. Then, significantly increased numbers of Z-linked genes and independent neo-W chromosome degeneration could accelerate the accumulation of genetic incompatibilities between populations and promote their divergence resulting in speciation.

Dissecting the Effects of Selection and Mutation on Genetic Diversity in Three Wood White (Leptidea) Butterfly Species

The relative role of natural selection and genetic drift in evolution is a major topic of debate in evolutionary biology. Most knowledge spring from a small group of organisms and originate from before it was possible to generate genome-wide data on genetic variation. Hence, it is necessary to extend to a larger number of taxonomic groups, descriptive and hypothesis-based research aiming at understanding the proximate and ultimate mechanisms underlying both levels of genetic polymorphism and the efficiency of natural selection. In this study, we used data from 60 whole-genome resequenced individuals of three cryptic butterfly species (Leptidea sp.), together with novel gene annotation information and population recombination data. We characterized the overall prevalence of natural selection and investigated the effects of mutation and linked selection on regional variation in nucleotide diversity. Our analyses showed that genome-wide diversity and rate of adaptive substitutions were comparatively low, whereas nonsynonymous to synonymous polymorphism and substitution levels were comparatively high in Leptidea, suggesting small long-term effective population sizes. Still, negative selection on linked sites (background selection) has resulted in reduced nucleotide diversity in regions with relatively high gene density and low recombination rate. We also found a significant effect of mutation rate variation on levels of polymorphism. Finally, there were considerable population differences in levels of genetic diversity and pervasiveness of selection against slightly deleterious alleles, in line with expectations from differences in estimated effective population sizes.

An updated checklist of the European Butterflies (Lepidoptera, Papilionoidea)

This paper presents an updated checklist of the butterflies of Europe, together with their original name combinations, and their occurrence status in each European country. According to this checklist, 496 species of the superfamily Papilionoidea occur in Europe. Changes in comparison with the last version (2.6.2) of Fauna Europaea are discussed. Compared to that version, 16 species are new additions, either due to cryptic species most of which have been discovered by molecular methods (13 cases) or due to discoveries of Asian species on the eastern border of the European territory in the Ural mountains (three cases). On the other hand, nine species had to be removed from the list, because they either do not occur in Europe or lost their species status due to new evidence. In addition, three species names had to be changed and 30 species changed their combination due to new evidence on phylogenetic relationships. Furthermore, minor corrections were applied to some authors’ names and years of publication. Finally, the namePolyommatusottomanusLefèbvre, 1831, which is threatened by its senior synonymLycaenalegeriFreyer, 1830, is declared anomen protectum, thereby conserving its name in the current combinationLycaenaottomana.

M. van Swaay, Vila R, Vliegenthart A, Wahlberg N, Verovnik R. An updated checklist of the European Butterflies (Lepidoptera, Papilionoidea)

This paper presents an updated checklist of the butterflies of Europe, together with their original name combinations, and their occurrence status in each European country. According to this checklist, 496 species of the superfamily Papilionoidea occur in Europe. Changes in comparison with the last version (2.6.2) of Fauna Europaea are discussed. Compared to that version, 16 species are new additions, either due to cryptic species most of which have been discovered by molecular methods (13 cases) or due to discoveries of Asian species on the eastern border of the European territory in the Ural mountains (three cases). On the other hand, nine species had to be removed from the list, because they either do not occur in Europe or lost their species status due to new evidence. In addition, three species names had to be changed and 30 species changed their combination due to new evidence on phylogenetic relationships. Furthermore, minor corrections were applied to some authors' names and years of publication. Finally, the name Polyommatusottomanus Lefèbvre, 1831, which is threatened by its senior synonym Lycaenalegeri Freyer, 1830, is declared a nomen protectum, thereby conserving its name in the current combination Lycaenaottomana.

Uniformity in premating reproductive isolation along an intraspecific cline

Premating reproductive isolation (RI) may reduce gene flow across populations that have differentiated in traits important for mate choice. Examining RI across genetic and phenotypic clines can inform the fundamental evolutionary processes that underlie population and lineage differentiation. We conducted female mate-choice studies across an intraspecific red-eyed treefrog cline in Costa Rica and Panama with 2 specific aims: (1) to characterize RI across the cline and examine the relationship between premating RI and genetic and phenotypic distance and (2) to evaluate our results within a broader evolutionary and taxonomic perspective through examination of other RI studies. We found that female red-eyed treefrogs prefer local males relative to non-local males, indicating that some premating RI has evolved in this system, but that preference strength is not associated with phenotypic or geographic distance. Our analysis of 65 other studies revealed no clear pattern between the strength of RI and geographic distribution (allopatry, parapatry, cline) or phenotypic distance, but revealed extreme variation and overlap in levels of intra- and interspecific levels of RI. This work contributes to a growing body of literature that examines intraspecific RI across a cline to understand the selective processes that shape evolutionary patterns at the earliest stages of divergence.

Versatility of multivalent orientation, inverted meiosis, and rescued fitness in holocentric chromosomal hybrids

Chromosomal rearrangements (e.g., fusions/fissions) have the potential to drive speciation. However, their accumulation in a population is generally viewed as unlikely, because chromosomal heterozygosity should lead to meiotic problems and aneuploid gametes. Canonical meiosis involves segregation of homologous chromosomes in meiosis I and sister chromatid segregation during meiosis II. In organisms with holocentric chromosomes, which are characterized by kinetic activity distributed along almost the entire chromosome length, this order may be inverted depending on their metaphase I orientation. Here we analyzed the evolutionary role of this intrinsic versatility of holocentric chromosomes, which is not available to monocentric ones, by studying F₁ to F₄ hybrids between two chromosomal races of the Wood White butterfly (Leptidea sinapis), separated by at least 24 chromosomal fusions/fissions. We found that these chromosomal rearrangements resulted in multiple meiotic multivalents, and, contrary to the theoretical prediction, the hybrids displayed relatively high reproductive fitness (42% of that of the control lines) and regular behavior of meiotic chromosomes. In the hybrids, we also discovered inverted meiosis, in which the first and critical stage of chromosome number reduction was replaced by the less risky stage of sister chromatid separation. We hypothesize that the ability to invert the order of the main meiotic events facilitates proper chromosome segregation and hence rescues fertility and viability in chromosomal hybrids, potentially promoting dynamic karyotype evolution and chromosomal speciation.

Chromosomal identification of cryptic species sharing their DNA barcodes: Polyommatus (Agrodiaetus) antidolus and P. (A.) morgani in Iran (Lepidoptera, Lycaenidae)

DNA barcoding has been suggested as a universal tool for molecular species identification; however, it cannot be applied in cases when morphologically similar species share their DNA barcodes due to the common ancestry or mitochondrial introgression. Here we analyze the karyotype of Polyommatus (Agrodiaetus) morgani (Le Cerf, 1909) from the region of its type locality in the southern Zagros Mountains in Iran, provide first chromosomal evidence for P. (A.) antidolus (Rebel, 1901) in Iran and demonstrate that these two species can be easily identified through analysis of their karyotypes whereas they share their mitochondrial barcodes.

Geographic isolation facilitates the evolution of reproductive isolation and morphological divergence

Geographic isolation is known to contribute to divergent evolution, resulting in unique phenotypes. Oftentimes morphologically distinct populations are found to be interfertile while reproductive isolation is found to exist within nominal morphological species revealing the existence of cryptic species. These disparities can be difficult to predict or explain especially when they do not reflect an inferred history of common ancestry which suggests that environmental factors affect the nature of ecological divergence. A series of laboratory experiments and observational studies were used to address what role biogeographic factors may play in the ecological divergence of Hyalella amphipods. It was found that geographic isolation plays a key role in the evolution of reproductive isolation and divergent morphology and that divergence cannot be explained by molecular genetic variation.

Rapid Increase in Genome Size as a Consequence of Transposable Element Hyperactivity in Wood-White (Leptidea) Butterflies

Characterizing and quantifying genome size variation among organisms and understanding if genome size evolves as a consequence of adaptive or stochastic processes have been long-standing goals in evolutionary biology. Here, we investigate genome size variation and association with transposable elements (TEs) across lepidopteran lineages using a novel genome assembly of the common wood-white (Leptidea sinapis) and population re-sequencing data from both L. sinapis and the closely related L. reali and L. juvernica together with 12 previously available lepidopteran genome assemblies. A phylogenetic analysis confirms established relationships among species, but identifies previously unknown intraspecific structure within Leptidea lineages. The genome assembly of L. sinapis is one of the largest of any lepidopteran taxon so far (643 Mb) and genome size is correlated with abundance of TEs, both in Lepidoptera in general and within Leptidea where L. juvernica from Kazakhstan has considerably larger genome size than any other Leptidea population. Specific TE subclasses have been active in different Lepidoptera lineages with a pronounced expansion of predominantly LINEs, DNA elements, and unclassified TEs in the Leptidea lineage after the split from other Pieridae. The rate of genome expansion in Leptidea in general has been in the range of four Mb/Million year (My), with an increase in a particular L. juvernica population to 72 Mb/My. The considerable differences in accumulation rates of specific TE classes in different lineages indicate that TE activity plays a major role in genome size evolution in butterflies and moths.

Evolutionary mechanisms of runaway chromosome number change in Agrodiaetus butterflies

Despite predictions of the classic, hybrid-sterility model of chromosomal speciation, some organisms demonstrate high rate of karyotype evolution. This rate is especially impressive in Agrodiaetus butterflies that rapidly evolved the greatest chromosome number diversity known in animal kingdom within a single subgenus. Here we analyzed karyotype evolution in Agrodiaetus using phylogenetic comparative methods. We found that chromosome numbers possess a strong phylogenetic signal. This disproves the chromosome megaevolution model that proposes multiple chromosome rearrangements to accumulate independently in each of closely related species. We found that Brownian motion gives a more adequate description of observed trait changes than Ornstein-Uhlenbeck model. This indicates that chromosome numbers evolve via random walk along branches of the phylogeny. We discovered a correlation between karyotype changes and phylogeny branch lengths. This gradual pattern is inconsistent with the hybrid-sterility model which, due to association of major chromosome changes with cladogenetic events, predicts a high degree of punctualism in karyotype evolution. Thus, low underdominace of chromosomal rearrangements and/or prevalence of the recombination-suppression model over the hybrid-sterility model of chromosome speciation are the most common engines of the runaway chromosome number change observed.

Panmixia and dispersal from the Mediterranean Basin to Macaronesian Islands of a macrolichen species

The Mediterranean region, comprising the Mediterranean Basin and the Macaronesian Islands, represents a center of diversification for many organisms. The genetic structure and connectivity of mainland and island microbial populations has been poorly explored, in particular in the case of symbiotic fungi. Here we investigated genetic diversity and spatial structure of the obligate outcrossing lichen-forming fungus Parmelina carporrhizans in the Mediterranean region. Using eight microsatellite and mating-type markers we showed that fungal populations are highly diverse but lack spatial structure. This is likely due to high connectivity and long distance dispersal of fungal spores. Consistent with low levels of linkage disequilibrium and lack of clonality, we detected both mating-type idiomorphs in all populations. Furthermore we showed that the Macaronesian Islands are the result of colonization from the Mediterranean Basin. The unidirectional gene flow, though, seemed not to be sufficient to counterbalance the effects of drift, resulting in comparatively allelic poor peripheral populations. Our study is the first to shed light on the high connectivity and lack of population structure in natural populations of a strictly sexual lichen fungus. Our data further support the view of the Macaronesian Islands as the end of the colonization road for this symbiotic ascomycete.

Karyosystematics and molecular taxonomy of the anomalous blue butterflies (Lepidoptera, Lycaenidae) from the Balkan Peninsula

The Balkan Peninsula represents one of the hottest biodiversity spots in Europe. However, the invertebrate fauna of this region is still insufficiently investigated, even in respect of such well-studied organisms as Lepidoptera. Here we use a combination of chromosomal, molecular and morphological markers to rearrange the group of so-called anomalous blue butterflies (also known as 'brown complex' of the subgenus Agrodiaetus Hübner, [1822] and as the Polyommatus (Agrodiaetus) admetus (Esper, 1783) species group) and to reveal its cryptic taxonomic structure. We demonstrate that Polyommatus aroaniensis (Brown, 1976) is not as widespread in the Balkans as was previously thought. In fact, it has a dot-like distribution range restricted to the Peloponnese Peninsula in South Greece. Polyommatus orphicus Kolev, 2005 is not as closely related to the Turkish species Polyommatus dantchenkoi (Lukhtanov & Wiemers, 2003) as was supposed earlier. Instead, it is a Balkan endemic represented by two subspecies: Polyommatus orphicus orphicus (Bulgaria) and Polyommatus orphicus eleniae Coutsis & De Prins, 2005 (Northern Greece). Polyommatus ripartii (Freyer, 1830) is represented in the Balkans by an endemic subspecies Polyommatus ripartii pelopi. The traditionally recognized Polyommatus admetus (Esper, 1783) is shown to be a heterogeneous complex and is divided into Polyommatus admetus sensu stricto (the Balkans and west Turkey) and Polyommatus yeranyani (Dantchenko & Lukhtanov, 2005) (east Turkey, Armenia, Azerbaijan and Iran). Polyommatus nephohiptamenos (Brown & Coutsis, 1978) is confirmed to be a species with a dot-like distribution range in Northern Greece. Finally, from Central Greece (Timfristos and Parnassos mountains) we describe Polyommatus timfristos Lukhtanov, Vishnevskaya & Shapoval, sp. n. which differs by its haploid chromosome number (n=38) from the closely related and morphologically similar Polyommatus aroaniensis (n=47-48) and Polyommatus orphicus (n=41-42). We provide chromosomal evidence for three separate south Balkan Pleistocene refugia (Peloponnesse, Central Greece and Northern Greece/South Bulgaria) and stress the biogeographic importance of Central Greece as a place of diversification. Then we argue that the data obtained have direct implications for butterfly karyology, taxonomy, biogeography and conservation.

Co-occurrence of ecologically equivalent cryptic species of spider wasps

Many cryptic species have been discovered in various taxonomic groups based on molecular phylogenetic analyses and mating experiments. Some sympatric cryptic species share equivalent resources, which contradicts the competitive exclusion principle. Two major theories have been proposed to explain the apparent lack of competitive exclusion, i.e. niche-based coexistence and neutral model, but a conclusive explanation is lacking. Here, we report the co-occurrence of cryptic spider wasp species appearing to be ecologically equivalent. Molecular phylogenetic analyses and mating experiments revealed that three phylogenetically closely related species are found sympatrically in Japan. These species share the same resources for larval food, and two of the species have the same niche for nesting sites, indicating a lack of competitive exclusion. This evidence may suggest that ecologically equivalent species can co-occur stably if their shared resources are sufficiently abundant that they cannot be over-exploited.

The sibling species Leptidea juvernica and L. sinapis (Lepidoptera, Pieridae) in the Balkan Peninsula: ecology, genetic structure, and morphological variation

The wood white butterfly Leptidea sinapis and its more recently discovered sibling species L. reali and L. juvernica have emerged as a model system for studying the speciation and evolution of cryptic species, as well as their ecological interactions in conditions of sympatry. Leptidea sinapis is widely distributed from Western Europe to Central Asia while the synmorphic L. juvernica and L. reali have allopatric distributions, both occurring in sympatry with L. sinapis and exhibiting an intricate, regionally variable ecological niche separation. Until now, the Balkan Peninsula remained one of the major unknowns in terms of distribution, genetic structure, and ecological preferences of the Leptidea triplet in Europe. We present new molecular and morphological data from a detailed survey of the region. Our DNA analysis suggests that the Balkan populations belong only to L. sinapis and L. juvernica, and that they are not genetically differentiated from other mainland populations. The distribution data reveal that L. sinapis is a widespread habitat generalist, while L. juvernica exhibits a transition from a habitat generalist in the western Balkans to a localised habitat specialist, confined to humid mountain habitats, in the east. The morphometric analysis of male genitalia and a comparison to data from other parts of the species' ranges suggest an interesting mosaic of regional morphological variation that is likely linked to environmental and ecological factors. We also demonstrate the usefulness of the genitalia scaling relations for analysing the morphological variation and discuss the implications for species identification.

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.

Genetic integrity of four species of Leptidea (Pieridae, Lepidoptera) as sampled in sympatry in West Siberia

In southern West Siberia, as many as four Leptidea Billberg, 1820 species are present sympatrically: Leptideaamurensis (Ménétriés, 1859), Leptideamorsei (Ménétriés, 1859), Leptideasinapis (Linnaeus, 1758) and Leptideajuvernica Williams, 1946. The two latter were recently recognised as nearly sibling species on morphological and molecular characters. Specimens intermediate as to their subtle diagnostic characters occurring in West Siberia and elsewhere were interpreted as resulted from limited introgression. This supposition was tested via populational morphological and molecular analysis of spring brood specimens of all the four species taken from a limited (4.5 × 0.2 km) area in the suburbs of Novosibirsk. The samples were analysed with respect to the genitalic morphology, external characters, three nuclear (CAD, H1 gene and ITS2) and one mitochondrial (COI) molecular markers, infection of the intracellular maternally inherited bacterial symbiont Wolbachia Hertig, 1836 and its wsp gene coding for a hypervariable surface protein. Interspecific variation of the nuclear CAD and ITS2 sequences and the mitochondrial COI gene in Leptideasinapis and Leptideajuvernica turned out concordant. The absence of molecular evidence of introgression suggests genetic integrity of these two species and allows their reliable identification by molecular characters. The genitalic (lengths of the saccus and valva) and external characters (wing pattern) of males overlap in Leptideasinapis and Leptideajuvernica, as identified by molecular markers and thus are not so helpful in actual species identification. Only the ductus bursae length showed no overlap and can be used for identification of females. The histone H1 gene appeared five times less variable over the four studied species than COI, and found to be identical in species Leptideasinapis and Leptideajuvernica. Wolbachia infection was found in all studied species. We identified three wsp variants of Wolbachia: 1) wsp-10 allele in Leptideaamurensis, Leptideasinapis, Leptideajuvernica; 2) a very similar wsp-687 allele in Leptideasinapis; and 3) wsp-688, highly divergent to the previous ones, in Leptideamorsei.

Dynamic karyotype evolution and unique sex determination systems in Leptidea wood white butterflies

Background

Chromosomal rearrangements have the potential to limit the rate and pattern of gene flow within and between species and thus play a direct role in promoting and maintaining speciation. Wood white butterflies of the genus Leptidea are excellent models to study the role of chromosome rearrangements in speciation because they show karyotype variability not only among but also within species. In this work, we investigated genome architecture of three cryptic Leptidea species (L. juvernica, L. sinapis and L. reali) by standard and molecular cytogenetic techniques in order to reveal causes of the karyotype variability.

Results

Chromosome numbers ranged from 2n = 85 to 91 in L. juvernica and 2n = 69 to 73 in L. sinapis (both from Czech populations) to 2n = 51 to 55 in L. reali (Spanish population). We observed significant differences in chromosome numbers and localization of cytogenetic markers (rDNA and H3 histone genes) within the offspring of individual females. Using FISH with the (TTAGG)_n telomeric probe we also documented the presence of multiple chromosome fusions and/or fissions and other complex rearrangements. Thus, the intraspecific karyotype variability is likely due to irregular chromosome segregation of multivalent meiotic configurations. The analysis of female meiotic chromosomes by GISH and CGH revealed multiple sex chromosomes: W₁W₂W₃Z₁Z₂Z₃Z₄ in L. juvernica, W₁W₂W₃Z₁Z₂Z₃ in L. sinapis and W₁W₂W₃W₄Z₁Z₂Z₃Z₄ in L. reali.

Conclusions

Our results suggest a dynamic karyotype evolution and point to the role of chromosomal rearrangements in the speciation of Leptidea butterflies. Moreover, our study revealed a curious sex determination system with 3–4 W and 3–4 Z chromosomes, which is unique in the Lepidoptera and which could also have played a role in the speciation process of the three Leptidea species.

Electronic supplementary material

The online version of this article (doi:10.1186/s12862-015-0375-4) contains supplementary material, which is available to authorized users.

Why do cryptic species tend not to co-occur? A case study on two cryptic pairs of butterflies

As cryptic diversity is being discovered, mostly thanks to advances in molecular techniques, it is becoming evident that many of these taxa display parapatric distributions in mainland and that they rarely coexist on islands. Genetic landscapes, haplotype networks and ecological niche modeling analyses were performed for two pairs of non-sister cryptic butterfly species, Aricia agestis-A. cramera and Polyommatus icarus-P. celina (Lycaenidae), to specifically assess non-coexistence on western Mediterranean islands, and to test potential causes producing such chequered distribution patterns. We show that the morphologically and ecologically equivalent pairs of species do not coexist on any of the studied islands, although nearly all islands are colonized by one of them. According to our models, the cryptic pairs displayed marked climatic preferences and 'precipitation during the driest quarter' was recovered as the most important climatic determinant. However, neither dispersal capacity, nor climatic or ecological factors fully explain the observed distributions across particular sea straits, and the existence of species interactions resulting in mutual exclusion is suggested as a necessary hypothesis. Given that the studied species are habitat generalists, feeding on virtually unlimited resources, we propose that reproductive interference, together with climatic preferences, sustain density-dependent mechanisms like "founder takes all" and impede coexistence on islands. Chequered distributions among cryptic taxa, both sister and non-sister, are common in butterflies, suggesting that the phenomenon revealed here could be important in determining biodiversity patterns.

CLIMBER: Climatic niche characteristics of the butterflies in Europe

Detailed information on species’ ecological niche characteristics that can be related to declines and extinctions is indispensable for a better understanding of the relationship between the occurrence and performance of wild species and their environment and, moreover, for an improved assessment of the impacts of global change. Knowledge on species characteristics such as habitat requirements is already available in the ecological literature for butterflies, but information about their climatic requirements is still lacking. Here we present a unique dataset on the climatic niche characteristics of 397 European butterflies representing 91% of the European species (see Appendix). These characteristics were obtained by combining detailed information on butterfly distributions in Europe (which also led to the ‘Distribution Atlas of Butterflies in Europe’) and the corresponding climatic conditions. The presented dataset comprises information for the position and breadth of the following climatic niche characteristics: mean annual temperature, range in annual temperature, growing degree days, annual precipitation sum, range in annual precipitation and soil water content. The climatic niche position is indicated by the median and mean value for each climate variable across a species’ range, accompanied by the 95% confidence interval for the mean and the number of grid cells used for calculations. Climatic niche breadth is indicated by the standard deviation and the minimum and maximum values for each climatic variable across a species’ range. Database compilation was based on high quality standards and the data are ready to use for a broad range of applications.

It is already evident that the information provided in this dataset is of great relevance for basic and applied ecology. Based on the species temperature index (STI, i.e. the mean temperature value per species), the community temperature index (CTI, i.e. the average STI value across the species in a community) was recently adopted as an indicator of climate change impact on biodiversity by the pan-European framework supporting the Convention on Biological Diversity (Streamlining European Biodiversity Indicators 2010) and has already been used in several scientific publications. The application potential of this database ranges from theoretical aspects such as assessments of past niche evolution or analyses of trait interdependencies to the very applied aspects of measuring, monitoring and projecting historical, ongoing and potential future responses to climate change using butterflies as an indicator.