Host-related fitness trade-offs in a presumed generalist parasitoid, Diaeretiella rapae (Hymenoptera: Aphidiidae)

Abstract.  1. Host ranges of parasitoid wasps are mediated by behavioural responses to hosts and their environment (infectivity), and development in hosts (virulence). Determinants of host range were measured in Diaeretiella rapae (Hymenoptera: Aphidiidae), which has been described as a generalist that attacks more than 60 species.
2. In northern Colorado, this wasp mainly attacks two hosts: cabbage aphid ( Brevicoryne brassicae ) and Russian wheat aphid ( Diuraphis noxia ). Here, laboratory experiments are described in which D. rapae originating from these two hosts were offered several hosts for oviposition. Both infectivity and virulence were measured.
3. Infectivity included host acceptance and handling time, while virulence was measured as productivity (number of progeny), survival of immatures within hosts, development time, and sex ratio.
4. Wasps had higher productivity and survival when attacking 'home' hosts than 'alternate' hosts, and trade-offs were found by quantitative genetic analyses to be genetically determined. Sex ratio and development times also showed trade-offs, but mainly related to the host environment in which females were reared.
5. In previous genetic studies in northern Colorado, populations were genetically subdivided on the scale of 1 km. The fitness differences described here could be strong enough to create populations adapted to different hosts, but it appears that gene flow is sufficient to prevent formation of separate lineages on the two hosts.
6. Rather than being a generalist with a broad host range, D. rapae is a serial specialist, attacking particular hosts according to availability in different seasons or in different geographical areas.     

Genetic variation in natural populations of the aphid Rhopalosiphum padi as revealed by maternally inherited markers

A survey on 148 clones of the aphid Rhopalosiphum padi from 11 widespread localities has been carried out to study the genetic structure of populations of this species as revealed by mitochondrial DNA restriction site and length polymorphisms as well as by restriction site analysis of a maternally inherited plasmid carried by the aphid eubacterial endosymbiont Buchnera aphidicola. Our results support the existence in the area under study of two main aphid maternal lineages strikingly coincidental with the two main reproductive categories displayed by this species. Those aphid clones possessing an incomplete life cycle that lacks the sexual phase (anholocyclic or androcyclic clones) show mitochondrial DNA (mtDNA) haplotype I and plasmid haplotype I, whereas those clones displaying the complete life cycle (holocyclic clones) posses some other distinct mtDNA haplotypes closely related to each other and plasmid haplotype II. While restriction-site analysis of maternally inherited markers points to a relatively ancient origin of anholocycly/androcycly (between 460 000 and 1 400 000 years) followed by interrupted gene flow with respect to the ancestral holocyclic population, mtDNA size variation also suggests that historical stochastic processes have a different effect on the evolution of both main aphid lineages. Evidence of occasional nuclear gene flow between lineages and its consequences on the correspondence between maternally inherited haplotypes and life cycle are also presented and discussed.     

Horizontal transfer of bacterial symbionts: heritability and fitness effects in a novel aphid host

Members of several bacterial lineages are known only as symbionts of insects and move among hosts through maternal transmission. Such vertical transfer promotes strong fidelity within these associations, favoring the evolution of microbially mediated effects that improve host fitness. However, phylogenetic evidence indicates occasional horizontal transfer among different insect species, suggesting that some microbial symbionts retain a generalized ability to infect multiple hosts. Here we examine the abilities of three vertically transmitted bacteria from the Gammaproteobacteria to infect and spread within a novel host species, the pea aphid, Acyrthosiphon pisum. Using microinjection, we transferred symbionts from three species of natural aphid hosts into a common host background, comparing transmission efficiencies between novel symbionts and those naturally infecting A. pisum. We also examined the fitness effects of two novel symbionts to determine whether they should persist under natural selection acting at the host level. Our results reveal that these heritable bacteria vary in their capacities to utilize A. pisum as a host. One of three novel symbionts failed to undergo efficient maternal transmission in A. pisum, and one of the two efficiently transmitted bacteria depressed aphid growth rates. Although these findings reveal that negative fitness effects and low transmission efficiency can prevent the establishment of a new infection following horizontal transmission, they also indicate that some symbionts can overcome these obstacles, accounting for their widespread distributions across aphids and related insects.     

Horizontal Transfer of Bacterial Symbionts: Heritability and Fitness Effects in a Novel Aphid Host

Members of several bacterial lineages are known only as symbionts of insects and move among hosts through maternal transmission. Such vertical transfer promotes strong fidelity within these associations, favoring the evolution of microbially mediated effects that improve host fitness. However, phylogenetic evidence indicates occasional horizontal transfer among different insect species, suggesting that some microbial symbionts retain a generalized ability to infect multiple hosts. Here we examine the abilities of three vertically transmitted bacteria from the Gammaproteobacteria to infect and spread within a novel host species, the pea aphid, Acyrthosiphon pisum. Using microinjection, we transferred symbionts from three species of natural aphid hosts into a common host background, comparing transmission efficiencies between novel symbionts and those naturally infecting A. pisum. We also examined the fitness effects of two novel symbionts to determine whether they should persist under natural selection acting at the host level. Our results reveal that these heritable bacteria vary in their capacities to utilize A. pisum as a host. One of three novel symbionts failed to undergo efficient maternal transmission in A. pisum, and one of the two efficiently transmitted bacteria depressed aphid growth rates. Although these findings reveal that negative fitness effects and low transmission efficiency can prevent the establishment of a new infection following horizontal transmission, they also indicate that some symbionts can overcome these obstacles, accounting for their widespread distributions across aphids and related insects.     

Host-range evolution in Aphidius parasitoids: fidelity, virulence and fitness trade-offs on an ancestral host

The diversity of parasitic insects remains one of the most conspicuous patterns on the planet. The principal factor thought to contribute to differentiation of populations and ultimately speciation is the intimate relationship parasites share with hosts and the potential for disruptive selection associated with using different host species. Traits that generate this diversity have been an intensely debated topic of central importance to the evolution of specialization and maintenance of ecological diversity. A fundamental hypothesis surrounding the evolution of specialization is that no single genotype is uniformly superior in all environments. This "trade-off" hypothesis suggests that negative fitness correlations can lead to specialization on different hosts as alternative stable strategies. In this study we demonstrate a trade-off in the ability of the parasitoid, Aphidius ervi, to maintain a high level of fitness on an ancestral and novel host, which suggests a genetic basis for host utilization that may limit host-range expansion in parasitoids. Furthermore, behavioral evidence suggests mechanisms that could promote specialization through induced host fidelity. Results are discussed in the context of host-affiliated ecological selection as a potential source driving diversification in parasitoid communities and the influence of host species heterogeneity on population differentiation and local adaptation.     

THE GENETIC STRUCTURE OF HOST PLANT ADAPTATION IN A SPATIAL PATCHWORK: DEMOGRAPHIC VARIABILITY AMONG RECIPROCALLY TRANSPLANTED PEA APHID CLONES

Populations of insect herbivores that feed on several host plant species may experience different selective forces on each host. When the hosts cooccur in a local area, herbivore populations can provide useful models for the study of evolutionary mechanisms in patchy environments. A first step in such a study involves determination of the genetic structure of host adaptation in the region: how is genetic variation for host use structured within and between subpopulations of herbivores on each host? The structure of genetic variation for host use reveals patterns of local adaptation, probable selective consequences of migration between hosts, and the potential for further evolution. To estimate the population structure of host adaptation in a patchwork, 7–11 pea aphid clones were collected at the beginning of the summer from each of two alfalfa and two red clover fields within a very localized area (about 15–20 km²). Using a reciprocal transplant in the field, replicates of these 35 clones were allowed to develop individually on each of the two crops. A complete life table was made for each replicate. Individual fitness was calculated from the life tables as the expected rate of population increase; longevity, age at first reproduction, and total fecundity were also measured for each clonal replicate. Currently, experimental estimates of genetic variation in complete life tables are virtually nonexistent for natural populations, even for single environments (Charlesworth, 1987); field studies are even less common. Because clones from each of two source crops were tested reciprocally on both hosts, variation in relative genotypic fitness on alfalfa and clover could be partitioned among clones within source crops, between fields of the same crop, and between source crops (alfalfa or red clover), providing a view of population structure. Significant clonal variation in relative performance on alfalfa and red clover was found: clones tended to have higher fitness on the crop from which they had been collected (the “home” crop) than they did on the “away” crop, suggesting local adaptation in response to patchy patterns of selection. Clonal variability within collections from the two crops suggests the potential for changes in the genetic constitution of these aphid populations within established fields as a result of clonal selection during the summer season. Significantly negative genetic correlations across crops were found for fitness and its major components. The possibility that these negative cross-environment correlations could act as evolutionary constraints on adaptation to the patchwork is considered.     

Genetic variation and founder effects in the parasitoid wasp, Diaeretiella rapae (M'intosh) (Hymenoptera: Braconidae: Aphidiidae), affecting its potential as a biological control agent

The effects of recent colonization on the aphid parasitoid, Diaeretiella rapae (M'Intosh) (Hymenoptera: Braconidae: Aphidiidae), in Western Australia were investigated. When compared with populations from the Old World, the results of a microsatellite analysis show that the insects have low allelic length and low allele frequency variation, revealing that these individuals experienced a significant founder effect. Marked genetic differentiation between populations was also revealed, which has potentially important implications for host utilization in this species when introduced to a new geographical area(s). Low genetic variation and gene flow in a founder population could limit evolutionary potential in Australia, including the ability of a population to mount a response to newly introduced hosts, such as the Russian wheat aphid, Diuraphis noxia (Mordvilko). Although the actual importance of genetic diversity in the success of biological control agents is unclear, current theory concerning the potential impact of genetic bottlenecks on additive genetic variance is discussed.     

Genetic diversity of the melon aphid <Emphasis Type="Italic">Aphis gossypii</Emphasis> (Hemiptera: Aphididae) in a diversified vegetable growing area

The melon aphid, Aphis gossypii Glover (Hemiptera: Aphididae), is a highly polyphagous species. To investigate its genetic diversity among seasonal host plants, we conducted a field study in a diversified vegetable growing area in Beijing, China. The molecular marker mtDNA COI (1563-bp long) was used to analyze the genetic diversity of aphid populations collected from eight plant species belonging to families Malvaceae, Cucurbitaceae and Rosaceae. A total of 33 haplotypes were identified, five of which were shared haplotypes, while the remaining 28 were unique haplotypes. At least one haplotype was shared by all eight A. gossypii populations. Aphid populations showed high levels of nucleotide and haplotype diversity. The genetic diversity indices were maximal on the hibiscus (the primary host), whereas minimal on cucumber and strawberry (the secondary hosts). The analysis of molecular variance showed that most of the variance was distributed within populations. Based on the genetic distances, the eight populations can be divided into two groups, associated with primary and the secondary host plant, respectively. The aphids possibly migrated from Hibiscus to watermelon, thereafter dispersed from watermelon to other secondary host plants. Watermelon was an important host probably due to its early growing season and large planting area. Our results highlighted the need to target the population dispersal for effective control.     

Population differentiation and genetic variation in performance on eight hosts in the pea aphid complex

Phytophagous insects frequently use multiple host-plant species leading to the evolution of specialized host-adapted populations and sometimes eventually to speciation. Some insects are confronted with a large number of host-plant species, which may provide complex routes of gene flow between host-adapted populations. The pea aphid (Acyrthosiphon pisum) attacks a broad range of plants in the Fabaceae and it is known that populations on Trifolium pratense and Medicago sativa can be highly specialized at exploiting these species. To find out whether adaptation to a broad range of co-occurring hosts has occurred, we tested the performance of pea aphid clones collected from eight host-plant genera on all of these plants in a reciprocal transfer experiment. We provide evidence for pervasive host-plant specialization. The high performance of all aphid clones on Vicia faba suggests that this host plant could be a site of gene flow between different populations that could limit further host-associated divergence. The genetic variance in host-plant usage was partitioned into within- and among-population components, which represent different levels of host adaptation. Little evidence of within-population trade-offs in performance on different plant species was found.     

Rapid evolution and specialization following host colonization in a yucca moth

Host shifts followed by specialization can result in sympatric genetic differentiation, and may have fuelled the diversification of phytophagous insects. This study examines a recent colonization of a non‐native host by Prodoxus quinquepunctellus (Lepidoptera: Prodoxidae). Allozyme differentiation was detected among different host feeding populations, yet was nearly absent among similar host feeding populations in sympatry. Geographical patterns of allozyme variation showed a much higher level of population structure among populations feeding on the derived host. Conversely, mtDNA haplotype frequencies were nearly homogeneous in the derived populations compared to the ancestral populations, suggesting a bottleneck and/or rapid fixation of haplotypes following host colonization. Moth emergence coincided with host plant flowering, and phenological differences between host species translated into allochronic isolation between populations feeding on different hosts. Derived moth populations also differed significantly in three ovipositor characters from ancestral populations. These findings suggest rapid host‐specific genetic differentiation, and specialization of moth emergence time and ovipositor morphology following host colonization.     

Genetic structure and host–parasite co‐divergence: evidence for trait‐specific local adaptation

Host–parasite co‐evolution can lead to genetic differentiation among isolated host–parasite populations and local adaptation between parasites and their hosts. However, tests of local adaptation rarely consider multiple fitness‐related traits although focus on a single component of fitness can be misleading. Here, we concomitantly examined genetic structure and co‐divergence patterns of the trematode Coitocaecum parvum and its crustacean host Paracalliope fluviatilis among isolated populations using the mitochondrial cytochrome oxidase I gene (COI). We then performed experimental cross‐infections between two genetically divergent host–parasite populations. Both hosts and parasites displayed genetic differentiation among populations, although genetic structure was less pronounced in the parasite. Data also supported a co‐divergence scenario between C. parvum and P. fluviatilis potentially related to local co‐adaptation. Results from cross‐infections indicated that some parasite lineages seemed to be locally adapted to their sympatric (home) hosts in which they achieved higher infection and survival rates than in allopatric (away) amphipods. However, local, intrinsic host and parasite characteristics (host behavioural or immunological resistance to infections, parasite infectivity or growth rate) also influenced patterns of host–parasite interactions. For example, overall host vulnerability to C. parvum varied between populations, regardless of parasite origin (local vs. foreign), potentially swamping apparent local co‐adaptation effects. Furthermore, local adaptation effects seemed trait specific; different components of parasite fitness (infection and survival rates, growth) responded differently to cross‐infections. Overall, data show that genetic differentiation is not inevitably coupled with local adaptation, and that the latter must be interpreted with caution in a multi‐trait context.     

Host Plant Specialization in the Sugarcane Aphid Melanaphis sacchari

Most aphids are highly specialized on one or two related plant species and generalist species often include sympatric populations adapted to different host plants. Our aim was to test the hypothesis of the existence of host specialized lineages of the aphid Melanaphis sacchari in Reunion Island. To this end, we investigated the genetic diversity of the aphid and its association with host plants by analyzing the effect of wild sorghum Sorghum bicolor subsp. verticilliflorum or sugarcane as host plants on the genetic structuring of populations and by performing laboratory host transfer experiments to detect trade-offs in host use. Genotyping of 31 samples with 10 microsatellite loci enabled identification of 13 multilocus genotypes (MLG). Three of these, Ms11, Ms16 and Ms15, were the most frequent ones. The genetic structure of the populations was linked to the host plants. Ms11 and Ms16 were significantly more frequently observed on sugarcane, while Ms15 was almost exclusively collected in colonies on wild sorghum. Laboratory transfer experiments demonstrated the existence of fitness trade-offs. An Ms11 isofemale lineage performed better on sugarcane than on sorghum, whereas an Ms15 lineage developed very poorly on sugarcane, and two Ms16 lineages showed no significant difference in performances between both hosts. Both field and laboratory results support the existence of host plant specialization in M. sacchari in Reunion Island, despite low genetic differentiation. This study illustrates the ability of asexual aphid lineages to rapidly undergo adaptive changes including shifting from one host plant to another.     

Influence of host phylogeographic patterns and incomplete lineage sorting on within-species genetic variability in Wigglesworthia species, obligate symbionts of tsetse flies

Vertical transmission of obligate symbionts generates a predictable evolutionary history of symbionts that reflects that of their hosts. In insects, evolutionary associations between symbionts and their hosts have been investigated primarily among species, leaving population-level processes largely unknown. In this study, we investigated the tsetse (Diptera: Glossinidae) bacterial symbiont, Wigglesworthia glossinidia, to determine whether observed codiversification of symbiont and tsetse host species extends to a single host species (Glossina fuscipes fuscipes) in Uganda. To explore symbiont genetic variation in G. f. fuscipes populations, we screened two variable loci (lon and lepA) from the Wigglesworthia glossinidia bacterium in the host species Glossina fuscipes fuscipes (W. g. fuscipes) and examined phylogeographic and demographic characteristics in multiple host populations. Symbiont genetic variation was apparent within and among populations. We identified two distinct symbiont lineages, in northern and southern Uganda. Incongruence length difference (ILD) tests indicated that the two lineages corresponded exactly to northern and southern G. f. fuscipes mitochondrial DNA (mtDNA) haplogroups (P = 1.0). Analysis of molecular variance (AMOVA) confirmed that most variation was partitioned between the northern and southern lineages defined by host mtDNA (85.44%). However, ILD tests rejected finer-scale congruence within the northern and southern populations (P = 0.009). This incongruence was potentially due to incomplete lineage sorting that resulted in novel combinations of symbiont genetic variants and host background. Identifying these novel combinations may have public health significance, since tsetse is the sole vector of sleeping sickness and Wigglesworthia is known to influence host vector competence. Thus, understanding the adaptive value of these host-symbiont combinations may afford opportunities to develop vector control methods.     

Facultative symbiont infections affect aphid reproduction

Some bacterial symbionts alter their hosts reproduction through various mechanisms that enhance their transmission in the host population. In addition to its obligatory symbiont Buchnera aphidicola, the pea aphid Acyrthosiphon pisum harbors several facultative symbionts influencing several aspects of host ecology. Aphids reproduce by cyclical parthenogenesis whereby clonal and sexual reproduction alternate within the annual life cycle. Many species, including the pea aphid, also show variation in their reproductive mode at the population level, with some lineages reproducing by cyclical parthenogenesis and others by permanent parthenogenesis. While the role of facultative symbionts has been well studied during the parthenogenetic phase of their aphid hosts, very little is known on their possible influence during the sexual phase. Here we investigated whether facultative symbionts modulate the capacity to produce sexual forms in various genetic backgrounds of the pea aphid with controlled symbiont composition and also in different aphid genotypes from natural populations with previously characterized infection status and reproductive mode. We found that most facultative symbionts exhibited detrimental effects on their hosts fitness under sex-inducing conditions in comparison with the reference lines. We also showed that the loss of sexual phase in permanently parthenogenetic lineages of A. pisum was not explained by facultative symbionts. Finally, we demonstrated that Spiroplasma infection annihilated the production of males in the host progeny by inducing a male-killing phenotype, an unexpected result for organisms such as aphids that reproduce primarily through clonal reproduction.     

The genetic architecture of a niche: variation and covariation in host use traits in the Colorado potato beetle

The genetic basis of host plant use by phytophagous insects can provide insight into the evolution of ecological niches, especially phenomena such as specialization and phylogenetic conservatism. We carried out a quantitative genetic analysis of multiple host use traits, estimated on five species of host plants, in the Colorado potato beetle, Leptinotarsa decemlineata (Coleoptera: Chrysomelidae). Mean values of all characters varied among host plants, providing evidence that adaptation to plants may require evolution of both behavioral (preference) and post-ingestive physiological (performance) characteristics. Significant additive genetic variation was detected for several characters on several hosts, but not in the capacity to use the two major hosts, a pattern that might be caused by directional selection. No negative genetic correlations across hosts were detected for any 'performance' traits, i.e. we found no evidence of trade-offs in fitness on different plants. Larval consumption was positively genetically correlated across host plants, suggesting that diet generalization might evolve as a distinct trait, rather than by independent evolution of feeding responses to each plant species, but several other traits did not show this pattern. We explored genetic correlations among traits expressed on a given plant species, in a first effort to shed light on the number of independent traits that may evolve in response to selection for host-plant utilization. Most traits were not correlated with each other, implying that adaptation to a novel potential host could be a complex, multidimensional 'character' that might constrain adaptation and contribute to the pronounced ecological specialization and the phylogenetic niche conservatism that characterize many clades of phytophagous insects.     

GEOGRAPHIC DISTRIBUTION OF MITOCHONDRIAL DNA VARIANTS AND THE HISTORICAL BIOGEOGRAPHY OF THE SPOTTED SALAMANDER,AMBYSTOMA MACULATUM

I analyzed geographic partitioning of mitochondrial DNA (mtDNA) restriction-site variants in the spotted salamander, Ambystoma maculatum. Two highly divergent and geographically separate genetic lineages were identified that differed by a minimum of 19 restriction sites (6% sequence divergence). One of the lineages has a disjunct distribution with very closely related haplotypes occurring in Missouri, Arkansas, North Carolina, and Virginia. The other lineage is found in Michigan, Illinois, and Alabama. The geographic separation of highly divergent mtDNA haplotypes, a pattern that was predicted based on the sedentary nature of these salamanders, is evidence for long-term barriers to gene flow. In contrast, the large-scale disjunction of very similar haplotypes suggests recent, long-distance gene flow and does not match the phylogeographic expectation for a small terrestrial vertebrate. I explain this potential contradiction in the level of importance assigned to gene flow by a scenario in which historical barriers to gene flow account for the two divergent mtDNA assemblages, but stochastic sorting of ancestral polymorphism is responsible for the large-scale geographic disjunction. Ten of 16 populations collected in the Ozark Highlands were fixed for the same haplotype. I attribute this lack of detectable variation to recent colonization of this area, a hypothesis that is supported by paleoecological data and demonstrates the potential benefits of combining data from paleobotany, geology, and other disciplines to reconstruct the historical biogeography of a species.     

Mitochondrial DNA analysis in Tibet: Implications for the origin of the Tibetan population and its adaptation to high altitude

Mitochondrial DNAs (mtDNAs) of 54 Tibetans residing at altitudes ranging from 3,000–4,500 m were amplified by polymerase chain reaction (PCR), examined by high-resolution restriction endonuclease analysis, and compared with those previously described in 10 other Asian and Siberian populations. This comparison revealed that more than 50% of Asian mtDNAs belong to a unique mtDNA lineage which is found only among Mongoloids, suggesting that this lineage most likely originated in Asia at an early stage of the human colonization of that continent. Within the Tibetan mtDNAs, sets of additional linked polymorphic sites defined seven minor lineages of related mtDNA haplotypes (haplogroups). The frequency and distribution of these haplogroups in modern Asian populations are supportive of previous genetic evidence that Tibetans, although located in southern Asia, share common ancestral origins with northern Mongoloid populations. This analysis of Tibetan mtDNAs also suggests that mtDNA mutations are unlikely to play a major role in the adaptation of Tibetans to high altitudes. © 1994 Wiley-Liss, Inc.     

Cost of host radiation in an RNA virus

Although host radiation allows a parasite to expand its ecological niche, traits governing the infection of multiple host types can decrease fitness in the original or alternate host environments. Reasons for this reduction in fitness include slower replication due to added genetic material or modifications, fitness trade-offs across host environments, and weaker selection resulting from simultaneous adaptation to multiple habitats. We examined the consequences of host radiation using vesicular stomatitis virus (VSV) and mammalian host cells in tissue culture. Replicate populations of VSV were allowed to evolve for 100 generations on the original host (BHK cells), on either of two novel hosts (HeLa and MDCK cells), or in environments where the availability of novel hosts fluctuated in a predictable or random way. As expected, each experimental population showed a substantial fitness gain in its own environment, but those evolved on new hosts (constant or fluctuating) suffered reduced competitiveness on the original host. However, whereas evolution on one novel host negatively correlated with performance on the unselected novel host, adaptation in fluctuating environments led to fitness improvements in both novel habitats.     

Genotype specificity among hosts,pathogens, and beneficial microbes influences the strength of symbiont‐mediated protection

The microbial symbionts of eukaryotes influence disease resistance in many host‐parasite systems. Symbionts show substantial variation in both genotype and phenotype, but it is unclear how natural selection maintains this variation. It is also unknown whether variable symbiont genotypes show specificity with the genotypes of hosts or parasites in natural populations. Genotype by genotype interactions are a necessary condition for coevolution between interacting species. Uncovering the patterns of genetic specificity among hosts, symbionts, and parasites is therefore critical for determining the role that symbionts play in host‐parasite coevolution. Here, we show that the strength of protection conferred against a fungal pathogen by a vertically transmitted symbiont of an aphid is influenced by both host‐symbiont and symbiont‐pathogen genotype by genotype interactions. Further, we show that certain symbiont phylogenetic clades have evolved to provide stronger protection against particular pathogen genotypes. However, we found no evidence of reciprocal adaptation of co‐occurring host and symbiont lineages. Our results suggest that genetic variation among symbiont strains may be maintained by antagonistic coevolution with their host and/or their host's parasites.