BACKGROUND
every plant and animal on this planet is engaged in beneficial and antagonistic relationships with microbes that profoundly influence their physiology, life history, and evolution.
the ubiquity and biological importance of long-term, stable, and beneficial host-microbe associations, or mutualisms, have been the focus of many investigations .
models of mutualisms possess common factors that contribute to their success, maintenance, and perpetuation.
cooperative association, fitness benefits for all parties, partner fidelity and preference are key factors that influence the establishment and maintenance of mutualisms  <cit> .
in this respect, the relationship between gram-negative xenorhabdus bacteria and their steinernema nematode hosts yields insights into processes that influence symbiont specificity, selection, and function.
furthermore, the inherent experimental tractability of these organisms and availability of genome sequences has spurred developments across a broad scientific front, such that this mutualism is now viewed as a model for basic research in ecology, evolution, biochemistry, and molecular genetics of symbiosis .
briefly, the xenorhabdus life cycle depends on two eukaryote hosts, steinernema nematodes and an insect  for their survival, dissemination, and reproduction.
the relationship between the bacteria and the nematodes is benign and mutually beneficial, whereas their relationship with the insect host is pathogenic  <cit> .
after invading an insect, nematodes  defecate xenorhabdus symbionts into the hemocoel where they multiply rapidly, neutralize the insect’s immune system, induce septicemia and toxemia, and begin digesting the fat-rich insect cadaver.
within the compromised insect, nematodes mature, and sexually reproduce resulting in many thousands of offspring using the metabolic resources liberated from the cadaver  <cit> .
once resources are exhausted, the nematodes develop into a non-feeding, environmentally-resistant “dauer” stage, known as the infective juvenile .
prior to emerging from the insect cadaver, ijs re-associate with xenorhabdus symbionts, which colonize the ij nematodes’ intestinal receptacle.
ijs carrying bacterial symbionts persist in the soil until a new insect host can be infected.
nutrients provisioned by the nematode ij are thought to support outgrowth and bacterial survival until the next insect host  <cit> .
each steinernema species associates with a single xenorhabdus species.
however, one xenorhabdus species can form associations with more than one steinernema host  <cit> .
the most promiscuous xenorhabdus species is xenorhabdus bovienii, which associates with numerous steinernema species comprising two different evolutionary clades in the steinernema phylogeny  <cit> .
previous studies revealed that the bacteria can contribute to virulence toward the insect host and support nematode reproduction  <cit> .
these benefits have also been demonstrated to have a strong phylogenetic correlation at the nematode intraspecific level  <cit>  with the driving selective interest of nematode symbiont choice likely differing between nematode species  <cit> .
complicating matters, in nature, insects can become co-infected with two steinernema species each harboring their respective bacterial symbiont  <cit> .
this requires that nematodes disassociate with their bacterial partner, reproduce and mature within a mixed community, and then re-associate with the correct symbiont to avoid the fitness costs associated with harboring the wrong bacterium.
this sets entomopathogenic nematode-bacterium symbioses apart from many other mutualistic systems where partners do not naturally disassociate under optimal environmental conditions.
symbiont switching can have detrimental effects to animal hosts.
in examples across invertebrate-microbe symbioses, symbionts serve important functions and provide adaptive features for their hosts by aiding in digestion, conferring pathogen/parasite resistance, and increasing host fitness  <cit> .
in fact, within a single species of nematode, association with a closely-related non-cognate symbiont results in reduced fitness  <cit> .
for this reason, partner fidelity and localization are especially important to host animals.
steinernema-xenorhabdus symbioses are excellent models for testing hypotheses concerning mutualism.
obligate mutualisms evolve and persist due to increased fitness of both parties through maintenance of partner fidelity and partner preference.
unlike many intracellular, vertically transmitted animal-microbe symbioses, entomopathogenic nematode-bacterial symbioses rely on partner re-association after each infection cycle.
due to the tractability of this system, we were able to test the fitness, fidelity, and preference of several nematode-bacterial partner pairs by switching their symbionts.
we show that across a number of x. bovienii-steinernema spp.
combinations, cognate partners display the hallmarks of mutualism including increased fitness in cooperation, partner fidelity, and partner preference.
additionally, we show that the most derived steinernema species tested s. puntauvense, is more promiscuous and is able to form relatively successful relationships with a wide variety of non-cognate xenorhabdus species.
contrastingly, the ancestral x. bovienii host, s. intermedium, is the most fastidious of the nematode species examined, performing poorly with all non-cognate bacteria.
xenorhabdus bovienii  is attractive to the largest variety of non-cognate nematode species.
taken together, these data show entomopathogenic nematode-bacterial mutualisms represent an important and useful model for understanding the evolution of mutualisms.
methods
nematode and bacteria isolates
three steinernema species , which are hosts of xenorhabdus bovienii, were used in this study to assess the cost/benefit of symbiont switching .
two of these species  were assessed in symbiont choice assays to determine symbiont preference .
additionally, steinernema feltiae and steinernema carpocapsae were included  in these assays given that these taxa have been extensively studied for the effects of symbiont switching events  <cit> .
details on nematode hosts and bacterial symbionts species/strain names, geographic origin, and sequence accession numbers are listed in table  <dig>  identity of nematodes and bacteria considered in this study was molecularly confirmed prior to the initiation of the experiments by sequencing of 28s and 12s ribosomal dna  genes for the nematode hosts and 16s genes for their bacterial symbionts.table  <dig>
steinernema and xenorhabdus species and strains considered in this study
s. carpocapsae
x. nematophila
s. feltiae
x. bovienii
s. feltiae
x. bovienii
s. feltiae
x. bovienii
s. intermedium
x. bovienii
s. kraussei
x. bovienii
s. kraussei
x. bovienii
s. jollieti
x. bovienii
s. oregonense
x. bovienii
s. puntauvense
x. bovienii
nematodes were reared in vivo with last instar wax moth larvae, galleria mellonella   following procedures described by kaya and stock  <cit> .
modified white traps were used to harvest emerging infective juvenile stages , which were stored in tissue culture flasks at a concentration of approximately  <dig> ijs/ml in 250 ml tissue culture flasks .
bacterial symbionts were extracted from each nematode host species/strain by sonication of naturally colonized ijs.
bacterial stocks were established by growing each bacterial species/strain to mid-log phase in luria-bertani  broth supplemented with  <dig> %  sodium pyruvate  <cit>  grown in a shaking, dark incubator at 28 °c following procedures described by akhurst  <cit> .
bacterial stocks were stored at −80 °c in lb broth supplemented with 20%  glycerol.
rearing of aposymbiotic nematodes
aposymbiotic  nematodes were obtained via alkaline axenization of eggs harvested from first generation gravid females as described by mcmullen and stock  <cit> .
eggs were seeded on liver-kidney agar  <cit> .
aposymbiotic ij nematodes were collected by removing the bottom portion with agar of the petri dish and placing it in a sterile modified white trap  <cit> .
ijs were stored in tissue culture flasks at 15 °c until use for experiments following the same procedures described above.
the aposymbiotic nature of the ijs was confirmed by grinding them during the first and second week of emergence with a motorized pestle  in lb broth.
the suspension was then plated onto differential nutrient agar supplemented with  <dig> %  bromothymol blue,  <dig> %  triphenyl-tetrazolium, and  <dig> %  sodium pyruvate   <cit>  to confirm presence/absence of bacterial symbionts.
host-symbiont interactions assays
host virulence
aposymbiotic ijs were surface sterilized according to procedures described by stock and goodrich-blair  <cit> .
ij suspension as diluted in grace’s medium to a concentration of  <dig> ijs in 10 μl.
bacteria cultures were grown overnight and then sub-cultured to a 1: <dig> dilution into lb broth and subsequently grown to od600nm  ~ <dig> .
bacterial suspension was then diluted in grace’s insect medium  to a concentration 200 cfu  in 10 Μl.
prior to injections, bacterial suspensions were plated to confirm that equal counts were used for all bacterial strains.
the experimental arena was a 12-well plate.
each g. mellonella larvae was co-injected with a nematode/bacterium treatment.
treatments consisted of co-injections of different steinernema/xenorhabdus species/strain pairs .
bacteria and nematodes were co-injected at the concentrations described above.
twelve g. mellonella larvae were used per treatment and each treatment were repeated three times independently.
insect mortality was monitored daily for 7 days post-injection.
g. mellonella larvae were considered dead when they did not respond to the gentle prodding of soft-grip forceps.
resulting cadavers and progeny were used for all subsequent experiments, except for the choice assays.table  <dig> cognate and non-cognate nematode-symbiont pairs considered in this study
grey boxes denote cognate symbiont-nematode pair
host reproductive fitness
reproductive fitness of s. intermedium, s. oregonense and s. puntauvense was measured considering two parameters: a) infection productivity , and b) emerging ij progeny .
for productive infection percentages, g. mellonella cadavers obtained from the virulence assays were considered.
each cadaver was thoroughly rinsed and transferred into a modified white trap at day seven post-injection.
cadavers were checked daily to record progeny emergence for up to 30 days.
cadavers producing at least  <dig> ijs were scored as ‘productive’, and productive infection percentage was calculated as the number of cadavers producing progeny out of the total number.
the emerging ij progeny was determined by collecting the nematode suspension from each modified white trap daily for 10 days of post-emergence and following volumetric dilution procedures per kaya and stock  <cit> .
symbiont carriage
ij symbiont carriage was determined following procedures described by goestch et al.
<cit> .
two technical replicates and three biological replicates were performed for each nematode-bacterium combination and averaged to obtain cfu ij− <dig>
f <dig> virulence
virulence of the emerging ij progeny  was tested considering a one-on-one assay as described by kaya and stock  <cit> .
the experimental arena was a 12-well plate, where each well was filled with 5 g of sterile sand.
an inoculum of  <dig> ijs delivered in 100 μl of distilled water was poured onto the sand.
one g. mellonella larvae was added into each well.
the lid was placed and the plate was incubated at 20 °c in the dark.
mortality was monitored every 8 h for a total of 96 h. experiments were repeated three times.
for data analysis, percent mortality registered at 48 h was used, as it showed the largest difference between combinations.
longevity of f <dig> progeny
ijs obtained from the host reproductive fitness assays were examined 6-months post-emergence at 30× magnification.
nematodes were considered dead when they adopted a straight shape and/or showed no movement, whereas nematodes were considered live when they were moving or had a j-shape position.
ijs were confirmed as dead by gentle prodding.
two technical replicates of at least  <dig> ijs were counted from each productive infection.
morphometric analysis
the effect of cognate/non-cognate bacterial symbiont on nematode size was assessed by measuring length and width of ij progeny  obtained from each treatment.
a total of  <dig> ijs/ treatment were randomly pooled.
ijs were heat killed in m <dig> buffer at 50-60 °c  <cit>  and subsequently fixed with triethanolamine formalin  following procedures described by goodey  <cit> .
ijs were mounted onto glass slides and measured using an olympus bx <dig> microscope equipped with differential interference contrast optics and olympus microsuite software .
bacteria choice assays
choice assays were performed to test whether different steinernema hosts of x. bovienii could discern between bacterial strains that are cognate vs. non-cognate.
for this purpose, three steinernema hosts of x. bovienii were considered: s. feltiae , s. oregonense , and s. puntauvense .
in addition, two other x. bovienii strains were included, xbffl  and xbi  to account for symbionts that have a range of fitness effects on nematode hosts.
s. carpocapsae  and its respective symbiont x. nematophila were considered for outgroup comparisons.
for each nematode species used, a no-choice assay was performed where all bacterial foci in our test arena  consisted of the symbiotic partner as a control.
aposymbiotic ijs were reared following procedures described by mcmullen and stock  <cit> .
prior to the experimental setup, nematodes were surface sterilized using the vacuum filtration method described by kaya and stock  <cit> .
due to low fitness of aposymbiotic s. intermedium nematodes, bacterial preference assays did not include this species as this would likely confound results.
bacterial strains were grown from glycerol stocks overnight  in 5 ml of lb broth with  <dig> % sodium pyruvate in a 15 ml centrifuge tube at 28 ± 1 °c in the dark with agitation.
50 μl of bacterial suspension was sub-cultured following the same conditions, but allowed to grow for 9 h prior to assay setup.
cultures were streaked onto nbta to confirm liquid cultures were phase i under the same conditions described above.
briefly, a 10 cm petri dish was filled with 20 ml of lipid agar  <cit>  supplemented with  <dig> % sodium pyruvate  <cit> .
a bacterium inoculum of 10 μl was inoculated onto the agar plate in four equidistant positions, approximately  <dig>  cm away from the edge.
for each plate, one cognate and one non-cognate bacterial strain were applied so that two inoculum positions of the same strain were not next to each other .
a 50 μl inoculum of  <dig> ijs was dispensed into the center of the plate.
there were five replicates for each nematode-bacteria combination.
controls consisted of four foci of cognate bacteria.
plates were incubated at 20-25 °c for 3 days in the dark.
this time period allowed for enough time for ijs to migrate to the preferred bacteria.
experiments were repeated twice, independently.
after 3 days, nematode migration to bacteria foci was evaluated.
to facilitate visualization and counting of nematodes, bacteria foci and nematode inoculum area  were cut and placed into a 5 cm petri dish containing approximately 10 ml of m <dig> buffer.
dishes were then placed on a shaker overnight to expedite release of the nematodes embedded in the agar and bacterial lawn.
the entire plate was then rinsed off to collect the ijs and stored in 50 ml centrifuge tubes at 4 °c until counted.
an attraction activity index  was calculated to assess nematode attraction for bacteria in the choice assays.
for each treatment, the number of nematodes found in each bacterial focus  was combined.
similar to zhang et al.
<cit> , the formula below was used to obtain the aai values.
aai ranges from  <dig> to − <dig> and the resulting values can be inferred as: 1 = all nematodes are attracted to non-cognate symbiont, 0 = equally prefer both bacteria, −1 = all nematode are repulsed by non-cognate strain.
the formula used was: aai =  –  / total number of nematodes in bacteria.
for no-choice assay controls, bacterial foci pairs were randomly assigned to follow data analysis procedure used for choice assays.
phylogenetic analyses
phylogenetic analyses were conducted to determine relative phylogenetic distance between cognate and non-cognate xenorhabdus strains.
mega v. <dig> was used to generate both nematode and bacterial phylogenies  <cit> .
for the xenorhabdus phylogeny, the following sequences were aligned using muscle  <cit>  and concatenated: 16s rdna, reca, dnan, gltx, gyrb, and infb .
maximum likelihood was carried out using a tn93 + g + i model of evolution with  <dig> bootstrap replicates.
the nematode phylogeny was generated using similar methods as the bacterial phylogeny with the mitochondrial and nuclear rrna genes: 12s and 28s, respectively .
however, a gtr + g model of evolution was used.
s. carpocapsae and its native symbiont x. nematophila were considered in nematode and bacterial phylogenies, respectively, as out group taxa.
statistical analyses
all statistical analyses were carried out using r, version  <dig> .
<dig>  <cit> .
nematode-bacterium virulence against g. mellonella were analyzed using a cox mixed effects model  from the package ‘coxme’  <cit> .
each nematode host tested was analyzed separately and their cognate symbiont was used as the reference strain to build the cox model.
the bacterial strain was considered as a fixed effect in the model, while a random effect was included to better explain the variance due to differences between experimental setups.
a post hoc tukey hsd test was conducted using the ‘multcomp’ package  <cit> .
final insect mortality at day seven and median mortality time  were calculated from data assessed with the ‘survival’ package  <cit> .
if a nematode-bacterium combination resulted in no deaths during the study period, a single death was added at day seven to allow the model to be stable.
either a one-way analysis of variance  or mixed effects anova  was considered to analyze most assays.
a one-way anova was considered to analyze productive infection and f <dig> progeny data.
in both instances, a single percentage was generated for each nematode-bacterium pair from each experimental setup.
the remaining measures  were assessed with a mixed effect model.
for all analyses, each nematode host was analyzed separately and the bacterial treatment was used as a fixed effect.
when a mixed effect anova was used, the experimental setup was used as a random effect to account for setup variation with the packages ‘lme4’  <cit>  and ‘car’  <cit> .
all residuals were visually inspected for normality and homoscedasticity.
a post hoc tukey hsd test was performed for pairwise comparisons between bacterial treatments with the package ‘lsmeans’  <cit> .
nematode choice data was analyzed by performing a two-sided one-sample t-test for each aai index with a μ =  <dig> .
p-values were adjusted using the benjamini-hochberg method to correct for false discovery rate in each nematode dataset.
linear regression analyses  were conducted on nematode and bacterial fitness and morphometric data.
nematode fitness was calculated as a factor of progeny count, progeny survival, and productive infection percentage relative to the cognate symbiont mean, while bacterial fitness was determined as a factor of nematode fitness and cfu ij− <dig> relative to the average cognate symbiont mean.
for morphometric analyses, the factor of length and width was analyzed.
each nematode host was analyzed separately.
RESULTS
all parameters measured across nematode-symbiont combinations were influenced by the pairing at hand with the exception of nematode size .
together the measured variables demonstrate the truly mutualistic nature of associations between steinernema nematodes and their xenorhabdus symbionts.
some non-cognate symbionts are sufficient to maintain virulence against insect hosts.
all steinernema spp.
tested demonstrated greater virulence with their cognate symbionts than with any other, non-cognate bacterium.
of them, symbiont switching in s. intermedium, a clade i species, had the most pronounced impact on virulence.
association with the cognate symbiont, xbi, with s. intermedium resulted in 50% insect mortality in 4 days and ultimately killed 100% of insect hosts within a week .
s. intermedium was also capable of moderate levels of virulence with both xbkbu and xbj resulting in  <dig> % and  <dig> % mortality by day seven, respectively .
association with xbffl and xbp did result in some insect mortality after 1 week, although lt <dig> was not observed within the 7-day observation period, and virulence when s. intermedium was colonized by these species was not statistically distinguishable from aposymbiotic s. intermedium or those species that were completely incompetent .table  <dig> non-cognate x. bovienii effect on nematode virulence
x. bovienii
s. intermedium a
s. oregonense a
s. puntauvense a
% mortalityb
<dig>  ± 0
4
a
<dig>  ±  <dig>
2
a
<dig>  ±  <dig>
2
a
acognate symbiont data are bolded
bpercent mortality after 1 week with standard error
clethal time to kill 50% of population reported as days
dletter-ranking generated from a cox mixed effects model and post hoc tukey test; see additional file 1: table s <dig> for summary of model
contrastingly, both assayed clade iii nematode species, s. puntauvense and s. oregonense, maintained some level of virulence with all tested bacteria and even without a bacterial partner .
as with s. intermedium, s. oregonense was most virulent when associated with their cognate symbiont; lt <dig> was shorter and 7-day insect mortality is higher in associations with xbo .
s. oregonense and s. puntauvense exhibited moderate virulence with all x. bovienii tested and with no associated symbiont .
interestingly, s. puntauvense had comparable lt <dig> times with many symbionts, though its cognate still produced the highest insect mortality at 7-days post-inoculation .
symbiont switching has negative impacts on fitness
despite many steinernema-xenorhabdus combinations producing successful infections in insect hosts, associations with non-cognate symbionts impacted nematode reproductive success.
the first measure of this impact was reflected in the proportion of productive infections from a given nematode-bacterium association.
we determined an infection to be productive if it resulted in the emergence of at least  <dig> ijs.
all species suffered negligible infection productivity when reared aposymbiotically and non-cognate associations were not sufficient to recover that loss in s. intermedium and s. oregonense .
s. intermedium had the least productive infections with non-cognate bacteria, with only xbkbu partially restoring infection productivity, and all combinations were not significantly different from the uncolonized aposymbiotic worms .
s. intermedium produced small numbers of ijs when partnered with xbj and xbkbu .
associations between non-cognates xbkca and xbp generated a comparable proportion of productive infections to the cognate symbiont in s. oregonense; however, the total number of progeny from these infections was significantly less than what was measured from xbo associations .
in contrast, no progeny was recovered when s. puntauvense associated with non-cognate bacterium, even though most combinations resulted in productive infections .
however, this pattern was not seen in combinations with xbi, xbj, xbkbu, nor xbo .fig.
<dig> productive infections resulting from steinernema-xenorhabdus experimental associations.
bars indicate the mean and se resulting from a one-way anova.
letters above bars show statistical ranking from post hoc tukey test.
n =  <dig> for each nematode-bacterium treatment.
these values were calculated as the percent of cadavers that produce progeny in each independent experimental setup.
cognate symbionts are highlighted with a heavy outline.
a. s. intermedium: f <dig>  20 =  <dig> , p <  <dig> , r
2 =  <dig> ; b. s. oregonense: f <dig>  20 =  <dig> , p <  <dig> , r
2 =  <dig> ; and c. s. puntauvense: f <dig>  20 =  <dig> , p  <dig> , r
2 =  <dig>
fig.
<dig>
steinernema progeny production with novel symbionts.
bars show the mean and se resulting from a mixed effects anova.
letters above bars demonstrate results from post hoc tukey test.
n =  <dig> for each nematode-bacterium treatment.
cognate symbionts are highlighted with a heavy outline.
a. s. intermedium: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig> ; b. s. oregonense: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig> ; and c. s. puntauvense: f <dig>  313 =  <dig> , p <  <dig> .
r
2
o =  <dig>
we also tested the longevity and virulence of the f <dig> ijs that emerged following association with non-cognate bacteria.
in general, the patterns observed in ij longevity followed the trends seen in progeny production.
nematode-symbiont pairings that yielded higher progeny production, also tended to have increased longevity.
six-months post-emergence, s. intermedium-xbj and -xbkbu pairings had ~20% ij survival, whereas s. oregonense-xbkca and -xbp pairings exhibited ~25% ij survival .
while these are both significantly less than the cognate symbiont, these pairings had higher survival than the other non-cognates that were capable of producing progeny .
all symbiont pairings with s. puntauvense, except xbi, had viable ijs 6-months post-emergence , though none of the non-cognate pairings survived as well as those associated with the cognate symbiont .fig.
<dig>
steinernema progeny survival with non-cognate symbionts.
bars indicate the mean and se resulting from a mixed effects anova.
letters above bars show results from post hoc tukey test.
n =  <dig> for each nematode-bacterium treatment.
data represents percentage of nematode alive in a pool of  <dig> ijs.
cognate symbionts are highlighted with a heavy outline.
a. s. intermedium: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig> ; b. s. oregonense: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig> ; and c. s. puntauvense: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig>
in addition to f <dig> ij longevity, virulence of this new generation is imperative for nematode success long term and establishment of novel associations.
strikingly, f <dig> ijs from all non-cognate associations had significantly less virulence than cognate associations .
in combinations where progeny were produced some virulence was observed no matter how slight .
consistent with other observations, s. intermedium was most susceptible to symbiont switching and resulted in a maximum of ~30% virulence in f <dig> ijs; whereas f1 s. puntauvense ijs had the most success with non-cognate symbionts and resulted in several combinations causing insect mortality >40% .
s. oregonense had intermediate success out of the three nematodes tested, with an observed insect mortality below 40% .fig.
<dig> f <dig> progeny virulence against g. mellonella.
bars indicate the mean and se resulting from a one-way anova.
letters above bars show results from post hoc tukey test.
n =  <dig> for each nematode-bacterium treatment.
values were calculated as the percent of insects that died by 48 h in each independent experimental setup.
cognate symbionts are highlighted with a heavy outline.
a. s. intermedium: f <dig>  18 =  <dig> , p =  <dig> , r
2 =  <dig> ; b. s. oregonense: f <dig>  18 =  <dig> , p <  <dig> , r
2 =  <dig> ; and c. s. puntauvense: f <dig>  18 =  <dig> , p  <dig> , r
2 =  <dig> .
post-
non-cognate bacteria are less successful at colonizing ijs
we also assayed the extent to which bacteria were able to establish themselves in cognate and non-cognate ijs by measuring bacterial carriage.
by and large, bacterial carriage follows the same patterns we saw in the other metrics discussed above.
s. intermedium carried the fewest non-cognate bacteria , s. oregonense can carry more species of non-cognate bacteria but at significantly fewer numbers , and s. puntauvense carried more non-cognate bacteria than the other two species though at lower concentrations than their cognate xbp .fig.
<dig>
steinernema symbiont carriage in f <dig> populations.
bars indicate the mean and se resulting from a mixed effects anova.
letters above bars show results from post hoc tukey test.
n =  <dig> for each nematode-bacterium treatment.
cognate symbionts are highlighted with a heavy outline.
a. s. intermedium: f <dig>  311 =  <dig> , p <  <dig> , r
2
o =  <dig> ; b. s. oregonense: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig> ; and c. s. puntauvense: f <dig>  313 =  <dig> , p <  <dig> , r
2
o =  <dig>
x. bovienii hosts discern between some non-cognates in bacterial choice assays
using choice assays, we determined that steinernema ijs were able to discern between some conspecific non-cognate bacterial strains and are repelled by heterospecific bacteria .
for example, s. carpocapsae was not attracted to any of the x. bovienii strains , with the exception of xbo .
the nematodes appeared to not be able to distinguish between this bacterium  and their cognate partner .
similarly, natural x. bovienii hosts  were not attracted to xn  .
steinernema species that naturally associate with x. bovienii exhibited strain-specific behavioral preferences within x. bovienii .
for instance, s. feltiae preferred xbffl and xbp equally as it did for its cognate symbiont.
this nematode was also attracted to xbi but repelled xbo .
s. oregonense was not attracted to xbffl and xbfsn, and did not demonstrate a differentiated preference for xbp and xbi and its cognate .
s. puntauvense did not display a preference for any of the tested conspecific bacteria .
interestingly, xbi was the only non-cognate bacterium that was always preferred  by the three x. bovienii nematode host species tested in this study .fig.
<dig> nematode preference for non-cognate symbionts.
a. s. feltiae, b. s. oregonense, c. s. puntauvense, and d. s. carpocapsae.
bars represent the mean and se resulting from two-sided one-sample t-test analyses with μ =  <dig>  p-values were adjusted with benjamini-hochberg method to correct for the false discovery rate and their results are indicated by the following symbols: ns = non-significant, * = p <  <dig> , ** = p <  <dig> , *** = p <  <dig> , **** = p <  <dig> .
n =  <dig> for each aai  value.
white bars indicate nematode attraction to non-cognate symbiont, while gray bars show repulsion.
cognate symbiont bars are highlighted with a heavy outline.
see additional file 1: table s <dig> for summary statistics
phylogeny as a predictor of fitness effects of non-cognate associations
a co-phylogenetic hypothesis for all bacteria-nematode pairs considered in this study was developed .
the overall fitness effects of all nematode-bacterial pairings were regressed against the phylogenetic distance between the test bacterium and the cognate bacterium .
in all cases, we show that fitness of both bacteria and nematodes declines as the phylogenetic distance between the cognate and non-cognate bacteria increases .
these results suggest that the more distant a bacterial strain is in relation to a nematode’s cognate symbiont, the more detrimental the fitness effects are on non-cognate nematode partner.
in particular, s. intermedium and s. oregonense are more sensitive to association with a phylogenetically distant symbiont, with  <dig> - <dig>  unit decrease in nematode and bacterial fitness for every  <dig>  nucleotide change in symbiont .
s. puntauvense is more tolerant to association with non-cognates and only suffers a  <dig> - <dig>  unit decrease in fitness for every  <dig>  bacterial nucleotide change .fig.
<dig> evolutionary relatedness between steinernema nematodes and their native symbionts and the fitness consequence of host-symbiont switching events.
a. co-phylogeny tanglegram of steinernema-xenorhabdus pairs used in experiments.
the steinernema phylogeny was generated with concatenated 12s and 28s rrna genes , while the xenorhabdus phylogeny shows the results from a concatenation of the following genes: 16s rrna, reca, dnan, gltx, gyrb, and infb .
internal nodes indicate bootstrap replicate values >50%.
panels b and c show relative nematode and bacterial fitness, respectively, by bacterial phylogenetic distance from native steinernema-xenorhabdus pairing.
for both graphs, data points indicate mean and se from a linear regression analysis; see additional file 1: table s <dig> for model summary.
s. intermedium = black triangles with black line, s. oregonense = gray squares with gray line, and s. puntauvense = hashed diamonds with dashed line
discussion
entomopathogenic nematodes and their bacterial partners are an important model system for studying mutualism.
the data presented here demonstrate that steinernema-xenorhabdus symbioses are mutually beneficial with regard to nematode and bacterial fitness across multiple species of nematode hosts.
association with cognate symbionts benefits the nematode host in terms of virulence against insect hosts, total progeny production, f <dig> progeny longevity, and virulence.
the bacterial partners also benefit by increased host colonization, and thus have an increased chance of being vectored into a new insect host, by their nematode host.
the final piece of this puzzle is partner fidelity.
not only did nematodes discern among different bacterial strains to varying degrees, but they tended to prefer bacteria which ultimately increase their virulence, longevity, and reproductive success.
generally, all cognate nematode-bacteria pairs outperformed all other combinations in terms of virulence, progeny production and longevity, bacterial carriage, and fitness of both parties.
with the exception of nematode size , all measured variables were influenced by symbiont partner associated with the given nematode host.
nematodes were also capable of discerning their cognate bacterium in a choice assay and tended to prefer their native symbiont to heterospecific strains, while their attraction toward non-cognate symbionts mirrored fitness and virulence trends measured in this study .
most importantly, the impacts for both nematode and bacterial fitness were correlated to the phylogenetic distance between the cognate bacterium and that colonizing partner .
these data solidify the argument that entomopathogenic nematodes and their bacterial symbionts represent a true mutualistic symbiosis.
while this trend is consistent across species, each tested nematode species had a distinct pattern of effects when associated with various x. bovienii strains.
for example, s. intermedium, the only species member of clade i tested, was very fastidious.
this nematode showed the highest virulence, progeny production, and f <dig> longevity with its native symbiont, xbi .
however, these same metrics were only minimally restored when the nematodes associated with xbkbu and xbj; no non-cognate symbiont restored even 50% virulence, ij production, or progeny success as the native xbi .
unfortunately, due to technical limitations, we were not able to determine bacterial preference in s. intermedium, but we did assay how other steinernema species responded to xbi.
interestingly, xbi was attractive to non-cognate nematode species such as s. feltiae, s. oregonense or s. puntauvense, despite conferring almost no fitness or virulence benefits .
given that its host, s. intermedium is in a basal  phylogenetic clade to these other x. bovienii hosts tested, xbi may exhibit some attractive trait which is generic to all x. bovienii, and as such is detectable by all x. bovienii host nematodes.
in this respect, it could be speculated an attraction trait, present in xbi, which more derived nematode hosts are cuing in on, may also be present in the cognate x. bovienii lineages.
it has been shown that differences in primary and secondary metabolite biosynthetic capabilities, amino acid metabolic capabilities, and a number of uncharacterized genes detected in the genomes of x. bovienii strains maybe be discernable by host nematodes  <cit> .
this was supported by the fact that s. carpocapsae, which does not carry x. bovienii, was not attracted to any x. bovienii bacteria assayed in this study, with the exception of xbo .
steinernema oregonense followed a similar pattern, as it performed best with its cognate symbiont in all metrics assessed.
for example, we showed that s. oregonense-xbp and -xbkca combinations had moderate success in terms of progeny production, longevity, and bacterial carriage but failed to restore performance to the level of associations with the cognate symbiont .
additionally, s. oregonense was attracted to xbp, although not significantly , but perhaps this reflects the nematode’s ability to detect some indication of its benefit over other less advantageous bacteria.
it could then be hypothesized that there are phenotypes of bacterial symbionts that the nematode can discriminate and react to behaviorally, such as chemical cues or signals they can perceive.
this phenomenon has been observed in associations of many taxa where host/partner finding is an important part of the symbiosis.
for example, leaf-cutting ants discriminate between concolonial, allocolonial, and allospecific fungi which differed measurably in the volatile, chemical signatures  <cit> .
beewolves recognize the colonization by opportunistic, non-symbiotic bacteria and do not pass those bacteria on to their offspring into brood cells  <cit> .
in choice assays, both ambrosia beetles and hermit crabs using olfactory-based chemical discrimination to locate their symbiotic partners, fungi and sea anemones respectively,  <cit> .
therefore, it is plausible that partner fidelity in entomopathogenic nematode-bacterium symbioses may be, in part, chemically mediated.
in contrast, our data showed that s. puntauvense was symbiotically promiscuous, demonstrating comparable infection success and progeny production with any x. bovienii feltiae bacterium  as it did with xbp .
we also showed that s. puntauvense was attracted to two non-cognate symbionts, xbfsn and xbi, in bacterial preference assays.
the preference of xbfsn may reflect phylogenetic relatedness of xbp to xbf* bacteria, suggesting that this clade may share common features which facilitate success in the first generation of cross-association but cannot be maintained in subsequent generations.
this was manifested in the observed limited virulence and longevity of f1 s. puntauvense-xbf* combinations , and may be also due to insufficient bacterial carriage in the ijs .
these findings are supported by the attraction assays in s. feltiae  and other measures of fitness and virulence  <cit> .
the demonstrated partner fidelity we see in these steinernema species coupled with the resulting multi-generational benefits of cognate associations are important for symbiotic stasis.
it is possible for an insect to become co-infected with more than one nematode species  simultaneously.
in doing so, nematodes must be able to identify their cognate symbiont or risk the fitness costs as manifested in the present study.
to complicate things further, x. bovienii strains exhibit differential virulence within an insect host  <cit> .
these observations suggest certain x. bovienii strains can directly outcompete other less virulent strains  <cit> .
additionally, it is possible that distinct strains of x. bovienii may have distinct nutritional requirements and/or differences in metabolic capability which could impact the nematode carrying them.
thus, the ability to discriminate native bacterial partners within a mixed infection is imperative for nematode success.
alternatively, nematode discrimination may be a mechanism to prevent co-infection events thereby reducing within host competition between different xenorhabdus strains.
previous evidence suggests that some steinernematid ijs respond to volatile chemical cues from insect hosts infected by both conspecifics and heterospecifics, and this chemical profile is influenced by the xenorhabdus symbionts  <cit> .
however, previous phylogenetic data leads to the hypothesis that host-symbiont switching events are rampant throughout the steinernema-xenorhabdus co-phylogeny which would suggest the ability of ijs to discern xenorhabdus partners that are best suited for nematode success.
steinernema-xenorhabdus symbioses exhibit all the hallmarks of mutualism increased fitness for both partners, divided but equally important roles in cooperative infections, and preference for native symbionts  <cit> .
additionally, the fitness costs of steinernema-xenorhabdus non-cognate associations are predictable based on phylogeny .
CONCLUSIONS
entomopathogenic nematodes and their bacterial symbionts provide an ideal model system for studying mutualism because of their manipulability, short generation times, and ease of rearing.
these stable associations carry measurable fitness benefits, have strong partner preference, and are phylogenetically predictable.
we have demonstrated the effects of partner-switching manifest across three species of steinernema nematodes which host different x. bovienii strains as symbionts.
each nematode species performed best in terms of reproductive fitness and virulence when associated with its cognate symbiont.
ijs colonized with native partners harbor more bacteria per nematode thus increasing the likelihood of the bacterium being passaged to offspring in a new insect host successfully.
finally, these associations were reinforced by strong partner preference evidence, shown in the symbiont choice assays, meaning x. bovienii strains produce reliable cues or signals that steinernema nematodes were able to detect and distinguish to varying degrees.
using the entomopathogenic nematode-bacterium association as a model for mutualism provides a practical avenue to study the ecological and evolutionary consequences of co-evolved symbioses.
additional file
additional file 1:
figure s <dig>  effect of non-cognate symbionts on nematode body size.
fig.
s <dig> schematic representation of setup for preference assays.
table s <dig> bacterial gene labels of sequences used for bacterial phylogeny.
table s <dig> summary results from cox mixed effects models.
table s <dig> summary statistics for choice assays.
table s <dig> summary results of linear regression models.
abbreviations
aaiattractions activity index
anovaanalysis of variance
cfucolony forming unit
dnadeoxyribonucleic acid
f1first generation progeny
ijinfective juvenile
lbluria-bertani
lt50median lethal time
nbtanutrient agar supplemented with  <dig> %  bromothymol blue,  <dig> %  triphenyl-tetrazolium, and  <dig> %  sodium pyruvate
odoptical density
rdnaribosomal dna
rnaribonucleic acid
rrnaribosomal rna
taftriethanolamine formalin
xbf*any x. bovienii strain associated with any s. feltiae strain
xbffl
x. bovienii associated with s. feltiae fl
xbfm
x. bovienii associated with s. feltiae moldova
xbfsn
x. bovienii associated with s. feltiae sn
xbi
x. bovienii associated with s. intermedium type
xbj
x. bovienii associated with s. jollieti monsanto
xbkbu
x. bovienii associated with s. kraussei nemasys-l
xbkca
x. bovienii associated with s. kraussei quebec
xbo
x. bovienii associated with s. oregonense os-10
xbp
x. bovienii associated with s. puntauvense li6
xn
x. nematophila associated with s. carpocapsae all

 every plant and animal on this planet is engaged in beneficial and antagonistic relationships with microbes that profoundly influence their physiology, life history, and evolution. the ubiquity and biological importance of long-term, stable, and beneficial host-microbe associations, or mutualisms, have been the focus of many investigations. models of mutualisms possess common factors that contribute to their success, maintenance, and perpetuation. cooperative association, fitness benefits for all parties, partner fidelity and preference are key factors that influence the establishment and maintenance of mutualisms. in this respect, the relationship between gram-negative xenorhabdus bacteria and their steinernema nematode hosts yields insights into processes that influence symbiont specificity, selection, and function. furthermore, the inherent experimental tractability of these organisms and availability of genome sequences has spurred developments across a broad scientific front, such that this mutualism is now viewed as a model for basic research in ecology, evolution, biochemistry, and molecular genetics of symbiosis. briefly, the xenorhabdus life cycle depends on two eukaryote hosts, steinernema nematodes and an insect for their survival, dissemination, and reproduction. the relationship between the bacteria and the nematodes is benign and mutually beneficial, whereas their relationship with the insect host is pathogenic. after invading an insect, nematodes defecate xenorhabdus symbionts into the hemocoel where they multiply rapidly, neutralize the insect’s immune system, induce septicemia and toxemia, and begin digesting the fat-rich insect cadaver. within the compromised insect, nematodes mature, and sexually reproduce resulting in many thousands of offspring using the metabolic resources liberated from the cadaver. once resources are exhausted, the nematodes develop into a non-feeding, environmentally-resistant “dauer” stage, known as the infective juvenile. prior to emerging from the insect cadaver, ijs re-associate with xenorhabdus symbionts, which colonize the ij nematodes’ intestinal receptacle. ijs carrying bacterial symbionts persist in the soil until a new insect host can be infected. nutrients provisioned by the nematode ij are thought to support outgrowth and bacterial survival until the next insect host. each steinernema species associates with a single xenorhabdus species. however, one xenorhabdus species can form associations with more than one steinernema host. the most promiscuous xenorhabdus species is xenorhabdus bovienii, which associates with numerous steinernema species comprising two different evolutionary clades in the steinernema phylogeny. previous studies revealed that the bacteria can contribute to virulence toward the insect host and support nematode reproduction. these benefits have also been demonstrated to have a strong phylogenetic correlation at the nematode intraspecific level with the driving selective interest of nematode symbiont choice likely differing between nematode species. complicating matters, in nature, insects can become co-infected with two steinernema species each harboring their respective bacterial symbiont. this requires that nematodes disassociate with their bacterial partner, reproduce and mature within a mixed community, and then re-associate with the correct symbiont to avoid the fitness costs associated with harboring the wrong bacterium. this sets entomopathogenic nematode-bacterium symbioses apart from many other mutualistic systems where partners do not naturally disassociate under optimal environmental conditions. symbiont switching can have detrimental effects to animal hosts. in examples across invertebrate-microbe symbioses, symbionts serve important functions and provide adaptive features for their hosts by aiding in digestion, conferring pathogen/parasite resistance, and increasing host fitness. in fact, within a single species of nematode, association with a closely-related non-cognate symbiont results in reduced fitness. for this reason, partner fidelity and localization are especially important to host animals. steinernema-xenorhabdus symbioses are excellent models for testing hypotheses concerning mutualism. obligate mutualisms evolve and persist due to increased fitness of both parties through maintenance of partner fidelity and partner preference. unlike many intracellular, vertically transmitted animal-microbe symbioses, entomopathogenic nematode-bacterial symbioses rely on partner re-association after each infection cycle. due to the tractability of this system, we were able to test the fitness, fidelity, and preference of several nematode-bacterial partner pairs by switching their symbionts. we show that across a number of x. bovienii-steinernema spp. combinations, cognate partners display the hallmarks of mutualism including increased fitness in cooperation, partner fidelity, and partner preference. additionally, we show that the most derived steinernema species tested s. puntauvense, is more promiscuous and is able to form relatively successful relationships with a wide variety of non-cognate xenorhabdus species. contrastingly, the ancestral x. bovienii host, s. intermedium, is the most fastidious of the nematode species examined, performing poorly with all non-cognate bacteria. xenorhabdus bovienii is attractive to the largest variety of non-cognate nematode species. taken together, these data show entomopathogenic nematode-bacterial mutualisms represent an important and useful model for understanding the evolution of mutualisms. methods nematode and bacteria isolates three steinernema species, which are hosts of xenorhabdus bovienii, were used in this study to assess the cost/benefit of symbiont switching. two of these species were assessed in symbiont choice assays to determine symbiont preference. additionally, steinernema feltiae and steinernema carpocapsae were included in these assays given that these taxa have been extensively studied for the effects of symbiont switching events. details on nematode hosts and bacterial symbionts species/strain names, geographic origin, and sequence accession numbers are listed in table identity of nematodes and bacteria considered in this study was molecularly confirmed prior to the initiation of the experiments by sequencing of 28s and 12s ribosomal dna genes for the nematode hosts and 16s genes for their bacterial symbionts.table steinernema and xenorhabdus species and strains considered in this study s. carpocapsae x. nematophila   s. feltiae x. bovienii   s. feltiae x. bovienii   s. feltiae x. bovienii s. intermedium x. bovienii   s. kraussei x. bovienii   s. kraussei x. bovienii   s. jollieti x. bovienii s. oregonense x. bovienii s. puntauvense x. bovienii nematodes were reared in vivo with last instar wax moth larvae, galleria mellonella following procedures described by kaya and stock. modified white traps were used to harvest emerging infective juvenile stages, which were stored in tissue culture flasks at a concentration of approximately ijs/ml in 250 ml tissue culture flasks. bacterial symbionts were extracted from each nematode host species/strain by sonication of naturally colonized ijs. bacterial stocks were established by growing each bacterial species/strain to mid-log phase in luria-bertani broth supplemented with % sodium pyruvate grown in a shaking, dark incubator at 28 °c following procedures described by akhurst. bacterial stocks were stored at −80 °c in lb broth supplemented with 20% glycerol. rearing of aposymbiotic nematodes aposymbiotic nematodes were obtained via alkaline axenization of eggs harvested from first generation gravid females as described by mcmullen and stock. eggs were seeded on liver-kidney agar. aposymbiotic ij nematodes were collected by removing the bottom portion with agar of the petri dish and placing it in a sterile modified white trap. ijs were stored in tissue culture flasks at 15 °c until use for experiments following the same procedures described above. the aposymbiotic nature of the ijs was confirmed by grinding them during the first and second week of emergence with a motorized pestle in lb broth. the suspension was then plated onto differential nutrient agar supplemented with % bromothymol blue, % triphenyl-tetrazolium, and % sodium pyruvate to confirm presence/absence of bacterial symbionts. host-symbiont interactions assays host virulence aposymbiotic ijs were surface sterilized according to procedures described by stock and goodrich-blair. ij suspension as diluted in grace’s medium to a concentration of ijs in 10 μl. bacteria cultures were grown overnight and then sub-cultured to a 1: dilution into lb broth and subsequently grown to od600nm ~. bacterial suspension was then diluted in grace’s insect medium to a concentration 200 cfu in 10 Μl. prior to injections, bacterial suspensions were plated to confirm that equal counts were used for all bacterial strains. the experimental arena was a 12-well plate. each g. mellonella larvae was co-injected with a nematode/bacterium treatment. treatments consisted of co-injections of different steinernema/xenorhabdus species/strain pairs. bacteria and nematodes were co-injected at the concentrations described above. twelve g. mellonella larvae were used per treatment and each treatment were repeated three times independently. insect mortality was monitored daily for 7 days post-injection. g. mellonella larvae were considered dead when they did not respond to the gentle prodding of soft-grip forceps. resulting cadavers and progeny were used for all subsequent experiments, except for the choice assays.table cognate and non-cognate nematode-symbiont pairs considered in this study grey boxes denote cognate symbiont-nematode pair host reproductive fitness reproductive fitness of s. intermedium, s. oregonense and s. puntauvense was measured considering two parameters: a) infection productivity, and b) emerging ij progeny. for productive infection percentages, g. mellonella cadavers obtained from the virulence assays were considered. each cadaver was thoroughly rinsed and transferred into a modified white trap at day seven post-injection. cadavers were checked daily to record progeny emergence for up to 30 days. cadavers producing at least ijs were scored as ‘productive’, and productive infection percentage was calculated as the number of cadavers producing progeny out of the total number. the emerging ij progeny was determined by collecting the nematode suspension from each modified white trap daily for 10 days of post-emergence and following volumetric dilution procedures per kaya and stock. symbiont carriage ij symbiont carriage was determined following procedures described by goestch et al.. two technical replicates and three biological replicates were performed for each nematode-bacterium combination and averaged to obtain cfu ij− f virulence virulence of the emerging ij progeny was tested considering a one-on-one assay as described by kaya and stock. the experimental arena was a 12-well plate, where each well was filled with 5 g of sterile sand. an inoculum of ijs delivered in 100 μl of distilled water was poured onto the sand. one g. mellonella larvae was added into each well. the lid was placed and the plate was incubated at 20 °c in the dark. mortality was monitored every 8 h for a total of 96 h. experiments were repeated three times. for data analysis, percent mortality registered at 48 h was used, as it showed the largest difference between combinations. longevity of f progeny ijs obtained from the host reproductive fitness assays were examined 6-months post-emergence at 30× magnification. nematodes were considered dead when they adopted a straight shape and/or showed no movement, whereas nematodes were considered live when they were moving or had a j-shape position. ijs were confirmed as dead by gentle prodding. two technical replicates of at least ijs were counted from each productive infection. morphometric analysis the effect of cognate/non-cognate bacterial symbiont on nematode size was assessed by measuring length and width of ij progeny obtained from each treatment. a total of ijs/ treatment were randomly pooled. ijs were heat killed in m buffer at 50-60 °c and subsequently fixed with triethanolamine formalin following procedures described by goodey. ijs were mounted onto glass slides and measured using an olympus bx microscope equipped with differential interference contrast optics and olympus microsuite software. bacteria choice assays choice assays were performed to test whether different steinernema hosts of x. bovienii could discern between bacterial strains that are cognate vs. non-cognate. for this purpose, three steinernema hosts of x. bovienii were considered: s. feltiae, s. oregonense, and s. puntauvense. in addition, two other x. bovienii strains were included, xbffl and xbi to account for symbionts that have a range of fitness effects on nematode hosts. s. carpocapsae and its respective symbiont x. nematophila were considered for outgroup comparisons. for each nematode species used, a no-choice assay was performed where all bacterial foci in our test arena consisted of the symbiotic partner as a control. aposymbiotic ijs were reared following procedures described by mcmullen and stock. prior to the experimental setup, nematodes were surface sterilized using the vacuum filtration method described by kaya and stock. due to low fitness of aposymbiotic s. intermedium nematodes, bacterial preference assays did not include this species as this would likely confound results. bacterial strains were grown from glycerol stocks overnight in 5 ml of lb broth with % sodium pyruvate in a 15 ml centrifuge tube at 28 ± 1 °c in the dark with agitation. 50 μl of bacterial suspension was sub-cultured following the same conditions, but allowed to grow for 9 h prior to assay setup. cultures were streaked onto nbta to confirm liquid cultures were phase i under the same conditions described above. briefly, a 10 cm petri dish was filled with 20 ml of lipid agar supplemented with % sodium pyruvate. a bacterium inoculum of 10 μl was inoculated onto the agar plate in four equidistant positions, approximately  cm away from the edge. for each plate, one cognate and one non-cognate bacterial strain were applied so that two inoculum positions of the same strain were not next to each other. a 50 μl inoculum of ijs was dispensed into the center of the plate. there were five replicates for each nematode-bacteria combination. controls consisted of four foci of cognate bacteria. plates were incubated at 20-25 °c for 3 days in the dark. this time period allowed for enough time for ijs to migrate to the preferred bacteria. experiments were repeated twice, independently. after 3 days, nematode migration to bacteria foci was evaluated. to facilitate visualization and counting of nematodes, bacteria foci and nematode inoculum area were cut and placed into a 5 cm petri dish containing approximately 10 ml of m buffer. dishes were then placed on a shaker overnight to expedite release of the nematodes embedded in the agar and bacterial lawn. the entire plate was then rinsed off to collect the ijs and stored in 50 ml centrifuge tubes at 4 °c until counted. an attraction activity index was calculated to assess nematode attraction for bacteria in the choice assays. for each treatment, the number of nematodes found in each bacterial focus was combined. similar to zhang et al., the formula below was used to obtain the aai values. aai ranges from to − and the resulting values can be inferred as: 1 = all nematodes are attracted to non-cognate symbiont, 0 = equally prefer both bacteria, −1 = all nematode are repulsed by non-cognate strain. the formula used was: aai =  – / total number of nematodes in bacteria. for no-choice assay controls, bacterial foci pairs were randomly assigned to follow data analysis procedure used for choice assays. phylogenetic analyses phylogenetic analyses were conducted to determine relative phylogenetic distance between cognate and non-cognate xenorhabdus strains. mega v. was used to generate both nematode and bacterial phylogenies. for the xenorhabdus phylogeny, the following sequences were aligned using muscle and concatenated: 16s rdna, reca, dnan, gltx, gyrb, and infb. maximum likelihood was carried out using a tn93 + g + i model of evolution with bootstrap replicates. the nematode phylogeny was generated using similar methods as the bacterial phylogeny with the mitochondrial and nuclear rrna genes: 12s and 28s, respectively. however, a gtr + g model of evolution was used. s. carpocapsae and its native symbiont x. nematophila were considered in nematode and bacterial phylogenies, respectively, as out group taxa. statistical analyses all statistical analyses were carried out using r, version.. nematode-bacterium virulence against g. mellonella were analyzed using a cox mixed effects model from the package ‘coxme’. each nematode host tested was analyzed separately and their cognate symbiont was used as the reference strain to build the cox model. the bacterial strain was considered as a fixed effect in the model, while a random effect was included to better explain the variance due to differences between experimental setups. a post hoc tukey hsd test was conducted using the ‘multcomp’ package. final insect mortality at day seven and median mortality time were calculated from data assessed with the ‘survival’ package. if a nematode-bacterium combination resulted in no deaths during the study period, a single death was added at day seven to allow the model to be stable. either a one-way analysis of variance or mixed effects anova was considered to analyze most assays. a one-way anova was considered to analyze productive infection and f progeny data. in both instances, a single percentage was generated for each nematode-bacterium pair from each experimental setup. the remaining measures were assessed with a mixed effect model. for all analyses, each nematode host was analyzed separately and the bacterial treatment was used as a fixed effect. when a mixed effect anova was used, the experimental setup was used as a random effect to account for setup variation with the packages ‘lme4’ and ‘car’. all residuals were visually inspected for normality and homoscedasticity. a post hoc tukey hsd test was performed for pairwise comparisons between bacterial treatments with the package ‘lsmeans’. nematode choice data was analyzed by performing a two-sided one-sample t-test for each aai index with a μ = . p-values were adjusted using the benjamini-hochberg method to correct for false discovery rate in each nematode dataset. linear regression analyses were conducted on nematode and bacterial fitness and morphometric data. nematode fitness was calculated as a factor of progeny count, progeny survival, and productive infection percentage relative to the cognate symbiont mean, while bacterial fitness was determined as a factor of nematode fitness and cfu ij− relative to the average cognate symbiont mean. for morphometric analyses, the factor of length and width was analyzed. each nematode host was analyzed separately.  all parameters measured across nematode-symbiont combinations were influenced by the pairing at hand with the exception of nematode size. together the measured variables demonstrate the truly mutualistic nature of associations between steinernema nematodes and their xenorhabdus symbionts. some non-cognate symbionts are sufficient to maintain virulence against insect hosts. all steinernema spp. tested demonstrated greater virulence with their cognate symbionts than with any other, non-cognate bacterium. of them, symbiont switching in s. intermedium, a clade i species, had the most pronounced impact on virulence. association with the cognate symbiont, xbi, with s. intermedium resulted in 50% insect mortality in 4 days and ultimately killed 100% of insect hosts within a week. s. intermedium was also capable of moderate levels of virulence with both xbkbu and xbj resulting in % and % mortality by day seven, respectively. association with xbffl and xbp did result in some insect mortality after 1 week, although lt was not observed within the 7-day observation period, and virulence when s. intermedium was colonized by these species was not statistically distinguishable from aposymbiotic s. intermedium or those species that were completely incompetent.table non-cognate x. bovienii effect on nematode virulence x. bovienii s. intermedium a s. oregonense a s. puntauvense a % mortalityb  ± 0 4 a  ±  2 a  ±  2 a acognate symbiont data are bolded bpercent mortality after 1 week with standard error clethal time to kill 50% of population reported as days dletter-ranking generated from a cox mixed effects model and post hoc tukey test; see additional file 1: table s for summary of model contrastingly, both assayed clade iii nematode species, s. puntauvense and s. oregonense, maintained some level of virulence with all tested bacteria and even without a bacterial partner. as with s. intermedium, s. oregonense was most virulent when associated with their cognate symbiont; lt was shorter and 7-day insect mortality is higher in associations with xbo. s. oregonense and s. puntauvense exhibited moderate virulence with all x. bovienii tested and with no associated symbiont. interestingly, s. puntauvense had comparable lt times with many symbionts, though its cognate still produced the highest insect mortality at 7-days post-inoculation. symbiont switching has negative impacts on fitness despite many steinernema-xenorhabdus combinations producing successful infections in insect hosts, associations with non-cognate symbionts impacted nematode reproductive success. the first measure of this impact was reflected in the proportion of productive infections from a given nematode-bacterium association. we determined an infection to be productive if it resulted in the emergence of at least ijs. all species suffered negligible infection productivity when reared aposymbiotically and non-cognate associations were not sufficient to recover that loss in s. intermedium and s. oregonense. s. intermedium had the least productive infections with non-cognate bacteria, with only xbkbu partially restoring infection productivity, and all combinations were not significantly different from the uncolonized aposymbiotic worms. s. intermedium produced small numbers of ijs when partnered with xbj and xbkbu. associations between non-cognates xbkca and xbp generated a comparable proportion of productive infections to the cognate symbiont in s. oregonense; however, the total number of progeny from these infections was significantly less than what was measured from xbo associations. in contrast, no progeny was recovered when s. puntauvense associated with non-cognate bacterium, even though most combinations resulted in productive infections. however, this pattern was not seen in combinations with xbi, xbj, xbkbu, nor xbo.fig. productive infections resulting from steinernema-xenorhabdus experimental associations. bars indicate the mean and se resulting from a one-way anova. letters above bars show statistical ranking from post hoc tukey test. n =  for each nematode-bacterium treatment. these values were calculated as the percent of cadavers that produce progeny in each independent experimental setup. cognate symbionts are highlighted with a heavy outline. a. s. intermedium: f 20 = , p < , r 2 = ; b. s. oregonense: f 20 = , p < , r 2 = ; and c. s. puntauvense: f 20 = , p , r 2 =  fig. steinernema progeny production with novel symbionts. bars show the mean and se resulting from a mixed effects anova. letters above bars demonstrate results from post hoc tukey test. n =  for each nematode-bacterium treatment. cognate symbionts are highlighted with a heavy outline. a. s. intermedium: f 313 = , p < , r 2 o = ; b. s. oregonense: f 313 = , p < , r 2 o = ; and c. s. puntauvense: f 313 = , p < . r 2 o =  we also tested the longevity and virulence of the f ijs that emerged following association with non-cognate bacteria. in general, the patterns observed in ij longevity followed the trends seen in progeny production. nematode-symbiont pairings that yielded higher progeny production, also tended to have increased longevity. six-months post-emergence, s. intermedium-xbj and -xbkbu pairings had ~20% ij survival, whereas s. oregonense-xbkca and -xbp pairings exhibited ~25% ij survival. while these are both significantly less than the cognate symbiont, these pairings had higher survival than the other non-cognates that were capable of producing progeny. all symbiont pairings with s. puntauvense, except xbi, had viable ijs 6-months post-emergence, though none of the non-cognate pairings survived as well as those associated with the cognate symbiont.fig. steinernema progeny survival with non-cognate symbionts. bars indicate the mean and se resulting from a mixed effects anova. letters above bars show results from post hoc tukey test. n =  for each nematode-bacterium treatment. data represents percentage of nematode alive in a pool of ijs. cognate symbionts are highlighted with a heavy outline. a. s. intermedium: f 313 = , p < , r 2 o = ; b. s. oregonense: f 313 = , p < , r 2 o = ; and c. s. puntauvense: f 313 = , p < , r 2 o =  in addition to f ij longevity, virulence of this new generation is imperative for nematode success long term and establishment of novel associations. strikingly, f ijs from all non-cognate associations had significantly less virulence than cognate associations. in combinations where progeny were produced some virulence was observed no matter how slight. consistent with other observations, s. intermedium was most susceptible to symbiont switching and resulted in a maximum of ~30% virulence in f ijs; whereas f1 s. puntauvense ijs had the most success with non-cognate symbionts and resulted in several combinations causing insect mortality >40%. s. oregonense had intermediate success out of the three nematodes tested, with an observed insect mortality below 40%.fig. f progeny virulence against g. mellonella. bars indicate the mean and se resulting from a one-way anova. letters above bars show results from post hoc tukey test. n =  for each nematode-bacterium treatment. values were calculated as the percent of insects that died by 48 h in each independent experimental setup. cognate symbionts are highlighted with a heavy outline. a. s. intermedium: f 18 = , p = , r 2 = ; b. s. oregonense: f 18 = , p < , r 2 = ; and c. s. puntauvense: f 18 = , p , r 2 = . post- non-cognate bacteria are less successful at colonizing ijs we also assayed the extent to which bacteria were able to establish themselves in cognate and non-cognate ijs by measuring bacterial carriage. by and large, bacterial carriage follows the same patterns we saw in the other metrics discussed above. s. intermedium carried the fewest non-cognate bacteria, s. oregonense can carry more species of non-cognate bacteria but at significantly fewer numbers, and s. puntauvense carried more non-cognate bacteria than the other two species though at lower concentrations than their cognate xbp.fig. steinernema symbiont carriage in f populations. bars indicate the mean and se resulting from a mixed effects anova. letters above bars show results from post hoc tukey test. n =  for each nematode-bacterium treatment. cognate symbionts are highlighted with a heavy outline. a. s. intermedium: f 311 = , p < , r 2 o = ; b. s. oregonense: f 313 = , p < , r 2 o = ; and c. s. puntauvense: f 313 = , p < , r 2 o =  x. bovienii hosts discern between some non-cognates in bacterial choice assays using choice assays, we determined that steinernema ijs were able to discern between some conspecific non-cognate bacterial strains and are repelled by heterospecific bacteria. for example, s. carpocapsae was not attracted to any of the x. bovienii strains, with the exception of xbo. the nematodes appeared to not be able to distinguish between this bacterium and their cognate partner. similarly, natural x. bovienii hosts were not attracted to xn. steinernema species that naturally associate with x. bovienii exhibited strain-specific behavioral preferences within x. bovienii. for instance, s. feltiae preferred xbffl and xbp equally as it did for its cognate symbiont. this nematode was also attracted to xbi but repelled xbo. s. oregonense was not attracted to xbffl and xbfsn, and did not demonstrate a differentiated preference for xbp and xbi and its cognate. s. puntauvense did not display a preference for any of the tested conspecific bacteria. interestingly, xbi was the only non-cognate bacterium that was always preferred by the three x. bovienii nematode host species tested in this study.fig. nematode preference for non-cognate symbionts. a. s. feltiae, b. s. oregonense, c. s. puntauvense, and d. s. carpocapsae. bars represent the mean and se resulting from two-sided one-sample t-test analyses with μ =  p-values were adjusted with benjamini-hochberg method to correct for the false discovery rate and their results are indicated by the following symbols: ns = non-significant, * = p < , ** = p < , *** = p < , **** = p < . n =  for each aai value. white bars indicate nematode attraction to non-cognate symbiont, while gray bars show repulsion. cognate symbiont bars are highlighted with a heavy outline. see additional file 1: table s for summary statistics phylogeny as a predictor of fitness effects of non-cognate associations a co-phylogenetic hypothesis for all bacteria-nematode pairs considered in this study was developed. the overall fitness effects of all nematode-bacterial pairings were regressed against the phylogenetic distance between the test bacterium and the cognate bacterium. in all cases, we show that fitness of both bacteria and nematodes declines as the phylogenetic distance between the cognate and non-cognate bacteria increases. these results suggest that the more distant a bacterial strain is in relation to a nematode’s cognate symbiont, the more detrimental the fitness effects are on non-cognate nematode partner. in particular, s. intermedium and s. oregonense are more sensitive to association with a phylogenetically distant symbiont, with - unit decrease in nematode and bacterial fitness for every nucleotide change in symbiont. s. puntauvense is more tolerant to association with non-cognates and only suffers a - unit decrease in fitness for every bacterial nucleotide change.fig. evolutionary relatedness between steinernema nematodes and their native symbionts and the fitness consequence of host-symbiont switching events. a. co-phylogeny tanglegram of steinernema-xenorhabdus pairs used in experiments. the steinernema phylogeny was generated with concatenated 12s and 28s rrna genes, while the xenorhabdus phylogeny shows the results from a concatenation of the following genes: 16s rrna, reca, dnan, gltx, gyrb, and infb. internal nodes indicate bootstrap replicate values >50%. panels b and c show relative nematode and bacterial fitness, respectively, by bacterial phylogenetic distance from native steinernema-xenorhabdus pairing. for both graphs, data points indicate mean and se from a linear regression analysis; see additional file 1: table s for model summary. s. intermedium = black triangles with black line, s. oregonense = gray squares with gray line, and s. puntauvense = hashed diamonds with dashed line discussion entomopathogenic nematodes and their bacterial partners are an important model system for studying mutualism. the data presented here demonstrate that steinernema-xenorhabdus symbioses are mutually beneficial with regard to nematode and bacterial fitness across multiple species of nematode hosts. association with cognate symbionts benefits the nematode host in terms of virulence against insect hosts, total progeny production, f progeny longevity, and virulence. the bacterial partners also benefit by increased host colonization, and thus have an increased chance of being vectored into a new insect host, by their nematode host. the final piece of this puzzle is partner fidelity. not only did nematodes discern among different bacterial strains to varying degrees, but they tended to prefer bacteria which ultimately increase their virulence, longevity, and reproductive success. generally, all cognate nematode-bacteria pairs outperformed all other combinations in terms of virulence, progeny production and longevity, bacterial carriage, and fitness of both parties. with the exception of nematode size, all measured variables were influenced by symbiont partner associated with the given nematode host. nematodes were also capable of discerning their cognate bacterium in a choice assay and tended to prefer their native symbiont to heterospecific strains, while their attraction toward non-cognate symbionts mirrored fitness and virulence trends measured in this study. most importantly, the impacts for both nematode and bacterial fitness were correlated to the phylogenetic distance between the cognate bacterium and that colonizing partner. these data solidify the argument that entomopathogenic nematodes and their bacterial symbionts represent a true mutualistic symbiosis. while this trend is consistent across species, each tested nematode species had a distinct pattern of effects when associated with various x. bovienii strains. for example, s. intermedium, the only species member of clade i tested, was very fastidious. this nematode showed the highest virulence, progeny production, and f longevity with its native symbiont, xbi. however, these same metrics were only minimally restored when the nematodes associated with xbkbu and xbj; no non-cognate symbiont restored even 50% virulence, ij production, or progeny success as the native xbi. unfortunately, due to technical limitations, we were not able to determine bacterial preference in s. intermedium, but we did assay how other steinernema species responded to xbi. interestingly, xbi was attractive to non-cognate nematode species such as s. feltiae, s. oregonense or s. puntauvense, despite conferring almost no fitness or virulence benefits. given that its host, s. intermedium is in a basal phylogenetic clade to these other x. bovienii hosts tested, xbi may exhibit some attractive trait which is generic to all x. bovienii, and as such is detectable by all x. bovienii host nematodes. in this respect, it could be speculated an attraction trait, present in xbi, which more derived nematode hosts are cuing in on, may also be present in the cognate x. bovienii lineages. it has been shown that differences in primary and secondary metabolite biosynthetic capabilities, amino acid metabolic capabilities, and a number of uncharacterized genes detected in the genomes of x. bovienii strains maybe be discernable by host nematodes. this was supported by the fact that s. carpocapsae, which does not carry x. bovienii, was not attracted to any x. bovienii bacteria assayed in this study, with the exception of xbo. steinernema oregonense followed a similar pattern, as it performed best with its cognate symbiont in all metrics assessed. for example, we showed that s. oregonense-xbp and -xbkca combinations had moderate success in terms of progeny production, longevity, and bacterial carriage but failed to restore performance to the level of associations with the cognate symbiont. additionally, s. oregonense was attracted to xbp, although not significantly, but perhaps this reflects the nematode’s ability to detect some indication of its benefit over other less advantageous bacteria. it could then be hypothesized that there are phenotypes of bacterial symbionts that the nematode can discriminate and react to behaviorally, such as chemical cues or signals they can perceive. this phenomenon has been observed in associations of many taxa where host/partner finding is an important part of the symbiosis. for example, leaf-cutting ants discriminate between concolonial, allocolonial, and allospecific fungi which differed measurably in the volatile, chemical signatures. beewolves recognize the colonization by opportunistic, non-symbiotic bacteria and do not pass those bacteria on to their offspring into brood cells. in choice assays, both ambrosia beetles and hermit crabs using olfactory-based chemical discrimination to locate their symbiotic partners, fungi and sea anemones respectively,. therefore, it is plausible that partner fidelity in entomopathogenic nematode-bacterium symbioses may be, in part, chemically mediated. in contrast, our data showed that s. puntauvense was symbiotically promiscuous, demonstrating comparable infection success and progeny production with any x. bovienii feltiae bacterium as it did with xbp. we also showed that s. puntauvense was attracted to two non-cognate symbionts, xbfsn and xbi, in bacterial preference assays. the preference of xbfsn may reflect phylogenetic relatedness of xbp to xbf* bacteria, suggesting that this clade may share common features which facilitate success in the first generation of cross-association but cannot be maintained in subsequent generations. this was manifested in the observed limited virulence and longevity of f1 s. puntauvense-xbf* combinations, and may be also due to insufficient bacterial carriage in the ijs. these findings are supported by the attraction assays in s. feltiae and other measures of fitness and virulence. the demonstrated partner fidelity we see in these steinernema species coupled with the resulting multi-generational benefits of cognate associations are important for symbiotic stasis. it is possible for an insect to become co-infected with more than one nematode species simultaneously. in doing so, nematodes must be able to identify their cognate symbiont or risk the fitness costs as manifested in the present study. to complicate things further, x. bovienii strains exhibit differential virulence within an insect host. these observations suggest certain x. bovienii strains can directly outcompete other less virulent strains. additionally, it is possible that distinct strains of x. bovienii may have distinct nutritional requirements and/or differences in metabolic capability which could impact the nematode carrying them. thus, the ability to discriminate native bacterial partners within a mixed infection is imperative for nematode success. alternatively, nematode discrimination may be a mechanism to prevent co-infection events thereby reducing within host competition between different xenorhabdus strains. previous evidence suggests that some steinernematid ijs respond to volatile chemical cues from insect hosts infected by both conspecifics and heterospecifics, and this chemical profile is influenced by the xenorhabdus symbionts. however, previous phylogenetic data leads to the hypothesis that host-symbiont switching events are rampant throughout the steinernema-xenorhabdus co-phylogeny which would suggest the ability of ijs to discern xenorhabdus partners that are best suited for nematode success. steinernema-xenorhabdus symbioses exhibit all the hallmarks of mutualism increased fitness for both partners, divided but equally important roles in cooperative infections, and preference for native symbionts. additionally, the fitness costs of steinernema-xenorhabdus non-cognate associations are predictable based on phylogeny.  entomopathogenic nematodes and their bacterial symbionts provide an ideal model system for studying mutualism because of their manipulability, short generation times, and ease of rearing. these stable associations carry measurable fitness benefits, have strong partner preference, and are phylogenetically predictable. we have demonstrated the effects of partner-switching manifest across three species of steinernema nematodes which host different x. bovienii strains as symbionts. each nematode species performed best in terms of reproductive fitness and virulence when associated with its cognate symbiont. ijs colonized with native partners harbor more bacteria per nematode thus increasing the likelihood of the bacterium being passaged to offspring in a new insect host successfully. finally, these associations were reinforced by strong partner preference evidence, shown in the symbiont choice assays, meaning x. bovienii strains produce reliable cues or signals that steinernema nematodes were able to detect and distinguish to varying degrees. using the entomopathogenic nematode-bacterium association as a model for mutualism provides a practical avenue to study the ecological and evolutionary consequences of co-evolved symbioses. additional file additional file 1: figure s effect of non-cognate symbionts on nematode body size. fig. s schematic representation of setup for preference assays. table s bacterial gene labels of sequences used for bacterial phylogeny. table s summary results from cox mixed effects models. table s summary statistics for choice assays. table s summary results of linear regression models. abbreviations aaiattractions activity index anovaanalysis of variance cfucolony forming unit dnadeoxyribonucleic acid f1first generation progeny ijinfective juvenile lbluria-bertani lt50median lethal time nbtanutrient agar supplemented with % bromothymol blue, % triphenyl-tetrazolium, and % sodium pyruvate odoptical density rdnaribosomal dna rnaribonucleic acid rrnaribosomal rna taftriethanolamine formalin xbf*any x. bovienii strain associated with any s. feltiae strain xbffl x. bovienii associated with s. feltiae fl xbfm x. bovienii associated with s. feltiae moldova xbfsn x. bovienii associated with s. feltiae sn xbi x. bovienii associated with s. intermedium type xbj x. bovienii associated with s. jollieti monsanto xbkbu x. bovienii associated with s. kraussei nemasys-l xbkca x. bovienii associated with s. kraussei quebec xbo x. bovienii associated with s. oregonense os-10 xbp x. bovienii associated with s. puntauvense li6 xn x. nematophila associated with s. carpocapsae all 
association with non-cognates resulted in decreased fitness, virulence, and bacterial carriage of the nematode-bacterial pairings. non-cognate nematode-bacteria pairings result in reductions for multiple measures of success, including total progeny production and virulence. entomopathogenic nematodes and their bacterial symbionts are a useful, tractable, and reliable model for testing hypotheses regarding the evolution, maintenance, persistence, and fate of mutualisms. entomopathogenic nematodessymbiosismutualismpartner choicenational sceince foundationios-0919565stock s. patricia national foundation for science ios-0920631graduate and professional student council grant n/anih iracda k-12n/apeterson brittany f. issue-copyright-statement© the author 2017 together steinernema nematodes and their bacterial symbionts successfully infect, kill, utilize, and exit their insect hosts. finally, we show that steinernema nematodes can distinguish heterospecific and some conspecific non-cognate symbionts in behavioral choice assays. steinernema-xenorhabdus symbioses are tightly governed by partner recognition and fidelity. additionally, nematode infective juveniles carry fewer bacterial cells when colonized by a non-cognate symbiont. this interaction can be complicated when two different nematodes co-infect an insect host.
 steinernematid nematodes form obligate symbioses with bacteria from the genus xenorhabdus. together steinernema nematodes and their bacterial symbionts successfully infect, kill, utilize, and exit their insect hosts. during this process the nematodes and bacteria disassociate requiring them to re-associate before emerging from the host. this interaction can be complicated when two different nematodes co-infect an insect host.  non-cognate nematode-bacteria pairings result in reductions for multiple measures of success, including total progeny production and virulence. additionally, nematode infective juveniles carry fewer bacterial cells when colonized by a non-cognate symbiont. finally, we show that steinernema nematodes can distinguish heterospecific and some conspecific non-cognate symbionts in behavioral choice assays.  steinernema-xenorhabdus symbioses are tightly governed by partner recognition and fidelity. association with non-cognates resulted in decreased fitness, virulence, and bacterial carriage of the nematode-bacterial pairings. entomopathogenic nematodes and their bacterial symbionts are a useful, tractable, and reliable model for testing hypotheses regarding the evolution, maintenance, persistence, and fate of mutualisms. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords entomopathogenic nematodessymbiosismutualismpartner choicenational sceince foundationios-0919565stock s. patricia national foundation for science ios-0920631graduate and professional student council grant n/anih iracda k-12n/apeterson brittany f. issue-copyright-statement© the author 2017 