Controversial patterns of Wolbachia infestation in the social parasitic Maculinea butterflies (Lepidoptera: Lycaenidae)

Evolutionary significant units (ESUs) and functional conservation units (FCUs) in the Phengaris (Maculinea) nausithous populations of East-Central Europe

The aim of this study was to analyse the level and structure of genetic and morphometric variation in the East-Central European populations of Phengaris (Maculinea) nausithous. On the basis of the life cycle of the species, we expected a low level of variation within the populations coupled with a relatively high level of differentiation among them. We also assumed that the differentiation among the populations has a regional pattern. In order to confirm these assumptions, we collected population samples from two regions within the Carpathian Basin (western Transdanubia and Transylvania) and from a region east of the Carpathian Mountains (Bukovina). The level of enzyme polymorphism and the amount of morphometric variation were investigated in the sampled populations. As the western (Transdanubia) and the eastern populations (Transylvania plus Bukovina) were suggested to belong to different subspecies (Rákosy et al. in Nota Lepidopterol 33:31-37, 2010), our working hypothesis was that these two groups of populations belong to different evolutionary significant units (ESUs). The results of all genetic analyses (PCA, Bayesian-clustering analysis and especially the UPGMA dendrogram) confirmed our assumption. The level of differentiation was higher between the western and eastern populations resulting in two monophyletic lineages of Phengaris (Maculinea) nausithous in East-Central Europe. This result suggests that these lineages can be considered as two different ESUs.

Phylogeography reveals the origin of the two phenological forms of large blue, Phengaris arion (Lepidoptera: Lycaenidae)

Our main goal was to investigate the phylogeography of the butterfly Phengaris arion to reveal the evolutionary origin of its ‘spring’ and ‘summer’ forms. Molecular analyses based on highly variable microsatellites, together with Wolbachia screening, were carried out on 34 populations in Europe. We found three well-defined genetic lineages of different origins: the Apennine, the central and the eastern. The highly distinct Apennine lineage is limited by the Alps and evaluated as an Evolutionary Significant Unit (ESU). Therefore, the taxon name ligurica, described from the Ligurian coast (Italy), should not be applied to denote the ‘summer form’ of the Pannonian region. The central lineage is limited by the Carpathians and the most eastern ranges of the Alps, and lacks major range fluctuations related to glaciations, although there is evidence for extra-Mediterranean refugia in the Carpathian Basin. The eastern clade could have had refugia in central Asia. Our results exclude the potential allopatric origin of the ‘spring’ and ‘summer’ arion, and support the hypothesis that the existence of the two forms could be a result of local adaptation to the distinctive phenology of host plant flowering which is manifested in the genetic differences between them. Wolbachia infection has been ruled out as a driver of sympatric speciation in P. arion.

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.

Incipient sympatric speciation via host race formation in Phengaris arion (Lepidoptera: Lycaenidae)

The plausibility of sympatric speciation is still debated despite increasing evidence, such as host races in insects. This speciation process may be occurring in the case of the two phenological forms of the obligatorily myrmecophilous Phengaris arion. The main goal of our research was to study the nature and causes of difference between these forms focusing primarily on the incipient speciation via host races. Molecular analyses based on highly variable microsatellites together with Wolbachia screening, male genitalia morphometrics and host ant studies were carried out on four syntopic sample pairs. Our results show that the two phenological forms of P. arion may meet the criteria for host plant races. They coexist in sympatry in certain parts of the species range which is allowed by the adaptation to the distinct phenology of the host plants. Negative selection acts against the intermediate individuals which are on the wing in the inappropriate time frame. Thus, disruptive selection affects and produces bimodal distributions of phenotypes. However, the phenology of food plants is not entirely distinct and fluctuates year by year. Therefore, the two forms can exchange genes occasionally depending on the length of the time slot when they can meet with each other. Consequently, the reproductive isolation could not be completed and the existence of the two arion forms may represent only an incipient stage of sympatric speciation. It is also clear that Wolbachia is likely not a driver of sympatric speciation in this case.

Assessment of the role of Wolbachia in mtDNA paraphyly and the evolution of unisexuality in Calligrapha (Coleoptera: Chrysomelidae)

Calligrapha is a New World leaf beetle genus that includes several unisexual species in northeastern North America. Each unisexual species had an independent hybrid origin involving different combinations of bisexual species. However, surprisingly, they all cluster in a single mtDNA clade and with some individuals of their parental species, which are in turn deeply polyphyletic for mtDNA. This pattern is suggestive of a selective sweep which, together with mtDNA taxonomic incongruence and occurrence of unisexuality in Calligrapha, led to hypothesize that Wolbachia might be responsible. I tested this hypothesis studying the correlation between diversity of Wolbachia and well-established mtDNA lineages in >500 specimens of two bisexual species of Calligrapha and their derived unisexual species. Wolbachia appears highly prevalent (83.4%), and fifteen new supergroup-A strains of the bacteria are characterized, belonging to three main classes: wCallA, occupying the whole species ranges, and wCallB and wCallC, narrowly parapatric, infecting beetles with highly divergent mtDNAs where they coexist. Most beetles (71.6%) carried double infections of wCallA with another sequence class. Bayesian inference of ancestral character states and association tests between bacterial diversity and the mtDNA genealogy show that each mtDNA lineage of Calligrapha has specific types of infection. Moreover, shifts can be explained by horizontal or vertical transfer from local populations to an expanding lineage and cytoplasmic incompatibility between wCallB and wCallC types, suggesting that the symbionts hitchhike with the host and are not responsible for selective mtDNA sweeps. Lack of evidence for sweeps and the fact that individuals in the unisexual clade are uninfected or infected by the widespread wCallA type indicate that Wolbachia does not induce unisexuality in Calligrapha, although they may manipulate host reproduction through cytoplasmic incompatibility.

Molecular diversity of Wolbachia in Lepidoptera: Prevalent allelic content and high recombination of MLST genes

Wolbachia are common endosymbiotic bacteria of Arthropoda and Nematoda that are ordinarily transmitted vertically in host lineages through the egg cytoplasm. Despite the great interest in the Wolbachia symbiont, many issues of its biology remain unclear, including its evolutionary history, routes of transfer among species, and the molecular mechanisms underlying the symbiont's effect on its host. In this report, we present data relating to Wolbachia infection in 120 species of 13 Lepidoptera families, mostly butterflies, from West Siberian localities based on Multilocus sequence typing (MLST) and the wsp locus and perform a comprehensive survey of the distribution of Wolbachia and its genetic diversity in Lepidoptera worldwide. We observed a high infection incidence in the studied region; this finding is probably also true for other temperate latitude regions because many studied species have broad Palearctic and even Holarctic distribution. Although 40 new MLST alleles and 31 new STs were described, there was no noticeable difference in the MLST allele content in butterflies and probably also in moths worldwide. A genetic analysis of Wolbachia strains revealed the MLST allele core in lepidopteran hosts worldwide, viz. the ST-41 allele content. The key finding of our study was the detection of rampant recombination among MLST haplotypes. High rates of homologous recombination between Wolbachia strains indicate a substantial contribution of genetic exchanges to the generation of new STs. This finding should be considered when discussing issues related to the reconstruction of Wolbachia evolution, divergence time, and the routes of Wolbachia transmission across arthropod hosts.

Host plant use drives genetic differentiation in syntopic populations of Maculinea alcon

The rare socially parasitic butterfly Maculinea alcon occurs in two forms, which are characteristic of hygric or xeric habitats and which exploit different host plants and host ants. The status of these two forms has been the subject of considerable controversy. Populations of the two forms are usually spatially distinct, but at Răscruci in Romania both forms occur on the same site (syntopically). We examined the genetic differentiation between the two forms using eight microsatellite markers, and compared with a nearby hygric site, Şardu. Our results showed that while the two forms are strongly differentiated at Răscruci, it is the xeric form there that is most similar to the hygric form at Şardu, and Bayesian clustering algorithms suggest that these two populations have exchanged genes relatively recently. We found strong evidence for population substructuring, caused by high within host ant nest relatedness, indicating very limited dispersal of most ovipositing females, but not association with particular host ant species. Our results are consistent with the results of larger scale phylogeographic studies that suggest that the two forms represent local ecotypes specialising on different host plants, each with a distinct flowering phenology, providing a temporal rather than spatial barrier to gene flow.