The maternally inherited intracellular bacteria Wolbachia can manipulate host reproduction in various ways that foster frequency increases within and among host populations. Manipulations involving cytoplasmic incompatibility (CI), where matings between infected males and uninfected females produce non-viable embryos, are common in arthropods and produce a reproductive advantage for infected females. CI was associated with the spread of Wolbachia variant wRi in Californian populations of Drosophila simulans, which was interpreted as a bistable wave, in which local infection frequencies tend to increase only once the infection becomes sufficiently common to offset imperfect maternal transmission and infection costs. However, maternally inherited Wolbachia are expected to evolve towards mutualism, and they are known to increase host fitness by protecting against infectious microbes or increasing fecundity. We describe the sequential spread over approximately 20 years in natural populations of D. simulans on the east coast of Australia of two Wolbachia variants (wAu and wRi), only one of which causes significant CI, with wRi displacing wAu since 2004. Wolbachia and mtDNA frequency data and analyses suggest that these dynamics, as well as the earlier spread in California, are best understood as Fisherian waves of favourable variants, in which local spread tends to occur from arbitrarily low frequencies. We discuss implications for Wolbachia-host dynamics and coevolution and for applications of Wolbachia to disease control.

Background

The interaction between the Blue Moon butterfly, Hypolimnas bolina, and Wolbachia has attracted interest because of the high prevalence of male-killing achieved within the species, the ecological consequences of this high prevalence, the intensity of selection on the host to suppress the infection, and the presence of multiple Wolbachia infections inducing different phenotypes. We examined diversity in the co-inherited marker, mtDNA, and the partitioning of this between individuals of different infection status, as a means to investigate the population biology and evolutionary history of the Wolbachia infections.

Results

Part of the mitochondrial COI gene was sequenced from 298 individuals of known infection status revealing ten different haplotypes. Despite very strong biological evidence that the sample represents a single species, the ten haplotypes did not fall within a monophyletic clade within the Hypolimnas genus, with one haplotype differing by 5% from the other nine. There were strong associations between infection status and mtDNA haplotype. The presence of wBol1 infection in association with strongly divergent haplotypes prompted closer examination of wBol1 genetic variation. This revealed the existence of two cryptic subtypes, wBol1a and wBol1b. The wBol1a infection, by far the most common, was in strict association with the single divergent mtDNA haplotype. The wBol1b infection was found with two haplotypes that were also observed in uninfected specimens. Finally, the wBol2 infection was associated with a large diversity of mtDNA haplotypes, most often shared with uninfected sympatric butterflies.

Conclusion

This data overall supports the hypothesis that high prevalence of male-killing Wolbachia (wBol1) in H. bolina is associated with very high transmission efficiency rather than regular horizontal transmission. It also suggests this infection has undergone a recent selective sweep and was introduced in this species through introgression. In contrast, the sharing of haplotypes between wBol2-infected and uninfected individuals indicates that this strain is not perfectly transmitted and/or shows a significant level of horizontal transmission.

Male-killing bacteria are widespread in arthropods, and can profoundly alter the reproductive biology of their host species. Here we detail the first case of complete suppression of a male killer. The nymphalid butterfly Hypolimnas bolina is infected with a strain of the bacterium Wolbachia, wBol1, which kills male host embryos in Polynesian populations, but does not do so in many areas of Southeast Asia, where both males and female adults are naturally infected, and wBol1-infected females produce a 1:1 sex ratio. We demonstrate that absence of male killing by wBol1 is associated with dominant zygotic suppression of the action of the male killer. Simulations demonstrate host suppressors of male-killer action can spread very rapidly, and historical data indicating the presence of male killing in Southeast Asia in the very recent past suggests suppressor spread has been a very recent occurrence. Thus, male killer/host interactions are much more dynamic than previously recognised, with rapid and dramatic loss of the phenotype. Our results also indicate that suppression can render male killers completely quiescent, leading to the conclusion that some species that do not currently express a male killer may have done so in the past, and thus that more species have had their biology affected by these parasites than previously believed.

New genetic diagnostic approaches have greatly aided efforts to document global biodiversity and improve biosecurity. This is especially true for organismal groups in which species diversity has been underestimated historically due to difficulties associated with sampling, the lack of clear morphological characteristics, and/or limited availability of taxonomic expertise. Among these methods, DNA sequence barcoding (also known as "DNA barcoding") and by extension, meta-barcoding for biological communities, has emerged as one of the most frequently utilized methods for DNA-based species identifications. Unfortunately, the use of DNA barcoding is limited by the availability of complete reference libraries (i.e., a collection of DNA sequences from morphologically identified species), and by the fact that the vast majority of species do not have sequences present in reference databases. Such conditions are critical especially in tropical locations that are simultaneously biodiversity rich and suffer from a lack of exploration and DNA characterization by trained taxonomic specialists. To facilitate efforts to document biodiversity in regions lacking complete reference libraries, we developed a novel statistical approach that categorizes unidentified species as being either likely native or likely nonnative based solely on measures of nucleotide diversity. We demonstrate the utility of this approach by categorizing a large sample of specimens of terrestrial insects and spiders (collected as part of the Moorea BioCode project) using a generalized linear mixed model (GLMM). Using a training data set of known endemic (n = 45) and known introduced species (n = 102), we then estimated the likely native/nonnative status for 4,663 specimens representing an estimated 1,288 species (412 identified species), including both those specimens that were either unidentified or whose endemic/introduced status was uncertain. Using this approach, we were able to increase the number of categorized specimens by a factor of 4.4 (from 794 to 3,497), and the number of categorized species by a factor of 4.8 from (147 to 707) at a rate much greater than chance (77.6% accuracy). The study identifies phylogenetic signatures of both native and nonnative species and suggests several practical applications for this approach including monitoring biodiversity and facilitating biosecurity.