Male-killing Wolbachia and mitochondrial selective sweep in a migratory African insect
© Graham and Wilson; licensee BioMed Central Ltd. 2012
Received: 15 May 2012
Accepted: 9 October 2012
Published: 15 October 2012
Skip to main content
© Graham and Wilson; licensee BioMed Central Ltd. 2012
Received: 15 May 2012
Accepted: 9 October 2012
Published: 15 October 2012
Numerous recent studies have shown that resident symbiotic microorganisms of insects play a fundamental role in host ecology and evolution. The lepidopteran pest, African armyworm (Spodoptera exempta), is a highly migratory and destructive species found throughout sub-Saharan Africa, that can experience eruptive outbreaks within the space of a single generation, making predicting population dynamics and pest control forecasting extremely difficult. Three strains of Wolbachia have recently been identified infecting this species in populations sampled from Tanzania. In this study, we examined the interaction between Wolbachia pipiensis infections and the co-inherited marker, mtDNA, within populations of armyworm, as a means to investigate the population biology and evolutionary history of Wolbachia and its host.
A Wolbachia-infected isofemale line was established in the laboratory. Phenotypic studies confirmed the strain wExe1 as a male-killer. Partial sequencing of the mitochondrial COI gene from 164 individual field-collected armyworm of known infection status revealed 17 different haplotypes. There was a strong association between Wolbachia infection status and mtDNA haplotype, with a single dominant haplotype, haplo1 (90.2% prevalence), harbouring the endosymbiont. All three Wolbachia strains were associated with this haplotype. This indicates that Wolbachia may be driving a selective sweep on armyworm haplotype diversity. Despite very strong biological and molecular evidence that the samples represent a single species (including from nuclear 28S gene markers), the 17 haplotypes did not fall into a monophyletic clade within the Spodoptera genus; with six haplotypes (2 each from 3 geographically separate populations) differing by >11% in their nucleotide sequence to the other eleven.
This study suggests that three strains of Wolbachia may be driving a selective sweep on armyworm haplotype diversity, and that based on COI sequence data, S. exempta is not a monophyletic group within the Spodoptera genus. This has clear implications for the use of mtDNA as neutral genetic markers in insects, and also demonstrates the impact of Wolbachia infections on host evolutionary genetics.
Recent studies into resident heritable symbiotic microorganisms have highlighted the central role they play in their insect host’s ecology and evolution [1, 2]. These symbionts can be classified as either obligate or facultative for host survival. Facultative symbionts are not essential for host development or reproduction, but their presence can impact upon host dynamics by manipulating host reproduction [3–5], as well as increasing host survival or fecundity [6–10], whereas removal of obligate symbionts results in the death of its host [11, 12]. One successful group of such symbiotic microorganisms belongs to the genus Wolbachia. These are intracellular, maternally-inherited bacteria belonging to the alpha-Proteobacteria group Rickettsia. They are among the most successful genus of bacteria, found in filarial nematodes, crustaceans, arachnids, and are estimated to occur in 20 - 70% of all insect species [13, 14].
Wolbachia have the ability to induce a number of reproductive manipulations of their hosts, such as feminisation of genetic males, induction of sperm-egg incompatibilities, thelytokous parthenogenesis, and male-killing [1, 15]. Through these processes, Wolbachia provides infected host-females with a relative reproductive advantage over uninfected females . Each of these mechanisms increases the symbiont-infection frequency by ensuring that the number of infected daughters produced by an infected female is greater than the average production of daughters per female. A single insect or population may be infected with more than one symbiotic microorganism and different populations may have different infection statuses [17, 18].
Both Wolbachia and host mitochondria are maternally transmitted and subsequently can be co-inherited by the offspring. The high mutation rate of mitochondrial DNA makes it a valuable evolutionary marker when compared to corresponding nuclear DNA . With this in mind, it also has the advantage of very low recombination, resulting in the whole mtDNA genome having the same genealogical history . It is well documented that symbionts such as Wolbachia can impact upon mtDNA diversity within host populations [21–25]. Indeed, the strong linkage between the two co-inherited markers makes them ideal candidates for the investigation of symbiont invasion-history and symbiont impact upon host genetics [21, 26, 27]. Previous studies have documented the role of Wolbachia in driving dramatic changes within host populations, due to the phenotypes induced by the symbionts . It is now widely accepted that endosymbiont screening and analysis should take place before any attempt to explain mtDNA patterns in terms of host ecology and evolution [20, 28].
The larval stage of the African armyworm, Spodoptera exempta (Lepidoptera: Noctuidae), is one of the most devastating crop pests in Africa, feeding upon many of the staple food crops such as maize, wheat, sorghum, millet, rice and pasture grasses. Most outbreaks occur on the eastern half of sub-Saharan Africa, as far north as Sudan and as far south as South Africa . The adult moths are highly migratory, often flying hundreds of kilometres over consecutive nights , with moth movements largely determined by the seasonal progression of the inter-tropical convergence zone (ITCZ). Therefore, it is believed that early-season armyworm outbreaks in central Tanzania essentially act as “source” populations for moths that will subsequently migrate to northern Tanzania and further northwards towards the horn of Africa .
A baculovirus, SpexNPV, is known to be present in field populations of armyworms causing larval mortality, which can have significant impact upon host population dynamics during sporadic natural epizootics. In a recent study, we found that armyworm harboured three strains of Wolbachia, designated wExe1, wExe2 and wExe3. Based on MLST classification, one of the strains, wExe1, was 100% identical to ST-125, which is a male-killing phenotype found in the nymphalid butterfly, Hypolimnas bolina. The other two strains of Wolbachia are new sequence-types, assigned ST-222 (wExe2 – clade B) and ST-223 (wExe3 – clade A). Interestingly, both laboratory and field data showed that infection with wExe1 strain increased host susceptibility to a baculovirus, which could have clear implications for the evolutionary history of the host and pathogen .
In the current study, our aim was to address the following questions: (i) How robust is the male-killing phenotype and efficiency of wExe1 vertical transmission? (ii) In the field, is global Wolbachia infection associated with particular host mitochondrial genotypes? (iii) Is there evidence of a selective sweep within the host population? To answer these questions, we established laboratory cultures of wExe1-infected armyworm to assess the infection phenotype. We obtained partial sequences of Wolbachia genes and the mitochondrial cytochrome oxidase I (COI) gene from field-collected S. exempta larvae and adults collected from pheromone traps. We examined the diversity of mitochondrial haplotypes, and analysed mtDNA variation to explore the potential association with different Wolbachia infection statuses.
Sex ratio and hatch ratio of insect lines reared under laboratory conditions, uninfected and infected with Wolbachia wExe1 strain
Egg hatch (%)
Egg hatch (%)
Egg hatch (%)
We explored the evolutionary history of Wolbachia pipientis infections in Tanzanian populations of the African armyworm, S. exempta. Three strains of Wolbachia have recently been detected at an infection prevalence of approximately 12% . Based on previous MLST analysis (using the five genes ftsZ, coxA, fbpA, hcpA and gatB) , one of these strains, wExe1, is 100% identical to ST-125, which is a male-killing phenotype found in the Blue Moon Butterfly, Hypolimnas bolina and, based on nucleotide identity for wsp and ftsZ genes, 100% identical to the Wolbachia isolate found in Tanzanian populations of the Common Acraea, Acraea encedon. This indicates that an infection-transmission event may have occurred between host species at some point in evolutionary history, possibly via a shared generalist parasitoid [34–37]. Indeed, a recent study has shown that plant-mediated horizontal transmission may also be possible in some insect-endosymbiont systems . The phenotypes of strains wExe2 and wExe3 are as yet undetermined. However, due to the absence of Wolbachia in male adult moths, it is possible that all three strains of Wolbachia display the male-killing phenotype. Male-killing Wolbachia have been previously identified in both Wolbachia Clade A  and Clade B e.g. [21, 32].
Male-killing is a widespread phenotype in lepidopteran-infecting Wolbachia and does appear to explain the sex ratio bias and egg-hatch rate in our infected lab lines. However, we cannot exclude the possibility of a more complicated scenario in the field, as Wolbachia is not the only cause of sex ratio distortion to be discovered in insects . Studies on a range of species have revealed a number of insect endosymbionts that can cause similar phenotypic effects to Wolbachia[41–43], including the generation of complex fitness interactions between infection-agent and host sex chromosomes, leading to sex-biased mortality [44, 45]; DGGE analysis suggests that no secondary endosymbionts occur in the armyworm , but this remains an interesting avenue for more detailed future study.
The lack of geographical structuring (of both Wolbachia infection status and host mitochondrial genetics) throughout Tanzania is most likely due to the highly migratory nature of the armyworm, thereby preventing genetic differentiations from establishing. This failure to detect any spatial structuring in the armyworm population is consistent with early genetic analyses using iso-enzymes over a much wider East African geographical range . The very large effective population size of this outbreaking species is also likely to be factor. The use of microsatellite markers is currently being examined as a technique allowing more sensitive exploration of the armyworm population structure .
Of particular interest was that all three strains of Wolbachia were found to be associated with the same host haplotype, haplo1. This haplotype was found in 90.2% of samples tested, suggestive of a selective sweep of this haplotype. No other host haplotype was observed to harbour Wolbachia, indicative of very low levels of horizontal symbiont transmission between individuals within a population. The high 28S nuclear gene variability, but selection for a single mtDNA haplotype, provides evidence that COI selection is being driven, and that this pattern is not due to demographic processes alone, such as a genetic bottleneck. We propose that this is due to the presence of multiple strains of Wolbachia (including the male-killer, wExe1). Two hypotheses have been previously proposed for why male-killers reduce mtDNA diversity, causing a sweep . Firstly, the initial invasion of Wolbachia, and the subsequent selection for the beneficial parasite genes will result in selective sweeps of the co-inherited mtDNA molecules. Secondly, once the Wolbachia invasion has reached equilibrium within the host population, the effective population size of mtDNA will be reduced because mtDNA mutations in uninfected females will tend to be lost . We propose that such a marked haplo1 sweep has occurred due to three separate Wolbachia-strain invasion events associated with mtDNA haplo1, each one enhancing the haplo1 sweep further. It is known that invasion events can cause selection for a particular haplotype, and therefore if three such events occurred in the same haplotype, a haplotype-skew towards 90% prevalence may not be unexpected. However, it must be stated that at this stage of our research, we can only speculate on this possibility.
So why is there such a large haplotype skew (90%) but only a moderate prevalence (12%) of Wolbachia? Several hypotheses have been proposed for maintaining Wolbachia prevalence equilibrium within host populations, including meta-population structure, mating preference, and environmentally-occurring antibiotics [39, 49]. All such spontaneous “cure” events would increase the ratio of uninfected individuals with the same haplotype to infected individuals. We propose two main hypotheses that may explain Wolbachia stability in our armyworm system. Firstly, inefficient transmission of male-killing Wolbachia. This occurs when an infected female produces offspring that are both infected and uninfected, resulting in individuals of both infection states having the same mtDNA profiles. Results from our laboratory cultures provide evidence for inefficient transmission in the S. exempta system, and this phenomenon would certainly help explain the skewed uninfected-haplo1 distribution. Wolbachia infection intensity within an individual may impact greatly upon the transmission efficiency between mother and offspring, whereby a reduction in infection intensity may lead to a decrease in vertical transmission efficiency and consequently loss of infection . In contrast, excessive infection intensity may result in pathology, resulting in negative effects upon host fitness, as seen in the Drosophilla melanogaster Wolbachia strain wMelpop. A number of previous studies have identified differences in endosymbiont intensity and infection-load within individual hosts [50, 52, 53]. We are currently developing quantitative assays to assess infection intensity and strain identification for the armyworm-Wolbachia system.
Secondly, Wolbachia-infected larvae may be more susceptible to an endemic armyworm baculovirus . This would keep a check on Wolbachia prevalence in populations by causing greater mortality in Wolbachia-infected larvae than non-infected haplo1 larvae, thereby further exaggerating the haplo1 uninfected:infected ratio. We propose that a combination of these two mechanisms occurs in the field and contributes to maintaining stability of Wolbachia and preserving the mtDNA haplo1 skew within populations of armyworm. In addition, inter-strain Wolbachia stability may exist if the three strains are adapted to different 'environments'. Armyworm experience huge environmental variations over the course of a season e.g. host density, temperature, food-plant changes, disease challenge, etc. The three Wolbachia strains may be differentially adapted to hosts experiencing different conditions e.g. those that do better under low or high density conditions, or they interact differently with baculovirus strains. Indeed, in field populations there was no correlation between wExe3 prevalence and viral mortality, indicating that wExe3 may not increase virus susceptibility in armyworms (unlike wExe1 and wExe2) .
S. exempta appears to be paraphyletic at the mtDNA level within the genus Spodoptera, with >11% nucleotide divergence amongst haplotypes within what observational and nuclear-gene 28S molecular studies indicate is one species. As previously discussed, mtDNA is a popular method of identifying species, and assessing biodiversity via barcoding protocols [54, 55]. Historically these have found intra-specific mtDNA diversity to be very low (typically <1%), and inter-specific diversity to be higher (>2%). However, an increasing number of studies are finding high levels of mtDNA diversity within classical species [25, 56], leading to the hypothesis that cryptic species and races are present in greater numbers than previously thought [57, 58]. In addition to this, endosymbionts are capable of driving mtDNA introgression from other neighbouring species, thereby confounding the effect of paraphyly within a host species [20, 59, 60], such as that observed previously in Drosophila, Acraea butterflies  and Stomoxys flies . This hypothesis may explain the paraphyly patterns observed in our S. exempta haplotype phylogeny.
All three strains of a Wolbachia infecting African armyworm S. exempta were associated with a single host haplotype, haplo1, which comprised 90.2% of the total samples tested. This study suggests that Wolbachia is driving a selective sweep for this particular haplotype, and that based on COI diversity, S. exempta is not a monophyletic group within the Spodoptera genus. This study supports previous research highlighting clear implications for the use of mtDNA as neutral genetic markers in insects.
Larval outbreaks were sampled in Tanzania between 2007 and 2011 as previously reported . Briefly, a minimum of thirty live larvae were collected in individual microtubes and stored in 100% ethanol for use in subsequent laboratory analysis. Adult males were caught throughout the season using a network of pheromone traps located across Tanzania. Trained trap operators collected specimens daily, and stored moths in 100% ethanol. The laboratory Spodoptera exempta culture was established from larvae collected in central Tanzania in January 2011, and high numbers were maintained at each generation to reduce inbreeding. All larval lines were maintained on a standard wheatgerm-based artificial diet  at a constant temperature of 25°C under a 12 hour light/dark cycle. Wolbachia-infected and Wolbachia-free insect lines were maintained in the same facility to ensure identical breeding conditions.
DNA extractions were performed using the AllPrep DNA/RNA Mini Kit (Qiagen, Crawley, UK) according to the manufacturer's protocol. Larvae were screened for Wolbachia infection by PCR amplifying the Wolbachia-specific wsp and ftsZ genes, using the primers wsp-81F/wsp-691R  and ftsZfl/ftsZrl, respectively. Reaction mixtures (50 μl) contained PCR buffers (10 mM Tris–HCl pH 8.3 at 25°C; 50 mM KCl; 1.5 mM MgCl; 0.001% gelatin), 5 mM each of dATP, dTTP, dCTP and dGTP, 10 mM of the relevant primers, 1 unit Taq polymerase (Qiagen, Crawley, UK) and approximately 10 ng DNA template. PCR was carried out in a Techne TC-512 thermal cycler (Bibby Scientific Ltd., Stone, UK), under the following reaction conditions: (i) 94°C for 5 min, 1 cycle; (ii) 94°C for 30s, 52°C for 30s, 7°C for 30s, 40 cycles; and (iii) 72°C for 5 min, 1 cycle. Wolbachia MLST analysis was undertaken using PCR protocols for amplification of the five reported Wolbachia MLST genes (ftsZ, coxA, fbpA, hcpA and gatB) as described elsewhere . PCR reaction-conditions were as above. All PCR products were run on a 1% agarose gel and visualized using ethidium bromide staining. PCR amplicons were purified (PCR Purification Kit, Qiagen), and sequenced (Source Bioscience, UK). The sequence data were analyzed against the Wolbachia MLST database (http://pubmlst.org/Wolbachia/).
To study armyworm haplotype diversity and distribution, amplification of the S. exempta mitochondrial cytochrome oxidase I (COI) gene was undertaken, using universal primers LCO-1490 and HCO-1298 , which yielded approximately 640bp amplicons. A sample comprising at least four larvae per 2007–08 population (n = 121), and then a random number of larvae and adults from the field season 2009–10 (n = 43) was investigated. Polymerase chain reaction (PCR) mixtures (50 μL) contained PCR buffer (10 mm Tris–HCl pH 8.3 at 25°C, 50 mm KCl, 1.5 mm MgCl, 0.001% gelatin), 10 μm each of dATP, dTTP, dCTP and dGTP, 0.1 μm of each primer, 1 unit Taq polymerase (Qiagen, Crawley, UK) and approximately 10 ng genomic DNA template. PCR reactions were as above. All PCR products were run on a 1% agarose gel and visualized using ethidium bromide staining. Nuclear DNA diversity was studied by amplifying a partial segment of the 28S ribosomal gene, using primers 28SFor 5' AAA GAT CGA ATG GGG AGA TTC ATC 3' and 28SRev 5' CGT CCT ACT AGG GGA GAA GTG CAC 3' to yield an approximately 150 bp product. PCR reactions were as above. All PCR products were run on a 1% agarose gel and visualized using ethidium bromide staining.
A proportion of all Wolbachia-infected armyworm lines used were treated with 0.03% tetracycline (10 mg/ml) to generate uninfected fly lines. This was achieved by feeding the armyworm larval stage, insuring ingestion of sufficient antibiotic to cure the Wolbachia infection, as previously described . Following the tetracycline treatment, armyworms were maintained for a generation to recover before experimental use. All tetracycline-treated lines were tested with Wolbachia-specific PCR to test for Wolbachia infection. Only clean uninfected lines were used in subsequent experiments.
PCR products were directly sequenced at Source BioScience UK Ltd. (Nottingham, UK). Sequences were viewed using BioEdit  and edited to remove universal primer regions. Preliminary identifications against previously published sequences were provided by BLAST . Sequence alignment was performed with ClustalW 1.8 using default parameters . Maximum-Likelihood phylogenetic analysis was performed having selected the most appropriate model based on lowest AIC score (MEGA 5.05 ). Sequences obtained in this study were deposited in GenBank [Genbank JQ315120-JQ315136]. Tajima’s D test, Fu and Li's D* and F* tests, haplotype diversity and nucleotide diversity were calculated using DnaSP vers. 5.10 . Median-joining haplotype networks were drawn using Network (version 126.96.36.199; ).
This work was supported by a grant from the Sustainable Agriculture Research for International Development (SARID) programme, funded by the United Kingdom’s Biotechnology and Biological Sciences Research Council (BBSRC) and the Department for International Development (DfID). The views expressed are not necessarily those of DfID or BBSRC. Additional support was provided by Lancaster University, U.K. We are extremely grateful for the support and assistance of the Tanzanian Ministry of Agriculture, Food and Cooperatives and to Ms Yamini Tummala and Mr Phill Nott for technical support. We are especially grateful to Mr Wilfred Mushobozi for coordinating field activities in Tanzania and Mr David Grzywacz, Dr Alan Shirras, Professor Jenny Cory, Dr Glenn Rhodes and Dr Kamal Ibrahim for useful discussion.
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.