Does urbanization explain differences in interactions between an insect herbivore and its natural enemies and mutualists?

Urbanization can alter the composition of arthropod communities. However, little is known about how urbanization affects ecological interactions. Using experimental colonies of the black bean aphid Aphis fabae Scopoli reared on Vicia faba L, we asked if patterns of predator-prey, host-parasitoid and ant-aphid mutualisms varied along an urbanization gradient across a large town in southern England. We recorded the presence of naturally occurring predators, parasitoid wasps and mutualistic ants together with aphid abundance. We examined how biotic (green areas and plant richness) and abiotic features (impervious surfaces and distance to town center) affected (1) aphid colony size, (2) the likelihood of finding predators, mutualistic ants and aphid mummies (indicating the presence of parasitoids), and (3) how the interplay among these factors affected patterns of parasitoid attack, predator abundance, mutualistic interactions and aphid abundance. Aphid abundance was best explained by the number of mutualistic ants attending the colonies. Aphid predators responded negatively to both the proportion of impervious surfaces and to the number of mutualistic ants farming the colonies, and positively to aphid population size, whereas parasitized aphids were found in colonies with higher numbers of aphids and ants. The number of mutualistic ants attending was positively associated with aphid colony size and negatively with the number of aphid predators. Our findings suggest that for insect-natural enemy interactions, urbanization may affect some groups, while not influencing others, and that local effects (mutualists, host plant presence) will also be key determinants of how urban ecological communities are formed.


Introduction
Urbanization is the defining feature of recent history; today over 50% (>90% in developed countries) of people live in urban environments (United Nations 2014). Urbanization is arguably the greatest anthropogenic transformation that ecological systems experience, and while most studies of urban ecology focus on changes to the diversity and abundance of species inhabiting towns and cities, attention has only started to focus on how assemblages of interacting species are formed in urban areas, and how this is affected by the intensity of urbanization (Bennett and Gratton 2012;Quispe and Fenoglio 2015;Pereira-Peixoto et al. 2016;Turrini et al. 2016). Fragmentation reduces populations of native plants (Benitez-Malvido 1998;Jules 1998;Williams et al. 2005), leads to decreased connectivity between vegetation patches and existing patches tend to be smaller (Medley et al. 1995;McKinney 2002) and therefore of reduced quality as habitat for many animal species (Bradley and Altizer 2007;Faeth et al. 2011;Turrini et al. 2016). There are also some dramatic physical changes from increased densities of roads, buildings and other sealed structures and microclimatic changes such as the urban heat island effect (Bradley and Altizer 2007;Faeth et al. 2011). Together, these changes affect the likelihood of encountering species at higher trophic levels (Faeth et al. 2005;Egerer et al. 2017). Understanding how such extreme anthropogenic habitat changes may affect patterns of ecological interactions is perhaps most tractable with arthropod model systems (McIntyre 2000;Bang and Faeth 2011), but experimental studies in urban ecosystems are scarce.
Urbanization has been shown to affect the structure of insect communities, resulting in dramatic changes in their abundance and richness (McIntyre 2000;Grimm et al. 2008;Raupp et al. 2010;Gardiner et al. 2014), most frequently leading to a loss of diversity (Kahn and Cornell 1989;Suarez et al. 1998;McKinney 2002;Shochat et al. 2004;Sadler et al. 2006;Clark et al. 2007;Magura et al. 2010;Uno et al. 2010;Bang and Faeth 2011;Bennett and Gratton 2012;Ramírez Restrepo and Halffter 2013). Few studies have considered how these changes influences the outcome of ecological interactions at multiple trophic levels (Shrewsbury and Raupp 2006;Bennett and Gratton 2012;Fenoglio et al. 2013;Pereira-Peixoto et al. 2016;Turrini et al. 2016). For example, abiotic environmental factors might interfere with biotic interactions, thereby modulating the strength of the trophic effects on food webs (Ritchie 2000;Preisser and Strong 2004;Turrini et al. 2016). Mooney et al. (2016) investigated if variation in light availability (shaded understory or open meadow) determines the abundance of the aphid Aphis helianthi feeding on the herb Ligusticum porteri. Aphid numbers were higher in open meadows than in shaded environments. This pattern was not due to the direct effects of light on aphid performance, plant quality or interactions with natural enemies, but instead was due to an indirect effect mediated by a mutualistic relationship with ants, which were more abundant in meadows. If, as expected, insects and other arthropods do respond to habitat structure, then we can predict that there will be differences not only in species assemblages, but also on trophic dynamics and species interactions as habitat configuration changes with urbanization.
If we consider that urban wildlife is subject to multiple changes in abiotic conditions simultaneously, it is not surprising that predicting the consequences of such changes for trophic processes and for direct and indirect species interactions is highly challenging (Turrini et al. 2016). However, a few trends have begun to appear. Urban areas are often characterized by reduced numbers of native vertebrate predators (McKinney 2002;Shochat 2004), an increased abundance of some urban-adapted species, which can potentially lead to increased competition and displacement (Hostetler and McIntyre 2001), altered behavior and phenology (Connor et al. 2002;Neil and Wu 2006), high densities of herbivorous arthropods (Dreistadt et al. 1990;Hanks and Denno 1993;Tooker and Hanks 2000), lower numbers of arthropod predators (Turrini et al. 2016) and lower numbers of parasitoids (Denys and Schmidt 1998;Bennett and Gratton 2012;Burks and Philpott 2017). All these changes can potentially lead to altered trophic structure, and we must recognize that trophic dynamics cannot be understood based only on our knowledge of species composition (Shochat et al. 2006). This way, evaluating empirically how trophic dynamics behave in urban environments may help us to make some broad and useful predictions regarding the effects that urbanization could have on multi-trophic interactions.
With direct trophic interactions such as predation, one species has a negative effect on the other species, but in indirect interactions one species can also positively affect another species through intermediate levels in a trophic cascade (Halaj and Wise 2001;Müller et al. 2005;Turrini et al. 2016). For example, the presence of some species of honeydewcollecting ants results in increased aphid numbers and also increased numbers of aphid parasitoids when protecting aphids from predators and incidentally also protecting parasitized aphids against predators and hyperparasitoids (Völkl 1992;Kaneko 2002). Nevertheless, the most recognized indirect trophic interactions are top-down trophic cascades in which predators influence plants by feeding on herbivores, thus reducing the consequences of herbivory (Schmitz et al. 2000;Shurin et al. 2002;Turrini et al. 2016).
Traditionally, research on trophic interactions and food webs mainly focus on direct interactions such as predation or parasitism, therefore the importance of non-trophic, indirect, and facilitative interactions has been rarely taken into consideration (Ohgushi 2008). Facilitative or positive interactions, like mutualisms, are rarely considered as potential factors affecting urban populations and communities (but see e.g. Thompson and McLachlan 2007;Gibb and Johansson 2010;Toby Kiers et al. 2010), and it is claimed that this type of positive interaction plays an important part in the structuring of some biological communities by providing refuge from predation or competition (Stachowicz 2001). Conversely, it is important to consider that mutualisms have formed over evolutionary time scales, and we do not know if mutualisms have evolved to be resilient enough to endure anthropogenic disturbances (Sachs and Simms 2006; Toby Kiers et al. 2010).
Host-parasitoid interactions are also likely to be considerably altered in urban ecosystems. Here, plant resources for herbivorous insects and their parasitoids are spatially subdivided and embedded in a matrix of built environment (Bennett and Gratton 2012;Fenoglio et al. 2013). These conditions are particularly prone to altering insect colonization and persistence, which may lead to altered trophic interactions (Fenoglio et al. 2013). Parasitoid insects are important biological control agents of herbivorous insect populations and have been found to be negatively affected by urbanization at both local and landscape spatial scales (Bennett and Gratton 2012;Fenoglio et al. 2013). Parasitoids are specialist organisms closely associated with their hosts (Kruess and Tscharntke 1994). Consequently, they might present higher sensitivity to environmental fluctuation and anthropogenic disturbance in comparison to less specialized species (Gibb and Hochuli 2002). Since some herbivore pest populations are limited by top-down control by parasitoids (Hawkins and Gross 1992), a decrease in parasitism or predation can favour pest outbreaks in these areas (Schmitz et al. 2000;Roslin et al. 2014).
Even less frequently considered is how these different ecological interactions (host-parasitoid, predator-prey, mutualisms) act together to affect the insect assemblages found in urban environments. Systems including different types of interactions and trophic groups have only recently started to be empirically examined (Halaj and Wise 2001;Lurgi et al. 2016). In this work we explore these interconnected biological interactions in an urban environment. We used a study system which consisted of experimental colonies of the herbivorous aphid Aphis fabae Scopoli reared on an herbaceous plant species (the dwarf broad bean Vicia faba L.) and their naturally occurring predators, parasitoid wasps and mutualistic ants along an urbanization gradient in a large town in southern England.
Study sites varied in the amount of impervious surfaces, green areas, plant species richness and position on the urban gradient. Distance from the town center is a variable frequently used as a proxy for urban gradients, as cities and towns frequently show gradients of urbanization from their centers to their edges, and that the biotic and abiotic factors that can potentially affect biological systems tend to follow and change as function of this gradient, resulting from variation in human population density and intensity of activity (Deichsel 2006;Clark et al. 2007;Bang and Faeth 2011). The extent of impervious cover (paved surfaces, structures such as buildings and roads) causes a variety of detrimental effects on arthropods (Morse et al. 2003;Sadler et al. 2006;Magura et al. 2008;Bennett and Gratton 2014), and it is a stronger predictor of urbanization gradients than broad classifications such as urban, suburban and rural areas (Ellis and Ramankutty 2008;Ramalho and Hobbs 2012;Savage et al. 2015). Variation in structure of green spaces within cities represents the availability of habitats for arthopods in gradients of urbanization, however, green spaces within cities that present complex structures and higher plant richness are thought to be of high quality as habitats for insects (Pauleit and Duhme 2000;Whitford et al. 2001;Turner et al. 2005).
Here, we report the results of a study asking a) if the relative performance of aphid colonies (i.e. aphid population numbers) was associated with urbanization; b) if the presence of natural enemies (insect predators, parasitoids) and mutualists (ants) found on colonies was determined by urbanization or aphid numbers; c) how biotic factors (the assemblage of natural enemies and mutualists, green areas, plant species diversity and aphid numbers) and abiotic factors (impervious surfaces, distance from urban centre) act in concert to determine herbivore population sizes and the occurrence of their mutualists and natural enemies.

Study sites and habitat variables
Study sites were located in Greater Reading, Berkshire (51°27'N, 0°58'W), a large town in southern England with a population of 290,000, which covers an area of ca. 72 km 2 . Twenty-eight experimental sites were selected, and sites were selected in order to capture an approximate gradient from very urbanized environments in the town center to suburban areas located on the south, covering areas of carparks, churchyards, parks, private and community gardens and woodlands. Each study site was at least 110 m apart (Fig. 1).
Land-use data for urbanization metrics were derived from the Ordnance Survey MasterMap® Topography layer, which represents topography at a scale of 1:1250. This is subdivided into a number of themes: administrative boundaries, buildings, heritage and antiquities, land, rail, roads, tracks and paths, structures, terrain and height and water. Using GIS techniques, 30-m-radius buffers were delimited from the sites where the experimental plants were located. This buffer size was chosen due to restrictions in access, and also due to limitations in sampling effort for the estimation of plant diversity.
Reclassification of urbanization metrics was performed to give the proportions of area represented by the following habitat types within those buffers: green areas, which was composed of gardens and lawns with ornamental plants, bushes, trees and shrubs; impervious surfaces, which comprised of buildings (any building or artificial structure made of concrete, brick or stone) and byways (roads, roadsides, tracks or paths made of impervious surfaces such as asphalt). This procedure was carried using QGIS 2.8.1 (QGIS Development Team 2015). In addition to these habitat variables, plant species richness within a 30-m radius of the study sites was estimated during the experimental period by visually counting all plant morphospecies within the area surrounding each experimental colony. This method is strongly correlated with species richness and it effectively captures variance between study areas, with the advantage of reduced sampling effort and increased effectiveness to achieve statistical power (Abadie et al. 2008;Schmiedel et al. 2016). Distance to the urban center (m) was calculated from each study site to a point at the central area of the town (Fig. 1).

Study system and summer recording
Black bean aphids Aphis fabae Scopoli were maintained in a monoclonal culture in the laboratory using plastic and mesh cages. Cultures were kept at a constant temperature of 20 ± 1°C and 16:8 h L:D light regime at ambient humidity on broad bean, Vicia faba L. (var. the Sutton dwarf). Three days before being allocated to the study sites, three adults were transferred from the culture and reared on 14-16 day old dwarf broad bean plants (18-22 cm in height) to allow new colonies to become established. These plants were sown in pots with traditional potting compost (Vitax Grower, Leicester, England), and watered as required. After three days, the established aphid colonies on broad bean plants were transferred to the 28 study sites.
Two days after experimental colonies were placed in the field, species and numbers of aphids, predators, ants and parasitized aphids (mummies) were recorded, and then subsequently every three days for five recording days providing a total of 17 days of sampling in the field. At the end of this sampling period the plant/aphid-colonies were removed and replaced by new ones in the field. Sampling was repeated four times in 2015 (sampling period one: May 16, 20, 24, 28 and June 1; period two: June 15, 19, 23, 27 and July 1; period 3: July 16, 20, 24, 28 and August 1; period four: August 14, 18, 22, 26 and 30).

Data analysis
All statistical analyses were carried out using R 3.1.2 (R Development Core Team 2014).
The dataset consisted of the cumulative numbers of predators, ants, aphids and aphid mummies of the five counting events on each of the four sampling periods. Some colonies were lost during the four sampling periods (three colonies on the first sampling period, eight colonies on the second sampling period, three colonies on the third sampling period and four colonies on the fourth sampling period), caused by poor plant health, herbivory of plants by snails and slugs, and also from damage or theft by the public. This resulted in 94 observations for analysis.
All counts of aphids, predators, ants, mummies and plant richness were either log-transformed or square root-transformed to deal with extreme values and to standardize and homogenize residuals (Crawley 2007;Zuur et al. 2009). To analyze aphid colony numbers we used a linear mixed model fitted by reduced maximum likelihood using package nlme (Pinheiro et al. 2016), and as fixed factors (explanatory variables) we used proportion of impervious surfaces, plant richness, distance to the town center, predator abundance, number of ants farming the colony, parasitized mummies and an interaction factor between ants and predator numbers. We accounted for repeated sampling of colonies in sites through time by adding period as a random effect. We removed the variable proportion of green areas from the set of explanatory variables since it was highly correlated to the proportion of impervious surfaces (r = −0.92).
To deal with the excess of zeros when modelling ants, predators and parasitized mummies as response variables, we transformed these variables as factors (presence or absence) and ran logistic regressions models Fig. 1 Study site location (n = 28) in Greater Reading, England. The x marks the point from which distance from the city center was calculated from each sampling location. Aerial image was obtained from the Ordnance Survey Edina MasterMap® with a binomial error distribution family (with canonical link logit) using the function 'glmer' of package lme4 (Bates et al. 2015), with period as a random effect and fitted by maximum likelihood (Crawley 2007). When modelling predators we used the proportion of impervious surfaces, plant richness, distance to the town center, aphid abundance, number of ants farming the colony, and number of parasitized mummies as explanatory factors. When modelling ants as response variable, we used the proportion of impervious surfaces, plant richness, distance to the town center, predator abundance, aphid numbers, and numbers of parasitized mummies. When analyzing parasitized mummies as the response variable, we removed the first sampling period from the dataset since no mummies were found on this period (leaving 69 observations in total), then we modelled this as a function of the proportion of impervious surfaces, plant richness, distance to the town center, predator abundance, aphid numbers, and number of ants.
Model selection was made by comparing all candidate models using Akaike's Information Criteria (Burnham and Anderson 2003) by developing a series of alternative mixed effect models that include different combinations of the explanatory variables (Zuur et al. 2009 , Table 1), by fitting the full model with the set of all possible explanatory variables and taking out the least significant term on each step (Crawley 2007). We then ranked the models according to AIC differences (Δ i = AIC i -AIC min , where AIC i is the model i value and AIC min is the best model value). Models with Δ i < 2 provide substantial support for a candidate model, whereas values of Δ i between 4 and 7 provide less support and Δ i > 10 indicates that the model is unlikely. We also calculated Akaike weights for all models, where these model weights are used to indicate the importance of a model, with increasing weights indicating the likelihood of a particular model as the overall best model (Burnham and Anderson 2003). Aikaike weights can also be used to calculate the relative importance of a variable by summing the Akaike weights of all models that include that variable (Burnham and Anderson 2003).
We checked if collinearity could be a potential issue in our models through variance inflation factors (VIF) which is used as an indicator of multicollinearity in multiple regression, with VIF values higher than 3 indicating that covariation between predictors may be a problem (Zuur et al. 2007). All our VIF values were in the range of 1.34-2.94. All response variables were checked for spatial autocorrelation through spline correlograms on package ncf (Bjornstad 2015), and we did not find any significant spatial structure in the response variables. We assessed the validity of all models by checking normality, independence and homogeneity of model residuals.

Study sites
Our study sites captured an urban gradient. The proportion of impervious surfaces was negatively correlated (r = −0.46), and plant diversity positively correlated (r = 0.61), with distance from town center. Plant richness varied from 14 to 100 species (mean ± SE: 35.86 ± 3.42), proportion of impervious surfaces varied from 0 to 0.862 (mean ± SE: 0.425 ± 0.051) and proportion of green areas around study sites varied from 0.138 to 1 (mean ± SE: 0.526 ± 0.050) (Fig. 1).

Aphid abundance
Model selection based on AIC differences revealed three model candidates (Δ i < 2) for explaining variance on aphid numbers, the first with predators, ants and parasitoids; the second with predators and ants and the third only with numbers of ants farming aphid colonies (Table 1, models 1, 2, 3). However, Akaike weights indicated that the first and third models are more likely to be the best models for explaining aphid numbers (Table 2), with ants farming the aphid colony the variable of the highest importance, being positively correlated with aphid increase (0.885 based on the sum of Akaike weights within models with Δ i < 2) (Fig. 2).

Aphid predators
Based on AIC differences two models were selected as candidates for explaining the presence of predators on aphid colonies (Table 1, models 4 and 5); however since model 5 is 1.76 times more likely to be the best model (evidence ratio = 0.467/0.265) we chose this model as the overall best model. As explanatory factors, this model included proportion of impervious surfaces, which negatively determined predator presence; number of aphids, positively determining predator presence; and number of ants farming the colony, which negatively influenced the presence of aphid predators (Table 2, Fig. 3).

Ants farming aphid colonies
Three candidate models were selected based on AIC differences for explaining the presence of ants farming aphid colonies (Table 1, models 6, 7, 8). Model 7 (Table 2), with number of aphids, predators and parasitized mummies as explanatory factors, seemingly to be the best model due to its higher Akaike weight (0.434). However, Fig. 3 only shows the logistic regression curves for number of aphids (positive; Fig. 4a) and number of predators (negative; Fig. 4b) as explanatory factors as the number of parasitized mummies was not significant at α = 0.05 (Table 2, model 7).

Parasitoid attack
Two candidate models were selected for explaining the presence of parasitized aphids on the colonies: first with the Only models with substantial support (Δ i < 2) are shown, and models highlighted in bold are considered the best model candidates and are further described on Table 2. Δ i = AIC differences, calculated by subtracting the model with lowest AIC value from other model AIC values. W = Akaike weights, with higher weights indicating increased model importance; ImpSurf = proportion of impervious surfaces within 30 m buffers; PlantRich = plant species richness within 30 m; DistCentre = distance to each study site to the town centre numbers of aphids, colony-farming ants and distance to the town center, and second with the first two variables but without distance to the town center (Table 1, models 9 and 10). Since distance to the town centre was not significant in model 9 (Table 2) we considered the model with only the numbers of aphids and colony-farming ants as the best overall model (model 10 in Table 2, Fig. 5). Both variables were positively correlated with the presence of parasitized aphids on the experimental colonies.

Discussion
Our aim was to investigate how urbanization may affect the intensity and outcome of interactions between species at various trophic levels, using the black bean aphid, its natural enemies and ant mutualists as a model system. Overall, we found that the presence of mutualistic ants, predators and parasitoids varied as a function of aphid numbers on the plants.
Predators were the only group affected by abiotic factors, with fewer predators found in areas with increased proportions of impervious surfaces. The presence of mutualistic ants was associated with an increase in both aphid and parasitoid numbers, and a decrease in numbers of aphid predators. In no case did local plant diversity or distance to the urban center affect the abundance of any of the interacting species. We found that Aphis fabae colony size was not affected by abiotic variables, something expected as each colony remained on their study site for a limited amount of time (~20 days for each sampling period), feeding on plants previously sown under identical conditions. This allows us to consider interactions at higher trophic levels without the confounding effects of plant and prey quality. The abundance of predators was significantly affected by aphid colony size, the number of ants farming aphid colonies and the proportion of impervious surfaces in the habitat. Density dependence in predation is a widely recognized factor (Sinclair and Pech 1996;Hixon and Carr 1997;Anderson 2001;Arditi et al. 2001;Holbrook and Schmitt 2002;Hixon and Jones 2005). In our experiment, ants attending aphid colonies greatly reduced predator numbers. Previous studies have reported that honeydew-collecting ants can alter predator abundance Fig. 3 Representation of logistic regression model of (a) proportion of impervious surfaces, (b) number of aphids and (c) number of ants farming aphid colonies in predicting the presence (1) or absence (0) of predators in the colonies. Note that a multivariate logistic mixed effects model was used (see Methods), but the trend line for a logistic regression model for just one explanatory variable on each panel was used to illustrate the direction of relationship between variables. Log scale used on x axis of panels (b) and (c) Fig. 2 Relationship between abundance of aphids and numbers of ants farming the aphid colonies throughout the four sampling periods. Although linear mixed-effects models were performed (see Methods), the linear model trend line is shown to illustrate the direction of the significant relationship between variables. Note log scale used on y and x axes (James et al. 1999;Wimp and Whitham 2001;Kaplan and Eubanks 2002). Neither of the above factors was unexpected. However, we also show that increased urbanization, measured as the proportion of impervious surfaces surrounding the field sites, was associated with a reduction in the numbers of predators recorded.
Insect predators are relatively generalist, and their abundance will be associated with the local population size of a range of prey species. Given the reduction in native plant diversity and abundance in urbanized areas (Dreistadt et al. 1990;Burton et al. 2005;Williams et al. 2005;Williams et al. 2008;Isaacs et al. 2009;Walker et al. 2009), it would be surprising if predators were not sensitive to urbanization (McKinney 2006;Jones and Leather 2012;Otoshi et al. 2015). Urban management techniques such as treading, bird feeding, mowing and pesticide application negatively impact predacious beetles and hemipterans (Morris and Rispin 1987;Helden and Leather 2004;Orros and Fellowes 2012;Jones and Leather 2012;Orros et al. 2015;Bennett and Lovell 2014;Smith et al. 2015). Human-induced extinctions and local extirpations are often biased towards higher trophic levels (Pauly et al. 1998;Jackson et al. 2001;Duffy 2002;Byrnes et al. 2005), and that losses of even one or two species that belong to higher trophic levels can cause cascading effects on species present on basal trophic levels (Paine 2002;Schmitz 2003) and critically affect ecosystem processes (Tilman et al. 1997;Byrnes et al. 2005;Hooper et al. 2005).
Urbanized environments might affect organisms at higher trophic levels more than their hosts or prey, particularly when they exhibit higher levels of resource specialization (Tscharntke et al. 1998;Bailey et al. 2005;Pereira-Peixoto et al. 2016). In our study system, this may apply to insect predators but does not appear to affect the likelihood of colonies suffering parasitoid attack. However, there was an indication that parasitized mummies were less frequently found on more urbanized sites of the gradient (closer to the town center, Table 2, model 9), but this factor was not statistically significant. There have been studies which found negative correlations between parasitism and urbanization in a landscape context (Gibb and Hochuli 2002;Bennett and Gratton 2012;Calegaro-Marques and Amato 2014), which was not our objective in this work. The presence of physical barriers and  Panel (b) shows the logistic regression curve for number of ants farming the colonies predicting number of parasitized aphids. Note that a multivariate logistic mixed effects model was used (see Methods); however, the trend line for a logistic regression model for just one explanatory variable on each panel was used to illustrate the direction of relationship between variables. Log scale used on x axis of both panels structures like buildings and roads in cities may make insect dispersal problematic, and present an obstacle for breeding and foraging (Wratten et al. 2003;Raupp et al. 2010;Peralta et al. 2011). On the other hand, vegetated areas bordering roads, pavements and streets may serve as biological corridors, particularly those that maintain higher plant diversity and density (Haddad et al. 2003;Peralta et al. 2011).
Although we did not directly measure the functional traits of each trophic guild and how this relates to variance in habitat, our results may suggest a certain degree of sensitivity to urbanization possibly associated with differential dispersal abilities and habitat requirements of predators and parasitoids. Predators were composed of several taxa (spiders, hoverfly larvae, aphid midges, flower bugs and ladybirds), presenting a range of dispersal abilities and with a variety of dietary breadths (Rotheray 1989). Aphid parasitoids are known to use a range of visual, acoustic or olfactory cues to locate potential host patches, including long-range olfactory cues originating from the host plant Vandermoten et al. 2012). This suggests that the differences in dispersal ability and habitat requirements may be responsible for the differences in vulnerability to urbanization between the two groups (predators and parasitoids) presented by our data. This possibility needs to be further tested.
We found that the mutualistic relationship between aphids and ants was responsible for a significant increase in aphid numbers. In our study, ant attendance at aphid colonies was not affected by habitat variables, and ant-attended colonies were present even on the most urbanized sites of the gradient. Mutualistic ants of aphids are known to protect aphid colonies from predator attack, to prevent mold growth when honeydew accumulates and to avoid aphid competition with other herbivores on the same resource (Stadler and Dixon 1998;Kaneko 2003;Yao 2014). The relationship between aphids and tending ants can then confer direct benefits to aphid survival, allowing highest feeding rates and nutrient uptake. At the same time, aphid-derived honeydew constitutes a nutrientrich food that may be essential for the survival and growth of ant colonies (Kaplan and Eubanks 2002;Tegelaar et al. 2013). Aphid parasitoids are less likely to be affected by the presence of ants on aphid colonies than predators. Although parasitoid wasps can sometimes be repelled by ants, once wasps successfully oviposit in aphids, these parasitized aphids also receive ant protection, which may in turn result in higher parasitoid emergence rates (Völkl 1992(Völkl , 1997Kaneko 2002;Yao 2014). Such patterns (a negative effect of ant presence on generalist predators, a positive effect on specialist enemies) was found by Wimp and Whitham (2001), who examined the mechanisms that determined arthropod community structure in a riparian zone dominated by cottonwood. Our urban ecosystem seemed to show similar trends.
Overall, only predators were affected by the features of urbanization measured on our study. This influence of environmental disturbance on multi-trophic interactions in urban habitats could result in important consequences for the assembly of local ecological communities, and also direct and practical implications for biocontrol services that natural enemies could provide in these habitats (Gibb and Hochuli 2002;Eubanks and Finke 2014;Calabuig et al. 2015;Philpott and Bichier 2017). For example, Turrini et al. (2016) investigated the effects of urbanization on trophic interactions and found that predators reduced aphid abundance less in urban than in agricultural ecosystems. This reduction in top-down regulation in urban areas resulted in urban plants having reduced biomass than plants in adjacent agricultural areas. Findings such as these emphasize that urbanization can influence not only interactions at higher trophic levels, but that these changes also affect plant communities through trophic cascades (Schemske et al. 1994;Brudvig et al. 2015). Our results highlight the negative effect of the main characteristic of cities, the increase amount of impervious surfaces, on an important trophic guild. Given that the amount of impervious surfaces is highly negatively correlated with proportion of green areas, our results reinforce the importance of maintaining and increasing the quality of urban green spaces as habitats for the conservation of biological diversity (Botkin and Beveridge 1997;Peralta et al. 2011), and consequently also on trophic dynamics.
One of the major challenges of ecology is to understand and predict the consequences of environmental changes for biodiversity and ecosystem functioning (van der Putten et al. 2004;Hooper et al. 2005). Variation in responses within and between trophic groups may cause restructuring of communities through changes in competitive, bottom-up and top-down control effects (Van der Putten et al. 2004). Any given species is affected by interactions with other species, therefore understanding how changes in species interactions potentially affect food web structure and function in urban habitats may help us to succeed when planning conservation strategies (Faeth et al. 2005;Faeth et al. 2011). To our knowledge, our work presents the first effort to address how interrelated multi-trophic interactions composed by herbivory, predation, parasitism and mutualism behave in urban habitats, with predation the most affected by the increase of urban features in the habitat. Our findings emphasize the need for careful consideration of how patterns of species interactions may be modified in urban settings, which is essential for conservation efforts that will promote ecosystem services and functioning in cities.