Wolbachia Prophage DNA Adenine Methyltransferase Genes in Different Drosophila-Wolbachia Associations

doi:10.1371/journal.pone.0019708

PLoS One. 2011; 6(5): e19708.

Published online 2011 May 6. doi: 10.1371/journal.pone.0019708

PMCID: PMC3089641

Copyright Saridaki et al. This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original author and source are credited.

Wolbachia Prophage DNA Adenine Methyltransferase Genes in Different Drosophila-Wolbachia Associations

Aggeliki Saridaki,¹ Panagiotis Sapountzis,¹ Harriet L. Harris,² Philip D. Batista,² Jennifer A. Biliske,² Harris Pavlikaki,^3,⁴ Stefan Oehler,^3,⁵ Charalambos Savakis,^3,^5,⁶ Henk R. Braig,⁷ and Kostas Bourtzis^1,^3,⁵^*

¹Department of Environmental and Natural Resources Management, University of Ioannina, Agrinio, Greece

²Department of Biological Sciences, University of Alberta, Edmonton, Canada

³Institute of Molecular Biology and Biotechnology, Foundation for Research and Technology-Hellas, Heraklion, Crete, Greece

⁴Technological Educational Institute of Kalamata, Kalamata, Greece

⁵Biomedical Sciences Research Center “Alexander Fleming”, Vari, Greece

⁶Medical School, University of Crete, Heraklion, Crete, Greece

⁷School of Biological Sciences, University of Bangor, Bangor, Gwynedd, United Kingdom

Michael Otto, Editor

National Institutes of Health, United States of America

* E-mail: kbourtz/at/uoi.gr

Conceived and designed the experiments: KB HRB HLH CS. Performed the experiments: AS PS PDB JAB HP SO. Analyzed the data: AS PS PDB HLH HRB KB. Contributed reagents/materials/analysis tools: KB HRB HLH CS. Wrote the paper: AS KB HRB HLH PDB CS SO.

Received February 16, 2011; Accepted April 11, 2011.

This article has been cited by other articles in PMC.

Abstract

Wolbachia is an obligatory intracellular bacterium which often manipulates the reproduction of its insect and isopod hosts. In contrast, Wolbachia is an essential symbiont in filarial nematodes. Lately, Wolbachia has been implicated in genomic imprinting of host DNA through cytosine methylation. The importance of DNA methylation in cell fate and biology calls for in depth studing of putative methylation-related genes. We present a molecular and phylogenetic analysis of a putative DNA adenine methyltransferase encoded by a prophage in the Wolbachia genome. Two slightly different copies of the gene, met1 and met2, exhibit a different distribution over various Wolbachia strains. The met2 gene is present in the majority of strains, in wAu, however, it contains a frameshift caused by a 2 bp deletion. Phylogenetic analysis of the met2 DNA sequences suggests a long association of the gene with the Wolbachia host strains. In addition, our analysis provides evidence for previously unnoticed multiple infections, the detection of which is critical for the molecular elucidation of modification and/or rescue mechanism of cytoplasmic incompatibility.

Introduction

Wolbachia pipientis is an obligate intracellular symbiont belonging to the α-proteobacteria. It is thought to be present in an estimated 66% of all insect species, including disease vectors of animals and plants [1]. It has also been found in terrestrial isopods, spiders, mites, springtails and nematodes. Some Wolbachia strains can modify host reproduction and distort sex ratio by inducing parthenogenesis, feminisation, male killing or cytoplasmic incompatibility (CI) (reviewed in [2]).

Wolbachia phage particles were first observed by Wright et al. [3]. Masui et al. [4] described prophage WO, a genetic element in Wolbachia strain wTai, containing about 26 open reading frames (ORFs) in 25 kb. The phage genome includes ORFs coding for capsid proteins, baseplate assembly proteins, integrase genes, several ankyrin-like proteins, as well as a potential methyltransferase. Further studies resulted in the isolation and characterization of the bacteriophage WO from the Wolbachia strain wCauB [5]–[7]. The genome of the purified phage is linear double-stranded DNA of about 43 kb, containing 47 ORFs (wCauB2) and 45 kb and 46 ORFs (wCauB3) [7]. Sequence analysis indicated that this phage genome includes ORFs coding for a DNA packaging protein, capsid proteins, baseplate assembly proteins, tail structural proteins and for several putative toxin-like secretory proteins [6], [7].

The Drosophila melanogaster wMel genome was released in 2004 [8]. Since then, three other Wolbachia genomes have become available [9]–[11] and a number of other sequencing projects are currently in progress. Two divergent prophage WO families, WO-A and WO-B have been identified. Family WO-B can be divided into three clades [4], [12]. Distribution surveys indicate that WO-B homologs occur in at least 89% of the two main lineages of Wolbachia that infect arthropods [12], [13]. Five WO-B-like prophage regions are present in the wPip genome, with some genes identical or highly similar between prophage copies, while other genes are unique. It seems likely that extensive recombination, duplication and insertion events have occurred between copies [10], [14]. In the highly recombining wRi genome, 4 prophage segments have been detected [11], while no prophage elements could be identified in the mutualistic wBm strain of the filarial nematode Brugia malayi [9]. Insertion sequences (IS) are frequently found in WO genomes and are considered to be a major factor driving phage recombination [8], [10], [11].

Since prophage regions have been found only in Wolbachia strains having a parasitic relationship with their hosts, it has been hypothesized that they contribute to the CI phenotype [4]. However, a study in different Culex pipiens strains detected no correlation between prophage orf7 gene type and CI [15]. Another survey showed that phage associated protein Gp15 is similar to a bacterial virulence factor. This gene was partially correlated with CI expression, suggesting that it could be linked to a CI gene [16]. However, sequence analyses found no phylogenetic clustering of phage genotypes congruent with the four major Wolbachia-induced sexual alterations [13]. In Nasonia vitripennis, 12% of Wolbachia cells were found to show lytic phage development. The density of the bacteriophage correlated inversely with the density of Wolbachia bacteria [17]. While density is one of the most critical determinants of penetrance of Wolbachia-induced phenotypes [18]–[20], only virion-free Wolbachia were observed to contact the host spermatids [17]. These observations led to the Phage Density Model hypothesis, which suggests that lytic phages negatively control Wolbachia densities and expression of symbiont functions (reviewed in [21]).

Both WO-A and WO-B prophages in wMel carry a gene that encodes a putative DNA adenine methyltransferase. Methyltransferase genes are carried by many bacteriophages [22]–[24], and modified bases are common in phage genomes [25], [26]. It is thus quite possible that a WO prophage methyltransferase has a phage-specific function or acts on the bacterial genome. Virulence and lysogeny of many pathogenic bacteria such as Escherichia, Salmonella, Yersinia, Vibrio, Haemophilus, Pasteurella, Aeromonas, Actinobacillus, Klebsiella, Brucella and Rickettsia are subject to control by adenine methylation [23], [27]–[29]. In some strains, methylation is essential for viability of the pathogen [30], [31]. Adenine methylation is also involved in cell-cycle regulation of bacteria [32]. In all these cases, methylation is performed by solitary or “orphan” methyltransferases, i.e. those that are not part of a restriction-modification (R-M) system. R-M systems protect bacteria against foreign DNA or act as selfish elements [33]. Negri et al. [34] provided the first evidence that a feminizing strain of Wolbachia interferes with host genetic imprinting through cytosine methylation of leafhopper DNA.

In the context of the fact that Wolbachia modifies the genome of its host insects, the presence of a DNA methyltransferase gene in its genome becomes even more intriguing. We therefore used molecular and phylogenetic approaches to detect and characterize the Wolbachia-phage embedded gene which encodes a putative DNA methyltransferase in a number of different Wolbachia strains of known CI properties. The ability of these strains to induce (mod status) and rescue CI (resc status) has previously been characterized [18], [35]–[38]. The potential involvement of this methyltransferase-like protein in host-Wolbachia symbiosis and in the induction of reproductive alterations is discussed.

Results

Sequence and distribution of the phage methyltransferase ORFs in different Wolbachia strains

Different Wolbachia strains of known CI phenotype (Table 1) were screened for the presence of the WO methyltransferase-like genes in order to test for a possible correlation with the CI phenotype. Degenerate and specific primers were designed for each of the two wMel methyltransferase genes (Table 2) and used to amplify DNA from a number of different Wolbachia-infected Drosophila strains. Wolbachia strains wMel, wMelCS, wAu and wHa contain both methyltransferase ORFs. In contrast, strains wRi, wNo, wYak, wTei and wSan only have one of the two ORFs (met2), while wMa and wMau do not contain any of the two methyltransferase ORFs. This was confirmed by Southern blot analysis using the two wMel methyltransferase ORFs as probes (Fig. 1). The wHa met ORFs were not further characterized.

Table 1

Insect lines and Wolbachia strains used in the present study.

Table 2

PCR primers used in the present study.

Figure 1

Southern blot analysis of met genes using genomic DNA from different Wolbachia infected species.

The met1 ORFs are more closely related to each other than to the respective met2 ORFs. ORF met2_wAu contains a 2 bp GC deletion after codon 92, leading to a truncated ORF lacking the C-teminal methyltransferase domain.

Further PCR analysis using met_102F/met_269R primers indicated the presence of three PCR products in wTei: While the first was of the expected size for a met2 product, the other two were about 500 and 900 bp larger (Fig. 2A). Sequencing of the largest product revealed the presence of ISWpi1, a Wolbachia-specific insertion sequence belonging to the IS5 family [39], in met2 gene. The intermediate sized PCR product could not be cloned in several attempts, probably because of the presence of unstable repeat sequences. Surprisingly, when PCR with the same primers was performed using DNA from D. simulans STCP flies transinfected with wTei [38], only two PCR products were detected, with a striking difference in the intensity ratio (Fig. 2A).

Figure 2

PCR analysis of met2 copies in Wolbachia strain wTei.

Cloning and sequencing of the methyltransferase ORFs from the closely related wSan, wTei and wYak strains also detected two different met2 ORFs; one closely related to the met2 sequences of supergroup A Wolbachia strains (met2_wTeiA for wTei), while the other is closely related to B supergroup strains (met2_wTeiB for wTei). Specific reverse primers were designed to anneal only to the A-group-like met2 (met_1024R; Fig. 2B) or the B-group-like met2 (TeiB_1024R; Fig. 2C). PCR reactions with these primers revealed that strains wSan, wTei and wYak bear a copy of the B-group-like met2; these genes were successfully transferred by microinjections into D. simulans STCP lines by Zabalou et al. [38] (Fig. 2C). The ISWpi1-disrupted met2 ORF of wTei was also transferred together with wTei into STCP, but not the intact copy of met2_wTeiA (Fig. 2B).

Phylogenetic Analysis

All methods used to reconstruct phylogenies yielded similar results. The three methods, distance, parsimony and maximum-likelihood (ML), make different evolutionary assumptions, thus their congruence provides strong support for the deduced phylogeny. We show only the tree derived by ML estimation (Fig. 3). The phylogenetic clustering of the met2 ORFs is similar to the currently accepted clustering of the respective Wolbachia strains: all met2 gene sequences coming from supergroup A strains cluster together (wMel, wMelCS, wRi, wAu, wTei, wYak and wSan strains). The met2 gene sequences from supergroup B strains also form a cluster and are distantly related to the supergroup A gene sequences. This suggests a long association of the methyltransferase genes, and consequently of the phages harbouring them, with the respective Wolbachia strains.

Figure 3

Phylogenetic tree of Wolbachia based on met gene sequences.

RT-PCR analysis

Transcriptional analysis of met2 genes was performed by RT-PCR on cDNA samples prepared from young adult male and female flies. Met2 transcripts were detected in all samples tested (Fig. 4). The met2 copy of wTei, which is disrupted by ISWpi1, was found to be transcriptionally silent (data not shown).

Figure 4

RT-PCR analysis of met2 genes.

Discussion

The release of the first Wolbachia genome (wMel strain) revealed that it contains two DNA methyltransferase genes met1 and met2, encoded by two prophages, WO-A and WO-B respectively [8]. This finding is intriguing in the light of the fact that Wolbachia-induced CI involves modification of the insect host chromosome [40]. The presence of phage-like particles in Wolbachia-infected hosts [3], [5], [17], [41] suggests an active role of the phage in Wolbachia biology. Thus, it is tempting to speculate that, beyond controlling lysogeny of the phage, the methyltransferases might be involved in triggering reproductive alterations imposed by Wolbachia on its host. We therefore undertook a survey of Wolbachia strains of known CI status for the presence of the methyltransferases, determined their sequences and reconstructed their phylogeny.

The met1 gene is only present in a few of the tested strains; there is no correlation with CI. Wolbachia strains wMel, wRi, wNo and wPip induce CI in permissive hosts [35], [42]–[44], and they all contain at least one functional copy of the met2 gene. In contrast, strain wAu [36] does not induce or rescue CI and its only met2 ORF is disrupted. Wolbachia strains wMau and wMa do not contain the met2 gene and also do not induce CI in their hosts [43], [45], [46]. On the other hand, wYak, wTei and wSan have recently been shown to fully rescue the wRi modification, while they are unable to induce CI in their native hosts [37]. wTei is, however, able to induce 100% CI after transfer into the permissive host D. simulans [38]. In this context, it is interesting to note that our data suggest a double or multiple infection of the original host of wTei, of which not all Wolbachia strain(s) were transferred upon transinfection of D. simulans (Fig. 2). Cloning and sequencing of PCR products, using different methyltransferase primer sets, supports the presence of a hidden double infection of D. teissieri. This could explain the phenotypic shift in CI properties that is observed between the natural host and the engineered strains.

When the rescue properties of all A group strains are examined, resc⁻ strains lack a functional A-group-like met2 gene, which is always present in resc⁺ strains (Table 3); a possible correlation between met2 and CI rescue should therefore be considered. B group strains (wNo, wMa, wMau) are all resc⁺, nevertheless they do not possess an A-group-like met2 gene. The genomes of these strains have not been sequenced and the presence or absence of prophage copies has not yet been documented. While wMa and wMau do not contain any met gene, wNo has a B-group-like met2 ORF; this could reflect a different mechanism regulating CI in B group Wolbachia strains.

Table 3

Distribution of phage methyltransferases in resc+ and resc− Wolbachia strains.

Transgenic expression of wMel met2 (WD0594) in D. melanogaster was recently reported using the UAS/GAL4 system [47]. This study revealed no modification of phenotype in flies expressing met2 ubiquitously and, similarly, when expressed specifically in the ovaries, no rescue phenotype was apparent in CI crosses. Although these data suggest that constitutive expression of the met2 gene does not alone drive the CI phenotype, it is still unclear what type of regulation met2 or any of the phage-related genes are subject to and how this affects the mechanism of CI.

Southern blot analysis indicates the presence of a met-like gene also in the Wolbachia strain wUni, which is known to induce parthenogenesis in the parasitic wasp Muscidifurax uniraptor (data not shown). The distribution of the met gene in parthenogenesis-inducing Wolbachia strains remains to be investigated. Interestingly, the mutualistic Wolbachia strain, which is present in the filarial nematode Brugia malayi, neither induces reproductive alterations nor carries a copy of the DNA methyltransferase genes.

Additionally, and important for any interpretation of the role of met2, we demonstrated expression of the gene in all Wolbachia strains with RT-PCR (Fig. 4). The Wolbachia phage DNA methyltransferase may be involved in the methylation of phage, bacterial, insect host genes or a combination of them. Although Drosophila had for a long time been considered to be free of DNA methylation, both the presence of methyltransferase genes in its genome [48], [49], and of 5-methylcytosine residues in the early stages of embryonic development [50], [51] have been demonstrated. Interestingly, a Dam-like methyltransferase has been implicated in male sterility in plants [52].

Base modification in bacterial genomes is performed by two classes of DNA methyltransferases: (i) those associated with restriction-modification systems, and (ii) solitary methyltransferases that do not have a restriction enzyme counterpart. Examples of the latter are the N6-adenine methyltransferases Dam and CcrM [53], [54]. In α-Proteobacteria, CcrM methylation regulates the cell cycle in Caulobacter crescentus, Rhizobium meliloti and Agrobacterium tumefaciens and plays a role in Brucella abortus infection (reviewed in [55]). Overexpression of CcrM in these bacteria results in the accumulation of multiple chromosomes, indicative of overinitiation of DNA replication [56], [57]. Wolbachia prophage methyltransferase could regulate several aspects of the symbiont's cell cycle by imposing a specific epigenetic signal.

In silico analysis of Wolbachia prophage methyltransferase has predicted an N-terminal ParB-like nuclease domain (data not shown) similar to the ParB of the parCBA operon in E. coli, which is important for plasmid stability and resolving dimeric or multimeric plasmids. ParB nucleases have also been reported in several other plasmid genomes. ParB nucleases are Ca++ dependent endonucleases with 5′ -3′ exonuclease activity [58], [59].

The methyltransferase genes, met1 and met2, are closely related (Fig. 5). The phylogenetic clustering of the methyltransferase genes, in particular of the met2 gene, is similar to the currently accepted clustering of the arthropod Wolbachia strains, both on the level of the major division of the Wolbachia strains into two supergroups, A and B, as well as on the lower level of clades and strains. Specifically, the met2-based tree is similar to the respective wsp-based tree (data not shown). This suggests a long association of the methyltransferases, and consequently of the phages carrying them, with the harbouring Wolbachia chromosomes (Fig. 3). However, translocation of met2 gene from the phage genome to Wolbachia chromosome cannot be excluded for any of the strains studied; such an event could explain why met2 phylogeny correlates with Wolbachia phylogeny.

Figure 5

Amino acid alignment of Met1 and Met2 proteins of Wolbachia strain wMel.

Wolbachia exhibits a fascinating array of host manipulations. The elucidation of the molecular basis of the host-symbiont interaction will allow insight in the regulation of fundamental cell biological processes. Future studies will address any potential direct or indirect effect of the methyltransferase(s) in the establishment of symbiosis and/or the induction of reproductive manipulations.

Materials and Methods

Insect lines

Insect lines and Wolbachia strains used in the present study are listed in Table 1. Flies were routinely grown at 25°C on standard cornmeal medium in uncrowded vials. The Drosophila simulans STCP lines were produced by Zabalou et al. [38] who transferred Wolbachia wYak, wTei, wSan strains into the same host background using embryonic cytoplasmic injections.

PCR analysis

DNA was extracted from adult flies using the NucleoSpin Tissue kit (Macherey-Nagel) or the CTAB protocol, as previously described [60]. The DNA was used as template for PCR reactions and Southern blot analysis. About 50 to 100 ng total DNA from adult flies were used as template in PCR reactions of Wolbachia targets. The presence of Wolbachia was initially determined using Wolbachia-specific 16S rDNA primers [61]. The primers used to amplify methyltransferase gene sequences are listed in Table 2. Standard PCR analysis was performed using GoTaq® Flexi DNA polymerase (Promega). To generate DNA templates for sequencing, PCR reactions were done using the Elongase Amplification System (Invitrogen, Glasgow, UK). The PCR products were A-tailed, cloned into the pGEM-T Easy vector (Promega, Wisconsin, USA) and transformed into competent E. coli XL1-Blue MRF cells (Stratagene, Amsterdam, The Netherlands). Plasmid DNA was extracted using the Nucleospin Plasmid kit (Macherey-Nagel). Sequencing reactions were performed with the SequiThermTM Excel Long ReadTM DNA sequencing Kit-LC (Epicentre Technologies) and the DNA sequence of the inserts was determined at the laboratory of Microchemistry, FoRTH, Heraklion (Greece) on a LiCor 4200 DNA sequencer. Three to six clones were sequenced from each individual. The methyltransferase gene sequences of this study have been deposited in the EMBL database under the accession numbers AJ851152 to AJ851164 and FR796473, and in the GenBank database under the accession number JF288559.

Southern blot analysis

Genomic DNA was prepared as reported previously [62] and digested with HaeIII. Agarose gel electrophoresis of DNA and blotting to nylon membranes were carried out using standard procedures [63]. DNA probes were prepared by random hexanucleotide priming [64]. Hybridization of ³²P-labelled probes to blotted DNA was performed using standard procedures [63].

RT-PCR analysis

For each of the tested Wolbachia-Drosophila associations, total RNA from young male and female adult flies was extracted using TRIzol (Invitrogen) and treated with RNase-free DNase (Invitrogen). First-strand cDNA was synthesized from 5 µg of total RNA using reverse transcriptase (SuperScript III; Invitrogen) and random primers (Promega) and the reactions were treated with RNase H. RNA integrity was assessed using the universal Wolbachia wsp 81F/691R primers [65]. All RNA samples were tested for genomic DNA contamination by performing PCR using wsp 81F/691R primers in DNaseI-treated RNA samples which were not reverse transcribed. Transcription of met2 was detected using GoTaq® Flexi DNA polymerase (Promega).

Alignment and Model Selection

Wolbachia adenine methyltransferase nucleotide and amino acid sequences were aligned using the ClustalW Multiple Alignment algorithm implemented in Geneious v.5.3.3. [66]. The appropriate evolutionary model JTT+ Γ was selected by the Akaike Information Criterion (AIC) using ProtTest v.2.4 [67]. Models of substitution for nucleotide alignments were selected using AIC in jModeltest v.0.1.1. [68]. The appropriate evolutionary model was TPM1uf+I+Γ.

Phylogenetic analysis

The evolutionary history was inferred by maximum likelihood criterion using PAUP* v.4.0b10 for the nucleotide alignment and PHYML for the protein alignment [69], [70]. Phylogenetic trees were generated using ML bootstrap analysis in PAUP (100 pseudoreplicates of heuristic search with 10 random sequence). The maximum likelihood method conducted using PHYML was performed using 100 bootstrap replicates, a fixed proportion of invariable sites, an estimated gamma distribution parameter and optimized topology, branch lengths and rate parameters. ML trees generated are midpoint rooted using Archaeopteryx v. 0.957b [71].

Methyltransferase nomenclature

The methyltransferase gene sequences were named based on the following system. Each methyltransferase gene name is composed of “met” (in italics) denoting methyltransferase, followed by the numbers 1 or 2, indicating their origin from phage WO-A or WO-B, respectively (phage nomenclature according to [8]) and concluded by the name of the Wolbachia strain harbouring the methyltransferase gene (eg wRi, wMel, etc). Using this nomenclature, the two methyltransferase genes present in the wMel strain are named met1_wMel and met2_wMel.

Acknowledgments

The authors wish to thank Richard Stouthamer for his support and encouragement as well as for critically reviewing an earlier version of the manuscript; Richard Stouthamer, Fabrice Vavre and Bill Ballard for insect samples and George Tsiamis for his help on the phylogenetic analysis.

Footnotes

Competing Interests: The authors have declared that no competing interests exist.

Funding: This research was supported in part by the European Union grant QLK3-2000-01079 (to KB, HRB, CS, and Richard Stouthamer), by the European Community's Seventh Framework Programme CSA-SA_REGPROT-2007-1 (to KB), by BBSRC, United Kingdom (to HRB), and by the NSERC, Canada grant (to HLH). The funders had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.

References

1. 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:215–220. [PMC free article] [PubMed]

2. Saridaki A, Bourtzis K. Wolbachia: more than just a bug in insects genitals. Curr Opin Microbiol. 2010;13:67–72. [PubMed]

3. Wright JD, Sjostrand FS, Portaro JK, Barr AR. Ultrastructure of Rickettsia-like microorganism Wolbachia pipientis and associated virus-like bodies in mosquito Culex pipiens. J Ultrastr Res. 1978;63:79–85. [PubMed]

4. Masui S, Kamoda S, Sasaki T, Ishikawa H. Distribution and evolution of bacteriophage WO in Wolbachia, the endosymbiont causing sexual alterations in arthropods. J Mol Evol. 2000;51:491–497. [PubMed]

5. Masui S, Kuroiwa H, Sasaki T, Inui M, Kuroiwa T, et al. Bacteriophage WO and virus-like particles in Wolbachia, an endosymbiont of arthropods. Biochem Biophys Res Comm. 2001;283:1099–1104. [PubMed]

6. Fujii Y, Kubo T, Ishikawa H, Sasaki T. Isolation and characterization of the bacteriophage WO from Wolbachia, an arthropod endosymbiont. Biochem Biophys Res Comm. 2004;317:1183–1188. [PubMed]

7. Tanaka K, Furukawa S, Nikoh N, Sasaki T, Fukatsu T. Complete WO phage sequences reveal their dynamic evolutionary trajectories and putative functional elements required for integration into the Wolbachia genome. Appl Environ Microbiol. 2009;75:5676–5686. [PMC free article] [PubMed]

8. Wu M, Sun LV, Vamathevan J, Riegler M, Deboy R, et al. Phylogenomics of the reproductive parasite Wolbachia pipientis wMel: A streamlined genome overrun by mobile genetic elements. PLoS Biol. 2004;2:1–15. [PMC free article] [PubMed]

9. Foster J, Ganatra M, Kamal I, Ware J, Makarova K, et al. The Wolbachia genome of Brugia malayi: endosymbiont evolution within a human pathogenic nematode. PLoS Biol. 2005;3:e121. [PMC free article] [PubMed]

10. Klasson L, Walker T, Sebaihia M, Sanders MJ, Quail MA, et al. Genome evolution of Wolbachia strain wPip from the Culex pipiens group. Mol Biol Evol. 2008;25:1877–1887. [PubMed]

11. Klasson L, Westberg J, Sapountzis P, Näslund K, Lutnaes Y, et al. The mosaic genome structure of the Wolbachia wRi strain infecting Drosophila simulans. Proc Natl Acad Sci U S A. 2009;106:5725–5730. [PubMed]

12. Bordenstein SR, Wernegreen JJ. Bacteriophage flux in endosymbionts (Wolbachia): infection frequency, lateral transfer, and recombination rates. Mol Biol Evol. 2004;21:1981–1991. [PubMed]

13. Gavotte L, Henri H, Stouthamer R, Charif D, Charlat S, et al. A survey of the bacteriophage WO in the endosymbiotic bacteria Wolbachia. Mol Biol Evol. 2007;24:427–435. [PubMed]

14. Ishmael N, Dunning Hotopp JC, Ioannidis P, Biber S, Sakamoto J, et al. Extensive genomic diversity of closely related Wolbachia strains. Microbiology. 2009;155:2211–2222. [PMC free article] [PubMed]

15. Sanogo YO, Eitam A, Dobson SL. No evidence for bacteriophage WO orf7 correlation with Wolbachia-induced cytoplasmic incompatibility in the Culex pipiens complex (Culicidae: Diptera). J Med Entomol. 2005;42:789–794. [PubMed]

16. Duron O, Bernard C, Unal S, Berthomieu A, Berticat C, et al. Tracking factors modulating cytoplasmic incompatibilities in the mosquito Culex pipiens. Mol Ecol. 2006;15:3061–3071. [PubMed]

17. Bordenstein SR, Marshall ML, Fry AJ, Kim U, Wernegreen JJ. The tripartite associations between bacteriophage, Wolbachia, and arthropods. PLoS Pathog. 2006;2:e43. [PMC free article] [PubMed]

18. Bourtzis K, Nirgianaki A, Markakis G, Savakis C. Wolbachia infection and cytoplasmic incompatibility in Drosophila species. Genetics. 1996;144:1063–1073. [PubMed]

19. Poinsot D, Bourtzis K, Markakis G, Savakis C, Merçot H. Wolbachia transfer from Drosophila melanogaster to D. simulans: host effect and cytoplasmic incompatibility relationships. Genetics. 1998;150:227–237. [PubMed]

20. Clark ME, Veneti Z, Bourtzis K, Karr TL. Wolbachia distribution and cytoplasmic incompatibility during sperm development: the cyst as the basic cellular unit of CI expression. Mech Dev. 2003;120:185–198. [PubMed]

21. Kent BN, Bordenstein SR. Phage WO of Wolbachia: lambda of the endosymbiont world. Trends Microbiol. 2010;18:173–181. [PMC free article] [PubMed]

22. Blaisdell BE, Campbell AM, Karlin S. Similarities and dissimilarities of phage genomes. Proc Natl Acad Sci U S A. 1996;93:5854–5859. [PubMed]

23. Murphy KC, Ritchie JM, Waldor MK, Løbner-Olesen A, Marinus MG. Dam methyltransferase is required for stable lysogeny of the Shiga toxin (Stx2)-encoding bacteriophage 933W of enterohemorrhagic Escherichia coli O157:H7. J Bacteriol. 2008;190:438–441. [PMC free article] [PubMed]

24. Lynch KH, Stothard P, Dennis JJ. Genomic analysis and relatedness of P2-like phages of the Burkholderia cepacia complex. BMC Genomics. 2010;11:599. [PMC free article] [PubMed]

25. Sternberg N, Coulby J. Cleavage of the bacteriophage P1 packaging site (pac) is regulated by adenine methylation. Proc Natl Acad Sci U S A. 1990;87:8070–8074. [PubMed]

26. Hattman S, Sun W. Escherichia coli OxyR modulation of bacteriophage Mu mom expression in dam⁺ cells can be attributed to its ability to bind P_mom promoter DNA. Nucleic Acids Res. 1997;25:4385–4388. [PMC free article] [PubMed]

27. Zhang JZ, Hao JF, Walker DH, Yu XJ. A mutation inactivating the methyltransferase gene in avirulent Madrid E strain of Rickettsia prowazekii reverted to wild type in the virulent revertant strain. Evir Vaccine. 2006;24:2317–2323. [PubMed]

28. Pouillot F, Fayolle C, Carniel E. A putative DNA adenine methyltransferase is involved in Yersinia pseudotuberculosis pathogenicity. Microbiology. 2007;153:2426–2434. [PubMed]

29. Low DA, Casadesús J. Clocks and switches: bacterial gene regulation by DNA adenine methylation. Curr Opin Microbiol. 2008;11:106–112. [PubMed]

30. Robertson GT, Reisenauer A, Wright R, Jensen RB, Jensen A, et al. The Brucella abortus CcrM DNA methyltransferase is essential for viability, and its overexpression attenuates intracellular replication in murine macrophages. J Bacteriol. 2000;182:3482–3489. [PMC free article] [PubMed]

31. Heusipp G, Fälker S, Schmidt MA. DNA adenine methylation and bacterial pathogenesis. Int J Med Microbiol. 2007;297:1–7. [PubMed]

32. Collier J, McAdams HH, Shapiro L. A DNA methylation ratchet governs progression through a bacterial cell cycle. Proc Natl Acad Sci U S A. 2007;104:17111–17116. [PubMed]

33. Murray NE. 2001 Fred Griffith review lecture. Immigration control of DNA in bacteria: self versus non-self. Microbiol. 2002;148:3–20. [PubMed]

34. Negri I, Franchini A, Gonella E, Daffonchio D, Mazzoglio PJ, et al. Unravelling the Wolbachia evolutionary role: the reprogramming of the host genomic imprinting. Proc Biol Sci. 2009;276:2485–2491. [PMC free article] [PubMed]

35. Merçot H, Llorente BM, Jacques M, Atlan A, Montchamp-Moreau C. Variability within the Seychelles cytoplasmic incompatibility system in Drosophila simulans. Genetics. 1995;141:1015–1023. [PubMed]

36. Hoffmann AA, Clancy DJ, Duncan J. Naturally-occurring Wolbachia infection in Drosophila simulans that does not cause cytoplasmic incompatibility. Heredity. 1996;76:1–8. [PubMed]

37. Zabalou S, Charlat S, Nirgianaki A, Lachaise D, Merçot H, et al. Natural Wolbachia infections in the Drosophila yakuba species complex do not induce cytoplasmic incompatibility but fully rescue the wRi modification. Genetics. 2004;167:827–834. [PubMed]

38. Zabalou S, Apostolaki A, Pattas S, Veneti Z, Paraskevopoulos C, et al. Multiple rescue factors within a Wolbachia strain. Genetics. 2008;178:2145–2160. [PubMed]

39. Cordaux R. ISWpi1 from Wolbachia pipientis defines a novel group of insertion sequences within the IS5 family. Gene. 2008;409:20–27. [PubMed]

40. Landmann F, Orsi GA, Loppin B, Sullivan W. Wolbachia-mediated cytoplasmic incompatibility is associated with impaired histone deposition in the male pronucleus. PLoS Pathog. 2009;5:e1000343. [PMC free article] [PubMed]

41. Gavotte L, Vavre F, Henri H, Ravallec M, Stouthamer R, et al. Diversity, distribution and specificity of WO phage infection in Wolbachia of four insect species. Insect Mol Biol. 2004;13:147–153. [PubMed]

42. Hoffmann AA, Turelli M, Immons GM. Unidirectional incompatibility between populations of Drosophila simulans. Evolution. 1986;40:692–701.

43. Bourtzis K, Dobson SL, Braig HR, O'Neill SL. Rescuing Wolbachia have been overlooked. Nature. 1998;391:852–853. [PubMed]

44. Sinkins SP, Walker T, Lynd AR, Steven AR, Makepeace BL, et al. Wolbachia variability and host effects on crossing type in Culex mosquitoes. Nature. 2005;436:257–260. [PubMed]

45. Giordano RS, O'Neill SL, Robertson HM. Wolbachia infections and the expression of cytoplasmic incompatibility in Drosophila sechellia and D. mauritiana. Genetics. 1995;140:1307–1317. [PubMed]

46. Merçot H, Poinsot D. Rescuing Wolbachia have been overlooked and discovered on Mount Kilimanjaro. Nature. 1998;391:853. [PubMed]

47. Yamada R, Iturbe-Ormaetze I, Brownlie JC, O'Neill SL. Functional test of the influence of Wolbachia genes on cytoplasmic incompatibility expression in Drosophila melanogaster. Insect Mol Biol. 2011;20:75–85. [PubMed]

48. Hung MS, Karthikeyan N, Huang B, Koo HC, Kiger J, et al. Drosophila proteins related to vertebrate DNA (5-cytosine) methyltansferases. Proc Natl Acad Sci U S A. 1999;96:11940–11945. [PubMed]

49. Lyko F, Whittaker AJ, Orr-Weaver TL, Jaenisch R. The putative Drosophila methyltransferase gene dDnmt2 is contained in a transposon-like element and is expressed specifically in ovaries. Mech Dev. 2000;95:215–217. [PubMed]

50. Lyko F, Ramsahoye BH, Jaenisch R. DNA methylation in Drosophila melanogaster. Nature. 2000;408:538–539. [PubMed]

51. Kunert N, Marhold J, Stanke J, Stach D, Lyko F. A Dnmt2-like protein mediates DNA methylation in Drosophila. Development. 2003;130:5083–5090. [PubMed]

52. Unger E, Betz S, Xu R, Cigan AM. Selection and orientation of adjacent genes influences DAM-mediated male sterility in transformed maize. Transg Res. 2001;10:409–422. [PubMed]

53. Low DA, Weyand NJ, Mahan MJ. Roles of DNA adenine methylation in regulating bacterial gene expression and virulence. Infect Immun. 2001;69:7197–7204. [PMC free article] [PubMed]

54. Løbner-Olesen A, Skovgaard O, Marinus M. Dam methylation: coordinating cellular processes. Curr Op Microbiol. 2005;8:154–160. [PubMed]

55. Wion D, Casadesús J. N6-methyl-adenine: an epigenetic signal for DNA-protein interactions. Nat Rev Microbiol. 2006;4:183–192. [PMC free article] [PubMed]

56. Wright R, Stephens C, Shapiro L. The CcrM DNA methyltransferase is widespread in the alpha subdivision of proteobacteria, and its essential functions are conserved in Rhizobium meliloti and Caulobacter crescentus. J Bacteriol. 1997;179:5869–5877. [PMC free article] [PubMed]

57. Kahng LS, Shapiro L. The CcrM DNA methyltransferase of Agrobacterium tumefaciens is essential, and its activity is cell cycle regulated. J Bacteriol. 2001;183:3065–3075. [PMC free article] [PubMed]

58. Grohmann E, Stanzer T, Schwab H. The ParB protein encoded by the RP4 par region is a Ca2+-dependent nuclease linearizing circular DNA substrates. Microbiol. 1997;143:3889–3898. [PubMed]

59. Johnson EP, Minser T, Schwab H, Burgin AB, Helinski DR. Plasmid RK2 ParB protein: Purification and nuclease properties. J Bacteriol. 1999;181:6010–6018. [PMC free article] [PubMed]

60. Navajas M, Lagnel J, Gutierrez J, Boursot P. Species wide homogeneity of nuclear ribosomal ITS2 sequences in the spider mite Tetranychus urticae contrasts with extensive mitochondrial COI polymorphism. Heredity. 1998;80:742–752. [PubMed]

61. O'Neill SL, Giordano R, Colbert AM, Karr TL, Robertson HM. 16S rRNA phylogenetic analysis of the bacterial endosymbionts associated with cytoplasmic incompatibility in insects. Proc Natl Acad Sci U S A. 1992;89:2699–2702. [PubMed]

62. Bourtzis K, Nirgianaki A, Onyango P, Savakis C. A prokaryotic dnaA sequence in D. melanogaster: Wolbachia infection and cytoplasmic incompatibility among laboratory strains. Insect Mol Biol. 1994;3:131–142. [PubMed]

63. Sambrook J, Fritsch EF, Maniatis T. Molecular Cloning: A Laboratory Manual. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 1989.

64. Feinberg AP, Vogelstein B. A technique for radiolabelling DNA restriction endonuclease fragments to high specific activity. Analyt Biochem. 1983;132:6–13. [PubMed]

65. Braig HR, Zhou W, Dobson SL, O'Neill SL. Cloning and characterization of a gene encoding the major surface protein of the bacterial endosymbiont Wolbachia pipientis. J Bacteriol. 1998;180:2373–2378. [PMC free article] [PubMed]

66. Drummond AJ, Ashton B, Buxton S, Cheung M, Cooper A, et al. Geneious, v5.3.3. 2010. Biomatters Ltd., Auckland, New Zealand. Available from http://www.geneious.com.

67. Abascal F, Zardoya R, Posada D. ProtTest: selection of best-fit models of protein evolution. Bioinformatics. 2005;21:2104–2105. [PubMed]

68. Posada D. jModelTest: phylogenetic model averaging. Mol Biol Evol. 2008;25:1253–1256. [PubMed]

69. Swofford DL. PAUP*. Phylogenetic Analysis Using Parsimony (*and Other Methods). 2003. Version 4. Sinauer Associates, Sunderland, Massachusetts.

70. Guindon S, Gascuel O. A simple, fast, and accurate algorithm to estimate large phylogenies by maximum likelihood. Syst Biol. 2003;52:696–704. [PubMed]

71. Han MV, Zmasek CM. PhyloXML: XML for evolutionary biology and comparative genomics. BMC Bioinformatics. 2009;10:356. [PMC free article] [PubMed]

72. Rousset F, Solignac M. Evolution of single and double Wolbachia symbioses during speciation in the Drosophila simulans complex. Proc Natl Acad Sci U S A. 1995;92:6389–6393. [PubMed]

73. Reynolds KT, Hoffmann AA. Male age, host effects and the weak expression or non-expression of cytoplasmic incompatibility in Drosophila strains infected by maternally-transmitted Wolbachia. Genet Res. 2002;80:79–87. [PubMed]

74. Min KT, Benzer S. Wolbachia, normally a symbiont of Drosophila, can be virulent, causing degeneration and early death. Proc Natl Acad Sci U S A. 1997;94:10792–10796. [PubMed]

75. James AC, Ballard JW. Expression of cytoplasmic incompatibility in Drosophila simulans and its impact on infection frequencies and distribution of Wolbachia pipientis. Evolution. 2000;54:1661–1672. [PubMed]

Articles from PLoS ONE are provided here courtesy of
Public Library of Science