Heterozygosity and fitness in a threatened songbird: blood parasite infection is explained by single-locus but not genome-wide effects

In non-pedigreed populations, insights into effects of inbreeding can be obtained by correlations between individual heterozygosity and fitness-related traits (HFCs). Using an information-theoretic approach, we explored whether heterozygosity of microsatellite markers, measured as internal relatedness (IR), is associated with infection by blood parasites (Plasmodium, Trypanosoma, or Leucocytozoon) in the threatened Aquatic Warbler (Acrocephalus paludicola). We also explored whether any of the markers is more influential than others, or than IR, in explaining blood parasitism (single-locus effects). While we observed that IR was a relatively important predictor of Plasmodium parasitism, we did not find strong evidence for IR to correlate with infection by the identified blood parasites, accounting for sex and population effects. Therefore, our data did not support negative inbreeding effects on blood parasite infection in the Aquatic Warbler. However, we found single-locus effects, such that individuals heterozygous at AW-03 and Ase19 had lower probability of infection by blood parasites pooled together and by Plasmodium, respectively. This indicates that these two markers are in linkage disequilibrium with unknown fitness loci which are related to resisting or clearing blood parasites, and which confer a heterozygote advantage in the Aquatic Warbler. Our results add to the growing evidence that single-locus effects contribute more to HFCs than formerly recognized and have implications for Aquatic Warbler conservation.


Introduction
Heterozygosity and fitness correlations (HFCs) relate individual genetic variation to traits associated with fitness, giving insight into inbreeding depression and population viability (Reed and Frankham 2003;Grueber et al. 2008). HFCs have been found for survival, clutch size, breeding success, growth rates, and parasite infestation (Coltman et al. 1999;Rowe et al. 1999;Slate et al. 2000;Townsend et al. 2009;White et al. 2015;Cézilly et al. 2016). HFC research allows for a better understanding of evolutionary ecology of a species and can support its conservation (Grueber et al. 2008).
HFCs have been studied using molecular markers, such as microsatellites. They arise due to genome-wide effects (the general effect hypothesis) and/or single-locus effects (the local effect hypothesis) (Hansson and Westerberg 2002). General effects are observed when the heterozygosities of functional and marker (i.e., neutral) loci are associated (identity disequilibrium, ID), as they both respond in a similar way to inbreeding. Thus, if fitness is impaired by the lowered heterozygosity of genes, there is a correlation between fitness and marker heterozygosity. For example, in the Soay Sheep (Ovis aries), microsatellite heterozygosity was inversely correlated to parasitism by gastrointestinal nematodes and mortality caused by these parasites (Coltman et al. 1999). Local effects are found when the heterozygosity of a marker locus is related to that of a fitness-associated locus (linkage disequilibrium, LD), for instance due to physical linkage (Szulkin et al. 2010). The neutral marker hitchhikes with the selected locus and influences the heterozygosity metric more than do other markers under study. For example, heterozygosity at a single microsatellite locus predicted whether tuberculosis was localized or spread in the body of Wild Boars (Sus scrofa) (Amos and Acevedo-Whitehouse 2009).
HFCs have been used to study inbreeding especially in animal populations, which lack information on pedigrees to calculate the individual coefficient of inbreeding (F). F defines the probability that two alleles at a locus are identical by descent. In HFC studies on non-pedigreed populations, the underlying assumption is that F is inversely related to heterozygosity at a number of molecular markers (multi-locus heterozygosity, MLH). This assumption has been subject to debate. On the one hand, the association between MLH and F was found to be weak when a low number of markers were used (Balloux et al. 2004;Slate et al. 2004). However, a more recent study (Forstmeier et al. 2012) made an important point that there is discrepancy between F calculated using a pedigree and the realized identity by descent, IBD (which determines the true rate of inbreeding), due to random events in Mendelian segregation. This discrepancy, rather than insufficient marker number, could explain the weak relationship between F and MLH. The strength of this relationship is also determined by variance in inbreeding, marker polymorphism, and the level of IBD (Hansson and Westerberg 2002;Csilléry et al. 2006;Ljungqvist et al. 2010;Taylor et al. 2010;Forstmeier et al. 2012). Overall, the literature shows that even at a modest number (e.g., up to 11), microsatellite markers can reliably detect HFCs that are attributable to inbreeding (Rowe et al. 1999;Slate et al. 2000;Hansson et al. 2004;Seddon et al. 2004;Forstmeier et al. 2012). To ensure legitimate inferences, it is important to test whether the markers of choice indicate inbreeding and whether there is variance in inbreeding, and to distinguish between the general and the local effects (Balloux et al. 2004;Pemberton 2004;David et al. 2007;Szulkin et al. 2010).
An important trait associated with fitness is resistance to parasites (Hamilton and Zuk 1982;Morand and Poulin 2000;Marzal 2012). In birds, some of the key parasites are vector-transmitted protozoans that invade the bloodstream. These include mainly the intra-erythrocytic hematozoans of the genera Plasmodium, Haemoproteus, and Leucocytozoon, and the extracellular kinetoplastid of the genus Trypanosoma (Atkinson and Van Riper III 1991;Zídková et al. 2012). The fitness costs of blood parasitism are mediated by extensive damage in a variety of tissues, including destruction of erythrocytes, and activation of the immune system (Molyneux et al. 1983;Atkinson and Van Riper III 1991;Valkiūnas 2005). Infection with Plasmodium, the agent of avian malaria, negatively correlates with hatching success, brood provisioning rates, and fledging success (Knowles et al. 2010), survival (Stone et al. 1971;Beier et al. 1981;Krama et al. 2015;Ilgūnas et al. 2016;Townsend et al. 2018b), predation risk (Møller and Nielsen 2007;Krama et al. 2015), song quality and output (Gilman et al. 2007), and fat score during migration (Shurulinkov et al. 2012). In passerines, experimental medication with anti-malarial drugs showed that hematozoan blood parasites negatively affect clutch size, hatching success, parental provisioning rates, nestling growth, fledging success, and survival (Merino et al. 2000;Marzal et al. 2005;Knowles et al. 2010;Martínez-de La Puente et al. 2010). Infections with trypanosomes have been associated with later arrival to breeding grounds in spring, siring fewer offspring, and less effective nest defence against a predator relative to uninfected individuals (Rätti et al. 1993;Hakkarainen et al. 1998;Dyrcz et al. 2005).
As in birds hemoparasite infection is an important fitness determinant, it is crucial to identify its underlying variables. One of the traits with a potential effect on infection occurrence is the genetic diversity of the host Frankham 2005). However, the link between genetic diversity and blood parasitism in birds is not well understood, because studies that addressed it have been scarce and yielded mixed results (MacDougall-Shackleton et al. 2005;Vallender et al. 2012;Boerner et al. 2013;Ferrer et al. 2014;Townsend et al. 2018a). In addition, only one study focused on a species with a small population (Ortego et al. 2007). In species with a small and fragmented population, studying HFCs is crucial, as it can benefit their conservation (Reed and Frankham 2003;Grueber et al. 2008).
Here, using an information-theoretic approach, we investigated whether the probability of infection with vector-transmitted blood parasites is predicted by (1) host MLH and/or, (2) single-locus effects of specific microsatellite loci in the Aquatic Warbler (Acrocephalus paludicola). This globally threatened passerine, a habitat specialist breeding in isolated wetlands, underwent a steep decline in population size (Briedis and Keišs 2016;Flade et al. 2018), which could have lowered its genetic diversity and increased inbreeding rates. The Aquatic Warbler is infected by blood parasites of the genera Leucocytozoon, Plasmodium, and Trypanosoma; in males, infection with trypanosomes negatively correlates with fitness-related traits (Dyrcz et al. 2005;Neto et al. 2015;Kubacka et al. 2019). These findings, converged with the previous results on haemosporidian parasites in other passerines (Merino et al. 2000;Marzal et al. 2005;Knowles et al. 2010;Martínez-de La Puente et al. 2010), indicate that blood parasitism likely affects the fitness status of the Aquatic Warbler. Therefore, this passerine comes as a good model to study HFCs related to blood parasite infection in species of conservation concern.

Study area and sampling
We sampled two core breeding populations of the Aquatic Warbler in Eastern Poland, differing by size and degree of fragmentation. The Biebrza Valley population (central position: 53°17′10.6′′ N, 22°33′41.9′′ E) is one of the largest (c. 2700 singing males), and occupies extensive and rather continuous mire habitat. The Polesie population (central position: 51°21′48.9′′ N, 23°17′21.4′′ E), located c. 250 km away, holds approximately 600 singing males and breeds at several small mires scattered across a large area.
In total, we obtained blood samples from 144 Aquatic Warblers between May and August 2014 (see Online Resource 1 for a list of samples). For the purpose of this analysis, five immature individuals were excluded, as their sex was unknown and because age affects blood parasite infection (Marzal et al. 2016). Hence, the sample sizes were 50 in the Biebrza Valley (17 females, 33 males) and 89 in Polesie (18 females, 71 males). The birds were sexed by the presence of cloacal protuberance (males) or brood-patch (females) and aged by plumage (Svensson 1992). We collected the blood by puncturing the brachial vein and stored the samples in 96% ethanol until laboratory analysis.
Both extraction methods that we applied produce yields of comparable DNA quality and quantity (Sameer et al. 2009;Mirmomeni et al. 2011). The pattern of null alleles in our sample (see Table 1) was similar to the one found in a sample collected from Aquatic Warblers in the Biebrza Valley and Polesie in another project and typed at the same microsatellite loci but with a different PCR protocol (Kubacka unpubl.). Thus, we had confidence that employing We used Genealex 6.5 Smouse 2006, 2012) to check the markers for deviation from the Hardy-Weinberg rule and the Genepop web interface (Raymond et al. 1995;Rousset 2008) to test for linkage disequilibrium, applying the sequential Bonferroni correction to the obtained p values. We assessed null allele frequency using the Microchecker and the Oosterhout estimator, which is well suited for nonequilibrium populations (Van Oosterhout et al. 2006).

Heterozygosity and identity disequilibrium
As a heterozygosity metric, we used the internal relatedness (IR). IR is a measure of parental relatedness that weighs heterozygosity by allele frequency. It is calculated with the formula: where H is the number of homozygous loci in an individual, N is the total number of loci, and f i is the frequency of the ith allele in the genotype. IR assumes values between -1 (all heterozygous) and 1 (all homozygous). Negative IR values indicate more outbred and positive IR values suggest more inbred individuals within the studied sample ). In the case of a species with fragmented population and low expected genetic diversity due to a steep recent decline, IR appears to be superior to homozygosity by loci (HL), which is better suited for open populations with relatively high immigration rates and heterozygosity (Aparicio et al. 2006). IR is a measure of heterozygosity that yields the strongest HFCs and is highly correlated with other measures, such as heterozygosity, standardized heterozygosity, and homozygosity by loci (Chapman et al. 2009;Forstmeier et al. 2012).
To test whether the marker loci showed ID, we tested for a heterozygosity-heterozygosity correlation (HHC) (Balloux et al. 2004). The set of marker loci was randomly divided in half to see whether the heterozygosity of the first group is correlated with that of the second group across individuals. This procedure was repeated 1000 times and the mean correlation coefficient was calculated with its confidence interval. If it is significantly greater than zero, a genome-wide effect is inferred that is likely due to inbreeding. In addition, we computed the g 2 statistic, which is an indicator of variance of inbreeding in the population. g 2 quantifies the covariance of heterozygosity at the studied loci standardized by their average heterozygosity (David et al. 2007).
We calculated IR, HHC, and g 2 in the R environment (v. 3.5.1; R Core Team 2018). IR and HHC were quantified using the package Rhh (v. 1.0.2; Alho et al. 2010) and the g 2 statistic was calculated with the package inbreedR (v. 0.3.2; Stoffel et al. 2016), running 1000 permutations.

Identification and diversity of blood parasites
The blood samples were molecularly screened for the presence of blood parasites in a different project (Kubacka et al. 2019). In brief, the DNA samples extracted following the phenol-chloroform protocol were used to PCR-amplify a 478 bp fragment of the mitochondrial cyt b gene to detect Leucocytozoon, Haemoproteus, and Plasmodium (Hellgren et al. 2004), and a 326 bp fragment of the small subunit ribosomal RNA (SSU rRNA) gene to identify Trypanosoma (Sehgal et al. 2001). PCR products of positive samples were purified and sequenced. The obtained sequences were aligned to sequences available in the MalAvi database (Bensch et al. 2009) for Leucocytozoon, Haemoproteus, and Plasmodium, and in GenBank for Trypanosoma (BLASTN, best hit match) (Kubacka et al. 2019).

Statistical analysis
We performed the analysis in the R environment (v. 3.5.1; R Core Team 2018). We used the Akaike information criterion corrected for small-sample size (AICc) and the informationtheoretic approach for inference (Burnham and Anderson 2002). In this approach, in contrast with the null hypothesis testing, the researcher creates an a priori set of biologically plausible candidate models, which are given relative AICc ranks corresponding to how well they explain the data. The best-ranking model (or a subset of best models) is then selected, with models within 2 AICc units relative to the top model (which has the lowest AICc) considered to be highly supported (Burnham and Anderson 2004). Each model is given a quantitative measure of relative support, by computing model likelihood (the relative likelihood of the model in the candidate set, given the data), the Akaike weight (ωAICc; the probability of the model for a given set of candidate models), and the evidence ratio (the ratio of the model's probability to the probability of the top model).
We ran two separate analyses, in both of which we used generalized linear models with binomial error structure and logit link, fitted using the 'glm' function from the stats package. First, we assessed the general heterozygosity effects on infection probability. We constructed three sets of candidate models, with the infection status by blood parasites pooled (candidate set 1), Plasmodium (candidate set 2), and Trypanosoma (candidate set 3) as a binary response variable, respectively. Because infection by Leucocytozoon was detected only in one individual (Kubacka et al. 2019), this genus was included in the pooled parasites models, but we did not run a separate analysis for it. In each of the candidate sets, the global (full) model included IR, as well as region (Biebrza or Polesie), sex (male or female) and all the two-way interactions. This allowed us to control for factors that affect infection rates and/or HFCs, such as population structure, sex, or habitat (Reid et al. 2007;Brekke et al. 2010;Szulkin et al. 2010;Sehgal 2015;Arct et al. 2017;Kubacka et al. 2019). Each candidate set contained a null model, which assumed that infection probability was constant in all individuals. By choosing this all-subset approach, we were also able to quantify the importance of IR in explaining infection, relative to sex and population. The importance of a variable is calculated by summing the Akaike weights across all the models in the set that contain this variable (Burnham and Anderson 2002). Second, to test for single-locus effects, we again constructed three sets of candidate models, with the response variables as above and the explanatory variables being the heterozygosities of the microsatellite loci, coded as 0 (homozygous) or 1 (heterozygous). Since six loci were discarded from the analyses of single-locus effects (see "Results"), there were 11 locus models. Each of the three candidate sets also included a model assuming that infection probability varied only with IR and a null model, yielding a total of 13 models in each candidate set. This enabled us to see whether any single locus explained the given infection status better than did other loci or IR, and better than the null model. Missing values (see Table 1) for a given locus were replaced by its mean heterozygosity (Szulkin et al. 2010). To ensure that missing values did not bias the results, we re-ran all the models taking standardized heterozygosity (SH) instead of IR. SH is best applied to a data set in which not all individuals are genotyped at all the loci. It is calculated by standardizing individual heterozygosity at the genotyped loci by their mean heterozygosity, thus making heterozygosity independent of the loci that are genotyped (Coltman et al. 1999).
In both analyses, the sum contrasts were used, IR was centred, and estimates were model-averaged across the 95% confidence set of the best models (i.e., a set including the top models the cumulative ωAICc of which equals or just exceeds 0.95) to account for model selection uncertainty (Burnham and Anderson 2002). We applied the package AICcmodavg

Heterozygosity and blood parasites: multi-locus effects
Probability of infection by blood parasites pooled was best explained by the constant infection model. The IR model, which was within the ≤ 2 ΔAICc subset, ranked third and had about twice lower probability relative to the null model ( Table 2). The odds ratio (OR) corresponding to the modelaveraged estimate of the IR effect (Table 3) suggested that probability of infection by any blood parasite increases about twofold with one unit increase in IR (Fig. 1a). However, the confidence interval of the IR estimate spanned zero, indicating that the effect of IR is not truly different from null. The relative importance of IR, sex, and population was 0.42, 0.46, and 0.36, respectively. The full model fitted the data well, as there was no significant difference between the fitted and the observed data (X 2 = 3.98, df = 8, p = 0.859). An analysis run on the data set excluding the individual parasitized by Leucocytozoon yielded similar results (not shown).
Infection by Plasmodium was also best explained by the constant infection model, with the IR model ranking second and having a similar probability ( Table 2). The estimate of the IR effect suggested an approximately twofold increase of Plasmodium infection probability with one unit increase in IR (Fig. 1b); however, its confidence intervals also overlapped zero (Table 3). The relative importance of IR, sex, and population in explaining Plasmodium infection was 0.60, 0.61, and 0.40, respectively. The global model fitted the data well (X 2 = 5.97, df = 8, p = 0.651). Infection by Trypanosoma was best explained by the population model, while the only model containing IR within ΔAICc < 2 ranked fifth, with a model probability of 0.08 (Table 2). The estimate of the IR effect was also not truly different from one (Table 3, Fig. 1c). The relative importance of IR, sex, and population in explaining Trypanosoma infection was 0.38, 0.51, and 0.73, respectively. The full model fitted the data well (X 2 = 4.60, df = 8, p = 0.799).
In all the analyses, the order and support for the models were nearly identical when they were run with standardized heterozygosity (SH) instead of IR (not shown).

Heterozygosity and blood parasites: single-locus effects
We found support that infection by blood parasites pooled was negatively affected by heterozygosity at locus AW-03.  Fig. 1 Relationship between internal relatedness (IR, centred) and the probability of infection by a at least one of the blood parasites identified in the study (Plasmodium, Trypanosoma, or Leucocytozoon), b Plasmodium, and c Trypanosoma. The lines were fitted based on the estimates from the top models that included IR (Table 2). Bands show confidence intervals and the circles represent the distribution of samples (0-uninfected; 1-infected). None of the associations was clearly supported by the data (see Tables 2 and 3) The AW-03 model was the only one within ΔAICc < 2; it was 9.8 times more probable than the IR model and 5.6 times more probable than the constant infection model (Table 4). The AW-03 estimate was -0.96 (95% CI -1.81 to -0.12) and its OR was 0.38 (95% CI 0.16-0.89), indicating that individuals heterozygous at the AW-03 locus were 0.38 times less likely to carry blood parasite infection relative to homozygous individuals. Infection by Plasmodium was best explained by the Ase19 and AW-03 models, which were the only ones within ΔAICc < 2. The top model, containing the Ase19 locus, was 4.3 times more probable than the IR model and 3.3 times more probable than the null model ( Table 4). The Ase19 estimate was -0.82 (95% CI -1.6 to -0.05), OR 0.44 (95% CI 0.20-0.95); the AW-03 estimate was -0.74 (95% CI -1.61 to 0.14), OR 0.48 (95% CI 0.20-1.15). Therefore, in individuals heterozygous at Ase19, the probability of infection with Plasmodium was 0.44 of that of homozygous individuals, while the effect of AW-03 was weakly supported, as the CIs of its OR spanned one.
There was considerable model selection uncertainty concerning Trypanosoma infection models, with none of the loci receiving clearly better support than the other loci or IR (Table 4). No estimates of the ΔAICc < 2 models showed an effect on infection probability, as all of their 95% CIs spanned zero (not shown). In all the three analyses, the order and support for models was also very similar when SH was used instead of IR (not shown).

Heterozygosity: genome-wide effects
We did not find strong evidence for an association between multi-locus heterozygosity and the incidence of infection by blood parasites in the Aquatic Warbler, a rare bird species that has gone through a dramatic decline in numbers and geographical range. Therefore, in our study, the relationship between MLH and infection status in this species cannot be convincingly explained by the general effects of heterozygosity. Still, heterozygosity appears to be a relatively important predictor of Plasmodium infection incidence in the Aquatic Warbler. Our results also indicate that heterozygosity has different importance as a predictor of infection by Plasmodium and by Trypanosoma, as it poorly explained infection by the latter.
Studies on birds that related infection by blood parasites to MLH have yielded mixed results. In Blue Tits (Cyanistes caeruleus), MLH at microsatellite markers showed a negative quadratic relationship with the probability of infection by Plasmodium relictum ). In the Common Buzzard (Buteo buteo), MLH was correlated with infection by Leucocytozoon buteonis through interaction with vole abundance and age of individuals (Boerner et al. 2013). In American Crows (Corvus brachyrhynchos), lower heterozygosity was associated with infection of offspring, but not adults, by Plasmodium spp. (Townsend et al. 2018a). However, other avian studies failed to find any relationship (Ortego et al. 2007;Vallender et al. 2012).
Our study was hampered by the lack of polymorphism of some of the microsatellite loci and presence of null alleles, which led to exclusion of 6 out of the 17 markers. The low final number of markers, together with their moderate polymorphism (Table 1), could explain why we did not find firm support for the effect of MLH on blood parasite infection status (Hansson and Westerberg 2002;Balloux et al. 2004;Slate et al. 2004;Ljungqvist et al. 2010;Forstmeier et al. 2012). On the other hand, it was also emphasized that finding an HFC is dependent on the variance in inbreeding rates (Hansson and Westerberg 2002;Csilléry et al. 2006;Ljungqvist et al. 2010;Taylor et al. 2010;Forstmeier et al. 2012). The isolation of breeding sites, the relatively high breeding philopatry, the skewed paternities and the recent steep decline of the Aquatic Warbler (Dyrcz et al. 2002(Dyrcz et al. , 2005Briedis and Keišs 2016;Bellebaum 2018;Flade et al. 2018) indicate that in this species, increased variation in inbreeding could be expected. Nevertheless, we found a low variance in inbreeding (g 2 ), which suggests low variation in relatedness in the studied areas.
It is also possible that MLH better explains quantitative infection rates, rather than infection probability, as the latter does not differentiate between acute and chronic infection. For example, in mountain White-crowned Sparrows (Zonotrichia leucophrys oriantha), the intensity of infection with hematozoan parasites was lower in more heterozygous individuals (MacDougall-Shackleton et al. 2005), and in Banded Mongooses (Mungos mungo), individuals with increased heterozygosity had lower intestinal parasite loads (Mitchell et al. 2017). Infection status might also be related to MLH only in a harsher environment (Boerner et al. 2013;Annavi et al. 2014). Finally, the genetic mechanism underlying resistance to blood parasites could be a trait determined by one or a few genes, making it less likely to observe correlation of genome-wide heterozygosity with blood parasite resistance (Szulkin et al. 2010;Bateson et al. 2016).

Heterozygosity: single-locus effects
We observed local effects of loci AW-03 and Ase19, which predicted a lower infection probability by the studied blood parasites taken together, and by Plasmodium, respectively, when in heterozygous state. This indicates that AW-03 and Ase19 are in linkage disequilibrium with unidentified fitness loci, which confer a heterozygote advantage in the Aquatic Warbler. Linkage disequilibrium between selected and apparently neutral loci can arise due to selection, population sub-structuring, and genetic bottleneck. However, unless recombination is low or suppressed, LD is disrupted with time (Hartl and Clark 2007). The severe decline experienced by the Aquatic Warbler during the last 100-150 years (Briedis and Keišs 2016; Flade et al. 2018) could have increased LD in its genome. Yet, our analyses of the same samples do not support a historical or recent genetic bottleneck in the study populations and no population sub-structuring was found between them (Kubacka unpubl.). Consequently, it is plausible that the association between heterozygosity of Ase19 and AW-03 and infection status in the Aquatic Warbler stems from their linkage with functional loci, such as those related to parasite resistance or infection clearance. Earlier research provides support for the association between Ase19 and immunity. In nestlings of the Bluethroat (Luscinia svecica), heterozygosity at this marker positively correlated with wing-web swelling following injection of phytohaemagglutinin (PHA) (Fossøy et al. 2009). In birds, response to PHA is associated with both innate and adaptive (T-cells) components of the immune system (Palacios et al. 2009;Vinkler et al. 2010;Salaberria et al. 2013). It is also positively related to the probability of recovery from infection with Haemoproteus sp. (Gonzalez et al. 1999).
In general, our results add to the increasing evidence that single-locus effects could underlie HFCs more than it was acknowledged earlier. Single-locus effects have been reported in previous HFC studies on wild-ranging animals (Hansson et al. 2004;Lieutenant-Gosselin and Bernatchez 2006;Ortego et al. 2007;Kaeuffer et al. 2008;Fossøy et al. 2009;Vangestel et al. 2011;García-Navas et al. 2014;Minias et al. 2015;Judson et al. 2018). However, only a few HFC studies related to disease documented local effects (Acevedo-Whitehouse et al. 2006Amos and Acevedo-Whitehouse 2009) and, to our knowledge, our study is the first to show them in a diseaserelated study on a bird.
The local effects of Ase19 and AW-03 could help indicate genomic regions that are under parasite-mediated selection, as pointed out by Acevedo-Whitehouse et al. (2009). In addition, as they suggest a heterozygote advantage at functional genes, they also imply that factors affecting genetic diversity in the Aquatic Warbler, such as the constant loss of habitats and the resulting genetic drift, could impact blood parasite resistance in this species. Finally, these local effects can be applied in conservation management of the species to identify individuals or populations resistant to blood parasites, for instance in translocation projects (Grueber et al. 2008).

Conclusion
While we found that internal relatedness at microsatellite loci appears to be an informative predictor of Plasmodium infection in the Aquatic Warbler, our results do not provide a solid support for a negative inbreeding effect on blood parasitism in this species. However, they demonstrate a case for single-locus effects on disease probability in birds, which, to our knowledge, has not been reported before. For the Aquatic Warbler, the observed local effects could have meaningful conservation implications. In general, they add to the accumulating evidence that local effects could underpin apparent genome-wide HFCs.
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://creat iveco mmons .org/licen ses/by/4.0/.