Skip to main content

Scale-dependent effects of habitat area on species interaction networks: invasive species alter relationships

Abstract

Background

The positive relationship between habitat area and species number is considered a fundamental rule in ecology. This relationship predicts that the link number of species interactions increases with habitat area, and structure is related to habitat area. Biological invasions can affect species interactions and area relationships. However, how these relationships change at different spatial scales has remained unexplored. We analysed understory plant–pollinator networks in seven temperate forest sites at 20 spatial scales (radius 120–2020 m) to clarify scale-associated relationships between forest area and plant–pollinator networks.

Results

The pooled data described interactions between 18 plant (including an exotic) and 89 pollinator (including an exotic) species. The total number of species and the number of interaction links between plant and pollinator species were negatively correlated with forest area, with the highest correlation coefficient at radii of 1520 and 1620 m, respectively. These results are not concordant with the pattern predicted by species–area relationships. However, when associations with exotic species were excluded, the total number of species and the number of interaction links were positively correlated with forest area (the highest correlation coefficient at a radius of 820 m). The network structure, i.e., connectance and nestedness, was also related to forest area (the highest correlation coefficients at radii of 720–820 m), when associations with exotics were excluded. In the study area, the exotic plant species Alliaria petiolata, which has invaded relatively small forest patches surrounded by agricultural fields, may have supported more native pollinator species than initially expected. Therefore, this invasive plant may have altered the original relationships between forest area and plant–pollinator networks.

Conclusions

Our results demonstrate scale-dependent effects of forest area on the size and structure of plant–pollinator networks. We also suggest that a single exotic plant species can impact plant–pollinator networks, even in temperate continental habitats.

Background

The relationship between species number and island area, namely, that species number increases with increasing island area, is a fundamental rule in ecology[14]. MacArthur and Wilson[1] hypothesized that this relationship results from a dynamic equilibrium between opposing immigration and extinction rates, which depend on island isolation and size, respectively. Because the same trend occurs in continental environments[13], the theory of island biogeography has been applied to the conservation of continental habitats[58].

Species–area relationships suggest that the number of interaction links among species, such as prey–predator and plant–pollinator species, increases with island or habitat area[9, 10]. However, surprisingly, only a few studies have tested this prediction, including one that examined plant–ant interactions on oceanic islands[9], and another that studied plant–pollinator interactions on continental habitats[11]. Additionally, the network structure of species interactions is empirically known to relate to the total number of species interacting[1214], which raises the expectation that there may be a close relationship between area and network structure[9]; however, only one study has tested this prediction at the landscape scale[15].

Species–area relationships differ among spatial scales, with the shapes and slopes of the relationships differing among local, regional, and global scales[3, 1618]. In addition, the effects of spatial scale on species–area relationships may depend on taxonomic group; for example, responses to habitat area differ between animals and plants[3, 19]. Given that networks of species interactions usually involve different taxonomic groups[20, 21], different spatial scales may affect the relationships between these networks and habitat area. However, how habitat area influences interaction networks at different spatial scales at the landscape level has remained unexplored.

Biological invasions impact the interactions among native species[22, 23], and island communities are more likely to be invaded and affected by exotic species than continental ones[2426]. However, exotic species may invade continental communities disturbed by human activities[14, 24]. Therefore, exotic species may impact the original relationship between habitat area and interaction networks even in continental environments.

Here we analysed plant–pollinator networks, both including and excluding exotic species, in temperate continental forests at different spatial scales to clarify scale-variable effects of habitat area on interaction networks. Plant–pollinator interactions provide excellent model systems for investigating the structure of species interactions[12, 2729]. Plant–pollinator interactions have recently been targeted for studying the network structure of plant–animal mutualistic interactions[2830]. For example, network metrics such as connectance and nestedness are often used to clarify the structure of mutualisitic interaction networks[1214, 31, 32]. As network metrics are correlated with the total number of interacting species[1214, 32, 33], and the number of species is related to habitat area[13], we expect that network structure is also related to habitat area. Reports of such area-related effects on network structure are rare[9], although the mathematical consequences of changing species numbers for network structure have been discussed[31, 34].

The main goal of this study was to test the following hypotheses: (1) the total number of species increases with forest area, (2) the network metrics of plant–pollinator interactions are related to the total number of species and forest area, (3) the relationships between forest area and plant–pollinator networks are scale-dependent, and (4) exotic species impact the relationships between forest area and plant–pollinator networks. Although the effects of island and habitat area on interaction networks have recently been reported[9, 11, 15], the scale-dependencies have never been examined.

Methods

Study system

We analysed a data set that was in part used by Taki and Kevan[35] to examine the effects of forest loss on the degree of specialization/generalization for plant and insect communities. Our analyses focused on the effects of forest area on the size and structure of the whole network in terms of scale-dependency and biological invasion. The study was conducted in Norfolk County, Ontario, Canada ( Additional file1; 42°37'–42°48' N, 80°25'–80°39' W; ca. 200 m a.s.l.). The study forests were located in southern Ontario’s deciduous Carolinian forest zone. The landscape is flat and characterized by distributed fragments of forest within intensively managed agricultural fields of crops such as corn, soybean, and tobacco. The forests were composed of deciduous trees such as oak and maple[36]. Small shrubs and herbs dominated the forest floors. The associations of flower-visiting insects with understory plants that flowered prior to leaf flush of the forest crown trees were examined (Figure 1) by observing the interactions between 18 flowering plants and 89 flower-visiting insect species ( Additional file2 and Additional file3). Although it is possible that not all of the flower-visitors were legitimate pollinators, we considered them to be pollinators to decipher the general patterns between habitat area and plant–pollinator networks. Garlic mustard Alliaria petiolata and honeybee Apis mellifera, both of which originated in Europe, are exotic species in this region ( Additional file2 and Additional file3). Both species are widely acknowledged to impact native communities[37, 38]. However, it is unlikely that A. mellifera has become established in the study forests, and the small numbers of individuals captured in this study likely escaped from domesticated colonies. Alliaria petiolata, on the other hand, has become established in these areas, and it invades forest edge and understory.

Figure 1
figure 1

Andrena bee visiting a flower of Claytonia virginica (left) and Trillium grandiflorum (right) in the forest understory.

Study design

The study was conducted in seven forest sites[35] ( Additional file1), which were selected using a Geographical Information System (GIS) with ArcView (Version 3.3, ESRI, Redlands, CA, USA). Seven geographical points that fell within the forest polygons in the study region were randomly chosen. The criteria for accepting a selected point included a minimum distance of 40 m from all edges of the forest polygons and a minimum separation distance of 4500 m from any other chosen geographical points. The geospatial data of forest coverage were produced using aerial photography (1:30000 and 1:50000) obtained from the Ontario Base Map Series in 2003 (Ontario Ministry of Natural Resources, Peterborough, ON, Canada).

Each of the seven selected sites consisted of a hexagonal transect with 20 m sides with the chosen geographical point marking the centre ( Additional file4). In constructing the hexagonal transects, the direction of the first radial arm of each transect was randomly chosen using a 1.5 m stick thrown into the air. The axes were marked with bamboo poles and a 120 m section of rope demarcated the perimeter. Although a better choice for transect shape would be a circle, hexagonal transects were more practical.

A belt transect method[39] was applied to the hexagonal transects. We sampled before canopy closure because most understory plant species bloom during this short season[4042]. Flower-visiting insects were sampled on sunny days when the temperature was at least 12.8°C. Sampling started at 11:30 and at 14:30, times at which flower visitations by insects are relatively frequent. Two sites were sampled each day in most cases. Four sampling cycles were conducted to ensure that sampling was carried out twice at 11:30 and twice at 14:30 for each of the seven sites (i.e., each site was sampled on four different days). All samplings were conducted from late April to late May. At each sampling, the same two researchers walked the perimeter of the hexagons five times at a slow pace. One person walked clockwise and the other walked counter-clockwise. Each of these samplings took 80–100 min. All insects visiting the flowers within 2 m on either side of the perimeter rope were sampled with insect nets and aspirators. The species of the flowering plants visited by each sampled insect were recorded. The species richness and abundance of all flowering plants in bloom at each site were recorded using a quadrate method with 1 x 1 m quadrates placed every meter along each side of the perimeter of each hexagon, creating 20 quadrates along the outer perimeter of each side and 19 along the inner row (234 quadrates in total for each recording). The cataloguing and census of plants took place four times for each site, coinciding with the first, second, third, and fourth insect sampling cycle (described above).

To quantify the amount of forest area at different spatial scales, circles of 20 radii from 120 to 2020 m (100–2000 m from the hexagonal transects) with 100 m intervals were created using ArcView on maps around each of the hexagonal transects. The scales followed previous studies on foraging ranges and scale-dependent effects in flower-visiting insects such as bees[4346]. The forest area at each of the scales was estimated by the amount of forest coverage (km2) within the circles. The environment surrounding the forests, i.e., intensively managed annual agricultural fields, was not habitat for most of native plant or pollinator species recorded in the study forests. Therefore, we can consider the forest coverage as habitat area for most plant and pollinator species (except for honeybees).

The degree of connectivity among forest patches has been considered one of the most important factors for evaluating the impacts of forest fragmentation on biodiversity [e.g.,[47]. In this study area, Taki et al.[46] showed that bee abundance and diversity increased with increasing forest area but not increasing forest edge. This suggests that connectivity has a lesser influence on plant-pollinator networks than the overall forested area in this study region. Therefore, our analyses did not focus on the effects of forest connectivity but rather the effects of forest area on plant-pollinator networks.

Data analysis

The numbers of species and interaction links were pooled at each forest site for the analyses. To determine the impacts of exotic species on the relationships between forest area and plant–pollinator networks, two data sets were analysed; the first included exotic species (i.e., contained all species) and the second excluded associations with exotic plant and pollinator species (i.e., contained only native species).

The numbers of plant species, pollinator species, total species, and occupancy by exotic garlic mustard A. petiolata were analysed in relation to forest area using Pearson’s correlation coefficient at each spatial scale. Prior to the analysis, forest area and species numbers were log transformed; occupancy of flowering A. petiolata (proportion of maximum numbers of quadrates out of 234) was square root-arcsine transformed.

To examine how the network metrics were related to the total number of species and forest area, the number of interaction links, connectance, and nestedness were calculated. On the basis of previous studies[9, 11], the number of observed interaction links between plant and pollinator species (I) can be predicted to increase with increases in the total number of species and habitat area.

Connectance (C), or the proportion of interactions actually observed amongst all possible interactions, usually represents the degree of redundancy in a system, with consequences for community stability[12, 48]. It was calculated as follows[12, 31, 48]:

C = I / ( P l × P o ) ,
(1)

where Pl and Po are the numbers of plant and pollinator species, respectively. On the basis of previous studies[12, 14, 32, 33], connectance can be predicted to decrease with increases in the total number of species and habitat area.

Nestedness is a pattern composed of asymmetrical interactions between generalists and specialists and symmetrical interactions among generalists. Nestedness is frequently detected within plant–pollinator networks[13, 29, 30, 49], although the concept of nestedness was originally used to analyse patterns of species occurrence and absence on a set of islands or in habitat fragments[50]. Nested networks are characterized by specialist species that tend to interact only with generalists; generalists that all interact with each other, forming a core of interacting species; and the absence of specialists that interact only with other specialists[28, 30]. We used the software ANINHADO 3.0.3[51, 52] to calculate nestedness (NODF; range 0–100). Relatively large values of NODF indicate a high degree of nestedness. Two null models (Models I and II provided by ANINHADO 3.0.3) were used to test the degree of nestedness expected from the basic network features. The first null model assumes that each randomly assigned plant and pollinator pair interacts with a constant probability, C (connectance). Therefore, it tests whether the observed NODF is higher than expected for random networks with a similar number of interactions. The second null model assumes that the probability of a plant, i, interacting with a pollinator, j, depends on the observed number of interactions of both species, such that

C ( r ij = 1 ) = k i P o + k j P l 1 2
(2)

where k is the observed number of interactions for the species, Pl is the number of plant species (rows) and Po is the number of pollinator species (columns). Therefore, null model II tests whether the observed NODF is higher than expected for random networks with similar heterogeneities of species interactions. Each community was compared to 1000 replicates generated by each null model. On the basis of the previous study[13], nestedness in mutualistic networks is predicted to increase with increases in the total number of species and habitat area.

The network metrics, link numbers, connectance, and nestedness were also analysed in relation to forest area using Pearson’s correlation coefficient at each spatial scale (JMP v. 7.0; SAS Institute, Cary, NC, USA). Before the analysis, C and NODF/100 were square root-arcsine transformed while link numbers and forest area were log transformed. Linear regression models were used to examine the relationships between the total number of species and network metrics (JMP v. 7.0). We also used linear regression models to investigate the relationships between forest area and network metrics. As the scale for linear regressions, we used the radius at which the highest correlation coefficient was detected.

Results

Relationships between forest area and species numbers

Among study sites, the numbers of plant species, pollinator species, and total species ranged from 3–12, 10–40, and 14–49, respectively (Table 1). The relationships between forest area and numbers of species were dependent on the spatial scale, and differed among plant and pollinator species (Figure 2a). The number of plant species was not related to forest area (at a radius of 820 m with the highest correlation coefficient; r2 = 0.36, F1,5 = 2.75, p = 0.158; Figures 2a,3a), while the number of pollinator species and total species were negatively correlated with forest area, with the highest correlation coefficient at radii of 1420 and 1520 m, respectively (Figure 2a,3c,3e). These results are not concordant with the patterns predicted by species–area relationships. Occupancy by exotic A. petiolata was negatively correlated with forest area, with the highest correlation coefficient at a radius of 1320 m (Figure 4). When exotic species, and the native species only observed interacting with exotics, were excluded, the numbers of plant species, pollinator species and total species were positively correlated with forest area, with the highest correlation coefficient at a radius of 820 m (Figures 2b,3b,3d,3f).

Table 1 Number of species and network metrics of plant–pollinator interactions in each study forest
Figure 2
figure 2

Correlation coefficients between forest area and numbers of species at each radius (range: 120–2020 m) among the study forests: (a) including exotic species, (b) excluding exotic species and the native species only observed interacting with exotics. Correlation coefficients in the shaded grey areas are statistically significant at p = 0.10 level.

Figure 3
figure 3

Relationship between forest area and (a, b) number of plant species, (c, d) number of pollinator species and (e, f) total number of species. Data for the radius with the highest correlation coefficients (Figure 2) regarding the effects of forest area are shown. Left (a, c, e) and right figures (b, d, f) are based on data including exotic species and excluding associations with exotics, respectively. Code numbers next to the circles indicate study site codes, which are provided in Table 1. Lines represent linear regressions: (b) number of plant species (PLS), log(PLS) = 1.8062 + 0.8348*log(FA), r2 = 0.49, F1,5 = 4.77, p = 0.081; (c) number of pollinator species (POS), log(POS) = 4.3847 – 1.3624*log(FA), r2 = 0.63, F1,5 = 8.60, p = 0.033; (d) log(POS) = 2.6504 + 1.3926*log(FA), r2 = 0.58, F1,5 = 6.88, p = 0.047; (e) total number of species (TS), log(TS) = 4.5598 – 1.1211*log(FA), r2 = 0.51, F1,5 = 5.29, p = 0.070; (f) log(TS) = 3.0195 + 1.1740*log(FA), r2 = 0.56, F1,5 = 6.30, p = 0.054.

Figure 4
figure 4

Relationship between forest area and occupancy by the exotic plant Alliaria petiolata. (a) Correlation coefficient between forest area and occupancy by A. petiolata at each radius (range: 120–2020 m) among the study forests. (b) Relationship between forest area at 1320 m radius and occupancy by A. petiolatai. Data for the radius with the highest correlation coefficients regarding the effects of forest area are shown. Code numbers next to the circles indicate study site codes, which are provided in Table 1. Correlation coefficients in the shaded grey areas are statistically significant at p = 0.10 level. Line represents a linear regression: occupancy by A. petiolata (AP, proportion of maximum numbers of quadrates out of 234), arcsin A P = 0.6452 – 0.5800*log(FA), r2 = 0.56, F1,5 = 6.46, p = 0.052.

Relationships between species numbers and interaction networks

The numbers of interaction links between plant and pollinator species ranged from 11–62 among the study sites (Table 1). Of these, the percentage involving exotic species ranged from 0–94.6% among the study sites (Table 1). The numbers of links were positively correlated with the total number of species (Figure 5a). This pattern did not change when associations with exotic species were excluded from the analysis (Figure 5b). Furthermore, connectance was not related to the total number of species when exotics were included (C, including exotics, 0.105–0.343; r2 = 0.10, F1,5 = 0.53, p = 0.50; Figure 5c). However, when associations with exotic species were excluded from the analysis, connectance was negatively correlated with the total number of species (C, excluding exotics, 0.103–0.500; Figure 5d). Most of the plant–pollinator networks were not nested (NODF, including exotics, 4.74–28.79; excluding exotics, 0–28.79; Table 1). The degree of nestedness was not related to the total number of species regardless of whether the associations with exotic species were excluded from the analysis (including exotics, r2 = 0.06, F1,5 = 0.33, p = 0.59; excluding exotics, r2 = 0.31, F1,5 = 2.21, p = 0.20; Figures 5e,5f).

Figure 5
figure 5

Relationship between total number of species and network metrics for (a, b) number of interaction links, (c, d) connectance, and (e, f) nestedness. Left (a, c, e) and right (b, d, f) panels are based on data including exotic species and excluding associations with exotics, respectively. Code numbers next to the circles indicate study site codes, which are provided in Table 1. Lines represent linear regressions: (a) number of interaction links (I), log(I) = −0.9601 + 1.2570*log(S), r2 = 0.96, F1,5 = 118.52, p = 0.0001; (b) log(I) = −1.0341 + 1.2774*log(S), r2 = 0.99, F1,5 = 334.95, p < 0.0001; (d) arcsin C = 1.0479 – 0.1894*log(S), r2 = 0.94, F1,5 = 74.16, p = 0.0003.

Relationships between forest area and interaction networks

The numbers of interaction links were negatively correlated with forest area, with the highest correlation coefficient at a radius of 1620 m (Figures 6a,7a), which is not concordant with the patterns predicted by species–area relationships. However, when associations with exotic species were excluded from the analysis, the number of interaction links was positively correlated with forest area, with the highest correlation coefficient at a radius of 820 m (Figures 6b,7b). Connectance was not related with forest area (at a radius of 820 m with the highest correlation coefficient; r2 = 0.34, F1,5 = 2.59, p = 0.169; Figures 6a,7c). However, when associations with exotics were excluded from the analysis, connectance was negatively correlated with forest area (highest correlation coefficient at the radius of 820 m; Table 1, Figures 6b,7d). Nestedness was positively correlated with forest area regardless of whether the associations with exotic species were excluded from the analysis, with the highest correlation coefficients at radii of 520 m including exotics and 620–720 m excluding exotics (Table 1, Figures 6,7e,7f).

Figure 6
figure 6

Correlation coefficients between forest area and network metrics (number of interaction links, connectance, and nestedness) at each radius (range: 120–2020 m) among the study forests: (a) including exotic species, (b) excluding associations with exotic species. Correlation coefficients in the shaded grey areas are statistically significant at p = 0.10 level.

Figure 7
figure 7

Relationship between forest area and network metrics for (a, b) number of interaction links, (c, d) connectance and (e, f) nestedness. Data for the radius with the highest correlation coefficients (Figure 6) regarding the effects of forest area are shown. Left (a, c, e) and right (b, d, f) figures are based on data including exotic species and excluding associations with exotics, respectively. Code numbers next to the circles indicate study site codes, which are provided in Table 1. Lines represent linear regressions: (a) number of interaction links (I), log(I) = 5.0660 – 1.52333*log(FA), r2 = 0.56, F1,5 = 6.32, p = 0.054; (b); log(I) = 2.8264 + 1.5556*log(FA), r2 = 0.59, F1,5 = 7.23, p = 0.043; (c) connectance (C), arcsin C = 0.4574 – 0.1231*log(FA), r2 = 0.34, F1,5 = 2.59, p = 0.169; (d) arcsin C = 0.4758 – 0.2245*log(FA), r2 = 0.53, F1,5 = 5.69, p = 0.063; (e) nestedness (NODF), arcsin N O D F / 100 = 0.5002 + 0.5195*log(FA), r2 = 0.53, F1,5 = 5.57, p = 0.065; (f), arcsin N O D F / 100 = 0.3963 + 0.2718*log(FA), r2 = 0.71, F1,5 = 12.1, p = 0.018.

Discussion

Some of our results did not support our hypotheses when exotic species were included (Figures 3e,5c,e,6a,7a,c). However, almost all of the hypotheses were verified when associations with exotic species were excluded from the networks (Figures 3f,5b,d,6b,7b,d,f). In particular, we determined that the relationship between forest area and plant–pollinator network depended on spatial scale (Figures 2,6). To the best of our knowledge, this is the first study to demonstrate scale-dependent effects of habitat area and exotic species on interaction networks.

Scale-dependent effects of forest area on plant–pollinator networks

The species–area relationships were vastly different between plant and pollinator species when exotic species were included in the analyses (Figures 3a,c). Additionally, the spatial scale at which the highest correlation coefficient was found differed between plant and pollinator species (Figure 2a). Responses to habitat area are generally different between plant and pollinator species[19]. However, the scale-dependent relationship was similar among the number of plant species, pollinator species, and interaction links, when exotic species as well as native species only observed interacting with exotics were excluded (Figure 2b6b). The similarities among the relationships between native plant and pollinator species are not likely to have been caused by our sampling scheme (i.e., insect collection on plants), because the results using a different sampling method (pan traps) also showed that native bee abundance and diversity increased with increasing forest area (with the highest correlation coefficients at 500–750 m radius[46]). The study[46] was conducted in the same region, but in different study forests (where A. petiolata did not flower). Therefore, the similar responses of native plant and native pollinator species to forest area caused a positive relationship between forest area and the numbers of links at the same spatial scale (i.e., 820 m; Figures 6b7b).

Previous studies have indicated that network metrics, connectance and nestedness are related to the total number of species involved in mutualistic interactions[1214, 32, 33]. Our results partly supported this pattern (Figures 5a,b,d; but Figure 5c,e,f). Although previous studies have used data sets composed of mutualistic networks at different sampling areas as well as different geographical regions[1214, 32, 33], all of our data were collected from within the same sampling area (234 m2) in the same region. Because the total number of species was related to the forest area (Figures 3e,f), so was the network structure (Figures 7a,b,d,e,f). More redundant networks with highly asymmetric interactions were found in relatively large forest areas (Figure 7d,e,f), suggesting that the stability of plant–pollinator networks might increase with forest area. However, these relationships may partly be explained by non-biological factors. For example, connectance may decrease with increasing numbers of possible interaction links, as the number of observations per species declines when the same absolute effort is made to sample networks of different sizes[31, 34]. In this study, however, connectance did not decrease with the total number of species when exotic species were included (Figure 5c). Also, nestedness was not related to the total number of species (Figures 5e,f), although it increased with increasing forest area (Figures 7e,f). These results suggest that the effects of forest area on connectance or nestedness were caused not only by changes in the total number of species, but also by other factors. For example, changes in the population densities of some species might affect the structure of interaction networks relative to habitat area because the population densities of various species are known to be related to habitat area[53]. Further studies are needed to clarify the mechanism driving the relationships between habitat area and network structure.

Impacts of exotic species

The scale-dependent relationships among forest area, total number of species, and interaction networks were different when interactions were considered with and without exotic species (Figures 23567), suggesting impacts of exotic species. Plant–pollinator interactions included only two exotic species, garlic mustard A. petiolata and the honeybee A. mellifera (Additional file2 and Additional file3). Although the contribution of A. mellifera to the interactions was not insignificant (Additional fileA3), A. petiolata was central to the interactions in at least two forests ( Additional file2). Although exotic plants are only rarely thought to invade temperate natural forests[24, 54], A. petiolata has frequently been reported to invade the forest edge and understory in North America[38]. Flowers of A. petiolata produce rich nectar, attracting a variety of native bees and flies[55]. Nectar-rich flowers of invasive plants can disturb native plant–pollinator interactions[56]. In the study areas of the present study, A. petiolata, which has invaded relatively small forest patches surrounded by agricultural fields (Figure 4b), may have supported more native pollinator species than initially expected ( Additional file3).

The presence of exotic species strongly influenced the scale-dependent relationships between forest area and the number of interaction links and connectance (Figures 7a–d). Thus, this invasive plant may have altered the original relationships between forest area and plant–pollinator networks and their scale-dependency. However, nestedness showed the same trend for both networks, with and without exotics (Figures 7e,f). Vilá et al.[23] hypothesised that invasive plants by being supergeneralists, both interacting with generalists and specialists, would increase the nestedness of the plant–pollinator network. Our results showed that excluding associations with exotic species increased the values of nestedness in some study sites (Figures 7e,f). However, this increase did not change the relationships between forest area and nestedness (Figures 7e,f).

Conclusions

Since the original publication of the equilibrium theory of island biogeography[1], species–area relationships have been extensively studied for various groups of organisms[3, 4]. Furthermore, the theory has been applied to the conservation of focal species in continental habitats[4, 5]. Sugiura[9] analysed the relationships between island area and plant–ant mutualistic interactions and suggested an extension of this basic species–area relationship to more specifically address “species interactions–area relationships”. Interaction networks on continental habitats as well as oceanic islands can be considered with a view to the species interactions–area relationships. Indeed, Sabatino et al.[11] indicated that the number of interaction links among flowering plants and their pollinators increases with habitat area in the continental environment. Valladares et al.[15] also indicated that the network structure of plant–leafminer–parasitoid webs is related to fragmented forest area. In this study, we determined that network metrics of plant–pollinator interactions are related to forest area and the relationships depend on spatial scale.

Each species has a particular spatial scale of habitat area that most strongly affects the abundance because body size and mobility differ among species[19, 44]. Therefore, the spatial scales of habitat area that most strongly affected the abundance generally differed among species[19]. In this study, we determined the spatial scales that most strongly affected the structure of plant–pollinator interactions, which suggests the presence of a spatial scale that most strongly influences the structuring and maintenance of the species interaction network.

Finally, exotic species appear to alter the relationship between habitat area and interaction network. Although many researchers have reported that exotic species impact interaction networks, particularly on oceanic islands[9, 14, 5759], we suggest that a single exotic plant species can impact this relationship, even in temperate continental habitats.

References

  1. MacArthur RH, Wilson EO: The theory of island biogeography. 1967, Princeton: Princeton University Press

    Google Scholar 

  2. Connor EF, McCoy ED: The statistics and biology of the species–area relationship. American Naturalist. 1979, 5: 397-411.

    Google Scholar 

  3. Rosenzweig ML: Species diversity in space and time. 1995, Cambridge: The University of Cambridge

    Book  Google Scholar 

  4. Lomolino MV: Ecology’s most general, yet protean pattern: the species–area relationship. Journal of Biogeography. 2000, 27: 17-26. 10.1046/j.1365-2699.2000.00377.x.

    Article  Google Scholar 

  5. Diamond JM: The island dilemma: lessons of modern biogeographic studies for the design of natural reserves. Biological Conservation. 1975, 7: 129-146. 10.1016/0006-3207(75)90052-X.

    Article  Google Scholar 

  6. Quammen D: The song of the dodo: Island biogeography in an age of extinctions. 1996, New York, London, Toront, Sydney: Scribner

    Google Scholar 

  7. Laurance WF: Beyond island biogeography theory: understanding habitat fragmentation in the real world. 2010, Princeton University Press: In The theory of island biogeography revisited. Edited by Losos JB, Ricklefs RE. Princeton, NJ, 214-236.

    Google Scholar 

  8. He F, Hubbell SP: Species–area relationships always overestimate extinction rates from habitat loss. Nature. 2011, 478: 368-371.

    Article  Google Scholar 

  9. Sugiura S: Species interactions–area relationships: biological invasions and network structure in relation to island area. Proceedings of the Royal Society B. 2010, 277: 1807-1815. 10.1098/rspb.2009.2086.

    Article  PubMed  PubMed Central  Google Scholar 

  10. Burkle LA, Alarcón R: The future of plant–pollinator diversity: understanding interaction networks across time, space, and global change. American Journal of Botany. 2011, 98: 528-538. 10.3732/ajb.1000391.

    Article  PubMed  Google Scholar 

  11. Sabatino M, Maceira N, Aizen MA: Direct effects of habitat area on interaction diversity in pollination webs. Ecological Applications. 2010, 20: 1491-1497. 10.1890/09-1626.1.

    Article  PubMed  Google Scholar 

  12. Jordano P: Patterns of mutualistic interactions in pollination and seed dispersal: connectance, dependence asymmetrics, and coevolution. American Naturalist. 1987, 129: 657-677. 10.1086/284665.

    Article  Google Scholar 

  13. Guimarães PRJ, Rico-Gray V, dos Reis SF, Thompson JN: Asymmetries in specialization in ant–plant mutualistic networks. Proceedings of the Royal Society B. 2006, 273: 2041-2047. 10.1098/rspb.2006.3548.

    Article  PubMed  PubMed Central  Google Scholar 

  14. Aizen MA, Morales CL, Morales JM: Invasive mutualists erode native pollination webs. PLoS Biology. 2008, 6: e31-10.1371/journal.pbio.0060031.

    Article  PubMed  PubMed Central  Google Scholar 

  15. Valladares G, Cagnolo L, Salvo A: Forest fragmentation leads to food web contraction. Oikos. 2012, 121: 299-305. 10.1111/j.1600-0706.2011.19671.x.

    Article  Google Scholar 

  16. Losos JB, Schluter D: Analysis of an evolutionary species–area relationship. Nature. 2000, 408: 847-850. 10.1038/35048558.

    Article  CAS  PubMed  Google Scholar 

  17. Crawley MJ, Harral JE: Scale dependence in plant biodiversity. Science. 2001, 291: 864-868. 10.1126/science.291.5505.864.

    Article  CAS  PubMed  Google Scholar 

  18. Hubbell SP: The unified neutral theory of biodiversity and biogeography. 2001, Princeton: Princeton University Press

    Google Scholar 

  19. Steffan-Dewenter I, Tscharntke T: Insect communities and biotic interactions on fragmented calcareous grasslands—a mini review. Biological Conservation. 2002, 104: 275-284. 10.1016/S0006-3207(01)00192-6.

    Article  Google Scholar 

  20. Montoya JM, Pimm SL, Sole RV: Ecological networks and their fragility. Nature. 2006, 442: 259-264. 10.1038/nature04927.

    Article  CAS  PubMed  Google Scholar 

  21. Ings TC, Montoya JM, Bascompte J, Blüthgen N, Brown L, Dormann CF, Edwards F, Figueroa D, Jacob U, Jones JI, Lauridsen RB, Ledger ME, Lewis HM, Olesen JM, van Veen FJ, Warren PH, Woodward G: Ecological networks – beyond food webs. Journal of Animal Ecology. 2008, 78: 253-69.

    Article  Google Scholar 

  22. Memmott J, Waser NM: Integration of alien plants into a native flower-pollinator visitation web. Proceedings of the Royal Society of London B. 2002, 269: 2395-2399. 10.1098/rspb.2002.2174.

    Article  Google Scholar 

  23. Vilà M, Bartomeus I, Dietzsch AC, Petanidou T, Steffan-Dewenter I, Stout JC, Tscheulinm T: Invasive plant integration into native plant-pollinator networks across Europe. Proceedings of the Royal Society B. 2009, 276: 3887-3893. 10.1098/rspb.2009.1076.

    Article  PubMed  PubMed Central  Google Scholar 

  24. Elton CS: The ecology of invasions by animals and plants. 1958, London: Methuen

    Book  Google Scholar 

  25. Simberloff D: Why do introduced species appear to devastate islands more than mainland areas?. Pacific Science. 1995, 49: 87-97.

    Google Scholar 

  26. Corlett RT: Invasive aliens on tropical East Asian islands. Biodiversity and Conservation. 2010, 19: 411-423. 10.1007/s10531-009-9624-4.

    Article  Google Scholar 

  27. Lewinsohn TM, Prado PI, Jordano P, Bascompte J, Olesen JM: Structure in plant–animal interaction assemblages. Oikos. 2006, 113: 174-184. 10.1111/j.0030-1299.2006.14583.x.

    Article  Google Scholar 

  28. Bascompte J, Jordano P: Plant–animal mutualistic networks: the architecture of biodiversity. Annual Review of Ecology, Evolution, and Systematics. 2007, 38: 567-593. 10.1146/annurev.ecolsys.38.091206.095818.

    Article  Google Scholar 

  29. Vázquez DP, Blüthgen N, Cagnolo L, Chacoff NP: Uniting pattern and process in plant–animal mutualistic networks: a review. Annals of Botany. 2009, 103: 1445-1457. 10.1093/aob/mcp057.

    Article  PubMed  PubMed Central  Google Scholar 

  30. Bascompte J, Jordano P, Melián CJ, Olesen JM: The nested assembly of plant–animal mutualistic networks. Proceeding of the National Academy of Sciences USA. 2003, 100: 9383-9387. 10.1073/pnas.1633576100.

    Article  CAS  Google Scholar 

  31. Blüthgen N, Fründ J, Vázquez DP, Menzel F: What do interaction network metrics tell us about specialization and biological traits?. Ecology. 2008, 89: 3387-3399. 10.1890/07-2121.1.

    Article  PubMed  Google Scholar 

  32. Olesen JM, Jordano P: Geographic patterns in plant–pollinator mutualistic networks. Ecology. 2002, 83: 2416-2424.

    Google Scholar 

  33. Blüthgen N, Menzel F, Hovestadt T, Fiala B, Blüthgen N: Specialization, constraints, and conflicting interests in mutualistic networks. Current Biology. 2007, 17: 341-346. 10.1016/j.cub.2006.12.039.

    Article  PubMed  Google Scholar 

  34. Kenny D, Loehle C: Are food webs randomly connected?. Ecology. 1991, 72: 1794-1799. 10.2307/1940978.

    Article  Google Scholar 

  35. Taki H, Kevan PG: Does habitat loss affect the communities of plants and insects equally in plant-pollinator interactions? Preliminary findings. Biodiversity and Conservation. 2007, 16: 3147-3161. 10.1007/s10531-007-9168-4.

    Article  Google Scholar 

  36. Waldron G: Trees of the Carolinian forest: a guide to species, their ecology and uses. 2003, Toronto: Boston Mills Press

    Google Scholar 

  37. Traveset A, Richardson DM: Biological invasions as disruptors of plant reproductive mutualisms. Trends in Ecology and Evolution. 2006, 21: 208-216. 10.1016/j.tree.2006.01.006.

    Article  PubMed  Google Scholar 

  38. Stinson K, Kaufman S, Durbin L, Lowenstein F: Impacts of Garlic Mustard invasion on a forest understory community. Northeastern Naturalist. 2007, 14: 73-88. 10.1656/1092-6194(2007)14[73:IOGMIO]2.0.CO;2.

    Article  Google Scholar 

  39. Memmott J: The structure of a plant-pollinator food web. Ecology Letters. 1999, 2: 276-280. 10.1046/j.1461-0248.1999.00087.x.

    Article  Google Scholar 

  40. Macior LW: Pollination ecology of vernal angiosperms. Oikos. 1978, 30: 452-460. 10.2307/3543340.

    Article  Google Scholar 

  41. Schemske DW, Wilson MF, Melampy MN, Miller LJ, Verner L, Schemske KM, Best LB: Flowering ecology of some spring woodland herbs. Ecology. 1978, 59: 351-366. 10.2307/1936379.

    Article  Google Scholar 

  42. Helenurm K: Barrett SCH: The reproductive biology of boreal forest herbs. II. Phenology of flowering and fruiting. Canadian Journal of Botany. 1987, 65: 2047-2056.

    Article  Google Scholar 

  43. Osborne JL, Clark SJ, Morris JR, Williams IH, Riley JR, Smith AD, Reynolds DR, Edwards AS: A landscape-scale study of bumble bee foraging range and constancy, using harmonic radar. Journal of Applied Ecology. 1999, 36: 519-533. 10.1046/j.1365-2664.1999.00428.x.

    Article  Google Scholar 

  44. Steffan-Dewenter I, Münzenberg U, Bürger C, Thies C, Tscharntke T: Scale-dependent effects of landscape context on three pollinator guilds. Ecology. 2002, 83: 1421-1432. 10.1890/0012-9658(2002)083[1421:SDEOLC]2.0.CO;2.

    Article  Google Scholar 

  45. Gathmann A, Tscharntke T: Foraging ranges of solitary bees. Journal of Animal Ecology. 2002, 71: 757-764. 10.1046/j.1365-2656.2002.00641.x.

    Article  Google Scholar 

  46. Taki H, Kevan PG, Ascher JS: Landscape effects of forest loss in a pollination system. Landscape Ecology. 2007, 22: 1575-1587. 10.1007/s10980-007-9153-z.

    Article  Google Scholar 

  47. Fahrig L: Effects of habitat fragmentation on biodiversity. Annual Review of Ecology, Evolution and Systematics. 2003, 34: 487-515. 10.1146/annurev.ecolsys.34.011802.132419.

    Article  Google Scholar 

  48. May RM: Stability and complexity in model ecosystems. 1973, Princeton, NJ: Princeton University Press

    Google Scholar 

  49. Ulrich W, Almeida-Neto M, Gotelli NJ: A consumer’s guide to nestedness analysis. Oikos. 2009, 118: 3-17. 10.1111/j.1600-0706.2008.17053.x.

    Article  Google Scholar 

  50. Atmar W, Patterson BD: The measure of order and disorder in the distribution of species in fragmented habitat. Oecologia. 1993, 96: 373-382. 10.1007/BF00317508.

    Article  Google Scholar 

  51. Guimarães PRJ, Guimarães P: Improving the analyses of nestedness for large sets of matrices. Environmental Modeling and Software. 2006, 21: 1512-1513. 10.1016/j.envsoft.2006.04.002.

    Article  Google Scholar 

  52. Almeida-Neto M, Guimarães PR, Guimarães PRJ, Loyola RD, Ulrich W: A consistent metric for nestedness analysis in ecological systems: reconciling concept and measurement. Oikos. 2008, 117: 1227-1239. 10.1111/j.0030-1299.2008.16644.x.

    Article  Google Scholar 

  53. Connor EF, Courtney AC, Yoder JM: Individual–area relationships: the relationship between animal population density and area. Ecology. 2000, 81: 734-748.

    Google Scholar 

  54. Von Holle B, Delcourt HR, Simberloff D: The importance of biological inertia in plant community resistance to invasion. Journal of Vegetation Science. 2003, 14: 425-432.

    Article  Google Scholar 

  55. Cruden RW, McClain AM, Shrivastava GP: Pollination biology and breeding system of Alliaria petiolata (Brassicaceae). Bulletin of the Torrey Botanical Club. 1996, 123: 273-280. 10.2307/2996775.

    Article  Google Scholar 

  56. Lopezaraiza-Mikel ME, Hayes RB, Whalley MR, Memmott J: The impact of an alien plant on a native plant–pollinator network: an experimental approach. Ecology Letters. 2007, 10: 539-550. 10.1111/j.1461-0248.2007.01055.x.

    Article  PubMed  Google Scholar 

  57. Olesen JM, Eskildsen LI, Venkatasamy S: Invasion of pollination networks on oceanic islands: importance of invader complexes and endemic super generalists. Diversity and Distributions. 2002, 8: 181-192. 10.1046/j.1472-4642.2002.00148.x.

    Article  Google Scholar 

  58. Sugiura S, Yamaura Y, Makihara H: Biological invasion into the nested assemblage of tree–beetle associations on the oceanic Ogasawara Islands. Biological Invasions. 2008, 10: 1061-1071. 10.1007/s10530-007-9184-z.

    Article  Google Scholar 

  59. Padrón B, Traveset A, Biedenweg T, Díaz D, Nogales M, Olesen JM: Impact of alien plant invaders on pollination networks in two archipelagos. PLoS ONE. 2009, 4: e6275-10.1371/journal.pone.0006275.

    Article  PubMed  PubMed Central  Google Scholar 

Download references

Acknowledgements

We thank Long Point Conservation Authority personnel, M. Armstrong, J. DeCloet, J. Knack, F. Saul and V. Sinnave for allowing us to access our study sites. We also thank P. Kevan, C. Fuss and A. and J. Boone for support and field assistance. We thank O. Honnay and anonymous reviewers for helpful comments on an earlier version of the manuscript. While conducting the field work, HT was supported by a scholarship from the Rotary Foundation. Financial support was provided by the Forestry and Forest Products Research Institute.

Author information

Authors and Affiliations

Authors

Corresponding author

Correspondence to Shinji Sugiura.

Additional information

Competing interests

The authors declare that they have no competing interests.

Authors’ contributions

SS and HT generated the original idea. HT designed the research. HT collected field data. SS and HT analysed data. SS wrote the paper. All authors read and approved the final manuscript.

Electronic supplementary material

12898_2012_209_MOESM1_ESM.eps

Additional file 1:Seven study sites with 20 radii ranging from 120 m to 2020 m (100-2000 m from the hexagonal transects) used to obtain forest area. All sites were in Norfolk Country, Ontario, Canada. Shaded areas represent forests, and the bar indicates a 2000 m scale. The geospatial data were obtained from the Ontario Base Map Series in 2003 (Ontario Ministry of Natural Resources, Peterborough, Ontario, Canada), Code numbers indicate study site codes, which are provided in Table 1. (EPS 5 MB)

Additional file 2:List of plant species at each forest site.(DOC 67 KB)

Additional file 3:List of pollinator species at each forest site.(DOC 202 KB)

12898_2012_209_MOESM4_ESM.eps

Additional file 4:Hexagonal transect with 20 m sides, with the chosen geographical point marking the centre (top). The axes were marked with bamboo poles and a 120 m section of rope demarcated the perimeter (bottom). (EPS 4 MB)

Authors’ original submitted files for images

Rights and permissions

This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Reprints and permissions

About this article

Cite this article

Sugiura, S., Taki, H. Scale-dependent effects of habitat area on species interaction networks: invasive species alter relationships. BMC Ecol 12, 11 (2012). https://doi.org/10.1186/1472-6785-12-11

Download citation

  • Received:

  • Accepted:

  • Published:

  • DOI: https://doi.org/10.1186/1472-6785-12-11

Keywords