Landscape structure and agricultural intensification are weak predictors of host range and parasitism rate of cereal aphids

Proportions of specialist and generalist primary parasitoids have been described by the resource breadth and the trade-off hypothesis. These alternative hypotheses predict either decreased or increased, respectively, parasitism rate of shared aphid species by specialist parasitoids. We tested both hypotheses and the confounding effects of landscape structure and agricultural intensification (AI) using extensive samplings of aphids and their parasitoids in Polish agricultural landscapes. Abundances, species composition of aphids, primary parasitoids, and parasitism rate of aphids by specialists and generalist parasitoids were analysed. Contrary to our expectations we found equally decreased parasitism rates by both types of primary parasitoids at higher aphid densities and thus proportion of specialists to generalists did not change with increasing host density. In line with the resource breadth hypothesis, specialist parasitoids had always lower abundances and parasitism rates than generalist parasitoids. Landscape diversity and agricultural intensification did not influence the host-parasitoid population dynamics. We speculate that these contrasting results could be caused by the additional density effects of secondary parasitoids. We conclude that simplistic two-trophic-level population models are not able to fully describe the complex dynamics of trophic networks. We also argue that agricultural intensification has lower effects on abundance and effectiveness of parasitoids than predicted by respective predator–prey models and empirical studies performed in controlled and artificial conditions.


Introduction
Aphids are important pests in agriculture and much effort has been made to reduce their negative impact on crop plants (Hassan 1994;Landis et al. 2000). One of the most diverse group of aphid enemies are primary parasitoids that have a high potential in limiting aphid numbers, especially under controlled conditions e.g. greenhouses (Yang et al. 2014). Many ecological factors have been shown to influence primary parasitoid effectiveness (e.g. predation, secondary parasitism, alternative prey). The host range of parasitoids (the degree of parasitoid host specialisation) may also affect parasitization rate (Monmany and Aide 2009;Rossinelli and Bacher 2015). For instance, at low densities of the target host (e.g. cereal aphids) the ability of generalist parasitoids to switch to other hosts might strongly reduce the parasitism rate of target aphid species allowing aphid populations to recover (Chow and Mackauer 1991).
Within primary parasitoids various degrees of host specialization can be found (Stilmant et al. 2008). Respective proportions of specialist and generalist parasitoids within a focal community have been described by two alternative hypotheses. Both hypotheses assume differences between generalists and specialists in efficiency of using shared hosts. According to the trade-off hypothesis the high performance of parasitoids on specific hosts comes to the cost of narrow host range (Poulin 1998). Therefore abundances and parasitism rates by specialist parasitoids will be higher in comparison to their generalist counterparts. The alternative resource breadth hypothesis assumes that generalists and specialists are equally efficient in exploiting shared hosts resulting in a predominance of generalist parasitoids due to the larger prey base that generalists can utilize (Gaston et al. 1997;Krasnov et al. 2004). However, recent tests of these contrasting hypotheses returned mixed results without giving a strong and cogent evidence in favour of one of these two views (Jaenike 1990;Scheirs et al. 2005;Agosta and Klemens 2009;Forister et al. 2012). In this respect, it does matter what type of parasitoid is the most effective in reducing pest numbers and which environmental factors favour this process.
Different proportions of specialist and generalist parasitoids stem from differences in the mode of parasitoid searching (Campan and Benrey 2004) and their ability to compete for the host (Brodeur and Rosenheim 2000;Sampaio et al. 2006). Apart from these factors, the landscape structure and agricultural intensification (AI) also trigger differences in abundance of these two types of parasitoids (Tscharntke et al. 2007).
High agricultural intensification is associated with environmental disturbances. It makes agro-ecosystems less stable in comparison to natural and semi-natural habitats (Kennedy and Storer 2000). According to quantitative modelling of Richmond et al. (2005), generalists prevail and contribute most to ecosystem functioning in unstable environments, while specialists excel under constant or slowly changing environmental conditions. Therefore the ratio of specialists to generalists may depend on the degree of ecosystem stability and may be different in various types of habitats. Hence, various regimes of agricultural intensification and landscape structure are potentially important and should be considered in testing of both hypotheses.
Field studies investigating the proportion of specialists to generalists are scarce (Straub et al. 2011;Rossinelli and Bacher 2015), and studies examining these two hypotheses in relation to landscape structure and AI basically do not exist.
In biocontrol programmes it is important to know which of these two types of parasitoids cause higher mortality on aphids and consequently performs better. The present work tries to answer the question, whether aphid abundance, landscape structure, and agricultural intensification affected the abundance of specialist and generalist primary parasitoids. Based on the resource breadth hypothesis (equal efficiency in exploiting cereal aphids) we predicted that generalist parasitoids will prevail and exert the stronger pressure on less intensively managed fields (less disturbed ecosystems) located in complex landscape, since such a landscape maintains larger populations of parasitoids with broader host range.

Materials and methods
In 2008 we studied four fields with high agricultural intensification (high AI fields) and four with low AI. These fields were selected from a regional pool according to number of management events (e.g. tillage), cereal yield, and the percentage of surrounding arable land. In 2009 we studied five different fields of each category, but subsequently ruled out one field from the analyses due to inappropriate management. In consequence eight and nine winter wheat [Triticum aestivum (L.)] fields in Central Poland were surveyed in 2008 and 2009, respectively (study design of the AGRIPOPES project, see Geiger et al. 2010a, b). Fields were above one hectare in area and located at least one kilometre apart. They were distributed within a square of about 30 · 30 km to minimize differences in the regional species pools among study farms.
As recommended by Thies et al. (2003), we quantified landscape structure as the proportion of individual habitat types (grassland, forest, arable land) within a circle of 500 m radius (area of about 78.5 ha included) going from the middle of each study field. The Shannon index of habitat types was calculated as a metric of landscape diversity. Abundances and species composition of aphids were estimated from five sampling points per field along a transect going from the edge to the centre of each field. At each point we sampled 20 randomly chosen winter wheat shoots resulting in a total of 100 tillers per field during each visit.
Abundances and species composition of aphid parasitoids were assessed from aphid mummies, collected at random in the whole field during 2-h surveys. Mummies were taken to the laboratory and kept individually in small vials. After emergence, adult parasitoids were identified to species level. We did not consider secondary parasitoids. Samples of both living aphids and those mummified by primary parasitoids were collected three and nine times in the season in 2008 and 2009, respectively, leading to a total of 105 samples.
We estimated the parasitism rate k by specialist and generalist parasitoids as the ratio of the numbers of a given type of mummies (aphids mummified by specialist or generalist parasitoids) to the all mummies (aphids parasitized by both types of primary parasitoids) collected during 2-h surveys during each visit plus living aphids that were recorded on 100 tillers of winter wheat. We used the difference Dk ¼ k specialists Àk generalists to quantify the proportions of both types. Additionally we checked the preferences of specialist and generalist primary parasitoids for particular aphid species, counting aphids parasitized by the focal parasitoid species. We tested for differences between the expected and the observed counts in different parasitoid categories (specialists vs generalists, and Aphidius vs Ephedrus vs Praon) using the v 2 -statistic.
We related abundances of generalist and specialist parasitoids and their aphid parasitism rate (dependent variables) to agricultural intensification (categorical predictor), and landscape diversity, percentage of arable land, aphids abundances, and the parasitism rate by the opposite category parasitoids (specialists and generalists) (metric predictors) using general linear modelling with Gaussian error structure and log-link functions. Due to the non-linear relationships squared parasitism rates entered the models as additional covariate. To assess the strength of influence of the predictors on the dependent variables we used partial g 2 values. Errors refer always to standard errors.
Neither landscape diversity nor percentage of arable land did detectably influence parasitism rates by spe-cialist and generalist parasitoids (Table 2) and their abundances (equivalent results as for parasitism rates, not shown). Total abundances of aphids and parasitoids were highly variable during the two study years (Table 1). Aphid abundances did not differ significantly (P(t 1,104 ) > 0.15) between the low (163 ± 24 ind.) and the high AI sites (111 ± 18 ind.). Similarly, the proportions of specialist and generalist parasitoids as quantified by Dk did not significantly depend on AI ( Fig. 1a; Table 2). AI alone explained only about 3% of the variability in Dk, and at most 3% of variability in parasitism rates (Table 2). Lower AI consistently increased abundances of generalist and specialist parasitoids (Fig. 1a). At higher AI, parasitism rates of specialists increased and those of generalists decreased (Fig. 1b), making total parasitism rate of aphids k independent of agricultural intensification (k(high AI) = 0.69; k(low AI) = 0.67). Rates of parasitism by generalist species were highest at intermediate parasitism by the specialist species (Fig. 2). Further, the proportions of both groups changed in dependence on host number (Fig. 3). Parasitism rates of aphids by specialists (Fig. 3a) and generalists (Fig. 3b) at both levels of agricultural intensification decreased with increasing host density. Proportions of both groups were highly variable at low and less variable at high aphid abundances (Fig. 3c). These effects were observed at both, high and low levels of AI. We did not find significant differences between both types of parasitoids in host association (3 · 2 contingency table P(v 2 ) > 0.1). Generally, the generalist species of parasitoids reached  Fig. 1a) and parasitism rates (Fig. 1b) than their specialist counterparts.

Discussion
Based on the general assumption that the type of land use differently influences host aphid density and their parasitism by specialist and generalist parasitoids (Monmany and Aide 2009), we predicted an association of generalist parasitoids with low AI and the more complex landscape. In line with the resource breadth hypothesis, we expected that less disturbed cereal fields and a diverse environment will promote generalist parasitoids and enhance the respective parasitism rate. Our results corroborates this hypothesis showing that generalists were always more numerous and had a higher aphid parasitism rate than specialists. However, we did not find any effect of environmental factors on parasitism rates. Similar results were obtained by Macfadyen et al. (2009), who also found a prevalence of the generalist parasitoid Ephedrus plagiator, irrespective of agricultural intensification. Gagic et al. (2012Gagic et al. ( , 2014 reported generalists to be associated with low AI early in the season while later being abundant regardless of the AI regime. We did not confirm the assumed relationship between the performance of generalist/specialist parasitoids and Table 2 General linear modelling linking the parasitism rates of aphids k by Aphidius uzbekistanicus, A. rhopalosiphi (specialists) and the generalist parasitoids A. ervi and Ephedrus spp. and Praon spp., as well as the difference in parasitism rate Dk to agricultural intensity AI, to the percentage of arable land, plant species richness, landscape diversity (metric variables), and to the presences of potential competitors. Given are partial g 2 values, and coefficients of determination R 2 of the whole model Parametric significance levels: * P <0.05 Fig. 1 The average abundances (a) and parasitism rates of aphids k (b) by specialist (dark grey bars) and generalist parasitoids (light grey bars), and the difference between both type of parasitoids Dk (white bars) on plots with high and low agricultural intensification (high AI, low AI). Error bars denote one standard error  Richmond et al. (2005) showed that generalists tolerate suboptimal conditions and prevail under highly variable environment.
Associating the environmental conditions with host availability, we expected that specialist parasitoids should benefit from situations when there is access to one specific host (cereal aphids) that is abundant, whereas generalists should have an advantage when such host is rare or unpredictable (Rossinelli and Bacher 2015). Instead, we found decreased parasitism by both types of primary parasitoids with increasing aphid numbers possibly reflecting a ''spreading of mortality risk'' strategy of primary parasitoids attacked by secondary parasitoids (Fig. 3a, b). If the death rate of primary parasitoid larvae increases with the clutch size of living and parasitized aphids due to secondary parasitism, it pays for primary parasitoids to escape from such a colony rather than to increase the number of ovipositions (Mackauer and Vo¨lkl 1993). This decrease was independent of host specialization resulting in equal proportions of parasitism by specialist and generalist parasitoids at higher densities of cereal aphids (Fig. 3c).
In this respect, Thies et al. (2005) pointed to trade-offs between higher aphid mortality due to parasitism and simultaneously increased colonization in complex landscapes resulting in decreased parasitism rates at high host density. Consequently, we argue that a higher parasitism rate of aphids should be expected at lower aphid density. In general, the AI and landscape effects we tested for, were weak and at the lower level of statistical detectability. However, we could not assess possible long-term effects on parasitoid community structure and abundance. In fact, we did not find strong evidence that agricultural intensification has a negative impact on biological control potential of parasitoids. Parasitism rates recorded on intensively and extensively managed fields were similar and not linked to our metrics of landscape complexity (Table 2). This is in opposite to results of some authors, who showed that the deterioration of biocontrol potential were primarily mediated by agricultural practices such as pesticide application (Jonsson et al. (2012) and simplification of landscape (Macfadyen et al. 2009;Thies et al. 2011;Veres et al. 2013;Rusch et al. 2016). In this respect, the important factor in the evaluation of landscape effects might be initial complexity of landscape structure. Jonsson et al. (2010) and Tscharntke et al. (2012) favoured bell shaped parasitism rate-complexity relationships with the consequence that we should expect detectable parasitoid biocontrol effects only at intermediate levels of landscape complexity. We note that agricultural intensification in Poland was comparably low during the last 20 years. Even the most intensively managed fields were less exploited than average fields in western Europe (Wretenberg 2006). Hence, the weak effects of AI and landscape structure on the observed pattern of parasitism rate might be due to the overall low degrees of AI and respective effects might have been masked by the lack of sufficiently homogenous conditions (Macfadyen et al. 2009;Jonsson et al. 2015). Probably for the same reasons Menalled et al. (2003) did not find any relationship between landscape structure and species richness.
These contradictions may also result from the complexity of ecological processes that involve multiple biotic-abiotic and biotic-biotic interactions, which are difficult to examine and control in the field. In this respect Schellhorn et al. (2015) stressed the importance of the spatial and temporal variability of aphids and their natural enemies for the final outcome of pest control. Host availability in previous seasons determine the rate of aphid parasitism in subsequent years (Holt and Lawton 1994). Our study also revealed significant year by year differences in community structure and dominance of aphids and their primary parasitoids (Tables 1,  2). Such a variability has already been reported by many authors without identification of the underlying factors (Askew and Shaw 1986;Jones and Weinzierl 1997;Thies et al. 2005;Gagic et al. 2014). Possibly the observed high species turnover is linked to environmental stress connected with AI (Warwick and Clarke 1993;Carpenter and Brock 2006) secondary parasitoid impact (Sullivan 1972), or climatic factors affecting the immigration, colonization, and development of aphids and their parasitoids (Pankanin-Franczyk and Ceryngier 1995;Pan-Fig. 3 Parasitism rate of aphids k by specialist (a) and generalist parasitoids (b) decreased with increasing aphid abundance (total counts). The respective difference Dk between both type of parasitoids (c) was independent of aphid abundance. Open dots low AI, black dots high AI. Regression lines in a and b drawn from all data points kanin-Franczyk and Sobota 1998; Langer and Hance 2004;Schellhorn et al. 2014;Hawro et al. 2015).
A major shortcoming of many studies relating land use intensification and community composition regards the temporal dynamics as these studies used short-term data only. This reduced number of replicates makes the identification of the factors and processes that trigger the long-term dynamics of aphid and parasitoid population structure challenging (Chaplin-Kramer et al. 2011). Future investigations on aphid-parasitoid dynamics in agricultural landscapes need to involve longer time scales (Menalled et al. 2003;Gagic et al. 2012Gagic et al. , 2014. Acknowledgements The study was part of the AGRIPOPES project (AGRIcultural POlicy-induced landscaPe changes: effects on biodiversity and Ecosystem Services; http://agripopes.net), funded through the European Science Foundation (EUROCORES programme), in collaboration with the Polish Ministry of Science and Higher Education. We would like to thank Ž eljko Tomanowic´for helping with identification of some Aphidius individuals. We would like also to express our special thanks to anonymous reviewers for valuable comments.
Open Access This article is distributed under the terms of the Creative Commons Attribution 4.0 International License ( http://creativecommons.org/licenses/by/4.0/), which permits unrestricted use, distribution, and reproduction in any medium, provided you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made.