Evidence for common horizontal transmission of Wolbachia among butterflies and moths
© Ahmed et al. 2016
Received: 20 October 2015
Accepted: 18 April 2016
Published: 27 May 2016
Wolbachia is one of the most widespread bacteria on Earth. Previous research on Wolbachia-host interactions indicates that the bacterium is typically transferred vertically, from mother to offspring, through the egg cytoplasm. Although horizontal transmission of Wolbachia from one species to another is reported to be common in arthropods, limited direct ecological evidence is available. In this study, we examine horizontal transmission of Wolbachia using a multilocus sequence typing (MLST) strains dataset and used Wolbachia and Lepidoptera genomes to search for evidence for lateral gene transfer (LGT) in Lepidoptera, one of the most diverse cosmopolitan insect orders. We constructed a phylogeny of arthropod-associated MLST Wolbachia strains and calibrated the age of Wolbachia strains associated with lepidopteran species.
Our results reveal inter-specific, inter-generic, inter-familial, and inter-ordinal horizontal transmission of Wolbachia strains, without discernible geographic patterns. We found at least seven probable cases of horizontal transmission among 31 species within Lepidoptera and between Lepidoptera and other arthropod hosts. The divergence time analysis revealed that Wolbachia is recently (22.6–4.7 mya, 95 % HPD) introduced in Lepidoptera. Analysis of nine Lepidoptera genomes (Bombyx mori, Danaus plexippus, Heliconius melpomene, Manduca sexta, Melitaea cinxia, Papilio glaucus, P. polytes, P. xuthus and Plutella xylostella) yielded one possible instance of Wolbachia LGT.
Our results provide evidence of high incidence of identical and multiple strains of Wolbachia among butterflies and moths, adding Lepidoptera to the growing body of evidence for common horizontal transmission of Wolbachia. This study demonstrates interesting dynamics of this remarkable and influential microorganism.
Offspring vertically inherit both nuclear and non-nuclear genetic material from their mothers . Among the non-nuclear material inherited are intracellular bacteria which are transferred vertically from mother to offspring and often live in symbioses with their hosts . These symbionts may be obligate (essential for host survival) or facultative, in which case they can increase or decrease host fitness [3, 4]. Obligate symbionts are found within specialized cells and typically share a long evolutionary history with their hosts , whereas facultative symbionts tend to have more recently formed host associations. Wolbachia (Alphaproteobacteria: Rickettsiales: Rickettsiaceae) is a genus of facultative endosymbiont common among arthropods that is estimated to have infected more than half of arthropod species , including two-thirds of all extant insect species . As with other facultative endosymbionts, Wolbachia has been thought to primarily undergo vertical transmission from mother to offspring with high fidelity . However, symbionts can also develop host associations via horizontal transmission between different host species [2, 4, 8]. Horizontal transmission is thought to be the most likely explanation for closely related symbionts occurring in phylogenetically distant insect lineages [2, 8–13]. There have been multiple phylogenetic and transinfection studies reporting evidence of Wolbachia transmission between both phylogenetically close and distant hosts [9, 14–18]; it is therefore probable that horizontal transmission of Wolbachia is occurring between some arthropod taxa .
Butterflies and moths (Lepidoptera) constitute one of the most diverse insect orders with nearly 158,000 described species . Lepidoptera play an important role in ecosystems and serve primarily as pollinators and herbivores, though some species feed on blood and other animal secretions [20–23]. The order includes many significant agricultural pests, and some species serve as models for many biological disciplines . Furthermore, lepidopteran larvae are hosts to other major insect radiations – the parasitic flies and wasps [25–27]. Despite the diversity of Lepidoptera and their many associations with other organisms, little is known about their bacterial community.
Wolbachia are some of the most widespread endosymbiotic microbes [6, 28–30]. In nematodes, Wolbachia interact mutually , and in arthropods, Wolbachia most commonly interact with their hosts via a parasitic manipulation of the reproductive system . Consequently, Wolbachia has been thought to undergo vertical transmission much more frequently than horizontal transmission . Wolbachia most commonly affects Lepidoptera via reproductive manipulation and can induce multiple phenotypes including feminization, male killing, and cytoplasmic incompatibility [31–33]. One strain of Wolbachia enhances the susceptibility of its lepidopteran host to baculovirus, rendering it a potential biological control agent against the agricultural pest Spodoptera exempta . It was recently estimated that approximately 80 % of Lepidoptera species are infected with Wolbachia , a prediction that is considerably higher than the 52 % estimated infection frequency across arthropods . However, the reported mean prevalence (27 %) in Lepidoptera  does not significantly differ from the estimated prevalence in arthropods (24 %) . The high incidence and low prevalence may reflect opportunities for substantial horizontal transfer of Wolbachia in Lepidoptera.
After Wolbachia undergoes stable horizontal transmission from natural to novel hosts, there are multiple possible phenotypic effects. We define “phenotype” as the set of observable characteristics of host result from its interaction with Wolbachia. The Wolbachia phenotype can become stronger, weaker, or remain the same, and in some cases, it can be changed to an unknown phenotype that is novel to the host . Additionally, once Wolbachia has successfully established a close relationship with its novel host, it may transfer a gene from its genome to the host genome over time . This is known as lateral gene transfer (LGT) [36, 37], and LGT is thought to be responsible for the presence of Wolbachia genes in 70 % of arthropod and nematode genomes [36, 38, 39]. A recent study showed evidence of ancient LGT of Enterococcus bacteria in Lepidoptera .
In this study, we 1) analyzed all published multilocus sequence typing strains (MLST) of Wolbachia including those from lepidopteran hosts in order to explore potential instances of horizontal transmission events, 2) surveyed transinfection experiments in Lepidoptera, to detail the factors underlying the host phenotype after horizontal transmission has occurred, and 3) searched for evidence of LGT between Wolbachia and Lepidoptera genomes. Our analyses reflect the complex dynamics of transmission between Wolbachia and their lepidopteran hosts.
We used multilocus sequence typing (MLST) strains based on five loci to identify and explore Wolbachia strain diversity. MLST provides a universal and unambiguous tool for strain typing, population genetics, and molecular evolutionary studies . MLST was developed as a universal genotyping tool for Wolbachia and was found effective for detecting diversity among strains within a single host species, as well as for identifying closely related strains found in different arthropod hosts . We downloaded and analyzed all 345 publically available strains of Wolbachia in arthropods and nematodes on March 31, 2014 from the PubMLST website (http://pubmlst.org/Wolbachia/) developed by Jolley and Maiden  (Additional file 1: Table S1). Approximately 26 % of these strains (90/345) were associated with lepidopteran hosts: 81 were strictly found in lepidopteran hosts whereas nine strains were found in both lepidopteran and non-lepidopteran arthropod hosts (Additional file 2: Table S2). Some of the strains from lepidopteran hosts (16/90) were unnamed and incomplete because not all five of the MLST loci were sequenced (gatB, coxA hcpA, ftsZ, and fbpA); these strains were designated as unassigned (UA) strains (Additional file 3: Table S3), and we included them in our analysis as such.
Sequence alignment and datasets
For ingroups, we included 345 MLST strains based on five MLST loci (gatB, coxA hcpA, ftsZ, and fbpA) of Wolbachia. For outgroups, we included bacteria closely related to Wolbachia: Anaplasma marginale (NCBI Genome accession no. NC_022760), Ehrlichia ruminantium (NC_006831) and Rickettsia slovaca (NC_017065), and extracted the five MLST loci from these genomes. These three outgroups and 345 ingroups were downloaded and aligned with the GINS-I algorithm in MAFFT . Geneious v8  was used to trim, align, and concatenated the five MLST loci. The best model and partitioning scheme were chosen using the Bayesian Information Criterion (BIC) in PartitionFinder v1.0.1  and resulted in two partitions (a combined first and second codon position [nt12]; and third codon positions only [nt3]).
Maximum likelihood (ML) phylogenetic analyses were conducted in RAxML v8  using a GTR + G model for each partition. To estimate the best ML tree in RAxML, we used the “–f a” option to estimate 1000 bootstraps and perform a likelihood search, as well as 200 “–f d” searches that started from a randomly generated parsimony tree, following the general methods of Kawahara et al. . We also estimated SH-like branch support  for the best topology in RAxML v8. We used the same method to construct a second ML tree for a smaller dataset of 51 strains found only in lepidopteran hosts, using three different outgroups: ID 37 from Supergroup D (host Brugia malayi, Nematoda), ID 505 from Supergroup C (host Onchocerca cervipedis, Nematoda) and ID 260 from Supergroup F (host Odontotermes horni, Isoptera).
A phylogeny of Wolbachia strains was also inferred with ClonalFrame v1.2  without outgroups. ClonalFrame uses information of substitution as well as recombination events and is therefore suitable to reconstruct bacterial evolution based on multilocus data . We performed ten separate runs, each with a burnin set to 250,000 generations and a sampling period of 750,000 generations, with a sampling frequency of 100. We chose the two runs with the highest mean log likelihood values and compared these to assess convergence of chains using the methods of Gelman and Rubin . Trees of the posterior samples of the converged runs were then combined to compute a majority rule consensus. We also calculated the ratio of nucleotides to point mutations (r/m).
Statistical parsimony network analysis has been shown to be useful for assessing species-level delimitation and to identify breaks in network connectivity [51–53]. Here we designated Wolbachia breaks in the network connectivity as identifying strains belonging to the Wolbachia species [54, 55]. In the present study, 90 strains were analyzed using a parsimony network approach  with TCS v.1.21  using a 95 % cut-off . The resulting networks identify both the relationships between the different haplotypes and the number of substitutions among connecting haplotypes .
A Mantel test was used to compute the Pearson correlation coefficient R using XLSTAT 2014 (http://www.xlstat.com). The test was performed on the pairwise node distance matrix of lepidopteran families from Regier et al.’s lepidopteran tree  and Wolbachia strains to test for significant association between matrices [60, 61].
Wolbachia strains from eight families of Lepidoptera were tested for codivergence.
We mapped the Wolbachia ClonalFrame tree onto the Lepidoptera phylogeny of Regier et al.  using JANE v4 . We reconstructed codivergence patterns with default cost values for cospeciation (0), duplication (1), duplication and host switch (2), loss (1), and failure to diverge (1). JANE analysis was performed using 500 generations and population sizes of 100. We selected an edge-based cost model and a node cost model; these models differ in counting the number events related to cospeciation, duplication and failure to diverge.
Divergence time estimation
To compare the age of Wolbachia divergence to previously published Lepidoptera divergence time estimations, we dated the splits of all Wolbachia strains found in lepidopteran species. Divergence time estimation analyses were performed in BEAST v2.1.3  and two independent calibrations were used to cross-validate our estimates . We applied the following calibration approaches: 1) using a recently published evolutionary rate of Wolbachia, estimated from Wolbachia genomes  and 2) using the age of a monophyletic set of strains shown to have strictly cospeciated with their hosts (bees) . We tested for the presence of a strict clock for nt12 and nt3 datasets using a likelihood ratio test (LRT)  in PAUP* v4.0 . Since the LRT test can be affected by recombination, we also used the relative-rate test (RRT) of Posada  in HyPHY , which can discriminate between strict and relaxed clock models in the presence of recombination. Because RRT requires that the outgroup taxa are recombination free, we used 3SEQ , implementing the full run mode for each gene to assure that the outgroup taxa did not have any recombinant genes. RRT analyses included taxa with unique sequences and no missing MLST loci and used two different outgroup MLST strains (13_Ekue_A_Ephestia_Pyralidae, 22_Aenc_B_Ugardan_Acraea_Nymphalidae). For the RRT, an alpha of ≤ 0.05 with a Bonferroni correction was treated as significant, and if any test was significant, then the strict clock is rejected . Since both the LRT and RRT rejected the strict clock, we estimated divergence times using a relaxed lognormal clock and applied one of the two calibrations to cross-validate estimates.
The first calibration scheme was based on the median rate (substitutions per site per generation) of the Wolbachia genome  reported in generations of Drosophila melanogaster, which was converted to year (10 generations per year) and scaled the rate to substitutions per site per million years (nt12 was 6.42× 10−3 [2.76 × 10−3 -1.29× 10−2, 95 % HPD] and nt3 was 6.87× 10−3 [2.88 × 10−3 -1.29× 10−2, 95 % HPD]). We set lognormal priors that spanned the 95 % HPD of the previous rate estimations (for nt12: lognormal M = 0.00642 and S = 0.45; for nt3: M = 0.00687 and S = 0.44). The second calibration scheme was based on the divergence time of MLST Wolbachia strains (wNLeu, wFla, wNPan) from Gerth et al. . The MRCA of these MLST strains is estimated at 1.7 mya (0.86–2.61, 95 % HPD) . We included these three strains in our divergence time analysis and calibrated the age of this group with a lognormal prior set to span the estimated HPD (M = 1.6 S = 0.33).
For each calibration scheme, we ran two BEAST analyses for a total of 4 runs using default settings for the remaining priors. We ran the MCMC chains for 150,000,000 generations, sampling every 1000th generation, and used Tracer  to ensure that the runs converged and had ESS values >200. For comparison with Wolbachia divergences, we applied the published divergence times of lepidopteran families [74, 75].
Evidence of LGT
MUMmer  was used to align Wolbachia and Lepidoptera genomes to search for evidence of LGT events. We used the following nine Wolbachia genomes: wBm (D) (host: Nematoda: Brugia malayi; AE017321) , wBol (B) (Lepidoptera: Nymphalidae: Hypolimnas bolina; CAOH01000001-CAOH0100014) , wMel (A) (Diptera: Drosophilidae: Drosophila melanogaster; NC_002978) , wPip (B) (Diptera: Culicidae: Culex quinquefasciatus; NC_010981) , wRi (A) (Diptera: Drosophilidae: Drosophila simulans; NC_012416) , wAlb (B) (Diptera: Culicidae: Aedes albopictus; CAGB01000001-CAGB01000165) , wVit (B) (Hymenoptera: Pteromalidae: Nasonia vitripennis; AERW00000000) , wHa (A) (Diptera: Drosophilidae: Drosophila simulans; CP003884) , and wNo (B) (Diptera: Drosophilidae: Drosophila simulans; CP003883) . At the time of this study, there were nine available Lepidoptera genomes that were used to search for possible LGT events: Bombyx mori , Danaus plexippus , Heliconius melpomene , Manduca sexta (http://agripestbase.org/manduca), Melitaea cinxia , Papilio glaucus , P. polytes, P. xuthus  and Plutella xylostella .
MLST strain diversity in Lepidoptera
All Wolbachia strains with known associated lepidopteran hosts were grouped in either Supergroup A or B (Additional file 2: Table S2). The majority of lepidopteran strains (76 total representing 32 unique MLST strains) belong to Supergroup B; the remaining (14 total strains representing 6 unique MLST) strains belonging to Supergroup A.
Phylogenetic analysis of MLST strains
In total, 345 Wolbachia strains were analyzed from insect hosts (Coleoptera, Diptera, Hemiptera, Hymenoptera, Isoptera, Orthoptera, Lepidoptera) and distantly related invertebrates (Arachnida, Crustacea, and Nematoda). The ML and ClonalFrame phylogenetic trees were divided into six major clades (Supergroups A-D, F, and H). The ClonalFrame tree also contained an additional clade with strains in Supergroups A, B, C and F; this likely represents sequences that underwent the most recombination. Supergroup A is closely related to Supergroup B (Fig. 1a, b). The strain wExe3, which has a lepidopteran host, was originally classified as A. However, it is basal to clade B with 98 % boostrap support in the ML tree, and it is denoted on Fig. 1 as “A*”. In addition, in the ML tree, strain wHyl, which has an arachnid host, was highly supported (bootstrap = 99) as being basal to the strain wExe3 (labeled “A**”, Fig. 1b). Supergroups A and B, along with A* and A**, were sister to a clade of strains previously classifed as Supergroup H, which further connects to Supergroup D and to Supergroup F. Supergroup C has high support (bootstrap = 85) as being a basal group near the outgroup (Fig. 1b). Most lepidopteran strains were classified in Supergroup B in both the ML and ClonalFrame trees (Fig. 1a). However, in the ClonalFrame tree, A* and A** were grouped in Supergroup A. In the ClonalFrame tree, Supergroup D has high support (bootstrap = 90) and is placed close to outgroups (Fig. 1a).
Gene network analyses of unique Wolbachia strains in Lepidoptera
Strains in Supergroup A were grouped into four networks. Two networks contained only one strain: Network 3 had the lepidopteran strain ST92 (from Ephestia kuehniella [Pyralidae]) and Network 4 had the lepidopteran strain ST223 (from Spodoptera exempta [Erebidae]). Network 2 contained two strains, both with lycanenid butterfly hosts, that were separated by two mutations: ST38 (Jamides alecto) and ST110 (Iraota rochana; Fig. 2b). Network 1 contained nine strains. ST19 was found in eight strains from eight host species: Ephestia kuehniella (Lepidoptera: Pyralidae), Ornipholidotos peucetia (Lepidoptera: Lycaneidae), Ceutorhynchus neglectus (Coleoptera: Curculionidae), and five ant species (Leptogenys sp., Leptomyrmex sp., Pheidole plagiara, P. planifrons, Technomyrmex albipes). The ninth strain, ST91 occurred on the nymphalid butterfly Hypolimnas bolina, and was separated by just one mutation from strain ST19.
Comparison of Wolbachia and Lepidoptera phylogenies
There was no strong congruence between the Wolbachia and lepidopteran phylogenies during mantel test. Analysis of the ML topologies for Wolbachia using JANE and the ML tree from Regier et al.’s  lepidopteran phylogeny at a p-value of 0.05 showed the reconstructions (cost = 92) with only 9 cospeciation events, 22 duplication, 19 duplication and host switching, 22 losses and 10 failure to diverge (Additional file 4: Figure S1).
The Mantel test analysis indicated that there were no significant correlations between the genetic distances of Wolbachia and host Lepidoptera (r = −0.072, P = 0.081 [indigenous]; r = 0.107, P = < 0.0001 [comparing the Wolbachia ClonalFrame tree with the ML tree of Regier et al. ]; r = 0.069, P = 0.019 [comparing the Wolbachia ML tree with the ML tree of Regier et al. ]).
Divergence time estimation
Geography of shared strains
Summary of previous transinfection studies in Lepidoptera
Results of published transinfection experiments of Wolbachia strains performed on lepidopteran hosts
Phenotype in Natural Host
Phenotype in Transinfected Host
1. Strain dependent Phenotype
E. kuehniella (−w)
E. kuehniella (Yokohama)
E. kuehniella (Tsuhiura) (−w)
2. Host dependent Phenotype
(a) Transferable multi potent strain
E. kuehniella (−w)
(b) Non-transferable strain
Bombyx mori (−w)
no stable infection
(c) Population dependent phenotypea
Hypolimnas bolina (Polynesia)
Hypolimnas bolina (South Asian) (+w)
3. Strain/host dependent Phenotype
E. kuehniella (−w)
O. scapulalis (−w)
Lateral gene transfer (LGT)
Comparisons of genomes of Wolbachia and Lepidoptera to test for traces of LGT
Total traces screened (mbp)
Wolbachia traces hitsa
Hits length (bp)
NCBI Blast matches
Previously, vector-mediated interspecific transmission was observed in Wolbachia through shared food sources [2, 102–105], ectoparasitic mites [106, 107], and parasitoids . Our study revealed that inter-specific, inter-familial, and inter-ordinal horizontal transmission is also common in Lepidoptera. Using phylogenetic, co-phylogenetic and network analyses, we found at least seven probable cases of horizontal transmission among 31 host species, both within Lepidoptera and between Lepidoptera and other arthropods. Three strains (ST3, ST40, ST146) were shared among three butterfly families (Lycaenidae, Nymphalidae, Pieridae). One strain (ST125) was shared between two butterfly families (Lycaenidae, Nymphalidae), and the distantly related moth family Noctuidae. Since the majority of lepidopteran larvae feed on plant tissue, and adults obtain nectar from flowers or tree sap, the close association of Lepidoptera with plants might lead to increased infection through host plant mediation . Strain 41 is the most widespread Wolbachia strain in butterflies; it was shared among eleven butterfly species in three families (Lycaenidae, Nymphalidae, Pieridae) and interestingly, it was also shared with one unidentified species of calyptrate fly. There are a number of known hymenopteran parasitoids that are found on both lepidopteran and dipteran hosts, and thus parasitoids may have mediated horizontal transfer .
Another strain, ST37, was found to be shared between the egg parasitoid Tetrastichus coeruleus, the social wasp Polistes dominula, and the lycaenid butterfly Athene emolus. Tetrastichus coeruleus is not known to parasitize lepidopterans. However, it parasitizes eggs of the common asparagus beetle, Crioceris asparagi , which shares a host plant with other Lepidoptera, such as the pest species Spodoptera exigua . Perhaps Wolbachia was transferred into a lepidopteran host through this shared host plant. Larvae of Polistes dominula are parasitoids of Chalcoela iphitalis (Lepidoptera: Pyralidae) , could serve as a possible route of Wolbachia transfer to a lepidopteran host. The Malaysian lycaenid butterfly, Athene emolus, is symbiotic with the ant species Oecophylla smaragdina. These ants guard A. emolus larvae and protect them from predators and parasites . We postulate that any one of these lepidopteran-hymenopteran interactions could potentially enable inter-ordinal transfer of ST37.
Strain ST19 also exhibits inter-ordinal transfer. It is shared among three different insect orders: Lepidoptera (the lycaenid butterfly Ornipholidotos peucetia, and the pyralid moth, Ephestia kuehniella), Hymenoptera (the ant species Leptogenys sp., Leptomyrmex sp., Pheidole planifrons, P. plagiara, and Technomyrmex albipes), and Coleoptera (the weevil Ceutorhynchus neglectus). Horizontal transmission of Wolbachia is also possible when an uninfected insect eats an infected one . Ceutorhynchus neglectus is parasitized by multiple wasps ; weevils also feed on flower pollen and nectar . It is thus possible that ST19 jumped across three insect orders either through shared host plants or via shared parasitoids.
The Mantel test revealed a weak correlation between genetic make-up of lepidopteran host and its endosymbiotic bacteria, Wolbachia, which further support horizontal transmission of Wolbachia within Lepidoptera. Co-phylogenetic analysis revealed common losses, duplication and host switches of Wolbachia strains within Lepidoptera.
We performed divergence time analyses on all available Wolbachia strains from Lepidoptera using two independent calibrations [65, 66]. Results from both calibrations cross-validate our divergence time estimates and suggest the conclusions are robust. Our analysis suggests that Wolbachia was recently introduced in Lepidoptera at a maximum age of ~23 mya. The Wolbachia divergence times, compared to the divergence times estimated by Wahlberg et al. , suggest lepidopteran families that are currently known to carry Wolbachia had already diversified before they became Wolbachia hosts. A recent study on insect evolution suggests the divergence between butterfly and moths and between Lepidoptera and other insect orders (Diptera, Coleoptera and Hymenoptera) took place between ~344-58 mya and the identical strains between them were acquired recently at a maximum of ~23 mya . Our divergence time analysis, in light of the most comprehensive Lepidoptera calibrated phylogeny, suggests that Wolbachia strains ST3, ST19, ST40, ST41, ST125 and ST146, are likely inter-familial horizontal transmissions, and ST125 and ST19 are inter-superfamilial horizontal transmissions [74, 75]. We also found that ST19, ST37, ST41 are clear cases of inter-ordinal horizontal transmission. The cospeciation events predicted in the co-phylogenetic analysis seems to be invalidated, given the lepidopteran estimated divergence times of Wahlberg et al. .
Facultative endosymbionts have already been shown to change host fitness or biology; pea aphids (Acyrthosiphon pisum) have facultative symbionts that protect their hosts against entomopathogenic fungi and parasitoid wasps, ameliorate the detrimental effects of heat, and influence host plant suitability [2, 116–118]. One main consequence of horizontal transmission is induction of unknown phenotypes of Wolbachia into the novel host . A recently discovered Wolbachia strain confers fitness benefits by increasing the resistance against natural pathogens in fruit flies . All previously published transinfection experiments in lepidopteran hosts arrived at similar conclusions that the phenotype induction after transinfection is determined by two factors strains and the host types [94–101]. It is necessary to investigate each strain’s genotype and phenotype in its natural host, as well as other possible hosts in which it may have been transferred through shared resources. In some cases, suppressors against phenotype can lead toward loss of phenotype . Therefore, some species that do not currently induce a phenotype may have done so in the past, implying that more species have had their biology affected by Wolbachia than previously estimated . In other cases, novel hosts can suppress the Wolbachia-mediated phenotype and enable the appearance of hidden phenotypes [100, 101]. Together, these studies suggest that Wolbachia strains possess the genetic makeup to induce multiple phenotypes .
The spread of endosymbionts in field populations by horizontal transmission have received little attention. The mechanisms driving horizontal transmission have mostly remained unclear; even the effects induced by common cases of horizontal transmission are currently unknown [2, 3]. Since there is no way to control horizontal transmission in the field, routes of transmission must be thoroughly studied in order to investigate the genotypes and phenotypes of strains in both natural and novel hosts.
Recently, a complete copy of the Wolbachia genome was found within the genome of Drosophila ananassae and large segments were found in seven other Drosophila species . During the original whole genome sequencing of the nematode, Brugia malayi, extensive levels of lateral gene transfer (LGT) were identified from its Wolbachia endosymbiont . LGT from the Wolbachia genome to the nuclear genome of its eukaryotic hosts is widespread [38, 39]. In a search of sequence data archives, about 70 % of arthropods and nematodes have evidence for LGT from Wolbachia [36, 38, 39]. We found one instance of possible Wolbachia LGT between strain wHa and the butterfly Melitaea cinxia. This result must be confirmed with PCR to rule out the possibility of a genome-sequencing error or contamination. We did not find any evidence of LGT from the Wolbachia genome to the other eight available genomes of Lepidoptera. Even Plutella xylostella, the only species known to have Wolbachia infection, did not yield any evidence of LGT in our analysis of its genome. For M. cinxia, the evidence we found of LGT transmission suggests it is or has been infected with Wolbachia. The method we used to search for possible LGT has previously been used effectively to trace LGT from Wolbachia  and from other bacterial species . The lack of evidence of LGT also supports our inference of a recent introduction of Wolbachia in Lepidoptera. Though these results are sound based on current available data, they are not conclusive; future studies should examine additional genomes and methods to trace LGT in Lepidoptera. The genome assemblies of eukaryotes often filter out bacterial sequences as contaminants and there might be possibility that Wolbachia genes may be present in the original sequencing reads, but not in the finished genome assemblies . We suggest future studies to examine the raw data read instead of assembled genomes to detect those genes, which might have filtered from the original sequencing reads.
Ahmed et al.  found geographic patterns in the infection status of Wolbachia, however, this survey did not find any such patterns in strain distribution. The study frequently found strains distributed across the continents, such as strains ST19, ST37, and ST41, which have been found in multiple hosts across Asia, Africa, Australia and North America. There is no generally accepted theory for how these strains were transferred between various hosts across continents, partially due to the difficulty in tracing the strains’ natural hosts. The comparison of phylogenies of Wolbachia and host Lepidoptera indicates that closely related strains have phylogenetically diverse hosts and vice versa. These examples of shared strains across distantly related families demonstrate that horizontal jumps might be result of recent acquisition of Wolbachia.
Currently, only eight families of Lepidoptera have published Wolbachia strain data. These include three moth families (Crambidae, Noctuidae, Pyralidae) five butterfly families (Hesperiidae, Lycaenidae, Nymphalidae, Papilionidae, Pieridae)that represent three Lepidoptera superfamilies (Noctuoidea, Pyraloidea, Papilionoidea), which contain about 50 % of described lepidopteran species . It would be interesting to explore the Wolbachia strains from other butterfly and moth families, in order to get a comprehensive estimate of the full extent of Wolbachia diversity and mode of transmission within this order.
We found evidence for several new instances of Wolbachia horizontal transmission in Lepidoptera. Our findings suggest that specific shared food sources and shared natural enemies are possible routes of horizontal transmission, but further studies are needed to conclusively determine these routes. We uncover evidence of Wolbachia inducing new phenotypes in novel hosts after horizontal transmission from natural hosts. However, Wolbachia-induced phenotypes have not been well studied for most natural hosts and potential novel hosts. Therefore, it is crucial to study additional Wolbachia-infected organisms in order to determine which species are natural hosts for each strain. It is also important to perform additional transinfection experiments to determine which species can sustain a stable infection. Data from these experiments will yield information about the phenotypes in both natural and novel hosts, revealing new insights into the mechanisms of Wolbachia-induced phenotypic change. Finally, further research into host genotypes should be conducted by analyzing additional genomes of potential hosts to search for the presence of inserted Wolbachia loci, in order to elucidate the function of these laterally transferred genes.
Consent to publish
Availability of data and materials
We provided the data at LabArchive (DOI: 10.6070/H48913W9).
David Plotkin provided useful suggestions; Andrei Sourakov helped with moth and butterfly figures.
This study was supported by the US National Science Foundation grants DEB-1354585, DEB-1541500, and the University of Florida Research Opportunity Seed Fund (ROSF) grant to AYK.
Open AccessThis article is distributed under the terms of the Creative Commons Attribution 4.0 International License (http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The Creative Commons Public Domain Dedication waiver (http://creativecommons.org/publicdomain/zero/1.0/) applies to the data made available in this article, unless otherwise stated.
- Cosmides LM, Tooby J. Cytoplasmic inheritance and intragenomic conflict. J Theor Biol. 1981;89(1):83–129.PubMedView ArticleGoogle Scholar
- Oliver KM, Degnan PH, Burke GR, Moran NA. Facultative symbionts in aphids and the horizontal transfer of ecologically important traits. Annu Rev Entomol. 2010;55:247–66.PubMedView ArticleGoogle Scholar
- Himler AG, Adachi-Hagimori T, Bergen JE, Kozuch A, Kelly SE, Tabashnik BE, Chiel E, Duckworth VE, Dennehy TJ, Zchori-Fein E, et al. Rapid spread of a bacterial symbiont in an invasive whitefly is driven by fitness benefits and female bias. Science. 2011;332(6026):254–6.Google Scholar
- Ahmed MZ, Li S-J, Xue X, Yin X-J, Ren S-X, Jiggins FM, Greeff JM, Qiu B-L. The intracellular bacterium Wolbachia uses parasitoid wasps as phoretic vectors for efficient horizontal transmission. PLoS Pathog. 2015;10(2):e1004672–2.Google Scholar
- Buchner P. Endosymbiosis of animals with plant microorganisms. New York: John Wiley and Sons Interscience; 1965. 332–338. Google Scholar
- Weinert LA, Araujo-Jnr EV, Ahmed MZ, Welch JJ. The incidence of bacterial endosymbionts in terrestrial arthropods. Proc R Soc B Biol Sci. 2015;282(1807):20150249.View ArticleGoogle Scholar
- Hilgenboecker K, Hammerstein P, Schlattmann P, Telschow A, Werren JH. How many species are infected with Wolbachia? - a statistical analysis of current data. FEMS Microbiol Lett. 2008;281(2):215–20.PubMedPubMed CentralView ArticleGoogle Scholar
- Russell JA, Latorre A, Sabater-Munoz B, Moya A, Moran NA. Side-stepping secondary symbionts: widespread horizontal transfer across and beyond the Aphidoidea. Mol Ecol. 2003;12(4):1061–75.PubMedView ArticleGoogle Scholar
- Vavre F, Fleury F, Lepetit D, Fouillet P, Bouletreau M. Phylogenetic evidence for horizontal transmission of Wolbachia in host-parasitoid associations. Mol Biol Evol. 1999;16(12):1711–23.PubMedView ArticleGoogle Scholar
- Noda H, Miyoshi T, Zhang Q, Watanabe K, Deng K, Hoshizaki S. Wolbachia infection shared among planthoppers (Homoptera : Delphacidae) and their endoparasite (Strepsiptera : Elenchidae): a probable case of interspecies transmission. Mol Ecol. 2001;10(8):2101–6.PubMedView ArticleGoogle Scholar
- Shoemaker DD, Machado CA, Molbo D, Werren JH, Windsor DM, Herre EA. The distribution of Wolbachia in fig wasps: correlations with host phylogeny, ecology and population structure. Proc R Soc B Biol Sci. 2002;269(1506):2257–67.View ArticleGoogle Scholar
- Baldo L, Ayoub NA, Hayashi CY, Russell JA, Stahlhut JK, Werren JH. Insight into the routes of Wolbachia invasion: high levels of horizontal transfer in the spider genus Agelenopsis revealed by Wolbachia strain and mitochondrial DNA diversity. Mol Ecol. 2008;17(2):557–69.PubMedView ArticleGoogle Scholar
- Ahmed MZ, De Barro PJ, Ren S-X, Greeff JM, Qiu B-L. Evidence for Horizontal Transmission of Secondary Endosymbionts in the Bemisia tabaci Cryptic Species Complex. PloS One. 2013;8(1):e53084.PubMedPubMed CentralView ArticleGoogle Scholar
- Boyle L, Oneill SL, Robertson HM, Karr TL. Interspecific and intraspecific horizontal transfer of Wolbachia in Drosophila. Science. 1993;260(5115):1796–9.PubMedView ArticleGoogle Scholar
- Heath BD, Butcher RDJ, Whitfield WGF, Hubbard SF. Horizontal transfer of Wolbachia between phylogenetically distant insect species by a naturally occurring mechanism. Curr Biol. 1999;9(6):313–6.PubMedView ArticleGoogle Scholar
- Kang L, Ma X, Cai L, Liao S, Sun L, Zhu H, Chen X, Shen D, Zhao S, Li C. Superinfection of Laodelphax striatellus with Wolbachia from Drosophila simulans. Heredity. 2003;90(1):71–6.Google Scholar
- Zabalou S, Riegler M, Theodorakopoulou M, Stauffer C, Savakis C, Bourtzis K. Wolbachia-induced cytoplasmic incompatibility as a means for insect pest population control. Proc Natl Acad Sci U S A. 2004;101(42):15042–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Riegler M, Charlat S, Stauffer C, Mercot H. Wolbachia transfer from Rhagoletis cerasi to Drosophila simulans: Investigating the outcomes of host-symbiont coevolution. Appl Environ Microbiol. 2004;70(1):273–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Van Nieukerken EJ, Kaila L, Kitching IJ, Kristensen NP, Lees DC, Minet J, Mitter C, Mutanen M, Regier JC, Simonsen TJ et al. Order Lepidoptera Linnaeus, 1758. Zootaxa. 2011;3148:212–21.Google Scholar
- Banziger H. Remarkable new cases of moths drinking human tears in Thailand (Lepidoptera: Thyatritridae, Sphingidae, Notodontidae). Nat Hist Bull Siam Soc. 1992;40(1):91–102.Google Scholar
- Pierce NE. Predatory and parasitic Lepidoptera: Carnivores living on plants. J Lepid Soc. 1995;49(4):412–53.Google Scholar
- Zaspel JM, Weller SJ, Branham MA. A comparative survey of proboscis morphology and associated structures in fruit-piercing, tear-feeding, and blood-feeding moths in Calpinae (Lepidoptera: Erebidae). Zoomorphology. 2011;130(3):203–25.View ArticleGoogle Scholar
- Plotkin D, Goddard J. Blood, sweat, and tears: a review of the hematophagous, sudophagous, and lachryphagous Lepidoptera. J Vector Ecol. 2013;38(2):289–94.PubMedView ArticleGoogle Scholar
- Roe AD, Weller SJ, Baixeras J, Brown J, Cummings MP, Davis DR, Kawahara AY, Parr CS, Regier JC, Rubinoff D, et al. Evolutionary framework for Lepidoptera model systems. 2010.Google Scholar
- Feener Jr DH, Brown BV. Diptera as parasitoids. Annu Rev Entomol. 1997;42:73–97.PubMedView ArticleGoogle Scholar
- Whitfield JB. Phylogeny and evolution of host-parasitoid interactions in hymenoptera. Annu Rev Entomol. 1998;43:129–51.PubMedView ArticleGoogle Scholar
- Pennacchio F, Strand MR. Evolution of developmental strategies in parasitic hymenoptera. Annu Rev Entomol. 2006;51:233–58.PubMedView ArticleGoogle Scholar
- Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008;6(10):741–51.PubMedView ArticleGoogle Scholar
- Ahmed MZ, Araujo-Jnr EV, Welch JJ, Kawahara AY. Wolbachia in butterflies and moths: geographic structure in infection frequency. Front Zool. 2015;12:16.PubMedPubMed CentralView ArticleGoogle Scholar
- Ahmed MZ, Greyvenstein OFC, Erasmus C, Welch JJ, Greeff JM. Consistently high incidence of Wolbachia in global fig wasp communities. Ecol Entomol. 2013;38(2):147–54.View ArticleGoogle Scholar
- Sasaki T, Ishikawa H. Wolbachia infections and cytoplasmic incompatibility in the almond moth and the mediterranean flour moth. Zool Sci. 1999;16(5):739–44.View ArticleGoogle Scholar
- Kageyama D, Nishimura G, Hoshizaki S, Ishikawa Y. Feminizing Wolbachia in an insect, Ostrinia furnacalis (Lepidoptera : Crambidae). Heredity. 2002;88:444–9.PubMedView ArticleGoogle Scholar
- Dyson EA, Hurst GDD. Persistence of an extreme sex-ratio bias in a natural population. Proc Natl Acad Sci U S A. 2004;101(17):6520–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Graham RI, Grzywacz D, Mushobozi WL, Wilson K. Wolbachia in a major African crop pest increases susceptibility to viral disease rather than protects. Ecol Lett. 2012;15(9):993–1000.PubMedView ArticleGoogle Scholar
- Hughes GL, Rasgon JL. Transinfection: a method to investigate Wolbachia-host interactions and control arthropod-borne disease. Insect Mol Biol. 2014;23(2):141–51.PubMedPubMed CentralView ArticleGoogle Scholar
- Hotopp JCD, Clark ME, Oliveira DCSG, Foster JM, Fischer P, Munoz Torres MC, Giebel JD, Kumar N, Ishmael N, Wang S et al. Widespread lateral gene transfer from intracellular bacteria to multicellular eukaryotes. Science. 2007;317(5845):1753–6.View ArticleGoogle Scholar
- Gyles C, Boerlin P. Horizontally transferred genetic elements and their role in pathogenesis of bacterial disease. Vet Pathol. 2014;51(2):328–40.PubMedView ArticleGoogle Scholar
- Hotopp JCD. Horizontal gene transfer between bacteria and animals. Trends Genet. 2011;27(4):157–63.View ArticleGoogle Scholar
- Robinson KM, Sieber KB, Hotopp JCD. A review of bacteria-animal lateral gene transfer may inform our understanding of diseases like cancer. PloS Genet. 2013;9(10):e1003877.PubMedPubMed CentralView ArticleGoogle Scholar
- Wheeler D, Redding AJ, Werren JH. Characterization of an ancient lepidopteran lateral gene transfer. PloS One. 2013;8(3):e59262.PubMedPubMed CentralView ArticleGoogle Scholar
- Baldo L, Hotopp JCD, Jolley KA, Bordenstein SR, Biber SA, Choudhury RR, Hayashi C, Maiden MCJ, Tettelin H, Werren JH. Multilocus sequence typing system for the endosymbiont Wolbachia pipientis. Appl Environ Microbiol. 2006;72(11):7098–110.PubMedPubMed CentralView ArticleGoogle Scholar
- Jolley KA, Maiden MCJ. BIGSdb: Scalable analysis of bacterial genome variation at the population level. BMC Bioinf. 2010;11:595.View ArticleGoogle Scholar
- Katoh K, Kuma K, Toh H, Miyata T. MAFFT version 5: improvement in accuracy of multiple sequence alignment. Nucleic Acids Res. 2005;33(2):511–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Kearse M, Moir R, Wilson A, Stones-Havas S, Cheung M, Sturrock S, Buxton S, Cooper A, Markowitz S, Duran C et al. Geneious Basic: An integrated and extendable desktop software platform for the organization and analysis of sequence data. Bioinformatics. 2012;28(12):1647–9.PubMedPubMed CentralView ArticleGoogle Scholar
- Lanfear R, Calcott B, Ho SYW, Guindon S. PartitionFinder: Combined selection of partitioning schemes and substitution models for phylogenetic analyses. Mol Biol Evol. 2012;29(6):1695–701.PubMedView ArticleGoogle Scholar
- Stamatakis A. RAxML version 8: a tool for phylogenetic analysis and post-analysis of large phylogenies. Bioinformatics. 2014;30(9):1312–3.PubMedPubMed CentralView ArticleGoogle Scholar
- Kawahara AY, Breinholt JW, Ponce FV, Haxaire J, Xiao L, Lamarre GPA, Rubinoff D, Kitching IJ. Evolution of Manduca sexta hornworms and relatives: Biogeographical analysis reveals an ancestral diversification in Central America. Mol Phylogenet Evol. 2013;68(3):381–6.PubMedView ArticleGoogle Scholar
- Guindon S, Dufayard J-F, Lefort V, Anisimova M, Hordijk W, Gascuel O. New Algorithms and Methods to Estimate Maximum-Likelihood Phylogenies: Assessing the Performance of PhyML 3.0. Syst Biol. 2010;59(3):307–21.PubMedView ArticleGoogle Scholar
- Didelot X, Falush D. Inference of bacterial microevolution using multilocus sequence data. Genetics. 2007;175(3):1251–66.PubMedPubMed CentralView ArticleGoogle Scholar
- Gelman A, DB R. Inference from iterative simulation using multiple sequences. Stat Sci. 1992;7:457–72.View ArticleGoogle Scholar
- Hart MW, Sunday J. Things fall apart: biological species form unconnected parsimony networks. Biol Lett. 2007;3(5):509–12.PubMedPubMed CentralView ArticleGoogle Scholar
- Chen H, Strand M, Norenburg JL, Sun S, Kajihara H, Chernyshev AV, Maslakova SA, Sundberg P. Statistical parsimony networks and species assemblages in Cephalotrichid Nemerteans (Nemertea). PloS One. 2010;5(9):e12885.PubMedPubMed CentralView ArticleGoogle Scholar
- De Barro P, Ahmed MZ. Genetic networking of the Bemisia tabaci cryptic species complex reveals pattern of biological invasions. PloS One. 2011;6(10):e25579.PubMedPubMed CentralView ArticleGoogle Scholar
- Pfarr K, Foster J, Slatko B, Hoerauf A, Eisen JA. On the taxonomic status of the intracellular bacterium Wolbachia pipientis: should this species name include the intracellular bacteria of filarial nematodes? Int J Syst Evol Microbiol. 2007;57:1677–8.PubMedPubMed CentralView ArticleGoogle Scholar
- Lo N, Paraskevopoulos C, Bourtzis K, O’Neill SL, Werren JH, Bordenstein SR, Bandi C. Taxonomic status of the intracellular bacterium Wolbachia pipientis. Int J Syst Evol Microbiol. 2007;57:654–7.PubMedView ArticleGoogle Scholar
- Posada D, Crandall KA. Intraspecific gene genealogies: trees grafting into networks. Trends Ecol Evol. 2001;16(1):37–45.PubMedView ArticleGoogle Scholar
- Clement M, Posada D, Crandall KA. TCS: a computer program to estimate gene genealogies. Mol Ecol. 2000;9(10):1657–9.PubMedView ArticleGoogle Scholar
- Templeton AR, Crandall KA, Sing CF. A Cladistic Analysis of phenotypic associations with haplotypes inferred from restriction endonuclease mapping and DNA sequence data. III. cladogram estimation. Genetics. 1992;132(2):619–33.PubMedPubMed CentralGoogle Scholar
- Regier JC, Mitter C, Zwick A, Bazinet AL, Cummings MP, Kawahara AY, Sohn J-C, Zwickl DJ, Cho S, Davis DR et al. A large-scale, higher-level, molecular phylogenetic study of the insect order lepidoptera (Moths and Butterflies). PLoS ONE. 2013;8(3):1–23.View ArticleGoogle Scholar
- Smouse PE, Long JC, Sokal RR. Multiple regression and correlation extensions of the mantel test of matrix correspondence. Syst Zool. 1986;35(4):627–32.View ArticleGoogle Scholar
- Maddison WP MD. Mesquite: A modular system for evolutionary analysis. Version 3.04. Available from http://mesquiteproject.org. 2015.
- Conow C, Fielder D, Ovadia Y, Libeskind-Hadas R. Jane: a new tool for the cophylogeny reconstruction problem. Algorithms Mol Biol. 2010;5:16.PubMedPubMed CentralView ArticleGoogle Scholar
- Bouckaert R, Heled J, Kuehnert D, Vaughan T, Wu C-H, Xie D, Suchard MA, Rambaut A, Drummond AJ. BEAST 2: A Software Platform for Bayesian Evolutionary Analysis. PloS Comput Biol. 2014;10(4):e1003537.PubMedPubMed CentralView ArticleGoogle Scholar
- Drummond AJ, Suchard MA, Xie D, Rambaut A. Bayesian Phylogenetics with BEAUti and the BEAST 1.7. Mol Biol Evol. 2012;29(8):1969–73.PubMedPubMed CentralView ArticleGoogle Scholar
- Richardson MF, Weinert LA, Welch JJ, Linheiro RS, Magwire MM, Jiggins FM, Bergman CM. Population genomics of the Wolbachia endosymbiont in Drosophila melanogaster. PLoS Genet. 2012;8(12):e1003129.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerth M, Roethe J, Bleidorn C. Tracing horizontal Wolbachia movements among bees (Anthophila): a combined approach using multilocus sequence typing data and host phylogeny. Mol Ecol. 2013;22(24):6149–62.PubMedView ArticleGoogle Scholar
- Felsenstein J. Evolutionary trees from DNA sequences: a maximum likelihood approach. J Mol Evol. 1981;17(6):368–76.PubMedView ArticleGoogle Scholar
- Swofford DL, Sullivan J. Phylogeny inference based on parsimony and other methods using PAUP*. 2009.View ArticleGoogle Scholar
- Posada D. Unveiling the molecular clock in the presence of recombination. Mol Biol Evol. 2001;18(10):1976–8.PubMedView ArticleGoogle Scholar
- Pond SLK, Frost SDW, Muse SV. HyPhy: hypothesis testing using phylogenies. Bioinformatics. 2005;21(5):676–9.PubMedView ArticleGoogle Scholar
- Boni MF, Posada D, Feldman MW. An exact nonparametric method for inferring mosaic structure in sequence triplets. Genetics. 2007;176(2):1035–47.PubMedPubMed CentralView ArticleGoogle Scholar
- Gerth M, Bleidorn C. A multilocus sequence typing (MLST) approach to diminish the problems that are associated with DNA barcoding: A reply to Stahlhut etal. 2012. Syst Biodivers. 2013;11(1):15–7.View ArticleGoogle Scholar
- Rambaut A, Suchard M, Xie D, Drummond A. Tracer v1.6, Available from http://beast.bio.ed.ac.uk/Tracer. 2014. Accessed 23 Apr. 2016.
- Wahlberg N, Wheat CW, Pena C. Timing and patterns in the taxonomic diversification of lepidoptera (Butterflies and Moths). PloS One. 2013;8(11):e80875.PubMedPubMed CentralView ArticleGoogle Scholar
- Misof B, Liu S, Meusemann K, Peters RS, Donath A, Mayer C, Frandsen PB, Ware J, Flouri T, Beutel RG et al. Phylogenomics resolves the timing and pattern of insect evolution. Science. 2014;346(6210):763–7.PubMedView ArticleGoogle Scholar
- Delcher AL, Phillippy A, Carlton J, Salzberg SL. Fast algorithms for large-scale genome alignment and comparison. Nucleic Acids Res. 2002;30(11):2478–83.PubMedPubMed CentralView ArticleGoogle Scholar
- Foster J, Ganatra M, Kamal I, Ware J, Makarova K, Ivanova N, Bhattacharyya A, Kapatral V, Kumar S, Posfai J et al. The Wolbachia genome of Brugia malayi: Endosymbiont evolution within a human pathogenic nematode. PloS Biol. 2005;3(4):599–614.View ArticleGoogle Scholar
- Duplouy A, Iturbe-Ormaetxe I, Beatson SA, Szubert JM, Brownlie JC, McMeniman CJ, McGraw EA, Hurst GDD, Charlat S, O'Neill SL et al. Draft genome sequence of the male-killing Wolbachia strain wBol1 reveals recent horizontal gene transfers from diverse sources. BMC Genomics. 2013;14:20.Google Scholar
- Wu M, Sun LV, Vamathevan J, Riegler M, Deboy R, Brownlie JC, McGraw EA, Martin W, Esser C, Ahmadinejad N Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: A streamlined genome overrun by mobile genetic elements. PloS Biol. 2004;2(3):327–41.View ArticleGoogle Scholar
- Klasson L, Walker T, Sebaihia M, Sanders MJ, Quail MA, Lord A, Sanders S, Earl J, O'Neill SL, Thomson N et al. Genome evolution of Wolbachia strain wPip from the Culex pipiens group. Mol Biol Evol. 2008;25(9):1877–87.Google Scholar
- Klasson L, Westberg J, Sapountzis P, Naesiund K, Lutnaes Y, Darby AC, Veneti Z, Chen L, Braig HR, Garrett R et al. The mosaic genome structure of the Wolbachia wRi strain infecting Drosophila simulans. Proc Natl Acad Sci U S A. 2009;106(14):5725–30.PubMedPubMed CentralView ArticleGoogle Scholar
- Mavingui P, Moro CV, Van T-V, Wisniewski-Dye F, Raquin V, Minard G, Tran F-H, Voronin D, Rouy Z, Bustos P et al. Whole-genome sequence of Wolbachia strain wAlbB, an endosymbiont of tiger mosquito vector Aedes albopictus. J Bacteriol. 2012;194(7):1840–0.Google Scholar
- Kent BN, Salichos L, Gibbons JG, Rokas A, Newton ILG, Clark ME, Bordenstein SR. Complete bacteriophage transfer in a bacterial endosymbiont (Wolbachia) determined by targeted genome capture. Genome Biol Evol. 2011;3:209–18.PubMedPubMed CentralView ArticleGoogle Scholar
- Ellegaard KM, Klasson L, Naslund K, Bourtzis K, Andersson SGE. Comparative genomics of Wolbachia and the bacterial species concept. PloS Genet. 2013;9(4):e1003381.PubMedPubMed CentralView ArticleGoogle Scholar
- Mita K, Kasahara M, Sasaki S, Nagayasu Y, Yamada T, Kanamori H, Namiki N, Kitagawa M, Yamashita H, Yasukochi Y et al. The genome sequence of silkworm, Bombyx mori. DNA Res. 2004;11(1):27–35.PubMedView ArticleGoogle Scholar
- Zhan S, Merlin C, Boore JL, Reppert SM. The monarch butterfly genome yields insights into long-distance migration. Cell. 2011;147(5):1171–85.PubMedPubMed CentralView ArticleGoogle Scholar
- Dasmahapatra KK, Walters JR, Briscoe AD, Davey JW, Whibley A, Nadeau NJ, Zimin AV, Hughes DST, Ferguson LC, Martin SH et al. Butterfly genome reveals promiscuous exchange of mimicry adaptations among species. Nature (London). 2012;487(7405):94–8.Google Scholar
- Ahola V, Lehtonen R, Somervuo P, Salmela L, Koskinen P, Rastas P, Välimäki N, Paulin L, Kvist J, Wahlberg N et al.: The Glanville fritillary genome retains an ancient karyotype and reveals selective chromosomal fusions in Lepidoptera. Nature Communications 2014:in press.Google Scholar
- Cong Q, Borek D, Otwinowski Z, Grishin NV. Tiger Swallowtail Genome Reveals Mechanisms for Speciation and Caterpillar Chemical Defense. Cell Rep. 2015;10(6):910–9.View ArticleGoogle Scholar
- Nishikawa H, Iijima T, Kajitani R, Yamaguchi J, Ando T, Suzuki Y, Sugano S, Fujiyama A, Kosugi S, Hirakawa H et al. A genetic mechanism for female-limited Batesian mimicry in Papilio butterfly. Nat Genet. 2015;47(4):405–U169.PubMedView ArticleGoogle Scholar
- Tang W, Yu L, He W, Yang G, Ke F, Baxter SW, You S, Douglas CJ, You M. DBM-DB: the diamondback moth genome database. Database. 2014;2014:bat087.PubMedPubMed CentralView ArticleGoogle Scholar
- Vos M, Didelot X. A comparison of homologous recombination rates in bacteria and archaea. ISME J. 2009;3(2):199–208.PubMedView ArticleGoogle Scholar
- Kawahara AY, Breinholt JW. Phylogenomics provides strong evidence for relationships of butterflies and moths. Proc R Soc B Biol Sci. 2014;281(1788):20140970.View ArticleGoogle Scholar
- Fujii Y, Kageyama D, Hoshizaki S, Ishikawa H, Sasaki T. Transfection of Wolbachia in Lepidoptera: the feminizer of the adzuki bean borer Ostrinia scapulalis causes male killing in the Mediterranean flour moth Ephestia kuehniella. Proc R Soc B Biol Sci. 2001;268(1469):855–9.View ArticleGoogle Scholar
- Kageyama D, Traut W. Opposite sex-specific effects of Wolbachia and interference with the sex determination of its host Ostrinia scapulalis. Proc R Soc B Biol Sci. 2004;271(1536):251–8.View ArticleGoogle Scholar
- Sasaki T, Ishikawa H. Transinfection of Wolbachia in the Mediterranean flour moth, Ephestia kuehniella, by embryonic microinjection. Heredity. 2000;85(2):130–5.PubMedView ArticleGoogle Scholar
- Sakamoto H, Ishikawa Y, Sasaki T, Kikuyama S, Tatsuki S, Hoshizaki S. Transinfection reveals the crucial importance of Wolbachia genotypes in determining the type of reproductive alteration in the host. Genet Res. 2005;85(3):205–10.PubMedView ArticleGoogle Scholar
- Sasaki T, Kubo T, Ishikawa H. Interspecific transfer of Wolbachia between two lepidopteran insects expressing cytoplasmic incompatibility: A Wolbachia variant naturally infecting Cadra cautella causes male killing in Ephestia kuehniella. Genetics. 2002;162(3):1313–9.PubMedPubMed CentralGoogle Scholar
- Kageyama D, Narita S, Noda H. Transfection of Feminizing Wolbachia Endosymbionts of the Butterfly, Eurema hecabe, into the Cell Culture and Various Immature Stages of the Silkmoth, Bombyx mori. Microb Ecol. 2008;56(4):733–41.PubMedView ArticleGoogle Scholar
- Hornett EA, Charlat S, Duplouy AMR, Davies N, Roderick GK, Wedell N, Hurst GDD. Evolution of male-killer suppression in a natural population. PLoS Biol. 2006;4(9):e283.PubMedPubMed CentralView ArticleGoogle Scholar
- Hornett EA, Duplouy AMR, Davies N, Roderick GK, Wedell N, Hurst GDD, Charlat S. You can’t keep a good parasite down: Evolution of a male-killer suppressor uncovers cytoplasmic incompatibility. Evolution. 2008;62(5):1258–63.PubMedView ArticleGoogle Scholar
- Huigens ME, Luck RF, Klaassen RHG, Maas M, Timmermans M, Stouthamer R. Infectious parthenogenesis. Nature. 2000;405(6783):178–9.PubMedView ArticleGoogle Scholar
- Huigens ME, de Almeida RP, Boons PAH, Luck RF, Stouthamer R. Natural interspecific and intraspecific horizontal transfer of parthenogenesis-inducing Wolbachia in Trichogramma wasps. Proc R Soc B Biol Sci. 2004;271(1538):509–15.View ArticleGoogle Scholar
- Rigaud T, Juchault P. Success and failure of horizontal transfers of feminizing Wolbachia endosymbionts in woodlice. J Evol Biol. 1995;8(2):249–55.View ArticleGoogle Scholar
- Sintupachee S, Milne JR, Poonchaisri S, Baimai V, Kittayapong P. Closely related Wolbachia strains within the pumpkin arthropod community and the potential for horizontal transmission via the plant. Microb Ecol. 2006;51(3):294–301.PubMedView ArticleGoogle Scholar
- Gehrer L, Vorburger C. Parasitoids as vectors of facultative bacterial endosymbionts in aphids. Biol Lett. 2012;8(4):613–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Jaenike J, Polak M, Fiskin A, Helou M, Minhas M. Interspecific transmission of endosymbiotic Spiroplasma by mites. Biol Lett. 2007;3(1):23–5.PubMedPubMed CentralView ArticleGoogle Scholar
- Apiwathnasorn C. Literature review of parasitoids of filth flies in Thailand: a list of species with brief notes on bionomics of common species. Southeast Asian J Trop Med Public Health. 2012;43(1):48–54.PubMedGoogle Scholar
- Poll JTK, Alphen JJM, Driessen GJJ. Biological control of the asparagus beetle (Crioceris asparagi) using Tetrastichus asparagi. Proc Sect Exp Appl Entomol Neth Entomol Soc. 1998;9:129–30.Google Scholar
- Jha RK, Tuan S, Chi H, Tang L. Life table and consumption capacity of corn earworm, Helicoverpa armigera, fed asparagus, Asparagus officinalis. J Insect Sci. 2014;14:34 (1 March 2014)-(2011 March 2014).PubMedPubMed CentralView ArticleGoogle Scholar
- Madden AA, Davis MM, Starks PT. First detailed report of brood parasitoidism in the invasive population of the paper wasp Polistes dominulus (Hymenoptera, Vespidae) in North America. Insect Soc. 2010;57(3):257–60.View ArticleGoogle Scholar
- Waldbauer G. What good are bugs? Insects in the web of life. 2003.Google Scholar
- Le Clec’h W, Chevalier FD, Genty L, Bertaux J, Bouchon D, Sicard M. Cannibalism and predation as paths for horizontal passage of Wolbachia between terrestrial isopods. PloS One. 2013;8(4):e60232.PubMedPubMed CentralView ArticleGoogle Scholar
- Mason PG, Miall JH, Bouchard P, Brauner A, Gillespie DR, Gibson GAP. The parasitoid communities associated with Ceutorhynchus species (Coleoptera: Curculionidae) in Ontario and Quebec, Canada. Can Entomol. 2014;146(2):224–35.View ArticleGoogle Scholar
- Juran I, Culjak TG, Grubisic D. Rape stem weevil (Ceutorhynchus napi Gyll. 1837) and cabbage stem weevil (Ceutorhynchus pallidactylus Marsh. 1802) (Coleoptera: Curculionidae) - Important oilseed rape pests. Agric Conspec Sci. 2011;76(2):93–100.Google Scholar
- Oliver KM, Moran NA, Hunter MS. Variation in resistance to parasitism in aphids is due to symbionts not host genotype. Proc Natl Acad Sci U S A. 2005;102(36):12795–800.PubMedPubMed CentralView ArticleGoogle Scholar
- Scarborough CL, Ferrari J, Godfray HCJ. Aphid protected from pathogen by endosymbiont. Science. 2005;310(5755):1781–1.Google Scholar
- Xue X, Li S-J, Ahmed MZ, De Barro PJ, Ren S-X, Qiu B-L. Inactivation of Wolbachia Reveals Its Biological Roles in Whitefly Host. Plos One. 2012;7(10):e48148.PubMedPubMed CentralView ArticleGoogle Scholar
- Teixeira L, Ferreira A, Ashburner M. The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. Plos Biology. 2008;6(12):2753–63.View ArticleGoogle Scholar
- Derks MFL, Smit S, Salis L, Schijlen E, Bossers A, Mateman C, Pijl AS, de Ridder D, Groenen MAM, Visser ME et al. The Genome of Winter Moth (Operophtera brumata) Provides a Genomic Perspective on Sexual Dimorphism and Phenology. Genome Biology and Evolution. 2015;7(8):2321–32.PubMedPubMed CentralView ArticleGoogle Scholar