The role of defensive symbionts in host–parasite coevolution

Understanding the coevolution of hosts and parasites is a long‐standing goal of evolutionary biology. There is a well‐developed theoretical framework to describe the evolution of host–parasite interactions under the assumption of direct, two‐species interactions, which can result in arms race dynamics or sustained genotype fluctuations driven by negative frequency dependence (Red Queen dynamics). However, many hosts rely on symbionts for defence against parasites. Whilst the ubiquity of defensive symbionts and their potential importance for disease control are increasingly recognized, there is still a gap in our understanding of how symbionts mediate or possibly take part in host–parasite coevolution. Herein we address this question by synthesizing information already available from theoretical and empirical studies. First, we briefly introduce current hypotheses on how defensive mutualisms evolved from more parasitic relationships and highlight exciting new experimental evidence showing that this can occur very rapidly. We go on to show that defensive symbionts influence virtually all important determinants of coevolutionary dynamics, namely the variation in host resistance available to selection by parasites, the specificity of host resistance, and the trade‐off structure between host resistance and other components of fitness. In light of these findings, we turn to the limited theory and experiments available for such three‐species interactions to assess the role of defensive symbionts in host–parasite coevolution. Specifically, we discuss under which conditions the defensive symbiont may take over from the host the reciprocal adaptation with parasites and undergo its own selection dynamics, thereby altering or relaxing selection on the hosts' own immune defences. Finally, we address potential effects of defensive symbionts on the evolution of parasite virulence. This is an important problem for which there is no single, clear‐cut prediction. The selection on parasite virulence resulting from the presence of defensive symbionts in their hosts will depend on the underlying mechanism of defence. We identify the evolutionary predictions for different functional categories of symbiont‐conferred resistance and we evaluate the empirical literature for supporting evidence. We end this review with outstanding questions and promising avenues for future research to improve our understanding of symbiont‐mediated coevolution between hosts and parasites.     

The natural occurrence of secondary bacterial symbionts in aphids

1. Many insects host secondary bacterial symbionts that are known to have wide‐ranging effects on their hosts, from host‐plant use to resistance against natural enemies. This has been most widely studied in aphids, which have become a model system to study insect–bacteria interactions. 2. While there is an increasing understanding of the role of symbionts in aphids from controlled laboratory studies, we are only beginning to explore the impact of hosting these symbionts on eco‐evolutionary dynamics in natural systems. To date, many research groups have identified bacterial symbionts from various aphid species, providing us with a bank of literature on aphid–symbiont associations in natural populations. 3. The role of secondary symbionts in aphids is discussed, and the taxonomic and geographical distribution of symbionts among aphids are summarised, and the potential reasons for the patterns observed. The need to test for multiple symbiont species (and co‐infections) across many individuals and the whole distribution range of an aphid is highlighted, including sampling on all known host‐plant species. 4. It is further important also to consider variation within the symbiont, the aphid‐host and the surrounding community, e.g. host‐plants or the natural enemies, to understand how these have the potential to mediate aphid–symbiont interactions. 5. Finally, the knowledge gained from experimental work should now be used to understand the role of aphid secondary symbionts in field systems, to fully understand the potentially far‐reaching consequences of aphid endosymbionts on community and ecosystem processes.     

Phylogeographic patterns in the Philippine archipelago influence symbiont diversity in the bobtail squid–Vibrio mutualism

Marine microbes encounter a myriad of biotic and abiotic factors that can impact fitness by limiting their range and capacity to move between habitats. This is especially true for environmentally transmitted bacteria that cycle between their hosts and the surrounding habitat. As geologic history, biogeography, and other factors such as water temperature, salinity, and physical barriers can inhibit bacterial movement to novel environments, we chose to examine the genetic architecture of Euprymna albatrossae (Mollusca: Cephalopoda) and their Vibrio fischeri symbionts in the Philippine archipelago using a combined phylogeographic approach. Eleven separate sites in the Philippine islands were examined using haplotype estimates that were examined via nested clade analysis to determine the relationship between E. albatrossae and V. fischeri populations and their geographic location. Identical analyses of molecular variance (AMOVA) were used to estimate variation within and between populations for host and symbiont genetic data. Host animals demonstrated a significant amount of variation within island groups, while symbiont variation was found within individual populations. Nested clade phylogenetic analysis revealed that hosts and symbionts may have colonized this area at different times, with a sudden change in habitat. Additionally, host data indicate restricted gene flow, whereas symbionts show range expansion, followed by periodic restriction to genetic flow. These differences between host and symbiont networks indicate that factors “outside the squid” influence distribution of Philippine V. fischeri. Our results shed light on how geography and changing environmental factors can impact marine symbiotic associations at both local and global scales.     

Detecting ancient codispersals and host shifts by double dating of host and parasite phylogenies: Application in proctophyllodid feather mites associated with passerine birds

Inferring cophylogeographic events requires matching the timing of these events on both host and symbiont (e.g., parasites) phylogenies because divergences of hosts and their symbionts may not temporally coincide, and host switches may occur. We investigate a large radiation of birds (Passeriformes) and their permanent symbionts, the proctophyllodid feather mites (117 species from 116 bird species; six genes, 11,468 nt aligned) using two time‐calibration strategies for mites: fossils only and host phylogeography only. Out of 10 putative cophylogeographic events 4 agree in timing for both symbiont and host events being synchronous co‐origins or codispersals; three were based on host shifts, but agree in timing being very close to the origin of modern hosts; two disagree; and one large basal mite split was seemingly independent from host phylogeography. Among these events was an ancient (21–25.3 Mya), synchronous codispersal from the Old World leading to the origin and diversifications of New World emberizoid passerids and their mites, the thraupis + quadratus species groups of Proctophyllodes. Our framework offers a more robust detection of host and symbiont cophylogeographic events (as compared to host‐symbiont reconciliation analysis and using host phylogeography for time‐calibration) and provides independent data for testing alternative hypotheses on timing of host diversification and dispersal.     

The more,the merrier? Obligate symbiont density changes over time under controlled environmental conditions,yet holds no clear fitness consequences

Symbiotic bacteria are highly diverse, play an important role in ecology and evolution, and are also of applied relevance because many pest insects rely on them for their success. However, the dynamics and regulation of symbiotic bacteria within hosts is complex and still poorly understood outside of a few model systems. One of the most intriguing symbiotic relationships is the obligate, tripartite nutritional mutualism in sap‐feeding, economically‐destructive mealybugs (Hemiptera: Sternorrhyncha: Pseudococcidae), which involves γ‐proteobacteria hosted within β‐proteobacteria hosted within the mealybugs. The present study examines whether there is population variation in symbiont density (i.e. infection intensity, or titre) in the citrus mealybug Planococcus citri (Risso) and how this impacts host life‐history. Symbiont density is found to differ significantly between populations when reared under controlled environmental conditions, indicating that the density of symbiont infections is influenced by host or symbiont genotype. However, symbiont density changes in populations over multiple generations, indicating that symbiont densities are dynamic. Surprisingly, given that the symbionts are essential nutritional mutualists, the density of the symbionts does not correlate significantly with either host fecundity or development. Higher levels of symbionts have no clear benefit to hosts and therefore appear to be superfluous, at least under constant, optimized environmental conditions. Excessive symbiont density may be an evolutionary artefact from a period of inefficient vertical transmission when the balance of conflict between host and symbiont was still being established.     

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.     

Landscape genomics of an obligate mutualism: Concordant and discordant population structures between the leafcutter ant Atta texana and its two main fungal symbiont types

To explore landscape genomics at the range limit of an obligate mutualism, we use genotyping‐by‐sequencing (ddRADseq) to quantify population structure and the effect of host–symbiont interactions between the northernmost fungus‐farming leafcutter ant Atta texana and its two main types of cultivated fungus. Genome‐wide differentiation between ants associated with either of the two fungal types is of the same order of magnitude as differentiation associated with temperature and precipitation across the ant's entire range, suggesting that specific ant–fungus genome–genome combinations may have been favoured by selection. For the ant hosts, we found a broad cline of genetic structure across the range, and a reduction of genetic diversity along the axis of range expansion towards the range margin. This population‐genetic structure was concordant between the ants and one cultivar type (M‐fungi, concordant clines) but discordant for the other cultivar type (T‐fungi). Discordance in population‐genetic structures between ant hosts and a fungal symbiont is surprising because the ant farmers codisperse with their vertically transmitted fungal symbionts. Discordance implies that (a) the fungi disperse also through between‐nest horizontal transfer or other unknown mechanisms, and (b) genetic drift and gene flow can differ in magnitude between each partner and between different ant–fungus combinations. Together, these findings imply that variation in the strength of drift and gene flow experienced by each mutualistic partner affects adaptation to environmental stress at the range margin, and genome–genome interactions between host and symbiont influence adaptive genetic differentiation of the host during range evolution in this obligate mutualism.     

Evolution of transmission mode in conditional mutualisms with spatial variation in symbiont quality

Many symbioses have costs and benefits to their hosts that vary with the environmental context, which itself may vary in space. The same symbiont may be a mutualist in one location and a parasite in another. Such spatially conditional mutualisms pose a dilemma for hosts, who might evolve (higher or lower) horizontal or vertical transmission to increase their chances of being infected only where the symbiont is beneficial. To determine how transmission in hosts might evolve, we modeled transmission evolution where the symbiont had a spatially conditional effect on either host lifespan or fecundity. We found that over ecological time, symbionts that affected lifespan but not fecundity led to high frequencies of infected hosts in areas where the symbiont was beneficial and low frequencies elsewhere. In response, hosts evolved increased horizontal transmission only when the symbiont affected lifespan. We also modeled transmission evolution in symbionts, which evolved high horizontal and vertical transmission, indicating a possible host–symbiont conflict over transmission mode. Our results suggest an eco‐evolutionary feedback where the component of host fitness affected by a conditionally mutualistic symbiont in turn determines its distribution in the population, and, through this, the transmission mode that evolves.     

Population dynamics can be more important than physiological limits for determining range shifts under climate change

Evidence is accumulating that species' responses to climate changes are best predicted by modelling the interaction of physiological limits, biotic processes and the effects of dispersal‐limitation. Using commercially harvested blacklip (Haliotis rubra) and greenlip abalone (Haliotis laevigata) as case studies, we determine the relative importance of accounting for interactions among physiology, metapopulation dynamics and exploitation in predictions of range (geographical occupancy) and abundance (spatially explicit density) under various climate change scenarios. Traditional correlative ecological niche models (ENM) predict that climate change will benefit the commercial exploitation of abalone by promoting increased abundances without any reduction in range size. However, models that account simultaneously for demographic processes and physiological responses to climate‐related factors result in future (and present) estimates of area of occupancy (AOO) and abundance that differ from those generated by ENMs alone. Range expansion and population growth are unlikely for blacklip abalone because of important interactions between climate‐dependent mortality and metapopulation processes; in contrast, greenlip abalone should increase in abundance despite a contraction in AOO. The strongly non‐linear relationship between abalone population size and AOO has important ramifications for the use of ENM predictions that rely on metrics describing change in habitat area as proxies for extinction risk. These results show that predicting species' responses to climate change often require physiological information to understand climatic range determinants, and a metapopulation model that can make full use of this data to more realistically account for processes such as local extirpation, demographic rescue, source‐sink dynamics and dispersal‐limitation.     

Does host outcrossing disrupt compatibility with heritable symbionts?

Vertically transmitted microbes are common in macro‐organisms and can enhance host defense against environmental stress. Because vertical transmission couples host and symbiont lineages, symbionts may become specialized to host species or genotypes. Specialization and contrasting reproductive modes of symbiotic partners could create incompatibilities between inherited symbionts and novel host genotypes when hosts outcross or hybridize. Such incompatibilities could manifest as failed colonization or poor symbiont growth in host offspring that are genetically dissimilar from their maternal host. Moreover, outcrossing between host species could influence both host and symbiont reproductive performance. We tested these hypotheses by manipulating outcrossing between populations and species of two grasses, Elymus virginicus and E. canadensis, that host vertically transmitted fungal endophytes (genus Epichloё). In both greenhouse and field settings, we found that host–symbiont compatibility was robust to variation in host genetic background, spanning within‐population, between‐population and between‐species crosses. Symbiont transmission into the F₁ generation was generally high and weakly affected by host outcrossing. Furthermore, endophytes grew equally well in planta regardless of host genetic background and transmitted at high frequencies into the F₂ generation. However, outcrossing, especially inter‐specific hybridization, reduced reproductive fitness of the host, and thereby the symbiont. Our results challenge the hypothesis that host genetic recombination, which typically exceeds that of symbionts, is a disruptive force in heritable symbioses. Instead, symbionts may be sufficiently generalized to tolerate ecologically realistic variation in host outcrossing.     

Quantifying the climatic niche of symbiont partners in a lichen symbiosis indicates mutualist‐mediated niche expansions

The large distributional areas and ecological niches of many lichenized fungi may in part be due to the plasticity in interactions between the fungus (mycobiont) and its algal or cyanobacterial partners (photobionts). On the one hand, broad‐scale phylogenetic analyses show that partner compatibility in lichens is rather constrained and shaped by reciprocal selection pressures and codiversification independent of ecological drivers. On the other hand, sub‐species‐level associations among lichen symbionts appear to be environmentally structured rather than phylogenetically constrained. In particular, switching between photobiont ecotypes with distinct environmental preferences has been hypothesized as an adaptive strategy for lichen‐forming fungi to broaden their ecological niche. The extent and direction of photobiont‐mediated range expansions in lichens, however, have not been examined comprehensively at a broad geographic scale. Here we investigate the population genetic structure of Lasallia pustulata symbionts at sub‐species‐level resolution across the mycobiont's Europe‐wide range, using fungal MCM7 and algal ITS rDNA sequence markers. We show that variance in occurrence probabilities in the geographic distribution of genetic diversity in mycobiont‐photobiont interactions is closely related to changes in climatic niches. Quantification of niche extent and overlap based on species distribution modeling and construction of Hutchinsonian climatic hypervolumes revealed that combinations of fungal–algal interactions change at the sub‐species level along latitudinal temperature gradients and in Mediterranean climate zones. Our study provides evidence for symbiont‐mediated niche expansion in lichens. We discuss our results in the light of symbiont polymorphism and partner switching as potential mechanisms of environmental adaptation and niche evolution in mutualisms.     

Cohabitation and roommate bias of symbiotic bacteria in insect hosts

Symbiotic interactions between insects and bacteria have long fascinated ecologists. Aphids have emerged as the model system on which to study the effect of endosymbiotic bacteria on their hosts. Aphid‐symbiont interactions are ecologically interesting as aphids host multiple secondary symbionts that can provide broad benefits, such as protection against heat stress or specialist natural enemies (parasitic wasps and entomopathogenic fungi). There are nine common aphid secondary symbionts and individual aphids host on average 1–2 symbionts. A cost‐benefit trade‐off for hosting symbionts is thought to explain why not all aphids host every possible symbiont in a population. Both positive and negative associations between various symbionts occur, and this could happen due to increased costs when cohosting certain combinations or as a consequence of competitive interactions between the symbionts within a host. In this issue of Molecular Ecology, Mathé‐Hubert, Kaech, Hertaeg, Jaenike, and Vorburger (2019) use data on the symbiont status of field‐collected aphids to inform a model on the evolution of symbiont co‐occurrence. They vary the effective female population size as well as the rate of horizontal and maternal transmission to infer the relative impact of symbiont‐symbiont interactions versus random drift. Additional data analysis revisits an association between two symbionts in a fruit fly species using a long‐term data set to highlight that such interactions are not limited to aphids.     

Symbiont‐mediated phenotypic variation without co‐evolution in an insect–fungus association

Recent studies have shown that symbionts can be a source of adaptive phenotypic variation for their hosts. It is assumed that co‐evolution between hosts and symbionts underlies these ecologically significant phenotypic traits. We tested this assumption in the ectosymbiotic fungal associate of the gall midge Asteromyia carbonifera. Phylogenetic analysis placed the fungal symbiont within a monophyletic clade formed by Botryosphaeria dothidea, a typically free‐living (i.e. not associated with an insect host) plant pathogen. Symbiont isolates from four divergent midge lineages demonstrated none of the patterns common to heritable microbial symbioses, including parallel diversification with their hosts, substitution rate acceleration, or A+T nucleotide bias. Amplified fragment length polymorphism genotyping of the symbiont revealed that within‐lineage genetic diversity was not clustered along host population lines. Culture‐based experiments demonstrated that the symbiont‐mediated variation in gall phenotype is not borne out in the absence of the midge. This study shows that symbionts can be important players in phenotypic variation for their hosts, even in the absence of a co‐evolutionary association.     

A symbiont's dispersal strategy: condition-dependent dispersal underlies predictable variation in direct transmission among hosts

Direct horizontal transmission of pathogenic and mutualistic symbionts has profound consequences for host and symbiont fitness alike. While the importance of contact rates for transmission is widely recognized, the processes that underlie variation in transmission during contact are rarely considered. Here, we took a symbiont''s perspective of transmission as a form of dispersal and adopted the concept of condition-dependent dispersal strategies from the study of free-living organisms to understand and predict variation in transmission in the cleaning symbiosis between crayfish and ectosymbiotic branchiobdellidan worms. Field study showed that symbiont reproductive success was correlated with host size and competition among worms for microhabitats. Laboratory experiments demonstrated high variability in transmission among host contacts. Moreover, symbionts were more likely to disperse when host size and competition for microhabitat created a fitness environment below a discrete minimum threshold. A predictive model based on a condition-dependent symbiont dispersal strategy correctly predicted transmission in 95% of experimental host encounters and the exact magnitude of transmission in 67%, both significantly better than predictions that assumed a fixed transmission rate. Our work provides a dispersal-based understanding of symbiont transmission and suggests adaptive symbiont dispersal strategies can explain variation in transmission dynamics and complex patterns of host infection.     

Vertical transmission as the key to the colonization of Madagascar by fungus-growing termites?

The mutualism between fungus-growing termites (Macrotermitinae) and their mutualistic fungi (Termitomyces) began in Africa. The fungus-growing termites have secondarily colonized Madagascar and only a subset of the genera found in Africa is found on this isolated island. Successful long-distance colonization may have been severely constrained by the obligate interaction of the termites with fungal symbionts and the need to acquire these symbionts secondarily from the environment for most species (horizontal symbiont transmission). Consistent with this hypothesis, we show that all extant species of fungus-growing termites of Madagascar are the result of a single colonization event of termites belonging to one of the only two groups with vertical symbiont transmission, and we date this event at approximately 13 Mya (Middle/Upper Miocene). Vertical symbiont transmission may therefore have facilitated long-distance dispersal since both partners disperse together. In contrast to their termite hosts, the fungal symbionts have colonized Madagascar multiple times, suggesting that the presence of fungus-growing termites may have facilitated secondary colonizations of the symbiont. Our findings indicate that the absence of the right symbionts in a new environment can prevent long-distance dispersal of symbioses relying on horizontal symbiont acquisition.     

Coevolutionary dynamics in one-to-many mutualistic systems

“One-to-many” mutualisms are often observed in nature. In this type of mutualism, each host individual can interact with many symbionts, whereas each individual symbiont can interact with only one host individual. Partner choice by the host is a potentially critical mechanism for maintaining such systems; however, its long-term effects on the coevolution between the hosts and symbionts have not been completely explored. In this study, I developed a simple mathematical model to describe the coevolutionary dynamics between hosts and symbionts in a one-to-many mutualism. I assumed that each host chooses a constant number of symbionts from a potential symbiont population, a fraction of which are chosen through preferential choice on the basis of the cooperativeness of the symbionts and the rest are chosen randomly. Using numerical calculations, I found that mutualism is maintained when the preferential choice is not very costly and the mutation rate of symbionts is large. I also found that symbionts that receive benefits from hosts without a return (cheater symbionts) and hosts that do not engage in preferential partner choice (indiscriminator hosts) can coexist with mutualist symbionts and discriminator hosts, respectively. The parameter domain of pure mutualism, i.e., free from cheater symbionts and indiscriminator hosts, can be narrower than the whole domain where the mutualism persists.     

Lucinidae/sulfur-oxidizing bacteria: ancestral heritage or opportunistic association? Further insights from the Bohol Sea (the Philippines)

The first studies of the 16S rRNA gene diversity of the bacterial symbionts found in lucinid clams did not clarify how symbiotic associations had evolved in this group. Indeed, although species-specific associations deriving from a putative ancestral symbiotic association have been described (coevolution scenario), associations between the same bacterial species and various host species (opportunistic scenario) have also been described. Here, we carried out a comparative molecular analysis of hosts, based on 18S and 28S rRNA gene sequences, and of symbionts, based on 16S rRNA gene sequences, to determine as to which evolutionary scenario led to modern lucinid/symbiont associations. For all sequences analyzed, we found only three bacterial symbiont species, two of which are harbored by lucinids colonizing mangrove swamps. The last symbiont is the most common and was found to be independent of biotope or depth. Another interesting feature is the similarity of ctenidial organization of lucinids from the Philippines to those described previously, with the exception that two bacterial morphotypes were observed in two different species (Gloverina rectangularis and Myrtea flabelliformis). Thus, there is apparently no specific association between Lucinidae and their symbionts, the association taking place according to which bacterial species is present in the environment.     

Environmental cues and symbiont microbe‐associated molecular patterns function in concert to drive the daily remodelling of the crypt‐cell brush border of the Euprymna scolopes light organ

Recent research has shown that the microbiota affects the biology of associated host epithelial tissues, including their circadian rhythms, although few data are available on how such influences shape the microarchitecture of the brush border. The squid‐vibrio system exhibits two modifications of the brush border that supports the symbionts: effacement and repolarization. Together these occur on a daily rhythm in adult animals, at the dawn expulsion of symbionts into the environment, and symbiont colonization of the juvenile host induces an increase in microvillar density. Here we sought to define how these processes are related and the roles of both symbiont colonization and environmental cues. Ultrastructural analyses showed that the juvenile‐organ brush borders also efface concomitantly with daily dawn‐cued expulsion of symbionts. Manipulation of the environmental light cue and juvenile symbiotic state demonstrated that this behaviour requires the light cue, but not colonization. In contrast, symbionts were required for the observed increase in microvillar density that accompanies post dawn brush‐border repolarization; this increase was induced solely by host exposure to phosphorylated lipid A of symbiont cells. These data demonstrate that a partnering of environmental and symbiont cues shapes the brush border and that microbe‐associated molecular patterns play a role in the regulation of brush‐border microarchitecture.     

A shift to parasitism in the jellyfish symbiont Symbiodinium microadriaticum

One of the outstanding and poorly understood examples of cooperation between species is found in corals, hydras and jellyfish that form symbioses with algae. These mutualistic algae are mostly acquired infectiously from the seawater and, according to models of virulence evolution, should be selected to parasitize their hosts. We altered algal transmission between jellyfish hosts in the laboratory to examine the potential for virulence evolution in this widespread symbiosis. In one experimental treatment, vertical transmission of algae (parent to offspring) selected for symbiont cooperation, because symbiont fitness was tied to host reproduction. In the other treatment, horizontal transmission (infectious spread) decoupled symbiont fitness from the host, potentially allowing parasitic symbionts to spread. Fitness estimates revealed a striking shift to parasitism in the horizontal treatment. The horizontally transmitted algae proliferated faster within hosts and had higher dispersal rates from hosts compared to the vertical treatment, while reducing host reproduction and growth. However, a trade-off was detected between harm caused to hosts and symbiont fitness. Virulence trade-offs have been modelled for pathogens and may be critical in stabilising 'infectious' symbioses. Our results demonstrate the dynamic nature of this symbiosis and illustrate the potential ease with which beneficial symbionts can evolve into parasites.     

Examination of the success rate of secondary symbiont manipulation by microinjection methods in the pea aphid system

Insects harbor a wide range of microbial symbionts, but their influence on host phenotypes is described in a limited number of biological models. One experimental approach to gain knowledge on the effects of symbionts to their hosts is to create insect lines with and without symbionts and examine their phenotypes. However, the success rate of symbiont elimination and introduction methods is dependent on several parameters that are scarcely tested or described. The pea aphid, Acyrthosiphon pisum Harris (Hemiptera: Aphididae), is a model insect of symbiosis studies. It harbors a primary symbiont that supplies the host with essential amino acids, and an array of secondary symbionts whose effects have been assessed by manipulating their presence/absence in the insect. Here, we describe the influence of key parameters on the success rate of symbiont manipulation using the pea aphid–secondary symbiont system. We compared two elimination methods differing in antibiotic treatment using several aphid–symbiont combinations. We also created new aphid host–symbiont combinations by secondary symbiont introduction and examined the effects of larval stage of recipient aphids on introduction success. Our study revealed that the aphid–symbiont combination has strong influence on both symbiont introduction and elimination success rates, and that the type of antibiotics and the larval stage of recipient aphids influence the elimination and introduction success rate, respectively.