Previous Article | Next Article ![]()
Applied and Environmental Microbiology, January 2006, p. 826-835, Vol. 72, No. 1
0099-2240/06/$08.00+0 doi:10.1128/AEM.72.1.826-835.2006
Copyright © 2006, American Society for Microbiology. All Rights Reserved.
Laboratories of Genome Dynamics, Center of Anatomy and Cell Biology, Medical University of Vienna, Vienna, Austria,1 School of Integrative Biology, University of Queensland, St. Lucia, Queensland 4072, Australia2
Received 18 August 2005/ Accepted 25 October 2005
|
|
|---|
|
|
|---|
In contrast to the well-studied Wolbachia associations in D. melanogaster and D. simulans, little is known about the occurrence of Wolbachia among American Neotropical Drosophila strains comprising two groups of species, the saltans group and the willistoni group (Fig. 1). There are presently two conflicting reports about the occurrence of Wolbachia in Neotropical Drosophila: Bourtzis et al. (4) screened a broad range of Drosophila species derived from various labs and from the Drosophila Species Stock Center (DSSC) in Bowling Green, Ohio (now held at the University of Arizona, Tucson). In their survey only two species out of the 41 stocks comprising 30 species were infected with Wolbachia. Interestingly, none of the analyzed DSSC fly lines was infected. The six Neotropical Drosophila species surveyed, including D. willistoni, D. prosaltans, and D. sturtevanti, were uninfected (4). The Neotropical samples surveyed originated from iso- or oligofemale lines kept at the DSSC since the 1950s. In contrast, Werren et al. (46) reported that a natural population of D. willistoni collected in the early 1990s in Panama was infected with Wolbachia. Its presence in D. willistoni was recently confirmed by discovering partial fragments of a Wolbachia genome in the Trace Archive of the D. willistoni genome sequencing project (37; J. Brownlie, personal communication). The genome sequence was derived from an isofemale line collected in the early 1990s in Guadeloupe (L. Ehrman, personal communication).
![]() View larger version (28K): [in a new window] |
FIG. 1. Phylogenetic relationship of the Sophophora radiation (41). The Wolbachia infection statuses of the Drosophila species shown were deduced from data published previously (melanogaster and obscura groups) and from this study (saltans and willistoni groups). In boldface are species that have been found to be infected with various Wolbachia strains based on wsp sequences (see Materials and Methods).
|
|
|
|---|
|
View this table: [in a new window] |
TABLE 1. Distribution of Wolbachia in natural populations and stocks of D. willistoni
|
|
View this table: [in a new window] |
TABLE 3. Distribution of Wolbachia in the saltans group
|
Phylogenetic analysis.
Multiple wsp sequence alignments, including the hypervariable regions (bases 217 to 252 and 520 to 582), were generated using the Clustal X program (40). Alignments were based on amino acid translations followed by manual modifications. A base substitution was included in the analysis if it occurred in two or more plasmid clones obtained from independent PCRs. Other substitutions were eliminated. The final alignment is available at ftp://ftp.ebi.ac.uk/pub/databases/embl/align/under accession number ALIGN_000917. Phylogenetic trees were constructed by applying PAUP* (39) in the absence of an available outgroup. Neighbor-joining analyses after midpoint rooting and unweighted-pair group method with arithmetic mean analyses yielded similar phylogenies, supporting the close relationship of wAu-like Wolbachia variants.
Single-fly Southern hybridization.
DNA extraction from individual 10-day-old female flies, restriction digestion with HindIII, vertical agarose gel separation, and membrane blotting were performed according to the protocol described by Junakovic (26). Nylon membranes were probed with the eluted wsp PCR fragment of wWil derived from the D. willistoni strain Pan 02 (Table 1) cloned into the pGEM-T Easy vector.
Semiquantitative genomic wsp PCRs.
The density of Wolbachia in D. willistoni was determined by semiquantitative wsp PCRs on 10 individual adult females of staged ages. After gel separation and SYBR Green I staining (Roche), the emission intensities of the obtained wsp fragments were determined and compared to wsp signal intensities derived from individual D. simulans (from a Coffs Harbor line) infected with wAu and D. simulans (from a Riverside line) infected with wRi.
Immunological studies.
Wolbachia density and tissue tropism of wWil in D. willistoni were determined using the polyclonal Wolbachia surface protein (WSP) antibody (11). WSP protein expression was analyzed via Western blotting of protein extracts derived from individual adult flies in independent replicates as well as whole-mount immunostainings on adult tissues and staged embryos (44). Rabbit anti-wsp antibody was used at a 1:500 dilution overnight at 4°C and detected after incubation with a 1:500 dilution of Alexa Fluor 488 goat anti-rabbit immunoglobulin G-labeled secondary antibody (Molecular Probes) at room temperature for 1 h. The total number of primordial germ line cells (PGCs) in stage 10 and later embryos of D. willistoni was determined using the pole cell-specific polyclonal rabbit anti-VASA antibody at a dilution of 1:1,000. Slides were stained for 3 min with 1 µg/ml DAPI (4',6'-diamidino-2-phenylindole) (Molecular Probes), rinsed, stained with 5 µg/ml propidium iodide (Molecular Probes) for 20 min, rinsed again, and mounted with ProLong antifade medium (Molecular Probes).
Fluorescence microscopy.
Immunostainings of embryos and ovaries were examined by using a Zeiss Axiomot 2 Epifluorescence microscope. Images were processed using Photoshop 6.0 (Adobe).
Nucleotide sequence accession numbers.
The wsp sequence data derived from Neotropical Wolbachia strains were deposited in GenBank under accession numbers AY620207 to AY620229 and DQ118779, as well as AY858801 for the respective sequence from D. ananassae collected in 2002 in Sao Tome. Sequences of the diagnostic VNTR-141 loci of D. simulans (Coffs Harbor) and D. willistoni were deposited in GenBank under accession numbers DQ118777 and DQ118778, respectively.
|
|
|---|
![]() View larger version (41K): [in a new window] |
FIG. 2. Intra- and interspecific distributions of Wolbachia in Neotropical Drosophila species. (A) Single-fly wsp PCR on eight strains of D. willistoni collected at different American locations in different years (see Table 1). For each D. willistoni strain tested, PCRs were performed separately on six individual 2-day-old female flies. Lines are as follows: 1, Pan 02; 2, wilC; 3, wilH; 4, wilB6; 5, Apa 5.1; 6, Para 4; 7, WIP4; 8, wilG1. (B) Genomic single-fly Southern blot hybridization probed with the wsp plasmid of wWil on individual 10-day-old females of D. melanogaster/wMel CS (lane 1), D. melanogaster/wMel ywc67 (lane 2), D. simulans/wRi (lane 3), D. simulans/wAu (lanes 4 and 5), D. willistoni/wWil treated with tetracycline (lane 6), and D. willistoni/wWil strain Pan 02 (lane 7). (C) Western immunoblotting using the anti-WSP antibody (1:1,000) on single-fly protein extracts derived from D. willistoni-T, the tetracycline-treated control line of JS 6.3 (lane 1), D. willistoni/wWil (lane 2), D. septentriosaltans/wSpt (lane 3), D. simulans/wAu (lane 4), D. simulans/wRi (lane 5), and D. melanogaster/wMel (lane 6). (D) VNTR-141 specific PCR on D. willistoni/wWil Pan 02 (lane 1), D. willistoni/wWil JS6.3 (lane 2), D. willistoni/wWil Para 4 (lane 3), D. simulans/wAu Coffs Harbor (lane 4), and D. simulans/wAu Yaounde 6 (lane 5). (E) Schematic comparison between the VNTR-141 loci (34) of wWil (top) and wAu (bottom). The basic unit is composed of a 15-bp repeat (stippled), a 23-bp hairpin (loop), an 18-bp insertion (hatched), and a 15-bp repeat (black). The size difference is caused by a 141-bp duplication in VNTR-141 of wAu.
|
Molecular characterization of the D. willistoni-specific Wolbachia strain wWil.
We sequenced fragments of two genes, wsp and ftsZ, from 12 Wolbachia-infected lines covering continental and Caribbean populations of D. willistoni in order to characterize the molecular structure and phylogenetic relationship of this Wolbachia association with other Wolbachia variants. Until recently these two diagnostic marker genes were regarded as the most informative for molecular Wolbachia variant classification (34). All isolated Wolbachia clones of D. willistoni were identical in their sequence. Below we refer to the strain as wWil. With respect to the wsp sequence of wWil obtained from the 12 infected lines (accession numbers AY620218 to AY620229, no sequence polymorphism could be detected. Moreover, all wsp and ftsZ sequences of D. willistoni were 100% identical to the respective wsp and ftsZ genes (accession numbers AF020067 and AY227739) of the Wolbachia variant wAu. As deduced from comparative Southern blots (Fig. 2B), the close relationship between wAu and wWil is corroborated by the conservation of the two HindIII restriction sites flanking the wsp locus.
In contrast to the identity of wWil and wAu at the wsp and ftsZ sequence level, comparative genomic single-fly Southern blots (Fig. 2B) and semiquantitative PCRs (data not shown) of infected individuals of D. willistoni and D. simulans showed clear quantitative differences. Strong signals comparable to those of wRi were obtained from wAu-infected D. simulans adults, and the intensity of wWil in similar-sized D. willistoni clearly showed a 70% reduction compared to that of wAu (Fig. 2B, lanes 4, 5, and 7). This quantitative effect was also detected at the WSP protein expression level by Western blots derived from single-fly protein extracts with the polyclonal anti-wsp antibody (Fig. 2C, lanes 2 and 4). The WSP proteins of wAu and wWil have the same molecular weight, whereas, for example, the homologues of two other Wolbachia variants that infect D. melanogaster and D. simulans (wMel and wRi, respectively) differ significantly (Fig. 2C, lanes 5 and 6).
In contrast to the wsp and ftsZ sequence identity between wWil and wAu, we were able to discriminate both strains at the genomic level by applying the recently isolated polymorphic marker VNTR-141 (36). This diagnostic marker covers the noncoding polymorphic VNTR-141 locus in wMel (positions 89003 to 90332 in the wMel chromosome). By performing VNTR-141-specific PCRs (Fig. 2D), we have obtained a 528-bp fragment from wAu (accession number DQ11877) and a 387-bp fragment from wWil (accession number DQ118778). The length difference is caused by a 141-bp duplication in wAu that is not present in wWil (Fig. 2E). Hence, wWil is closely related but not identical to wAu of D. simulans.
Extreme pole cell tropism of wWil in D. willistoni embryos.
Whole-mount immunostainings were performed on early embryos and ovaries of both fly species, using the anti-WSP antibody. In early embryos of D. simulans, wAu bacteria were detected in somatic and germ line tissues during all stages of embryonic development (Fig. 3A). Nuclei of earlier blastodermal stages were infected with Wolbachia, with some significant enrichment in the posterior pole cell region in both D. simulans and D. willistoni. Such posterior accumulations of wAu in D. simulans blastodermal embryos were reported recently (44). In contrast to this, during stages 9 and 10 of embryonic development of D. willistoni, wWil specifically targets the germ line (Fig. 3B). In early gastrulating embryos, shortly after pole cell invagination, somatic tissues of D. simulans were heavily infected by wAu (Fig. 3C). At this developmental stage, wWil bacteria in D. willistoni are selectively targeting a small number of primordial germ line cells, whereas somatic tissues are devoid of bacteria (Fig. 3D). Later, during stages 12 to 14, in the course of germ band retraction, only one lateral pair of five or six PGC nuclei was infected by wWil (Fig. 3E to H). Control immunostainings with the Drosophila germ line-specific VASA antibody (28) showed that, in contrast to D. melanogaster, the Neotropical species D. willistoni harbors a reduced number of PGCs which perfectly colocalize with wWil (data not shown). Based on the tight temporal and spatial association between the host-encoded VASA protein and WSP-expressing wWil, we assume that this intracellular parasite possesses a molecular association with the host-expressed, pole cell-specific vasa RNA or with its encoded protein.
![]() View larger version (59K): [in a new window] |
FIG. 3. Distribution of Wolbachia in Drosophila embryos. (A and B) Whole-mount immunostainings with rabbit anti-WSP antibody (green) on early-stage embryos of D. simulans infected with wAu (A) and D. willistoni JS 6.3 infected with wWil (B). (C and D) Stage 9 to 10 embryos of D. simulans (C) and D. willistoni (D) are shown in detail. Whereas wAu in D. simulans uniformly infects both somatic and germ line cells, wWil selectively targets a very limited number of primordial germ cells. In D. willistoni, wWil is not detectable at the immunological level in somatic cells at embryonic stage 9 or later. (E to H) Lateral views of a stage 12 embryo of D. willistoni infected with wWil (E and G) and their enhanced magnifications show a row of five heavily infected primordial germ cell nuclei on both lateral sides of the embryo (F and H).
|
|
View this table: [in a new window] |
TABLE 2. Distribution of Wolbachia in the willistoni group
|
![]() View larger version (18K): [in a new window] |
FIG. 4. Unweighted-pair group method with arithmetic mean cladogram based on wsp sequence alignment, including the hypervariable region of the Wolbachia strains (50) derived from Neotropical Drosophila species (boldface) and from earlier reported host species (lightface). Host species found to harbor polymorphic Wolbachia variants are indicated by asterisks.
|
Six wSpt wsp sequences were isolated from three different D. septentriosaltans lines collected in Panama between 1998 and 2002 (Table 3). At least four different wSPT subtypes can be distinguished according to their wsp sequences (Fig. 4): wSPT BCI1 (accession number AY620209) is identical to wCer2 (accession number AF418557) of the cherry fruit fly Rhagoletis cerasi (33) and to wTei (accession number AY291347) and wYak (accession number AY291348) of D. teissieri and D. yakuba, respectively (8). The variant wSpt PLR1 (accession number AY620211) clusters with wSpt PLR2 (accession number AY620212), BCI2 (accession number AY620210), and PNM1 (accession number AY620213). The latter three wsp clones are identical at the sequence level but stem from three different Panamanian D. septentriosaltans populations (Table 3). The fourth subtype, wSpt PNM2 (accession number AY620214), is the most divergent variant positioned between wMel (accession number AF020072) of D. melanogaster and the wAu-like Wolbachia clade (Table 4). All lines of D. septentriosaltans tested are multiply infected with wsp variants of wSpt. For example, individual flies from the PNM strain from Panama City harbor at least two different types of wsp sequences. Each wsp variant sequenced seems to be part of an intact open reading frame encoding a 196-amino-acid (aa) section of the WSP protein. The observed wsp sequence polymorphism of wSpt variants within D. septentriosaltans is manifested even at the protein level (Table 4). With respect to the WSP consensus sequence wBCL1 (accession number AY620209), two amino acid substitutions are found, i.e., in the sequence of wSpt PLR1, PLR2, BCI2, PNM1, and PNM2 at consensus position aa 24 (Tyr to His) and in the variant PLR1 (accession number AY620211) at position 126 (Asp to Gly).
|
View this table: [in a new window] |
TABLE 4. Variable nucleotide and amino acid sites in the wsp sequence of the closely related wAu-like Wolbachia strains of Drosophila
|
wStv Wolbachia in D. sturtevanti.
Our survey yielded another new Wolbachia variant, wStv, which was isolated from D. sturtevanti, a member of the sturtevanti subgroup (Fig. 1). The distribution pattern of the wStv infection within its host species is patchy; e.g., wStv is present in the isofemale line Pan 6 (accession number AY620216) but is absent from Pan 12 (Table 3). As deduced from wsp sequence data wStv belongs to A-group Wolbachia but is distantly related to the wAu-like variants (Fig. 4). Three closely related but distinctive variants of wStv were isolated as singly occurring infections from three Panamanian populations (accession numbers AY620215, AY620216, and AY620217) (Fig. 4). Interestingly, the wsp sequence of wStv MI (accession number AY620215) collected in Maria Eugenia, Panama, is identical to that of wWhi (accession number AF237886) isolated from the phlebotomine sand fly Lutzomyia shannoni in Colombia (31). Those authors proposed, based on an extensive data set showing that other non-American populations of L. shannoni are free of Wolbachia, that L. shannoni probably acquired wWhi recently from another host in America.
|
|
|---|
The complete absence of wAu-like Wolbachia in the related willistoni group species tested (Table 2) supports the idea that D. willistoni was infected after its speciation. Since all recently collected continental samples of D. willistoni are infected with wWil, we assume that this Wolbachia strain reached fixation in continental populations of D. willistoni. wWil's pole cell tropism and its 100% transmission rate, seen in lab lines, were probably crucial factors. A recent Wolbachia-driven process should also be detected in the biogeographic distribution of mitochondrial variation, but this has not yet been looked at in the context of Wolbachia infections. A departure from an expected ratio of mitochondrial versus nuclear DNA polymorphism has been reported when comparing different populations of D. willistoni, and a selective mitochondrial sweep has been suggested as one plausible reason (38; J. Silva and M. Kidwell, personal communication).
wAu-like Wolbachia originated in saltans group species.
We found polymorphic but closely related wsp sequences of wPro and wSpt in the host species D. prosaltans and D. septentriosaltans, respectively. This implies that these Wolbachia variants are an outcome of old associations with Neotropical Drosophila species. Independent multiple horizontal transfers with closely related Wolbachia strains are less likely. The progenitor of wPro and wSpt presumably infected the common ancestor of both host species before speciation and subsequently diverged at the wsp sequence level in the course of long-term vertical transmission. Host-specific diagnostic sites within wsp correspond with our hypothesis (Table 4). Therefore, we suggest that wAu-like variants evolved in the American Neotropical saltans group species and are potential donors for the horizontal transmission to D. willistoni. A similar event has been suggested for Wolbachia associations among the Old World sibling species D. simulans and D. sechellia, where original Wolbachia infections in an original species have not yet yielded a sequence divergence in wsp in the sibling species (6).
Recent horizontal transfer into D. similans: origin of the wAu infection.
Non-CI-inducing wAu of D. simulans (17) is found worldwide, including in Australia, Madagascar, Cameroon, parts of Europe and Japan, Ecuador, Jamaica, and the southern United States (2, 3, 7, 23). The overlapping geographic distribution of populations of D. simulans, D. willistoni, and other Neotropical Drosophila species in Central America, together with wsp and fstZ sequence identity of the two Wolbachia variants wAu and wWil, strongly suggests a recent horizontal transfer of Wolbachia from an original native Neotropical Drosophila-Wolbachia guild to the immigrating Old World species D. simulans. To date D. willistoni can be regarded as the most likely donor species of this transfer. Recent transfers of transposable elements between D. willistoni and another immigrating Old World Drosophila species, D. melanogaster, have been shown for the canonical P transposon (10, 18), and for the retrotransposon copia (13, 25). Furthermore, the male-killing bacterium Spiroplasma poulsonii of the D. willistoni group species D. nebulosa has recently infected immigrating D. melanogaster populations in Brazil (30). Extensive phylogenetic studies of hosts and their parasites suggest horizontal transmission of Wolbachia variants between distantly related insect species (5, 16, 43, 46). Furthermore, it has been experimentally demonstrated that Wolbachia can be shuffled horizontally within and between Trichogramma parasitoid species (19, 20).
In agreement with the hypothesis of an American origin of wAu and opposed to an African origin (7) is the extensive analysis of mitochondrial variation in D. simulans. wAu is globally associated with the mitochondrial-siII haplotype of D. simulans (23). However, some African populations of D. simulans harboring the siII haplotype are uninfected. Ballard proposed recently that uninfected flies migrated to Ecuador and acquired wAu in a horizontal transmission event from an unknown host source (2). Subsequently, wAu spread throughout natural populations of D. simulans worldwide. The infection model outlined by Ballard, based on mitochondrial haplotypes and geographic distribution of wAu-infected D. simulans, is in line with our hypothesis that a Neotropical species such as D. willistoni could be the donor species of wAu.
In summary, we suggest a potential evolutionary scenario: wAu-like variants evolved in the guild of the Neotropical saltans group, being vertically transmitted and/or horizontally shuffled between related host species over a long period of time. More recently, a proto-wAu-like strain, the ancestor of wWil, infected horizontally a locally isolated population of D. willistoni, most likely in Central America. In this population, wWil evolved perfect maternal transmission through an extreme tissue tropism towards the germ line of D. willistoni. Within the last 300 years, immigrating D. simulans flies from Africa may have become infected by wWil or by another wAu-like strain from infected Neotropical Drosophila species through vectors such as parasitoid wasps (20, 43). wAu-infected D. simulans has then spread worldwide (2). An alternative source for wAu is an acquisition from outside the closely related Wolbachia pool of Neotropical Drosophila species, but if so, the fact that Neotropical Drosophila species are infected with closely related Wolbachia strains will need to be explained. It is unclear how wWil and wAu drove themselves through host populations. Presently, neither wAu in D. simulans nor wWil in D. willistoni is able to induce measurable levels of CI (17; W. J. Miller, unpublished data). The possibility that they did so in the past cannot be excluded. As reported by Ballard and coworkers, wAu seems to induce weak levels of CI in some infected populations of D. simulans from Florida (3, 23, 29). Alternatively, the driving force for the spreading of wAu-like strains could be a positive fitness contribution to their hosts that remains to be elucidated. The phenotypes of the Neotropical Wolbachia strains still need to be elucidated. The wAu-like strains are nested within the Mel cluster (50) of closely related Wolbachia strains that have a variety of phenotypic effects in other host species. Based on the wsp sequence, the variant wSpt BCI1 of D. septentriosaltans is identical to the infection of the African species D. yakuba and D. teissieri (27) and to wCer2 of R. cerasi (33). Whereas wCer2 causes strong CI in R. cerasi and in transinfected Ceratitis capitata (49) and intermediate CI in transinfected D. simulans (35), wTei and wYak do not induce CI but are able to fully rescue the wRi mod function in their original host (48). Wolbachia infections of D. melanogaster (32, 50) and the quinaria group member D. recens induce CI (45). The strain wInnA causes male killing in the related D. innubila (22), where it is regarded as an ancestral infection (12).
The present paper shows the complexity of evolutionary dynamics of Wolbachia in Neotropical Drosophila species and its success in colonizing the Old World species D. simulans. Both wWil and wAu successfully colonized natural populations of D. willistoni in America and of D. simulans globally. The detailed understanding of the evolutionary "jump-and-go" dynamics of Wolbachia will have important implications for practical applications of this symbiont as a vector system and in biological pest control management.
We are very much obliged to Petra Zwirn, Edeltraud Kehrer, Ingrid Gerstl, and Sabine Lagger for excellent technical support and assistance. We highly appreciate the valuable material contributions by Margaret Kidwell, Andrew Holyake, Egon Bartel, Kim van der Linde, Peter Chabora, and Francesco Ayala, who donated numerous Neotropical Drosophila lines, and by Paul Lasko and Kostas Bourtzis, who provided the anti-VASA and anti-WSP antibodies, respectively. Moreover, we thank Margaret Kidwell, Joana Silva, Andrew Holyoake, and Filipa Vala for stimulating discussions and Filipa Vala, Jeremy Brownlie, and two anonymous reviewers for valuable comments on earlier versions of the manuscript.
|
|
|---|
This article has been cited by other articles:
| |||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||||
Copyright © 2009 by the American Society for Microbiology. For an alternate route to Journals.ASM.org, visit: http://intl-journals.asm.org | More Info»