BACKGROUND
recent advances have shown that the community of microbial species within a host has substantial and direct effects on host fitness and trophic interactions in the host community  <cit> .
when multiple species of microparasites infect the same host, the interactions among parasites can exacerbate the effects of the microbial communities on host fitness  <cit>  and select for mechanisms that alter the epidemiology of the parasites .
investigations into the pattern and mechanisms involved in the interactions among microbes within a host are important as the overwhelming majority of vertebrates in nature are infected by multiple microbial species simultaneously, often with multiple lineages of the same species  <cit> .
strains of the same species are likely to occupy the same ecological niche and are therefore predicted to have more intense competitive interactions affecting their respective fitness; in addition strains can also facilitate each other, maximizing overall infection and the transmission  <cit> .
here we use borrelia burgdorferi, the bacterium responsible for human lyme disease, as a model to investigate patterns and mechanisms of ecological interactions among microbial strains of the same species and to investigate the effects of these interactions on microbial fitness.
ixodid ticks, the genus that carries and transmits b. burgdorferi among vertebrate hosts, can acquire multiple b. burgdorferi strains by feeding on a single infected host  <cit> .
in geographic regions where b. burgdorferi is prevalent, individual white footed mice - peromyscus leucopus, the main host of the bacterium in the wild - are often infected with multiple strains of b. burgdorferi due to either the multiplicity of infection within individual feeding ticks  <cit>  or to the large number of infectious tick bites each individual vertebrate incurs over its lifetime  <cit> .
thus, multiple b. burgdorferi strains commonly infect individual vertebrate hosts and individual tick vectors, providing the opportunity for both direct competition  and indirect interactions .
strains of b. burgdorferi are transmitted to ticks feeding on peromyscus leucopus mice at much higher rates when the mouse is only infected with a single strain than when the mouse is infected by two or more strains  <cit> .
these results are consistent in both experimentally and naturally infected mice, indicating that the fitness of each strain is compromised by co-infections with other strains.
while these data suggest that b. burgdorferi strains interact competitively within a host, this hypothesis has not been experimentally investigated and the types of competitive interactions among strains within hosts are not known.
any of the three competition types  could be operating in this system.
for example, a primary strain could actively prevent the dissemination of a secondary strain in adequate sites, reducing overall transmission success of the secondary strains to feeding ticks.
alternatively, a strain could exploit and deplete a resource at the expense of others thanks to a more efficient exploitation or early segregation of the resource.
lastly, a strain could prime the host immune response to target an incoming strain in a process resembling apparent competition  <cit> .
these processes are not mutually exclusive and each mechanism could affect the competitive outcome among infecting strains.
in this study we looked for empirical evidence of competition among three coinfecting strains of b. burgdorferi, assessed the ecological patterns of the competition, and investigated a potential molecular mechanism affecting competition.
the wild bacteria b. burgdorferi is particularly variable at the genetic locus coding for the outer surface protein c   <cit> ; ospc is critical for the establishment of the bacteria in early stages of mammalian infection  <cit>  even though the exact mechanism of action is still debated  <cit>  and is known to trigger a strong immune response  <cit> .
some strains are more commonly found in certain host species than others  <cit> .
we used the ospc strains a, k, and n, which are infectious to humans  <cit> .
mice received two sequential inoculations to mimic the timing in which hosts can be exposed to b. burgdorferi-infected ticks in nature  <cit> .
the fully crossed-design created three groups of mice exposed to homologous infections , and six groups of mice exposed to heterologous infections .
we investigated the fitness of each strain as the rate of transmission of that strain to feeding ticks four times during the 90-day experiment .
additionally, we investigated if strains interfered with the colonization of disseminated sites within hosts by testing for the presence of each strain in four tissues of each host.
these data allow us to assess the competitive hierarchy among strains, the effect of a secondary infection on the fitness of a disseminated strain, and the ability of each strain to colonize a previously infected host.
to test if the observed competition was mediated by the adaptive immune system of the host, we also measured the immunoglobulin g produced against strain-specific antigen , against the flagellin protein  common to all b. burgdorferi strains, and against all antigens .
methods
borrelia burgdorferi and mice
we used three borrelia burgdorferi sensu stricto isolates that differed at their ospc locus to experimentally infect  <dig> three-month old female c3h/hej mice .
the b. burgdorferi isolates 97– <dig> , 97– <dig> , and  <dig>  were obtained through the center for disease control as part of the standard patient care, which did not require ethical approval.
each strain  was grown in bsk-ii complete medium  <cit>  at 35°c, washed three times with phosphate buffer saline  and concentrations were equalized after cell counting by hemocytometer, prior to subcutaneous injection of 5* <dig> cells in 200 μl of pbs.
mice were randomly assigned to a co-infection treatment  after acclimation to the animal facility for one week.
mice were randomly paired and kept two per cage.
two individuals in groups kk and nn died for unknown reasons during the 90-day course of the experiment.
mice in the homologous co-infection treatment groups were infected with same strain on days  <dig> and day  <dig>  of the experiment .
mice in the heterologous co-infection treatments groups were infected by two different strains during the experiment .
at day  <dig>  mice were sacrificed by co <dig> inhalation and blood was collected by cardiac puncture.
tissue samples from the ear, the bladder, the heart, and the mammary glands were extracted and frozen for subsequent analysis of tissue tropism of each strain .
the animal experimentation protocol was approved by the university of pennsylvania iacuc .figure  <dig>
infection in mouse tissues at day  <dig> was dominated by the primary strain.
few tissue samples tested positive for the secondary strains.
ear tissue was less frequently infected than the other organs.
n =  <dig> mice per group except in groups kk and nn where n =  <dig> mice.
infection status
the infection status of each mouse was determined using nested pcr of the ospc gene of mouse tissue samples and xenodiagnostic ticks.
xenodiagnosis, allowing uninfected larval ticks to feed on a host, is a traditional method of diagnosis for vector-borne diseases and can be used to estimate the transmission efficiency of b. burgdorferi from an infected mouse to naïve larval ticks.
xenodiagnostic larval ticks  were placed on each mouse on days  <dig>   <dig>   <dig>  and  <dig> and allowed to feed to repletion  <cit> .
engorged larval ticks were kept individually in microtubes containing a 1-cm <dig> square of moistened blotting paper until they molted to the nymphal stage.
dna from  <dig> whole nymphal ticks  or 25 μg of each mouse tissue was extracted using the dneasy blood and tissue kit .
the presence of b. burgdorferi in each sample was detected by nested pcr of the ospc gene using two sets of specifically designed primers that targeted conserved regions for all our strains .
the first pcr reaction amplified a 1087 bp fragment using primers oc- <dig>  and oc <dig> .
the second pcr reaction amplified a 657 bp fragment using primers oc <dig>  and oc <dig> .
the first pcr reaction contained 4 μl of the dna extraction,  <dig>  mm mgcl <dig>   <dig>  mm dntp,  <dig>  μm of each primer, 1x buffer, and  <dig> unit of recombinant taq  in a total volume of 50 μl.
after an initial denaturation at 95°c for 1 min, the mixture was run for 30 cycles at 95°c for 40 sec, 54°c for 35 sec, and 72°c for 60 sec.
the conditions for the second pcr reaction were the same as above except that 2 μl of product from the first pcr reaction was used as template, annealing temperature was 53°c, and the reaction was run for 40 cycles.
the resulting products were run on a 1% ethidium bromide agarose gel to assess the presence of b. burgdorferi.
all samples were tested by pcr two times to protect against false positive and false negative assignments.
samples with multiple strains
pcr products were digested with allele-specific restriction enzymes to identify the ospc alleles present in each sample.
the ospc alleles a, k, and n can be cleaved by restriction enzymes bgl ii , pst i, and sac i , respectively.
digestion products were run on 2% ethidium bromide agarose gel in parallel with control pcr products .
immune response
the total titer of circulating immunoglobulin-g  and the titer of igg specific for b. burgdorferi flagellin, ospc-a, ospc-k, or ospc-n was measured by elisa.
blood was collected by submandibular puncture from each mouse on days  <dig>   <dig>   <dig>  and  <dig>  centrifuged at 2500 rpm for 10 min, and the sera from each sample was collected and frozen at −20°c.
total igg titers were measured using the mouse total igg elisa kit ready-set-go!
and following the protocol of the manufacturer.
the titers of antibodies against individual b. burgdorferi proteins were measured by coating nunc maxisorbent 96-well plates with the b. burgdorferi protein overnight.
recombinant b. burgdorferi ospc proteins were produced as described in  <cit>  and were diluted at 10 ng/μl in pbs.
recombinant flagellin was purchased  and diluted to 2 ng/μl in 5 m urea.
after incubation, each well was blocked with 2% bsa for two hours and incubated with serum diluted 1: <dig> in 2% bsa for each ospc protein or 1: <dig> in pbs for flagellin.
each well was then incubated for 45 min with the secondary antibody  diluted at 1:5000).
all wells were washed three times with pbs-tween  <dig> % between each step.
the color reagent p-nitrophenyl phosphate  was added and the kinetic phosphatase reaction was monitored at 405 nm for 30 minutes .
general statistical methods
differences in the infection of mouse tissues and the mouse-to-tick transmission were modeled using generalized linear mixed effects models  with binomial errors while differences in antibody titers were modeled using glmms with normal error terms.
the models used to test the effect of the order of inoculation or strain identity on the response variables included strain identity , the order of inoculation , and the date of xenodiagnoses  as fixed factors, all possible 2-way and 3-way interaction terms, and mouse identity  as a random factor.
the effect of the order of inoculation on mouse-to-tick transmission during the last two transmission events when both strains were present in the mouse  also included ticks nested within mouse as a random factor.
the statistical significance of fixed factors and interactions among fixed factors for each model were determined using likelihood ratio tests.
additionally, the most parsimonious models were selected using the aic model selection approach.
analyses were done using the lme <dig> package  <cit>  in r .
final models are exposed in additional file  <dig>
RESULTS
primary strain infection and transmission
evidence of infection by the primary strain was present in both the tissue samples and the xenodiagnostic ticks  from all mice.
after sacrifice, the primary strain was found in at least one tissue from all mice although most mice showed evidence of the primary strain in multiple tissues .
evidence of active infection was least prevalent in ear tissue  but common in the heart, bladder and mammary gland .
the identity of the strain affected the probability of detecting infection in tissues  as strain a was more commonly detected than strains k or n .
the primary strain was also frequently detected in xenodiagnostic ticks that had taken their larval blood meal from the experimental mice  .
of the  <dig> ticks that tested positive for b. burgdorferi, the primary strain was present in  <dig>  the identity of the strain affected the rate at which the strain was transmitted to feeding ticks , with strain k being transmitted at the highest rate followed by strain a .
the rate of transmission of the primary strain to feeding ticks varied among the four xenodiagnostic time points , although the pattern of temporal variation differed significantly among strains .
while strain k was transmitted to feeding ticks at a consistently high rate throughout the experiment, the transmission rate of strain a decreased over time, and the transmission rate of strain n peaked in the middle time points.figure  <dig>
average mouse-to-tick transmission rates  of the primary strain.
strain a - blue circles; strain k - red squares; strain n - green triangles.
the mouse-to-tick transmission rate of the secondary strain was pooled for the three strain types .
inf  <dig> - day of inoculation with primary strain; inf  <dig> - day of inoculation with secondary strain.
d <dig>  d <dig>  d <dig>  d <dig> refer to the days of blood sampling and xenodiagnosis.
sac = sacrifice for biopsy of organs .
secondary strain infection and transmission
the secondary strain – the second strain the mice were exposed to – was detected in only seven of the  <dig> mice that were challenged with two different strains.
of these seven mice, only three had evidence of infection with the secondary strain in tissue samples; strain a was found in the ear, heart, and mammary tissues of one mouse and in the bladder of another mouse, and strain n was found in the heart of a kn mouse .
this mouse, like four others, showed evidence of secondary strain infection only by one to two positive xenodiagnostic ticks.
strain k was transmitted from three mice, strain a from one mouse, and strain n from one mouse.
neither strain identity, time of xenodiagnosis, nor any interaction term significantly affected the mouse-to-tick transmission rate of the secondary strain.
interestingly, the identity of the secondary strain did not influence the mouse-to-tick transmission of the primary strain  .
since only seven ticks carried the secondary strain , we lacked power to statistically test the effect of the primary strain on the transmission of the secondary strain, but we observed transmission of the secondary strain with each of our primary strains.
for mice inoculated with multiple b. burgdorferi strains, mouse-to-tick transmission of the primary strain was an order of magnitude higher than that of the secondary strain during the last two transmission events when both strains were present in the mouse .
thus, the order of inoculation of the strains had a significant effect on the probability of transmission to mice .
the log-odds ratio of the contrast between the secondary and the primary strain was − <dig>  .
that is, the probability of transmission to feeding larvae was 31-fold greater when a strain is inoculated first than when it was inoculated second.table  <dig>
mouse-to-tick transmission was more successful for the primary strain than the secondary strain
d46
d65
strain
primary
secondary
primary
secondary
transmission success was calculated as the proportion of ticks infected with strain  <dig>  relative to the total number of infected ticks.
the number of ticks infected with each strain are in brackets.
immune response
the primary strain induced a strong immune response as measured by total immunoglobulin g  production as well as a strong specific immune response against flagellin and against the ospc genotype expressed by the primary strain .
the anti-ospc response induced against the primary strain had low cross-reactivity with the ospc antigen of the secondary strains.
in other words, the anti-ospc antibodies induced by strain a infection had low affinity for ospc type k nor type n proteins and vice versa.
the inoculation of the secondary strain did not affect the total igg level, the anti-fla response, nor the specific anti-ospc response targeting either the primary or the secondary strain, even in the mice with conclusive evidence of infection by the secondary strain.figure  <dig>
the primary strain induced a strong strain-specific adaptive immune response.
the secondary strain did not alter immune response.
each panel shows one of nine co-infection treatments.
antibody profiles are similar within columns, which share the primary strain, but differ among rows, which share the secondary strain.
this pattern indicates that the primary strain dominated the antibody response and that antibody profiles differed among primary strains .
arrows indicate the timing of the primary and secondary inoculation.
shown are the mean antibody titers  for total igg , anti-flagellin igg , anti-ospca igg , anti-ospck , and anti-ospcn .
to facilitate viewing, absolute differences in scale among the five antibody variables were removed by standardizing them to z-scores .
we used principal component analysis  to characterize the antibody profiles of the mice as the five antibody response variables  were strongly and positively correlated .
the first principal component  captured  <dig> % of the total variance in the data  and represented the overall igg response as the five antibody variables were all positively correlated with this component.
the second principal component  explained  <dig> % of the total variance  and represented a contrast between high levels of total igg and ospck-igg and low levels of fla-igg and ospca-igg.
thus pc <dig> can be considered as a measure of specificity because individuals that have a strong anti-ospca response have a weak anti-ospck response and vice versa.
the ability of our pca to capture the relevant variation in the antibody profiles is shown by the scatterplot of pc <dig> versus pc <dig>  which causes mice to cluster according to the primary strain and not according to the secondary strain .
the overall igg response depended on the primary strain , and the day of blood sampling , but was independent of the secondary strain .
similarly, only the primary strain affected the ospc-specific immune response  while the day of blood sampling and the secondary strain had no effect .
discussion
inhibitory priority effect
interactions among co-infecting microbes in a host can alter the within-host microbial community and ultimately influence the diversity and population dynamics of microbes in nature.
we used experimental infections to determine the competitive interactions among strains of borrelia burgdorferi, the pathogenic bacterium that causes human lyme disease.
we detected the presence of strong, inhibitory competitive interactions among strains of b. burgdorferi within a host.
the order in which strains were inoculated into the mouse was the dominant factor in the competitive outcome among strains.
in all cases, the first strain to infect had a strong fitness advantage primarily because the colonization of mouse tissues by subsequent strains was inhibited.
although the strains differed in their rates of transmission in the absence of competition, indicating a potential competitive hierarchy, all strains were excluded when secondary in infection order regardless of the primary and secondary strain identity, suggesting the absence of a competitive hierarchy in sequential infections.
the primary strain evoked a strong and specific adaptive immune response, whereas there was no change in antibody level following the secondary inoculation contrary to expectation of activation of a memory immune response.
overall, the data do not support a major role of apparent competition due to the adaptive immune response mediated by igg in the observed priority effect.
evidence of infection by the secondary strain was rare  suggesting that the primary infection either reduced the probability of infection by the secondary strain or limited the population size of the secondary strain within the host.
the data presented are consistent with the latter hypothesis as strains with small population sizes are difficult to detect in host tissues by pcr and would rarely be transmitted to feeding ticks  <cit> .
indirect interactions among strains mediated by the host immune system were not readily apparent from the data.
all mice mounted a vigorous immune response against the initial strain, which could have been protective against subsequent b. burgdorferi infections  <cit> .
however, the hypothesis that initial infections reduce the probability of infection by subsequent strains, similar to expectations from vaccination, is not congruous.
it seems unlikely that previously induced antibodies against other borrelia specific antigens, such as the antigenic proteins outer surface protein a  or ospc, can explain the observed priority effect data for four reasons.
first, the exclusion of the secondary strain occurred without increasing antibody titers following the inoculation of the secondary strain as would be expected if memory antibodies were involved.
if memory antibodies were activated, they would trigger the proliferation of b cells that would then produce large amount of antibodies  <cit> .
in contrast, the constant level of antibodies in our data suggests an immunosuppressive effect by the primary strain  <cit> .
second, at least 29% of mice were infected by the secondary strain but at very low levels.
vaccination effect would result in less mice infected with the second strain and would affect the transmission of the bacteria to ticks of all strains independently of the order of infection, which is not what we observed in mice successfully infected by the secondary strain  <cit> .
third, bacteria in culture express different phenotypes, some expressing ospc, others expressing other proteins  <cit> .
if the rejection of the secondary strain was due to a specific antigen recognition, a proportion of the inoculated bacteria would be immune to the specific antibodies and would be able to disperse.
lastly, it is unlikely that the titer of antibodies from the first immune response was at saturation as mice primed with the n strain had significantly lower antibody titers than those initially infected by strains a or k yet still showed no increase when challenged by a second strain .
thus, the mechanism leading to the observed priority effect pattern is likely caused by direct or indirect exclusion  of the second strain from colonizing disseminated sites in the host.
future investigations using serial sampling of disseminated sites are necessary to observe direct interactions among strains.
ecological consequences of priority effect
the strength of the observed priority effect could, repeated over several years in a homogenous environment, reduce the diversity of b. burgdorferi strains in nature and rapidly lead to a uniform population.
however, the diversity of b. burgdorferi strains remains high in many geographic areas and even within wild mice  <cit> .
high strain diversity within wild mice has been observed in geographic regions where both tick vectors and b. burgdorferi are highly prevalent and mice are fed upon by tens of infected ticks, many carrying multiple strains, within a period of 30 days  <cit> .
thus, the wild mice in these studies were challenged by multiple strains nearly simultaneously, which would exclude the possibility of detecting a priority effect.
competitive exclusion due to the priority effect mechanism may have an impact in geographic areas where b. burgdorferi and ticks are present but at low prevalence, such as areas where lyme disease is spreading  <cit> .
in these areas, the time between infectious tick bites can be long resulting in competitive exclusion of the latter strains and dramatic fitness benefits to the most prevalent strains.
this effect has the potential to inhibit or prevent the invasion of subsequent strains into a geographic region.
additionally, in many areas the season in which ticks are active in extends over several months such that strains in late-season feeding ticks that take their blood meal at the end of the season will have lower success than those from early-season feeding ticks, even in highly endemic regions  <cit> .
this scenario suggests that bacteria face strong selective pressures for the emergence of mechanisms to secure their transmission such as vector behavior manipulation for early or synchronous host seeking.
CONCLUSIONS
the strong inhibitory priority effect is the dominant mechanism controlling the competitive outcome among b. burgdorferi strains in a sequentially challenged vertebrate host.
this inhibitory priority effect occurs regardless of the fitness differences among strains in the absence of competition.
although the molecular mechanism leading to the observed priority effect could not be precisely determined with the current experiments, it is effective within 10 days of the secondary inoculation and results in the inhibition of dissemination or colonization of the host.
the priority effect has the potential to limit or reduce the diversity of strains in nature, especially in low prevalence areas where sequential challenges are likely.
expansion of b. burgdorferi into new geographic areas or changes in tick questing behaviors in response to climate change may increase the impact of the priority effect on competitive interactions among b. burgdorferi strains and cause the most common b. burgdorferi lineages, which are also the most infectious in humans  <cit> , to become even more prevalent.
additional file
additional file 1: figure s <dig>  chart summarizing the diagnostic of infection by first and secondary strain for each mouse.
figure s <dig>  the immune response depends on the primary strain.
table s <dig>  organs were infected by the primary inoculated strain, though strains did not infect at the same rate and some organs were more susceptible to infection, without interactions.
table s <dig>  mouse-to-tick transmission of primary strain depended on the identity of the primary strain, the day of xenodiagnoses , and the strain:day interaction.
table s <dig>  among heterologous mice at d <dig> and d <dig>  mouse-to-tick transmission depended on the order of inoculation , strain identity , the day of xenodiagnoses , and the strain:day interaction.
table s <dig>  the titers of igg were correlated.
table s <dig>  the results from the principal component analysis .
table s <dig>  the strength of the immune response, summarized by pc <dig>  depended on the day of blood sampling and the primary strain, but was independent of the secondary strain.
table s <dig>  the specific adaptive immune response, summarized by pc <dig>  depended on the primary inoculated strain and the date considered, regardless of the secondary strain.
competing interests
the authors declare that they have no competing interests.
authors’ contributions
gd and db designed the experiment and wrote the manuscript.
gd conducted the experiment with the technical assistance of td, cg and sm.
gd conducted the data analysis.
all authors contributed to the final manuscript.
all authors read and approved the final manuscript.
authors’ information
godefroy devevey was a snsf post-doctoral associate in the brisson lab.
he is currently working with amy pedersen at the university of edinburgh, uk.
he is interested in parasite communities, the nature and consequences of parasite-parasite interactions, and the differential roles of hosts to generalist parasites.
christopher j. graves contributed to this work as a research technician in the brisson lab, and is currently a phd candidate in the department of ecology and evolutionary biology at brown university.
trang dang was a research technician in the brisson lab.
she is currently a graduate student with dr. andrew blaustein at oregon state university, studying disease ecology in amphibians with special interests in pathogen diversity and virulence.
sarah murray is an undergraduate studying physics at the university of pennsylvania.
she distinguished by an excellent contribution to this work during her junior year.
dustin brisson is the head of the evolution and ecology of disease systems lab at the university of pennsylvania.
many thanks to maarten j. voordouw for his important help with the statistical analysis and to andrea graham for constructive discussion on the results.
thank you to the two reviewers of peerage of science who largely helped to improve this manuscript.
the study was supported by grant ai <dig> from the national institute of health and grant ck <dig> from the center for diseases control and prevention.
gd was partially funded by the swiss science foundation fellowship for prospective researcher pblap3-127724/ <dig>

 recent advances have shown that the community of microbial species within a host has substantial and direct effects on host fitness and trophic interactions in the host community. when multiple species of microparasites infect the same host, the interactions among parasites can exacerbate the effects of the microbial communities on host fitness and select for mechanisms that alter the epidemiology of the parasites. investigations into the pattern and mechanisms involved in the interactions among microbes within a host are important as the overwhelming majority of vertebrates in nature are infected by multiple microbial species simultaneously, often with multiple lineages of the same species. strains of the same species are likely to occupy the same ecological niche and are therefore predicted to have more intense competitive interactions affecting their respective fitness; in addition strains can also facilitate each other, maximizing overall infection and the transmission. here we use borrelia burgdorferi, the bacterium responsible for human lyme disease, as a model to investigate patterns and mechanisms of ecological interactions among microbial strains of the same species and to investigate the effects of these interactions on microbial fitness. ixodid ticks, the genus that carries and transmits b. burgdorferi among vertebrate hosts, can acquire multiple b. burgdorferi strains by feeding on a single infected host. in geographic regions where b. burgdorferi is prevalent, individual white footed mice - peromyscus leucopus, the main host of the bacterium in the wild - are often infected with multiple strains of b. burgdorferi due to either the multiplicity of infection within individual feeding ticks or to the large number of infectious tick bites each individual vertebrate incurs over its lifetime. thus, multiple b. burgdorferi strains commonly infect individual vertebrate hosts and individual tick vectors, providing the opportunity for both direct competition and indirect interactions. strains of b. burgdorferi are transmitted to ticks feeding on peromyscus leucopus mice at much higher rates when the mouse is only infected with a single strain than when the mouse is infected by two or more strains. these results are consistent in both experimentally and naturally infected mice, indicating that the fitness of each strain is compromised by co-infections with other strains. while these data suggest that b. burgdorferi strains interact competitively within a host, this hypothesis has not been experimentally investigated and the types of competitive interactions among strains within hosts are not known. any of the three competition types could be operating in this system. for example, a primary strain could actively prevent the dissemination of a secondary strain in adequate sites, reducing overall transmission success of the secondary strains to feeding ticks. alternatively, a strain could exploit and deplete a resource at the expense of others thanks to a more efficient exploitation or early segregation of the resource. lastly, a strain could prime the host immune response to target an incoming strain in a process resembling apparent competition. these processes are not mutually exclusive and each mechanism could affect the competitive outcome among infecting strains. in this study we looked for empirical evidence of competition among three coinfecting strains of b. burgdorferi, assessed the ecological patterns of the competition, and investigated a potential molecular mechanism affecting competition. the wild bacteria b. burgdorferi is particularly variable at the genetic locus coding for the outer surface protein c; ospc is critical for the establishment of the bacteria in early stages of mammalian infection even though the exact mechanism of action is still debated and is known to trigger a strong immune response. some strains are more commonly found in certain host species than others. we used the ospc strains a, k, and n, which are infectious to humans. mice received two sequential inoculations to mimic the timing in which hosts can be exposed to b. burgdorferi-infected ticks in nature. the fully crossed-design created three groups of mice exposed to homologous infections, and six groups of mice exposed to heterologous infections. we investigated the fitness of each strain as the rate of transmission of that strain to feeding ticks four times during the 90-day experiment. additionally, we investigated if strains interfered with the colonization of disseminated sites within hosts by testing for the presence of each strain in four tissues of each host. these data allow us to assess the competitive hierarchy among strains, the effect of a secondary infection on the fitness of a disseminated strain, and the ability of each strain to colonize a previously infected host. to test if the observed competition was mediated by the adaptive immune system of the host, we also measured the immunoglobulin g produced against strain-specific antigen, against the flagellin protein common to all b. burgdorferi strains, and against all antigens. methods borrelia burgdorferi and mice we used three borrelia burgdorferi sensu stricto isolates that differed at their ospc locus to experimentally infect three-month old female c3h/hej mice. the b. burgdorferi isolates 97–, 97–, and were obtained through the center for disease control as part of the standard patient care, which did not require ethical approval. each strain was grown in bsk-ii complete medium at 35°c, washed three times with phosphate buffer saline and concentrations were equalized after cell counting by hemocytometer, prior to subcutaneous injection of 5* cells in 200 μl of pbs. mice were randomly assigned to a co-infection treatment after acclimation to the animal facility for one week. mice were randomly paired and kept two per cage. two individuals in groups kk and nn died for unknown reasons during the 90-day course of the experiment. mice in the homologous co-infection treatment groups were infected with same strain on days and day of the experiment. mice in the heterologous co-infection treatments groups were infected by two different strains during the experiment. at day mice were sacrificed by co inhalation and blood was collected by cardiac puncture. tissue samples from the ear, the bladder, the heart, and the mammary glands were extracted and frozen for subsequent analysis of tissue tropism of each strain. the animal experimentation protocol was approved by the university of pennsylvania iacuc.figure infection in mouse tissues at day was dominated by the primary strain. few tissue samples tested positive for the secondary strains. ear tissue was less frequently infected than the other organs. n =  mice per group except in groups kk and nn where n =  mice. infection status the infection status of each mouse was determined using nested pcr of the ospc gene of mouse tissue samples and xenodiagnostic ticks. xenodiagnosis, allowing uninfected larval ticks to feed on a host, is a traditional method of diagnosis for vector-borne diseases and can be used to estimate the transmission efficiency of b. burgdorferi from an infected mouse to naïve larval ticks. xenodiagnostic larval ticks were placed on each mouse on days and and allowed to feed to repletion. engorged larval ticks were kept individually in microtubes containing a 1-cm square of moistened blotting paper until they molted to the nymphal stage. dna from whole nymphal ticks or 25 μg of each mouse tissue was extracted using the dneasy blood and tissue kit. the presence of b. burgdorferi in each sample was detected by nested pcr of the ospc gene using two sets of specifically designed primers that targeted conserved regions for all our strains. the first pcr reaction amplified a 1087 bp fragment using primers oc- and oc. the second pcr reaction amplified a 657 bp fragment using primers oc and oc. the first pcr reaction contained 4 μl of the dna extraction,  mm mgcl  mm dntp,  μm of each primer, 1x buffer, and unit of recombinant taq in a total volume of 50 μl. after an initial denaturation at 95°c for 1 min, the mixture was run for 30 cycles at 95°c for 40 sec, 54°c for 35 sec, and 72°c for 60 sec. the conditions for the second pcr reaction were the same as above except that 2 μl of product from the first pcr reaction was used as template, annealing temperature was 53°c, and the reaction was run for 40 cycles. the resulting products were run on a 1% ethidium bromide agarose gel to assess the presence of b. burgdorferi. all samples were tested by pcr two times to protect against false positive and false negative assignments. samples with multiple strains pcr products were digested with allele-specific restriction enzymes to identify the ospc alleles present in each sample. the ospc alleles a, k, and n can be cleaved by restriction enzymes bgl ii, pst i, and sac i, respectively. digestion products were run on 2% ethidium bromide agarose gel in parallel with control pcr products. immune response the total titer of circulating immunoglobulin-g and the titer of igg specific for b. burgdorferi flagellin, ospc-a, ospc-k, or ospc-n was measured by elisa. blood was collected by submandibular puncture from each mouse on days and centrifuged at 2500 rpm for 10 min, and the sera from each sample was collected and frozen at −20°c. total igg titers were measured using the mouse total igg elisa kit ready-set-go! and following the protocol of the manufacturer. the titers of antibodies against individual b. burgdorferi proteins were measured by coating nunc maxisorbent 96-well plates with the b. burgdorferi protein overnight. recombinant b. burgdorferi ospc proteins were produced as described in and were diluted at 10 ng/μl in pbs. recombinant flagellin was purchased and diluted to 2 ng/μl in 5 m urea. after incubation, each well was blocked with 2% bsa for two hours and incubated with serum diluted 1: in 2% bsa for each ospc protein or 1: in pbs for flagellin. each well was then incubated for 45 min with the secondary antibody diluted at 1:5000). all wells were washed three times with pbs-tween % between each step. the color reagent p-nitrophenyl phosphate was added and the kinetic phosphatase reaction was monitored at 405 nm for 30 minutes. general statistical methods differences in the infection of mouse tissues and the mouse-to-tick transmission were modeled using generalized linear mixed effects models with binomial errors while differences in antibody titers were modeled using glmms with normal error terms. the models used to test the effect of the order of inoculation or strain identity on the response variables included strain identity, the order of inoculation, and the date of xenodiagnoses as fixed factors, all possible 2-way and 3-way interaction terms, and mouse identity as a random factor. the effect of the order of inoculation on mouse-to-tick transmission during the last two transmission events when both strains were present in the mouse also included ticks nested within mouse as a random factor. the statistical significance of fixed factors and interactions among fixed factors for each model were determined using likelihood ratio tests. additionally, the most parsimonious models were selected using the aic model selection approach. analyses were done using the lme package in r. final models are exposed in additional file  primary strain infection and transmission evidence of infection by the primary strain was present in both the tissue samples and the xenodiagnostic ticks from all mice. after sacrifice, the primary strain was found in at least one tissue from all mice although most mice showed evidence of the primary strain in multiple tissues. evidence of active infection was least prevalent in ear tissue but common in the heart, bladder and mammary gland. the identity of the strain affected the probability of detecting infection in tissues as strain a was more commonly detected than strains k or n. the primary strain was also frequently detected in xenodiagnostic ticks that had taken their larval blood meal from the experimental mice. of the ticks that tested positive for b. burgdorferi, the primary strain was present in the identity of the strain affected the rate at which the strain was transmitted to feeding ticks, with strain k being transmitted at the highest rate followed by strain a. the rate of transmission of the primary strain to feeding ticks varied among the four xenodiagnostic time points, although the pattern of temporal variation differed significantly among strains. while strain k was transmitted to feeding ticks at a consistently high rate throughout the experiment, the transmission rate of strain a decreased over time, and the transmission rate of strain n peaked in the middle time points.figure average mouse-to-tick transmission rates of the primary strain. strain a - blue circles; strain k - red squares; strain n - green triangles. the mouse-to-tick transmission rate of the secondary strain was pooled for the three strain types. inf - day of inoculation with primary strain; inf - day of inoculation with secondary strain. d d d d refer to the days of blood sampling and xenodiagnosis. sac = sacrifice for biopsy of organs. secondary strain infection and transmission the secondary strain – the second strain the mice were exposed to – was detected in only seven of the mice that were challenged with two different strains. of these seven mice, only three had evidence of infection with the secondary strain in tissue samples; strain a was found in the ear, heart, and mammary tissues of one mouse and in the bladder of another mouse, and strain n was found in the heart of a kn mouse. this mouse, like four others, showed evidence of secondary strain infection only by one to two positive xenodiagnostic ticks. strain k was transmitted from three mice, strain a from one mouse, and strain n from one mouse. neither strain identity, time of xenodiagnosis, nor any interaction term significantly affected the mouse-to-tick transmission rate of the secondary strain. interestingly, the identity of the secondary strain did not influence the mouse-to-tick transmission of the primary strain. since only seven ticks carried the secondary strain, we lacked power to statistically test the effect of the primary strain on the transmission of the secondary strain, but we observed transmission of the secondary strain with each of our primary strains. for mice inoculated with multiple b. burgdorferi strains, mouse-to-tick transmission of the primary strain was an order of magnitude higher than that of the secondary strain during the last two transmission events when both strains were present in the mouse. thus, the order of inoculation of the strains had a significant effect on the probability of transmission to mice. the log-odds ratio of the contrast between the secondary and the primary strain was −. that is, the probability of transmission to feeding larvae was 31-fold greater when a strain is inoculated first than when it was inoculated second.table mouse-to-tick transmission was more successful for the primary strain than the secondary strain d46 d65 strain primary secondary primary secondary transmission success was calculated as the proportion of ticks infected with strain relative to the total number of infected ticks. the number of ticks infected with each strain are in brackets. immune response the primary strain induced a strong immune response as measured by total immunoglobulin g production as well as a strong specific immune response against flagellin and against the ospc genotype expressed by the primary strain. the anti-ospc response induced against the primary strain had low cross-reactivity with the ospc antigen of the secondary strains. in other words, the anti-ospc antibodies induced by strain a infection had low affinity for ospc type k nor type n proteins and vice versa. the inoculation of the secondary strain did not affect the total igg level, the anti-fla response, nor the specific anti-ospc response targeting either the primary or the secondary strain, even in the mice with conclusive evidence of infection by the secondary strain.figure the primary strain induced a strong strain-specific adaptive immune response. the secondary strain did not alter immune response. each panel shows one of nine co-infection treatments. antibody profiles are similar within columns, which share the primary strain, but differ among rows, which share the secondary strain. this pattern indicates that the primary strain dominated the antibody response and that antibody profiles differed among primary strains. arrows indicate the timing of the primary and secondary inoculation. shown are the mean antibody titers for total igg, anti-flagellin igg, anti-ospca igg, anti-ospck, and anti-ospcn. to facilitate viewing, absolute differences in scale among the five antibody variables were removed by standardizing them to z-scores. we used principal component analysis to characterize the antibody profiles of the mice as the five antibody response variables were strongly and positively correlated. the first principal component captured % of the total variance in the data and represented the overall igg response as the five antibody variables were all positively correlated with this component. the second principal component explained % of the total variance and represented a contrast between high levels of total igg and ospck-igg and low levels of fla-igg and ospca-igg. thus pc can be considered as a measure of specificity because individuals that have a strong anti-ospca response have a weak anti-ospck response and vice versa. the ability of our pca to capture the relevant variation in the antibody profiles is shown by the scatterplot of pc versus pc which causes mice to cluster according to the primary strain and not according to the secondary strain. the overall igg response depended on the primary strain, and the day of blood sampling, but was independent of the secondary strain. similarly, only the primary strain affected the ospc-specific immune response while the day of blood sampling and the secondary strain had no effect. discussion inhibitory priority effect interactions among co-infecting microbes in a host can alter the within-host microbial community and ultimately influence the diversity and population dynamics of microbes in nature. we used experimental infections to determine the competitive interactions among strains of borrelia burgdorferi, the pathogenic bacterium that causes human lyme disease. we detected the presence of strong, inhibitory competitive interactions among strains of b. burgdorferi within a host. the order in which strains were inoculated into the mouse was the dominant factor in the competitive outcome among strains. in all cases, the first strain to infect had a strong fitness advantage primarily because the colonization of mouse tissues by subsequent strains was inhibited. although the strains differed in their rates of transmission in the absence of competition, indicating a potential competitive hierarchy, all strains were excluded when secondary in infection order regardless of the primary and secondary strain identity, suggesting the absence of a competitive hierarchy in sequential infections. the primary strain evoked a strong and specific adaptive immune response, whereas there was no change in antibody level following the secondary inoculation contrary to expectation of activation of a memory immune response. overall, the data do not support a major role of apparent competition due to the adaptive immune response mediated by igg in the observed priority effect. evidence of infection by the secondary strain was rare suggesting that the primary infection either reduced the probability of infection by the secondary strain or limited the population size of the secondary strain within the host. the data presented are consistent with the latter hypothesis as strains with small population sizes are difficult to detect in host tissues by pcr and would rarely be transmitted to feeding ticks. indirect interactions among strains mediated by the host immune system were not readily apparent from the data. all mice mounted a vigorous immune response against the initial strain, which could have been protective against subsequent b. burgdorferi infections. however, the hypothesis that initial infections reduce the probability of infection by subsequent strains, similar to expectations from vaccination, is not congruous. it seems unlikely that previously induced antibodies against other borrelia specific antigens, such as the antigenic proteins outer surface protein a or ospc, can explain the observed priority effect data for four reasons. first, the exclusion of the secondary strain occurred without increasing antibody titers following the inoculation of the secondary strain as would be expected if memory antibodies were involved. if memory antibodies were activated, they would trigger the proliferation of b cells that would then produce large amount of antibodies. in contrast, the constant level of antibodies in our data suggests an immunosuppressive effect by the primary strain. second, at least 29% of mice were infected by the secondary strain but at very low levels. vaccination effect would result in less mice infected with the second strain and would affect the transmission of the bacteria to ticks of all strains independently of the order of infection, which is not what we observed in mice successfully infected by the secondary strain. third, bacteria in culture express different phenotypes, some expressing ospc, others expressing other proteins. if the rejection of the secondary strain was due to a specific antigen recognition, a proportion of the inoculated bacteria would be immune to the specific antibodies and would be able to disperse. lastly, it is unlikely that the titer of antibodies from the first immune response was at saturation as mice primed with the n strain had significantly lower antibody titers than those initially infected by strains a or k yet still showed no increase when challenged by a second strain. thus, the mechanism leading to the observed priority effect pattern is likely caused by direct or indirect exclusion of the second strain from colonizing disseminated sites in the host. future investigations using serial sampling of disseminated sites are necessary to observe direct interactions among strains. ecological consequences of priority effect the strength of the observed priority effect could, repeated over several years in a homogenous environment, reduce the diversity of b. burgdorferi strains in nature and rapidly lead to a uniform population. however, the diversity of b. burgdorferi strains remains high in many geographic areas and even within wild mice. high strain diversity within wild mice has been observed in geographic regions where both tick vectors and b. burgdorferi are highly prevalent and mice are fed upon by tens of infected ticks, many carrying multiple strains, within a period of 30 days. thus, the wild mice in these studies were challenged by multiple strains nearly simultaneously, which would exclude the possibility of detecting a priority effect. competitive exclusion due to the priority effect mechanism may have an impact in geographic areas where b. burgdorferi and ticks are present but at low prevalence, such as areas where lyme disease is spreading. in these areas, the time between infectious tick bites can be long resulting in competitive exclusion of the latter strains and dramatic fitness benefits to the most prevalent strains. this effect has the potential to inhibit or prevent the invasion of subsequent strains into a geographic region. additionally, in many areas the season in which ticks are active in extends over several months such that strains in late-season feeding ticks that take their blood meal at the end of the season will have lower success than those from early-season feeding ticks, even in highly endemic regions. this scenario suggests that bacteria face strong selective pressures for the emergence of mechanisms to secure their transmission such as vector behavior manipulation for early or synchronous host seeking.  the strong inhibitory priority effect is the dominant mechanism controlling the competitive outcome among b. burgdorferi strains in a sequentially challenged vertebrate host. this inhibitory priority effect occurs regardless of the fitness differences among strains in the absence of competition. although the molecular mechanism leading to the observed priority effect could not be precisely determined with the current experiments, it is effective within 10 days of the secondary inoculation and results in the inhibition of dissemination or colonization of the host. the priority effect has the potential to limit or reduce the diversity of strains in nature, especially in low prevalence areas where sequential challenges are likely. expansion of b. burgdorferi into new geographic areas or changes in tick questing behaviors in response to climate change may increase the impact of the priority effect on competitive interactions among b. burgdorferi strains and cause the most common b. burgdorferi lineages, which are also the most infectious in humans, to become even more prevalent. additional file additional file 1: figure s chart summarizing the diagnostic of infection by first and secondary strain for each mouse. figure s the immune response depends on the primary strain. table s organs were infected by the primary inoculated strain, though strains did not infect at the same rate and some organs were more susceptible to infection, without interactions. table s mouse-to-tick transmission of primary strain depended on the identity of the primary strain, the day of xenodiagnoses, and the strain:day interaction. table s among heterologous mice at d and d mouse-to-tick transmission depended on the order of inoculation, strain identity, the day of xenodiagnoses, and the strain:day interaction. table s the titers of igg were correlated. table s the results from the principal component analysis. table s the strength of the immune response, summarized by pc depended on the day of blood sampling and the primary strain, but was independent of the secondary strain. table s the specific adaptive immune response, summarized by pc depended on the primary inoculated strain and the date considered, regardless of the secondary strain. competing interests the authors declare that they have no competing interests. authors’ contributions gd and db designed the experiment and wrote the manuscript. gd conducted the experiment with the technical assistance of td, cg and sm. gd conducted the data analysis. all authors contributed to the final manuscript. all authors read and approved the final manuscript. authors’ information godefroy devevey was a snsf post-doctoral associate in the brisson lab. he is currently working with amy pedersen at the university of edinburgh, uk. he is interested in parasite communities, the nature and consequences of parasite-parasite interactions, and the differential roles of hosts to generalist parasites. christopher j. graves contributed to this work as a research technician in the brisson lab, and is currently a phd candidate in the department of ecology and evolutionary biology at brown university. trang dang was a research technician in the brisson lab. she is currently a graduate student with dr. andrew blaustein at oregon state university, studying disease ecology in amphibians with special interests in pathogen diversity and virulence. sarah murray is an undergraduate studying physics at the university of pennsylvania. she distinguished by an excellent contribution to this work during her junior year. dustin brisson is the head of the evolution and ecology of disease systems lab at the university of pennsylvania. many thanks to maarten j. voordouw for his important help with the statistical analysis and to andrea graham for constructive discussion on the results. thank you to the two reviewers of peerage of science who largely helped to improve this manuscript. the study was supported by grant ai from the national institute of health and grant ck from the center for diseases control and prevention. gd was partially funded by the swiss science foundation fellowship for prospective researcher pblap3-127724/ 
the strong inhibitory priority effect is a dominant mechanism underlying competition for transmission between coinfecting b. burgdorferi strains, most likely through resource exploitation. we investigated the within-host interactions among strains of borrelia burgdorferi, the causative agent of lyme disease, using experimental infections in mice. we used data from these experimental infections to assess the effect of coinfection on bacterial dissemination and fitness as well as the effect of coinfection on host immune response compared to single infection. we used a fully crossed-design with three distinct strains, each group of hosts receiving two sequential inoculations. contrary to expectation, the strong and specific adaptive immune response evoked against the primary strain was not followed by production of immunoglobulins after the inoculation of the secondary strain, neither against strain-specific antigen nor against antigens common to all strains. in all cases, the primary strain a mouse was infected with had an absolute fitness advantage primarily since it was transmitted an order of magnitude more than the secondary strain. the mechanism of exclusion of the secondary strain is an inhibition of the colonization of mouse tissues, even though 29% of mice showed some evidence of infection by secondary strain. hence, the data do not support a major role of the immune response in the observed priority effect. the microbial community associated with a host is indeed influenced by a complex network of direct and indirect interactions between the host and the lineages of microbes it harbors, but the mechanisms are rarely established. the infection and transmission data strongly indicate a competitive interaction among b. burgdorferi strains within a host in which the order of appearance of the strain is the main determinant of the competitive outcome. this pattern is well described by the classic priority effect in the ecological literature. the observed priority effect could shape bacterial diversity in nature, with consequences in epidemiology and evolution of the disease. borrelia burdgorferiospcpriority effecttransmissioncoinfectionstrain competitionapparent competitionissue-copyright-statement© the author 2015
 within-host microbial communities and interactions among microbes are increasingly recognized as important factors influencing host health and pathogen transmission. the microbial community associated with a host is indeed influenced by a complex network of direct and indirect interactions between the host and the lineages of microbes it harbors, but the mechanisms are rarely established. we investigated the within-host interactions among strains of borrelia burgdorferi, the causative agent of lyme disease, using experimental infections in mice. we used a fully crossed-design with three distinct strains, each group of hosts receiving two sequential inoculations. we used data from these experimental infections to assess the effect of coinfection on bacterial dissemination and fitness as well as the effect of coinfection on host immune response compared to single infection.  the infection and transmission data strongly indicate a competitive interaction among b. burgdorferi strains within a host in which the order of appearance of the strain is the main determinant of the competitive outcome. this pattern is well described by the classic priority effect in the ecological literature. in all cases, the primary strain a mouse was infected with had an absolute fitness advantage primarily since it was transmitted an order of magnitude more than the secondary strain. the mechanism of exclusion of the secondary strain is an inhibition of the colonization of mouse tissues, even though 29% of mice showed some evidence of infection by secondary strain. contrary to expectation, the strong and specific adaptive immune response evoked against the primary strain was not followed by production of immunoglobulins after the inoculation of the secondary strain, neither against strain-specific antigen nor against antigens common to all strains. hence, the data do not support a major role of the immune response in the observed priority effect.  the strong inhibitory priority effect is a dominant mechanism underlying competition for transmission between coinfecting b. burgdorferi strains, most likely through resource exploitation. the observed priority effect could shape bacterial diversity in nature, with consequences in epidemiology and evolution of the disease. electronic supplementary material the online version of this article contains supplementary material, which is available to authorized users. keywords borrelia burdgorferiospcpriority effecttransmissioncoinfectionstrain competitionapparent competitionissue-copyright-statement© the author 2015 