# Levels of pathogen virulence and host resistance both shape the antibody response to an emerging bacterial disease

Daisy Gates, Molly Staley, Luc Tardy, Mathieu Giraudeau, Geoffrey Hill, Kevin Mcgraw, Camille Bonneaud

► To cite this version:

Daisy Gates, Molly Staley, Luc Tardy, Mathieu Giraudeau, Geoffrey Hill, et al.. Levels of pathogen virulence and host resistance both shape the antibody response to an emerging bacterial disease. Scientific Reports, 2021, 11 (1), 10.1038/s41598-021-87464-9. hal-03418447

HAL Id: hal-03418447

https://hal.science/hal-03418447

Submitted on 7 Nov 2021

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Check for updates

OPEN

# Levels of pathogen virulence and host resistance both shape the antibody response to an emerging bacterial disease

Daisy E. Gates<fnref n="1" />, Molly Staley²,⁴, Luc Tardy<fnref n="1" />, Mathieu Giraudeau<fnref n="1" />,³,⁵, Geoffrey E. Hill², Kevin J. McGraw³ & Camille Bonneaud<fnref n="1" />✉

Quantifying variation in the ability to fight infection among free-living hosts is challenging and often constrained to one or a few measures of immune activity. While such measures are typically taken to reflect host resistance, they can also be shaped by pathogen effects, for example, if more virulent strains trigger more robust immune responses. Here, we test the extent to which pathogen-specific antibody levels, a commonly used measure of immunocompetence, reflect variation in host resistance versus pathogen virulence, and whether these antibodies effectively clear infection. House finches (*Haemorhous mexicanus*) from resistant and susceptible populations were inoculated with > 50 isolates of their novel *Mycoplasma gallisepticum* pathogen collected over a 20-year period during which virulence increased. Serum antibody levels were higher in finches from resistant populations and increased with year of pathogen sampling. Higher antibody levels, however, did not subsequently give rise to greater reductions in pathogen load. Our results show that antibody responses can be shaped by levels of host resistance and pathogen virulence, and do not necessarily signal immune clearance ability. While the generality of this novel finding remains unclear, particularly outside of mycoplasmas, it cautions against using antibody levels as implicit proxies for immunocompetence and/or host resistance.

Within any population of host, individuals typically vary in their ability to fight infection. Although understanding this variation is key to predicting host-pathogen interactions, from epidemiological dynamics to their co-evolutionary trajectories, quantifying it in free-living host populations presents a serious challenge<fnref n="1" />⁻⁴. Often, it is only possible to take one or a few measurements of immune molecules or cells that are then used as proxies of an individual's ability to clear infection². There are reasons to question, however, whether such proxies are always good predictors of resistance. For instance, variation in the immune responses measured may, instead, reflect the level of virulence of the pathogen strain, with more virulent strains eliciting more robust responses⁵⁻⁸. Furthermore, across infected individuals, the immune responses measured may not signal clearance ability if their primary role is to prevent infection establishment, or if different pathogen strains vary in their ability to evade immunity⁹. While differences in the immune response to a given pathogen clone are likely to reflect among-host variation in resistance, any natural variation in the pathogen, including in its virulence, may complicate the relationship between immune activity and resistance. Determining the effectiveness and adaptive benefits of immune activity will, as a result, require functional tests that are conducted across a range of different pathogen strains.

Circulating levels of antibodies (immunoglobulins) are often used to estimate immunocompetence in natural populations because they can be relatively easily measured¹⁰. Antibodies are produced by B lymphocytes in response to specific antigens and can control invading pathogens through various mechanisms, including neutralization, prevention of adherence, enhancement of phagocytosis through opsonisation or activation of the complement cascade¹¹. In natural populations of Soay sheep (*Ovis aries*), for example, total nematode-specific antibody levels in blood plasma were found to negatively correlate with parasite egg counts in the faeces of lambs over a 25-year period, suggesting that antibody levels are a useful metric of resistance in this system¹²,¹³.

<footnote n="1">¹Centre for Ecology and Conservation, University of Exeter, Penryn, Cornwall TR10 9FE, UK. ²Department Biological Science, Auburn University, Auburn, Alabama 36849-5414, USA. ³School of Life Sciences, Arizona State University, Tempe, AZ 85287-4501, USA. ⁴Present address: Biology Department, Loyola University Chicago, Chicago, IL 60660-1537, USA. ⁵Present address: Centre for Ecological and Evolutionary Research On Cancer, UMR CNRS/IRD/UM 5290 MIVEGEC, 34394 Montpellier, France. ✉email: c.bonneaud@exeter.ac.uk</footnote>
However, antibody levels against the fungal disease *Pseudogymnoascus destructans* were found to be lower in two European *Myotis* bat species that have coevolved with the pathogen than in newly exposed North American little brown bats (*Myotis lucifugus*), suggesting that, in this case, the antibody response is not consistently an accurate predictor of survival to infection¹⁴. Together these studies highlight the complex relationship between antibody responses and their protective benefits to hosts.

We can use the recently established host-pathogen system consisting of the house finch (*Haemorhous mexicanus*) from North America and its conjunctivitis-causing bacterium, *Mycoplasma gallisepticum*, to test the role of host resistance versus pathogen virulence in shaping among-host variation in immune activity. *M. gallisepticum* emerged in house finches in 1994 following a single host shift from poultry¹⁵⁻¹⁷, causing a devastating epidemic that spread rapidly throughout the eastern United States (U.S.)¹⁸,¹⁹. In disease-exposed populations, the intensity of pathogen-induced selection led to the spread of resistance from standing genetic variation within 12 years of epidemic outbreak²⁰,²¹, although some western populations have remained unexposed to present day²². The evolution of host resistance, in turn, favoured the evolution of increasing pathogen virulence²¹, with late-epidemic isolates causing greater putative host mortality and more severe signs of conjunctivitis than earlier ones²³. As a result, this system presents the necessary variation in host resistance and pathogen virulence required for testing the contribution of each in explaining variation in the immune response to infection.

Here, we conducted a test on the systemic antibody (immunoglobulin Y-IgY) response of house finches to *M. gallisepticum*, which is commonly used to determine past or current infections¹⁶,²⁴. IgYs are the most abundant circulating antibody isotype in birds, reptiles and amphibians, and the functional equivalent of mammalian IgG²⁵. This antibody class immobilises pathogens through binding via agglutination or coating pathogen surfaces (referred to as opsonisation), and can also activate the complement cascade²⁶. We measured serum levels of *M. gallisepticum*-specific IgYs in samples obtained from an infection experiment in which house finches from *M. gallisepticum*-exposed populations that have evolved resistance, as well as from unexposed populations that are susceptible, were inoculated with *M. gallisepticum* isolates of varying virulence and collected over a 20-year period from outbreak²¹. A previous study showed that, following inoculation, finches from resistant and susceptible populations were just as likely to become infected and displayed similar pathogen loads 8 days post-inoculation, but that finches from resistant populations subsequently cleared the infection²⁷.

In our experiment, > 50 pathogen isolates sampled in 9 separate years (1-9 isolates collected per year over a 20-year period) were each inoculated into two house finches, one from resistant and one from susceptible populations. Doing so allowed us to test the extent to which host resistance and pathogen virulence shape the antibody response using a regression-based approach. We favoured a regression-based approach (i.e. maximising the number of isolates) over a replicate-based approach (i.e. maximising the number of replicates per isolate), because we wished to encapsulate sufficient natural variation in both host resistance and pathogen virulence to quantify statistically the relative impacts of each in shaping antibody levels. This aim is at odds with the alternative replicate-based approach for two reasons. First, while the replicate-based approach is preferred for clarifying differences among a limited set of isolates (e.g. A vs. B vs. C), it is not appropriate for quantifying the relative explanatory power of two or more predictors. Second, since in naturally evolving systems one can only replicate the pathogen isolate, there is a danger of spuriously increasing the relative importance of pathogen traits relative to host traits, as replicates will necessarily reduce natural variation and measurement error. Nevertheless, we maximised the power of our regression-based approach using a paired-experimental design in which each isolate was inoculated into two house finches of the same age (i.e. ~ 6 months old) and that had never been exposed to the pathogen previously, but with one originating from resistant and the other from susceptible populations (see "Methods" section).

We then tested the following hypotheses. (1) We tested whether infection intensity (i.e. pathogen load) was associated with serum levels of *M. gallisepticum*-specific antibodies. Because pathogen load is an outcome of the host-pathogen interaction, we then tested whether any association between pathogen load and antibody levels was driven by differences in host resistance and/or pathogen virulence. Antibody levels would reflect resistance if finches from resistant populations displayed higher levels of specific IgYs than finches from susceptible populations. By contrast, antibody levels would reflect variation in pathogen virulence if late-epidemic isolates (i.e. those that are more virulent) elicited higher antibody responses than early-epidemic ones. (2) We then tested the extent to which antibody responses were associated with a change in pathogen load over the course of the infection. Here antibody levels would be indicative of immune clearance ability if individuals displaying higher antibody responses subsequently exhibited the greatest decline in pathogen load.

## Results

**Antibody levels and associating factors.** Mean antibody level over the course of the experiment was 1.85 ± 0.26 ELISA Units (EU)/ml (range 1.49–3.10). Antibody levels were highest at 28 dpi and more comparable on 14 and 35 dpi (linear mixed model; $\chi^2=43.2, p<0.0001$; estimate (14 vs. 28 dpi) ± SE = 0.09 ± 0.01; estimate (14 vs. 35 dpi) ± SE = 0.02 ± 0.01; Figs. 1 and S1). Despite changes in antibody levels between 14 and 35 dpi, antibody levels remained weakly but significantly repeatable within individuals (R = 0.38 ± 0.06; likelihood ratio test = 201.7, *p* < 0.0001; birds from susceptible populations only: R = 0.38 ± 0.09; birds from resistant populations only: R = 0.38 ± 0.09). Although we found no effect of body mass on antibody levels on any of the day measured, sex significantly associated with antibody levels at 28 and 35 dpi, with males displaying significantly higher antibody levels on those days than females (Table 1). As a consequence, all statistical results are provided with sex included as a cofactor.

We tested whether infection intensity (i.e. pathogen load) predicted subsequent antibody levels, and whether any effect of pathogen load was driven by differences in host resistance and/or pathogen virulence. We found that bacterial load at 8 dpi significantly predicted antibody levels at 14 dpi (linear mixed model;
<figure><img src="image_3.png" /><figcaption>**Figure 1.** Variation in antibody levels in both sexes over the course of the infection. We measured antibody levels (in ELISA Units (EU)/ml) at 14, 28 and 35 days post-inoculation (dpi). Boxplots show the median and inter-quartile range, whiskers show maximum values and points represent raw data (dots = females, triangles = males).</figcaption></figure>

estimate ± SE = 0.0007 ± 0.0002, $\chi^2$ = 7.7, *p* = 0.005). To investigate the contribution of variation in host resistance and/or pathogen load in driving the association between bacterial load and antibody levels, we fitted whether the hosts were from susceptible or resistant populations²⁷ (hereafter: host resistance status) and year of pathogen sampling²³ as our explanatory variables, since all other metrics of infection (e.g. clinical symptoms severity) will reflect the outcome of the interaction between host resistance and pathogen virulence. The association between bacterial load at 8 dpi and antibody levels at 14 dpi was explained by a significant effect of both host resistance

<table><caption>**Table 1.** Effects of sex and body mass on antibody levels at 14, 28 and 35 dpi. We ran linear models with antibody levels at 14 dpi or sqrt(antibody levels at 28 or 35 dpi) as the response term and with sex and body mass at the start of the experiment as explanatory terms.</caption><thead><tr><th>Response term</th><th>Explanatory term</th><th>Estimate ± SE</th><th>t statistics</th><th>p value</th></tr></thead><tbody><tr><td rowspan="2">Antibodies at 14 dpi</td><td>Sex</td><td>0.05 ± 0.04</td><td>t<sub>2.96</sub> = 1.3</td><td>0.19</td></tr><tr><td>Initial body mass</td><td>0.02 ± 0.015</td><td>t<sub>2.96</sub> = 1.4</td><td>0.16</td></tr><tr><td rowspan="2">Antibodies at 28 dpi</td><td>Sex</td><td>0.06 ± 0.02</td><td>t<sub>2.97</sub> = 2.7</td><td>0.008*</td></tr><tr><td>Initial body mass</td><td>-0.004 ± 0.009</td><td>t<sub>2.97</sub> = -0.5</td><td>0.64</td></tr><tr><td rowspan="2">Antibodies at 35 dpi</td><td>Sex</td><td>0.04 ± 0.02</td><td>t<sub>2.97</sub> = 2.5</td><td>0.014*</td></tr><tr><td>Initial body mass</td><td>-0.009 ± 0.007</td><td>t<sub>2.97</sub> = -1.2</td><td>0.23</td></tr></tbody></table>

status and year of pathogen sampling (Table 2). Indeed, antibody levels at 14 dpi were significantly higher in birds from resistant relative to susceptible populations (resistant = 1.8 ± 0.2 EU/ml; susceptible = 1.7 ± 1.4 EU/ml), and they increased significantly as the pathogen was sampled progressively later in the epidemic and hence became more virulent (Fig. 2; Table 2).

These results were largely mirrored at subsequent time points. While bacterial load at 21 dpi did not predict antibody levels on day 28 dpi (linear mixed model; estimate ± SE = 0.0002 ± 0.0002, $\chi^2$ = 1.6, *p* = 0.20), bacterial load at 28 dpi significantly and positively predicted antibody levels at 35 dpi (linear mixed model; estimate ± SE = 0.0008 ± 0.0001, $\chi^2$ = 26.2, *p* < 0.0001). In addition, antibody levels at 35 dpi were significantly associated with the year of pathogen sampling but not host resistance status (Table 2), with a similar effect of year of pathogen sampling in resistant and susceptible host populations (linear mixed model; year of pathogen sampling × host resistance status interaction effect: estimate ± SE = 0.002 ± 0.002, $\chi^2$ = 1.4, *p* = 0.23). Taken together, our results indicate that bacterial load shaped subsequent antibody levels, and that the host resistance status and the year of pathogen sampling both drove this association. Indeed, although individuals from both resistant and susceptible populations were able to mount an antibody response, the response of birds from resistant populations was greater than that of birds from susceptible populations. Further, antibody levels increased as the pathogen isolate was sampled later in the epidemic and was, hence, more virulent.

**Antibodies levels and immune clearance ability.** We tested whether antibody responses indicated immune clearance ability and gave rise to a subsequent decline in pathogen load. A main effect of antibody levels at 14 dpi on the change in bacterial load between 14 and 21 dpi would not suggest causation because of
<table><caption>**Table 2.** Effects of host resistance and pathogen virulence on antibody levels. We ran linear mixed models with either antibody levels at 14 dpi or sqrt(antibody levels at 28 or 35 dpi) as the response term, and with pathogen load at 8, 21 or 28 dpi, host resistance status (resistant vs. susceptible populations) and the year of pathogen sampling (1994–2015) as explanatory terms; sex was included as a cofactor and pathogen isolate identity as a random term.</caption><thead><tr><th>Response term</th><th>Explanatory term</th><th>Estimate ± SE</th><th>χ<sup>2</sup></th><th>p value</th></tr></thead><tbody><tr><td rowspan="4">Antibody levels at 14 dpi</td><td>Bacterial load at 8 dpi</td><td>0.0004 ± 0.0003</td><td>2.1</td><td>0.015</td></tr><tr><td>Host resistance status</td><td>-0.09 ± 0.03</td><td>8.4</td><td>0.004*</td></tr><tr><td>Pathogen year</td><td>0.007 ± 0.003</td><td>8.2</td><td>0.004*</td></tr><tr><td>Sex</td><td>0.03 ± 0.03</td><td>1.0</td><td>0.3</td></tr><tr><td rowspan="4">Antibody levels at 35 dpi</td><td>Bacterial load at 28 dpi</td><td>0.0007 ± 0.0002</td><td>20.3</td><td>&lt;0.0001*</td></tr><tr><td>Host resistance status</td><td>-0.01 ± 0.01</td><td>1.0</td><td>0.33</td></tr><tr><td>Pathogen year</td><td>0.003 ± 0.001</td><td>8.3</td><td>0.004*</td></tr><tr><td>Sex</td><td>0.04 ± 0.01</td><td>8.3</td><td>0.004*</td></tr></tbody></table>

the correlations that exist between antibody levels at 14 dpi and bacterial load at 8 dpi (see above). By contrast, causation would be suggested if individuals exhibiting the same pathogen load, but differing antibody levels at 14 dpi, exhibited contrasting changes in bacterial load between 14 and 21 dpi. For example, antibody levels would suggest immune clearance ability if individuals with high bacterial loads and high antibody levels at 14 dpi showed a reduction in bacterial load between 14 and 21 dpi that was greater than those with high bacterial loads and low antibody levels at 14 dpi. However, we found no evidence that the antibody response at 14 dpi was indicative of the change in bacterial load between 14 and 21 dpi (Table 3). Similarly, we found no evidence that the antibody response at 14 dpi was indicative of the change in bacterial load between 14 and 28 dpi (Table 3). Our measures of antibody levels therefore do not appear reflective of the ability of house finches to clear *M. gallisepticum* infection.

## Discussion

By experimentally inoculating house finches from coevolved (resistant) and unexposed (susceptible) populations with pathogen isolates sampled over 20-year of the epidemic during which virulence is known to have increased²³, we were able to test the roles of host resistance and pathogen virulence in shaping the antibody response. We found that all finches mounted an antibody response and that finches from resistant populations mounted a significantly higher antibody response than those from susceptible populations. In addition, antibody levels increased positively with year of pathogen sampling. In other words, although antibody levels should reflect levels of resistance when all birds are exposed to the same isolate, the association between antibody levels and resistance will be obscured when multiple pathogen isolates of differencing virulence are involved. Finally, we found that antibody levels were, however, not indicative of immune clearance ability. Indeed, there was no evidence to suggest that antibody levels on our measurement days shaped pathogen load because, for a given level of pathogen load, individuals with higher antibody levels were not subsequently better able to clear the infection. Taken together, our results suggest that antibody production is reflective of both the levels of host resistance and pathogen virulence, and cannot, in our case, be used as a proxy for immune clearance ability.

Our results are unlikely to be confounded by differences in the birds' prior exposure to the pathogen over their short life before capture, or by genetic differences between the isolates of *M. gallisepticum* used. First, regarding the potential issue of differences in prior exposure, we ensured that (a) none of the birds used in this study was infected with *M. gallisepticum* at the time of capture; and that (b) none had *M. gallisepticum*-specific antibodies before experimental onset. Together, these indicate that none of the birds had been previously exposed and that the antibody response to experimental inoculation was a primary, not a secondary, response to infection. Second, differences in the *M. gallisepticum*-specific IgY levels detected using our method are unlikely to be shaped by genetic differences between isolates for three reasons. (1) Coating the plate microwells with inactivated whole *M. gallisepticum* of poultry origin provides numerous antigens with which house finch antibodies can interact, thereby reducing the likelihood that variation in antibody levels will reflect differences in the ability to cross-react and bind the coated antigens. (2) Using the same methodology, we previously showed that house finches inoculated with chicken and house finch *M. gallisepticum* both mounted detectable antibody responses, but the response to the house finch isolate was greater—a pattern which is opposite to the one expected if the assay had been preferentially biased in favour of antibodies directed against chicken *M. gallisepticum*²⁸. (3) If antibody levels were a function of the genetic distance between the house finch *M. gallisepticum* isolate used for inoculation and the poultry isolate used for coating the microwells, then we should see a decrease in antibody levels with years since the jump from poultry into house finches, as *M. gallisepticum* diverged away from its poultry ancestor and adapted to its novel host. However, this was not the case and, in fact, antibody levels were instead found to increase as the isolate used for inoculation was sampled progressively later in the epidemic (see “Results” section).

That *M. gallisepticum*-specific antibodies are associated not only with the level of resistance of the host population, but also with the level of pathogen virulence, sheds light on previous findings in this system. By inoculating house finches from resistant populations with four *M. gallisepticum* isolates, Grodio and colleagues²⁹ found that the highest antigen-specific serum and lachrymal antibody levels were produced in response to the isolate that
<figure><img src="image_4.png" /></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 2.** Association between antibody levels (in EU/ml), year of pathogen sampling (1994–2017) and host population (resistant and susceptible) at (A) 14 dpi, and (B) 35 dpi (statistics provided in Table 2). Points represent raw values; lines are predicted from the model (solid = resistant host populations; dashed = susceptible host populations), with standard errors represented by ribbons. Note that points were jittered for clarity.</figcaption></figure>

caused the most severe signs of conjunctivitis, and lowest following exposure to the isolate that caused the least. Our results confirm that variation in pathogen virulence can shape levels of *M. gallisepticum*-specific antibodies, although the mechanism by which it does so remains to be determined.
<table><caption>**Table 3.** Effects of antibody levels on changes in pathogen load. We ran linear mixed models with the difference in pathogen load between 21 and 14 dpi or between 28 and 14 dpi as the response term, and with antibody levels at 14 dpi, pathogen load at 14 dpi and their interaction as explanatory terms; sex was included as a cofactor and pathogen isolate identity as the random term.</caption><thead><tr><th>Response term</th><th>Explanatory term</th><th>Estimate ± SE</th><th>χ<sup>2</sup></th><th>p value</th></tr></thead><tbody><tr><td rowspan="4">Change in bacterial load between 14 and 21 dpi</td><td>Antibody levels at 14 dpi</td><td>-13.1 ± 31.1</td><td>1.3</td><td>0.26</td></tr><tr><td>Bacterial load at 14 dpi</td><td>0.03 ± 0.61</td><td>89.4</td><td>&lt;0.0001*</td></tr><tr><td>Antibody levels at 14 dpi × bacterial load at 14 dpi</td><td>-0.35 ± 0.32</td><td>1.3</td><td>0.26</td></tr><tr><td>Sex</td><td>-10.6 ± 9.5</td><td>1.1</td><td>0.29</td></tr><tr><td rowspan="4">Change in bacterial load between 14 and 28 dpi</td><td>Antibody levels at 14 dpi</td><td>1.3 ± 29.9</td><td>0.03</td><td>0.86</td></tr><tr><td>Bacterial load at 14 dpi</td><td>-0.58 ± 0.58</td><td>12.8</td><td>&lt;0.0001*</td></tr><tr><td>Antibody levels at 14 dpi × bacterial load at 14 dpi</td><td>-0.11 ± 0.3</td><td>0.1</td><td>0.71</td></tr><tr><td>Sex</td><td>-6.4 ± 10.2</td><td>0.34</td><td>0.56</td></tr></tbody></table>

A greater antibody response has been found to be associated with increased resistance to *M. gallisepticum* infection in poultry. For example, vaccination of chickens with an avirulent strain of *M. gallisepticum* (GT5) conferred secondary immunity to a subsequent infection with the virulent strain R<sub>low</sub>, with vaccinated chickens displaying lower tracheal lesion scores than sham-immunised controls, as well as a thickening of the trachea due to the infiltration of leukocyte cells<sup>30</sup>. This reduced lesion and inflammatory response of GT5-vaccinated chicken was associated with a *M. gallisepticum*-specific serum IgY response that was both significantly more rapid (as early as 4 days post-challenge) and reached higher levels than that of sham-immunised controls. In our study, none of the finches used had previously been exposed to the pathogen, which means that they were all mounting primary, rather than secondary, immune responses to infection. Nevertheless, house finches from populations that have evolved resistance through immune clearance mounted similarly higher serum *M. gallisepticum*-specific IgY responses at 14, 28 and 35 days post-inoculation than finches from susceptible populations. This suggest that the antibody response contributes to shaping resistance upon first and secondary exposure to *M. gallisepticum*. However, our lack of evidence for an association between antibody levels and subsequent declines in pathogen load also suggests that antibodies do not play a central role in the immune clearance process.

That antibodies play no or a limited role in pathogen clearance has previously been suggested for other mycoplasma species<sup>31</sup>. Antigen-specific antibody levels may not mediate clearance if antibodies serve other protective functions or if the pathogen is able to avoid antibody detection or binding. There are three lines of evidence regarding why such processes could occur during *M. gallisepticum* infection. First, in poultry, antigen-specific antibodies may be more important in preventing the establishment of *M. gallisepticum* infection than in mediating its subsequent clearance. Antigen-specific IgYs present at the infection site (i.e. the respiratory mucosa) and either produced in situ or originating from serum, are indeed thought to mediate immunity to *M. gallisepticum* by preventing bacterial attachment to host tissues<sup>32, 33</sup>. Second, *M. gallisepticum* may be able to avoid antibody detection. Antibodies are thought to have limited protective benefits against intracellular pathogens, with infections by such pathogens thought to be better controlled through cell-mediated immune responses<sup>34</sup>. In accordance, *M. gallisepticum* has been shown capable of invading avian cells and persisting in the intracellular environment<sup>35</sup>, and the jump from poultry into house finches was associated with an increased ability to do both<sup>36</sup>. Third, *M. gallisepticum* may be able to prevent antibody binding. *Mycoplasma* spp. have been shown to possess a specific two-protein system (Mycoplasma Ig binding protein and Mycoplasma Ig protease; MIB-MIP) that captures and cleaves IgYs<sup>37</sup>. Homologs of the MIB-MIP system were found in *M. gallisepticum*<sup>37</sup>, but whether this system allows the pathogen to evade detection by house finch antibodies remains to be tested. Further work is required to understand the precise function of antigen-specific serum antibodies in mediating the resistance that house finches have evolved to *M. gallisepticum*.

In conclusion, we show that antibody production in response to a bacterial pathogen increases positively with levels of both host resistance and pathogen virulence, but fails to be associated with subsequent declines in pathogen load. Our results stress the need for caution when using a limited number of measures of immune activity as estimates of host resistance, even when such measures correlate with variation in disease severity. Further studies are needed to test the degree to which antibody levels confer protection against infection or reflect infection severity in natural populations of hosts, as well as to determine the extent to which our results can be generalised to other pathogens, particularly non-mycoplasma ones.

## Methods

**Capture, housing and experimental inoculation.** In summer 2015, wild hatch-year house finches were trapped from urban and suburban sites in an area of Arizona where finches are known to have never been exposed to *M. gallisepticum* (i.e. susceptible populations)<sup>38</sup>. Birds were weighed, banded, and immediately brought back to Arizona State University, where they were housed in an indoor aviary for the remainder of the experiment<sup>21</sup>. During the same time period, we also caught hatch-year house finches from urban and suburban sites in Alabama. House finches in Alabama have been exposed to *M. gallisepticum* since initial disease outbreak and have evolved resistance. Birds were immediately brought back to aviaries at Auburn University<sup>21</sup>. In all birds from Arizona and Alabama, the lack of past and current infection with *M. gallisepticum* was verified by testing for the presence of: (1) *M. gallisepticum*-specific antibodies using the serum plate agglutination
assay¹⁶, and (2) *M. gallisepticum* bacteria using PCR amplification of *M. gallisepticum* DNA in choanal swabs³⁹. As expected, none of the birds tested positive for either test in Arizona²². Birds positive for either test in Alabama were released immediately; this ensured that: (1) none of the birds had any *M. gallisepticum*-specific IgYs in circulation before experimental onset; (2) antibody responses measured in this study were primary (not secondary) responses to infection; and (3) individuals remaining were not a biased sample of the population (i.e. those that were resistant to the infection). All birds from Arizona (N = 59; 29 females, 30 males) and all birds from Alabama that tested negative for both tests (N = 53; 29 females, 24 males) then underwent a >30-day quarantine period. During this time, and although none of the birds displayed any sign of infection with other diseases, they were all treated prophylactically for *Trichomonas gallinae* with carnidazole (Spartrix, Janssen/Elanco) and *Isospora* spp. with sulfadimethoxine. At the end of this quarantine period, all birds from Alabama were transported in an air-conditioned vehicle to the aviary at Arizona State University. Care was used to minimise travel time (< 30 h), movement and stress to the birds; food and water was provided ad libitum throughout the trip and the birds were regularly checked for any signs of distress or injury. The results obtained in this study are unlikely to be impacted by automobile travel, as the population-difference in the response to infection measured in this experiment was equivalent to the one detected in an earlier study conducted at Auburn University (Alabama), and for which finches caught in Arizona had been driven the same distance in the opposite direction to Alabama²⁰. Following the arrival of birds from Alabama at Arizona State University, all birds were allowed to acclimate for one additional month prior to experimental onset and provided with ad libitum food and water throughout. All birds were weighed (± 0.01 g) at the start and end of the experiment using a top-pan balance.

Each house finch was inoculated with one of 55 *M. gallisepticum* isolates via eye drops²¹. Isolates were originally obtained from naturally infected, wild-caught house finches and collected over a 20-year period across 8 different states in the eastern U.S. (mainly from Alabama)²¹. None of the isolates had been passed in culture more than 3 times⁴⁰ and were stored at -80 °C. Isolates were administered via 20 µL of culture containing 1 × 10⁴ to 1 × 10⁶ colour changing units/mL of *M. gallisepticum* in both eyes, which corresponds to concentrations ranging from 4.1 × 10⁵ to 3.0 × 10⁶ bacterial cells/µL (average ± SE = 1.4 × 10⁶ ± 0.6 × 10⁶ bacterial cells/µL)²³. We verified that variation in the number of bacterial cells inoculated (i.e. inoculation dose) did not affect antibody levels on any of the three days they were measured (linear models; antibodies at 14 days post-inoculation (dpi): estimate ± SE = -2.0 × 10⁻⁹ ± 1.6 × 10⁻⁹, t₁,₉₃ = -1.4, p = 0.18; antibodies at 28 dpi: estimate ± SE = -2.7 × 10⁻⁹ ± 2.8 × 10⁻⁹, t₁,₉₄ = 1.0, p = 0.34; antibodies at 35 dpi: estimate ± SE = 3.7 × 10⁻¹⁰ ± 2.0 × 10⁻⁹, t₁,₉₄ = 0.2, p = 0.86). The experiment was stopped at 35 dpi and all birds were euthanized as stipulated by home office licensing. Protocols were approved by Institutional Animal Care and Use Committees (IACUC) at Auburn University (permit # PRN 2015-2721) and Arizona State University (permit #15-1438R), by Institutional Biological Use Authorizations to Auburn University (# BUA 500), and by the University of Exeter's ethics committee. The experiment was performed in strict accordance with relevant guidelines and regulations, including the USA federal regulations of the Animal Welfare Act (AWA, 7 USC 2131), the Institutional Animal Care & Use Committees, and the Animals (Scientific Procedures) Act 1986 (UK). This study is reported in accordance with ARRIVE guidelines (https://arriveguidelines.org).

**Pathogen load.** Pathogen load was measured by quantitative amplification of *M. gallisepticum* DNA from conjunctival and tracheal swabs obtained at 8, 14, 21 and 28 dpi²³. DNA was extracted using a QIAGEN DNeasy Blood and Tissue Kit according to manufacturer's protocols. Multiplex quantitative PCRs for *mgc2* and *rag1* in each sample was then conducted using an Applied Biosystems StepOnePlus Real-Time PCR system²³. Reactions were run in 20 µL volume containing: 2 µL of sample genomic DNA template, 1 µL each of 10 µM *mgc110-F/R* and *rag1-102-F/R* primers (total 4ul), 0.5 µL each of 10 µM *Mgc110-JOE* and *Rag1-102-6FAM* fluorescent hydrolysis probes (total 1 µL), 10 µL of 2x *qPCRBIO Probe Mix HI-ROX (PCR BIOSYSTEMS)* and 5 µL Nuclease-free water (Ambion). Reactions were then run at 95 °C for 3 min, followed by 45 cycles of 95 °C for 1 s and 60 °C for 20 s. Samples were run in duplicate alongside serial dilutions of plasmid-based standards. Amplification data was exported to LinRegPCR v.2017.1 for calculation of individual reaction efficiencies and quantification of low-amplification samples⁴¹, ⁴²; between-run variation was normalised using Factor qPCR v.2016.0⁴³, with plasmid standard serial dilutions used for factor correction. Pathogen load was then determined as the number of *M. gallisepticum* cells divided by the number of house finch cells to control for variation in sampling efficiency⁴⁴.

**Antibody levels.** At 14, 28 and 35 dpi, ~75 µL of whole blood was taken by brachial venipuncture and serum separated through centrifugation at 2000XG for 15 min at 4 °C before being stored at -80 °C. Goat Polyclonal Passerine IgY-heavy and light chain Antibody (Bethyl Laboratories Inc., TX, USA) was diluted to 1:1000 dilution according to manufacturer protocol. The prepared antibody was then conjugated to a visible proprietary HRP (horse radish peroxidase) ligand for visualization in enzyme-linked immunosorbent assay (ELISA), according to instructions accompanying the Lighting-Link HRP Conjugation Kit (Innova Biosciences Ltd., Cambridge, UK) for an antibody:HRP molar ratio of 1:4. *M. gallisepticum*-specific antibody detection was then performed with ELISA, with methods adapted from¹⁰, ²⁸. Finch serum samples were diluted to 1:2400 in 1x sample conjugate diluent (Affinitech LTD, USA) and 100 µL of each sample was added to plate microwells coated with inactivated whole *M. gallisepticum* of poultry origin (Affinitech LTD). A standard curve of pooled sample was run on each plate in a twofold dilution series from 1:200 to 1:6400. Samples were run in duplicate, and care was taken to ensure that duplicates were distant from one-another on the plate to minimise variation caused by evaporation and edge effects. Plates were washed three times in 1x wash buffer (Bethyl Laboratories Inc.), before adding 100 µL of the prepared HRP-conjugated Passerine IgY-heavy and light chain antibody, diluted to a 1:10,000 concentration with sample conjugate diluent (Bethyl Laboratories Inc). After incubation for one hour in darkness,
the plate was washed three times and 100 µL of TMB one component HRP microwell substrate (Bethyl Laboratories Inc.) was added to each well and incubated for 15 min in darkness. The reaction was stopped with 100 µL of 0.18 M H₂SO₄ and the plate was read on a plate reader at 450 nm within 30 min of stopping the reaction.

**Statistical analyses.** All statistical analyses were conducted in R 3.3.2⁴⁵. We used the package lme4⁴⁶ to perform linear mixed models. The significance of explanatory terms were assessed by comparing models with and without each term using likelihood ratio tests. Figures were produced using the package ggplot2⁴⁷. Out of the 112 birds experimentally-inoculated, we removed 4 that died during the course of the experiment due to incomplete data. To create a fully balanced paired design, we then considered only the isolates that were each inoculated into one bird from resistant populations and one from susceptible populations, and hence removed a further 8 birds from the analyses. We thus analysed data from a total of 50 pairs of birds (resistant populations: 26 females and 24 males; susceptible populations: 25 females and 25 males), each inoculated with one of 47 pathogen isolates (note that 3 isolates—2 from 1995 and one from 2007—were inoculated into 2 pairs of birds each).

Antibody levels on all days combined were natural log transformed, and antibody levels at 28 and 35 dpi were square root transformed, to fit a normal distribution. To determine whether antibody levels varied over the course of the infection, we ran a linear mixed model with ln(antibody level) as the response term, day post-inoculation (dpi) as the explanatory term, and pathogen isolate identity as the random term. Tests of the repeatability of antibody levels within individuals were conducted using the rpt function⁴⁸, with ln(antibody levels) as the response term and house finch identity as the explanatory term. We tested the associations between antibody levels and sex or body mass using 3 linear models with antibody levels at 14 dpi or sqrt(antibody levels at 28 or 35 dpi) as the response term and with sex and body mass at the start of the experiment as explanatory terms (note: inclusion of isolate identity as a random term did not qualitatively affect the outcome of those models).

To test whether pathogen load predicted antibody levels, we ran linear mixed models with either antibody levels at 14 dpi or sqrt(antibody levels at 28 or 35 dpi) as the response term, and with pathogen load at 8, 21 or 28 dpi as an explanatory term; sex was included as a cofactor and pathogen isolate identity as a random term to maintain the paired design in the model. To then test the contribution of variation in host resistance and/or pathogen load in driving any association between bacterial load and antibody levels, we re-ran those models including host resistance status (resistant vs. susceptible populations) and the year of pathogen sampling (1994–2015) as additional explanatory terms. We used host resistance status to test for an effect of host resistance, and year of pathogen sampling for an effect of pathogen virulence, rather than metrics taken during the infection (e.g. clinical symptoms severity, pathogen load), because such metrics are unavoidably shaped by both host resistance and pathogen virulence. Finally, to determine whether antibody levels shaped changes in pathogen load, we ran linear mixed models with the difference in pathogen load between 21 and 14 dpi or between 28 and 14 dpi as the response term, and with antibody levels at 14 dpi, pathogen load at 14 dpi and their interaction as explanatory terms; sex was included as a cofactor and pathogen isolate identity as the random term.

**Data availability**

Data reported in this paper have been deposited in Dryad Digital Repository https://doi.org/10.5061/dryad.x69p8czht.

Received: 1 December 2020; Accepted: 24 March 2021

Published online: 15 April 2021

**References**

1. Biard, C., Monceau, K., Motreuil, S. & Moreau, J. Interpreting immunological indices: the importance of taking parasite community into account. An example in blackbirds *Turdus merula*. *Methods Ecol. Evol.* **6**, 960–972. https://doi.org/10.1111/2041-210x.12371 (2015).

2. Boughton, R. K., Joop, G. & Armitage, S. A. O. Outdoor immunology: methodological considerations for ecologists. *Funct. Ecol.* **25**, 81–100. https://doi.org/10.1111/j.1365-2435.2010.01817.x (2011).

3. Maizels, R. M. & Nussey, D. H. Into the wild: digging at immunology’s evolutionary roots. *Nat. Immunol.* **14**, 879–883. https://doi.org/10.1038/ni.2643 (2013).

4. Martin, L. B., Weil, Z. M. & Nelson, R. J. Refining approaches and diversifying directions in ecoimmunology. *Integr. Comp. Biol.* **46**, 1030–1039. https://doi.org/10.1093/icb/icl039 (2006).

5. Johnson, W. et al. Pathogenic and humoral immune responses to porcine reproductive and respiratory syndrome virus (PRRSV) are related to viral load in acute infection. *Vet. Immunol. Immunopathol.* **102**, 233–247. https://doi.org/10.1016/j.vetimm.2004.09.010 (2004).

6. Ortiz, R. H. et al. Differences in virulence and immune response induced in a murine model by isolates of *Mycobacterium ulcerans* from different geographic areas. *Clin. Exp. Immunol.* **157**, 271–281. https://doi.org/10.1111/j.1365-2249.2009.03941.x (2009).

7. Sela, U., Euler, C. W., da Rosa, J. C. & Fischetti, V. A. Strains of bacterial species induce a greatly varied acute adaptive immune response: the contribution of the accessory genome. *PLoS Pathog.* https://doi.org/10.1371/journal.ppat.1006726 (2018).

8. Skjesol, A. et al. IPN with high and low virulence: host immune responses and viral mutations during infection. *Virol. J.* https://doi.org/10.1186/1743-422x-8-396 (2011).

9. Hornef, M. W., Wick, M. J., Rhen, M. & Normark, S. Bacterial strategies for overcoming host innate and adaptive immune responses. *Nat. Immunol.* **3**, 1033–1040. https://doi.org/10.1038/ni1102-1033 (2002).

10. Fassbinder-Orth, C. A. et al. Immunoglobulin detection inwild birds: effectiveness of three secondary anti-avian IgY antibodies in direct ELISAs in 41 avian species. *Methods Ecol. Evol.* **7**, 1174–1181. https://doi.org/10.1111/2041-210x.12583 (2016).

11. Janeway, C. *Immunobiology: The Immune System in Health and Disease* (Garland Science, 2005).

12. Coltman, D. W., Pilkington, J., Kruuk, L. E. B., Wilson, K. & Pemberton, J. M. Positive genetic correlation between parasite resistance and body size in a free-living ungulate population. *Evolution* **55**, 2116–2125 (2001).

13. Hayward, A. D. et al. Natural selection on individual variation in tolerance of gastrointestinal nematode infection. *PLoS Biol.* https://doi.org/10.1371/journal.pbio.1001917 (2014).
14. Johnson, J. S. et al. Antibodies to *Pseudogymnoascus destructans* are not sufficient for protection against white-nose syndrome. *Ecol. Evol.* **5**, 2203–2214. https://doi.org/10.1002/ece3.1502 (2015).

15. Fischer, J. R., Stallknecht, D. E., Luttrell, M. P., Dhondt, A. A. & Converse, K. A. Mycoplasmal conjunctivitis in wild songbirds: the spread of a new contagious disease in a mobile host population. *Emerg. Infect. Dis.* **3**, 69–72. https://doi.org/10.3201/eid0301.970110 (1997).

16. Luttrell, M. P., Fischer, J. R., Stallknecht, D. E. & Kleven, S. H. Field investigation of *Mycoplasma gallisepticum* infections in house finches (*Carpodacus mexicanus*) from Maryland and Georgia. *Avian Dis.* **40**, 335–341 (1996).

17. Delaney, N. F. et al. Ultrafast evolution and loss of CRISPRs following a host shift in a novel wildlife pathogen, *Mycoplasma gallisepticum*. *Plos Genet.* https://doi.org/10.1371/journal.pgen.1002511 (2012).

18. Dhondt, A. A., Tessaglia, D. L. & Slothower, R. L. Epidemic mycoplasmal conjunctivitis in house finches from Eastern North America. *J. Wildl. Dis.* **34**, 265–280 (1998).

19. Nolan, P. M., Hill, G. E. & Stoehr, A. M. Sex, size, and plumage redness predict house finch survival in an epidemic. *Proc. R. Soc. Lond. Ser. B Biol. Sci.* **265**, 961–965 (1998).

20. Bonneaud, C. et al. Rapid evolution of disease resistance is accompanied by functional changes in gene expression in a wild bird. *Proc. Natl. Acad. Sci. U.S.A.* **108**, 7866–7871 (2011).

21. Bonneaud, C. et al. Rapid antagonistic coevolution in an emerging pathogen and its vertebrate host. *Curr. Biol.* **28**, 2978–2983 (2018).

22. Staley, M., Bonneaud, C., McGraw, K. J., Vleck, C. M. & Hill, G. E. Detection of *Mycoplasma gallisepticum* in House Finches (*Haemorhous mexicanus*) from Arizona. *Avian Dis.* **62**, 14–17 (2018).

23. Tardy, L., Giraudeau, M., Hill, G. E., McGraw, K. J. & Bonneaud, C. Contrasting evolution of virulence and replication rate in an emerging bacterial pathogen. *Proc. Natl. Acad. Sci. U.S.A.* **116**, 16927–16932. https://doi.org/10.1073/pnas.1901556116 (2019).

24. Grodio, J. L., Buckles, E. L. & Schat, K. A. Production of house finch (*Carpodacus mexicanus*) IgA specific anti-sera and its application in immunohistochemistry and in ELISA for detection of *Mycoplasma gallisepticum*-specific IgA. *Vet. Immunol. Immunopathol.* **132**, 288–294 (2009).

25. Warr, G. W., Magor, K. E. & Higgins, D. A. IgY—clues to the origins of modern antibodies. *Immunol. Today* **16**, 392–398. https://doi.org/10.1016/0167-5699(95)80008-5 (1995).

26. Diebolder, C. A. et al. Complement is activated by IgG hexamers assembled at the cell surface. *Science* **343**, 1260–1263. https://doi.org/10.1126/science.1248943 (2014).

27. Bonneaud, C. et al. Evolution of both host resistance and tolerance to an emerging bacterial pathogen. *Evol. Lett.* **3**, 544–554. https://doi.org/10.1002/evl3.133 (2019).

28. Staley, M., Hill, G. E., Josefson, C. C., Armbruster, J. W. & Bonneaud, C. Bacterial pathogen emergence requires more than direct contact with a novel passerine host. *Infect. Immun.* **86**, 9. https://doi.org/10.1128/iai.00863-17 (2018).

29. Grodio, J. L. et al. Pathogenicity and immunogenicity of three *Mycoplasma gallisepticum* isolates in house finches (*Carpodacus mexicanus*). *Vet. Microbiol.* **155**, 53–61. https://doi.org/10.1016/j.vetmic.2011.08.003 (2012).

30. Javed, M. A. et al. Correlates of immune protection in chickens vaccinated with *Mycoplasma gallisepticum* strain GT5 following challenge with pathogenic *M-gallisepticum* strain R-low. *Infect. Immun.* **73**, 5410–5419 (2005).

31. Dumke, R. & Jacobs, E. Antibody response to *Mycoplasma pneumoniae*: Protection of host and influence on outbreaks?. *Front. Microbiol.* **7**, 7. https://doi.org/10.3389/fmicb.2016.00039 (2016).

32. Avakian, A. P. & Ley, D. H. Protective immune-response to *Mycoplasma gallisepticum* demonstrated in respiratory-tract washings from *M-gallisepticum*-infected chickens. *Avian Dis.* **37**, 697–705. https://doi.org/10.2307/1592017 (1993).

33. Yagihashi, T. & Tajima, M. Antibody-responses in sera and respiratory secretions from chickens infected with *Mycoplasma gallisepticum*. *Avian Dis.* **30**, 543–550. https://doi.org/10.2307/1590419 (1986).

34. Glatman-Freedman, A. & Casadevall, A. Serum therapy for tuberculosis revisited: reappraisal of the role of antibody-mediated immunity against *Mycobacterium tuberculosis*. *Clin. Microbiol. Rev.* **11**, 514. https://doi.org/10.1128/cmr.11.3.514 (1998).

35. Vogl, G. et al. *Mycoplasma gallisepticum* invades chicken erythrocytes during infection. *Infect. Immun.* **76**, 71–77. https://doi.org/10.1128/iai.00871-07 (2008).

36. Dowling, A. J., Hill, G. E. & Bonneaud, C. Multiple differences in pathogen-host cell interactions following a bacterial host shift. *Sci. Rep.* **10**, 6779 (2020).

37. Arfi, Y. et al. MIB-MIP is a mycoplasma system that captures and cleaves immunoglobulin G. *Proc. Natl. Acad. Sci. U.S.A.* **113**, 5406–5411. https://doi.org/10.1073/pnas.1600546113 (2016).

38. Staley, M., Bonneaud, C., McGraw, K. J., Vleck, C. M. & Hill, G. E. Detection of *Mycoplasma gallisepticum* in House Finches (*Haemorhous mexicanus*) from Arizona. *Avian Dis.* https://doi.org/10.1637/11610-021317-RegR (2018).

39. Roberts, S. R., Nolan, P. M., Lauerman, L. H., Li, L. Q. & Hill, G. E. Characterization of the mycoplasmal conjunctivitis epizootic in a house finch population in the southeastern USA. *J. Wildl. Dis.* **37**, 82–88 (2001).

40. Papazisi, L. et al. GapA and CrmA coexpression is essential for *Mycoplasma gallisepticum* cytadherence and virulence. *Infect. Immun.* **70**, 6839–6845 (2002).

41. Ruijter, J. M. et al. Amplification efficiency: linking baseline and bias in the analysis of quantitative PCR data. *Nucleic Acids Res.* **37**, 12. https://doi.org/10.1093/nar/gkp045 (2009).

42. Tuomi, J. M., Voorbraak, F., Jones, D. L. & Ruijter, J. M. Bias in the C-q value observed with hydrolysis probe based quantitative PCR can be corrected with the estimated PCR efficiency value. *Methods* **50**, 313–322. https://doi.org/10.1016/j.ymeth.2010.02.003 (2010).

43. Ruijter, J., Villalba, A., Hellemans, J., Untergasser, A. & van den Hoff, M. Removal of between-run variation in a multi-plate qPCR experiment. *Biomol. Detect. Quantif.* **5**, 10–14 (2015).

44. Grodio, J. L., Dhondt, K. V., O'Connell, P. H. & Schat, K. A. Detection and quantification of *Mycoplasma gallisepticum* genome load in conjunctival samples of experimentally infected house finches (*Carpodacus mexicanus*) using real-time polymerase chain reaction. *Avian Pathol.* **37**, 385–391. https://doi.org/10.1080/03079450802216629 (2008).

45. R: A language and environment for statistical computing (R Foundation for Statistical Computing, Vienna, 2016).

46. Bates, D., Maechler, M., Bolker, B. & Walker, S. Fitting linear mixed-effects models using lme4. *J. Stat. Softw.* **67**, 1–48. https://doi.org/10.18637/jss.v067.i01 (2015).

47. ggplot2: Elegant Graphics for Data Analysis (Springer, 2009).

48. Repeatability Estimation for Gaussian and Non-Gaussian Data v. 0.9.21 (2018).

## Acknowledgements

This research was supported by a Natural Environment Research Council standard grant to C.B. (NE/M00256X) and a Genetics Society Fieldwork Heredity grant to D.E.G. We thank A. Russell, P. Sirand-Pugnet and Y. Arfi for helpful discussion and/or constructive comments on the manuscript. We thank M. Cook for assisting with bird captures in Arizona, A. Santos and undergraduates in the lab of W. R. Hood for assisting with bird captures in Alabama, and A. K. Ziegler for assisting with the experiment in Arizona.
**Author contributions**

C.B. conceived and designed the study. G.E.H., K.J.M., M.G., M.S., C.B. obtained the animals and/or bacterial isolates. M.G., K.J.M., L.T. conducted the experiment. D.E.G., M.S. quantified antibody levels. L.T. quantified bacterial load. C.B., D.E.G. analysed the data and wrote the paper. All authors provided comments and approved the final submission of the manuscript.

**Competing interests**

The authors declare no competing interests.

**Additional information**

**Supplementary Information** The online version contains supplementary material available at https://doi.org/10.1038/s41598-021-87464-9.

**Correspondence** and requests for materials should be addressed to C.B.

**Reprints and permissions information** is available at www.nature.com/reprints.

**Publisher's note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_6.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2021