Host-Symbiont Conflicts: Positive Selection on an Outer Membrane Protein of Parasitic but not Mutualistic Rickettsiaceae

Francis M. Jiggins*, Gregory D. D. Hurst{dagger} and Ziheng Yang{dagger}

*Department of Genetics, University of Cambridge;
{dagger}Department of Biology, University College London


    Abstract
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
The Rickettsiaceae is a family of intracellular bacterial symbionts that includes both vertically transmitted parasites that spread by manipulating the reproduction of their host (Wolbachia in arthropods) and horizontally transmitted parasites (represented by Cowdria ruminantium), and mutualists (Wolbachia pipientis in nematode worms). We have investigated the nature of natural selection acting on an outer membrane protein, the wsp gene in Wolbachia and its homologue map1 in Cowdria, thought likely to be involved in host-parasite interactions in these bacteria. The ratio of nonsynonymous to synonymous substitution rates (dN/dS) at individual amino acid sites or at lineages within the gene's phylogeny was estimated using maximum likelihood models of codon substitution. The first hypothesis we tested was that this protein is under positive selection in the parasitic but not in the mutualistic Rickettsiaceae. This hypothesis was supported as positive selection and was detected in Cowdria and arthropod Wolbachia sequence evolution but not in the evolution of Wolbachia sequences from nematodes. Furthermore, this selection was concentrated outside the transmembrane region of the protein and, therefore, in the regions of the protein that may interact with the host. The second hypothesis tested was that positive selection would be stronger in the strains of arthropod Wolbachia that distort the host sex ratio than in those that induce cytoplasmic incompatibility. However, we found no support for this hypothesis. In conclusion, our results are consistent with the hypothesis that antagonistic coevolution causes faster evolution of surface protein sequences in parasites than in mutualists. Confirmation of this conclusion awaits the replication of these results both in additional genes and across more bacterial taxa. The regions of the wsp and map1 genes we identified as likely to be involved in host-parasite arms races should be examined in future studies of parasite virulence and host immune responses, and during the design of vaccines.


    Introduction
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Most animals and plants have intimate associations with a diverse range of microbial symbionts. These symbionts may be classified into mutualists, whose presence increases host fitness, and parasites, which are deleterious. When the association is parasitic, coevolution is expected to take the form of an arms race because host and parasite meet adaptation with counteradaptation (Ehrlich and Raven 1964Citation ). However, when the association is mutualistic, arms races driven by the host evolving to eliminate the pathogen are expected to be absent. This scenario generates the prediction that genes involved in host-symbiont interactions will evolve rapidly because of selection continually favoring novel traits, whereas those involved in mutualistic interactions will tend to be under stabilizing selection. The aim of this study was to test this hypothesis by directly estimating whether natural selection has caused the diversification of an outer membrane protein in parasites and mutualists.

The symbionts studied were intracellular bacterial symbionts of the family Rickettsiaceae. Most of the species in this family are arthropod-vectored parasites of vertebrates, which are represented in our study by the species Cowdria ruminantium, a tick-borne pathogen of ruminants. The second host-symbiont interaction seen in the Rickettsiaceae is reproductive parasitism, which is exhibited by Wolbachia bacteria infecting a range of arthropod hosts. These bacteria are primarily vertically transmitted from mother to offspring, and they spread by manipulating the reproduction of their hosts to enhance their own transmission. The main manipulations are distorting the host sex ratio toward females and inducing cytoplasmic incompatibility (CI) (note that conflict may be restricted to sex-ratio lineages—see subsequently) (Stouthamer, Breeuwer, and Hurst 1999Citation ). The third class of association is mutualism, which is exemplified by the interaction between Wolbachia and their filarial nematode worm hosts. Unlike the parasitic Rickettsiaceae, if the Wolbachia infection of nematode worms is cured with tetracycline, the hosts themselves are harmed (Bandi et al. 1999Citation ). The phylogenies of the bacteria and their nematode hosts are congruent, suggesting that these associations have persisted for far greater periods of time than did the parasitisms (Bandi, Anderson, and Blaxter 1998Citation ).

The arthropod Wolbachia mostly fall into two phenotypic classes, which are intermingled on phylogenies of bacterial genes (Zhou, Rousset, and O'Neill 1998Citation ). In the first type are the sex-ratio distorters, in which the bacterium increases the production of daughters at the expense of sons. This benefits the bacteria, which are transmitted only through females, but is detrimental to the host, whose fitness is usually maximized by investing equally in sons and daughters (Cosmides and Tooby 1981Citation ). Selection, therefore, acts to promote host genes that prevent parasite action and transmission. In the second class of Wolbachia are those that induce CI. These strains of Wolbachia invade the host population by reducing the reproductive success of uninfected females, thereby increasing the frequency of the infected cytoplasm in the population. In diploid hosts they cause a reduction in the viability of eggs when infected males mate with females that do not carry the same strain of Wolbachia. This reduces the fitness of uninfected females in the population, thereby indirectly increasing the relative fitness of the infected females and leading to the spread of the infection through the population (Caspari and Watson 1959Citation ; Turelli and Hoffmann 1991Citation ). The extent of host-parasite conflict is more limited in these associations. Selection on both the bacterium and the female hosts (the transmitting sex) tends to act to increase both the transmission efficiency of the bacterium and the level of compatibility between infected males and females (Turelli 1994Citation ). Therefore, there may be little conflict between CI Wolbachia and their hosts.

Within the Rickettsiaceae, therefore, we would predict that the level of host-symbiont conflict would vary between the different types of symbiosis. In mutualistic interactions there is little conflict, and there is unlikely to be any selection on the host either to suppress the bacterial phenotype or to eliminate the infection. Therefore, comparing the major taxonomic groups, we expect host parasite arms races in the arthropod Wolbachia and C. ruminantium but not in the nematode Wolbachia. At a finer taxonomic scale, comparing the strains of arthropod Wolbachia, we expect a greater amount of host-parasite conflict in sex ratio–distorting strains than in those that induce CI.

These hypotheses can be tested by investigating the nature of selection acting on a bacterial gene that is likely to be involved in host-parasite interactions. We used the outer membrane proteins map1 in C. ruminantium (Vanvliet et al. 1994Citation ) and its homologue wsp in Wolbachia (Braig et al. 1998Citation ). The map1 gene is a member of a multigene family the members of which are closely related to genes found in rickettsial parasites of humans and canines (Sulsona, Mahan, and Barbet 1999Citation ). Phylogenetic analysis of these genes suggests that they have diverged after the split from Wolbachia, although this result relies on the uncertain rooting of the tree (Ohasi et al. 1998Citation ). Outer membrane proteins are likely to be involved in host-parasite interactions, and this has been confirmed in the case of map1, which is known to act as an antigen (Vanvliet et al. 1994Citation ; Perez et al. 1998Citation ). To investigate whether natural selection has favored diversification in these proteins along specific lineages of the bacterial phylogeny, or at specific amino acid sites, we used maximum likelihood analysis of a model of codon substitution to reconstruct the ratio of nonsynonymous to synonymous nucleotide substitutions (the dN/dS ratio) (Yang and Bielawski 2000Citation ). Neutral evolution will produce a dN/dS value of 1, whereas dN/dS < 1 indicates purifying selection, and dN/dS > 1 indicates positive selection. Studies of the wsp gene have shown that it evolves faster in arthropods than in nematodes (Bazzocchi et al. 2000bCitation ), but positive selection was not detected in the arthropod sequences (Zhou, Rousset, and O'Neill 1998Citation ). More recent analyses of a larger data set also found that in the majority of pairwise comparisons dN/dS < 1 (Schulenburg et al. 2000Citation ). However, the pairwise estimation of the average dN/dS ratio along the entire length of a protein is a highly conservative criterion for detecting positive selection.

There were three specific a priori hypotheses tested for these genes. First, the surface protein wsp in C. ruminantium and arthropod Wolbachia is positively selected, whereas the protein in nematode Wolbachia is not. Second, positive selection is stronger on sex ratio–distorting strains than on CI strains. Third, these differences are in regions of the protein that are likely to interact with the host.


    Methods
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Sequence Data
We used published sequences of the wsp gene from Wolbachia and the map1 gene from C. ruminantium. The sequences were used to create four separate alignments (alignment done by eye) containing the nematode Wolbachia (10 sequences), the arthropod Wolbachia (52 sequences), the C. ruminantium (14 sequences), and a fourth alignment containing all the Wolbachia (62 sequences).

In the case of the arthropod Wolbachia, only sequences for which the bacterium's probable phenotype was known were included. In this alignment the codons were numbered according to the reference sequence from strain wRi (accession number AF020070). Because of the missing data, codons 1–36 and 223–229 were excluded from the alignment; of these, codons 1–24 are a signal sequence that is cleaved from the mature peptide (Braig et al. 1998Citation ). In addition, codons 59–65 and 204–213 were excluded from this alignment because they contained numerous insertions and deletions, making the alignment of homologous sites uncertain.

The alignment of the C. ruminantium sequences included all the amino acid sites in the mature peptide, and the codons are numbered according to the reference sequence Welgevonden (accession number U49843). By homology with the Wolbachia sequence, codons 1–25 in the alignment are the signal peptide. The codons of Wolbachia from filarial nematodes were numbered according to the reference sequence from Brugia malayi (accession number AJ252061). The beginning and the end of the full protein sequence were missing from the alignment (codons 1–54 and 203–240). Of these, codons 1–24 are the signal peptide. When both the arthropod and nematode sequences were aligned, it was necessary to exclude slightly different regions from those omitted in the previous alignments. This alignment, numbered according to reference sequence wRi, excluded codons 1–65 and 202–230.

The four alignments have been submitted to the EMBL alignment database under the accession numbers ALIGN_000198–ALIGN_000201. These alignments contain the accession numbers of the individual sequences used.

Tree Reconstruction
The phylogeny of the sequences in each of the four alignments was reconstructed by maximum likelihood using the program PAUP* v.4.0b8 (Swofford 1998Citation ). First, the model of sequence evolution was selected for each alignment by comparing the likelihood scores using the program Modeltest v.3.04 (Posada and Crandall 1998Citation ). The models used were HKY85+G for all the Wolbachia alignments and GTR+I for the C. ruminantium alignment. For each alignment, model parameters were estimated from a maximum parsimony tree and then fixed in the maximum likelihood tree search using the nearest neighbor interchanges. These trees were then used for the estimation of dN/dS by maximum likelihood.

To assess the sensitivity of the results to the tree topology, a maximum parsimony tree was also generated, and all the branches that were supported by less than 50% of 100 bootstrap replicates were collapsed into multifurcations. These trees were used to replicate the results based on the maximum likelihood topologies.

The dN/dS Ratio Across Sites
The first hypothesis that was tested was that sites in the arthropod Wolbachia and C. ruminantium sequences are under positive selection, whereas positive selection is absent in the nematode Wolbachia. A maximum likelihood model of codon substitution was used to identify amino acid sites under positive selection (Nielsen and Yang 1998Citation ; Yang et al. 2000Citation ). The alignments of sequences from arthropod Wolbachia, nematode Wolbachia, and C. ruminantium were analyzed using the codeml program in the PAML package (v3.0a [Yang 1997Citation ]). The distribution of the dN/dS ratio ({omega}) across sites was estimated for the alignments of sequences from arthropod Wolbachia, nematode Wolbachia, and C. ruminantium. Four different models of codon substitution were used (Yang et al. 2000Citation ). M0 (one ratio) assumes that all sites have the same value of {omega}. M3 (discrete) assumes three different classes of sites with different {omega} ratios. M7 (beta) allows sites to have 10 different values of {omega}, calculated from the beta distribution with parameters p and q. The beta distribution is bounded between 0 and 1 and, thus, constitutes a null model for testing positive selection. M8 (beta + {omega}) is similar to M7 but with an additional {omega} category that can exceed 1.

The hypothesis that some amino acid sites are under positive selection can be tested by comparing two nested models with a likelihood ratio test, where these models differ in whether or not they allow some sites to have values of {omega} greater than 1. In this case M7 is compared with M8, and M3 with M0. The former comparison (M7 and M8) is the most stringent test of positive selection (Anisimova, Bielawski, and Yang 2001Citation ). The null distribution of the likelihood ratio test statistic (2{Delta}l, where {Delta}l is the difference between the log-likelihood scores of the two models) can be approximated using the {chi}2 distribution with the degree of freedom (df) being the difference in the number of free parameters between the two models (df = 2 when comparing M7 and M8, and df = 4 when comparing M0 and M3). Analyses of simulated data sets indicate that this test is conservative (Anisimova, Bielawski, and Yang 2001Citation ).

When the likelihood ratio tests suggest the presence of sites under positive selection, an empirical Bayes method is used to calculate the posterior probabilities that each site falls into the {omega} classes (Yang et al. 2000Citation ). Sites with high probabilities of belonging to the class with {omega} > 1 are likely to be under positive selection.

Prediction of Transmembrane Regions
The hypothesis that positive selection results from host-parasite interactions predicts that the positively selected regions of the protein will mostly be in the outer membrane region of the protein. To test this, the transmembrane regions of the protein were predicted from the sequence by using three different methods, DAS (Cserzo et al. 1997Citation ), MemSat (Jones, Taylor, and Thorton 1994Citation ), and TMPred (Hofmann and Stoffel 1993Citation ). To ensure the results, robust predictions were made using all three programs for four arthropod Wolbachia sequences and four C. ruminantium sequences. These sequences were chosen because they covered all, or nearly all, the mature protein.

The dN/dS Ratio Along Lineages
The hypothesis that the {omega} ratio differs between arthropod and nematode sequences was tested directly using a likelihood ratio test comparing models allowing different {omega} ratios along branches within the phylogeny (Yang 1998Citation ; Yang and Nielsen 1998Citation ). The {omega} ratio was estimated using the combined Wolbachia alignment either by forcing all branches to have the same {omega} or by allowing the nematode and arthropod clades to have different {omega}s. The likelihood scores were then compared using a likelihood ratio test, assuming that the test statistic follows a {chi}2 distribution with df = 1.

The final hypothesis tested was that, within the arthropod Wolbachia, the selection on the wsp gene is stronger in sex ratio–distorting strains than in CI strains. The two phenotypes are scattered across the phylogeny of the gene. First, the phenotype of internal branches was reconstructed by maximum parsimony. The value of {omega} was then estimated either by forcing each branch to have the same value of {omega} or by allowing sex-ratio and CI branches to have two different values of {omega}. The hypothesis that {omega} is different between CI and sex-ratio strains was tested by comparing the likelihood scores of the one-ratio and the two-ratio models using a likelihood ratio test, with the test statistic compared with a {chi}2 distribution with df = 1.


    Results
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
Positive Selection in Parasitic and Mutualistic Rickettsiaceae
To illustrate the relationships of the different bacterial strains, the phylogenies of map1 and wsp are shown in figures 1 and 2 . In accordance with our predictions we detected positive selection in the arthropod Wolbachia and the C. ruminantium sequences but not in the nematode Wolbachia. The likelihood scores of the different codon substitution models and sites under selection are shown in tables 1–3 . In the arthropod Wolbachia, likelihood ratio tests indicated that models that included positive selection significantly increased the likelihood scores when compared with models with no positive selection (M0 vs. M3: 2{Delta}l = 565.32, df = 4, P < 0.001; M7 vs. M8: 2{Delta}l = 91.62, df = 2, P < 0.001). Similarly, these tests supported the presence of positive selection in the C. ruminantium sequences (M0 vs. M3: 2{Delta}l = 185.5, df = 4, P < 0.001; M7 vs. M8: 2{Delta}l = 7.7, df = 2, P < 0.05). These results contrast with those from the nematode Wolbachia, in which no sites under positive selection were detected (table 2 ). The difference in the dN/dS ratio at homologous sites in the wsp gene from arthropod and nematode Wolbachia is shown graphically in figure 3 .



View larger version (26K):
[in this window]
[in a new window]
 
Fig. 1.—The maximum likelihood phylogeny of the wsp gene of Wolbachia bacteria from both arthropod and nematode hosts. The taxa are labeled with the name of the host organism or the bacterial strain, and the scale refers to the mean number of substitutions per nucleotide site

 


View larger version (16K):
[in this window]
[in a new window]
 
Fig. 2.—The maximum likelihood phylogeny of the map1 gene from Cowdria ruminantium. The taxa are labeled with the names of the strains, and the scale refers to the mean number of substitutions per nucleotide site

 

View this table:
[in this window]
[in a new window]
 
Table 1 Positively Selected Sites, Log-Likelihood Scores, and Parameter Estimates for the wsp Gene from Arthopod Wolbachia

 

View this table:
[in this window]
[in a new window]
 
Table 2 Positively Selected Sites, Log-Likelihood Scores, and Parameter Estimates for the wsp Gene from Nematode Wolbachiaa

 


View larger version (54K):
[in this window]
[in a new window]
 
Fig. 3.—The predicted dN/dS ratio at homologous amino acid sites in the wsp gene of the parasitic arthropod Wolbachia (top) and the mutualistic nematode Wolbachia (bottom). This is calculated as the average of the {omega} ratios over the 11 classes from model M8, with the posterior probabilities used as weights. This is an approximation of the posterior mean of {omega} for the site. Homologous amino acid sites were aligned for the figure using the program SIM (Huang and Miller 1991Citation ). Also marked are the putative transmembrane regions of the arthropod protein (estimated by DAS) and the regions of the sequence that were missing from the analysis, both of which are described in the text

 
We directly tested whether the value of {omega} averaged across sites was greater in the arthropod than in the nematode sequences by comparing the likelihood score of a model in which nematode and arthropod lineages had the same {omega} against a model in which {omega} differed between the two clades. The two-ratio model estimated {omega} to be greater in arthropods than in nematodes (table 4 ), and this model significantly increased the likelihood score relative to the one-ratio model (2{Delta}l = 15.86, df = 1, P < 0.001).


View this table:
[in this window]
[in a new window]
 
Table 4 Differences in the dN/dS Ratio of the wsp Gene in Nematode Wolbachia (mutualists) and Arthropod Wolbachia (parasites) Sequences

 
Finally, we investigated the robustness and consistency of these results. First, we analyzed the sequences from the A and B groups (see fig. 1 ) of arthropod Wolbachia separately. In both cases positive selection was detected both by comparing M7 versus M8 and M0 versus M3 (P < 0.001 in all cases; data not shown). Furthermore, similar regions of the protein were under positive selection in both groups; the eight positively selected codons in the A group were also detected in the B group. A smaller number of positively selected sites was detected in the A group (8) compared with the B group (20) or with the combined data (20). This is probably because of the high estimate of {omega} in the A-group analysis ({omega} = 9), meaning that fewer sites were assigned to this category. Second, we investigated the sensitivity of our results to the tree topology by replicating the analysis using a maximum parsimony bootstrap consensus tree. The results under this tree topology were virtually identical to those obtained using the maximum likelihood tree. Other factors such as recombination or concerted evolution may further confound the phylogeny, and therefore we replicated the analysis for a third time using a star phylogeny. Again, positive selection was detected in the parasitic but not in the mutualistic data sets. Third, different regions of the protein were excluded from the analysis of the arthropod and nematode wsp sequences (fig. 1 ). We repeated the analysis including only homologous codons in both data sets. Again, this had little effect on the results. Finally, the nematode Wolbachia wsp data set had fewer taxa and a shorter tree length. Therefore, we repeated the analysis on four alignments of between seven and ten arthropod Wolbachia sequences with a tree length similar to that of the nematode data set. In all four replicates, positive selection was detected.

Prediction of Transmembrane Regions
Two transmembrane regions were consistently predicted for each of the four arthropod wsp sequences. In the reference sequence (wRi) these were inferred to be at codons 117–125 and 140–147 by DAS, at 111–127 and 138–150 by MemSat, and at 111–130 and 136–153 by TMPred. In the other three sequences examined, at least two of the three programs predicted these two transmembrane regions. Furthermore, these regions were similar to those predicted in an earlier study (Braig et al. 1998Citation ). The methods failed to produce a consistent prediction as to the orientation of the transmembrane regions. Consistent with our hypothesis, the transmembrane region is not subject to positive selection (fig. 1 ). The predicted transmembrane regions for the four C. ruminantium sequences used were not consistent across methods and sequences, and these results were, therefore, discarded.

Selection Pressure on Sex Ratio–Distorting and CI Strains
The second hypothesis tested was whether positive selection was stronger along sex-ratio Wolbachia lineages than along CI Wolbachia lineages. First, internal branches were assigned to either phenotype by maximum parsimony. The value of {omega} was then estimated using a two-ratio model in which {omega} differed between sex-ratio and CI branches. Contrary to our prediction, {omega} was estimated to be slightly higher in the CI branches ({omega}CI = 0.37 vs. {omega}SR = 0.33; table 5 ). When this model is compared with the one-ratio model, there is no evidence of any difference in the dN/dS ratio along CI and sex-ratio lineages (2{Delta}l = 0.32, df = 1, not significant).


View this table:
[in this window]
[in a new window]
 
Table 5 Differences in the dN/dS Ratio of the wsp Gene of Sex Ratio–Distorting and Cytoplasmic incompatibility (CI) Wolbachia

 

    Discussion
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
We have tested two a priori hypotheses. The predictions of the first, that host-symbiont conflicts are stronger in the parasitic than in the mutualistic Rickettsiaceae, were supported. However, the predictions of the second, that conflict is greater in sex ratio–distorting strains than in CI strains, were not supported.

Parasites and Mutualists
The first hypothesis was that antagonistic coevolution between the host and the parasitic Rickettsiaceae should exert stronger selection on the symbionts when compared with mutualistic relationships. As predicted, the outer membrane protein studied was under positive selection in the parasites but not in the mutualists. Furthermore, sites predicted to be under selection fell outside the predicted transmembrane regions and, hence, possibly in regions of the protein exposed to the host. The importance of the map1 gene in host-parasite interactions has been confirmed experimentally in C. ruminantium, where variation in the sequence of this protein has been shown to correlate with variable antigenic properties (Perez et al. 1998Citation ).

An important limitation of this analysis is that we have analyzed data only from a single gene in two parasitic and in one mutualistic bacterial clades. This clearly raises the possibility that some other factors have confounded our analysis, such as a change in protein function or changes in the ecology of the bacteria other than between parasitism and mutualism. Therefore, the conclusion that parasitism but not mutualism selects for accelerated change in protein sequence needs to be replicated both across additional genes and for further bacterial symbioses. The conclusion that the positive selection pressure is from the host immune system will need to be confirmed by functional studies because convincing support for this idea is available only for C. ruminantium (Perez et al. 1998Citation ). The possible interactions between the arthropod host and Wolbachia remain a matter for speculation. However, surface proteins may play a role in recognition by the host and also potentially as targets of the host in attempts to resist parasite proliferation or the entry of the bacteria into the egg.

It has been suggested that map1 could be a candidate for vaccine development (Perez et al. 1998Citation ). Our analysis has identified sites likely to be under positive selection. One simple explanation is that those sites are involved in antigenic variation in map1 and, thus, evolve at fast rates because of the selective pressure from the host immune response. Identification of those sites may thus provide insights into future efforts of vaccine development. Similarly, the positively selected sites in the wsp gene (the map1 homologue) of arthropod Wolbachia are also likely to be important in host-parasite interactions. Our structural predictions located transmembrane regions of the protein, which are predicted to be under stabilizing or purifying selection in the maximum likelihood analysis. The structural prediction and the maximum likelihood comparison are thus consistent as we expect the transmembrane regions to be conserved.

The filarial nematode worms infected with mutualistic Wolbachia are themselves parasites of vertebrates. However, because the bacteria are within the cells of the worm, it is perhaps less likely (although not impossible) that the vertebrate immune system will exert any selection pressure on the outer membrane proteins of Wolbachia. Despite this, cats infected with filarial nematodes produced antibodies to the wsp protein (Bazzocchi et al. 2000aCitation ). There may be an immune response to Wolbachia outside the worm that are probably released from dead nematodes. Indeed, treatment of filarial nematode infections is known to cause an inflammatory response in mice through the release of Wolbachia from dead worms (Taylor, Cross, and Bilo 2000Citation ). The absence of positive selection on the wsp gene of these bacteria is consistent with the hypothesis that this vertebrate immune response does not result in the elimination of the nematode worms. Testing of this hypothesis is timely.

Sex-Ratio Distortion and CI
The second hypothesis was that, within the arthropod Wolbachia, the selection on sex-ratio strains would be greater than on strains that induce CI. This is expected because, all else being equal, infection with a sex-ratio distorter is deleterious. In contrast, there may be a cost for being an uninfected female in a population infected with CI Wolbachia because such females cannot mate successfully with infected males. However, this hypothesis was not supported by our analysis because we could not detect any difference in the dN/dS ratio in sex-ratio and CI strains. A possible factor limiting the power of our test is the observation that strains may switch fairly frequently between sex-ratio distortion and CI (F. Jiggins, unpublished data). If this is the case, then our assignment of internal and even terminal branches of the tree to a particular phenotype may have been inaccurate.

However, if selection on CI and sex-ratio strains is similar, then positive selection on the arthropod wsp sequences may result from aspects of the parasitic lifestyle that are held in common by both the phenotypic classes. In particular, parasitism is associated with frequent shifts between host species, whereas mutualism leads to long-term stable relationships involving cocladogenesis of host and bacterium (Bandi, Anderson, and Blaxter 1998Citation ; Zhou, Rousset, and O'Neill 1998Citation ). Therefore, continual adaptation to novel hosts may, in part, account for the positive selection on wsp in arthropod Wolbachia. It is noteworthy that the host range of C. ruminantium is probably intermediate between the two Wolbachia data sets.

The validity of all these conclusions depends on our ability to accurately estimate the dN/dS ratio. One point of concern is that the phylogeny of the wsp gene is known not to be very robust (Schulenburg et al. 2000Citation ), and factors such as recombination within the gene may confound the reconstruction of the tree (Jiggins et al. 2001Citation ; Werren and Bartos 2001Citation ). However, this problem does not appear serious because the method used is known to be insensitive to the tree topology (Yang et al. 2000Citation ; Swanson et al. 2001Citation ). This was confirmed for our data set, where different trees all produced similar results.


View this table:
[in this window]
[in a new window]
 
Table 3 Positively Selected Sites, Log-Likelihood Scores, and Parameter Estimates for the map1 Gene from Cowdria ruminantium

 

    Acknowledgements
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 
We thank two anonymous reviewers and the associate editor for their comments on this manuscript. F.M.J. is funded by a Research Fellowship from Emmanuel College, Cambridge, and a Natural Environment Research Council small grant. G.D.D.H. is supported by a Biotechnology and Biological Sciences Research Council D. Phillips Fellowship. Z.Y. is supported by a Biotechnology and Biological Sciences Research Council grant and an EU Human Frontier Science Programme grant.


    Footnotes
 
Pekka Pamilo, Reviewing Editor

Abbreviation: CI, cytoplasmic incompatibility. Back

Keywords: coevolution Cowdria map1 gene positive selection Rickettsiaceae Wolbachia wsp gene Back

Address for correspondence and reprints: Francis M. Jiggins, Department of Genetics, University of Cambridge, Downing Street, Cambridge CB2 3EH, United Kingdom. fmj1001{at}mole.bio.cam.ac.uk Back


    References
 TOP
 Abstract
 Introduction
 Methods
 Results
 Discussion
 Acknowledgements
 References
 

    Anisimova M., J. P. Bielawski, Z. Yang, 2001 Accuracy and power of the likelihood ratio test in detecting adaptive molecular evolution Mol. Biol. Evol 18:1585-1592[Abstract/Free Full Text]

    Bandi C., C. G. Anderson, M. L. Blaxter, 1998 Phylogeny of Wolbachia in filarial nematodes Proc. R. Soc. Lond. B 265:2407-2413[ISI][Medline]

    Bandi C., J. W. Mccall, C. Genchi, S. Corona, L. Venco, L. Sacchi, 1999 Effects of tetracycline on the filarial worms Brugia pahangi and Dirofilaria immitis and their bacterial endosymbionts Wolbachia Int. J. Parasitol 29:357-364[ISI][Medline]

    Bazzocchi C., F. Ceciliani, J. W. Mccall, I. Ricci, C. Genchi, C. Bandi, 2000a. Antigenic role of the endosymbionts of filarial nematodes: IgG response against the Wolbachia surface protein in cats infected with Dirofilaria immitis Proc. R. Soc. Lond. Ser. B Biol. Sci 267:2511-2516[ISI][Medline]

    Bazzocchi C., W. Jamnongluk, S. L. O'neill, T. J. C. Anderson, C. Genchi, C. Bandi, 2000b. wsp Gene sequences from the Wolbachia of filarial nematodes Curr. Microbiol 41:96-100[ISI][Medline]

    Braig H. R., W. Zhou, S. L. Dobson, S. L. O'neill, 1998 Cloning and characterization of a gene encoding the major surface protein of the bacterial endosymbiont Wolbachia pipientis J. Bacteriol 180:2373-2378[Abstract/Free Full Text]

    Caspari E., G. Watson, 1959 On the evolutionary importance of cytoplasmic sterility in mosquitoes Evolution 13:568-570[ISI]

    Cosmides L. M., J. Tooby, 1981 Cytoplasmic inheritance and intragenomic conflict J. Theor. Biol 89:83-129[ISI][Medline]

    Cserzo M., E. Wallin, I. Simon, G. Vonheijne, A. Elofsson, 1997 Prediction of transmembrane alpha-helices in prokaryotic membrane proteins: the dense alignment surface method Protein Eng 10:673-676[Abstract]

    Ehrlich P. R., P. H. Raven, 1964 Butterflies and plants: a study in coevolution Evolution 18:586-608[ISI]

    Hofmann K., W. Stoffel, 1993 TMbase—a database of membrane spanning proteins segments Biol. Chem. Hoppe-Seyler 374:166.

    Huang X. Q., W. Miller, 1991 A time-efficient, linear-space local similarity algorithm Adv. Appl. Math 12:337-357[ISI]

    Jiggins F. M., J. H. G.V. D. Schulenburg, G. D. D. Hurst, M. E. N. Majerus, 2001 Recombination confounds interpretations of Wolbachia evolution Proc. R. Soc. Lond. B 268:1423-1427[ISI][Medline]

    Jones D. T., W. R. Taylor, J. M. Thorton, 1994 A model recognition approach to the prediction of all-helical membrane-protein structure and topology Biochemistry 33:3038-3049[ISI][Medline]

    Nielsen R., Z. Yang, 1998 Likelihood models for detecting positively selected amino acid sites and applications to the HIV-1 envelope gene Genetics 148:929-936[Abstract/Free Full Text]

    Ohasi N., A. Unver, N. Zhi, Y. Rikihisa, 1998 Cloning and characterization of multigenes encoding the immunodominant 30-kilodalton major outer membrane proteins of Ehlichia canis and application of the recombinant protein for serodiagnosis J. Clin. Microbiol 36:2671-2680[Abstract/Free Full Text]

    Perez J. M., D. Martinez, C. Sheikboudou, F. Jongejan, A. Bensaid, 1998 Characterization of variable immunodominant antigens of Cowdria ruminantium by ELISA and immunoblots Parasite Immunol 20:613-622[ISI][Medline]

    Posada D., K. Crandall, 1998 Modeltest: testing the model of DNA substitution Bioinformatics 14:817-818[Abstract]

    Schulenburg J. H.G. V. D., G. D. D. Hurst, T. M. E. Huigens, M. M. M. Van Meer, F. M. Jiggins, M. E. N. Majerus, 2000 Molecular evolution and phylogentic utility of WolbachiaftsZ and wsp gene sequences with special reference to the origin of male-killing Mol. Biol. Evol 17:584-600[Abstract/Free Full Text]

    Stouthamer R., J. A. J. Breeuwer, G. D. D. Hurst, 1999 Wolbachia pipientis: microbial manipulator of arthropod reproduction Annu. Rev. Microbiol 53:71-102[ISI][Medline]

    Sulsona C. R., S. M. Mahan, A. F. Barbet, 1999 The map1 gene of Cowdria ruminantium is a member of a multigene family containing both conserved and variable genes Biochem. Biophys. Res. Commun 257:300-305[ISI][Medline]

    Swanson W. J., Z. Yang, M. F. Wolfner, C. F. Aquadro, 2001 Positive Darwinian selection drives the evolution of several female reproductive proteins in mammals Proc. Natl. Acad. Sci. USA 98:2509-2514[Abstract/Free Full Text]

    Swofford D. L., 1998 PAUP*: phylogenetic analysis using parsimony (* and other methods). Version 4 Sinauer Associates, Sunderland, Mass

    Taylor M. J., H. F. Cross, K. Bilo, 2000 Inflammatory responses induced by the filarial nematode Brugia malayi are mediated by lipopolysaccharide-like activity from endosymbiotic Wolbachia bacteria J. Exp. Med 191:1429-1435[Abstract/Free Full Text]

    Turelli M., 1994 Evolution of incompatibility-inducing microbes and their hosts Evolution 48:1500-1513[ISI]

    Turelli M., A. A. Hoffmann, 1991 Rapid spread of an inherited incompatibility factor in California Drosophila Nature 353:440-442[ISI][Medline]

    Vanvliet A. H. M., F. Jongejan, M. Vankleef, B. A. M. Vanderzeijst, 1994 Molecular-cloning, sequence-analysis, and expression of the gene encoding the immunodominant 32-kilodalton protein of Cowdria-ruminantium Infect. Immun 62:1451-1456.[Abstract]

    Werren J. H., J. D. Bartos, 2001 Recombination in Wolbachia Curr. Biol 11:431-435[ISI][Medline]

    Yang Z., 1997 PAML: a program package for phylogenetic analysis by maximum likelihood Comput. Appl. Biosci 13:555-556[Medline]

    ———. 1998 Likelihood ratio tests for detecting positive selection and application to primate lysozyme evolution Mol. Biol. Evol 15:568-573[Abstract]

    Yang Z. H., J. P. Bielawski, 2000 Statistical methods for detecting molecular adaptation Trends Ecol. Evol 15:496-503[ISI][Medline]

    Yang Z., R. Nielsen, 1998 Synonymous and nonsynonymous rate variation in nuclear genes of mammals J. Mol. Evol 46:409-418[ISI][Medline]

    Yang Z. H., R. Nielsen, N. Goldman, A. M. K. Pedersen, 2000 Codon-substitution models for heterogeneous selection pressure at amino acid sites Genetics 155:431-449[Abstract/Free Full Text]

    Zhou W. F., F. Rousset, S. O'neill, 1998 Phylogeny and PCR based classification of Wolbachia strains using wsp gene sequences Proc. R. Soc. Lond. B 265:509-515[ISI][Medline]

Accepted for publication April 12, 2002.