Wolbachia in butterflies and moths: geographic structure in infection frequency
© Ahmed et al. 2015
Received: 24 April 2015
Accepted: 18 May 2015
Published: 16 July 2015
Butterflies and moths (Lepidoptera) constitute one of the most diverse insect orders, and play an important role in ecosystem function. However, little is known in terms of their bacterial communities. Wolbachia, perhaps the most common and widespread intracellular bacterium on Earth, can manipulate the physiology and reproduction of its hosts, and is transmitted vertically from mother to offspring, or sometimes horizontally between species. While its role in some hosts has been studied extensively, its incidence across Lepidoptera is poorly understood. A recent analysis using a beta-binomial model to infer the between-species distribution of prevalence estimated that approximately 40 % of arthropod species are infected with Wolbachia, but particular taxonomic groups and ecological niches seem to display substantially higher or lower incidences. In this study, we took an initial step and applied a similar, maximum likelihood approach to 300 species of Lepidoptera (7604 individuals from 660 populations) belonging to 17 families and 10 superfamilies, and sampled from 36 countries, representing all continents excluding Antarctica.
Approximately a quarter to a third of individuals appear to be infected with Wolbachia, and around 80 % of Lepidoptera species are infected at a non-negligible frequency. This incidence estimate is very high compared to arthropods in general. Wolbachia infection in Lepidoptera is shown to vary between families, but there is no evidence for closely related groups to show similar infection levels. True butterflies (Papilionoidea) are overrepresented in our data, however, our estimates show this group can be taken as a representative for the other major lepidopteran superfamilies. We also show substantial variation in infection level according to geography – closer locations tend to show similar infection levels. We further show that variation in geography is due to a latitudinal gradient in Wolbachia infection, with lower frequencies towards higher latitudes.
Our comprehensive survey of Wolbachia infection in Lepidoptera suggests that infection incidence is very high, and provides evidence that climate and geography are strong predictors of infection frequency.
KeywordsBacteria Butterfly Latitudinal gradient Moth
Butterflies and moths (Lepidoptera) constitute one of the most diverse insect orders with more than 157,000 described species . Lepidoptera play an important role in ecosystem function primarily as pollinators and herbivores, though some species feed on blood and other animal secretions [2–5]. The order includes many significant agricultural pests, and some species serve as models across biological disciplines . Furthermore, lepidopteran larvae are hosts to another major insect radiation – the parasitic flies and wasps [7–9]. Despite the diversity of Lepidoptera and their many associations with other organisms, little is known about the bacterial community associated with the order.
Wolbachia (Alphaproteobacteria: Rickettsiales: Rickettsiaceae) is a genus of intracellular bacteria that infects many arthropods and nematodes , and is considered the most widespread endosymbiotic bacterium on Earth [10–13]. Although there is some evidence of mutualism, it is well known that Wolbachia can be involved in the parasitic manipulation of its host’s reproductive system . Wolbachia is usually transmitted vertically from mother to offspring; horizontal transmission is also possible [10, 14–16]. Wolbachia infection frequency, both within species (prevalence) and among species (incidence), is a crucial parameter for understanding the dynamic processes behind host-endosymbiont interactions. Infection frequencies vary with transmission efficiencies, and the phenotypic effects of infection . Some Wolbachia phenotypes increase the infection frequency in the host population whereas others decrease it . Our study lays the groundwork for further biological investigations of the effects of Wolbachia on host Lepidoptera.
A pioneering study by Werren et al.  concluded that 17 % of neotropical insect species are infected with Wolbachia. A subsequent, expanded study by Werren and Windsor  concluded that at least 20 % of all insect species are infected. These studies reported sample incidence (i.e., the proportion of samples infected), but many of the samples included a small number of individuals, making low prevalence infections hard to detect. To solve this problem, Hilgenboecker et al.  used beta-binomial modeling to infer the between-species distribution of prevalence from the screen data. Using 20 PCR screens for Wolbachia, their best-fitting distribution of prevalence had peaks at very low and very high prevalence values – suggesting that many low prevalence infections were missed by previous studies. They subsequently calculated that 66 % of arthropod species were infected. Zug and Hammerstein  recently applied the same model to data from Duron et al.  and calculated a 40 % incidence in arthropods.
Wolbachia infection frequencies can have substantial variation both within and among particular taxonomic groups or ecological niches [19, 20, 22]. For example, Ahmed et al.  applied the beta-binomial approach to screening data from 172 species of fig wasps, and estimated a Wolbachia prevalence distribution with peaks at the two extreme values. They further estimated >80 % incidence for fig wasps, considerably higher than all previous estimates for Wolbachia infection frequency across arthropods . Species of Lepidoptera also exhibit extreme variation in Wolbachia prevalence, such as the butterfly Hypolimnas bolina (50–100 %) [23, 24] and the moth Plutella xylostella (0–40 %) . However, despite these studies, there is little understanding of why infection levels vary. Although much work has been conducted on Wolbachia-Lepidoptera interactions [19, 20, 23–29], broad patterns of Wolbachia prevalence and incidence in the order have not been thoroughly characterised. In this study, we synthesize prior work on Wolbachia in Lepidoptera and perform a new analysis of the combined data set using the beta-binomial approach. In addition, we test for variation in infection patterns between taxonomic groups and geographic regions.
Materials and methods
Summary of Wolbachia infection percentages in Lepidoptera
No. of populations
No. of species
No. of individuals
A. Host Taxonomy (Superfamily)
B. Host Geography (Continents)
Estimators of Wolbachia prevalence and incidence
Test for predictors of Wolbachia infection
We performed Mantel tests to assess whether phylogenetic distance among lepidopteran families or geography is correlated with the absolute differences in our moment-based estimates of their mean prevalence (μ) and/or incidence (x c ) of Wolbachia infection. To avoid pseudoreplication, we used randomization to calculate significance. These tests are distribution free, and so have reduced power compared to a full model-based parametric test. However, standard Brownian motion models are not appropriate for incidence and prevalence estimates, and no obvious model-based alternative exists. These tests used the “mantel.rtest” function in the “ade4” package in R version 2.13.0 .
To calculate geographical distance between locations, we used coordinates from the estimated midpoints of the countries where the samples were collected, and these were converted into distances using the Meeus method, which assumes an ellipsoid shape for the planet, in the R “geosphere” package . To calculate evolutionary distances, we used patristic distances as found in the phylogeny of Regier et al. .
Test for sample size bias
A major problem with comparative analyses is the sampling bias that arises from a focus on populations or species that are already known to contain infection . To test for potential bias in our data, we used two approaches. First, we binned screens by sample size (0–10, 10–20, 20–30 etc.) and used Spearman’s rank correlation to test for an association between sample size and the estimated incidence values. Second, we used standard binomial regression to ask whether the probability an individual being infected correlated positively with the population sample size.
We assembled a data set from the literature, containing screens of 7604 lepidopteran individuals from 660 populations covering a broad taxonomic and geographic range (Table 1). Initial tests for sampling bias showed no tendency for well-sampled populations to have higher levels of infection prevalence in these data (Spearman’s rank correlation, p = 0.67; Slope of best-fit binomial regression = −0.00081).
Wolbachia incidence and prevalence in Lepidoptera
Samples from 129 of 300 Lepidoptera species were infected with Wolbachia, resulting in a sample incidence of 43 % (Table 1, Additional file 2: Table S2). However, as with previous studies, these estimates are not very thorough because they ignore low prevalence infections, severely underestimating the true incidence .
Maximum likelihood estimates of levels of Wolbachia infection in Lepidoptera
Mean prevalence, μ
Incidence, x 0.001
Complete database; beta distribution (eq. 3)
0.28 (0.25, 0.31)
0.84 (0.77, 0.90)
Complete database; doubly inflated distribution (eq. 4)
0.27 (0.24, 0.31)
0.77 (0.63, 1.00)
Complete database; beta distribution
0.26 (0.23, 0.29)
0.81 (0.72, 0.88)
Complete database; doubly inflated distribution
0.25 (0.22, 0.29)
0.88 (0.65, 0.95)
Multiple primers; beta distribution
0.32 (0.25, 0.40)
0.86 (0.72, 0.98)
These estimates are open to criticism in several respects. First, the beta distribution might not be flexible enough to fit the true distribution of prevalences in nature. Indeed, the raw data (black bars in Fig. 1), suggest a trimodal distribution, which cannot be accommodated by the beta distribution (eq. (3)). To address this criticism, we followed Weinert et al.  and fit the more flexible doubly-inflated beta distribution (eq. (4)), which can fit trimodal distributions. This more complex distribution was favoured by model selection criterion (Akaike Information Criterion: beta 1439.79, doubly-inflated 1432.14). Nevertheless, the estimate of mean prevalence was little changed, although the estimate of incidence was reduced from 84 % to 77 % (Table 2).
Another serious criticism of our estimates is the highly biased sampling of our database, and the quality of the individual data points. As our database is dominated by true butterflies (Superfamily Papilionoidea; Table 1), we first carried out separate estimates for this Superfamily (Table 2). Again, we estimated a high incidence, which in this case was increased fitting the favoured doubly-inflated distribution (Table 2; Akaike Information Criterion: beta 1031.83, doubly-inflated 1025.271). Furthermore, this estimate was almost unchanged if we restricted our analysis to the subset of the data obtained with multiple primer sets, which are less likely to be affected by false positives due to amplification of contaminants (; Table 2). Therefore, we can be relatively confident of our estimates for the true butterflies, although results reported below suggest that they may also be broadly representative of the other major lepidopteran superfamilies in our database (see below).
Wolbachia infection in Lepidoptera is predictable from host geography
No phylogenetic signal for Wolbachia infection in Lepidoptera
To determine whether these estimates showed phylogenetic signal (i.e., whether there was a tendency for closely-related lepidopteran groups to show similar levels of incidence and prevalence), we used a Mantel test on the family-level data, comparing differences in levels of Wolbachia infection to the phylogenetic distance between the families. While correlations were positive, we found no significant relationship in either prevalence or incidence case (n =10; μ: r = 0.12, P = 0.20; n = 9; x 0.001 : r = 0.054, P = 0.43; x 0.0001 : r = 0.067, P = 0.41) (Fig. 3).
Lepidoptera are among the most diverse orders of insects, with more than 157,000 described species . Lepidoptera are also among the most widely distributed, inhabiting all terrestrial biomes except Antarctica. Their larvae are predominantly herbivores, and their adults play a key role in ecological systems as pollinators . Some species rank among the most destructive agricultural pests, causing significant damage to crops, stored products and natural forests . There is evidence that Wolbachia in Lepidoptera have both parasitic and mutualistic relationships [17, 23, 37–40]. A well-known effect of Wolbachia in Lepidoptera is reproductive manipulation, including feminization, androcide, and cytoplasmic incompatibility [23, 37, 38]. One species of Wolbachia enhances the susceptibility of its lepidopteran host to baculovirus, rendering it a potential biological control agent against the agricultural pest Spodoptera exempta . In the current study, we demonstrate that Wolbachia infects a high proportion of Lepidoptera individuals and species, reflecting the significance of Wolbachia’s role in moth and butterfly biology .
Our global survey of published screenings of over 300 lepidopteran species found that 43 % of samples were infected with Wolbachia (Table 1), a substantially higher estimate than most previous localized reports: 16.2 % (n = 43 species ) from Panama, 35.2 % (n = 34 ) from the UK, 14.3 % (n = 21 ) from the US, 17 % (n = 24 ) from Uganda, 45 % (n = 49 ) from Japan, 43 % (n = 7 ) from Western Europe, and 52 % from India (n = 56 ).
We also expanded on previous studies by using a model-based analysis to estimate infection frequency across the order. Using a beta-binomial approach to estimate the distribution of prevalence across species , we estimated that the vast majority of Lepidoptera are infected with Wolbachia (around 80 %), a much greater frequency than has been estimated across arthropods as a whole [11, 12]. However, the mean prevalence in Lepidoptera reported here (27 % in the preferred model) does not significantly differ from the estimated prevalence in arthropods: the value reported by Hilgenboecker et al.  (25.3 %) falls within our confidence interval boundaries (24-31 %). The very high incidence that we estimate may reflect the opportunities for substantial horizontal transfer of Wolbachia in Lepidoptera. It must be pointed out that our results come from a data base consisting largely of true butterflies (Papilionoidea, comprising 83 % of the 660 Lepidoptera populations used). However, the results suggest that true butterflies might not be unrepresentative of the other major lepidopteran superfamilies, from which they do not differ significantly in our estimates (Fig. 3). This homogeneity in high frequency infection across groups makes sense in light of the fact that the majority of lepidopteran larvae feed on plant tissue and that adults obtain nectar from flowers or tree sap , if these food sources are possible ways to mediate cross-infection in different populations and species .
In arthropods, bacterial infection frequencies can be influenced by abiotic factors, such as geographical location and climate condition [43–46], or by biotic factors such as host genetic variation and competition with other endosymbionts on the same host [47, 48]. We found no correlation between Wolbachia infection frequency and phylogenetic relatedness of lepidopteran host groups. However, there was a significant correlation between infection frequency and host geography. Wolbachia infection of lepidopteran species has been known to differ between geographical regions , though not substantially, and other arthropod groups have shown no evidence of geographical variation [20, 22]. Our study revealed that Wolbachia infection frequencies in Lepidoptera are higher at lower absolute latitudes, suggesting that infection is greater in warmer climates (Fig. 2). This result is in agreement with the findings of Toju and Fukatsu , who observed high Wolbachia infection frequencies in weevils that were found in climates with higher temperature. However, Liu et al.  and Morag et al.  found that the higher infection frequencies occur in regions with moderate climatic conditions, as opposed to geographic regions with extreme climates. Furthermore, Sumi et al.  found no seasonal effect on Wolbachia infection in the butterfly Pseudozizeeria maha (Lycaenidae). Our results seem difficult to reconcile with experimental demonstrations that increasing temperatures reduce Wolbachia infection frequencies [50–53]. However, latitudinal gradients are also correlated with additional factors, such as species densities and their interactions , so there are other plausible explanations for the correlation between Wolbachia infection and absolute latitude, such as climatic factors. This study serves as a foundation for future research that can provide more insight into the factors that impact Wolbachia and its interactions with Lepidoptera.
We collated data from what we believe are all published sources of Wolbachia infection in Lepidoptera. We estimate that around 80 % of species and from a quarter to a third of individuals are infected at non-negligible frequency. We found no evidence that closely related taxonomic groups show similar infection levels, but we did find geographical variation, with closer locations tending to show similar infection levels. We also show that latitudinal gradient appears to be an important factor in infection level, with lower frequencies towards higher latitudes. Our study is the first to show such a latitudinal gradient in Wolbachia infection at such a broad taxonomic and geographic scale, and suggests that symbiont infection might be predictable from ecological variables.
We thank the McGuire Center for Lepidoptera and Biodiversity, Florida Museum of Natural History for their continued support. David Plokin and Jesse Breinholt provided useful suggestions. This study was supported by the University of Florida Research Opportunity Seed Fund (ROSF) and the National Science Foundation grant number DEB-1354585 to AYK.
- van Nieukerken EJ, Kaila L, Kitching IJ, Kristensen NP, Lees DC, Minet J, 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. Journal of the Lepidopterists’ Society. 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. doi:10.1007/s00435-011-0132-1.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 JW, Cummings MP et al. Evolutionary framework for Lepidoptera model systems. Genetics and Molecular Biology of Lepidoptera, Contemporary Topics in Entomology, eds Goldsmith M and Marec F. Gainesville: CRC Press; 2009.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. doi:10.1146/annurev.ento.51.110104.151029.PubMedView ArticleGoogle Scholar
- Werren JH, Baldo L, Clark ME. Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008;6(10):741–51. doi:10.1038/nrmicro1969.PubMedView 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. doi:10.1111/j.1574-6968.2008.01110.x.PubMed CentralPubMedView ArticleGoogle Scholar
- Zug R, Hammerstein P. Still a host of hosts for Wolbachia: Analysis of recent data suggests that 40 % of terrestrial arthropod species are infected. Plos One. 2012;7(6). doi:10.1371/journal.pone.0038544.Google Scholar
- Weinert L, Araujo-Jnr E, Ahmed M, Welch JJ. The incidence of bacterial endosymbionts in terrestrial arthropods. Proceedings of the Royal Society B: Biological Sciences. 2015; 282(1807). doi:10.1098/rspb.2015.0249.Google Scholar
- Stahlhut JK, Desjardins CA, Clark ME, Baldo L, Russell JA, Werren JH, et al. The mushroom habitat as an ecological arena for global exchange of Wolbachia. Mol Ecol. 2010;19(9):1940–52. doi:10.1111/j.1365-294X.2010.04572.x.PubMedView ArticleGoogle Scholar
- Bennett GM, Pantoja NA, O’Grady PM. Diversity and phylogenetic relationships of Wolbachia in Drosophila and other native Hawaiian insects. Fly. 2012;6(4):273–83. doi:10.4161/fly.21161.PubMed CentralPubMedView ArticleGoogle Scholar
- Nikoh N, Tanaka K, Shibata F, Kondo N, Hizume M, Shimada M, et al. Wolbachia genome integrated in an insect chromosome: evolution and fate of laterally transferred endosymbiont genes. Genome Res. 2008;18(2):272–80. doi:10.1101/gr.7144908.PubMed CentralPubMedView ArticleGoogle Scholar
- Narita S, Shimajiri Y, Nomura M. Strong cytoplasmic incompatibility and high vertical transmission rate can explain the high frequencies of Wolbachia infection in Japanese populations of Colias erate poliographus (Lepidoptera: Pieridae). Bull Entomol Res. 2009;99(4):385–91. doi:10.1017/s0007485308006469.PubMedView ArticleGoogle Scholar
- O'Neill S, Hoffmann A, Werren J. Influential passengers: Inherited microorganisms and arthropod reproduction. Influential passengers: Inherited microorganisms and arthropod reproduction. 1997.Google Scholar
- Werren JH, Windsor D, Guo LR. Distribution of Wolbachia among Neotropical Arthropods. Proceedings of the Royal Society B-Biological Sciences. 1995;1364:197–204.View ArticleGoogle Scholar
- Werren JH, Windsor DM. Wolbachia infection frequencies in insects: evidence of a global equilibrium? Proceedings of the Royal Society B-Biological Sciences. 2000;267(1450):1277–85. doi:10.1098/rspb.2000.1139.PubMed CentralView ArticleGoogle Scholar
- Duron O, Bouchon D, Boutin S, Bellamy L, Zhou L, Engelstaedter J et al. The diversity of reproductive parasites among arthropods: Wolbachia do not walk alone. Bmc Biology. 2008;6. doi:10.1186/1741-7007-6-27.Google 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. doi:10.1111/een.12002.View 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. doi:10.1073/pnas.0304068101.PubMed CentralPubMedView ArticleGoogle Scholar
- Charlat S, Hornett EA, Dyson EA, Ho PPY, Loc NT, Schilthuizen M, et al. Prevalence and penetrance variation of male-killing Wolbachia across Indo-Pacific populations of the butterfly Hypolimnas bolina. Mol Ecol. 2005;14(11):3525–30. doi:10.1111/j.1365-294X.2005.02678.x.PubMedView ArticleGoogle Scholar
- Delgado AM, Cook JM. Effects of a sex-ratio distorting endosymbiont on mtDNA variation in a global insect pest. BMC Evolutionary Biology. 2009;9:10 pp.-pp.Google Scholar
- West SA, Cook JM, Werren JH, Godfray HCJ. Wolbachia in two insect host-parasitoid communities. Mol Ecol. 1998;7(11):1457–65. doi:10.1046/j.1365-294x.1998.00467.x.PubMedView ArticleGoogle Scholar
- Jiggins FM, Hurst GDD, Schulenburg J, Majerus MEN. Two male-killing Wolbachia strains coexist within a population of the butterfly Acraea encedon. Heredity. 2001;86:161–6. doi:10.1046/j.1365-2540.2001.00804.x.PubMedView ArticleGoogle Scholar
- Tagami Y, Miura K. Distribution and prevalence of Wolbachia in Japanese populations of Lepidoptera. Insect Mol Biol. 2004;13(4):359–64. doi:10.1111/j.0962-1075.2004.00492.x.PubMedView ArticleGoogle Scholar
- Salunke BK, Salunkhe RC, Dhotre DP, Walujkar SA, Khandagale AB, Chaudhari R, et al. Determination of Wolbachia diversity in butterflies from western Ghats, India, by a multigene approach. Appl Environ Microbiol. 2012;78(12):4458–67. doi:10.1128/aem.07298-11.PubMed CentralPubMedView ArticleGoogle Scholar
- Edwards A. Likelihood expanded ed. Baltimore: Johns Hopkins University Press; 1992.Google Scholar
- Chessel D, Duffour A, Dray S. Ade4: Analysis of Ecological Data: Exploratory and Euclidean Methods in Environmental Sciences. http://cran.r-project.org/web/packages/ade4/ 2011.
- Hijmans R, Williams E, Vennes C. Geosphere: Spherical Trigonometry. http://cran.r-project.org/web/packages/geosphere/. 2014((version 1.3-8)).
- Regier JC, Mitter C, Zwick A, Bazinet AL, Cummings MP, Kawahara AY, 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
- Simoes PM, Mialdea G, Reiss D, Sagot MF, Charlat S. Wolbachia detection: an assessment of standard PCR protocols. Mol Ecol Resour. 2011;11(3):567–72. doi:10.1111/j.1755-0998.2010.02955.x.PubMedView ArticleGoogle Scholar
- Bloch D, Werdenberg N, Erhardt A. Pollination crisis in the butterfly-pollinated wild carnation dianthus carthusianorum? New Phytol. 2006;169(4):699–706. doi:10.1111/j.1469-8137.2005.01653.x.PubMedView ArticleGoogle Scholar
- Bloem KA, Bloem S, Carpenter JE. Impact of moth suppression/eradication programmes using the sterile insect technique or inherited sterility. Principles and Practice in Area-Wide Integrated Pest Management: Sterile Insect Technique; 2005.View ArticleGoogle Scholar
- Sasaki T, Ishikawa H. Wolbachia infections and cytoplasmic incompatibility in the almond moth and the mediterranean flour moth. Zoolog Sci. 1999;16(5):739–44. doi:10.2108/zsj.16.739.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. doi:10.1038/sj/hdy/6800077.PubMedView ArticleGoogle Scholar
- Kaiser W, Huguet E, Casas J, Commin C, Giron D. Plant green-island phenotype induced by leaf-miners is mediated by bacterial symbionts. Proceedings of the Royal Society B-Biological Sciences. 2010;277(1692):2311–9. doi:10.1098/rspb.2010.0214.PubMed 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. doi:10.1111/j.1461-0248.2012.01820.x.PubMedView ArticleGoogle Scholar
- Chapman R. Foraging and food choice in phytophagous insects. In Chemical Ecology (ed by H JD) Eolss Publishers, Oxford. 2009.Google Scholar
- Caspi-Fluger A, Inbar M, Mozes-Daube N, Katzir N, Portnoy V, Belausov E, et al. Horizontal transmission of the insect symbiont Rickettsia is plant-mediated. Proc R Soc B Biol Sci. 2012;279(1734):1791–6. doi:10.1098/rspb.2011.2095.View ArticleGoogle Scholar
- Tsuchida T, Koga R, Shibao H, Matsumoto T, Fukatsu T. Diversity and geographic distribution of secondary endosymbiotic bacteria in natural populations of the pea aphid. Acyrthosiphon pisum Molecular Ecology. 2002;11(10):2123–35. doi:10.1046/j.1365-294X.2002.01606.x.PubMedView ArticleGoogle Scholar
- Liu Y, Miao H, Hong XY. Distribution of the endosymbiotic bacterium Cardinium in Chinese populations of the carmine spider mite Tetranychus cinnabarinus (Acari : Tetranychidae). J Appl Entomol. 2006;130(9–10):523–9. doi:10.1111/j.1439-0418.2006.01112.x.View ArticleGoogle Scholar
- Toju H, Fukatsu T. Diversity and infection prevalence of endosymbionts in natural populations of the chestnut weevil: relevance of local climate and host plants. Mol Ecol. 2011;20(4):853–68. doi:10.1111/j.1365-294X.2010.04980.x.PubMedView ArticleGoogle Scholar
- Morag N, Klement E, Saroya Y, Lensky I, Gottlieb Y. Prevalence of the symbiont Cardinium in Culicoides (Diptera: Ceratopogonidae) vector species is associated with land surface temperature. Am Fasten J. 2012;26(10):4025–34. doi:10.1096/fj.12-210419.Google Scholar
- Chiel E, Gottlieb Y, Zchori-Fein E, Mozes-Daube N, Katzir N, Inbar M, et al. Biotype-dependent secondary symbiont communities in sympatric populations of Bemisia tabaci. Bull Entomol Res. 2007;97(4):407–13. doi:10.1017/s0007485307005159.PubMedView ArticleGoogle Scholar
- Gueguen G, Vavre F, Gnankine O, Peterschmitt M, Charif D, Chiel E, et al. Endosymbiont metacommunities, mtDNA diversity and the evolution of the Bemisia tabaci (Hemiptera: Aleyrodidae) species complex. Mol Ecol. 2010;19(19):4365–78. doi:10.1111/j.1365-294X.2010.04775.x.PubMedView ArticleGoogle Scholar
- Sumi T, Miura K, Miyatake T. No seasonal trend in infection of the pale grass blue butterfly, Zizeeria maha (Lepidoptera: Lycaenidae), by Wolbachia. Appl Entomol Zool. 2013;48(1):35–8. doi:10.1007/s13355-012-0152-4.View ArticleGoogle Scholar
- Van Opijnen T, Breeuwer JAJ. High temperatures eliminate Wolbachia, a cytoplasmic incompatibility inducing endosymbiont, from the two-spotted spider mite. Exp Appl Acarol. 1999;23(11):871–81. doi:10.1023/a:1006363604916.PubMedView ArticleGoogle Scholar
- Hurst GDD, Jiggins FM, Robinson SJW. What causes inefficient transmission of male-killing Wolbachia in Drosophila? Heredity. 2001;87:220–6. doi:10.1046/j.1365-2540.2001.00917.x.PubMedView ArticleGoogle Scholar
- Jia F-X, Yang M-S, Yang W-J, Wang J-J. Influence of Continuous High Temperature Conditions on Wolbachia Infection Frequency and the Fitness of Liposcelis tricolor (Psocoptera: Liposcelididae). Environ Entomol. 2009;38(5):1365–72.PubMedView ArticleGoogle Scholar
- Guruprasad NM, Mouton L, Puttaraju HP. Effect of Wolbachia infection and temperature variations on the fecundity of the Uzifly Exorista sorbillans (Diptera: Tachinidae). Symbiosis. 2011;54(3):151–8. doi:10.1007/s13199-011-0138-y.View ArticleGoogle Scholar
- Schemske DW, Mittelbach GG, Cornell HV, Sobel JM, Roy K. Is There a Latitudinal Gradient in the Importance of Biotic Interactions? Annual Review of Ecology Evolution and Systematics. Annual Review of Ecology Evolution and Systematics, 2009. p. 245–69.Google Scholar
This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly credited. 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.