Genetic structure of sibling butterfly species affected by Wolbachia infection sweep: evolutionary and biogeographical implications

Arthropod mtDNA paraphyly: a case study of introgressive origin

Mitochondrial paraphyly between arthropod species is not uncommon and has been speculated to largely be the result of incomplete lineage sorting (ILS) of ancestral variation within the common ancestor of both species, with hybridization playing only a minor role. However, in the absence of comparable nuclear genetic data, the relative roles of ILS and hybridization in explaining mitochondrial DNA (mtDNA) paraphyly remain unclear. Hybridization itself is a multifaceted gateway to mtDNA paraphyly, which may lead to paraphyly across both the nuclear and mitochondrial genomes, or paraphyly that is largely restricted to the mitochondrial genome. These different outcomes will depend upon the frequency of hybridization, its demographic context, and the extent to which mtDNA is subject to direct selection, indirect selection, or neutral processes. Here, we describe extensive mtDNA paraphyly between two species of iron-clad beetle (Zopheridae) and evaluate competing explanations for its origin. We first test between hypotheses of ILS and hybridization, revealing strong nuclear genetic differentiation between species, but with the complete replacement of Tarphius simplex mtDNA through the introgression of at least 5 mtDNA haplotypes from T. canariensis. We then contrast explanations of direct selection, indirect selection, or genetic drift for observed patterns of mtDNA introgression. Our results highlight how introgression can lead to complex patterns of mtDNA paraphyly across arthropod species, while simultaneously revealing the challenges for understanding the selective or neutral drivers that underpin such patterns.

Beyond gene flow: (non)-parallelism of secondary contact in a pair of highly differentiated sibling species

Replicated secondary contact zones can provide insights into the barriers to gene flow that are important during speciation and can reveal to which degree secondary contact may result in similar evolutionary outcomes. Here, we studied two secondary contact zones between highly differentiated Alpine butterflies of the genus Erebia using whole-genome resequencing data. We assessed the genomic relationships between populations and species and found hybridization to be rare, with no to little current or historical introgression in either contact zone. There are large similarities between contact zones, consistent with an allopatric origin of interspecific differentiation, with no indications for ongoing reinforcing selection. Consistent with expected reduced effective population size, we further find that scaffolds related to the Z-chromosome show increased differentiation compared to the already high levels across the entire genome, which could also hint towards a contribution of the Z chromosome to species divergence in this system. Finally, we detected the presence of the endosymbiont Wolbachia, which can cause reproductive isolation between its hosts, in all E. cassioides, while it appears to be fully or largely absent in contact zone populations of E. tyndarus. We discuss how this rare pattern may have arisen and how it may have affected the dynamics of speciation upon secondary contact.

Validating a Mitochondrial Sweep Accompanying the Rapid Spread of a Maternally Inherited Symbiont

Maternally inherited symbiotic bacteria that interfere with the reproduction of their hosts can contribute to selective sweeps of mitochondrial haplotypes through hitch-hiking or coordinate inheritance of cytoplasmic bacteria and host mitochondria. The sweep will be manifested by genetic variations of mitochondrial genomic DNA of symbiont-infected hosts relative to their uninfected counterparts. In particular, at the population level, infected specimens will show a reduced mitochondrial DNA polymorphism compared to that in the nuclear DNA. This may challenge the use of mitochondrial DNA sequences as neutral genetic markers, as the mitochondrial patterns will reflect the evolutionary history of parasitism, rather than the sole evolutionary history of the host. Here, I describe a detailed step-by-step procedure to infer the occurrence and timing of symbiont-induced mitochondrial sweeps in host species.

Diversity of Wolbachia infection and its influence on mitochondrial DNA variation in the diamondback moth, Plutella xylostella

Plutella xylostella is a pest that severely damages cruciferous vegetables worldwide and has been shown to be infected with the maternally inherited bacteria Wolbachia, with the main infected strain was plutWB1. In this study, we performed a large-scale global sampling of P. xylostella and amplified 3 mtDNA genes of P. xylostella and 6 Wolbachia genes to analyze the infection status, diversity of Wolbachia in P. xylostella, and its effect on mtDNA variation in P. xylostella. This study provides a conservative estimate of Wolbachia infection rates in P. xylostella, which was found to be 7% (104/1440). The ST 108 (plutWB1) was shared among butterfly species and the moth species P. xylostella, revealing that Wolbachia strain plutWB1 acquisition in P. xylostella may be through horizontal transmission. The Parafit analyses indicated a significant association between Wolbachia and Wolbachia-infected P. xylostella individuals, and individuals infected with plutWB1 tended to cluster in the basal positions of the phylogenetic tree based on the mtDNA data. Additionally, Wolbachia infections were associated with increased mtDNA polymorphism in the infected P. xylostella population. These data suggest that Wolbachia endosymbionts may have a potential effect on mtDNA variation of P. xylostella.

A continental-scale survey of Wolbachia infections in blue butterflies reveals evidence of interspecific transfer and invasion dynamics

Infections by maternally inherited bacterial endosymbionts, especially Wolbachia, are common in insects and other invertebrates but infection dynamics across species ranges are largely under studied. Specifically, we lack a broad understanding of the origin of Wolbachia infections in novel hosts, and the historical and geographical dynamics of infections that are critical for identifying the factors governing their spread. We used Genotype-by-Sequencing data from previous population genomics studies for range-wide surveys of Wolbachia presence and genetic diversity in North American butterflies of the genus Lycaeides. As few as one sequence read identified by assembly to a Wolbachia reference genome provided high accuracy in detecting infections in host butterflies as determined by confirmatory PCR tests, and maximum accuracy was achieved with a threshold of only 5 sequence reads per host individual. Using this threshold, we detected Wolbachia in all but 2 of the 107 sampling localities spanning the continent, with infection frequencies within populations ranging from 0% to 100% of individuals, but with most localities having high infection frequencies (mean = 91% infection rate). Three major lineages of Wolbachia were identified as separate strains that appear to represent 3 separate invasions of Lycaeides butterflies by Wolbachia. Overall, we found extensive evidence for acquisition of Wolbachia through interspecific transfer between host lineages. Strain wLycC was confined to a single butterfly taxon, hybrid lineages derived from it, and closely adjacent populations in other taxa. While the other 2 strains were detected throughout the rest of the continent, strain wLycB almost always co-occurred with wLycA. Our demographic modeling suggests wLycB is a recent invasion. Within strain wLycA, the 2 most frequent haplotypes are confined almost exclusively to separate butterfly taxa with haplotype A1 observed largely in Lycaeides melissa and haplotype A2 observed most often in Lycaeides idas localities, consistent with either cladogenic mode of infection acquisition from a common ancestor or by hybridization and accompanying mutation. More than 1 major Wolbachia strain was observed in 15 localities. These results demonstrate the utility of using resequencing data from hosts to quantify Wolbachia genetic variation and infection frequency and provide evidence of multiple colonizations of novel hosts through hybridization between butterfly lineages and complex dynamics between Wolbachia strains.

Wolbachia in Black Spiny Whiteflies and Their New Parasitoid Wasp in Japan: Evidence of the Distinct Infection Status on Aleurocanthus camelliae Cryptic Species Complex

Wolbachia, an alphaproteobacterial reproductive parasite, can cause profound mitochondrial divergence in insects, which might eventually be a part of cryptic speciation. Aleurocanthus camelliae is a cryptic species complex consisting of several morphospecies and/or haplotypes that are genetically different but morphologically indistinctive. However, little is known about the Wolbachia infection status in these tea and Citrus pests. Thus, this study aimed to profile the diversity and phenotypic characteristics of Wolbachia natural infections in the A. camelliae cryptic species complex. A monophyletic strain of Wolbachia that infected the A. camelliae cryptic species complex (wAlec) with different patterns was discovered. Whiteflies that are morphologically identical to Aleurocanthus spiniferus (Aleurocanthus cf. A. spiniferus in Eurya japonica and A. spiniferus in Citrus) were grouped into uninfected populations, whereas the fixed infection was detected in A. camelliae B1 from Theaceae. The rapid evolution of wAlec was also found to occur through a high recombination event, which produced subgroups A and B in wAlec. It may also be associated with the non-cytoplasmic incompatibility (CI) phenotype of wAlec due to undetectable CI-related genes from phage WO (WOAlec). The current discovery of a novel cryptic species of A. camelliae led to a discussion about the oscillation hypothesis, which may provide insights on cryptic speciation, particularly on how specialization and host expansion have been recorded among these species. This study also identified a parasitoid wasp belonging to the genus Eretmocerus in A. camelliae, for the first time in Japan.

Wolbachia infection and genetic diversity of Italian populations of Philaenus spumarius, the main vector of Xylella fastidiosa in Europe

Philaenus spumarius is a cosmopolitan species that has become a major threat to European agriculture being recognized as the main vector of the introduced plant pathogen Xylella fastidiosa, the agent of the “olive quick decline syndrome”, a disease which is devastating olive orchards in southern Italy. Wolbachia are bacterial symbionts of many insects, frequently as reproductive parasites, sometime by establishing mutualistic relationships, able to spread within host populations. Philaenus spumarius harbors Wolbachia, but the role played by this symbiont is unknown and data on the infection prevalence within host populations are limited. Here, the Wolbachia infection rate was analyzed in relation to the geographic distribution and the genetic diversity of the Italian populations of P. spumarius. Analysis of the COI gene sequences revealed a geographically structured distribution of the three main mitochondrial lineages of P. spumarius. Wolbachia was detected in half of the populations sampled in northern Italy where most individuals belonged to the western-Mediterranean lineage. All populations sampled in southern and central Italy, where the individuals of the eastern-Mediterranean lineage were largely prevalent, were uninfected. Individuals of the north-eastern lineage were found only in populations from the Alps in the northernmost part of Italy, at high altitudes. In this area, Wolbachia infection reached the highest prevalence, with no difference between north-eastern and western-Mediterranean lineage. Analysis of molecular diversity of COI sequences suggested no significant effect of Wolbachia on population genetics of P. spumarius. Using the MLST approach, six new Wolbachia sequence types were identified. Using FISH, Wolbachia were observed within the host’s reproductive tissues and salivary glands. Results obtained led us to discuss the role of Wolbachia in P. spumarius, the factors influencing the geographic distribution of the infection, and the exploitation of Wolbachia for the control of the vector insect to reduce the spread of X. fastidiosa.

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

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

Incongruent molecular and morphological variation in the crab spider Synemaglobosum (Araneae, Thomisidae) in Europe

Establishing species boundaries is one of the challenges taxonomists around the world have been tackling for centuries. The relation between intraspecific and interspecific variability is still under discussion and in many taxa it remains understudied. Here the hypothesis of single versus multiple species of the crab spider Synemaglobosum (Fabricius) is tested. The wide distribution range as well as its high morphological variability makes this species an interesting candidate for re-evaluation using an integrative approach. This study combines information from barcoding, phylogenetic reconstruction based on mitochondrial CO1 and ITS2 of more than 60 specimens collected over a wide range of European localities, and morphology. The findings show deep clades with up to 6% mean pairwise distance in the CO1 barcode without any biogeographical pattern. The nuclear ITS2 gene did not support the CO1 clades. Morphological assessment of somatic and genital characters in males and females and a morphometric analysis of the male palp uncovered high intraspecific variation that does not match the CO1 or ITS2 phylogenies or biogeography either. Screening for endosymbiotic Wolbachia bacteria was conducted and only a single infected specimen was found. Several scenarios might explain these inconsistent patterns. While the deep divergences in the barcoding marker might suggest cryptic or ongoing speciation or geographical isolation in the past, the lack of congruent variation in the nuclear ITS2 gene or the studied morphological character systems, especially the male palp, indicates that S.globosum might simply be highly polymorphic both in terms of its mtDNA and morphology. Therefore, more data on ecology and behaviour and full genome sequences are necessary to ultimately resolve this taxonomically intriguing case.

Wolbachia-driven selective sweep in a range expanding insect species

Background

Evolutionary processes can cause strong spatial genetic signatures, such as local loss of genetic diversity, or conflicting histories from mitochondrial versus nuclear markers. Investigating these genetic patterns is important, as they may reveal obscured processes and players. The maternally inherited bacterium Wolbachia is among the most widespread symbionts in insects. Wolbachia typically spreads within host species by conferring direct fitness benefits, and/or by manipulating its host reproduction to favour infected over uninfected females. Under sufficient selective advantage, the mitochondrial haplotype associated with the favoured maternally-inherited symbiotic strains will spread (i.e. hitchhike), resulting in low mitochondrial genetic variation across the host species range.

Method

The common bluetail damselfly (Ischnura elegans: van der Linden, 1820) has recently emerged as a model organism for genetics and genomic signatures of range expansion during climate change. Although there is accumulating data on the consequences of such expansion on the genetics of I. elegans, no study has screened for Wolbachia in the damselfly genus Ischnura. Here, we present the biogeographic variation in Wolbachia prevalence and penetrance across Europe and Japan (including samples from 17 populations), and from close relatives in the Mediterranean area (i.e. I. genei: Rambur, 1842; and I. saharensis: Aguesse, 1958).

Results

Our data reveal (a) multiple Wolbachia-strains, (b) potential transfer of the symbiont through hybridization, (c) higher infection rates at higher latitudes, and (d) reduced mitochondrial diversity in the north-west populations, indicative of hitchhiking associated with the selective sweep of the most common strain. We found low mitochondrial haplotype diversity in the Wolbachia-infected north-western European populations (Sweden, Scotland, the Netherlands, Belgium, France and Italy) of I. elegans, and, conversely, higher mitochondrial diversity in populations with low penetrance of Wolbachia (Ukraine, Greece, Montenegro and Cyprus). The timing of the selective sweep associated with infected lineages was estimated between 20,000 and 44,000 years before present, which is consistent with the end of the last glacial period about 20,000 years.

Conclusions

Our findings provide an example of how endosymbiont infections can shape spatial variation in their host evolutionary genetics during postglacial expansion. These results also challenge population genetic studies that do not consider the prevalence of symbionts in many insects, which we show can impact geographic patterns of mitochondrial genetic diversity.

Sympatric occurrence of deeply diverged mitochondrial DNA lineages in Siberian geometrid moths (Lepidoptera: Geometridae): cryptic speciation, mitochondrial introgression, secondary admixture or effect of Wolbachia?

The divergent sympatric mitochondrial lineages within traditionally recognized species present a challenge regularly faced by taxonomists and evolutionary biologists. We encountered this problem when studying the Siberian geometrid moths, Alcis deversata and Thalera chlorosaria. Within each of these species we found two deeply diverged mitochondrial lineages that demonstrated a level of genetic differentiation exceeding the standard interspecific DNA barcode threshold. Using analyses of nuclear genes, morphology, ecological preferences and Wolbachia endosymbionts, we tested five hypotheses that might explain the mitochondrial pattern observed: cryptic speciation, ancestral polymorphism, interspecific mitochondrial introgression, secondary admixture of allopatrically evolved populations and an effect of intracellular Wolbachia endosymbionts. We demonstrate that in A. deversata and Th. chlorosaria the mitochondrial differences are not correlated with differences in nuclear genes, morphology, ecology and Wolbachia infection status, thus not supporting the hypothesis of cryptic species and an effect of Wolbachia. Mitochondrial introgression can lead to a situation in which one species has both its own mitochondrial lineage and the lineage obtained from another species. We found this situation in the species pair Alcis repandata and Alcis extinctaria. We conclude that the mitochondrial heterogeneity in A. deversata and Th. chlorosaria is most likely to be attributable to the secondary admixture of allopatrically evolved populations.

Population structure, adaptation and divergence of the meadow spittlebug, Philaenus spumarius (Hemiptera, Aphrophoridae), revealed by genomic and morphological data

Understanding patterns of population differentiation and gene flow in insect vectors of plant diseases is crucial for the implementation of management programs of disease. We investigated morphological and genome-wide variation across the distribution range of the spittlebug Philaenus spumarius (Linnaeus, 1758) (Hemiptera, Auchenorrhyncha, Aphrophoridae), presently the most important vector of the plant pathogenic bacterium Xylella fastidiosa Wells et al., 1987 in Europe. We found genome-wide divergence between P. spumarius and a very closely related species, P. tesselatus Melichar, 1899, at RAD sequencing markers. The two species may be identified by the morphology of male genitalia but are not differentiated at mitochondrial COI, making DNA barcoding with this gene ineffective. This highlights the importance of using integrative approaches in taxonomy. We detected admixture between P. tesselatus from Morocco and P. spumarius from the Iberian Peninsula, suggesting gene-flow between them. Within P. spumarius, we found a pattern of isolation-by-distance in European populations, likely acting alongside other factors restricting gene flow. Varying levels of co-occurrence of different lineages, showing heterogeneous levels of admixture, suggest other isolation mechanisms. The transatlantic populations of North America and Azores were genetically closer to the British population analyzed here, suggesting an origin from North-Western Europe, as already detected with mitochondrial DNA. Nevertheless, these may have been produced through different colonization events. We detected SNPs with signatures of positive selection associated with environmental variables, especially related to extremes and range variation in temperature and precipitation. The population genomics approach provided new insights into the patterns of divergence, gene flow and adaptation in these spittlebugs and led to several hypotheses that require further local investigation.

Diversification, selective sweep, and body size in the invasive Palearctic alfalfa weevil infected with Wolbachia

The alfalfa weevil Hypera postica, native to the Western Palearctic, is an invasive legume pest with two divergent mitochondrial clades in its invading regions, the Western clade and the Eastern/Egyptian clade. However, knowledge regarding the native populations is limited. The Western clade is infected with the endosymbiotic bacteria Wolbachia that cause cytoplasmic incompatibility in host weevils. Our aim was to elucidate the spatial genetic structure of this insect and the effect of Wolbachia on its population diversity. We analyzed two mitochondrial and two nuclear genes of the weevil from its native ranges. The Western clade was distributed in western/central Europe, whereas the Eastern/Egyptian clade was distributed from the Mediterranean basin to central Asia. Intermediate mitotypes were found from the Balkans to central Asia. Most Western clade individuals in western Europe were infected with an identical Wolbachia strain. Mitochondrial genetic diversity of the infected individuals was minimal. The infected clades demonstrated a higher nonsynonymous/synonymous substitution rate ratio than the uninfected clades, suggesting a higher fixation of nonsynonymous mutations due to a selective sweep by Wolbachia. Trans-Mediterranean and within-European dispersal routes were supported. We suggest that the ancestral populations diversified by geographic isolation due to glaciations and that the diversity was reduced in the west by a recent Wolbachia-driven sweep(s). The intermediate clade exhibited a body size and host plant that differed from the other clades. Pros and cons of the possible use of infected-clade males to control uninfected populations are discussed.

Wolbachia affects mitochondrial population structure in two systems of closely related Palaearctic blue butterflies

The bacterium Wolbachia infects many insect species and spreads by diverse vertical and horizontal means. As co-inherited organisms, these bacteria often cause problems in mitochondrial phylogeny inference. The phylogenetic relationships of many closely related Palaearctic blue butterflies (Lepidoptera: Lycaenidae: Polyommatinae) are ambiguous. We considered the patterns of Wolbachia infection and mitochondrial diversity in two systems: Aricia agestis/Aricia artaxerxes and the Pseudophilotes baton species complex. We sampled butterflies across their distribution ranges and sequenced one butterfly mitochondrial gene and two Wolbachia genes. Both butterfly systems had uninfected and infected populations, and harboured several Wolbachia strains. Wolbachia was highly prevalent in A. artaxerxes and the host's mitochondrial structure was shallow, in contrast to A. agestis. Similar bacterial alleles infected both Aricia species from nearby sites, pointing to a possible horizontal transfer. Mitochondrial history of the P. baton species complex mirrored its Wolbachia infection and not the taxonomical division. Pseudophilotes baton and P. vicrama formed a hybrid zone in Europe. Wolbachia could obscure mitochondrial history, but knowledge on the infection helps us to understand the observed patterns. Testing for Wolbachia should be routine in mitochondrial DNA studies.

Host plants and Wolbachia shape the population genetics of sympatric herbivore populations

Changing climate and land-use practices have the potential to bring previously isolated populations of pest insects into new sympatry. This heightens the need to better understand how differing patterns of host-plant association, and unique endosymbionts, serve to promote genetic isolation or integration. We addressed these factors in populations of potato psyllid, Bactericera cockerelli (Šulc), a generalist herbivore that vectors a bacterial pathogen (Candidatus Liberibacter solanacearum, causal pathogen of zebra chip disease) of potato (Solanum tuberosum L.). Genome-wide SNP data revealed two major genetic clusters-psyllids collected from potato crops were genetically similar to psyllids found on a common weed, Lycium spp., but dissimilar from those found on another common non-crop host, Solanum dulcamara L. Most psyllids found on Lycium spp. and potato represented a single mitochondrial cytochrome oxidase I (COI) haplotype that has been suggested to not be native to the region, and whose arrival may have been concurrent with zebra chip disease first emerging. The putatively introduced COI haplotype usually co-occurred with endosymbiotic Wolbachia, while the putatively resident COI haplotype generally did not. Genetic intermediates between the two genetic populations of insects were rare, consistent with recent sympatry or reproductive isolation, although admixture patterns of apparent hybrids were consistent with introgression of genes from introduced into resident populations. Our results suggest that both host-plant associations and endosymbionts are shaping the population genetic structure of sympatric psyllid populations associated with different non-crop hosts. It is of future interest to explicitly examine vectorial capacity of the two populations and their potential hybrids, as population structure and hybridization might alter regional vector capacity and disease outbreaks.

Wolbachia have made it twice: Hybrid introgression between two sister species of Eurema butterflies

Wolbachia, cytoplasmically inherited endosymbionts of arthropods, are known to hijack their host reproduction in various ways to increase their own vertical transmission. This may lead to the selective sweep of associated mitochondria, which can have a large impact on the evolution of mitochondrial lineages. In Japan, two different Wolbacahia strains (wCI and wFem) are found in two sister species of pierid butterflies, Eurema mandarina and Eurema hecabe. In both species, females infected with wCI (C females) produce offspring with a nearly 1:1 sex ratio, while females infected with both wCI and wFem (CF females) produce all-female offspring. Previous studies have suggested the historical occurrence of hybrid introgression in C individuals between the two species. Furthermore, hybrid introgression in CF individuals is suggested by the distinct mitochondrial lineages between C females and CF females of E. mandarina. In this study, we performed phylogenetic analyses based on nuclear DNA and mitochondrial DNA markers of E. hecabe with previously published data on E. mandarina. We found that the nuclear DNA of this species significantly diverged from that of E. mandarina. By contrast, mitochondrial DNA haplotypes comprised two clades, mostly reflecting Wolbachia infection status rather than the individual species. Collectively, our results support the previously suggested occurrence of two independent historical events wherein the cytoplasms of CF females and C females moved between E. hecabe and E. mandarina through hybrid introgression.

How consistent is RAD-seq divergence with DNA-barcode based clustering in insects?

Promoted by the barcoding approach, mitochondrial DNA is more than ever used as a molecular marker to identify species boundaries. Yet, it has been repeatedly argued that it may be poorly suited for this purpose, especially in insects where mitochondria are often associated with invasive intracellular bacteria that may promote their introgression. Here, we inform this debate by assessing how divergent nuclear genomes can be when mitochondrial barcodes indicate very high proximity. To this end, we obtained RAD-seq data from 92 barcode-based species-like units (operational taxonomic units [OTUs]) spanning four insect orders. In 100% of the cases, the observed median nuclear divergence was lower than 2%, a value that was recently estimated as one below which nuclear gene flow is not uncommon. These results suggest that although mitochondria may occasionally leak between species, this process is rare enough in insects to make DNA barcoding a reliable tool for clustering specimens into species-like units.

Multiplex PCR for identification of two butterfly sister species: Eurema mandarina and Eurema hecabe

Objective

In insects, closely related species are often difficult or impossible to distinguish solely by morphological traits. Mitochondrial DNA (mtDNA) markers are often useful and reliable for distinguishing closely related species. However, useful mtDNA markers can be unavailable, particularly when such species pairs experienced hybrid introgression in the past. Although polymorphic nuclear DNA markers would be necessary to distinguish such species pairs, recombination, multiple copies, and slower mutation rates of the nuclear DNA compared with those of mtDNA often make it challenging. The objective of this study was to develop a multiplex polymerase chain reaction that can reliably amplify and distinguish the Tpi sequences of Eurema mandarina and Eurema hecabe.

Results

We successfully analyzed the nucleotide sequences of the Z chromosome-linked triose phosphate isomerase (Tpi) gene to develop a multiplex polymerase chain reaction (PCR) that amplified ca. 120-bp products for E. mandarina and ca. 375-bp products for E. hecabe. We suggest that multiplex PCR using Tpi with appropriately designed primers can be used to accurately and reliably distinguish between other closely related Lepidoptera species.

Genetic Diversity and Wolbachia Infection Patterns in a Globally Distributed Invasive Ant

Understanding the phylogeographic history of an invasive species may facilitate reconstructing the history and routes of its invasion. The longhorn crazy ant, Paratrechina longicornis, is a ubiquitous agricultural and household pest throughout much of the tropics and subtropics, but little is known about the history of its spread. Here, we examine worldwide genetic variation in P. longicornis and its associated Wolbachia bacterial symbionts. Analyses of mtDNA sequences of 248 P. longicornis workers (one per colony) from 13 geographic regions reveal two highly diverged mtDNA clades that co-occur in most of the geographic regions. These two mtDNA clades are associated with different Wolbachia infection patterns, but are not congruent with patterns of nDNA (microsatellite) variation. Multilocus sequence typing reveals two distinct Wolbachia strains in P. longicornis, namely, wLonA and wLonF. The evolutionary histories of these two strains differ; wLonA appears to be primarily transmitted maternally, and patterns of mtDNA and nDNA variation and wLonA infection status are consistent with a relatively recent Wolbachia-induced selective sweep. In contrast, the observed patterns of mtDNA variation and wLonF infections suggest frequent horizontal transfer and losses of wLonF infections. The lack of nDNA structure among sampled geographic regions coupled with the finding that numerous mtDNA haplotypes are shared among regions implies that inadvertent long-distance movement through human commerce is common in P. longicornis and has shaped the genetic structure of this invasive ant worldwide.

Global genetic diversity, lineage distribution, and Wolbachia infection of the alfalfa weevil Hypera postica (Coleoptera: Curculionidae)

The alfalfa weevil (Hypera postica) is a well-known example of a worldwide-distributed pest with high genetic variation. Based on the mitochondrial genes, the alfalfa weevil clusters into two main mitochondrial lineages. However, there is no clear picture of the global diversity and distribution of these lineages; neither the drivers of its diversification are known. However, it appears likely that historic demographic events including founder effects played a role. In addition, Wolbachia, a widespread intracellular parasite/symbiont, likely played an important role in the evolution of the species. Wolbachia infection so far was only detected in the Western lineage of H. postica with no information on the infecting strain, its frequency, and its consequences on the genetic diversity of the host. We here used a combination of mitochondrial and nuclear sequences of the host and sequence information on Wolbachia to document the distribution of strains and the degree of infection. The Eastern lineage has a higher genetic diversity and is found in the Mediterranean, the Middle East, Eastern Europe, and eastern America, whereas the less diverse Western lineage is found in Central Europe and the western America. Both lineages are infected with the same common strain of Wolbachia belonging to Supergroup B. Based on neutrality tests, selection tests, and the current distribution and diversification of Wolbachia in H. postica, we suggested the Wolbachia infection did not shape genetic diversity of the host. The introduced populations in the United States are generally genetically less diverse, which is in line with founder effects.

Dynamics of an Ongoing Wolbachia Spread in the European Cherry Fruit Fly, Rhagoletis cerasi (Diptera: Tephritidae)

Numerous terrestrial arthropods are infected with the alphaproteobacterium Wolbachia. This endosymbiont is usually transmitted vertically from infected females to their offspring and can alter the reproduction of hosts through various manipulations, like cytoplasmic incompatibility (CI), enhancing its spread in new host populations. Studies on the spatial and temporal dynamics of Wolbachia under natural conditions are scarce. Here, we analyzed Wolbachia infection frequencies in populations of the European cherry fruit fly, Rhagoletis cerasi (L.), in central Germany—an area of an ongoing spread of the CI-inducing strain wCer2. In total, 295 individuals from 19 populations were PCR-screened for the presence of wCer2 and their mitochondrial haplotype. Results were compared with historic data to understand the infection dynamics of the ongoing wCer2 invasion. An overall wCer2 infection frequency of about 30% was found, ranging from 0% to 100% per population. In contrast to an expected smooth transition from wCer2-infected to completely wCer2-uninfected populations, a relatively scattered infection pattern across geography was observed. Moreover, a strong Wolbachia-haplotype association was detected, with only a few rare misassociations. Our results show a complex dynamic of an ongoing Wolbachia spread in natural field populations of R. cerasi.

Game of clones: Is Wolbachia inducing speciation in a weevil with a mixed reproductive mode?

Parthenogenesis is widely distributed in Metazoa but it is especially frequent in weevils (Coleoptera, Curculionidae) with one fifth of all known cases. Previous studies have shown that in the tribe Naupactini parthenogenetic reproduction most likely originated with an infection of the endoparasitic bacterium Wolbachia pipientis. In particular, Pantomorus postfasciatus possess a mixed reproductive mode: some populations have males while in others they are absent, and females produce clones by thelytoky. To better understand this scenario, we studied the population structure and infection status in 64 individuals of P. postfasciatus from Argentina and Brazil. We sequenced two mitochondrial (COI and COII) and one nuclear (ITS-1) fragments and obtained two very divergent haplogroups, one corresponding to the sexual populations uninfected with Wolbachia, and another conforming a monophyletic parthenogenetic (or presumptively parthenogenetic) and infected clade. Each of these haplogroups was identified as an independently evolutionary unit by all species delimitation analyses accomplished: multilocus *BEAST and BP&P, and single locus GMYC and K/θ rule. Additionally, present evidence suggests that Wolbachia infection occurred at least twice in all-female populations of P. postfasciatus with two different bacterial strains. Speciation mediated by Wolbachia is a recently described phenomenon and the case of P. postfasciatus is the first known case in a diplo-diploid insect. A model that describes how thelytoky-inducing phenotypes of Wolbachia could generate new lineages is discussed.

Influence of Wolbachia infection on mitochondrial DNA variation in the genus Polytremis (Lepidoptera: Hesperiidae)

The maternally inherited obligate bacteria Wolbachia is known for infecting the reproductive tissues of a wide range of arthropods and can contribute to phylogenetically discordant patterns between mtDNA and nDNA. In this study, we tested for an association between mito-nuclear discordance in Polytremis and Wolbachia infection. Six of the 17 species of Polytremis were found to be infected with Wolbachia. Overall, 34% (70/204) of Polytremis specimens were Wolbachia positive and three strains of Wolbachia identified using a wsp marker were further characterized as six strains based on MLST markers. Wolbachia acquisition in Polytremis appears to occur mainly through horizontal transmission rather than codivergence based on comparison of the divergence times of Wolbachia and Polytremis species. At the intraspecific level, one of the Wolbachia infections (wNas1) is associated with reduced mtDNA polymorphism in the infected Polytremis population. At the interspecific level, there is one case of mito-nuclear discordance likely caused by introgression of P. fukia mtDNA into P. nascens driven by another Wolbachia strain (wNas3). Based on an absence of infected males, we suspect that one Wolbachia strain (wNas2) affects sex ratio, but the phenotypic effects of the other strains are unclear. These data reveal a dynamic interaction between Polytremis and Wolbachia endosymbionts affecting patterns of mtDNA variation.

Evolutionary Mechanisms of Varying Chromosome Numbers in the Radiation of Erebia Butterflies

The evolution of intrinsic barriers to gene flow is a crucial step in the process of speciation. Chromosomal changes caused by fusion and fission events are one such barrier and are common in several groups of Lepidoptera. However, it remains unclear if and how chromosomal changes have contributed to speciation in this group. I tested for a phylogenetic signal of varying chromosome numbers in Erebia butterflies by combining existing sequence data with karyological information. I also compared different models of trait evolution in order to infer the underlying evolutionary mechanisms. Overall, I found significant phylogenetic signals that are consistent with non-neutral trait evolution only when parts of the mitochondrial genome were included, suggesting cytonuclear discordances. The adaptive evolutionary model tested in this study consistently outperformed the neutral model of trait evolution. Taken together, these results suggest that, unlike other Lepidoptera groups, changes in chromosome numbers may have played a role in the diversification of Erebia butterflies.

Wolbachia-induced meiotic drive and feminization is associated with an independent occurrence of selective mitochondrial sweep in a butterfly

Maternally inherited Wolbachia endosymbionts manipulate arthropod reproduction in various ways. In the butterfly Eurema mandarina, a cytoplasmic incompatibility-inducing Wolbachia strain wCI and the associated mtDNA haplotypes are known to originate from the sister species Eurema hecabe, which offered a good case study for microbe-mediated hybrid introgression. Besides wCI, some females with the Z0 karyotype harbour a distinct Wolbachia strain wFem, which causes all-female production by meiotic drive and feminization. We report that a considerable proportion of E. mandarina females (65.7%) were infected with both wCI and wFem (CF) on Tanegashima Island. While females singly infected with wCI (C) produced offspring at a 1 : 1 sex ratio, CF females produced only females. Although Z-linked sequence polymorphism showed no signs of divergence between C and CF females, mtDNA split into two discrete clades; one consisted of C females and the other CF females, both of which formed a clade with E. hecabe but not with uninfected E. mandarina This suggests that CF matrilines also, but independently, experienced a selective sweep after hybrid introgression from E. hecabe Distinct evolutionary forces were suggested to have caused C and CF matrilines to diverge, which would be irreversible because of the particular phenotype of wFem.

Feminizing Wolbachia endosymbiont disrupts maternal sex chromosome inheritance in a butterfly species

Wolbachia is a maternally inherited ubiquitous endosymbiotic bacterium of arthropods that displays a diverse repertoire of host reproductive manipulations. For the first time, we demonstrate that Wolbachia manipulates sex chromosome inheritance in a sexually reproducing insect. Eurema mandarina butterfly females on Tanegashima Island, Japan, are infected with the wFem Wolbachia strain and produce all‐female offspring, while antibiotic treatment results in male offspring. Fluorescence in situ hybridization (FISH) revealed that wFem‐positive and wFem‐negative females have Z0 and WZ sex chromosome sets, respectively, demonstrating the predicted absence of the W chromosome in wFem‐infected lineages. Genomic quantitative polymerase chain reaction (qPCR) analysis showed that wFem‐positive females lay only Z0 eggs that carry a paternal Z, whereas females from lineages that are naturally wFem‐negative lay both WZ and ZZ eggs. In contrast, antibiotic treatment of adult wFem females resulted in the production of Z0 and ZZ eggs, suggesting that this Wolbachia strain can disrupt the maternal inheritance of Z chromosomes. Moreover, most male offspring produced by antibiotic‐treated wFem females had a ZZ karyotype, implying reduced survival of Z0 individuals in the absence of feminizing effects of Wolbachia. Antibiotic treatment of wFem‐infected larvae induced male‐specific splicing of the doublesex (dsx) gene transcript, causing an intersex phenotype. Thus, the absence of the female‐determining W chromosome in Z0 individuals is functionally compensated by Wolbachia‐mediated conversion of sex determination. We discuss how Wolbachia may manipulate the host chromosome inheritance and that Wolbachia may have acquired this coordinated dual mode of reproductive manipulation first by the evolution of female‐determining function and then cytoplasmically induced disruption of sex chromosome inheritance.

Biogeography and systematics of endemic island damselflies: The Nesobasis and Melanesobasis (Odonata: Zygoptera) of Fiji

The study of island fauna has greatly informed our understanding of the evolution of diversity. We here examine the phylogenetics, biogeography, and diversification of the damselfly genera Nesobasis and Melanesobasis, endemic to the Fiji Islands, to explore mechanisms of speciation in these highly speciose groups. Using mitochondrial (COI, 12S) and nuclear (ITS) replicons, we recovered garli-part maximum likelihood and mrbayes Bayesian phylogenetic hypotheses for 26 species of Nesobasis and eight species/subspecies of Melanesobasis. Biogeographical patterns were explored using lagrange and bayes-lagrange and interpreted through beast relaxed clock dating analyses. We found that Nesobasis and Melanesobasis have radiated throughout Fiji, but are not sister groups. For Nesobasis, while the two largest islands of the archipelago-Viti Levu and Vanua Levu-currently host two distinct species assemblages, they do not represent phylogenetic clades; of the three major groupings each contains some Viti Levu and some Vanua Levu species, suggesting independent colonization events across the archipelago. Our beast analysis suggests a high level of species diversification around 2-6 Ma. Our ancestral area reconstruction (rasp-lagrange) suggests that both dispersal and vicariance events contributed to the evolution of diversity. We thus conclude that the evolutionary history of Nesobasis and Melanesobasis is complex; while inter-island dispersal followed by speciation (i.e., peripatry) has contributed to diversity, speciation within islands appears to have taken place a number of times as well. This speciation has taken place relatively recently and appears to be driven more by reproductive isolation than by ecological differentiation: while species in Nesobasis are morphologically distinct from one another, they are ecologically very similar, and currently are found to exist sympatrically throughout the islands on which they are distributed. We consider the potential for allopatric speciation within islands, as well as the influence of parasitic endosymbionts, to explain the high rates of speciation in these damselflies.

The hitchhiker's guide to Europe: the infection dynamics of an ongoing Wolbachia invasion and mitochondrial selective sweep in Rhagoletis cerasi

Wolbachia is a maternally inherited and ubiquitous endosymbiont of insects. It can hijack host reproduction by manipulations such as cytoplasmic incompatibility (CI) to enhance vertical transmission. Horizontal transmission of Wolbachia can also result in the colonization of new mitochondrial lineages. In this study, we present a 15‐year‐long survey of Wolbachia in the cherry fruit fly Rhagoletis cerasi across Europe and the spatiotemporal distribution of two prevalent strains, wCer1 and wCer2, and associated mitochondrial haplotypes in Germany. Across most of Europe, populations consisted of either 100% singly (wCer1) infected individuals with haplotype HT1, or 100% doubly (wCer1&2) infected individuals with haplotype HT2, differentiated only by a single nucleotide polymorphism. In central Germany, singly infected populations were surrounded by transitional populations, consisting of both singly and doubly infected individuals, sandwiched between populations fixed for wCer1&2. Populations with fixed infection status showed perfect association of infection and mitochondria, suggesting a recent CI‐driven selective sweep of wCer2 linked with HT2. Spatial analysis revealed a range expansion for wCer2 and a large transition zone in which wCer2 splashes appeared to coalesce into doubly infected populations. Unexpectedly, the transition zone contained a large proportion (22%) of wCer1&2 individuals with HT1, suggesting frequent intraspecific horizontal transmission. However, this horizontal transmission did not break the strict association between infection types and haplotypes in populations outside the transition zone, suggesting that this horizontally acquired Wolbachia infection may be transient. Our study provides new insights into the rarely studied Wolbachia invasion dynamics in field populations.

Taxonomic status and molecular phylogeography of two sibling species of Polytremis (Lepidoptera: Hesperiidae)

The skipper Polytremis theca species complex is widely distributed in the south of the Qinling Mountains in China. A recent study of the Polytremis genus suggested that this species might encompass two differentiated lineages. We tested this hypothesis, by carrying out a phylogenetic study of this agricultural pest based on nationwide sampling and the evaluation of mitochondrial and nuclear DNA markers. We show that this species is actually an amalgamation of two sibling taxa (P. t. theca and P. t. fukia), which displayed levels of genetic divergence as great as those generally found between sister species in the Polytremis genus, suggesting that they actually correspond to two distinct species. The Divergence time estimates suggest that an active period of speciation within Polytremis occurred within the Pleistocene eras. Based on its distinct phylogenetic placement and geographical isolation, we suggest that the subspecies should be elevated to full species status under the phylogenetic species concept, which has significant management implications.

Reassessment of Species Diversity of the Subfamily Denticollinae (Coleoptera: Elateridae) through DNA Barcoding

The subfamily Denticollinae is a taxonomically diverse group in the family Elateridae. Denticollinae includes many morphologically similar species and crop pests, as well as many undescribed species at each local fauna. To construct a rapid and reliable identification system for this subfamily, the effectiveness of molecular species identification was assessed based on 421 cytochrome c oxidase subunit I (COI) sequences of 84 morphologically identified species. Among the 84 morphospecies, molecular species identification of 60 species (71.4%) was consistent with their morphological identifications. Six cryptic and/or pseudocryptic species with large genetic divergence (>5%) were confirmed by their sympatric or allopatric distributions. However, 18 species, including a subspecies, had ambiguous genetic distances and shared overlapping intra- and interspecific genetic distances (range: 2.12%–3.67%) suggesting incomplete lineage sorting, introgression of mitochondrial genome, or affection by endosymbionts, such as Wolbachia infection, between species and simple genetic variation within species. In this study, we propose a conservative threshold of 3.6% for convenient molecular operational taxonomic unit (MOTU) identification in the subfamily Denticollinae based on the results of pairwise genetic distances analyses using neighbor-joining, mothur, Automatic Barcode Gap Discovery analysis, and tree-based species delimitation by Poisson Tree Processes analysis. Using the 3.6% threshold, we identified 87 MOTUs and found 8 MOTUs in the interval between 2.5% to 3.5%. Evaluation of MOTUs identified in this range requires integrative species delimitation, including review of morphological and ecological differences as well as sensitive genetic markers. From this study, we confirmed that COI sequence is useful for reassessing species diversity for polymorphic and polytypic species occurring in sympatric and allopatric distributions, and for a single species having an extensively large habitat.

Multiple endosymbiont infections and reproductive manipulations in a linyphiid spider population

In many arthropods, maternally inherited endosymbiotic bacteria can increase infection frequency by manipulating host reproduction. Multiple infections of different bacteria in a single host population are common, yet few studies have documented concurrent endosymbiont phenotypes or explored their potential interactions. We hypothesized that spiders might be a particularly useful taxon for investigating endosymbiont interactions, because they are host to a plethora of endosymbiotic bacteria and frequently exhibit multiple infections. We established two matrilines from the same population of the linyphiid spider Mermessus fradeorum and then used antibiotic curing and controlled mating assays to demonstrate that each matriline was subject to a distinct endosymbiotic reproductive manipulation. One matriline was co-infected with Rickettsia and Wolbachia and produced offspring with a radical female bias. Antibiotic treatment eliminated both endosymbionts and restored an even sex ratio to subsequent generations. Chromosomal and fecundity observations suggest a feminization mechanism. In the other matriline, a separate factorial mating assay of cured and infected spiders demonstrated strong cytoplasmic incompatibility (CI) induced by a different strain of Wolbachia. However, males with this Wolbachia induced only mild CI when mated with the Rickettsia-Wolbachia females. In a subsequent survey of a field population of M. fradeorum, we detected these same three endosymbionts infecting 55% of the spiders in almost all possible combinations, with nearly half of the infected spiders exhibiting multiple infection. Our results suggest that a dynamic network of endosymbionts may interact both within multiply infected hosts and within a population subject to multiple strong reproductive manipulations.

Population genetic structure and Wolbachia infection in an endangered butterfly, Zizina emelina (Lepidoptera, Lycaenidae), in Japan

Zizina emelina (de l'Orza) is listed on Japan's Red Data List as an endangered species because of loss of its principal food plant and habitat. We compared parts of the mitochondrial and nuclear genes of this species to investigate the level of genetic differentiation among the 14 extant populations. We also examined infection of the butterfly with the bacterium Wolbachia to clarify the bacterium's effects on the host population's genetic structure. Mitochondrial and nuclear DNA analyses revealed that haplotype composition differed significantly among most of the populations, and the fixation index F ST was positively correlated with geographic distance. In addition, we found three strains of Wolbachia, one of which was a male killer; these strains were prevalent in several populations. There was linkage between some host mitochondrial haplotypes and the three Wolbachia strains, although no significant differences were found in a comparison of host mitochondrial genetic diversity with nuclear genetic diversity in Wolbachia-infected or -uninfected populations. These genetic analyses and Wolbachia infection findings show that Z. emelina has little migratory activity and that little gene flow occurs among the current populations.

Modulation of host immunity and reproduction by horizontally acquired Wolbachia

The Wolbachia are symbiotic bacteria vertically transmitted from one host generation to another. However, a growing amount of data shows that horizontal transfers of Wolbachia also frequently occur within and between host species. The consequences of the arrival of new symbionts on host physiology can be studied by their experimental introduction in asymbiotic hosts. After experimental transfers of the eight major isopod Wolbachia strains in the isopod Porcellio dilatatus only two of them (wCon and wDil) were found to (1) have no pathogenic effect on the host and (2) be able to pass vertically to the host offspring. In the present work, we studied the influence of these two strains, able to complete an horizontal transfer, on immunity and reproduction of P. dilatatus at two stages of the transfer: (1) in recipient hosts that encounter the symbionts: to test the influence of symbiont when acquired during host life and (2) in vertically infected offspring: to test the influence of a symbiotic interaction occurring all lifelong. The impact of Wolbachia varied depending on the stage: there were clearer effects in vertically infected individuals than in those that acquired the symbionts during their lives. Moreover, the two Wolbachia strains showed contrasted effects: the strain wCon tended to reduce the reproductive investment but to maintain or increase immune parameters whilst wDil had positive effects on reproductive investment but decreased the investment in some immune parameters. These results suggest that horizontally acquisition of Wolbachia can influence the balance between host immune and reproductive traits.

Wolbachia infection status and genetic structure in natural populations of Polytremis nascens (Lepidoptera: Hesperiidae)

The maternally inherited obligate bacteria Wolbachia is known for infecting the reproductive tissues of a wide range of arthropods. In this study, we surveyed Wolbachia infections in Polytremis nascens (Lepidoptera: Hesperiidae) from 14 locations in China by amplifying the 16S rRNA gene with a nested PCR method and revealed the effect of Wolbachia on host mitochondrial DNA. The results show that 31% (21/67) are Wolbachia positive among all specimens and mainly prevails in southern populations in China. No significant difference in the prevalence is found between the sexes. Notably, the nucleotide diversity of Wolbachia infected butterflies is smaller compared to that of uninfected butterflies. The mitochondrial DNA of infected group appear to be not evolving neutrally (Tajima's D value=-2.3303 and Fu's F values=-3.7068). The analysis of molecular variance shows significant differentiation of mitochondrial haplotypes between infected and uninfected specimens (FST=0.6064). The mismatch analysis speculated the different expansion pattern in Wolbachia infected specimens and all P. nascens specimens. These results suggest that the populations of P. nascens may have recently been subjected to a Wolbachia-induced sweep. Additionally, phylogenetic analysis differentiated the mitochondrial haplotypes of P. nascens into three major clades. The clades are in perfect agreement with the pattern of Wolbachia infection. One of the clades grouped with the butterflies infected with Wolbachia. The remaining two clades grouped with uninfected butterflies from the central-west of China populations and Eastern and Southern China populations respectively, which are isolated mainly by the Yangtze River. The analysis of haplotype networks, geographic distribution and population size change shows that Haplotype 1 in central-west of China is the ancestral haplotype and the populations of P. nascens are expanded.

Gene-flow in a mosaic hybrid zone: is local introgression adaptive?

Genome-wide scans of genetic differentiation between hybridizing taxa can identify genome regions with unusual rates of introgression. Regions of high differentiation might represent barriers to gene flow, while regions of low differentiation might indicate adaptive introgression-the spread of selectively beneficial alleles between reproductively isolated genetic backgrounds. Here we conduct a scan for unusual patterns of differentiation in a mosaic hybrid zone between two mussel species, Mytilus edulis and M. galloprovincialis. One outlying locus, mac-1, showed a characteristic footprint of local introgression, with abnormally high frequency of edulis-derived alleles in a patch of M. galloprovincialis enclosed within the mosaic zone, but low frequencies outside of the zone. Further analysis of DNA sequences showed that almost all of the edulis allelic diversity had introgressed into the M. galloprovincialis background in this patch. We then used a variety of approaches to test the hypothesis that there had been adaptive introgression at mac-1. Simulations and model fitting with maximum-likelihood and approximate Bayesian computation approaches suggested that adaptive introgression could generate a "soft sweep," which was qualitatively consistent with our data. Although the migration rate required was high, it was compatible with the functioning of an effective barrier to gene flow as revealed by demographic inferences. As such, adaptive introgression could explain both the reduced intraspecific differentiation around mac-1 and the high diversity of introgressed alleles, although a localized change in barrier strength may also be invoked. Together, our results emphasize the need to account for the complex history of secondary contacts in interpreting outlier loci.

Wolbachia effects in natural populations of Chorthippus parallelus from the Pyrenean hybrid zone

We evaluate for the first time the effect of Wolbachia infection, involving two different supergroups, on the structure and dynamics of the hybrid zone between two subspecies of Chorthippus parallelus (Orthoptera) in the Pyrenees. Wolbachia infection showed no effects on female fecundity or a slight increment in females infected by F supergroup, although in the last case it has to be well established. Cytoplasmic incompatibility (CI) is confirmed in crosses carried out in the field between individuals from a natural hybrid population. This CI, registered as the relative reduction in embryo production (sh ), was of sh = 0.355 and sh = 0.286 in unidirectional crosses involving B and F supergroups, respectively. CI also occurred in bidirectional crosses (sh = 0.147) but with a weaker intensity. The transmission rates of the two Wolbachia strains (B and F) were estimated by the optimization of a theoretical model to reach the infection frequencies observed in certain population. To fit this scenario, both supergroups should present transmission rates close to 1. Further, we have simulated the infection dynamics, and hence, the capacity of Wolbachia to structure the population of the host insects and to affect to reproduction and genetic introgression in the hybrid zone. This represents a first example of the influence of Wolbachia in an insect natural hybrid zone.

Draft Genome Sequence of Commensalibacter papalotli MX01, a Symbiont Identified from the Guts of Overwintering Monarch Butterflies

We report the draft genome sequence of Commensalibacter papalotli strain MX01, isolated from the intestines of an overwintering monarch butterfly. The 2,332,652-bp AT-biased genome of C. papalotli MX01 is the smallest genome for a member of the Acetobacteraceae family and provides the first evidence of plasmids in Commensalibacter.

Metamorphosis of a butterfly-associated bacterial community

Butterflies are charismatic insects that have long been a focus of biological research. They are also habitats for microorganisms, yet these microbial symbionts are little-studied, despite their likely importance to butterfly ecology and evolution. In particular, the diversity and composition of the microbial communities inhabiting adult butterflies remain uncharacterized, and it is unknown how the larval (caterpillar) and adult microbiota compare. To address these knowledge gaps, we used Illumina sequencing of 16S rRNA genes from internal bacterial communities associated with multiple life stages of the neotropical butterfly Heliconius erato. We found that the leaf-chewing larvae and nectar- and pollen-feeding adults of H. erato contain markedly distinct bacterial communities, a pattern presumably rooted in their distinct diets. Larvae and adult butterflies host relatively small and similar numbers of bacterial phylotypes, but few are common to both stages. The larval microbiota clearly simplifies and reorganizes during metamorphosis; thus, structural changes in a butterfly's bacterial community parallel those in its own morphology. We furthermore identify specific bacterial taxa that may mediate larval and adult feeding biology in Heliconius and other butterflies. Although male and female Heliconius adults differ in reproductive physiology and degree of pollen feeding, bacterial communities associated with H. erato are not sexually dimorphic. Lastly, we show that captive and wild individuals host different microbiota, a finding that may have important implications for the relevance of experimental studies using captive butterflies.

Wolbachia infections mimic cryptic speciation in two parasitic butterfly species, Phengaris teleius and P. nausithous (Lepidoptera: Lycaenidae)

Deep mitochondrial divergence within species may result from cryptic speciation, from phylogeographic isolation or from endosymbiotic bacteria like Wolbachia that manipulate host reproduction. Phengaris butterflies are social parasites that spend most of their life in close relationship with ants. Previously, cryptic speciation has been hypothesised for two Phengaris species based on divergent mtDNA sequences. Since Phengaris species are highly endangered, the existence of cryptic species would have drastic consequences for conservation and management. We tested for cryptic speciation and alternative scenarios in P. teleius and P. nausithous based on a comprehensive sample across their Palaearctic ranges using COI gene sequences, nuclear microsatellites and tests for Wolbachia. In both species a deep mitochondrial split occurring 0.65-1.97 myrs ago was observed that did not correspond with microsatellite data but was concordant with Wolbachia infection. Haplotypes previously attributed to cryptic species were part of the Wolbachia-infected clades. In both species remaining phylogeographic structure was largely consistent between mitochondrial and nuclear genomes. In P. teleius several mitochondrial and nuclear groups were observed in East Asia while a single haplogroup and nuclear cluster prevailed across continental Eurasia. Neutrality tests suggested rapid demographic expansion into that area. In contrast, P. nausithous had several mitochondrial and nuclear groups in Europe, suggesting a complex phylogeographic history in the western part of the species range. We conclude that deep intraspecific divergences found in DNA barcode studies do not necessarily need to represent cryptic speciation but instead can be due to both infection by Wolbachia and phylogeographic structure.

Speciation in the asexual realm: is the parthenogenetic weevil Naupactus cervinus a complex of species in statu nascendi?

Population genetic theory shows that asexual organisms may evolve into species, which behave as independent evolutionary units. As a result, they form genotypic clusters separated by deep gaps due to geographic isolation and/or divergent selection. Identification of several genetically divergent groups of weevils embodied in the nominal species Naupactus cervinus deserves further study, in order to test if these lineages are evolving independently. In the present paper we tested if the parthenogenetic weevil N. cervinus, native to South America and broadly distributed throughout the world, contains more than one evolutionary unit. For this purpose, we applied three different approaches, a multilocus phylogenetic analysis, the GMYC approach and the K/θ method. We accomplished these analyses through a survey of mitochondrial (COI and COII genes) and nuclear (ITS1 sequence) genetic variation and morphometric analysis in a sample which included individuals from different locations within the native geographic range of N. cervinus. In addition, we compared the divergence accumulated in this species with that in another weevil of the same tribe (Naupactini) showing identical reproductive mode to see if similar levels of morphological variation matches similar levels of genetic divergence. We report the presence of two independent evolutionary units living in sympatry in forest areas. The incongruence between mitochondrial and nuclear datasets analyzed herein reflects incomplete lineage sorting of the nuclear marker and different evolutionary rates between genomes. Ecological divergence driven by natural selection (sympatry) or secondary contact after geographic isolation (allopatry) might explain the deep gaps in mitochondrial phylogenies. Instead, Wolbachia infection was ruled out as a causal factor for such differentiation. We conclude that N. cervinus is probably a species complex with at least two well differentiated lineages that would represent a cluster of species in statu nascendi.

Linkage map of the peppered moth, Biston betularia (Lepidoptera, Geometridae): a model of industrial melanism

We have constructed a linkage map for the peppered moth (Biston betularia), the classical ecological genetics model of industrial melanism, aimed both at localizing the network of loci controlling melanism and making inferences about chromosome dynamics. The linkage map, which is based primarily on amplified fragment length polymorphisms (AFLPs) and genes, consists of 31 linkage groups (LGs; consistent with the karyotype). Comparison with the evolutionarily distant Bombyx mori suggests that the gene content of chromosomes is highly conserved. Gene order is conserved on the autosomes, but noticeably less so on the Z chromosome, as confirmed by physical mapping using bacterial artificial chromosome fluorescence in situ hybridization (BAC-FISH). Synteny mapping identified three pairs of B. betularia LGs (11/29, 23/30 and 24/31) as being orthologous to three B. mori chromosomes (11, 23 and 24, respectively). A similar finding in an outgroup moth (Plutella xylostella) indicates that the B. mori karyotype (n=28) is a phylogenetically derived state resulting from three chromosome fusions. As with other Lepidoptera, the B. betularia W chromosome consists largely of repetitive sequence, but exceptionally we found a W homolog of a Z-linked gene (laminin A), possibly resulting from ectopic recombination between the sex chromosomes. The B. betularia linkage map, featuring the network of known melanization genes, serves as a resource for melanism research in Lepidoptera. Moreover, its close resemblance to the ancestral lepidopteran karyotype (n=31) makes it a useful reference point for reconstructing chromosome dynamic events and ancestral genome architectures. Our study highlights the unusual evolutionary stability of lepidopteran autosomes; in contrast, higher rates of intrachromosomal rearrangements support a special role of the Z chromosome in adaptive evolution and speciation.

Population genetics of two asexually and sexually reproducing psocids species inferred by the analysis of mitochondrial and nuclear DNA sequences

Background

The psocids Liposcelis bostrychophila and L. entomophila (Psocoptera: Liposcelididae) are found throughout the world and are often associated with humans, food stores and habitations. These insects have developed high levels of resistance to various insecticides in grain storage systems. However, the population genetic structure and gene flow of psocids has not been well categorized, which is helpful to plan appropriate strategies for the control of these pests.

Methodology/Principal Findings

The two species were sampled from 15 localities in China and analyzed for polymorphisms at the mitochondrial DNA (Cytb) and ITS (ITS1-5.8S-ITS2) regions. In total, 177 individual L. bostrychophila and 272 individual L. entomophila were analysed. Both Cytb and ITS sequences showed high genetic diversity for the two species with haplotype diversities ranged from 0.154±0.126 to 1.000±0.045, and significant population differentiation (mean F_ST = 0.358 for L. bostrychophila; mean F_ST = 0.336 for L. entomophila) was also detected among populations investigated. A Mantel test indicated that for both species there was no evidence for isolation-by-distance (IBD). The neutrality test and mismatch distribution statistics revealed that the two species might have undergone population expansions in the past.

Conclusion

Both L. bostrychophila and L. entomophila displayed high genetic diversity and widespread population genetic differentiation within and between populations. The significant population differentiation detected for both psocids may be mainly due to other factors, such as genetic drift, inbreeding or control practices, and less by geographic distance since an IBD effect was not found.