- Open Access
Influence of Wolbachiaon host gene expression in an obligatory symbiosis
© Kremer et al; licensee BioMed Central Ltd. 2012
- Published: 18 January 2012
Wolbachia are intracellular bacteria known to be facultative reproductive parasites of numerous arthropod hosts. Apart from these reproductive manipulations, recent findings indicate that Wolbachia may also modify the host’s physiology, notably its immune function. In the parasitoid wasp, Asobara tabida, Wolbachia is necessary for oogenesis completion, and aposymbiotic females are unable to produce viable offspring. The absence of egg production is also associated with an increase in programmed cell death in the ovaries of aposymbiotic females, suggesting that a mechanism that ensures the maintenance of Wolbachia in the wasp could also be responsible for this dependence. In order to decipher the general mechanisms underlying host-Wolbachia interactions and the origin of the dependence, we developed transcriptomic approaches to compare gene expression in symbiotic and aposymbiotic individuals.
As no genetic data were available on A. tabida, we constructed several Expressed Sequence Tags (EST) libraries, and obtained 12,551 unigenes from this species. Gene expression was compared between symbiotic and aposymbiotic ovaries through in silico analysis and in vitro subtraction (SSH). As pleiotropic functions involved in immunity and development could play a major role in the establishment of dependence, the expression of genes involved in oogenesis, programmed cell death (PCD) and immunity (broad sense) was analyzed by quantitative RT-PCR. We showed that Wolbachia might interfere with these numerous biological processes, in particular some related to oxidative stress regulation. We also showed that Wolbachia may interact with immune gene expression to ensure its persistence within the host.
This study allowed us to constitute the first major dataset of the transcriptome of A. tabida, a species that is a model system for both host/Wolbachia and host/parasitoid interactions. More specifically, our results highlighted that symbiont infection may interfere with numerous pivotal processes at the individual level, suggesting that the impact of Wolbachia should also be investigated beyond reproductive manipulations.
- Gene Ontology
- Suppression Subtractive Hybridization
- Parasitoid Wasp
- Cytoplasmic Incompatibility
- Wolbachia Infection
Symbiotic communities of eukaryotic organisms are known to influence host developmental programs  and also to shape immune response against pathogens . Interestingly, some genes/pathways (e.g. programmed cell death) have a pleiotropic role in immunity and development, and could play a major role in the maintenance of a specific bacterial community. For instance, the homeobox gene Caudal is involved in the formation of the antero-posterior body axis of Drosophila, but also in the regulation of the commensal gut microbiota . In the squid-vibrio association, it has recently been shown that the regulation of a peptidoglycan recognition protein (PGRP), classically involved in innate immunity, plays a role in the activation of the apoptotic process initiating the morphogenetic changes of the symbiont-harboring organ . The generality of the interplay between immunity and development during symbiosis is currently unknown.
Wolbachia (Anaplasmataceae) is among the most abundant intracellular bacteria. It infects both arthropods and nematodes, and is known to be a master manipulator of host biology . Wolbachia is generally a facultative reproductive parasite in arthropods, and invades the host population by inducing cytoplasmic incompatibility, male-killing, feminization or thelytokous parthenogenesis .
Another extended phenotype due to the presence of Wolbachia is observed in the parasitoid wasp Asobara tabida (Hymenoptera: Braconidae), in which aposymbiotic females exhibit a strong developmental defect. Surprisingly, Wolbachia has become necessary for egg production in this wasp, since aposymbiotic females are unable to produce viable offspring . Interestingly, A. tabida is the only member of the genus Asobara to be dependent on Wolbachia for oogenesis, which suggests that the dependence has evolved recently, and makes it possible to study the molecular mechanisms underlying this transition. In addition, polymorphism of the ovarian phenotype is observed in natural populations after the elimination of Wolbachia: some aposymbiotic females do not produce eggs, whereas others produce a few eggs that die prematurely [7, 8]. This polymorphism constitutes a tool to better understand the influence of these molecular mechanisms on the severity of the ovarian phenotype and on the evolution of dependence.
At a mechanistic level, cytological analysis of the ovarian phenotype has begun to shed light on the mechanisms underlying dependence in A. tabida. Indeed, eliminating Wolbachia triggers programmed cell death (PCD) in the egg chambers within the ovaries of A. tabida females . As egg production is tightly controlled by two main apoptotic checkpoints during oogenesis , the deregulation of PCD in aposymbiotic wasps must result in female inability to complete oogenesis. Because PCD is frequently involved in infection processes by bacterial pathogens , it has been hypothesised that a mechanism underlying the maintenance of Wolbachia at the individual level may have given rise to the evolution of dependence through its pleiotropic role in immunity and development .
This hypothesis is supported by recent findings showing that consequences of Wolbachia infection in insects may extend far beyond the classical effect on reproduction, by impacting host physiology and immunity. Wolbachia could play a role as a nutritional mutualist, by influencing iron utilization by its Drosophila hosts [13, 14]. Wolbachia infection has also been shown to generate oxidative stress in one Aedes aegypti cell line, which reacts by the over-expression of host antioxidant genes . Interestingly, Reactive Oxygen Species (ROS) are known to play a major role in immunity as a first line of defence  but also as a mechanism insuring microbe homeostasis . Finally, Wolbachia is known to confer resistance against RNA viral infection in D. melanogaster and D. simulans [18, 19], and against various pathogens in the mosquito A. aegypti, notably by priming the innate immune system [20, 21].
To summarize, increasing evidence is emerging on the phenotypic effects of Wolbachia infection on host physiology and immunity [18, 19, 22]. However, few studies have attempted to describe the molecular mechanisms underlying these phenotypic effects in natural systems [20, 21, 23, 24]. The objective of this paper is to clarify the effect of Wolbachia on gene expression in a particular symbiotic association in which Wolbachia affects developmental processes, through its effect on wasp oogenesis. For that purpose, we used both global and dedicated transcriptomic approaches.
Even though A. tabida is a model system in host/parasitoid and host/Wolbachia interactions, no genetic data were available for this parasitoid wasp. Thus, the first aim of this study was to build a reference transcriptome based on several tissues (ovaries, whole females) and physiological conditions (symbiosis, immune challenge). By sequencing 10 cDNA libraries (one of which is a normalized library), we provide here the first large-scale, genetic information on this wasp. The second aim of the study was to better understand how dependence arose in this particular species by deciphering the molecular mechanisms underlying this evolutionary transition. An overview of functions that could be differentially expressed in response to symbiosis was outlined through in silico analyses on ovaries EST libraries (Gene Ontology-based bioinformatics) and in vitro subtractions (Suppressive Subtraction Hybridizations). Then, we focused on candidate genes involved in immunity (broad sense), programmed cell death and oogenesis; functions which could play a major role in the control of ovarian phenotype through pleiotropy. Using quantitative real-time PCR, we thus characterized the effect of symbiosis on host gene expression in both males and females, in two populations exhibiting extreme ovarian phenotypes.
Asobara tabida (Hymenoptera: Braconidae) is a solitary endoparasitoid laying its eggs into the first or second instar larvae of Drosophila species. After Drosophila pupation, the parasitoid becomes an ectoparasite, and consumes its host before it itself pupates prior to emerging.
A. tabida is naturally infected by three strains of the intracellular bacterium Wolbachia (wAtab1, wAtab2 and wAtab3): wAtab1 and wAtab2 induce cytoplasmic incompatibility, and only wAtab3 is required for oogenesis completion [6, 25].
Polymorphism of ovarian phenotype in populations
After Wolbachia removal, the ovarian phenotype displays a high level of intra-species variation: whereas uninfected females of the Pi strain (Pierrefeu, France) produce no eggs, uninfected females of the NA strain (Saanich, Canada) produce a small number of aborting eggs . In this study, we used the NA strain and a Pi-derived strain (Pi3). Pi3 was obtained by moderate antibiotic treatment, and contains only the obligatory Wolbachia strain wAtab3 . The lines are stable, and have been maintained by regular sib-matings without antibiotic treatment for about 100 generations. When comparing NA and Pi3 aposymbiotic individuals, three things must be kept in mind. (i) Any differences observed may be explained by the host genotype, whether they are directly linked to the ovarian phenotype or not. (ii) Because NA is triply infected whereas Pi3 is singly infected, differences could also be due to the presence or absence of wAtab1 and wAtab2. (iii) NA and Pi3 symbiotic individuals have differing bacterial community compositions due to the moderate antibiotic treatment of Pi3 .
Wasps were allowed to parasite Wolbachia-free D. melanogaster. Insects were reared on axenic medium  and maintained under controlled conditions (climate chambers at 21°C, 70% relative humidity and cycle LD 12:12). Young adults (0-1 day old) were collected and anesthetized on ice before being dissected in a drop of PBS and/or stored until use at -80°C.
Because we were interested in determining the effect of symbiosis, we performed antibiotic treatments to produce Wolbachia-free (i.e. aposymbiotic) wasps. Even though antibiotics could also affect host gene expression directly (e.g. cytotoxicity, modification of mitochondrial metabolism) or indirectly (e.g. change in gut microflora), antibiotic treatment is the only efficient method to eliminate Wolbachia from A. tabida. Aposymbiotic females are sterile, and so it is impossible to establish and maintain aposymbiotic lines. Hence, antibiotic treatments had to be administered just before the experiment to obtain aposymbiotic wasps, as described in . Briefly, rifampicin 2% (Hoechst, Germany) was added to the axenic nutritive medium to reach a final concentration of 2 mg/g of standard diet. Seventy D. melanogaster eggs were deposited in this medium, and allowed to be parasitized by three female wasps. The developing Drosophila thus transferred the antibiotic to each of the endoparasitoid wasp larvae, rendering them aposymbiotic. As a control, the same procedure was performed without the antibiotic treatment.
Because we were interested in identifying immunity-related genes, we performed a challenge by the intracellular bacteria Salmonella typhimurium (strain 12023G, Grenoble) to enhance the immune response of A. tabida (Pi3 strain). Bacteria were prepared from a 2 h-culture initially started with a 1/10 dilution of an overnight culture (LB + ampicillin, 37°C, 190 rpm). Bacteria were rinsed twice and concentrated in 1 mL of fresh LB medium. Immune challenge was performed by injecting 13.2 nL of the mother solution (corresponding to 1.8x105 bacteria) in the thorax of young (0-1 day old) females (Nanoject II injector, Drummond, Broomall, PA). As a control, 13.2 nL of fresh LB medium was injected as described above. Individuals were collected 3h, 6h and 12h after challenge (or LB injection), and stored until use at -80°C.
Constitution of a reference transcriptome and comparison of gene expression profiles between EST libraries
Preparation of a normalized library for cDNA sequences acquisition
In order to build a transcriptome of reference for A.tabida, we constructed a normalized library (N) based on both whole females (mix of complex tissues) and ovaries (organ of interest), in various physiological conditions (with or without symbionts/pathogens). To limit host genetic variability, only the Pi3 strain was used for the library preparation. The normalized library was constructed by Evrogen (Moscow, Russia) from an equimolar proportion of total RNA prepared from aposymbiotic ovaries, symbiotic ovaries, and 3h-, 6h-, 12h-challenged symbiotic females. Total RNA samples were used for ds cDNA synthesis using the SMART approach . SMART-prepared, amplified cDNA was then normalized using the DSN normalization method . Normalization included cDNA denaturing/re-association, treatment by duplex-specific nuclease (DSN)  and amplification of normalized fraction by PCR. Normalized cDNA was purified using QIAquick PCR Purification Kit (Qiagen, Alameda, CA), digested with restriction enzyme Sfi1, purified (BD Chroma Spin - 1000 column), and ligated into pAL 17.3 vector (Evrogen) for Escherichia coli transformation.
Preparation of EST libraries for in silicocomparisons between symbiotic and aposymbiotic ovaries
In order to increase the number of transcripts from the ovaries and to determine the influence of symbiosis on host gene expression, we constructed EST libraries on aposymbiotic (OA1 and OA2, the quality of the OA2 library being slightly lower) and symbiotic (OS) ovaries (Pi strain). Total RNA was extracted from a large number of ovaries (nOA=196, nOS=120) as described in , and treated with DNAse (TurboDNase, Ambion, Applied Biosystems, Austin, TX), following the Manufacturer’s instructions. Tissue libraries were prepared using Creator SMART cDNA Library Construction kit (Clontech/BD biosciences, PaloAlto, CA), following the Manufacturer’s instructions. cDNA was digested by Sfi1, purified (BD Chroma Spin – 400 column), and ligated into pDNRlib vector for E. coli transformation.
Preparation of Suppression Subtractive Hybridizations (SSH) libraries for in vitrocomparisons
Because in silico comparisons of EST libraries can be limited by the depth coverage, we also used a complementary technique to compare gene expression by directly screening differentially-expressed transcripts through SSH.
In order to better understand the influence of ovarian phenotype, we performed SSHs between aposymbiotic (A) and symbiotic (S) ovaries in two populations exhibiting extreme phenotypes (Pi3: no eggs in aposymbiotic ovaries, NA: few abnormal eggs in aposymbiotic ovaries). Total RNA was extracted from a large number of ovaries [nA=373 and nS=458 for SSHs-1 A-S (Pi strain, distal part of ovaries), nA=nS=200 for SSHs-2 A-S (NA strain, whole ovaries)] and treated with DNAse (TurboDNase, Ambion, Applied Biosystems, Austin, TX), following the Manufacturer’s instructions. Amplified ds cDNA was prepared using a SMART approach . SMART Oligo II oligonucleotide (Clontech/BD biosciences, PaloAlto, CA) and CDS primer were used for first-strand cDNA synthesis. SMART-amplified cDNA samples were further digested by RsaI endonuclease. Subtractive hybridizations were performed using the SSH method in both directions (Aposymbiotic vs. Symbiotic A/S and vice-versa S/A) as described in [32, 33] using the PCR-Select cDNA Subtraction Kit (Clontech/BD biosciences, PaloAlto, CA). In order to reduce the number of false-positive clones in the SSH-generated libraries, the MOS procedure (Mirror Orientation Selection) was performed by Evrogen (Moscow, Russia) for SSH2s A-S, as described in . Purified subtracted cDNAs from SSH1s A-S were cloned into the PCR 2.1 TOPO vector (Invitrogen, Cergy-Pontoise, France) and used for E.coli transformation. 137 and 72 clones (SSH1-A/S and SSH1-S/A), respectively, were selected for further confirmation. Purified cDNA from SSH2s A-S were cloned into the pAL16 vector (Evrogen) and used for E. coli transformation. 480 clones for each subtraction were selected for further confirmation. PCR-amplified inserts from clones representing differentially-expressed gene products were confirmed by differential hybridization using either DIG-labeled (SSH1s A-S; DIG high prime DNA labeling and detection starter kit, Roche, Meylan, France) or P-32-labeled (SSH2s A-S), subtracted cDNA probes.
Finally, in order to characterize genes responding to bacterial challenge, we performed SSHs between extracts from whole females, challenged or not challenged by S. typhimurium (SSHs C-NC, nC=nNC=40 females), see above for bacterial challenge procedure. The preparation of these SSHs has been performed by Evrogen (Moscow, Russia) with the same procedure as for SSH2s A-S.
EST sequencing, data processing and analysis
Clustering and assembly of the ESTs were performed with TGICL  to obtain putative unique transcripts (unigenes) composed of contiguous ESTs (contigs) and unique ESTs (singletons). To do this, a pairwise comparison was first performed using a modified version of megablast (minimum similarity 94%). Clustering was done with tclust, which proceeds by a transitive approach (minimum overlap: 60 bp at 20 bp maximum of the end of the sequence). Assembly was done with CAP3 (minimum similarity 94%).
To detect unigene similarities with other species, several blasts (with high cut-off e-values) were performed against the following databases: NCBI nr (blastx (release: 1 March 2011); e-value < 5, HSP length > 33aa), Refseq genomic database (blastn, e-value < 10), Unigene division Arthropods (tblastx, #8 Aedes aegypti, #37 Anopheles gambiae, #3 Apis mellifera, #3 Bombyx mori, #53 Drosophila melanogaster, #9 Tribolium castaneum; e-value < 5), Nasonia vitripennis Nvit OGS_v1.0 (CDS predicted by Gnomon (NCBI)) and Wolbachia sequences from Genbank (blastn (release 164); e-value < e-20). Gene Ontology annotation was carried out using Blast2go software . During the first step (mapping), a pool of candidate GO terms was obtained for each unigene by retrieving GO terms associated with the hits obtained after a blastx search against NCBI nr. During the second step (annotation), reliable GO terms were selected from the pool of candidate GO terms by applying the Score Function (SF) of Blast2go with permissive annotation parameters (EC_weight=1, e-value_filter=0.1, GO_weight=5, HSP/hit coverage cut-off=0%). In the third step of the annotation procedure, the pool of GO terms selected during the annotation step was merged with GO terms associated with Interpro domain (Interpro predictions based on the longest ORF). Finally, the Annex augmentation step was run to modulate the annotation by adding GO terms derived from implicit relationships between GO terms .
In order to extract the biological processes and molecular functions statistically over-represented in aposymbiotic libraries, we performed a hyper-geometrical test between GO terms from the aposymbiotic libraries (OA1 and OA2) and those from the OS library, which corresponds to natural physiological conditions. The p-values were then adjusted using Bonferroni’s correction. To perform a functional enrichment analysis of the unigenes extracted from the SSH, we used the FatiGO web tool  on the OS library. With respect to the GO analysis, levels 3 and 6 were chosen to describe biological processes, and level 4 was chosen to describe molecular functions.
Gene expression measurement by quantitative RT-PCR (qRT-PCR)
We sought to determine the effect of symbiosis on the expression of a set of candidate genes involved in immunity, programmed cell death and oogenesis. For that purpose, we first compared gene expression between symbiotic and aposymbiotic samples, in ovaries (to characterize the dependence phenotype induced by Wolbachia) and then in males (to provide additional information concerning the specificity of the process). In order to limit the influence of the presence of eggs in symbiotic vs. aposymbiotic ovaries of the Pi3 strain, only the distal part that does not contain eggs (DPOv) was dissected in that strain. Because the dependence phenotype is determined by the host genotype , we compared gene expression between two populations exhibiting extreme ovarian phenotypes.
Total RNA was extracted from 5 replicates of 10 males or 10 full (NA)/partial (Pi) ovaries, as described in . Total RNA was purified from potential DNA contamination by DNase treatment (Turbo DNAse, Ambion, Applied Biosystems, Austin, TX). First-strand cDNA synthesis was performed from 500 ng of total RNA using the Superscript III enzyme (Invitrogen, Cergy-Pontoise, France) and oligodT primers, according to the Manufacturer’s instructions. For each biological sample, 4 ng of cDNA was spotted in duplicate in a 96-well plate (Microlab star, Hamilton, Bonaduz, Switzerland). Quantitative PCR was performed using LightCycler LC480 system (Roche, Meylan, France) as follows: 5 min at 95°C, 35 times [15 s at 95°C, 10s at 58°C, 20 s at 72°C], 20 s at 70°C. A melting curve was recorded at the end of the PCR amplification to confirm that a unique transcript product had been amplified. The reaction mixture consisted of 0.5 µM of each primer, 5 µL of Fast SYBR-Green Master Mix (Roche, Meylan, France), and 4 µL of diluted cDNA (corresponding to 4 ng of cDNA). Primers used for quantitative PCR are summarized in Additional File 1. In order to calculate PCR efficiencies, standard curves were plotted using seven dilutions (10–107 copies) of a previously amplified PCR product purified using Nucleospin Extract II kit (Macherey-Nagel, Hoerdt, France). Expression data were estimated by calculating E−Cp, where E corresponds to the efficiency of the PCR reaction, and Cp to the crossing point . Candidate gene expression was normalized by the geometric mean of the expression level of three housekeeping genes (Ribosomal L6, β-tubulin, and Elongation factor 1γ), and analyzed by Wilcoxon’s test. The p-values were then adjusted using false discovery rate’s correction (FDR, R software, version 2.12).
More than 12,000 unigenes sequenced in cDNA libraries
However, most unigenes were obtained from the normalized library and the ovary libraries (Fig. 2C). In addition, the overlap between libraries was low, suggesting that the sampling effort should be increased to perform a transcriptomic analysis at the gene level. Indeed, 60% of the unigenes were defined by a single EST (Fig. 2B). Furthermore, the two aposymbiotic libraries (OA1 and OA2) only partially overlapped (Fig. 2C), sharing 345 unigenes, corresponding to 16% of OA1 and 26% of OA2, respectively.
Functional annotation was performed on the 12,511 unigenes using Blast against various databases and using the Gene Ontology procedure (method summarized in Fig. 1B, results in Fig. 2D). Gene Ontology (GO) is a structured vocabulary describing gene products with terms taken from different ontologies, such as the molecular function or biological process resulting from the coordinate action of molecular functions. We chose high e-value cut-offs because of the ancient divergence between A. tabida and the closest sequenced genomes. In addition, divergence can be very high for fast-evolving genes like immune effectors. The principal database sources for the GO annotation were UniprotKB (55%), Flybase (21%) and Mouse Genome Informatics (19%). Around 70% of the unigenes had Blast similarities, mainly against N. vitripennis (15 %), Apis mellifera (13%), Harpegnathos saltator (11%), Camponotus floridanus (11%), Solenopsis invicta (8%) and Tribolium castaneum (2%), with an e-value lower than e-20 for more than 55% of the unigenes. Undetectable similarity could correspond to the UTR part of the cDNA, or to species-specific genes. Around 40% of unigenes were annotated after the Blast2go annotation procedure for High Scoring Pair (HSP) over a hit length coverage cut-off of 0%. We used permissive annotation parameters since our goal was to keep the maximum functional annotation even if it involves only a very short portion of the unigene (e.g. a domain). Adding Interproscan prediction and running the Annex augmentation procedure increased the number of unigenes annotated. While we kept the unigenes/GO datatset corresponding to the minimum HSP coverage percentage, the mean number of GO terms assigned per unigene was 1.66 GO (Fig. 2E).
Functional analysis of the symbiotic interaction
Functions under-represented in wasp ovaries in response to Wolbachia infection
OA1, level 3
(n = 95)
interspecies interaction between organisms
stem cell maintenance
mRNA splice site selection
reproductive developmental process
generation of precursor metabolites and energy
cellular component organization and biogenesis
ensheathment of neurons
OA2, level 3
(n = 16)
stem cell maintenance
OA1, level 4
glutathione peroxidase activity
cell adhesion molecule binding
oxidoreductase activity, acting on superoxide radicals as acceptor
transferase activity, transferring alkyl or aryl (other than methyl) groups
flavine monoNucleotid binding
protein transmembrane transporter activity
oxidoreductase activity, acting on peroxide as acceptor
transferase activity, transferring sulfur-containing groups
lipoic acid binding
substrate-specific transmembrane transporter activity
translation elongation factor activity
heme-copper terminal oxidase activity
oxidoreductase activity, acting on heme group of donors
oxidoreductase activity, acting on NADH or NADPH
hydrolase activity, acting on acid anhydrides
active transmembrane transporter activity
electron carrier activity
fatty acid binding
OA2, level 4
(n = 92)
oxidoreductase activity, acting on paired donors, with incorporation or reduction of molecular oxygen
transferase activity, transferring sulfur-containing groups
thiamin pyrophosphate binding
fatty acid binding
Functional enrichment analysis
SSH2A vs. OS
cellular component organization & biogenesis
ribonucleoprotein complex biogenesis & assembly
organelle organization & biogenesis
ribosome biogenesis & assembly
structural constituent of ribosome
SSH1A vs. OS
ferric iron binding
SSH2S vs. OS
no significant terms
SSH1S vs. OS
no significant terms
Suppressive Subtraction Hybridizations were performed between wasps challenged with S. typhimurium and unchallenged wasps (SSHs C-NC) in order to detect immune genes. However, the SSH-C was saturated with the antimicrobial peptide Hymenoptaecin, and so was not informative.
Expression of genes related to immunity (broad sense), programmed cell death, and oogenesis
Overall analysis of differential gene expression in response to Wolbachia removal
If we focus on genes involved in immunity (Toll, Imd, JNK, JAK-STAT, RNAi pathways), expression patterns were relatively clear in males. Indeed, NA males showed a very limited response to Wolbachia removal, whereas Pi3 males mainly showed an under-expression of these immune genes in symbiotic individuals. On the opposite, immune genes were mainly over-expressed in symbiotic ovaries of both strains, with however a higher differential expression in Pi3 ovaries. This difference could be attributable to the ovarian phenotype, but also to other phenotypic traits controlled by the female genotype. Furthermore, numerous genes involved in immune functions (e.g. Toll, Cactus, Dorsal, Basket) may also play an important role during the development. Since their transcripts may accumulate during oogenesis, expression results associated with these genes have to be interpreted with caution in aposymbiotic females whose oogenetic process is markedly affected. Curiously, in most of these immune pathways, but particularly the Toll and JAK-STAT pathways, expression profiles depended on the gene being investigated. Indeed, genes upstream in the pathways were mainly over-expressed in symbiotic individuals, whereas downstream effectors, such as anti-microbial peptides and TEPs, were mainly down-regulated in response to symbiosis. It is also interesting to note that gene expression was generally much lower in ovaries than in males, suggesting that this tissue may display limited immuno-competency.
In order to study immunity in its broad sense, we also took into account processes involved in the stress response and programmed cell death, as they can also be involved in limiting bacterial infection. Unfortunately, very few genes involved in canonical pathways of apoptosis and autophagy were detected among the libraries, which limited the scope of our investigation. Expression patterns were once again very different in NA males and Pi3 males. In Pi3 males, genes involved in stress and programmed cell death were mainly under-expressed in response to symbiosis. It is difficult to interpret the response of NA males to symbiosis, since the very few genes that were differentially regulated were either up or down-regulated within a given pathway. In the ovaries, where cytological analyses have highlighted apoptotic and autophagic processes in aposymbiotic ovaries [,Rancès, pers. com.], processes associated with PCD were either unchanged in response to symbiosis (NA strain) or, surprisingly, over-expressed in symbiotic ovaries (Pi3 strain). In Pi3 and NA ovaries, genes involved in the stress response (detoxification, folding) were mainly under-expressed in response to symbiosis, which confirms the trend highlighted by the analyses of EST libraries.
Wolbachia is known to play a role in oogenesis completion in A. tabida , and to restore fertility to the Sxlf4 D. melanogaster mutant . Therefore, we studied the expression of genes known to be involved in sex determination in Drosophila (Sxl, Ix) and also in oogenesis and embryogenesis. Expression of Sxl and Ix was not limited to one sex, as shown by , and varied in response to symbiosis in all the populations investigated. However, the function of Sxl in sex determination seems to be restricted to Drosophila and could have another role (not female specific) in other insect species . Genes involved in oogenesis and embryogenesis were all over-expressed in symbiotic ovaries, and more significantly so in the Pi ovaries. These findings are thus congruent with the ovarian phenotype of aposymbiotic females (without eggs in the Pi3 strain, and with a few eggs in the NA strain). Patterns in gene expression could be explained by the ovarian phenotype’s being related either to a direct role in oogenesis or to mRNA storage in the eggs for subsequent embryo development.
Phenotypic effects of Wolbachia on host biology are being increasingly reported in arthropod species . Furthermore, growing numbers of Wolbachia genomes have now been sequenced from strains inducing various phenotypic effects [45–49], which provides essential information about the biology and evolution of the symbiont. However, very few studies have focused on the overall response of the host to the presence of Wolbachia in natural associations [20, 21, 23, 24]. Most studies have focused on host response after stable [20, 21] or transient infection by Wolbachia , or in cell cultures [23, 51].
The first goal of this work was to generate a first reference transcriptome of A. tabida, a model system both for host/Wolbachia  and host/parasitoid interactions [52, 53]. The 12,511 unigenes we isolated from the wasp A. tabida constitute a valuable resource for further genetic studies of these interactions. For example, the host transcriptional response to parasitoid attack has been studied in D. melanogaster using microarrays , but large-scale analyses in wasps are currently lacking. The genetic information provided here may help to fill this gap.
The second objective was to detect differentially-represented functions in response to symbiosis. Direct analysis of the libraries was limited by the sequencing depth at the gene level, and thus required an analysis based on the GO term level. Several genes associated with candidate functions were extracted from the current ESTs dataset, and were thoroughly studied through qRT-PCR. The current transcriptomic map can now be used as a backbone for high-throughput sequencing (e.g. Illumina) to provide an accurate global analysis of genes that are differentially expressed in response to symbiosis.
Through different approaches, we identified various biological processes that were transcriptionally affected by Wolbachia removal. Indeed, almost all the genes we studied using qRT-PCR were differently regulated in male and/or females at least in one population. The difference in gene expression was generally less than 2-fold, and could not have been detected by microarray analyses. The influence of Wolbachia removal on gene expression was expected in the ovaries, where the absence of Wolbachia dramatically alters the ovarian structure. In males however, the current absence of extended phenotype in response to Wolbachia removal would not have suggested changes in gene expression.
As has already been pointed out, these results must be treated with caution. In aposymbiotic individuals, antibiotic treatment could indeed have directly influenced mitochondrial metabolism  and gene expression because of its general cytotoxic effect. Antibiotics could also have indirectly influenced gene expression through the elimination of other bacteria (e.g. present in the gut community ). We are confident that the variations observed must have been due (or at least largely due) to Wolbachia infection. Indeed, we would expect the direct effects of antibiotics to affect both strains similarly. However, we found that (1) direct effects of the antibiotic treatment may be very limited, as very few genes were differentially regulated in NA males, (2) no gene (except Transferrin) was differentially expressed in all comparisons, and (3) as expected, the Pi3 strain was more sensitive to Wolbachia removal than the NA strain. These results suggest either that changes in gene expression are due to the host genotype in response to Wolbachia removal, or that the potential antibiotic effect impacts the expression of genes also involved in the ovarian phenotype.
As variation in dependence phenotype is determined by the host nuclear genotype , we studied transcriptional response to symbiosis in two populations with extreme ovarian phenotypes. However, the comparison between Pi3 and NA populations could have been obscured by their different evolutionary histories and symbiotic status regarding Wolbachia strains and other bacteria. To discard this hypothesis, we subsequently measured the expression of some genes in two strains originating from a same population (Saintte Foy-lès-Lyon, France), but exhibiting different ovarian phenotypes . These strains were genetically related and both triply-infected, and similar patterns were observed as in the comparison between Pi3 and NA ovaries . Hence, variation in gene expression in response to symbiosis must be driven by the genetic background associated with the dependence phenotype.
Growing evidence shows that the presence of a symbiont can dramatically affect host immunity . For instance, Wigglesworthia reduces susceptibility of the tsetse fly to infection by Trypanosoma by modulating PGRP-LB [58, 59], and the male-killer Spiroplasma weakens antimicrobial expression in D. melanogaster . Immuno-modulation by a symbiont could thus be a way of circumventing the host’s immune system and/or to increase host fitness and ability to cope with common pathogens, thus ensuring that the symbiont is maintained within the host. Although Wolbachia is hidden in a host-derived vacuole, the transcriptomic analyses presented here suggest that the host organism detects its presence, and that Wolbachia may not only adopt an ‘immune-escape’ strategy. Indeed, Wolbachia seems to influence host immune system in its broad sense, including both canonical pathways and the stress response to external stimuli. These observations led us to wonder how Wolbachia is detected within the cell, how Wolbachia evades the host immune system, and what are the consequences of these manipulations on host cell physiology.
In the present study, most of the canonical immune PGRP receptors were differentially-regulated in the presence of Wolbachia, probably through lipoprotein or polysaccharide binding, and the outcome of the interaction tended towards under-expression of immune effectors of the Toll, Imd and JAK-STAT pathways. Even when the regulation cascade was too complex to analyze, the expression patterns of most immune genes were modified in response to symbiosis, suggesting that Wolbachia may adopt an active strategy of immune evasion in A. tabida. However, as few immune genes from the Toll signaling pathway are also known to play a role in development, expression data have to be interpreted with caution with respect to the important development defect of ovaries in aposymbiotic females. The regulation appeared to be tissue or sex-specific, immune genes being expressed to a greater extent in males than in ovarian tissues. Wolbachia is mainly concentrated in the ovaries of females, whereas they are spread more widely throughout the male body . Hence, modulation of immune pathways could be both gene- and tissue-specific, as shown in the differential immune regulation of bacteriocytes vs. whole body in Sitophilus zeamais . The immune response to Wolbachia also seems to be host strain-specific, with the Pi3 strain generally exhibiting a more pronounced pattern than the NA strain. Finally, the immune response to Wolbachia seems to be host-specific, as Drosophila simulans did not repress or induce antimicrobial peptides production , whereas the D. melanogaster cell line over-expressed antimicrobial peptides in response to Wolbachia infection . Similarly, the presence of Wolbachia tends to increase immune gene expression in the mosquito hosts when stably introduced [20, 21, 50].
By comparing aposymbiotic and symbiotic tissues of A. tabida, we also highlighted the influence of Wolbachia on host immunity in its broad sense, and especially on the regulation of cell homeostasis and the oxidative environment, which are known to play a key role in physiological responses to invasion by pathogens. Indeed, processes involved in the control of the oxidative environment were highlighted both in in silico and in vitro subtractions, and confirmed by qRT-PCR. Given these observations, we further demonstrated the influence of Wolbachia on iron homeostasis and oxidative stress regulation in A. tabida [8, 14]. We confirmed the differential expression of Ferritin, a protein involved in iron storage and transport, in males, females and ovaries from the Pi strain . Since the control of iron homeostasis pleiotropically affects apoptotic and oogenetic processes, its perturbation may have played a role in the evolution of dependence . Interestingly, ROS also interfere with oogenesis in mosquitoes  and Drosophila , probably by controlling apoptotic checkpoints . The influence of Wolbachia on iron homeostasis was not restricted to A. tabida, since we demonstrated a similar effect in D. simulans and in an A. aegypti cell line . Hence, processes highlighted in an association in which Wolbachia induces an extreme phenotype also shed light on more general processes in host/Wolbachia interactions. In the present study, the stress response was not restricted to iron regulation, as other chaperones and enzymes involved in detoxification were also differentially expressed in response to Wolbachia symbiosis, in both males and females. These results suggest a general regulation of the oxidative environment, not solely restricted to the ovaries where the phenotype is observed. Genes involved in the stress response were generally over-expressed in aposymbiotic individuals, suggesting either that Wolbachia has a protective effect on host physiology/immunity or that host compensatory mechanisms have been developed to reduce the harmful impact of the presence of Wolbachia . Interestingly, we observed a differential response in Pi3 vs. NA strains through quantitative RT-PCR, which was confirmed in another population with similar phenotypes . These results suggest that host gene expression has evolved to tolerate the presence of Wolbachia, and that the Pi3 genotype is more sensitive to its presence.
Finally, some striking similarities emerge when these results are compared with two other models that have been used in similar studies, but which have radically different extended phenotypes and types of relationships (i.e. Armadillidium vulgare/Wolbachia and Sitophilus orizae/SOPE) [66, 67]. Functions such as oxidative stress regulation [8, 14] and classical immune pathways  have already been highlighted, and appear again as being shared between symbiotic associations. Apoptosis has previously been highlighted in A. tabida, owing to the strong cellular phenotype induced by the removal of Wolbachia , but also appears to be shared by the other associations. Finally, new functions, such as autophagy, have been detected in all three associations, raising the possibility that this pathway also plays a central role in symbiotic interactions. All these functions are also shared in host-pathogen interactions, suggesting the existence of a common language between bacteria and their hosts, whatever the form their interaction takes. However, a detailed analysis of these pathways revealed that they may be under- or over-regulated, depending on the symbiotic association. These differences in gene regulation may reflect different co-evolutionary dynamics (e.g. an arms race or cooperation between the partners), and/or different selective pressures due to symbiont location. When symbionts are not restricted to specific tissues, deleterious side effects on other traits, e.g. the response to pathogens or developmental processes modulated by the pleiotropic action of genes, may indeed limit or shape the expression of these pathways.
In this study, we identified 12,511 unigenes from the parasitoid wasp A. tabida, which can now facilitate future genetic studies on host/Wolbachia and host/parasitoid interactions. We also highlighted that Wolbachia might interfere with the expression of genes involved in development, PCD and immunity, especially through the regulation of oxidative stress. These results confirm that Wolbachia does not only impact its host reproduction, but may also influence more globally the biology and physiology of its hosts with potential unprecedented effects on the evolution of their life history.
We would like to thank two anonymous reviewers for their helpful comments on the manuscript, and Suzanne Peyer for reviewing the English text. We would like to express our sincere thanks to Christine Oger (DTAMB, IFR 41, Université de Lyon) for her help in using the Microlabstar Hamilton. A. tabida sequences were obtained within the framework of the “Functional Genomics and Immune Signaling in Invertebrate Endosymbiosis” program, conducted in collaboration with the Centre National de Séquençage, Genoscope (Evry, France). This work was supported by funding from UMR CNRS 5558, IFR 41 and GDR 2153, a grant from the Agence Nationale de la Recherche (ANR-06-BLANC-0316 "EndoSymbArt"), and a grant from the Fondation Innovations en Infectiologie (FINOVI 005).
This article has been published as part of BMC Microbiology Volume 11 Supplement 1, 2012: Arthropod symbioses: from fundamental studies to pest and disease mangement. The full contents of the supplement are available online at http://0-www.biomedcentral.com.brum.beds.ac.uk/1471-2180/12?issue=S1.
- McFall-Ngai MJ: Unseen forces: the influence of bacteria on animal development. Dev Biol. 2002, 242: 1-14. 10.1006/dbio.2001.0522.View ArticlePubMedGoogle Scholar
- Ivanov II, Littman DR: Modulation of immune homeostasis by commensal bacteria. Curr Opin Microbiol. 2011, 14: 106-114. 10.1016/j.mib.2010.12.003.PubMed CentralView ArticlePubMedGoogle Scholar
- Ryu J-H, Kim S-H, Lee H-Y, Bai JY, Nam YD, Bae JW, Lee DG, Shin SC, Ha E-M, Lee W-J: Innate immune homeostasis by the homeobox gene caudal and commensal-gut mutualism in Drosophila. Science. 2008, 319: 777-782. 10.1126/science.1149357.View ArticlePubMedGoogle Scholar
- Troll JV, Adin DM, Wier AM, Paquette N, Silverman N, Goldman WE, Stadermann FJ, Stabb EV, McFall-Ngai MJ: Peptidoglycan induces loss of a nuclear peptidoglycan recognition protein during host tissue development in a beneficial animal-bacterial symbiosis. Cell Microbiol. 2009, 11: 1114-1127. 10.1111/j.1462-5822.2009.01315.x.PubMed CentralView ArticlePubMedGoogle Scholar
- Werren JH, Baldo L, Clark ME: Wolbachia: master manipulators of invertebrate biology. Nat Rev Microbiol. 2008, 6: 741-751. 10.1038/nrmicro1969.View ArticlePubMedGoogle Scholar
- Dedeine F, Vavre F, Fleury F, Loppin B, Hochberg ME, Bouletreau M: Removing symbiotic Wolbachia bacteria specifically inhibits oogenesis in a parasitic wasp. Proc Natl Acad Sci U S A. 2001, 98: 6247-6252. 10.1073/pnas.101304298.PubMed CentralView ArticlePubMedGoogle Scholar
- Dedeine F, Boulétreau M, Vavre F: Wolbachia requirement for oogenesis: occurrence within the genus Asobara (Hymenoptera, Braconidae) and evidence for intraspecific variation in A. tabida. Heredity. 2005, 95: 394-400. 10.1038/sj.hdy.6800739.View ArticlePubMedGoogle Scholar
- Kremer N, Dedeine F, Charif D, Finet C, Allemand R, Vavre F: Do variable compensatory mechanisms explain the polymorphism of the dependence phenotype in the Asobara tabida-Wolbachia association?. Evolution. 2010, 64: 2969-2979.PubMedGoogle Scholar
- Pannebakker BA, Loppin B, Elemans CPH, Humblot L, Vavre F: Parasitic inhibition of cell death facilitates symbiosis. Proc Natl Acad Sci U S A. 2007, 104: 213-215. 10.1073/pnas.0607845104.PubMed CentralView ArticlePubMedGoogle Scholar
- McCall K: Eggs over easy: cell death in the Drosophila ovary. Dev Biol. 2004, 274: 3-14. 10.1016/j.ydbio.2004.07.017.View ArticlePubMedGoogle Scholar
- Böhme L, Rudel T: Host cell death machinery as a target for bacterial pathogens. Microbes Infect. 2009, 11: 1063-1070. 10.1016/j.micinf.2009.08.014.View ArticlePubMedGoogle Scholar
- Vavre F, Kremer N, Pannebakker BA, Loppin B, Mavingui P: Is symbiosis evolution influenced by the pleiotropic role of programmed cell death in immunity and development?. Insect Symbiosis. Edited by: K. Bourtzis and T. A. Miller. 2008, CRC Press, Boca Raton, 3: 57-75.Google Scholar
- Brownlie JC, Cass BN, Riegler M, Witsenburg JJ, Iturbe-Ormaetxe I, McGraw EA, O'Neill SL: Evidence for metabolic provisioning by a common invertebrate endosymbiont, Wolbachia pipientis, during periods of nutritional stress. PLoS Pathog. 2009, 5: e1000368-10.1371/journal.ppat.1000368.PubMed CentralView ArticlePubMedGoogle Scholar
- Kremer N, Voronin D, Charif D, Mavingui P, Mollereau B, Vavre F: Wolbachia interferes with ferritin expression and iron metabolism in insects. PLoS Pathog. 2009, 5: e1000630-10.1371/journal.ppat.1000630.PubMed CentralView ArticlePubMedGoogle Scholar
- Brennan LJ, Keddie BA, Braig HR, Harris HL: The endosymbiont Wolbachia pipientis induces the expression of host antioxidant proteins in an Aedes albopictus cell line. PLoS One. 2008, 3: e2083-10.1371/journal.pone.0002083.PubMed CentralView ArticlePubMedGoogle Scholar
- Molina-Cruz A, DeJong RJ, Charles B, Gupta L, Kumar S, Jaramillo-Gutierrez G, Barillas-Mury C: Reactive oxygen species modulate Anopheles gambiae immunity against bacteria and Plasmodium. J Biol Chem. 2008, 283: 3217-3223.View ArticlePubMedGoogle Scholar
- Ryu J-H, Ha E-M, Lee W-J: Innate immunity and gut-microbe mutualism in Drosophila. Dev Comp Immunol. 2010, 34: 369-376. 10.1016/j.dci.2009.11.010.View ArticlePubMedGoogle Scholar
- Teixeira L, Ferreira A, Ashburner M: The bacterial symbiont Wolbachia induces resistance to RNA viral infections in Drosophila melanogaster. PLoS Biol. 2008, 6: e2-View ArticlePubMedGoogle Scholar
- Osborne SE, Leong YS, O'Neill SL, Johnson KN: Variation in antiviral protection mediated by different Wolbachia strains in Drosophila simulans. PLoS Pathog. 2009, 5: e1000656-10.1371/journal.ppat.1000656.PubMed CentralView ArticlePubMedGoogle Scholar
- Moreira LA, Iturbe-Ormaetxe I, Jeffery JA, Lu G, Pyke AT, Hedges LM, Rocha BC, Hall-Mendelin S, Day A, Riegler M, Hugo LE, Johnson KN, Kay BH, McGraw EA, van den Hurk AF, Ryan PA, O'Neill SL: A Wolbachia symbiont in Aedes aegypti limits infection with dengue, Chikungunya, and Plasmodium. Cell. 2009, 139: 1268-1278. 10.1016/j.cell.2009.11.042.View ArticlePubMedGoogle Scholar
- Bian G, Xu Y, Lu P, Xie Y, Xi Z: The endosymbiotic bacterium Wolbachia induces resistance to dengue virus in Aedes aegypti. PLoS Pathog. 2010, 6: e1000833-10.1371/journal.ppat.1000833.PubMed CentralView ArticlePubMedGoogle Scholar
- Iturbe-Ormaetxe I, O'Neill SL: Wolbachia-host interactions: connecting phenotype to genotype. Curr Opin Microbiol. 2007, 10: 221-224. 10.1016/j.mib.2007.05.002.View ArticlePubMedGoogle Scholar
- Xi Z, Gavotte L, Xie Y, Dobson SL: Genome-wide analysis of the interaction between the endosymbiotic bacterium Wolbachia and its Drosophila host. BMC Genomics. 2008, 9: 1-10.1186/1471-2164-9-1.PubMed CentralView ArticlePubMedGoogle Scholar
- Ikeya T, Broughton S, Alic N, Grandison R, Partridge L: The endosymbiont Wolbachia increases insulin/IGF-like signalling in Drosophila. Proc Biol Sci. 2009, 276: 3799-37807. 10.1098/rspb.2009.0778.PubMed CentralView ArticlePubMedGoogle Scholar
- Dedeine F, Vavre F, Shoemaker DD, Boulétreau M: Intra-individual coexistence of a Wolbachia strain required for host oogenesis with two strains inducing cytoplasmic incompatibility in the wasp Asobara tabida. Evolution. 2004, 58: 2167-2174.View ArticlePubMedGoogle Scholar
- Zouache K, Voronin D, Tran-Van V, Mavingui P: Composition of bacterial communities associated with natural and laboratory populations of Asobara tabida infected with Wolbachia. Appl Environ Microbiol. 2009, 75: 3755-3764. 10.1128/AEM.02964-08.PubMed CentralView ArticlePubMedGoogle Scholar
- David J: A new medium for rearing Drosophila in axenic conditions. Drosophila Info Serv. 1962, 36: 128-Google Scholar
- Zhu YY, Machleder EM, Chenchik A, Li R, Siebert PD: Reverse transcriptase template switching: a SMART approach for full-length cDNA library construction. Biotechniques. 2001, 30: 892-897.PubMedGoogle Scholar
- Zhulidov PA, Bogdanova EA, Shcheglov AS, Vagner LL, Khaspekov GL, Kozhemyako VB, Matz MV, Meleshkevitch E, Moroz LL, Lukyanov SA, Shagin DA: Simple cDNA normalization using kamchatka crab duplex-specific nuclease. Nucleic Acids Res. 2004, 32: e37-10.1093/nar/gnh031.PubMed CentralView ArticlePubMedGoogle Scholar
- Shagin DA, Rebrikov DV, Kozhemyako VB, Altshuler IM, Shcheglov AS, Zhulidov PA, Bogdanova EA, Staroverov DB, Rasskazov VA, Lukyanov S: A novel method for SNP detection using a new duplex-specific nuclease from crab hepatopancreas. Genome Res. 2002, 12: 1935-1942. 10.1101/gr.547002.PubMed CentralView ArticlePubMedGoogle Scholar
- Chomczynski P, Sacchi N: Single-step method of RNA isolation by acid guanidinium thiocyanate-phenol-chloroform extraction. Anal Biochem. 1987, 162: 156-159.View ArticlePubMedGoogle Scholar
- Diatchenko L, Lukyanov S, Lau YF, Siebert PD: Suppression subtractive hybridization: a versatile method for identifying differentially expressed genes. Methods Enzymol. 1999, 303: 349-380.View ArticlePubMedGoogle Scholar
- Diatchenko L, Lau YF, Campbell AP, Chenchik A, Moqadam F, Huang B, Lukyanov S, Lukyanov K, Gurskaya N, Sverdlov ED, Siebert PD: Suppression subtractive hybridization: a method for generating differentially regulated or tissue-specific cDNA probes and libraries. Proc Natl Acad Sci U S A. 1996, 93: 6025-6030. 10.1073/pnas.93.12.6025.PubMed CentralView ArticlePubMedGoogle Scholar
- Rebrikov DV, Britanova OV, Gurskaya NG, Lukyanov KA, Tarabykin VS, Lukyanov SA: Mirror orientation selection (MOS): a method for eliminating false positive clones from libraries generated by suppression subtractive hybridization. Nucleic Acids Res. 2000, 28: E90-10.1093/nar/28.20.e90.PubMed CentralView ArticlePubMedGoogle Scholar
- Ewing B, Hillier L, Wendl MC, Green P: Base-calling of automated sequencer traces using phred. I. Accuracy assessment. Genome Res. 1998, 8: 175-185.View ArticlePubMedGoogle Scholar
- Ewing B, Green P: Base-calling of automated sequencer traces using phred. II. Error probabilities. Genome Res. 1998, 8: 186-94.View ArticlePubMedGoogle Scholar
- Pertea G, Huang X, Liang F, Antonescu V, Sultana R, Karamycheva S, Lee Y, White J, Cheung F, Parvizi B, Tsai J, Quackenbush J: TIGR Gene Indices clustering tools (TGICL): a software system for fast clustering of large EST datasets. Bioinformatics. 2003, 19: 651-652. 10.1093/bioinformatics/btg034.View ArticlePubMedGoogle Scholar
- Conesa A, Götz S, García-Gómez JM, Terol J, Talón M, Robles M: Blast2GO: a universal tool for annotation, visualization and analysis in functional genomics research. Bioinformatics. 2005, 21: 3674-3676. 10.1093/bioinformatics/bti610.View ArticlePubMedGoogle Scholar
- Götz S, García-Gómez JM, Terol J, Williams TD, Nagaraj SH, Nueda MJ, Robles M, Talón M, Dopazo J, Conesa A: High-throughput functional annotation and data mining with the Blast2GO suite. Nucleic Acids Res. 2008, 36: 3420-3435. 10.1093/nar/gkn176.PubMed CentralView ArticlePubMedGoogle Scholar
- Al-Shahrour F, Díaz-Uriarte R, Dopazo J: FatiGO: a web tool for finding significant associations of Gene Ontology terms with groups of genes. Bioinformatics. 2004, 20: 578-580. 10.1093/bioinformatics/btg455.View ArticlePubMedGoogle Scholar
- Pfaffl MW: A new mathematical model for relative quantification in real-time RT-PCR. Nucleic Acids Res. 2001, 29: e45-10.1093/nar/29.9.e45.PubMed CentralView ArticlePubMedGoogle Scholar
- Starr DJ, Cline TW: A host parasite interaction rescues Drosophila oogenesis defects. Nature. 2002, 418: 76-79. 10.1038/nature00843.View ArticlePubMedGoogle Scholar
- Tarone AM, Nasser YM, Nuzhdin SV: Genetic variation for expression of the sex determination pathway genes in Drosophila melanogaster. Genet Res. 2005, 86: 31-40. 10.1017/S0016672305007706.View ArticlePubMedGoogle Scholar
- Cline TW, Dorsett M, Sun S, Harrison MM, Dines J, Sefton L, Megna L: Evolution of the Drosophila feminizing switch gene Sex-lethal. Genetics. 2010, 186: 1321-1336. 10.1534/genetics.110.121202.PubMed CentralView ArticlePubMedGoogle Scholar
- Wu M, Sun LV, Vamathevan J, Riegler M, Deboy R, Brownlie JC, McGraw EA, Martin W, Esser C, Ahmadinejad N, Wiegand C, Madupu R, Beanan MJ, Brinkac LM, Daugherty SC, Durkin AS, Kolonay JF, Nelson WC, Mohamoud Y, Lee P, Berry K, Young MB, Utterback T, Weidman J, Nierman WC, Paulsen IT, Nelson KE, Tettelin H, O'Neill SL, Eisen JA: Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: a streamlined genome overrun by mobile genetic elements. PLoS Biol. 2004, 2: E69-10.1371/journal.pbio.0020069.PubMed CentralView ArticlePubMedGoogle Scholar
- Foster J, Ganatra M, Kamal I, Ware J, Makarova K, Ivanova N, Bhattacharyya A, Kapatral V, Kumar S, Posfai J, Vincze T, Ingram J, Moran L, Lapidus A, Omelchenko M, Kyrpides N, Ghedin E, Wang S, Goltsman E, Joukov V, Ostrovskaya O, Tsukerman K, Mazur M, Comb D, Koonin E, Slatko B: The Wolbachia genome of Brugia malayi: endosymbiont evolution within a human pathogenic nematode. PLoS Biol. 2005, 3: e121-10.1371/journal.pbio.0030121.PubMed CentralView ArticlePubMedGoogle Scholar
- Klasson L, Walker T, Sebaihia M, Sanders MJ, Quail MA, Lord A, Sanders S, Earl J, O'Neill SL, Thomson N, Sinkins SP, Parkhill J: Genome evolution of Wolbachia strain wPip from the Culex pipiens group. Mol Biol Evol. 2008, 25: 1877-1887. 10.1093/molbev/msn133.PubMed CentralView ArticlePubMedGoogle Scholar
- Klasson L, Westberg J, Sapountzis P, Näslund K, Lutnaes Y, Darby AC, Veneti Z, Chen L, Braig HR, Garrett R, Bourtzis K, Andersson SGE: The mosaic genome structure of the Wolbachia wRi strain infecting Drosophila simulans. Proc Natl Acad Sci U S A. 2009, 106: 5725-5730. 10.1073/pnas.0810753106.PubMed CentralView ArticlePubMedGoogle Scholar
- Salzberg SL, Puiu D, Sommer DD, Nene V, Lee NH: Genome sequence of the Wolbachia endosymbiont of Culex quinquefasciatus JHB. J Bacteriol. 2009, 191: 1725-10.1128/JB.01731-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Kambris Z, Blagborough AM, Pinto SB, Blagrove MSC, Godfray HCJ, Sinden RE, Sinkins SP: Wolbachia stimulates immune gene expression and inhibits plasmodium development in Anopheles gambiae. PLoS Pathog. 2010, 6: e1001143-10.1371/journal.ppat.1001143.PubMed CentralView ArticlePubMedGoogle Scholar
- Hughes GL, Ren X, Ramirez JL, Sakamoto JM, Bailey JA, Jedlicka AE, Rasgon JL: Wolbachia Infections in Anopheles gambiae Cells: Transcriptomic Characterization of a Novel Host-Symbiont Interaction. PLoS Pathog. 2011, 7: e1001296-10.1371/journal.ppat.1001296.PubMed CentralView ArticlePubMedGoogle Scholar
- Eslin P, Prévost G, Havard S, Doury G: Immune resistance of Drosophila hosts against Asobara parasitoids: cellular aspects. Adv Parasitol. 2009, 70: 189-215.View ArticlePubMedGoogle Scholar
- Fleury F, Gibert P, Ris N, Allemand R: Ecology and life history evolution of frugivorous Drosophila parasitoids. Adv Parasitol. 2009, 70: 3-44.View ArticlePubMedGoogle Scholar
- Wertheim B, Kraaijeveld AR, Schuster E, Blanc E, Hopkins M, Pletcher SD, Strand MR, Partridge L, Godfray : Genome-wide gene expression in response to parasitoid attack in Drosophila. Genome Biology. 2005, 6: R94-10.1186/gb-2005-6-11-r94.PubMed CentralView ArticlePubMedGoogle Scholar
- Ballard JWO, Melvin RG: Tetracycline treatment influences mitochondrial metabolism and mtDNA density two generations after treatment in Drosophila. Insect Mol Biol. 2007, 16: 799-802. 10.1111/j.1365-2583.2007.00760.x.View ArticlePubMedGoogle Scholar
- Lee W-J: Bacterial-modulated host immunity and stem cell activation for gut homeostasis. Genes Dev. 2009, 23: 2260-2265. 10.1101/gad.1858709.PubMed CentralView ArticlePubMedGoogle Scholar
- Gross R, Vavre F, Heddi A, Hurst GDD, Zchori-Fein E, Bourtzis K: Immunity and symbiosis. Mol Microbiol. 2009, 73: 751-759. 10.1111/j.1365-2958.2009.06820.x.View ArticlePubMedGoogle Scholar
- Pais R, Lohs C, Wu Y, Wang J, Aksoy S: The obligate mutualist Wigglesworthia glossinidia influences reproduction, digestion, and immunity processes of its host, the tsetse fly. Appl Environ Microbiol. 2008, 74: 5965-5974. 10.1128/AEM.00741-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Wang J, Wu Y, Yang G, Aksoy S: Interactions between mutualist Wigglesworthia and tsetse peptidoglycan recognition protein (PGRP-LB) influence trypanosome transmission. Proc Natl Acad Sci USA. 2009, 106: 12133-8. 10.1073/pnas.0901226106.PubMed CentralView ArticlePubMedGoogle Scholar
- Anbutsu H, Fukatsu T: Evasion, suppression and tolerance of Drosophila innate immunity by a male-killing Spiroplasma endosymbiont. Insect Mol Biol. 2010, 19: 481-488.PubMedGoogle Scholar
- Mouton L, Dedeine F, Henri H, Boulétreau M, Profizi N, Vavre F: Virulence, multiple infections and regulation of symbiotic population in the Wolbachia-Asobara tabida symbiosis. Genetics. 2004, 168: 181-189. 10.1534/genetics.104.026716.PubMed CentralView ArticlePubMedGoogle Scholar
- Anselme C, Pérez-Brocal V, Vallier A, Vincent-Monegat C, Charif D, Latorre A, Moya A, Heddi A: Identification of the weevil immune genes and their expression in the bacteriome tissue. BMC Biol. 2008, 6: 43-10.1186/1741-7007-6-43.PubMed CentralView ArticlePubMedGoogle Scholar
- Bourtzis K, Pettigrew MM, O'Neill SL: Wolbachia neither induces nor suppresses transcripts encoding antimicrobial peptides. Insect Mol Biol. 2000, 9: 635-639. 10.1046/j.1365-2583.2000.00224.x.View ArticlePubMedGoogle Scholar
- DeJong RJ, Miller LM, Molina-Cruz A, Gupta L, Kumar S, Barillas-Mury C: Reactive oxygen species detoxification by catalase is a major determinant of fecundity in the mosquito Anopheles gambiae. Proc Natl Acad Sci U S A. 2007, 104: 2121-2126. 10.1073/pnas.0608407104.PubMed CentralView ArticlePubMedGoogle Scholar
- Parkes TL, Kirby K, Phillips JP, Hilliker AJ: Transgenic analysis of the cSOD-null phenotypic syndrome in Drosophila. Genome. 1998, 41: 642-651.View ArticlePubMedGoogle Scholar
- Chevalier F, Herbinière-Gaboreau J, Charif D, Mitta G, Gavory F, Wincker P, Grève P, Braquart-Varnier C, Bouchon D: Feminizing Wolbachia: a transcriptomics approach with insights on the immune response genes in Armadillidium vulgare. BMC Microbiol. 2012, 12 (Suppl 1): S1-10.1186/1471-2180-12-S1-S1.PubMed CentralView ArticlePubMedGoogle Scholar
- Vigneron A, Charif D, Vincent-Monegat C, Vallier A, Gavory F, Wincker P, Heddi A: Host gene response to endosymbiont and pathogen in the cereal weevil Sitophilus oryzae. BMC Microbiol. 2012, 12 (Suppl 1): S14-10.1186/1471-2180-12-S1-S14.PubMed CentralView ArticlePubMedGoogle Scholar
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.