Leukocyte profile variation in Dupont’s Lark (Chersophilus duponti) in Spain and Morocco

Stress in birds has been widely studied through the measurement of heterophil-to-lymphocyte ratio (H/L ratio). In this study we aimed to assess for the first time the potential variation of stress, measured as H/L ratio, associated to geography (between-country variation) and seasonality (between seasons and within the breeding season), as well as the leukocyte profiles, in the threatened Dupont’s Lark (Chersophilus duponti), using samples from Spain and Morocco. Furthermore, we tested whether variation in H/L ratio was associated with variables such as population density, presence of blood parasites and individual body condition. We found that H/L ratio did not vary between countries, but individuals captured during the breeding season showed higher values of H/L compared to non-breeding ones. Neither male density, nor date within the breeding season had an effect on the H/L ratio. In Spain, individuals with higher body condition showed lower H/L ratio regardless of whether they were malaria-infected. In Morocco, malaria-infected individuals showed higher values of H/L ratio than the non-infected birds. Moreover, we found that our average values of H/L ratio in Morocco were within the ranges of other passerines, but not in Spain. Individuals with higher H/L ratios may be more stressed or present higher capability to face stressful situations. Although H/L ratio is a useful and relatively easy way to obtain measure of stress, the impact that the environment might have on stress and the way it is explained by H/L ratio must be addressed carefully. This study provides new insight for this species’ biology and provides useful reference information to test the status and survival of other populations.


Introduction
In many species, the physiological stress of individuals is often the first compensatory response to intrinsic or extrinsic stressors-including those mediated by human activitiesaimed to preserve homeostasis and improve the chances of survival and reproduction (Dantzer et al. 2014). Living in risky or stressful environments might chronically activate the physiological stress response of individuals. Therefore, studying and quantifying the physiological responses of individuals to stressors can serve to understand whether individual-level responses impact the viability of populations (Wikelski and Cooke 2006;Cooke et al. 2013) and stress metrics can help to assess the status and future fate of natural populations (Dantzer et al. 2014;McCormick and Romero 2017;Szott et al. 2020).
Levels of health status or stress in birds is often measured through blood parameters, like hormone concentration or leukocyte cell counts. Hormones (e.g. corticosterone) have been widely used to assess physiological stress (Davis et al. 2008;Cockrem 2013;Skwarska 2018). However, their concentrations in blood might change rapidly (1-3 min) following disturbance (short-term stress indicator), which makes its use challenging (Scanes et al. 2020) when working with individuals in the wild, whose control and manipulation is not immediate. Alternatively, the use of haematological parameters, such as the heterophil-to-lymphocyte ratio (hereafter H/L ratio) quantified in blood smears, is considered a reliable measure of long-term physiological stress on birds (Davis and Maney 2018;Gormally and Romero 2020), as they vary after some time (1-4 h depending on the species and the stressor's nature) since exposure to the stressor and they tend to endure longer than corticosterone levels in blood (Vleck et al. 2000;Davis et al. 2008;Müller et al. 2011;Lebigre et al. 2012;Goessling et al. 2015; Davis and Maney. 2018;Skwarska 2018;Minias 2019;Latimer et al. 2020).
Heterophils are phagocytic leukocytes involved in responses to inflammation, infection, and stress (also known as innate defence). Lymphocytes, on the other side, are the largest cell group of the immune system and are involved in immunoglobulin production and control of the immune defence (adaptive immune defence) (Davis et al. 2008;Banbura et al. 2013;Skwarska 2018). In case of exposure to an external factor, causing an injury or infection, a hormonemediated alteration of the number of leukocytes in blood is produced, so that the number of heterophils increases to fight the infection (Müller et al. 2011), and the number of lymphocytes decreases to avoid an excessive expenditure of energy (Goessling et al. 2015; Granthon and Williams 2017;Minias 2019). Thus, high values of H/L result from an increase in the production of heterophils and/or a concomitant reduction in lymphocytes as a response to stressors. Among the common stressors identified as having an impact on the H/L ratio of birds are: fear (Campo et al. 2005), handling procedures (Cīrule et al. 2012;Lebigre et al. 2012;Scanes et al. 2020), reproduction costs (Vleck et al. 2000;Ochs and Dawson 2008;Jakubas et al. 2011;Lebigre et al. 2012;Skawarska 2018) or harsh environmental conditions (e.g. low food availability, water deprivation, bad weather conditions; Ruiz et al. 2002;Ochs and Dawson 2008;Plischke et al. 2010;Skawarska 2018). Therefore, differences in H/L ratio may be expected due to population geographic location, year of sampling or time of the year (Norte et al. 2009a). Variation of H/L ratio has been also related to the stress induced by parasites. Parasites usually drain host physiological resources and are therefore expected to influence host leukocyte profiles. The presence of blood parasites, for example, can trigger changes in the levels of leukocytes in birds (Dunn et al. 2013). Several factors can influence the change in the H/L ratio of an individual, such as the degree of virulence of a given parasite, the degree of parasitemia as well as its interaction with the age, sex and body condition (Norte et al. 2009b;Becker et al. 2019). However, the relationship between H/L ratio and blood parasites is yet to be clarified, as non-consistent responses of H/L ratio under parasite infections have been found (Norte et al. 2009b;Skwarska 2018;Becker et al. 2019). For example, some studies have reported that avian malaria parasites (Haemoproteus, Plasmodium) increase H/L ratios (Norte et al. 2009b), whereas others found decreased H/L ratios (Dunn et al. 2013) or even no evidence of such variations (Cornelius et al. 2014;Granthon and Williams 2017).
Interpreting results from H/L ratios can be challenging because baseline values of the species studied should be considered to interpret stress correctly (Davis and Maney 2018). Furthermore, it is not always clear whether a high H/L ratio is related to the effect of a potential stressor or a sign of good health status (see Davis et al. 2008). In the context of blood parasite infections, for example, high H/L ratios have been interpreted both as an indicator of good ability to control and repel parasite infections (Minias 2019), and as evidence of active infections (Davis et al. 2004;Becker et al. 2019). Contradictory results have been also found on the relationship between H/L ratio and birds' body condition. Ochs and Dawson (2008) showed that female Tree swallows (Tachycineta bicolor) in poorer body condition had higher H/L ratio due to higher levels of stress. On the contrary, Norte et al. (2009c) and Powell et al. (2013) found no significant relationship between the H/L ratio and individual body condition of Great Tits (Parus major) and Noisy Miners (Manorina melanocephala), respectively. Therefore, further research assessing the relationship between the H/L ratio, parasites and the body condition of individuals on different study species under varying conditions (sex and age of individuals, temporality, etc.) is needed to improve our understanding of this topic.
Dupont's Lark (Chersophilus duponti) is a mediumsized passerine (c. 40 g) classified by the IUCN as Vulnerable globally (BirdLife International 2020), and in Europe (BirdLife International 2020). Its distribution is restricted to natural steppes of Spain and North of Africa (García et al. 2008a;Suárez 2010;García-Antón et al. 2019). In Europe, the species has suffered an annual population decline rate of 3.9% between (Gómez-Catasús et al. 2018. In North Africa, this species also seems to be experiencing a severe decline (BirdLife International 2020). However, there is no updated information on this respect. Moroccan and Iberian populations are long-term isolated (García et al. 2008b), live in different ecological context and show morphological differences (García et al. 2008b;Vögeli et al. 2017;García-Anton et al. 2019). These two populations currently seem to differ in size and conservation status, being the Moroccan population seven times larger than the Iberian one (García et al. 2008a). Habitat loss and fragmentation are the main threats for the species (Tella et al. 2005;Suárez 2010; Pérez-Granados and López-Iborra 2013), which partly explain the concentration of individuals at large connected habitat patches with high food availability (Vögeli et al. 2010;Pérez-Granados and López-Iborra 2013;Gómez-Catasús et al. 2019). Despite the current species' decline, there is only one study, based on feather corticosterone and stable isotopes values, assessing the physiological variation of the species in relation to environmental conditions (Fairhurst et al. 2013).
The aim of this study is to investigate whether the H/L ratio in Dupont's Lark is sensitive to both external and internal stressors. First, we examined the potential variation of the H/L ratio between countries, comparing samples from Spain and Morocco to assess whether the countries could be considered as independent units of analysis, and between seasons within each country, to test the hypothesis that stress values are higher during the breeding season, compared to the non-breeding season (Wingfield et al. 1997;Norte et al. 2009a;Pap et al. 2010;Skwarska 2018). During the breeding season, individuals are considered to face higher exposure to environmental stress factors, like the emergence of seasonal parasites, increased energetic expenditure, strong trade-offs between immune function and reproduction, or high conspecific competition. The specific features of Mediterranean shrubsteppes, where the study species inhabits, impose particular sources of stress linked to high predation rates and severe time constraints, such as decreased food availability and increased temperature as the breeding season progresses (see for example Calero-Riestra et al. 2013). Second, we evaluated whether the H/L ratio varied with population density, under the hypothesis that stress increases due to increased social contacts (Jakubas et al. 2011). Finally, we tested whether the H/L ratio was associated with individual's body condition per se, considering the interaction body condition-malaria infection status, expecting that malariainfected individuals and/or individuals with poorer body condition would show higher H/L ratio (Skwarska 2018;Becker et al. 2019).

Study sites
The study was carried out in three different periods: 2005-2006, 2012-2013 and 2017-2018 (Table 1), as part of a long-term monitoring programme started in 2005, on the ecology and conservation of the species in Spain and Morocco.
In Spain, we sampled 19 sites in three different regions (same as in García-Antón et al. 2019)  the Spanish distribution of the species and corresponding to the main breeding range of the species in the Iberian Peninsula, with continental Mediterranean climate (locations 3-18, Fig. 1). For a more detailed information see Table S1.

Field and laboratory procedures
Birds were captured during both the breeding and the non-breeding seasons (see Table 1) using tape luring and spring-traps baited with mealworms (Tenebrio molitor larvae). Breeding season lapsed from March 1st to July 10th in Spain (Herranz et al. 1994;Suárez 2010;Pérez-Granados et al. 2017) and from February 1st to June 19th in Morocco (Thévenot et al. 2003, Garcia et al. 2008a. Outside these periods, captures were considered as belonging to the nonbreeding season. Because the capture method used is inherently biased toward males (Laiolo et al. 2007), we excluded the few females (N = 12) and juveniles (N = 2) captured from analyses. A total of 300 Dupont's Lark males were captured (Table 1) Table 1).
All birds were ringed to avoid collecting repeated samples over the study period. Blood samples were immediately collected after capture of each individual in the field (see below) and afterwards, the following morphological measurements were taken: wing length, length of the 8th primary (P8) and tail length to the nearest 1 mm with a rule, tarsus length (bent toes) using a digital calliper (accuracy 0.01 mm), and body mass using a digital balance to the nearest 0.1 g. Handling time was reduced to the minimum and never exceeded 30 min.
Blood samples were collected from the jugular or brachial vein, using 0.5 ml insulin syringes and capillary tubes. One drop of blood was used to make a standard smear on a microscope slide, which was fixed in 99% ethylic alcohol for 180 s and air-dried. Afterwards, they were stained with Wright Giemsa stain, and stored at room temperature until microscope inspection (Davis 2005;Banbura et al. 2013). We followed procedures for leukocyte counting described in Campbell (1995) and Carr and Rodak (2010). Slides were visually scanned using an optical microscope (OLYMPUS BX41) under 100 × magnification (using an oil-immersion objective), in a zig-zag pattern, to count the types of leukocytes (heterophils, lymphocytes, eosinophils, basophils and monocytes). A total of 100 leukocytes were identified per slide, noting the proportions of each leukocyte type, and the H/L ratio was calculated by dividing the number of heterophils in these 100 leukocytes by the number of lymphocytes. Only fields of view with relatively uniform distribution of cells were examined. All smears were examined by people with the same training and experience (CPG, SM).
Additionally, we investigated the presence of avian malaria parasites (Plasmodium and Haemoproteus) in samples from both countries using the 2005-2006 sample set. Total genomic DNA was extracted from blood using ammonium acetate and ethanol to precipitate proteins, and purified DNA was diluted to a working concentration of 25 ng/µl. The samples were screened for the presence of Plasmodium and Haemoproteus using a nested polymerase chain reactions (PCR) protocol (Waldenström et al. 2004), designed to amplify a 479-bp fragment of the mitochondrial cytochrome b gene of each parasite genus. We evaluated 2.5 µl of each final reaction on 2% agarose gels stained with ethidium bromide and using 1xTAE buffer. We repeated the protocol three times to confirm negative results.
All birds were captured and handled with the corresponding permissions of both regional and national Spanish and Moroccan authorities, and procedures were approved by the Local Ethical Committee for Animal Experiments of the Autónoma University of Madrid (CEI80-1468-A229, see Ethical Statement).

Statistical analysis
Since both environmental conditions and stressors can vary widely geographically and over time, particularly with respect to time of breeding, we firstly investigated whether both factors had an effect on the H/L ratio (dependent variable) using linear mixed-models (LMM). We used country, season (reproductive/non-reproductive period) and their interaction as fixed factors, and the random effect of sampling year was introduced to account for the potentially high correlation of the data from the same year. Pairwise comparisons among countries and seasons were explored by post-hoc Tukey's test using the package emmeans (Length 2021) in R software.
After detecting significant variation of H/L ratio considering the interaction between season and country (see below), and as it has been shown how populations differed both morphologically and genetically between Spain and Morocco because of their isolation (García et al 2008b;García-Antón et al. 2018), we focussed further analysis within each country. Furthermore, as we obtained population density estimates only for the breeding season, we restricted further analyses to this period within each country. A linear regression model (LM) was performed to investigate whether the variation in H/L ratio (response variable) was related to potential stressors in birds' environment during the breeding season, using population density as a proxy for intraspecific competition (Pérez-Granados and López-Iborra 2015) and date (linear and quadratic term) as potential predictors. Population density estimates in each study site (males per 10 ha) were obtained from García et al. (2008a), Suárez (2010) and Traba et al (2019) (see Table S1). Due to the existing variation in timing of reproduction between countries (Herranz et al. 1994), the variable date was relativized to the start of the season in each country (1st of March for Spain and 1st of February for Morocco), as described in Suárez (2010) and Thévenot et al. (2003). The variable year was included as a fixed factor to control for potential inter-annual variability, since the number of years sampled within each country was not high enough to be considered as a random factor. Lastly, we also included longitude (x) and latitude (y) (UTM coordinates of locations where the animals were captured), as well as the longitude*latitude interaction term to account for potential spatial effects.
We finally explored the relationship between the H/L ratio during the breeding season and body condition and infection status of captured birds using data from 2005 and 2006. We ran a linear regression model using H/L ratio as response variable. The explanatory variables for this model were malaria infection (infected/non-infected), body condition, the interaction between them, and date (linear and quadratic term). As in the previous model, we also included the variable year and longitude (x) and latitude (y), as well as the longitude*latitude interaction term to account for potential spatial correlation. We used the residuals of the regression between body mass and structural size as an index of individual body condition, following the procedure described by Peig and Green (2009). We used wing length as a measure of structural body size because it was the morphological measure that showed the highest correlation with body mass in both countries (Pearson's test, Spain: R 2 = 0.38; p < 0.001; Morocco: R 2 = 0.49; p = 0.001).
Final models were chosen based on the deviance statistic for model comparison. Using drop1 function in R (Chambers 1992) we removed those variables that showed the highest p-value, running afterwards a new reduced model and comparing it with the previous one using an Anova (Log-Likelihood ratio Chi-squared test) for the significance of the model. This procedure was carried out until all remaining variables showed a significant effect, and a significant difference between consecutive models tested using Anova function was found.
Lastly, we assessed the variance explained by each model. For the linear mixed model, to know the variance explained by the fixed factors (marginal R 2 ) and from the full model (conditional R 2 , which considers the variance explained by both the fixed and random factors) we used the R2GLMM function from the package "MuMIn" (Barton 2020) following the procedure described by Nakagawa and Schielzeth (2013). On the other hand, the variance explained by the linear regression models was measured as the adjusted R 2 .
Statistical analyses were performed using R software (version 3.6.0) (R Core Team 2019). Mixed model was run using the package "lme4" (Bates et al. 2015). Normality of model residuals was visually checked for all models using QQ plots. The H/L ratio was log-transformed to attain normality for the analyses. Continuous predictors were standardized (mean = 0 and SD = 1) within each country. Standarization was obtained using the function decostand of the R package "vegan" (Oksanen et el. 2020).

Spatio-temporal variation in H/L ratio
The H/L ratio was significantly higher during the breeding season than during the non-breeding season (χ 2 = 10.23; df = 1; p = 0.001) (Fig. 2), and the interaction between season and country had only a marginal effect on H/L ratio (χ 2 = 3.07; df = 1; p = 0.079). Tukey's test showed that differences in H/L ratio between both seasons were significant in Spain, but not in Morocco (Table S2) (Fig. 2). There were no overall differences in H/L ration between both countries (χ 2 = 0.41; df = 1; p = 0.519). The full model explained 21.2% of variation in H/L ratio, though the variance of H/L ratio explained by the fixed factors was 10.9% (R 2 GLMM (m) = 0.109).

Local variation in H/L ratios
In Spain, the H/L ratio varied significantly with the geographic position, as well as with the year of sampling (adjusted R 2 = 0.14; Table 3). H/L ratio values decreased with the interaction between longitude and latitude, being the values lower to the north-east (estimate ± SE = − 2.235e-11 ± 1.043e-11; t = − 2.14; p < 0.05). Population density and date were removed from the initital model, as they were not significant in any country. In Morocco, the final model showed that H/L ratio values varied significantly only with longitude (adjusted R 2 = 0.14), being the values lower to the east (estimate ± SE = − 1.504e-06 ± 5.598e-07; t = − 2.69; p < 0.05; Table 3). The influence of malaria infection, which showed a similar prevalence in both countries, and body condition on H/L ratio differed between countries. In Spain, H/L ratio marginally varied with the interaction body condition*malaria (F 1,32 = 3.16; p = 0.08), showing a significant variation with longitude (F 1,32 = 4.90; p < 0.05) and latitude (F 1,32 = 8.01; p < 0.05) (adjusted R 2 = 0.32). Non-infected individuals showed an overall independence of H/L values with respect to the body condition, whereas infected individuals showed Table 2 H/L ratio and proportion of the different types of leukocytes (heterophils, lymphocytes, eosinophils, monocytes and basophils) (mean ± SD) on Dupont's Lark males at each study site sampled in Spain and Morocco The number between brackets is referred to the code used for each location at the map of the main text ( Fig. 1  a decline in H/L values with increasing values of body condition. The infected individuals with high body condition showed lower H/L ratios than non-infected ones (Table 4, Fig. 3), while infected poor-conditon individuals were predicted to show higher H/L ratios than the group of noninfected individuals. Date was removed from the initial model as no significant variation of H/L ratio was detected. Moreover, as we only had data about malaria infection for birds trapped during one year, the variable year was not considered in this case. We did not detect a statistically significant relationship between avian malaria and body condition in Morocco; H/L values in infected individuals did not vary with body condition and were only slightly higher (marginally non-significant: F 1,32 = 3.24; p = 0.08) than in uninfected individuals. H/L ratio also increased significantly to the west (F 1,32 = 4.77; p < 0.05) (adjusted R 2 = 0.16) ( Table 4).

Discussion
Our study is the first to provide leukocyte profiles and to examine variation in H/L ratio in the threatened Dupont's Lark, as well as the link between H/L ratio and different environmental stressors in this species. We found values of H/L ratio similar to those found in other passerines, such as Great tits (Parus major  We observed a temporal variation in this physiological measure of stress. We also observed different local patterns of response in relation with infections status and body condition. The limitations in the knowledge about the relationship between H/L ratio and health status and the unbalanced nature of our dataset implies some caution about the conclusions, but we believe that our results may inform about the health status of populations of this threatened lark. H/L ratio levels were higher during the breeding season compared to the non-breeding season, especially in Spain. We found no relationship between H/L ratio and population densityFinally, when analysing the variation of H/L ratio in relation with individuals' condition, we did not observe a relationship between avian malaria andbody condition in Morocco. However, in Spain we did show how the presence of malaria affects the H/L ratio of individuals and how this effect could be mitigated with good body condition.

Spatio-temporal variation
Previous studies have demonstrated how Spanish and Moroccan populations of Duponts' Lark have been genetically disconnected for about 24,000 years (García et al. 2008b). That isolation has resulted not only in a genetical differentiation between countries (García et al. 2008b), but also in morphological variation (García et al. 2008b;García-Antón et al. 2018). Furthermore, habitat structure and composition differ between both study areas (García et al. 2008a;Suárez 2010;García-Antón et al. 2019). All these differences, however, were not substantiated in significant differences in H/L ratio between countries, which suggests similar H/L baseline accounting both in Spain and Morocco.
For instance, as regards temporal variation, we found the same pattern of seasonal change of the H/L ratio in both countries, whith higher H/L ratio during the breeding season than during the non-breeding period. However, this variation was significant only in Spain. In this case, the higher H/L ratio during the breeding season could be related to the territorial behaviour of the species, when there is high competition for finding and defending mates and territories (Garza et al. 2005;Pérez-Granados and López-Iborra 2015). However, this result does not match the lack of relationship between H/L ratio and population density (but see below for further explanations), so it could be also explained by energy demanding constraints during breedingother than intrasexual competition. The relationship between breeding activity and the H/L ratio was already shown by Lebigre et al. (2012), who found that the H/L ratio on Black grouse (Tetrao tetrix) increased during the mating season because of the competition between males and higher demands of energy for reproduction. Our findings are also consistent with those by Pap et al. (2010), who found that the H/L ratio of the Great tit increased during the breeding season (probably caused by an increase of energy demand for parental care and/or a higher presence of parasites, and thus higher levels of infestation), to decrease afterwards (less energetically stressful season and/or absence of parasites in winter). During the non-breeding season, Dupont's Lark is less territorial and some individuals perform dispersal movements (Suárez et al. 2006). Therefore, competition among males during the nonbreeding period might be reduced, compared to the breeding season, which may partly explain the lower H/L ratio found during the non-breeding season. Moreover, we also found that the variance explained by our model was higher when considering the random variable year, which could be due to the potential influence of yearly weather variablility and related food supply on the H/L ratio values (Norte et al. 2009a;Banbura et al. 2013;Fairhurst et al. 2013). Finally, we cannot exclude the role of the high nest predation rates reported in this species (over 80% of nests in some Spanish populations, Herranz et al. 1994;Pérez-Granados et al. 2017;Gómez-Catasús et al. 2021) and other berian groundnesting shrubsteppe passerines (> 70%, Suarez et al. 1993, Yanes andSuárez 1995) in explaining the observed H/L ratio seasonal differences.

Local variation
In terms of potential environmental stressors, the H/L ratio did not vary significantly in relation to the density of males, contrary to our expectations. High social stress is recognized as a possible cost of reproduction (Jakubas et al. 2011). In several avian species nesting in high density areas, individuals are affected by strong social stress and, subsequently, show an increase in the overall H/L ratio induced by exposure to external stressors (i.e. risk of injury, agonistic interactions, etc.). For example, Vleck et al. (2000) found that H/L ratio in Adélie penguins (Pygoscelis adeliae) decreased after hatching, which suggested a previous stressful situation due to intraspecific competition within the colony and the pressure to defend the clutch against predators. In addition, increased population density has been related to higher infection rate by horizontally transmitted parasites and pathogens, and colonial species often show higher parasitic pressure in large groups (Tella 2002;Brown and Brown 2004) that can affect the stress status of individuals (Raouf et al. 2006). Therefore, it is unlikely that our results actually reflect that a higher male density does not imply a higher overall stress. Conversely, it may indicate (i) a bias in the prevalence/pathogenicity of parasites between populations, (ii) that our density estimates are not the most adequate to reflect social stress, or (iii) that birds from dense areas are actively coping with stressors. First, we found that infections affected the H/L ratio, so the potential differences between populations in this causal factor (or other unmeasured confounders, such as local food availability or local vegetation community Fairhurst et al. 2013) might potentially prevent the detection of more subtle relationships between population density and stress. Second, it is posible that mere the number of males in a certain population does not reflect the true population density. Besides, the spatial scale used to estimate population density might not be the most appropriate. Lastly, it is possible that the populations studied did not reach high enough density values as to have a detectable effect on stress levels.
Once we controlled for the potential spatial (location) and temporal (year) dependency of the data, we analysed the potential influence of infection on the H/L ratio. In relation to the potential variation of H/L ratio with internal stressors we showed how the H/L ratio associated to malaria infection differed between countriesIn Spain, our data suggest a link between H/L ratio and the negative covariation between body condition and infection. This result was similar to the one obtained by Becker et al. (2019) in Dark-eyed juncos (Junco hyemalis), where only infected birds with low body condition showed higher values of H/L ratio. Stressful environments (e.g. resource-poor sites or sites with more pathogenic strains) might reveal negative body conditioninfection relationships that are not apparent in more benign environments for hosts. Low values of H/L ratio on individuals with high body condition could be explained by their major resistance against pathogen infections (Jiménez-Peñuela et al. 2019), or the difference in number of parasites needed to cause a significant infection (Marzal et al. 2008). In Morocco, H/L ratio showed no relationship with individual condition. The marginal association with malaria infection could suggest an immune reaction mediated by parasitism (Norte et al. 2009b). Although the different relation of malaria infection,along with body condition,on the H/L ratio is similar to those reported on several bird species (e.g. Cornelius et al. 2014;Granthon and Williams 2017;Skwarska 2018;Becker et al. 2019), results should be interpreted carefully, as the number of infected individuals was quite restricted, which could influence our model performance. Moreover, the potential differences in parasite diversity and pathogenecity within each country (not tested) could also explain these diverging results. In addition, the impossibility of ageing Dupont's Lark and the small number of females trapped prevented testing the effect of age and sex on H/L ratio variation (Norte et al. 2009b;Jakubas et al. 2011;Isaksson et al. 2013;Dantzer et al. 2014).

Conclusions
In this study, we provide the first values of H/L ratio for populations of the threatened Dupont's Lark and assessed the potential impact of different stressors on H/L ratio found in the species. Our results highlight season (breeding vs. non-breeding) as an important source of variation of H/L ratio in the species. In addition, we demonstrate the importance of considering individual-level variables, and their interaction, when studying H/L ratio. Our results are based only on adult males, mainly sampled during the breeding season. Therefore, further research considering not only juveniles and females, but also a more evened data for both seasons within each country and for several years. will contribute to a better understanding of how environmental and individual stressors may affect Dupont's Lark responses to environmental stress (Norte et al. 2009a;Skwarska 2018). Nevertheless, our results might serve as starting point for comparison in future studies aiming to elucidate the health status of threatened Dupont's Lark populations, as well as for future conservation programmes.

Supplementary Information
The online version contains supplementary material available at https:// doi. org/ 10. 1007/ s10336-021-01958-x. whose helpful comments contributed to improve the paper. This work was approved by the Local Ethical Committee for Animal Experiments of the Autónoma University of Madrid (CEI80-1468-A229). All participants accepted to be part of the study.
Funding Open Access funding provided thanks to the CRUE-CSIC agreement with Springer Nature.
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/.