Effects of forest structure on the interaction between avian hosts, dipteran vectors and haemosporidian parasites
BMC Ecology volume 20, Article number: 47 (2020)
Forest habitats are important biodiversity refuges for a wide variety of bird species. Parasitism may modulate host species presence and abundance, and parasite effects can change according to forest management practices. Such processes are not well studied in vector-borne avian haemosporidians. We analyzed the effects of forest management on bird-dipteran-haemosporidian interactions, using seven common bird species in managed and unmanaged beech forest habitats in northeastern Germany. We assumed that forest structural heterogeneity affects parasite population parameters in avian hosts (i.e., prevalence and parasitemia), through its effect on the condition of the avian host but also through varying vector abundances.
Parasite prevalence was high (about 80%) and homogeneous across different beech forest categories (i.e., young, old, unmanaged) and for all bird species, except Erithacus rubecula (35%). Parasitemia varied across bird species but not across forest categories within each avian species (lowest parasitemia were found in E. rubecula, Turdus merula, and Turdus philomelos). In our study system, we found that vector abundance was not the main driver of parasite dynamics. We found that forest structure affects parasite infection probability directly and potentially host condition via available resources that have to be used either to combat infections (i.e., high parasitemia) or to maintain a good body condition.
The effects of each of the predictors were bird species-specific, and we found that Diptera vectors were not the foremost influence in our host-vector-parasite system. Effects of forest habitat variables indicated that for most bird species in this study, habitat regulation of infection probability was more likely (i.e., E. rubecula, Fringilla coelebs, Sylvia atricapilla), whereas for Parus major habitat characteristics impacted first individuals' body condition and subsequently the probability of infection. Our findings emphasize the need of species-specific analyses and to use continuous forest structural parameters (e.g., the proportion of gap, south facing aspect) to better understand habitat and land use effects on host-vector-parasite dynamics.
Parasites are an essential component of natural systems and ecosystem processes through their effects on host individuals and thus populations . Parasites typically have negative effects on hosts and may also affect hosts' inter- and intraspecific interactions . Host physiological changes due to an infection can potentially change host behavior , fitness, reproductive success, or survival . Chronic infections produce recurring disease symptoms that may adversely affect hosts throughout their life , regulating host populations and modifying interactions within communities . Bottom-up and top-down processes are likely in host-parasite systems, specifically when the environment either influences first the parasites (top-down regulation) or the host physical condition (bottom-up regulation) . So far, bottom-up and top-down processes have been suggested  but not demonstrated for host-vector-parasite interactions in a varying environment. Therefore, our understanding of the functioning of natural communities will increase by studying the relative importance of both regulatory forces within the same community.
Three-quarters of all bird species use forest habitats . Increasingly, habitats are subjected to human alteration, resulting in habitat degradation, fragmentation and diversity loss through land use intensification . Forest modifications are particularly driven by forestry and management regimes , where local species richness is affected through forest structural heterogeneity and/or forest type . Forest habitats provide resources for avian hosts, such as nesting sites, shelter, foraging opportunities and food . However, not all forest types are optimal for all bird species, and there are species-specific preferences of habitat resources .
Habitat conditions and resource availability may have direct effects on host fitness . The effects on fitness are reflected in body growth and condition , and also in the birds' ability to invest in immune responses to fight-off pathogens (e.g. ). Outcomes of this trade-off can be observed in bilateral asymmetries, such as that of tarsus growth, which reflect early developmental challenges that depend on available resources of the place the young birds inhabit . The life history and habitat preferences of the avian hosts play an essential role in determining the likelihood of an infection [6, 18]. Hence, habitat change alters ecosystem processes and species interactions, leading for example to higher infection rates in more intensively used habitats [19, 20].
Haemosporidian parasites are intracellular parasites that affect a wide range of vertebrate species. These parasites are transmitted by dipterans of different families in which the sexual reproduction takes place. Vertebrate hosts are then required for the asexual phase of the parasite's reproduction . The effects of an infection vary greatly among host species and can lead to chronic lifelong infections (e.g. ). Many effects on individuals are non-lethal, such as weight loss and anemia , decreased reproduction , changes in body asymmetries , or inhibited feather growth . In addition to the fitness effects, parasite infections acting in synergy with anthropogenic modifications of the environment can cause increased infection risk and changes in local abundance of hosts and vectors (e.g. [25, 26]).
Studies incorporating all three components in the interaction between blood parasites, avian hosts and their insect vectors are scarce (e.g. ). Similarly to the avian hosts, vectors may be affected by structural differences in forest habitats or forest management . This ultimately affects vector distribution and abundance, subsequently influencing exposure to parasites and transmission risk [28, 29], changing the spatial variation of prevalence in hosts . Forest habitat structure and management may therefore affect host-vector-parasite interactions in a predictable way, as a function of both avian and insect abundances. In a previous study, the forest habitat was included explicitly as forest categories (i.e., forest stands dominated by one tree species, such as beech or spruce) and forest structure (i.e., continuous variables assessing the interior structural heterogeneity, such as amount of gaps, canopy heterogeneity and canopy height), but the vectors have not been explicitly included in the models so far .
The original hypotheses (c.f. ) can be interpreted as a 'top-down'-hypothesis and a 'bottom-up'-hypothesis. In the top-down-hypothesis, forest structure (as a proxy for management) a priori determines bird infection risk due to the habitat suitability for vectors, in turn triggering an immune response in birds, which then consumes resources fighting-off the infection that would otherwise be available for body growth (measured by asymmetries in tarsus and feather growth). We regard this hypothesis as top-down regulation from the perspective of the avian host, since the parasite effects determine the bird's condition. In the bottom-up-hypothesis, development of the avian host is primarily affected by the forest structure through the availability of resources in the given habitat . Suboptimal habitats provide fewer resources to invest in immune response and defending against parasite infections; the host needs to invest more energy in foraging to compensate for the detrimental effects of habitat quality. Here, there is bottom-up regulation from the perspective of the avian hosts, since the hosts' body condition determines infection risk. In this study, we extend the original hypotheses in two ways:
First, we extend the number of avian host species from two to seven species, which will improve our understanding of the relationship between host traits on parasitism (e.g., the prevalence of avian haemosporidians has been shown to be positively related to the abundance of hosts species; ).
Second, because vectors are known to be affected by forest habitat structure , we explicitly add dipteran vectors as an important mechanistic driver for the variation and spread of avian haemosporidians [28, 33]. Based on the importance of vector abundance for the transmission of avian haemosporidians , we finally hypothesize that the amended models would better explain the host-vector-parasite interactions, compared to the previous study in which the models only included bird hosts and forest structure .
We approach the analysis of our data using two different statistical methods: through model selection and Structural Equation Models (SEMs). The usage of SEMs is novel approach in the field of host-vector-parasite interaction [19, 34], and allows the analysis of multiple pathways and associations among variables, thus offering a way to improve our understanding of host-vector-parasite interactions. We explain our statistical approach in more detail in “Methods” section.
We captured 498 individuals of seven avian host species over two years. Parasite status was analyzed for 441 individuals using microscopy, and 403 individuals through Polymerase Chain Reaction (PCR). Either method resulted in a parasite prevalence of 65% (respectively 288 out of 441 and 260 out of 403 individuals). For 398 individuals we could determine the status through both methods, resulting in 80% infection prevalence over the seven avian host species (Table 1). In 71% of all infected individuals both methods yielded a positive result for the same individual, indicating 29% discrepancy between microscopy and PCR methods. The total prevalence over all avian host species is higher using a combination of both microscopy and the PCR methods (80%).
E. rubecula had a significantly lower mean prevalence than the other avian species (35%; Table 1; Fig. 1a,). We did not find an effect of forest category on parasite prevalence for each avian species separately (Fig. 1a). Median parasitemia varied significantly between avian host species (Mood's median, χ2 = 85.68, df = 6, P < 0.01; Fig. 1b). We did not find significant differences in relation to parasitemia within an avian host species in relation to the forest categories (Fig. 1b). Our model selection approach did not show any effect of forest category on either parasite prevalence or parasitemia; instead, both asymmetry parameters were important factors to be included in a model explaining haemosporidian prevalence (Table 2). The analysis showed that tarsus asymmetry was negatively related to parasitemia and prevalence (Table 3, Additional file 1: Fig. S1, S2).
We captured 4623 female vectors and the dataset included 15 species of Ceratopogonidae (biting midges), 12 species of Culicidae (mosquitoes) and 2 species of Simuliidae. The most observed species in terms of female abundance per vector family were: Aedes cantans (Culidiae), Culicoides impunctatus (Ceratopogonidae) and Simulium lundstromi (Simuliidae). The highest number of species for families Culicidae and Ceratopogonidae were observed in the old age forest category. We only caught black fly individuals in the age class young category (Additional file 1: Table S1). The model selection analysis showed that adding female vector abundance did not improve model fit in relation to parasite prevalence and parasitemia (Table 2).
Accordingly, best-fitted SEMs (Structural Equation Models) indicated that female vector abundance has limited relevance in explaining parasite dynamics for our bird-dipteran-haemosporidian systems (Fig. 2; Additional file 1: Tables S2–S5). Incorporating the measure of vector abundance did not increase model fit, nor did model fit increase when we applied parasite genera and female vector specific relationships (i.e., Ceratopogonidae-Haemoproteus, Culicidae-Plasmodium and Simuliidae-Leucocytozoon (see the test statistics for the '1b' and '2b' models in Additional file 1: Tables S2–S5).
The overall SEM models, in which we analyzed all avian species together but with either primary 3 or tarsus asymmetry, conform to the idea that habitat management and structure (e.g., south facing slope, proportion of gaps) directly affect parasite infection probability and parasitemia first (i.e., top-down regulation; Fig. 2; Additional file 1: Fig. S2a).
However, there were differences in the outcome for each bird species. In the case of P. major and T. merula, forest structural variables first affected bird condition (i.e., bird asymmetries), which subsequently affected the probability of acquiring an infection as well as the severity of the infection (i.e., bottom-up regulation; Additional file 1: Tables S2–S5; Figs. S3, S4). In contrast, S. atricapilla, E. rubecula and F. coelebs were fitted better by the general pattern, where habitat structure directly affects parasite prevalence and parasitemia first (i.e., top-down regulation; Additional file 1: Tables S2–S5; Figs. S3, S4). Finally, in both overall SEM models, the correlation between parasite prevalence and parasitemia was significantly positive. (Fig. 2; Additional file 1: Tables S2–S5).
We analyzed how host-vector-parasite interactions respond to local forest management (GLMM) and forest structure (SEM). The general pattern indicated (i.e., the overall models) parasite regulation (i.e., top-down) from the perspective of the avian hosts. Habitat structure (e.g., south facing slope, proportion of gaps, entropy) first affects the probability of acquiring an infection (Fig. 2, Additional file 1: Figs. S3, S4), which then develops into a severe disease for the individual (measured by parasitemia), and subsequently forces birds to use limited resources to fight-off infections instead of using these resources for maintenance and physical development. For the birds, it implies a negative indirect effect on avian development, given that the birds would have to use resources to fight off infections (i.e., via immune responses; increased leucocyte numbers) instead of using them for growth and maintenance. We partially identified that the proportion of south facing forest (aspect) and the proportion of gaps affect avian haemosporidian transmission dynamics, but effects can vary across different bird species. Parasite transmission dynamics therefore occurred in different ways among different avian species (e.g. ). Importantly, vector abundance did not have an important effect on the transmission of avian haemosporidians in our study system [see also 20].
In one of our avian host species (S. atricapilla) both parasite prevalence and parasitemia were related to increased tarsus asymmetry, which might be an indirect effect of parasite infection. Haemosporidians are reported to affect tail feather growth [35, 36], and similarly other studies report increased asymmetry as a result of ectoparasite infections such as mites and cimicids, affecting the length of feathers, wings, and tarsus . Body asymmetries can be a result of increased stress due to changes in habitat . Asymmetries can be caused by an increased susceptibility and exposure to parasites  and negatively affect host fitness . We found different infection dynamics for different bird species; for some species habitat structure seems to affect first their condition (bottom-up regulation), and for other species forest structure affects first the probability of an haemosporidian infection (top-down regulation). In order to tease apart if parasite infection leads to suboptimal bird conditions or vice versa, an experimental approach would be useful. Most likely, both scenarios, or mechanisms, occur in nature and will vary as a function of interacting species.
The method of parasite detection also played an important role; the total prevalence is higher when using a combination of both microscopy and molecular methods. Total parasite prevalence would therefore have been underestimated using only one method . We also found differences in prevalence and parasitemia between avian species; thus, grouping parasite parameters over all avian host species would disregard subtle effects on actual infections. Our study confirms that parasite prevalence is variable across bird species [e.g. 18, 41]. Grouping host species may lead to a lower overall prevalence, such as when combining infections of 27 European species that lead to a mean prevalence of 26% . Thus, the approach to better understand host-vector-parasite systems has to follow a reductionist framework, studying the same system across different geographical locations, gradually increasing the complexity (i.e., adding functional traits and other interacting organisms) in such a way that both contingent and general outcomes can be easily identified .
Contrary to our predictions, we did not find an effect of female vector abundance on haemosporidian parasitological parameters in any of our avian host species. We hypothesized that forest structure should indirectly affect parasite infection risk because of its known effects on vector abundance [16, 19]. Variation of parasite prevalence in birds has been observed as within-year seasonal variation, reflected by the seasonal variation of vector abundance, distribution or exposure . To illustrate how the presence of vectors influences parasite dynamics, some studies have shown that parasites can be completely absent in specific sea bird populations or it can be as high as 80% in Afro-tropical water birds, which seems to be an effect of inhabiting either marine or freshwater environments (in combination with both phylogenetic and life-history traits) and possibly as a result of the abundance and presence of the Dipteran vectors [44,45,46]. However, our results suggest that vector abundance and seasonal variation are not always important drivers of infection risk . Other factors and processes beyond vector abundance could be more influential, for example host-selection by the vector  or host-compatibility and specialization of the parasites . Vectors of avian haemosporidians may play a major role in areas where they are extremely abundant , or when their influence is more pronounced at times when birds are prone to infections because of suppression of the immune system, e.g., after long distance migration  or during the breeding season . Moreover, parasite induced phenotypic change in the host is a widespread phenomenon that also influences transmission dynamics. Some studies conducted on avian malaria have shown enhanced vector attraction toward infected hosts (indirect manipulation hypothesis; (e.g. ), implying a relation between infection status and vector behavior. Finally, our data is grouping the abundance of many different species within Diptera families, which can dilute subtle effects on transmission dynamics by particular vector species; thus, future studies should consider to tackle this issue whenever Diptera sample sizes of individual species allows (e.g. ).
We did not find any effect of the a priori determined forest categories on parasite infection status, which indicates limited effect of coarse forest management categories on infection dynamics in our system. However, in similar forest types and ecosystems, it was shown that parasite infection parameters of F. coelebs were significantly higher in beech forests compared to spruce and mixed forests, and for S. atricapilla there was significantly higher prevalence in spruce dominated forest stands than in mixed and beech forests [16, 19]. Such differences were interpreted as a result of increased intra-specific competition for suitable habitat (i.e., beech forest for F. coelebs), or as living in suboptimal habitats (i.e., S. atricapilla in spruce forests); in both cases there was a negative impact of forest structure parameters on body condition and immune status, increasing parasite infection risk [16, 19]. In our study, we analyzed host-vector-parasite interactions within beech forests only; beech forest age classes did not per se represent the intricate structural differences as a result of management, as for example spruce forest does. However, other studies show that parasite prevalence varies across bird species, even within similar habitat types [16, 18, 19]. In this study, using quantitative forest structural variables (i.e., LiDAR) allowed a better understanding of ecosystem processes or species distributions , and has value for disentangling the complex relationships between the effects of land use (management) and forest structure on host-parasite interactions . We here confirm the relevance of fine-scale forest structural parameters (e.g., gaps, south facing slope, entropy) on ecosystem processes. Therefore, we recommend that future analyses of the effects of land use on host-vector-parasite systems include continuous structural parameters whenever possible, instead of relying on coarse habitat classes.
We found that parasite prevalence and parasitemia in seven avian hosts are affected by differences in fine-scale habitat structural parameters. In some cases, fine-scale habitat translates into negative indirect impacts on bird condition (i.e., physical asymmetry). Land use, such as forest fragmentation, deforestation, and urbanization have large effects on parasite dynamics at local geographical scales (e.g. [25, 55]). Other studies have shown that land use management influences all components of the host-vector-parasite system: avian hosts , vector abundance (e.g. ), and parasite prevalence . Our study confirms some parts of previous research but adding new insights and confirm that even subtle changes in forest structure, can have significant effects on host-vector-parasite interactions. Our study results have relevance to other vector-borne diseases and ecosystem processes, and it emphasizes the need to consider species-specific habitat effects on parasite dynamics for each host species.
Our study was part of the DFG Biodiversity Exploratories, a long-term and large-scale biodiversity research scheme in Germany . The study site is located in the northeast of Germany, within the Schorfheide-Chorin Biosphere Reserve (52°58′ N, 13°45′ E). The Biosphere Reserve is a young lowland glacial landscape, covering about 1,300 km2, and altitude ranges between 3 and 140 m.a.s.l., with mean annual precipitation of 500–600 mm . The study site contains 50 experimental forest plots, from which we sampled a total of 21 plots. All plots are 100 × 100 m with at least a 30 m buffer of the same forest structure per plot and within the buffer. We focused on forest stands dominated (i.e., > 70% canopy cover) by European beech (Fagus sylvatica) and selected three forest categories based on median stand age. The three forest categories were: young (age class; < 136 years), old (age class; > 136 years) and unmanaged (for at least 60 years). We selected seven replicate forest plots for each of the three forest categories, resulting in 21 plots (Additional file 1: Table S6).
We further categorized forest structure and heterogeneity at the plot level using LiDAR (Light Detection and Ranging) to derive continuous forest structural parameters . LiDAR is an active remote sensing system that can be used to map structure including vegetation height, density and other characteristics over large regions (e.g. ecosystem or regional scale). LiDAR variables were selected based on their relevance to the research topic of bird-vector-parasite interactions: forests entropy (local vertical variation in canopy structure or canopy surface heterogeneity), gap (proportion of canopy characterized as gap), and south facing aspect (proportion of south facing canopy) [c.f. 19, 20, 58].
We captured potential vectors on the same day when we captured avian hosts. We used one BG Sentinel Trap per forest plot, baited with CO2 (outflow of 200 g per day), Biogents Sweet-scent and an ultraviolet module to improve catch rate (Biogents AG, Regensburg, Germany), traps were placed in the center of the forest plot on the forest floor. All traps ran for 18 h, starting typically around 1 p.m. until the next morning. Insects were placed in a freezer at − 20 °C before immersion in 80% analytical grade ethanol. Insects were classified to the family level  and sorted to sex by morphology using a stereo microscope, species were determined by PCR and sequencing [20, 60]. Vector abundance was calculated as the sum of all female individuals per forest plot for the families Culicidae, Ceratopogonidae, and Simuliidae that respectively transmit avian haemosporidians of the genera Plasmodium, Haemoproteus and Leucocytozoon .
Host model species
Seven avian host species were selected for our study based on capture events and their habitat preferences: Eurasian blue tit (Cyanistes caeruleus), European robin (Erithacus rubecula), common chaffinch (Fringilla coelebs), great tit (Parus major), Eurasian blackcap (Sylvia atricapilla), common blackbird (Turdus merula) and song thrush (Turdus philomelos). We merged all captures of a given avian host species within each forest category.
All seven avian host species are common in forests, but with varying preferences for forest structure and forest types. C. caeruleus prefers all kinds of tree species, avoiding conifers and open stands . E. rubecula prefers medium dense forests, both broad leaved and mixed woodlands with cool and moist places while avoiding open spaces . F. coelebs is an arboreal species with a preference for mixed deciduous forests during the breeding season but can be found also in other mixed and conifer forests . P. major is a lowland species, avoiding conifer forest stands, and prefers mixed deciduous over pure deciduous forest stands, also prefers high structure of the undergrowth for cover and foraging . S. atricapilla forages and sings higher in the canopy, requires not too dense undergrowth for nesting and particularly avoids pure conifer forest stands . Both T. merula and T. philomelos tolerate wet and humid conditions, especially for invertebrate food items, and both can be found in diverse habitats, from dense woodlands to open landscape. T. merula is least frequent in the canopy and forages in all other layers and has no preference for conifer versus broadleaved forest stands. T. philomelos can be found where tree species provide enough shade and undergrowth, like beech forests .
Bird sampling and blood smear preparation
Birds were sampled three times in 2014 and two times in 2015 on each plot from April until June—with at least 14 days in-between sampling visits of the same plot . We captured birds from April until June using eight mist nets of 9 m × 2.5 m (mesh size 16 mm, nylon, Ecotone), placed at four locations within each experimental plot. Nets were opened 30 min after sunrise and left open for five consecutive hours aiming at the bird activity peak. Nets were checked with intervals of approximately 20 min. We used playback stations near the nets, playing territorial songs of the seven avian host species for improved capture success. We took measurements on both sides of each individual using a digital caliper (tarsus and third primary feather from the outside, i.e., primary 3) . Asymmetry was calculated as a ratio between right and left measurement, and as an absolute deviation from 1 (a ratio of 1 implies symmetry) to avoid using a measure for directed asymmetry. Birds were ringed and any recaptures within the same sampling period were excluded from the analysis.
About 30 ɥl of blood were taken from the brachial vein using a micro-capillary tube. Blood smears were prepared shortly after drawing the blood and another part was saved for molecular diagnostics . Blood smears were air dried, fixed in absolute methanol and stained with Giemsa for microscopic analyses (a solution consisting of 1.0 g Na2HPO4 + 0.7 g KH2PO4 + 1 L of double distilled H2O).
Hematological parameters (parasite prevalence, parasitemia, leucocyte count and H/L-ratio)
Blood smears were scanned with a Leica microscope (DM 5500B) and analyzed in oil immersion at 1000× magnification. Images of blood slides were screened manually to count the number of blood parasites (per 10,000 erythrocytes) as a measure of parasite load (i.e., parasitemia). We counted the leucocytes and determined heterophil to lymphocyte ratio (H/L-ratio) to assess birds' immune response . The H/L-ratio is a measure of stress in birds and is positively related to the magnitude of the stressor . We used a nested PCR protocol which allows simultaneous screening of Haemosporidia of the genera Plasmodium, Haemoproteus and Leucocytozoon , but cannot make a distinction between the first two genera.
For the analysis, we used a combined prevalence measure based on both PCR and microscopy. This means that when either technique showed a positive infection, we regarded this avian host as infected. If we detected an infection only through PCR, we assumed a very low parasitemia (a parasite count of 1) not detected by using the microscope . However, in such cases we repeated each PCR at least two times to avoid false positives.
Data analysis and statistics
All data was log-transformed (log10 x + 1) except for parasite prevalence (binomial), and were analyzed using R, version 4.0.1 . Mood's median tests were carried out if the data did not follow normality (i.e., with parasitemia) using the R-package "RVAideMemoire"  using the approach as outlined in Rózsa et al. [73, 74]. Generalized Linear Mixed Models (GLMM) were applied using the "lme4" package  and test results were obtained using the "lmerTest" and "piecewiseSEM" packages [76, 77], the factor bird species was used as a random effect in order to capture species-specific responses [18, 19]. Results from tests were considered statistically significant if the P-value was less than 0.05.
Initially, we tested the effect of forest categories on parasite prevalence or parasitemia per avian host species. The effect of forest category (factor, 3-levels) on parasite prevalence was tested by fitting one GLM (binomial data structure) and Mood's median tests for the effect on parasitemia (data not normally distributed). Subsequently, we used a model selection approach in combination with Generalized Linear Mixed Models (GLMM) to analyze the data in a broader context to understand mean associations between parasite infections and forest categories, avian host species, vector abundance, body asymmetry and immune response. We tested the effect of forest category, female vector abundance, H/L-ratio, leucocyte count and both tarsus and primary-3 asymmetry on both parasite prevalence and parasitemia. Except for 'forest category', all predictors were continuous. We stepwise removed one predictor variable at a time (based on the highest P-value) and selected a best-fitted model based on the largest difference of the Akaike Information Criterion corrected for small sample sizes (ΔAICc) to the full model, considering a better model-fit with decreasing AICc from the initial model [78, 79]. The results of the model selection formed the basis of the subsequent approach were model estimates were acquired using R-package "MuMIn" 
In order to have a more thorough understanding of the ecological dynamics in our host-vector-parasite systems, we applied Structural Equation Modeling (SEM), which allows analyzing alternative pathways of associations among variables at the same time, helping to determine the strength of association (covariance) between variables [81, 82]. We tested the original models, as well as the amended models that included a variable for vector abundance. In the model selection approach, we grouped parasite genera and vector families together, but since the parasite genera of Plasmodium, Haemoproteus and Leucocytozoon are each transmitted by dipteran vectors of the different dipteran families (i.e., Plasmodium-Culicidae, Haemoproteus-Ceratopogonidae and Leucocytozoon-Simuliidae), we made further additional parasite-vector specific analysis. However, Plasmodium and Haemoproteus infections cannot be separated with the PCR protocol of Hellgren et al.  and thus their respective vectors were grouped . We included parasite genus and vector family explicitly in the models (Plasmodium-Culicidae and Haemoproteus-Ceratopogonidae vector abundance). We applied all models to all bird species merged and separately for each bird species (i.e., overall and species-specific model; Additional file 1: Figs. S2, S3; Tables S2–S5) and separated the models by each tarsus and primary 3 asymmetries. Each model was analyzed separately with either the full dataset, or with sub-selections in case of the species-specific models. Models with a P-value (based on the chi-square value, χ2) < 0.05 were rejected, and we subsequently retained models with both Goodness-of-Fit Index (GFI) and Comparative Fit Index (CFI) of 0.9 or higher . From all models that met the criteria, the model with the lowest AIC (Akaike Information Criterion) value was selected as the best-fitted model . In the figures, path coefficients are indicated as non-standardized parameter estimates and the value in parenthesis indicate standard error values [83, 84]. SEMs were created using the R-package "lavaan", version 0.5–23.1097 .
Availability of data and materials
All data generated or analyzed during this study are included in this published article and its Additional files.
generalized linear model
generalized linear mixed model
heterophil to lymphocyte ratio
light detection and ranging
polymerase chain reaction
structural equation model
Lafferty KD, Allesina S, Arim M, Briggs CJ, De Leo G, Dobson AP, Dunne JA, Johnson PT, Kuris AM, Marcogliese DJ, et al. Parasites in food webs: the ultimate missing links. Ecol Lett. 2008;11(6):533–46.
Hatcher MJ, Dick JTA, Dunn AM. Diverse effects of parasites in ecosystems: linking interdependent processes. Front Ecol Environ. 2012;10(4):186–94.
Garcia-Longoria L, Møller AP, Balbontín J, de Lope F, Marzal A. Do malaria parasites manipulate the escape behaviour of their avian hosts? An experimental study. Parasitol Res. 2015;114(12):4493–501.
Sinclair ARE, Fryxell JM, Caughley G. Wildlife ecology, conservation, and management. New York: Blackwell Publishing; 2006.
Hammers M, Komdeur J, Kingma SA, Hutchings K, Fairfield EA, Gilroy DL, Richardson DS. Age-specific haemosporidian infection dynamics and survival in Seychelles warblers. Sci Rep. 2016;6:29720.
Santiago-Alarcon D, Carbó-Ramírez P, Macgregor-Fors I, Chávez-Zichinelli CA, Yeh PJ. The prevalence of avian haemosporidian parasites in an invasive bird is lower in urban than non-urban environments. Ibis. 2018. https://doi.org/10.1111/ibi.12699.
Hudson PJ, Dobson AP, Lafferty KD. Is a healthy ecosystem one that is rich in parasites? Trends Ecol Evol. 2006;21(7):381–5.
Birds occur in all major habitat types, with forest being particularly important. https://datazone.birdlife.org/sowb/casestudy/birds-occur-in-all-major-habitat-types-with-forest-being-particularly-important
Sala OE, Chapin FS III, Armesto JJ, Berlow E, Bloomfield J, Rodolfo D, Huber-Sanwald E, Huenneke LF, Jackson RB, Kinzig A, et al. Global biodiversity scenarios for the year 2100. Science. 2000;287:1770–4.
Schall P, Gossner MM, Heinrichs S, Fischer M, Boch S, Prati D, Jung K, Baumgartner V, Blaser S, Böhm S, et al. The impact of even-aged and uneven-aged forest management on regional biodiversity of multiple taxa in European beech forests. J Appl Ecol. 2018;55(1):267–78.
Wells K, Böhm SM, Boch S, Fischer M, Kalko EKV. Local and landscape-scale forest attributes differ in their impact on bird assemblages across years in forest production landscapes. Basic Appl Ecol. 2011;12(2):97–106.
Renner SC, Baur S, Possler A, Winkler J, Kalko EKV, Bates PJJ, Mello MAR. Food preferences of winter bird communities in different forest types. PLoS ONE. 2012;7(12):1–10.
Gregory RD, Gaston KJ. Explanations of commonness and rarity in British breeding birds: separating resource use and resource availability. Oikos. 1996;2000(88):515–26.
Carrascal LM, Senar JC, Mozetich I, Uribe F, Domenech J. Interactions among environmental stress, body condition, nutritional status, and dominance in Great Tits. Auk. 1998;115(3):727–38.
Carbonell R, Perez-Tris J, Tellería JL. Effects of habitat heterogeneity and local adaptation on the body condition of a forest passerine at the edge of its distribitional range. Biol J Lin Soc. 2003;78:479–88.
Lüdtke B, Moser I, Santiago-Alarcon D, Fischer M, Kalko EKV, Schaefer HM, Suarez-Rubio M, Tschapka M, Renner SC. Associations of forest type, parasitism and body condition of two European passerines, Fringilla coelebs and Sylvia atricapilla. PLoS ONE. 2013;8(12):e81395–e8139581395.
Marzal A, Asghar M, Rodríguez L, Reviriego M, Hermosell IG, Balbontín J, Garcia-Longoria L, de Lope F, Bensch S. Co-infections by malaria parasites decrease feather growth but not feather quality in house martin. J Avian Biol. 2013;44:1.
Santiago-Alarcon D, MacGregor-Fors I, Falfan I, Ludtke B, Segelbacher G, Schaefer HM, Renner S. Parasites in space and time: a case study of haemosporidian spatiotemporal prevalence in urban birds. Int J Parasitol. 2019;49(3–4):235–46.
Renner SC, Lüdtke B, Kaiser S, Kienle J, Schaefer HM, Segelbacher G, Tschapka M, Santiago-Alarcon D. Forests of opportunities and mischief: disentangling the interactions between forests, parasites and immune responses. Int J Parasitol. 2016;46(9):571–9.
van Hoesel W, Marzal A, Magallanes S, Santiago-Alarcon D, Ibáñez-Bernal S, Renner SC. Management of ecosystems alters vector dynamics and haemosporidian infections. Sci Rep. 2019;9(1):8779.
Valkiūnas G. Avian malaria parasites and other haemosporidia. New York: CRC Press; 2005.
Palinauskas V, Valkiūnas G, Bolshakov CV, Bensch S. Plasmodium relictum (lineage P-SGS1): effects on experimentally infected passerine birds. Exp Parasitol. 2008;120(4):372–80.
Marzal A, de Lope F, Navarro C, Moller AP. Malarial parasites decrease reproductive success: an experimental study in a passerine bird. Oecologia. 2005;142(4):541–5.
Brown CR, Brown MB. Ectoparasites cause increased bilateral asymmetry of naturally selected traits in a colonial bird. J Evol Biol. 2002;15(6):1067–75.
Hernández-Lara C, González-García F, Santiago-Alarcon D. Spatial and seasonal variation of avian malaria infections in five different land use types within a Neotropical montane forest matrix. Landscape Urban Planning. 2017;157:151–60.
Abella-Medrano CA, Ibáñez-Bernal S, MacGregor-Fors I, Santiago-Alarcon D. Spatiotemporal variation of mosquito diversity (Diptera: Culicidae) at places with different land-use types within a neotropical montane cloud forest matrix. Parasi Vectors. 2015;8(1):487–487.
Ferraguti M, Martínez-de la Puente J, Bensch S, Roiz D, Ruiz S, Viana DS, Soriguer RC, Figuerola J. Ecological determinants of avian malaria infections: an integrative analysis at landscape, mosquito and vertebrate community levels. J Anim Ecol. 2018;87(3):727–40.
Abella-Medrano CA, Ibanez-Bernal S, Carbo-Ramirez P, Santiago-Alarcon D. Blood-meal preferences and avian malaria detection in mosquitoes (Diptera: Culicidae) captured at different land use types within a neotropical montane cloud forest matrix. Parasitol Int. 2018;67(3):313–20.
Lachish S, Knowles SCL, Alves R, Sepil I, Davies A, Lee S, Wood MJ, Sheldon BC. Spatial determinants of infection risk in a multi-species avian malaria system. Ecography. 2012;36(5):587–98.
Martínez-Abraín A, Esparza B, Oro D. Lack of blood parasites in bird species: does absence of blood parasite vectors explain it all? Ardeola. 2004;51(1):225–32.
Lens L, Van Dongen S, Norris K, Githiru M, Matthysen E. Avian persistence in fragmented rainforest. Science. 2002;298(5596):1236–8.
Ellis VA, Medeiros MCI, Collins MD, Sari EHR, Coffey ED, Dickerson RC, Lugarini C, Stratford JA, Henry DR, Merrill L, et al. Prevalence of avian haemosporidian parasites is positively related to the abundance of host species at multiple sites within a region. Parasitol Res. 2017;3:73–80.
Medeiros MCI, Ricklefs RE, Brawn JD, Hamer GL. Plasmodium prevalence across avian host species is positively associated with exposure to mosquito vectors. Parasitology. 2015;142:1612.
Santiago-Alarcon D, Ricklefs RE, Parker PG. Parasitism in the endemic Galápagos Dove (Zenaida galapagoensis) and its relation to host genetic diversity and immune response. Stud Avian Biol. 2012;42:31–41.
Marzal A, Reviriego M, Hermosell IG, Balbontin J, Bensch S, Relinque C, Rodriguez L, Garcia-Longoria L, de Lope F. Malaria infection and feather growth rate predict reproductive success in house martins. Oecologia. 2013;171(4):853–61.
Coon CA, Garcia-Longoria L, Martin LB, Magallanes S, de Lope F, Marzal A. Malaria infection negatively affects feather growth rate in the house sparrow Passer domesticus. J Avian Biol. 2016;47(6):779–87.
Anciaes M, Marini M. The effects of fragmentation on fluctuating asymmetry in passerine birds of Brazilian tropical forests. J Appl Ecol. 2000;37(6):1013–28.
Møller AP. Developmental stability and fitness: a review. Am Nat. 1997;149(5):916–32.
Møller AP. Parasitism and developmental instability of hosts: a review. Oikos. 1996;77:189–96.
Valkiūnas G, Iezhova TA, Križanauskienė A, Palinauskas V, Sehgal RN, Bensch S. A comparative analysis of microscopy and PCR-based detection methods for blood parasites. J Parasitol. 2008;94(6):1395–401.
Santiago-Alarcon D, MacGregor-Fors I, Kühnert K, Segelbacher G, Schaefer HM. Avian haemosporidian parasites in an urban forest and their relationship to bird size and abundance. Urban Ecosyst. 2016;19(1):331–46.
Scheuerlein A, Ricklefs RE. Prevalence of blood parasites in European passeriform birds. Proc Biol Sci R Soc. 2004;271(1546):1363–70.
Schmitz OJ. Resolving ecosystem complexity (MPB-47). Princeton: Princeton University Press; 2010.
Lutz HL, Hochachka WM, Engel JI, Bell JA, Tkach VV, Bates JM, Hackett SJ, Weckstein JD. Parasite prevalence corresponds to host life history in a diverse assemblage of afrotropical birds and haemosporidian parasites. PLoS ONE. 2015;10(4):1–24.
Quillfeldt P, Arriero E, Martínez J, Masello JF, Merino S. Prevalence of blood parasites in seabirds-a review. Front Zool. 2011;8(1):26.
Mendes L, Piersma T, Lecoq M, Spaans B, Ricklefs RE. Disease-limited distributions? Contrasts in the prevalence of avian malaria in shorebirds species using marine and freshwater habitats. Oikos. 2004;2005(109):396–404.
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(19):8040–57.
Pettersson E, Bensch S, Ander M, Chirico J, Sigvald R, Ignell R. Molecular identification of bloodmeals and species composition in Culicoides biting midges. Med Vet Entomol. 2013;27(1):104–12.
Ellis VA, Collins MD, Medeiros MC, Sari EH, Coffey ED, Dickerson RC, Lugarini C, Stratford JA, Henry DR, Merrill L, et al. Local host specialization, host-switching, and dispersal shape the regional distributions of avian haemosporidian parasites. Proc Natl Acad Sci USA. 2015;112(36):11294–9.
Ishtiaq F, Renner S. Bird migration and vector-borne parasite transmission. In: Santiago-Alarcon D, Marzal A, editors. Avian malaria and related parasites in the tropics. Berlin: Springer International Publishing; 2020.
Knowles SC, Nakagawa S, Sheldon BC. Elevated reproductive effort increases blood parasitaemia and decreases immune function in birds: a meta-regression approach. Funct Ecol. 2009;23(2):405–15.
Yan J, Martínez-de la Puente J, Gangoso L, Gutiérrez-López R, Soriguer R, Figuerola J. Avian malaria infection intensity influences mosquito feeding patterns. Int J Parasitol. 2018;48(3–4):257–64.
Martínez-De La Puente J, Ferraguti M, Ruiz S, Roiz D, Soriguer RC, Figuerola J. Culex pipiens forms and urbanization: effects on blood feeding sources and transmission of avian Plasmodium. Malar J. 2016;15(1):589.
Bradbury RB, Hill RA, Mason DC, Hinsley SA, Wilson JD, Balzter H, Anderson GQA, Whittingham MJ, Davenport IJ, Bellamy PE. Modelling relationships between birds and vegetation structure using airborne LiDAR data: a review with case studies from agricultural and woodland environments. Ibis. 2005;147(3):443–52.
Sehgal RNM. Deforestation and avian infectious diseases. J Exp Biol. 2010;213(6):955–60.
Fischer M, Bossdorf O, Gockel S, Hänsel F, Hemp A, Hessenmöller D, Korte G, Nieschulze J, Pfeiffer S, Prati D, et al. Implementing large-scale and long-term functional biodiversity research: the biodiversity exploratories. Basic Appl Ecol. 2010;11(6):473–85.
Schall P, Ammer C. EP-summary of Lidar derived forest structure parameter v1.1.4. Repository: https://www.bexis.uni-jena.de/.DatasetId=17066. 2016.
Renner SC, Suarez-Rubio M, Kaiser S, Nieschulze J, Kalko EKV, Tschapka M, Jung K. Divergent response to forest structure of two mobile vertebrate groups. For Ecol Manage. 2018;415–416:129–38.
Oosterbroek P. The European families of the Diptera - Identification, diagnosis, biology. Utrecht: KNNV Publishing; 2006.
Martínez-de la Puente J, Martínez J, Ferraguti M, Morales-del Nuez A, Castro N, Figuerola J. Genetic characterization and molecular identification of the bloodmeal sources of the potential bluetongue vector Culicoides obsoletus in the Canary Islands. Spain. Parasit Vectors. 2012;5:147–147.
Cramp S. Handbook of the birds of Europe, the Middle East and North Africa. The birds of the western Palearctic, vol. VII. Flycatchers to Shrikes. Oxford: Oxford University Press; 1993.
Cramp S. Handbook of the birds of Europe, the Middle East and North Africa. The birds of the Western Palearctic, vol. V. Tyrant Flycatchers to Thrushes; 1988.
Cramp S. Handbook of the birds of Europe, the Middle East and North Africa. The birds of the western Palearctic, vol. VIII. Crows to Finches. Oxford: Oxford University Press; 1994.
Cramp S. Handbook of the birds of Europe, the Middle East and North Africa. The birds of the Western Palearctic, vol. VI. Warblers. Oxford: Oxford University Press; 1992.
Renner SC, van Hoesel W. Ecological and Functional Traits in 99 Bird Species over a Large-Scale Gradient in Germany. Data. 2017;1:13.
Svensson L. Identification Guide to European Passerines. Thetford: British Trust for Ornithology; 1992.
Campbell TW, Ellis CK. Avian and exotic animal hematology and cytology. New York: Wiley; 2013.
Davis A, Maney D, Maerz J. The use of leukocyte profiles to measure stress in vertebrates: a review for ecologists. Funct Ecol. 2008;22(5):760–72.
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(4):797–802.
Krams I, Suraka V, Cīrule D, Hukkanen M, Tummeleht L, Mierauskas P, Rytkönen S, Rantala MJ, Vrublevska J, Orell M. A comparison of microscopy and PCR diagnostics for low intensity infections of haemosporidian parasites in the Siberian tit Poecile cinctus. Ann Zool Fennici. 2012;1:331–40.
R Core Team. R: A language and environment for statistical computing. Vienna: R Foundation for Statistical Computing; 2020.
Hervé M. RVAideMemoire: testing and plotting procedures for biostatistics. R package version 0.9–68. 2017.
Rózsa L, Reiczigel J, Majoros G. Quantifying parasites in samples of hosts. J Parasitol. 2000;86(2):228–32.
Reiczigel J, Marozzi M, Fábián I, Rózsa L. Biostatistics for parasitologists–a primer to quantitative parasitology. Trends Parasitol. 2019;35(4):277–81.
Bates D, Mächler M, Bolker B, Walker S. Fitting linear mixed-effects models using lme4. arXiv preprint arXiv:14065823. 2014.
Kuznetsova A, Bruun Brockhoff P, Haubo Bojesen Christensen R. lmerTest: Tests in Linear Mixed Effects Models 2.0–33. 2016.
Lefcheck JS. piecewiseSEM: Piecewise structural equation modelling in r for ecology, evolution, and systematics. Methods Ecol Evol. 2016;7(5):573–9.
Zuur AF, Ieno EN, Walker NJ, Saveliev AA, Smith GM. Mixed effects models and extensions in ecology with R; 2009.
Burnham KP, Anderson DR. Model Selection and Multimodel Inference. New York: Springer; 2004.
Bartoń K. Package ‘MuMIn’. 2016.
Grace JB. Structural Equation Modeling for Observational Studies. J Wildl Manag. 2008;72(1):14–22.
Santiago-Alarcon D, Palinauskas V, Schaefer HM. Diptera vectors of avian Haemosporidian parasites: untangling parasite life cycles and their taxonomy. Biol Rev Camb Philos Soc. 2012;87(4):928–64.
Fan Y, Chen J, Shirkey G, John R, Wu SR, Park H, Shao C. Applications of structural equation modeling (SEM) in ecological studies: an updated review. Ecol Process. 2016;5(1):19–19.
Rosseel Y. lavaan: An R Package for Structural Equation Modeling. J Stat Softw. 2012;48(2):1–36.
We thank the local management teams of the Biodiversity Exploratories for their incredible support in the field. We thank F. Fischer, W.W. Weisser, K.E. Linsenmair, and F. Buscot for their role in setting up the Biodiversity Exploratories. We thank P. Schall and C. Ammer for the available dataset with regards to LiDAR derived forest structure parameters. We thank J. Kienle, C. Leutgeb, B. Lüdtke, I. Moser, A. Gil-Pastor, M. Salzmann, B. Springer and M. Wojta for additional support.
The work has been funded by the DFG (http://www.dfg.de) Priority Program 1374 ‘Biodiversity-Exploratories’ (Re1733/6-1; SCHA1008/8-1). D.S.-A. was funded by the Alexander von Humboldt Foundation (post-doctoral grant) and by Consejo Nacional de Ciencia y Tecnología (CONACYT, project number CB-2011-01-168524). AM was funded by the Spanish Ministry of Economy and Competitiveness (CGL2015-64650P) and Junta de Extremadura (IB16121). The funding organizations or permit organizations had no role in study design, data collection and analysis, decision to publish, or preparation of the manuscript.
Ethics approval and consent to participate
All bird handling procedures, particularly capture and drawing blood, were performed in compliance with laws and regulations of the German federal and Brandenburg state regulations (Bundesnaturschutzgesetz, BNatSchG, Federal Act for the Protection of Nature). Permits were granted by the Landesamt für Umwelt, Gesundheit und Verbraucherschutz, Potsdam (animal protection permit number: Para-Aves, 2347-3-2014 and capture permit for birds RO/4610/SOB-0907D).
Consent for publication
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.
About this article
Cite this article
van Hoesel, W., Santiago-Alarcon, D., Marzal, A. et al. Effects of forest structure on the interaction between avian hosts, dipteran vectors and haemosporidian parasites. BMC Ecol 20, 47 (2020). https://doi.org/10.1186/s12898-020-00315-5