- Research article
- Open Access
Do avian blood parasites influence hypoxia physiology in a high elevation environment?
BMC Ecologyvolume 18, Article number: 15 (2018)
Montane birds which engage in elevational movements have evolved to cope with fluctuations in environmental hypoxia, through changes in physiological parameters associated with blood oxygen-carrying capacity such as haemoglobin concentration (Hb) and haematocrit (Hct). In particular, elevational migrants which winter at low elevations, encounter varying intensities of avian haemosporidian parasites as they traverse heterogeneous environments. Whilst high intensity parasite infections lead to anaemia, one can expect that the ability to cope with haemosporidian infections should be a key trait for elevational migrants that must be balanced against reducing the oxygen-carrying capacity of blood in response to high elevation. In this study, we explored the links between environmental hypoxia, migration, and disease ecology by examining natural variation in infections status and intensity of avian haemoporidians across a suite of Himalayan birds with different migratory strategies while controlling for host phylogeny.
We found predictably large variation in haemoglobin levels across the elevational gradient and this pattern was strongly influenced by season and whether birds are elevational migrants. The overall malaria infection intensity declined with elevation whereas Hb and Hct decreased with increase in parasite intensity, suggesting an important role of malaria parasites on hypoxia stressed birds in high elevation environments.
Our results provide a key insight into how physiological measures and sub-clinical infections might affect dynamics of high-elevation bird populations. We suggest a potential impact of avian elevational migration on disease dynamics and exposure to high intensity infections with disease spread in the face of climate change, which will exacerbate hypoxic stress and negative effects of chronic avian malaria infection on survival and reproductive success in wild birds. Future work on chronic parasite infections must consider parasite intensity, rather than relying on infection status alone.
Environmental hypoxia—the decreased partial pressure of oxygen—and cold temperature are key drivers in the evolution of high elevation adaptations in montane organisms . In particular, birds that undertake elevational migrations are exposed to considerable seasonal fluctuations in hypoxia levels and temperature across the year. Therefore elevational migrants must evolve physiological strategies to regulate haemoglobin (hereafter Hb) concentration, an important measure of blood oxygen-carrying capacity , as they move between elevations. High intensity infections with avian ectoparasites and haemoparasites have been shown to lower the oxygen-carrying capacity of their hosts by reducing Hb levels through direct consumption or destruction of red blood cells, often leading to anaemia . Thus the capacity to regulate Hb concentration and cope with hypoxia can be significantly compromised by haemoparasite infection in elevational migrants which are at increased risk of infection due to the following several factors. Elevational migrants tend to encounter more parasites as they are exposed to high prevalence areas at low elevations with optimal climatic conditions for parasite transmission (e.g. ), have increased physiological demands incurred due to changes in immune conditions during migration ( but see ), have a relapse of chronic infections prior to or during the journey , and have increased risk of exposure due to host aggregation during migration or at the wintering sites . Given the physiological strategies and life-history characteristics, high intensity parasite infections could be especially costly in hypoxic environments for montane birds. Recent evidence suggests that Hb is a more reliable indicator of host response to parasitism than lowered haematocrit (Hct, total red blood cell volume), due to the commonly observed occurrence of erythropoiesis in parasitized hosts . In birds, regenerative or haemolytic anaemia leads to increased erythropoiesis that is rapid production of immature red blood cells stimulated by hypoxia and yet showing normal Hct, despite a substantial parasite load (intensity). To date, only a few studies have compared association between Hb levels and ectoparasite load [9,10,11] and haemoparasites , and these have only been in juvenile altricial birds .
Avian haemosporidians (Apicomplexa: Haemosporida; Plasmodium and other related genera Haemoproteus spp. and Leucocytozoon spp.; henceforth referred to as ‘avian malaria’) are vector-borne blood parasites that place significant selective pressure on their host populations by reducing reproductive success , immune response [14, 15], and chronic-infection driven fitness . Increased parasite intensities can have negative effects on bird populations by reducing growth  and causing higher mortality and/or lower birth rates .
The temporal and spatial dynamics, and epidemiology of avian malaria are strongly governed by ecological (season, habitat quality, elevation) [19, 20], demographical (host and vector density) , and environmental factors (temperature, precipitation) . Seasonal variation in ambient temperature, influences host condition through an increase in energy demands , leading to stress, and thereby increasing susceptibility to infection . In addition, host and vector abundance govern the frequency-dependent transmission of the pathogens [22, 25,26,27,28]. Cosgrove et al.  showed a marked seasonality in Plasmodium transmission in temperate environments, which appears to result in zero prevalence estimates in winters. This is in stark contrast to the tropics where there is a lack of seasonal variation in Plasmodium spp. prevalence in birds despite a strong correlation between mosquito abundance and temperature . Given these seasonal effects, the role of avian haemosporidians on their host has mainly been examined in the breeding season where the spring–summer emergence of parasitemia in avian hosts coincides with a high allocation of energy in reproduction, and high levels of corticosterone in the blood. However, the impact of parasite infection in the non-breeding season has largely remained neglected but see . Furthermore, climate change is likely to facilitate expanding geographic range of Plasmodium species with an increase in the transmission window in habitats that are currently too cool to sustain vector populations and parasite development.
The western Himalayan bird community provides a unique opportunity to understand the interplay between hypoxia and blood parasite intensity. In this montane system, birds exhibit two migration strategies; species are either year-round high elevation residents (sedentary) or seasonal elevational migrants. Elevational migrants winter at low elevations or in the plains (≤ 1500 m above sea level; a.s.l.) and move to breeding grounds at higher elevations (2600–4000 m a.s.l. or even higher) during the summer season . In the wintering grounds, resident birds may act as reservoirs for blood parasites, increasing the risk of infection for migrants. Elevational migrants thus, encounter a more diverse fauna of parasites compared with sedentary species. Barve et al.  showed that high elevation residents regulate Hb by having high mean cell haemoglobin concentration (MCHC), while elevational migrants adjust their Hb through fluctuation in Hct as they move between high and low elevations. The degree to which elevation and blood parasites influence Hb levels across elevational migrants and resident birds has remained unexplored and can have important implications for the ecology of parasites and montane birds. Given that thermal conditions and suitable vectors are present to transmit and maintain the infection, such migrants can form an effective bridge for parasites between the wintering and breeding grounds. Hence, this increases the risk of infection to naïve resident birds at high elevations which might not have evolved to cope with decreased Hb levels as a result of parasite infection. Additionally, elevational migrants are exposed to hypoxia for a relatively short period compared to their high elevation resident counterparts. Therefore, an examination of Hb concentration and haemosporidians infections in resident and migrant species can reveal whether physiological stress and haemolytic anaemia might be more detrimental to high elevation residents, with more haemoglobin lost with the haemolysis of each red blood cell, exacerbating physiological stress across a large elevational gradient.
In this study we examined the infection status and infection intensity of three parasite genera Plasmodium, Haemoproteus and Leucocytozoon in 18 Himalayan avian species which differ in their migration strategies, and their erythropoietic response to hypoxia. We predicted that: (i) seasonal variation in probability of infection and infection intensity will be higher in elevational migrants than their resident counterparts; (ii) infection status and infection intensity will show a negative relationship with Hct and Hb levels; and (iii) high infection intensity and low Hb levels will be associated with elevational migrants as they are likely to experience more chronic hypoxia in the breeding season than high elevation residents.
Study site, avian blood sampling and physiological parameters
The bird blood sampling was conducted in the non-breeding season (January–March) and breeding season (April–May) in Kedarnath Wildlife Division, Uttarakhand, India. We selected 18 species of passerine birds based on their abundance and diversity across seven elevations: 1000, 1500 (low), 2100 (mid), 2650, 2800, 3000 and 3200 m (high) a.s.l. as mentioned in Barve et al.  (Fig. 1a, Additional file 1: Table S1).
At each sampling site, 6–12 mist nets were set up in high bird activity locations. Each bird caught was identified and ringed with a numbered metal ring. Bird blood from the sub-brachial wing vein (never exceeding 1% of the individual’s body weight) was collected in SET Buffer (20–40 μL in 500 μL buffer 0.15 M NaCl, 0.05 M Tris, M EDTA, pH 8.0) or FTA cards (Whatman) for molecular analyses. Captured individuals were released at the site immediately after processing. In addition, for each blood sample, the concentration of Hb was determined using a portable HemoCue Hb 201+ photometer (HemoCue Hb, Ängelholm, Sweden). In this study we measured Hb for 18 bird species and obtained 286 samples in the breeding season and 287 in the non-breeding season (Additional file 2: Table S2 and Fig. S1 for summary of bird species with infection status and Hb values).
Whenever a sufficient sample (20–30 μL) was available, we also measured Hct (n = 245; Additional file 2: Table S2) for the same samples, using a Zipocrit Portable Centrifuge (LW Scientific Inc., Lawrenceville, GA, USA) with a spin time of 5 min. Both Hb and Hct measurements were highly repeatable within this study .
In addition, for the 2014–2015 bird sampling in the breeding season (April–May), thin blood smears were prepared on glass slides to identify whether infections exhibited parasite gametocytes, which would indicate that elevational migrants act as ‘true’ hosts for that particular parasite strain . The smears were air-dried, fixed in 100% methanol and stained with Giemsa.
Infection status and intensity estimation
DNA extractions were conducted using the phenol chloroform extraction method . The presence of Plasmodium, Haemoproteus and Leucocytozoon was assessed using parasite-specific primers designed to amplify partial fragments the cytochrome b (cyt-b) gene using nested PCR protocol . Each plate run included a parasite positive control and also a negative control (water instead of template DNA) to examine for any potential contamination. We screened all parasite negative samples for bird DNA (for cyt-b gene) following Dumbacher et al. . The resulting PCR products were then sequenced in both directions. Sequences were assembled, aligned and edited using SEQUENCHER version 5.2 (Genecodes Corp., Ann Arbor, MI, USA). We then identified sequences to genus using their closest sequence matches in the GenBank or MalAvi databases . For each positive PCR product (without blood smear in non-breeding season; January–March), parasite intensities were calculated as relative quantification values (RQ) as follows. We performed two quantitative PCRs (qPCR): one targeting the 18sRNA gene of avian haemosporidians following Ishtiaq et al. , and the other targeting the 18sRNA gene of the bird (as described in ). RQ can be interpreted as the fold-amount of the target gene (parasite 18s rDNA) with respect to the amount of the reference gene (Bird18s rDNA) and was calculated as 2−(Ct18s parasite−Ct18s Bird) where Ct represents the number of PCR cycles at which fluorescence was first detected as statistically significant above the baseline, which is inversely correlated with the initial amount of DNA in a sample. Each DNA sample was run in duplicate and the average values were used for further analysis. For samples run on different plates we used one reference sample common to all plates. We ran all samples, including the inter run calibrator (IRC), in duplicate. Runs were validated only if the non-template control (NTC) and the negative control did not exhibit fluorescence curves that crossed the threshold line, and the positive control gave a fluorescence curve that crossed the threshold line within 38 cycles (Ct ≤ 38). Following each run, we examined the results to ensure that replicates had similar Ct values and melting peaks. We validated intensity data against microscopy by including positive samples which had corresponding blood smear to calculate intensity. These results showed a good correspondence between parasitemia estimated by qPCR and microscopy  which allowed us to compare intensity results using qPCR and traditional microscopy techniques. Similar results have been reported in other studies (e.g. ). For convenience, RQ values were log transformed into intensity categories so it can be combined with the data with intensities scored on smears (see below). For example, assuming that in one microscopic field is 300–400 erythrocytes, then: ++++ should be more than approx. 3% (in percentage; 3 of 10 erythrocytes are infected); +++ should be about 0.3–1.3%, ++ should be about 0.03–0.2%, + should be about 0.003–0.02%. So, qPCR values (host/parasite ratio) of 1.796519336 were considered as +++; 0.011073119 + 0.064422307 ++. This way we generated intensity data for all samples ranging from submicroscopic, low (+), medium (++), high (+++) to very high (++++) intensity.
Blood smears were examined for the presence of infected blood cells and gametocytes following Godfrey et al. . Briefly, all slides were first examined at low magnification (500×) for approximately 100 fields and then at least 100 fields were studied at high magnification (1000×) using an oil immersion lens (taking into account the 10× contribution from the ocular combined with objectives of the magnification reported). The intensity of infection was recorded as per the following criteria: low (+); 1–10 parasites per 100 thin film fields, medium (++); 11–100 parasites per 100 thin film fields, high (+++); 1–10 parasites in one thin film field and very high (++++); more than 10 parasites in one thin film field.
Modelling disease risk in Hb compromised species
For comparison of infection status across three parasite genera, we carried out four distinct analyses. We classified each sampled species as either an elevational migrant (EM status = 1) where almost all individuals migrate to significantly lower elevations or as sedentary (EM status = 0) were species that are consistently found in their breeding range in the winter in our study site. Information on elevational movement and distribution was obtained from a field survey  in our study site and Rasmussen and Anderton . The first analysis consisted of generalized linear mixed models (GLMMs, function glmer in lme4; Bates et al. ) to assess whether individual infection status with Plasmodium, Haemoproteus and Leucocytozoon was influenced by elevation, bird migration strategy and season, as fixed effects, and bird taxonomy (species nested in genus and genus nested in family) as a random effect. All models were specified with a binomial error distribution and logit link function. The significance of fixed effects was evaluated with Wald’s χ2-tests .
In the second analysis, we used a linear mixed effect models to assess differences in infection intensity of the three parasite genera (Plasmodium, Haemoproteus and Leucocytozoon) with elevation, season and migration status as fixed effects, and bird species as a random effect. For the intensity data, we used only positive samples with microscopy and qPCR methods. We further analysed the relationship between infection status and parasite intensity with Hct sampled across species. Owing to limited data, Hct was not analysed for migration strategies and season. Finally, we sought to tease apart the effects of migration status, parasite intensity, elevation, and season, on Hb levels. In order to investigate factors relating to Hb (at a variety of complexities/scales), we constructed predictive regression models using the linear mixed-effects kinship model fit by maximum likelihood (lmekin) within the “coxme” package . For each model, we controlled for phylogenetic effects by including a phylogenetic correlation matrix as a random effect using the “ape” package . The phylogenetic correlation matrix was derived from a phylogeny constructed from information by Jetz et al. . All data was checked for normality and corrected for over-dispersion as required.
Prior to running predictive models, correlations between predictor variables were assessed to test for collinearity using correlation matrices and simple linear models generated for each pairwise combination of predictors. There was a significant positive relationship between Hb and elevation in both elevational migrants (R2 = 0.24, P < 0.001) and residents (R2 = 0.16, P < 0.001; ). We used intensity as an interactive term in all models since infection intensity can vary with time, seasons, age and sex [20, 47] and showed a negative association with Hb levels (see below). All other pairwise comparisons yielded absolute correlation values and variance explained by linear models (adjusted R2) of < 0.5, indicating no issues of collinearity, so all predictors were retained.
Using data for 15 species in summer and nine species in winter (18 species in total), we fitted 18 linear mixed models separately for Leucocytozoon, Plasmodium, Haemoproteus, mixed infections and total infected birds. We used AIC-based multi-model inference to identify well-supported statistical models that describe the relationships between Hb and biological parameters relevant to elevational distribution, such as (i) infection intensity, (ii) elevation (low, mid, high), (iii) EM status (1 = elevational migrants, 0 = resident species), (iv) season (breeding = 1, nonbreeding = 0), as fixed effects. In all cases, bird species was included as a random effect. The models were written based on variables that were ecologically relevant to hypoxia physiology (Additional file 3: Table S3) and model selection was done using the MuMIn package .
We used the Akaike Information Criteria (AIC; see ), to select the best-fit model. We generated a set of 18 candidate models (including global model and single predictor models) which were ranked by AICc (Additional file 3: Table S3). Models with a difference (ΔAICc) of ≤ 2 are as parsimonious as the best-fit model (lowest AICc), and is considered strong evidence that the quality of the candidate models differ from one another . Relative importance of the traits was assessed and model-averaged estimates were derived to account for model uncertainty . We used a model-averaging approach to check the validity of the top-ranking model in each case, only including models with ΔAICc < 2 . The relative importance (RI) of each parameter after model-averaging, was calculated by summing wi across all models in which the parameter was present. Analyses were conducted in R v. 3.0.1 .
Of the 573 birds sampled, 44.8% (95% CI 40.7–48.9%) were positive for avian haemosporidians infections. Leucocytozoon spp. showed significantly higher infection prevalence (39.7%; 95% CI 35.9–43.9%) than Haemoproteus spp. (10.2%; 95% CI 7.8–12.7%) and Plasmodium spp. (3.1%; 95% CI 1.7–4.5%) (Wald’s χ2 = 232.88, df = 2, P < 0.0001). Only 7.3% (95% CI 5.1–9.4%) were co-infections.
Seasonal variation in the probability of infection and parasite intensity
In the first analysis, across the three parasite genera, there was no significant difference in the infection status with elevation, season and migration status of species (Table 1). However, the probability of Plasmodium infection was high for low elevations, whereas Leucocytozoon infection risk increased with elevation, albeit as submicroscopic infections. In the second analysis the infection intensities of Leucocytozoon (Wald’s χ2 = 15.70.12, df = 1, P < 0.0001) was higher in non-breeding than in the breeding season, but with no effect of bird migration status (Table 2). Total parasite intensity showed a significant decrease with elevation (Wald’s χ2 = 7.12, df = 2, P < 0.01; Fig. 1b) in the breeding season. Birds had a higher infection intensity in the non-breeding season across mid-elevation (Wald’s χ2 = 6.06, df = 2, P < 0.04; Table 2).
Effect of parasite intensity on Hct and Hb
In the third analysis, there was a significant negative relationship between Hct and overall avian malaria intensity (b = − 0.01, t = − 2.45, P < 0.02; Fig. 1c). The probability of infection with Plasmodium spp. showed significant increases with Hct (Wald’s χ2 = 4.54, df = 1, P < 0.03; Fig. 2). However, neither Haemoproteus spp. (Wald’s χ2 = 1.69, df = 1, P = 0.19) nor Leucocytozoon spp. (Wald’s χ2 = 0.11, df = 1, P = 0.73) showed significant variation in infection status with Hct. The total parasite intensity as a single predictor showed a negative relationship with Hb (b = − 0.22, t = − 2.68, P < 0.001).
In the final analysis, for total parasite intensity, model-averaged estimates derived from 90% model set agreed with the best-approximating model with high elevation, elevational migrant and interaction of breeding season and intensity. High Hb levels were influenced by high elevation, elevational migrant and breeding season. Although, for Leucocytozoon infected birds, the model-averaged estimates indicated that high elevation influenced the high Hb levels, however, parasite intensity appeared to have negative effect on Hb. For Haemoproteus infected birds, none of the parameter showed any influence on Hb. However, for Plasmodium infected birds, Hb was positively influenced by elevational migrants. Mixed infections showed the negative influence of parasite intensity on Hb levels as the top model (b = − 0.54, t = − 2.78, P < 0.001; Fig. 1d; Table 3).
This is the first study which links interspecific variation in migratory strategies, variation in erythropoietic responses to hypoxia and parasite infection and intensity across an elevational gradient. Our current understanding of the seasonality and epidemiology of the avian haemosporidian community has been largely derived from studies conducted in the breeding season in temperate regions, when parasite infections peak in spring and summer with an increase in vector activity [52, 53]. Our study on western Himalayan birds shows that parasites persist in winters and infections are higher in intensity in the non-breeding season. We also demonstrate that high parasite intensity has a negative influence on both Hb and Hct suggesting a role of haemosporidian infections in impacting hypoxic stress in high elevation Himalayan birds.
Seasonal variation in parasite intensity and infection status
The presence of parasitemia throughout our study period (January–May; 2014–2015) suggests various underlying mechanisms; (i) the continued presence of all three parasite genera in medium to high intensity during the non-breeding season, highlighting the environmental stress in winters which likely induces a relapse of the previous year’s infections [52, 54]; (ii) scarcity of food in highly seasonal habitats like the Himalayas (e.g. ), which directly affects body condition [56, 57], and is often reflected in high levels of infection intensity in winters. Migrant birds generally have poorer condition in winters compared to summer  and body condition can deteriorate with high intensity infections . In addition, Hb declines with haemosporidian infections , and a lower Hb with high intensity infections mainly driven by Leucocytozoon spp. in the non-breeding season, would reflect a seasonal decline in trophic conditions. High intensity infections in elevational migrants in winter is thus an important finding, which has important implications for seasonal disease dynamics with the potential for carryover effects to the breeding season. Many elevational migrants sampled at high elevations during the breeding season were infected with submicroscopic infections.
Dunn et al.  showed high incidence of Haemoproteus spp. in the non-breeding season, however, at low intensity. In our study, total parasite intensity was high in the non-breeding season. Plasmodium spp. intensity showed no variation across seasons which is likely to be chronic stages of infections, whereas Leucocytozoon spp. showed high intensity in the non-breeding season regardless of bird migration status, which can be classified either as chronic or relapses. These findings have important implications for over-winter ecology of Leucocytozoon spp. in elevational migrants. Nonetheless, the total parasite intensity showed a significant decrease with elevation despite no change in parasite infection status with elevation. Temperature is a key driver of haematophagous arthropod vectors and parasite life-history traits that combine to determine transmission intensity. Even though the probability of infection with three parasite genera at these sites is variable, the sampling sites were selected primarily to illustrate the contrasting thermal environments that can exist across relatively small spatial scales within a region. Our sampling sites vary in elevation (Fig. 1a) and have marked differences in mean temperatures and diurnal temperature ranges, and thus have divergent effects on parasite development and persistence at a population level . The daily temperatures in the low elevation sites during the non-breeding season were warmer (10–25 °C) than mid to high elevation sites (4–10 °C; FI unpub. data) which suggests that low elevations have optimal conditions and lack any thermal constraints for the parasites and this probably allows the parasites to be present in the blood stream during non-breeding season. This would imply an ongoing transmission at the low elevation during winter as long as potential vector species are present, whereas at colder higher elevation sites, parasites withdraw from the blood stream until the ideal conditions appear.
Effect of parasite intensity on Hct and Hb
In western Himalayan birds, high Hb levels were strongly associated high elevation, and the breeding season where elevational migrants are likely to experience chronic hypoxia (potentially compounded by increased energetic demands) . In our linear mixed model analyses to examine factors influencing Hb levels in parasite-specific infected bird models, parasite intensity was selected as an important predictor in mixed infection birds, however, parasite intensity showed negative significant effect in Leucocytozoon infected birds with confidence intervals for the parameter estimates included zero, so there was little evidence that Leucocytozoon parasite had any effect on Hb levels. These results are consistent with other studies that report negative effects of parasite intensity on haemoglobin concentration in birds [3, 18, 59]; whereas some studies report no effect of parasite intensity on Hb levels . Based on a recent experimental study, even low intensity chronic infection influences physiology and fitness of hosts . Therefore, measurement of parasite intensity is important rather than relying only on infection status in evaluating the role of infections on host fitness and physiology.
Whilst the prevalence estimates of haemosporidians are snapshots in time and space [62, 63], our findings are important from physiological as well as ecological perspectives. Previous studies have shown that there is a correlation between parasite intensity and the severity of human malarial infections-high parasite intensity is correlated with high mortality rates . It is quite possible that high intensity blood parasite infections are pathogenic to elevational migrants and lead to excessive haemolysis of parasitized erythrocytes in malaria infection, and may lead to anaemia (e.g. [65, 66]). This is further supported by the low frequency of very high intensity parasite infections in high elevation sites. Many elevational migrants were detected with submicroscopic infection (detected by PCR but not by microscopy and this implies absence of infective parasite stage—gametocytes) in high elevation environments in peak breeding season despite being found in medium or high intensity at low elevations. These findings are in stark contrast with other studies in temperate regions where the majority of the birds with latent infections (birds infected during the prior breeding season) return to the breeding grounds and experience an increase in parasite blood stages—relapse [54, 67]. Parasite intensity in Plasmodium falciparum often shows a clear trend of decreasing with decreasing transmission intensity  and our finding implies that elevation can be a proxy for transmission intensity in avian haemosporidians, highlighting that studies should not only rely on prevalence based models as intensity drives the transmission intensity in high elevation sites and the combination of molecular tools and microscopy is critical in assessing parasite load. While we quantified relative intensity across Himalayan birds, one of the limitations of using mix intensity data—by converting qPCR values into intensity categories as scored on smears—meant loss of quantitative information which would have provided better resolution. Our estimates could be compromised for genus level intensity in few co-infected birds as qPCR primers amplify all three parasite genera successfully . Nonetheless, our parasite intensity estimates from qPCR and microscopy techniques were comparable as we screened all samples using nested PCR protocols [38, 39].
Similarly, Hct, an important predictor of blood oxygen transport in elevational migrants, showed a negative effect with overall parasite intensity suggesting that the infections with haemosporidians lowers haematocrit. A recent field study of wild red-winged blackbirds population demonstrated that birds with higher Plasmodium and Haemoproteus parasitemia have lower hematocrit and a higher rate of red blood cell production . The probability of Plasmodium infection showed a slight increase with Hct, probably as a result of upregulation of erythrocyte production, although, this warrants further investigation. Plasmodium remained absent in high elevation environments in resident birds and showed distinctly low prevalence in migrant populations. Given that all infections in our study are identified at the parasite genus level and we have not been able to detect differences in pathogenicity within a genus, the pathogenicity could be related to specific parasite lineages which could in turn have a variable effect . In general, Haemoproteus has lower pathogenicity on its host than Plasmodium [70,71,72]. It is possible that Plasmodium infections are pathogenic to high elevation residents, which is reflected by low incidence at high elevations (but see below). Finally, some host species may be more susceptible to infection than others, and the species chosen for this study may have co-evolved with Leucocytozoon spp. parasites becoming resistant to them , a factor that probably explains the generally high prevalence of this parasite genus in western Himalayan birds. The high prevalence of Leucocytozoon spp. throughout the elevational gradient is in contrast with other studies (e.g. ) where Leucocytozoon showed an association with high elevation rather than in adjacent lowlands. In our western Himalayan system, there was no difference in infection status but intensity showed a sharp decline across the elevational gradient which points towards the role of high intensity infections lowering Hct and Hb levels at high elevations.
From an ecological perspective, climate change increases the risk of malaria in many birds . In addition to rising mean temperatures, diurnal fluctuations in temperature have also been shown to affect the rate of parasite development  and range expansion, or allows a longer breeding season for arthropod vectors, both of which can alter transmission rates of malaria parasites. The relatively low prevalence of Plasmodium in high elevation environments could be attributed to micro-climatic differences governing the vector populations between low and high elevations and therefore influencing the spread of the disease across elevational gradients. The high elevation Himalayan bird community may not generally be exposed to vectors and parasites throughout the year—Plasmodium transmission is restricted by thermal constraints and distribution of potential vectors , and Plasmodium spp. was not found in the resident birds in the high elevation sites and was limited to the elevational migrants. With rising temperature predicted by imminent climate warming , and the subsequent range expansion of mosquito species and Plasmodium spp. parasites to high elevations, this might pose significant risks to naïve bird populations (e.g. Hawaii islands; ).
Our study explored the links between seasonal migration, hypoxia physiology and disease dynamics. We found a higher parasite intensity in the non-breeding season regardless of migratory strategy, which is in line with previous studies suggesting that winter induces the relapse of the previous year’s infections. Our results show that high intensity in mixed haemosporidian infections have the predicted negative effect on oxygen-carrying capacity of blood after controlling for host species. Future work on chronic parasite infection must consider parasite intensity, rather than only relying on infection status. These results imply a role of elevational migration in disease dynamics. The potentially increased exposure to high intensity infections with upslope disease spread in the face of climate change will exacerbate hypoxic stress experienced by high elevation birds.
above sea level
generalized linear mixed models
polymerase chain reaction
Ramirez J-M, Folkow LP, Blix AS. Hypoxia tolerance in mammals and birds: from the wilderness to the clinic. Annu Rev Physiol. 2007;69:113–43.
Campbell TW. Avian hematology and cytology. Ames: Iowa State University Press; 1995.
O’Brien EL, Morrison BL, Johnson LS. Assessing the effects of haematophagous ectoparasites on the health of nestling birds: haematocrit vs haemoglobin levels in house wrens parasitized by blow fly larvae. J Avian Biol. 2001;32:73–6.
Waldenström J, Bensch S, Kiboi S, Hasselquist D, Ottosson U. Cross-species infection of blood parasites between resident and migratory songbirds in Africa. Mol Ecol. 2002;11:1545–54.
Owen J, Moore F. Seasonal differences in immunological condition of three species of thrushes. Condor. 2006;108:389–98.
Hahn S, Buer S, Dimitrov D, Emmenegger T, Ivanova K, Zehtindjiev P, Buttemer WA. Low intensity blood parasite infections do not reduce the aerobic performance of migratory birds. Proc R Soc B. 2017;285:2307.
Gylfe A, Bergström S, Lundström J, Olsen B. Reactivation of borrelia infection in birds. Nature. 2000;403:724–5.
Krauss S, Stallknecht DE, Negovetich NJ, Niles LJ, Webby RJ, Webster RG. Coincident ruddy turnstone migration and horseshoe crab spawning creates an ecological “hot spot” for influenza viruses. Proc R Sci Lond B. 2010;277:3373–9.
Powlesland RG. Effects of haematophagous mite Ornithonyssus bursa on nestling starlings in New Zealand. NZ J Zool. 1977;4:85–94.
Chapman BR, George JE. The effects of ectoparasites on cliff swallow growth and survival. In: Loye JE, Zuk M, editors. Bird–parasite interactions. Ecology, evolution and behaviour. Oxford: Oxford University Press; 1991. p. 69–92.
Whitworth TW, Bennett GF. Pathogenicity of larval Protocalliphora (Diptera: Calliphoridae) parasitizing nestling birds. Can J Zool. 1992;70:2184–91.
Krams IA, Suraka V, Rantala MJ, Sepp T, Mierauskas P, Vrublevska J, Krama T. Acute infection of avian malaria impairs concentration of haemoglobin and survival in juvenile altricial birds. J Zool. 2013;291:34–41.
Marzal A, de Lope F, Navarro C, Møller AP. Malarial parasites decrease reproductive success: an experimental study in a passerine bird. Oecologia. 2005;142:541–5.
Merino S, Moreno J, Sanz JJ, Arriero E. Are avian blood parasites pathogenic in the wild? A medication experiment in blue tits (Parus caeruleus). Proc R Sci Lond B. 2000;267:2507–10.
Tomás G, Merino S, Moreno J, Morales J, Martínez-de la Puente J. Impact of blood parasites on immunoglobulin level and parental effort: a medication field experiment on a wild passerine. Funct Ecol. 2007;21:125–33.
Asghar M, Hasselquist D, Hansson D, Zehtindjiev P, Westerdahl H, Bensch S. Hidden costs of infection: chronic malaria accelerates telomere degradation and senescence in wild birds. Science. 2015;347:436–8.
Valkiūnas G, Zickus T, Shapoval AP, Iezhova TA. Effect of Haemoproteus belopolskyi (Haemosporida: Haemoproteidae) on body mass of the blackcap Sylvia atricapilla. J Parasitol. 2006;92:1123–5.
Dawson RD, Bortolotti GR. Effects of hematozoan parasites on condition and return rates of American Kestrels. Auk. 2000;117:373–80.
Rooyen J, Lalubin F, Glaziot O, Christe P. Altitudinal variation in haemosporidian parasite distribution in great tit populations. Parasit Vectors. 2013;6:139.
Wood MJ, Cosgrove CL, Wilkin TA, Knowles SC, Day KP, Sheldon BC. Within-population variation in prevalence and lineage distribution of avian malaria in blue tits, Cyanistes caeruleus. Mol Ecol. 2007;16:3263–73.
Lachish S, Knowles SC, Alves R, Wood MJ, Sheldon BC. Infection dynamics of endemic malaria in a wild bird population: parasite species-dependent drivers of spatial and temporal variation in transmission rates. J Anim Ecol. 2011;80:1207–16.
LaPointe DA, Goff ML, Atkinson CT. Comparative susceptibility of introduced forest-dwelling mosquitoes in Hawai’i to avian malaria, Plasmodium relictum. J Parasitol. 2005;91:843–9.
Hasselquist D, Nilsson J-A. Physiological mechanisms mediating costs of immune responses: what can we learn from studies of birds? Anim Behav. 2012;83:1303–12.
Sheldon BC, Verhulst S. Ecological immunology: costly parasite defences and trade-offs in evolutionary ecology. TREE. 1996;11:317–21.
Hay SI, Cox J, Rogers DJ, Randolph SE, Stern DI, Shanks GD, et al. Climate change and the resurgence of malaria in the East African highlands. Nature. 2002;415:905–9.
Okanga S, Cumming GS, Hockney PAR. Avian malaria prevalence and mosquito abundance in the western Cape, South Africa. Malaria J. 2013;12:370.
Medeiros MC, Ricklefs RE, Brawn JD, Hamer GL. Plasmodium prevalence across avian host species is positively associated with exposure to mosquito vectors. Parasitology. 2015;142:1612–20.
Ishtiaq F, Bowden CGR, Jhala YV. Seasonal dynamics in mosquito abundance and temperature do not influence avian malaria prevalence in the Himalayan foothills. Ecol Evol. 2017;7:8040–57.
Cosgrove CL, Wood MJ, Sheldon BC. Seasonal variation in Plasmodium prevalence in a population of blue tits Cyanistes caeruleus. J Anim Ecol. 2008;77:540–8.
Dunn J, Goodman SJ, Benton TG, Hamer KC. Active blood parasite infection during the non-breeding season: an overlooked issue in declining populations? BMC Ecol. 2013;13:30.
Dixit S, Joshi V, Barve S. Bird diversity of the Amrutganga Valley, Kedarnath, Uttarakhand, India with an emphasis on the elevational distribution of species. Check List. 1874;2016:12.
Barve S, Dhondt AA, Mathur VB, Ishtiaq F, Cheviron ZA. Life history characteristics influence physiological strategies to cope with hypoxia in Himalayan birds. Proc R Sci Lond B. 2016;283:20162201.
Valkiūnas G, Iezhova TA, Loiseau C, Sehgal RN. Nested cytochrome b polymerase chain reaction diagnostics detect sporozoites of hemosporidian parasites in peripheral blood of naturally infected birds. J Parasitol. 2009;95:1512–5.
Sambrook J, Fritsch EF, Maniatis T. Molecular cloning, a laboratory manual. 2nd ed. Cold Spring Harbor: Cold Spring Harbor Laboratory Press; 1989.
Hellgren O, Waldenström J, Bensch S. A new PCR assay for simultaneous studies of Leucocytozoon, Plasmodium, and Haemoproteus from avian blood. J Parasitol. 2004;90:797–802.
Dumbacher JP, Pratt TK, Fleischer RC. Phylogeny of the owlet-nightjars (Aves: Aegothelidae) based on mitochondrial DNA sequence. Mol Phyl Evol. 2003;29:540–9.
Bensch S, Hellgen O, Perez-Tris J. MalAvi: a public database of malaria parasites and related haemosporidians in avian hosts based on mitochondrial cytochrome b lineages. Mol Ecol Resour. 2009;9:1353–8. http://mbioserv4.mbioekol.lu.se/avianmalaria/index.html.
Ishtiaq F, Rao M, Huang X, Bensch S. Estimating prevalence of avian haemosporidians in natural populations—a comparative study on screening protocols. Parasit Vectors. 2017. https://doi.org/10.1186/s13071-017-2066-z.
Asghar M, Hasselquist D, Bensch S. Are chronic avian haemosporidian infections costly in wild birds? J Avian Biol. 2011;42:530–7.
Godfrey RD, Fedynich AM, Pence DB. Quantification of haematozoa in blood smears. J Wildl Dis. 1987;23:558–65.
Rasmussen PC, Anderton JC. Birds of South Asia. The ripley guide. Washington and Barcelona: Smithsonian Institution and Lynx Edicions; 2005.
Bates D, Maechler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. J Stat Softw. 2015;67:1–48.
Bolker BM, Brooks ME, Clark CJ, Geange SW, Poulsen JR, Stevens MH, White JS. Generalized linear mixed models: a practical guide for ecology and evolution. TREE. 2009;24:127–35.
Therneau T. Coxme: mixed effects Cox models. 2012. R package version 2.2-3. https://cran.r-project.org/web/packages/coxme/coxme.pdf. Accessed 27 Feb 2018.
Paradis E, Claude J, Strimmer K. APE: analyses of phylogenetics and evolution in R language. Bioinformatics. 2004;20:289–90.
Jetz W, Thomas GH, Joy JB, Hartmann K, Mooers AO. The global diversity of birds in space and time. Nature. 2012;491:444–8.
Bensch S, Waldenström J, Jonzén N, Westerdahl H, Hansson B, Sejberg D, Hasselquist D. Temporal dynamics and diversity of avian malaria parasites in a single host species. J Anim Ecol. 2007;76:112–22.
Barton K. MuMIn: multi-model inference. R package version 1.0. 0. 2012. http://cran.r-project.org/package=MuMin. Accessed 27 Feb 2018.
Burnham KP, Anderson DR, Huyvaert KP. AIC model selection and multimodel inference in behavioral ecology: some background, observations, and comparisons. Behav Ecol Sociobiol. 2011;65:23–35.
Burnham KP, Anderson DR. Model selection and multimodel inference: a practical information—theoretic approach. 2nd ed. New York: Springer; 2002.
R Development Core Team. R: a language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2015.
Beaudoin RL, Applegate JE, David DE, McLean RG. A model for the ecology of avian malaria. J Wildl Dis. 1971;7:5–13.
Cranston PS, Ramsdale CD, Snow KP, White GB. Adults, larvae and pupae of British mosquitoes (Culicidae). Ambleside: Freshwater Biological Association; 1987.
Applegate JE, Beaudoin RL. Mechanism of spring relapse avian malaria: effect of gonadotrophin and corticosterone. J Wildl Dis. 1970;6:443–7.
Sherry TW, Johnson MD, Strong AM. Does winter food limit populations of migratory birds? In: Greenberg R, Marra PP, editors. Birds of two worlds the ecology and evolution of migration. Baltimore: Johns Hopkins University Press; 2005. p. 414–25.
Strong AM, Sherry TW. Habitat-specific effects of food abundance on the condition of ovenbirds wintering in Jamaica. J Anim Ecol. 2000;69:883–95.
Brown DR, Sherry TW. Food supply controls the body condition of a migrant bird wintering in the tropics. Oecologia. 2006;149:22–32.
Blanford JI, Blanford S, Crane RG, Mann ME, Paaijmans KP, Schreiber KV, Thomas MB. Implications of temperature variations for malaria parasite development across Africa. Sci Rep. 2013;3:1300.
Dudaniec R, Kleindorfer S, Fessl B. Effects of the introduced ectoparasite Philornis downsi on haemoglobin level and nestling survival in Darwin’s small groundfinch (Geospiza fuliginosa). Austral Ecol. 2006;31:88–94.
Colombelli-Négrel D, Kleindorfer S. In superb fairy wrens (Malurus cyaneus), nuptial males have more blood parasites and higher haemoglobin concentration than eclipsed males. J Zool. 2008;56:117–21.
Schoenle LA, Kernbach M, Haussmann MF, Bonier F, Moore IT. An experimental test of the physiological consequences of avian malaria infection. J Anim Ecol. 2017;86:1483–96.
Beadell JS, Gering E, Austin J, Dumbacher JP, Peirce MA, Pratt TK, et al. Prevalence and differential host-specificity of two avian blood parasite genera in the Australo–Papuan region. Mol Ecol. 2004;13:3829–44.
Loiseau C, Harrigan RJ, Bichet C, Julliard R, Garnier S, Lendvai AZ, et al. Predictions of avian Plasmodium expansion under climate change. Sci Rep. 2013;3:1126.
Tangpukdee N, Krudsood S, Kano S, Wilairatana P. Falciparum malaria parasitemia index for predicting severe malaria. Int J Lab Hematol. 2012;34:320–7.
Koncan RM, Clark DT. Anemia in ducks infected with Leucocytozoon simondi. J Protozool. 1966;13:465–8.
Maley GJM, Desser SS. Anemia in Leucocytozoon simondii infections. I. Quantification of anemia, gametocytemia, and osmotic fragility of erythrocytes in naturally infected Pekin ducklings. Can J Zool. 1977;55:355–8.
Valkiũnas G. Avian malaria parasites and other haemosporidia. Boca Raton: CRC Press; 2005.
Okell LC, Bousema T, Griffin JT, Ouédrago AL, Ghani AC, Drakeley CJ. Factors determining the occurrence of submicroscopic malaria infections and their relevance for control. Nat Commun. 2012;3:1237.
Bennett GF, Peirce MA, Ashford RW. Avian haematozoa: mortality and pathogenicity. J Nat Hist. 1993;27:993–1001.
van Riper IIIC, van Riper SG, Goff ML, Laird M. The epizootiology and ecological significance of malaria in Hawaiian land birds. Ecol Mon. 1986;56:327–44.
Atkinson CT, Dusek RJ, Woods KL, Iko WM. Pathogenicity of avian malaria in experimentally-infected Hawaii Amakihi. J Wildl Dis. 2000;36:197–204.
Belo NO, Rodríguez-Ferraro A, Braga EM, Ricklefs RE. Diversity of avian haemosporidians in arid zones of northern Venezuela. Parasitology. 2012;139:1021–8.
Palinauskas V, Valkiūnas GN, Bolshakov CV, Bensch S. Plasmodium relictum (lineage P-SGS1): effects on experimentally infected passerine birds. Exp Parasitol. 2008;120:372–80.
Galen SC, Witt CC. Diverse avian malaria and other haemosporidian parasites in Andean house wrens: evidence for regional co-diversification by host-switching. J Avian Biol. 2014;45:374–86.
Garamszegi LZ. Climate change increases the risk of malaria in birds. Glob Change Biol. 2011;17:1751–9.
Paaijmans KP, Blanford S, Bell AS, Blanford JI, Read AF, Thomas MB. Influence of climate on malaria transmission depends on daily temperature variation. Proc Nat Acad Sci USA. 2010;107:15135–9.
Shrestha UB, Gautam S, Bawa KS. Widespread climate change in the himalayas and associated changes in local ecosystems. PLoS ONE. 2012;7(5):e36741.
del Hoyo J, Elliot A, Sartagal J, Christie DA, De Juana E. Handbook of the birds of the world alive. Barcelona: Lynx Edicions. 2014. http://www.hbw.com. Accessed 5 Mar 2018.
FI conceived the idea. FI, and SB designed the experiment; FI and SB conducted the fieldwork. FI lead the laboratory experiments; FI analysed the data and wrote the manuscript, and SB reviewed the manuscript. Both authors read and approved the final manuscript.
This study was funded by Wellcome Trust/DBT India Alliance Fellowship (IA/I(S)/12/2/500629) to FI, and the Athena Fund at the Cornell Lab of Ornithology and the Explorer’s Club Exploration Fund to SB. The field experiments comply with the current laws of the India where the study was performed. We thank Uttarakhand Forest Department for ethical approval and permission for the collection of avian blood samples. We thank numerous research assistants who helped us during the fieldwork. We thank H. Maithani, M. Rao, P. Jeganathan and C. Bowden for help with fieldwork and K. Narayanan, P. Dey for help with lab work. We also thank Dr. Staffan Bensch, the editor and two anonymous reviewers for constructive comments that improved the manuscript.
The authors declare that they have no competing interests.
Consent for publication
Attached in additional material.
Ethics approval and consent to participate
The field experiments comply with the current laws of the India where the study was performed. We thank Uttarakhand Forest Department for ethical approval and permission for the collection of avian blood samples.
This study was funded by Wellcome Trust/DBT India Alliance Fellowship (IA/I(S)/12/2/500629) to FI, and the Athena Fund at the Cornell Lab of Ornithology and the Explorer’s Club Exploration Fund to SB.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.