DOI: 10.13102/sociobiology.v62i1.56-64Sociobiology 62(1): 56-64 (March, 2015) Open access journal: http://periodicos.uefs.br/ojs/index.php/sociobiology ISSN: 0361-6525 Topological Structure of Plant-bee Networks in Four Mexican Environments Introduction One of many goals in ecology is to understand how different mechanisms shape ecological interactions within natural communities (Loreau et al., 2001; Dyer et al., 2010). Using a network approach several studies have described complex, hyper- diverse communities (Bascompte et al., 2003; Lewinsohn et al., 2006; Hagen et al., 2012; Bascompte & Jordano, 2014). Most studies have primarily focused on mutualistic interactions between plants and animals (mainly pollinators and seed dispersers) (e.g. Vázquez & Aizen 2004; Olesen et al., 2008; Dupont et al., 2009 Donatti et al., 2011; Mello et al., 2011). Furthermore, in tropical regions, reproduction (pollination, fruit and seed dispersal) of more than 90% of shrubs and trees depends on animals (e.g., insects, birds, bats) (Bascompte & Jordano, 2008). Abstract We analyzed the topological structure (e.g., links per species, connectance, core-periphery analyses, specialization, nestedness, modularity) of plant-bee interactions of four areas in Mexico. We used qualitative data (binary networks). Mexico exhibits a complex geography and community variation that can affect bee networks. Network architecture is variable within large spatial scales, thus our results should vary according to site characteristics (La Mancha and Totula in Veracruz, Carrillo Puerto in Quintana Roo, and the Tehuacan- Cuicatlan valley, in Puebla), type of vegetation, endemism, altitude, size of area sampled. Network topology varied among sites, and the presence of nested or modular patterns was analyzed for robustness to simulated species extinctions. The lowest species richness was recorded for the Quintana Roo site (15 plant, 25 bee species), and the highest for the Tehuacan-Cuicatlan valley site (88 plant, 231 bee species). There was a tendency to have more connected networks when species richness was low and networks with greater species richness had a higher number of interactions. The distribution of interactions differed between environments but not due to network size and all were significantly nested. The robustness to cumulative extinctions showed a different pattern at each site; the most robust network was at Carrillo Puerto, which also was the site with less species. Sites with more endemic species (e.g. Tehuacan) had more specialized interactions, and were more susceptible to extinction. Sociobiology An international journal on social insects VA Ramírez-Flores1, R Villanueva-Gutiérrez2, DW Roubik3, CH Vergara4, N Lara-Rodríguez5, W Dáttilo1, MEB Ferrer6, V Rico-Gray1 Article History Edited by Gilberto M. M. Santos, UEFS, Brazil Received 19 January 2015 Initial acceptance 20 February 2015 Final acceptance 22 February 2015 Keywords Bee, plant, network, endemism, Mexico. Corresponding author Dr. Victor Rico-Gray Instituto de Neuroetologia Universidad Veracruzana Calle Dr. Luis Castelazo s/n, Col. Industrial- Animas, Xalapa, Veracruz 91190, Mexico E-Mail: vricogray@yahoo.com Traditionally, plant-pollinator systems have been studied as small groups of species (Stephens et al., 2009; Santos et al., 2012). However, modern network tools allow simultaneous analysis of species over different spatial and/or temporal scales (Stephens et al., 2009; Lange et al., 2013). Such analysis allows the study of complex systems and mega-diverse communities (Jordano et al., 2009; Hagen et al., 2012). Although several communities defy analysis based on few, incomplete samples, as the behavior of the system is often very different from the sum of its parts (Jordano et al., 2009). For instance, when most interactions are not reached by a study, the ‘system’ remains, essentially, unknown. Nonetheless, the conservation and management of the studied communities, as well as of pollination services, can receive useful information through these efforts (Burkle &Alarcón, 2011). RESEARCH ARTICLE - BEES 1 - Universidad Veracruzana, Xalapa, Veracruz, Mexico; 2 - Colegio de la Frontera Sur, ECOSUR-Unidad Chetumal, Mexico; 3 - Smithsonian Institution, Washington, DC., U.S.A.; 4 - Universidad de las Américas, Puebla, Mexico; 5 - Universidad de Alicante, San Vicente del Raspeig (Alicante), Spain; 6 - Universidad Complutense, Madrid, Spain Sociobiology 62(1): 56-64 (March, 2015) 57 Understanding the ecological and evolutionary causes and consequences of spatial and/or temporal variation in plant/ floral visitor systems is very important when formulating basic and applied questions on the structure and function of a community, evolution of floral characteristics (e.g. Hernández-Yáñez et al., 2013), and development of optimal conservation strategies (e.g. Burkle & Alarcon, 2011). Although landscape heterogeneity may augment the number of interactive species, these may not constantly share the same habitat (Jordano et al., 2009; Clemente et al., 2013). For instance, Olesen and Jordano (2002) analyzed plant-pollinator networks from different geographical areas and found that interaction patterns co-varied with factors such as latitude, altitude and habitat isolation (islands). Also, plant species are more generalized (i.e. have more interactions) at higher latitude and in lowland habitats. In oceanic islands, for instance, networks are often specialized (i.e. have less interactions), pollinators are scarce, and a low proportion of animals exists, relative to plant species, which may add to the idea of higher specificity in the tropics that at higher latitudes (Olesen & Jordano 2002). Based on high species numbers and the abundance of certain species, bees are probably the insect group best-adapted to floral biology and flower visits. This group is essential for pollination and sexual reproduction of most flowering plants (including many cultivars) (Roubik, 1992; Michener, 2000). The latter is very important in terms of securing future food supplies (Deguines et al., 2014). Specific studies (reviewed by Hagen et al., 2012) of plant-pollinator networks consider all floral visitors of different taxonomic groups (e.g. bees, flies, butterflies, beetles, wasps, birds, and bats), and have evaluated how different factors (biotic, abiotic) structure networks. However, despite the general growth of knowledge on network topology, little is known about how plant- floral visitor networks are structured and how they vary over space when only one taxonomic group in considered. Mexico exhibits a complex geography and community variation, both of which should affect plant-bee networks. Thus, network architecture (e.g. links per species, connectivity, core-periphery analyses, specialization, nestedness) is likely change within the country so that results should vary according to site characteristics (e.g. type of vegetation, endemism, altitude, size of area sampled). Here we describe the topological structure of plant-bee interactions in four distinct areas in Mexico. In particular, we searched for the presence of modular patterns and analyzed the pattern of secondary extinctions and stability of the networks in the range sampled. Material and Methods Study area and field work Field work was accomplished in four Mexican sites: 1) Tehuacan-Cuicatlan valley, Puebla; 2) surroundings of Carrillo Puerto, Quintana Roo; 3) La Mancha, Veracruz; and 4) a traditional management coffee plantation in the montane forest in Totutla, Veracruz (Table 1). Due to the number of researchers involved, field work differed among sites, although entomological nets were used throughout. The surveys in the Tehuacan and Totutla sites were guided by C. Vergara, the survey in Quintana Roo was guided by Villanueva and Roubick, and the survey at the La Mancha site was guided by N. Lara and V. Rico-Gray. The number of different areas per site was given by the vegetation changes (see Table 1). All accounts of visits by potential bee pollinators were determined by observations of bees on flowers and not from previous knowledge of the interactions. In the Tehuacan-Cuicatlan valley the vegetation is a xerophytic shrubland (Rzedowski 1978), and Zavala-Hurtado (1982) described four vegetation types: shrubs with spines (‘matorral espinoso’), tetechera, cardonal and izotal (Yucca spp.). Even though some of the vegetation types are in Spanish, the latter names are based on dominant plant species (see also Table 1). The biota of the valley is characterized by a high number of endemic species, the highest density known in Mexico (Méndez- Larios et al., 2004). Sampling was accomplished with monthly observations (3-5 days per visit) between July 1996 and June 1997. Five areas with similar plant composition were repeatedly surveyed (total area of site was 10,000 km2), which included the main vegetation associations and the border of agricultural fields. Plants with flowers were selected and bees were collected, killed and mounted the same day, and plants collected (see Vergara- Briceño, 1999). Bees and plant vouchers are deposited in the collection at Universidad de Las Americas. In Carrillo Puerto, Quintana Roo field work was carried out in areas near the Carrillo Puerto-Valladolid highway and the road from Carrillo Puerto to Vigia Chico, which is within the Sian ka’an Biosphere Reserve, an area of a half million hectares (e.g. CIQRO, 1983). Five collecting sites were selected along a 48 km transect, which were visited every two months between June 2008 and September 2009. Bees were collected in each site within a 1 km radius, then killed, mounted and identified. Bees and plant vouchers are deposited in ECOSUR. Two sites were sampled in the state of Veracruz: Centro de Investigaciones Costeras La Mancha (CICOLMA- La Mancha) and in the vicinity of Totutla. At La Mancha, bees were collected monthly between March 2007 and March 2008 (five days per visit) along established trails that included different environmental types (see Rico-Gray, 1993). Bees and plants were collected, identified using the collections of Instituto de Ecologia, A.C. (IEXA, XAL) and deposited in those collections. The ranch “Finca El Mirador” is located in the vicinity of Totutla. Bees were collected every month in October-November 2005 and January through September 2006 in three sites within the ranch: El Desengaño, Loma Araucarias and Loma Ventura. A part of this site has been converted to an organic coffee plantation. Bees and plants were collected and plants were identified in XAL herbarium of Instituto de Ecologia, A.C. and bees were identified by C. Vergara (UDLA-P), with the help of Terry Griswold (USDA Bee Laboratory, Logan UT, U.S.A.) and Michael Engel (Biodiversity Institute, University of Kansas, Lawrence, KS, U.S.A.). Plants and bees are deposited in C. Vergara’s collection. VA Ramírez-Flores et al. - Bee-plant Networks 58 Si te A lti tu de ( m ) C lim at e Te m pe ra tu re o C Pr ec ip ita tio n (m m ) V eg et at io n L a M an ch a 19 º 3 6’ N –9 6º 2 2’ W 20 W ar m s ub hu m id w ith d ry p er io d, w et in s um m er , a nd a c ol d fr on t s ea so n in w in te r 24 -2 6 15 00 D un e pi on ee rs , s an d du ne s cr ub , t ro pi ca l d ry an d de ci du ou s fo re st s, m an gr ov e an d fr es h w at er s w am ps C ar ri llo P ue rt o 19 o 42 ’- 1 9o 2 5’ N 87 o 52 ’ – 8 8o 0 6’ W 30 W ar m s ub hu m id w ith s um m er r ai n 25 .8 11 58 Su bd ec id uo us tr op ic al f or es t a nd s ec on da ry ve ge ta tio n Te hu ac an -C ui ca tla n 18 o 14 ’ – 1 8o 2 2’ N 97 o 26 ’ – 9 3o 3 4’ W 60 0- 30 00 W ar m s ub hu m id , s um m er r ai n, m ar ke d dr y pe ri od d ur in g su m m er r ai n 21 40 0 Sr ub s w ith s pi ne s, c ol um na r ca ct i a ss oc ia - tio ns , a nd Y uc ca s p. a ss oc ia tio ns To tu tla 19 º 1 2’ N –9 0º 5 3’ W 10 10 Te m pe ra te , h um id w ith s um m er r ai n, ha rd ly a ny th er m ic o sc ill at io n 20 18 25 R us tic c of fe e pl an ta tio n, m on ta ne f or es t T ab le 1 . C ha ra ct er is tic s pe r si te ( G ar cí a 19 81 ; I N E G I 19 88 ; V al ie nt e B an ue t 1 99 8; D ía z- C as te la zo 2 00 4; C as til lo -C am po s et a l. 20 06 ; D el oy a et a l. 20 07 ). A ls o, s ee th e m et ho ds s ec tio n (“ st ud y si te s an d fie ld w or k” ). Network analyses The manner in which networks were sampled, only allows us to build binary (qualitative) networks. We defined each sampling site as an adjacency matrix C, in which Cij represents the number of interactions between bee species I and plant species j (Bascompte et al., 2003). We then calculated the following metrics derived from the network theory to describe the bee-plant networks in each of the four sites (Tehuacan-Cuicatlan valley, Quintana Roo, La Mancha, and Totutla): links per species, connectance, network specialization, nestedness, modularity, and robustness to species extinction. We used these metrics mainly because they provide a way to describe the organization of our networks in a way that allows direct comparison with previous studies on plant-pollinator networks and among our four sites. We calculated the average degree for bees and plants using the arithmetic means of degrees (number of interactions in which each species was involved) of all species. We calculated the (C) connectance, defined as the fraction of registered interactions with respect to the total number of possible interactions: C= I (B x P) where, I is the total number of interactions observed, B is the number of bees species, and P is the number of plant species (Jordano, 1987). We calculated the degree of specialization (H2’) for each network (Blüthgen et al., 2006). The index ranges from zero (extreme generalization) to one (extreme specialization). The index is derived from the Shannon index (H’ is not dependent on sample size), rare species are not equally presented by the index and is robust to changes in sampling intensity and the number of interacting species in the network (Poole, 1974; Blüthgen et al., 2006). We used the package bipartite (Dormann et al., 2009) in “R” to perform network-specialization analysis. We used the NODF metric (Nestedness metric based on Overlap and Decreasing Fill) (Almeida-Neto et al., 2008) to estimate nestedness in each network. To calculate nestedness we used ANINHADO 3.0 (Guimarães & Guimarães, 2006). Under perfect nestedness, generalist species interact with each other, while the specialist species interact with generalist species (Bascompte et al., 2003). We tested the significance of nestedness for each network through 1,000 simulated networks generated by the Null Model II (CE) (Bascompte et al., 2003), which allowed us to determine if the degree of nestedness was higher than those expected when compared to random patterns of interactions. In this null model, the probability that a bee-plant interaction occurs is proportional to the mean number of interactions of both species (Bascompte et al., 2003). Due to differences of sampling, and species richness, in our study sites, we also calculated nestedness controlling potential differences in richness and heterogeneity of interactions among our four sites by using a z-score, which normalizes the deviation from a mean value (Almeida-Neto et al., 2008). We also tested whether there were groups of bee species that interacted more frequently with a particular group of plant species (Mello et al., 2011). We assessed modularity (M) within each network based on simulated annealing (SA) (Guimerà et al., 2004; Guimerà & Amaral 2005) using the Sociobiology 62(1): 56-64 (March, 2015) 59 software MODULAR (Marquitti et al., 2013). The index ranges from zero to one, and high values of M indicate that different bees and plants species are grouped by modules that are semi-independent, within the network (Olesen et al., 2007). The significance of M also was estimated with Null Model II, described above. If a network was significantly modular, we also tested within each module for a nested pattern. We calculated network robustness for each site, and to species extinction based on the cumulative random removal of species (sensu Burgos et al., 2007). We also used the robustness index (R) to estimate the area under the extinction curve (Burgos et al., 2007; Dáttilo, 2012). This index ranges from zero to one, and values of R close to one come from robust network curves, because their extinction curve decreases slowly (Burgos et al., 2007). When compared to other curves, a robustness index is not as sensitive to curve shape (Burgos et al., 2007). The generalist core of species per network (Gc), and peripheral species, were categorized using core-periphery components (Dáttilo et al., 2013). When Gc >1 the species are those with the greatest number of interactions when compared with other species in the same trophic level. They constitute the generalist species core. In contrast, when Gc <1 the species exhibit a low number of interactions related to other species in the same trophic level, and form the periphery of the network (Dáttilo et al., 2013). Results Our results show that network topology among sites varied widely (Fig 1). The lowest species richness was recorded for the Quintana Roo site (15 plant and 25 bee species), and the highest was Tehuacan-Cuicatlan (88 plant and 231 bee species). The networks of La Mancha and Totutla had intermediate species richness (La Mancha: 65 plant and 43 bee species; Totutla: 56 plant and 156 bee species, Table 2). Networks with lower species richness were more connected: the Quintana Roo site exhibited C= 0.26 and Tehuacan- Cuicatlan C= 0.03 (Table 2). Greater species richness was linked to a higher number of interactions: Tehuacan-Cuicatlan, 764 bee- plant interactions; Totutla, 615 bee-plant interactions; La Mancha 134 bee-plant interactions; and the Quintana Roo site 98 bee-plant interactions. However, the distribution of interactions was different among the four environments, apparently unrelated to network size. Totutla exhibited the highest mean number of interactions per species (3.18), followed by the Quintana Roo site (2.45) and Tehuacan-Cuicatlan (2.39). The network from La Mancha had the least number of interactions per species (1.24) (Table 2). The Totutla network had the lowest specialization value (H2’= 0.39); the other specialization values were: the Quintana Roo site (H2’= 0.40), La Mancha (H2’= 0.69), and Tehuacan- Cuicatlan (H2’= 0.58). Moreover, the four networks were significantly nested (P < 0.05) (Z-scoreNODF values: Quintana Roo site= 1.68, La Mancha= 5.74, Tehuacan-Cuicatlan = 15.28, and Totutla= 24.33). Additionaly, the Tehuacan-Cuicatlan network exhibited significant modularity (M= 0.469) (P= 0.02), with only four of the resulting eight modules significantly nested (Table 2). Table 2. Network descriptors for bee-plant networks studied in four sites of Mexico. Network descriptor Tehuacan- Cuicatlan valley La Mancha Carrillo Puerto Totutla Plant richness 88 65 15 56 Bee richness 231 43 25 156 Connectance 0.03 0.04 0.26 0.08 Network specialization (H2’) 0.58 0.69 0.40 0.39 Links per species 2.39 1.24 2.45 3.18 No. of interactions 764 134 98 616 Nestedness (NODF) 11.76* 14.62* 38.39* 48.38* Nestedness (Z-ScoreNODF) 15.28 5.74 1.68 24.33 Modularity 0.46* 0.61 0.31 0.34 Robustness to cumulative plant extinctions (R) 0.32 0.49 0.86 0.49 Robustness to cumulative bee extinctions (R) 0.49 0.25 0.65 0.26 No. of species in the generalist core (% of total) 15 (17.04%) 10 (15.38%) 4 (26.66%) 6 (10.71%) No. of species in the generalist core (% of total) 24 (10.38%) 5 (11.62%) 4 (12%) 16 (10.25%) * represents P-values < 0.05 Fig 1. Ecological network involving bees (right) and plants (left) sampled in four different sites of Mexico: (A) Totutla, Veracruz (B) La Mancha, Veracruz (C) Carrillo Puerto, Quintana Roo and (D) Tehuacan- Cuicatlan valley, Puebla. Each node or rectangle represents one plant or bee species, and lines represent bee–plant interactions. Rectangle size represents the number of interactions of a certain species, for instance, a larger rectangle is equal to more interactions by that species. The robustness to cumulative extinctions showed a different pattern at each site. The most robust network was the Quintana Roo site, where plants exhibited a robustness of R= 0.86 and bees R= 0.65. The bees at La Mancha were the most vulnerable (R= 0.25), while plants (R= 0.49) exhibited a similar value that for Tehuacan-Cuicatlan. The Totutla VA Ramírez-Flores et al. - Bee-plant Networks 60 Discussion We found that in the four areas bees exhibit a nested pattern of foraging, independent of variation in local and landscape environmental factors or species composition. Plant– pollinator networks have been evaluated in different types of habitats around the world, and in most cases they exhibit a highly nested pattern (Olesen & Jordano, 2002; Bascompte et al., 2003; Jordano et al., 2003). Despite differences in habitat structure among areas, we observed that our networks were significantly nested, suggesting that different populations of plants and bees interact in a similar way independent of variations in local and landscape environmental factors, and are highly stable as predicted by previous studies (Bascompte et al., 2003; Vázquez et al., 2009; Hagen et al., 2012). Biologically, the nested pattern decreases interspecific competition within mutualistic networks and increases species coexistence, generating a stable pattern of plant-animal interactions (Bastolla et al., 2009). Based on this perspective, we show that our nested networks have great stability to simulated species loss when compared with structurally different networks (e.g. highly compartmentalized) (Mello et al., 2011; Dáttilo, 2012; Clemente et al., 2013). However, when we evaluated network robustness to simulated species extinction, we observed that areas with greater specialization (e.g. Tehuacan-Cuicatlan valley and La Mancha) tend to be less robust. This probably occurred because environments with a high number of specialist interactions tend to have lower functional redundancy, and therefore, when a species is extinct in a trophic level only a few species from another trophic level can “dampen” the system (Mouchet et al., 2010; Joner et al., 2011). Due to the high specialization of bee-plant networks in the Tehuacan-Cuicatlan valley and La Mancha generated by different mechanisms, it is expected that these areas are functionally less redundant and resilient to disturbances (e.g. Sánchez-Galván et al., 2012). This was confirmed when species were removed from these networks, the extinction curves declined faster compared with other networks, because deletion of one species necessarily causes the deletion of other partner species. Our results demonstrate that environments with many specialized interactions are highly susceptible to loss of species and this pattern should occur in other systems of ecological networks with a high degree of specialization (wasp-plant and ant-myrmecophyte) (Dáttilo et al., 2012; Clemente et al., 2013). It has been noted that as species richness increases in the community, the number of registered interactions also increases; however, connectance decreases exponentially (Jordano, 1987; Blüthgen, 2010). Our results corroborate these studies, since the most species-rich networks (Tehuacan-Cuicatlan and La Mancha) were the least connected networks. Our results suggest that this occurred because bees and plants within these networks had fewer interactions per species, when compared with the other study sites. This result indicates that most interactions in less connected networks tend to be pairwise (Fonseca & Ganade, network exhibited very similar values to those of La Mancha (bees R= 0.26, plants R= 0.49). Finally, the plants in Tehuacan- Cuicatlan were most vulnerable to species extinctions (R= 0.32), while bees exhibited higher robustness to species loss (R=0.49) (Table 2). The Quintana Roo site network exhibited the most robust species core (26.66% of total plant species in the network), while only three bee species (12% the total, see Appendix) were included in the network core. Compared to the other three sites, the number of plant species in the core of the Totutla network was the smallest (10.71%) (n= 6 species, see Appendix), with a similar result for core bee species (10.25%) (n= 16 species) (Figure 2). The La Mancha plant species core was the second smallest (15.38%), however, more species were included (n= 10 species, see Appendix), the La Mancha bee core exhibited an intermediate value (11.62%) (n= 5 species, see Appendix). Core bees for Tehuacan-Cuicatlan were fewer (10.38%), although composed by 24 species (see Appendix). The plant core exhibited intermediate size (17.04%) and was composed by 15 species (see Appendix). Overlap is meaningless among core species (plants, bees): 1 plant species between sites 1 and 2, and 1 plant species between sites 2 and 3; 1 bee species between sites 1 and 4, and 2 bee species between sites 1 and 2. Fig 2. Robustness to cumulative species removal of plants and bees in (A-B) Totutla, Veracruz (C-D) La Mancha, Veracruz (E-F) Carrillo Puerto, Quintana Roo and (G-H) Tehuacan-Cuicatlan valley, Puebla. The images on the left (A-C-E-G) simulate the extinction of plant species, and the images on the right (B-D-F-H) simulate the extinction of bee species. The x-axis axis represents the proportion of primary extinctions of a trophic level, and the y-axis represents the proportion of species in the other trophic level still alive. Sociobiology 62(1): 56-64 (March, 2015) 61 1996), which could generate a high degree of specialization within such networks, as found in our networks. Why were the Tehuacan-Cuicatlan and La Mancha networks the most specialized? We believe that different mechanisms could be responsible for the high degree of specialization found in these areas. In La Mancha, sampling was conducted in locations with different environmental complexities, and therefore, so that certain groups of plants and bees could be habitat specialists (i.e. collected in few sites) and cause more specialized interaction profiles (McPeek, 1996; Warren et al., 2001). The high degree of specialization in Tehuacan-Cuicatlan possibly occurred due to the extreme endemism that occurs in the area. For instance, areas with large numbers of endemic species tend to have more specialists and rare interactions (Olesen et al., 2002; Hansen et al., 2002), which is due to the great geographical isolation over evolutionary time and may lead to specialized processes among species (Thompson, 2005). Interestingly, when data on floristic richness and endemism is compared with that of other natural areas in Mexico, it becomes evident that Tehuacan-Cuicatlan is highest for endemic species (Méndez- Larios et al., 2004). Therefore, based on the differences between Tehuacan-Cuicatlan and La Mancha, it is important to consider that structural characteristics of the vegetation and endemism may be important factors shaping the ecological interactions involving bees and plants. Several studies have demonstrated that pollination networks can be highly modular (Olesen et al., 2007; Ramos-Jiliberto et al., 2010; Hernández-Yáñez et al., 2013). These structures only refer to the presence of groups of species (plants, animals) that interact more intimately with other species in their group rather than with other species in the network (Olesen et al., 2007). However, these studies have evaluated networks considering different taxonomic groups together. Here we only considered bees, and we should expect a scarce presence of modules, because having only one type of organism a certain historical evolution is shared. In fact, the presence of a modular pattern depends on several factors, from evolutionary to spatial and temporal (Olesen et al., 2007; Mello et al., 2011). Our results only show a modular pattern for the plant-bee network collected in Tehuacan-Cuicatlan valley. It is possible that through evolutionary time, the high endemism found in this area contributed to the convergence of traits in phylogenetically related bee groups, which currently may visit a group of common plants. Finally, part of the Totutla forest has been converted to an organic coffee plantation, altering the interactions naturally found there (Steffan-Dewenter et al., 2001; Tylianakis et al., 2007). Our results suggest that the abundance of Bidens pilosa (Compositae, Asteraceae), a common secondary plant species, could be explained due to the type of management used. In summary, our results show that structural characteristics of the environments over space are extremely important factors that shape interactions involving plants and bees. Moreover, we found that areas with endemic species tend to have more specialized interactions, and that the nature of this characteristic can also influence the patterns of interactions among organisms and the susceptibility to extinction of species. Therefore, our results reinforce the role of environmental variation in species maintenance and biodiversity in hyper- diverse systems. Acknowledgments We wish to thank the following organizations that allowed us to either accomplish field work or support the analyses and writing of this manuscript: Universidad de Alicante (NLR), Universidad Veracruzana (VRG, WD), and the Smithsonian Institution (DWR, RVG). For Tehuacan-Cuicatlan, funding was provided by the Comision Nacional para el Conocimiento y Uso de la Biodiversidad (CONABIO, H278 Apoidea (Hymenoptera) del Valle de Zapotitlan de las Salinas, Puebla), and by Consejo Nacional de Ciencia y Tecnologia (CONACyT, 2635-PN, “Diversidad, fenologia y flora visitada por apoideos del Valle de Zapotitlan de las Salinas, (Puebla)”. The study in Totutla was supported by a grant from Mexico’s Environmental Ministry (SEMARNAT- CONACyT 2002-C01-0194) to CV. The Dean office of Instituto de Investigacion y Posgrado (INIP) at Universidad de las Americas- Puebla, provided funding and logistic support. We thank Wilberto Colli-Ucán and Margarito Tuz-Novelo for their help in field work in Carrillo Puerto region. For the Totula study we wish to thank Jorge A. Mueller Grohmann for allowing us to work in his property. For the Tehuacan-Cuicatlan study Rhonda and Terry Griswold and Robert W. Brooks participated in the collection of specimens; Mariana Cuautle helped with field work for the entire collecting season of 1996-1997. We thank Terry W. Griswold (Bee Biology and Systematics Laboratory, USDA/ARS- University of Utah, Logan, Utah), Wallace LaBerge (Illinois Natural History Survey), Michael S. Engel and Robert W. Brooks (Snow Entomological Division, Nat ural History Museum, University of Kansas, Lawrence, Kansas), and Virginia Meléndez (Universidad Autónoma de Yucatán) for bee identification. References Almeida-Neto, M., Guimarães-Jr, P.R., Guimarães, P., Loyola, R.D. & Urlich, W. (2008). A consistent metric for nestedness analysis in ecological systems: reconciling concept and measurement. Oikos, 117: 1227-1239. Bascompte, J., Jordano, P., Melián, C.J. & Olesen, J.M. (2003). The nested assembly of plant-animal mutualistic networks. Proceedings of the National Academy of Sciences USA, 100: 9383-87. Bascompte, J. & Jordano, P. (2008). Redes mutualistas de especies. Revista de Investigacion y Ciencia, 384:50-59. Bascompte, J. & Jordano, P. (2014). Mutualistic Networks. Monographs in Population Biology No. 53. Princeton University Press, Princeton. Bastolla, U., Fortuna. M.A., Pascual-García, A., Ferrera, VA Ramírez-Flores et al. - Bee-plant Networks 62 A., Luque, B. & Bascompte, J. (2009). The architecture of mutualistic networks minimizes competition and increases biodiversity. Nature, 458: 1018-1020. Blüthgen, N. (2010). Why network analysis is often disconnected from community ecology: A critique and an ecologist’s guide. Basic and Applied Ecology, 11:185-195. Blüthgen, N., Menzel, F. & Blüthgen, N. (2006). Measuring specialization in species interaction networks. BMC Ecology, 6: 9. Burgos, E., Ceva, H., Perazzo, R.P.J., Devoto, M., Medan, D., Zimmermann, M. & Delbue, A.M. (2007). Why nestedness in mutualistic networks? Journal of Theoretical Biology, 249: 307-313. Burkle, L.A. & Alarcón, R. (2011). The future of plant – pollinator diversity: understanding interaction networks across time, space and global change. American Journal of Botany, 98:1-11. Castillo-Campos, G., Bello, A.C.T. (2006). Entornos Veracruzanos: La Costa de La Mancha, Instituto de Ecología, A.C., Xalapa, México. CIQRO (Centro de Investigaciones de Quintana Roo) (1983). Sian Ka’an: estudios preliminares de una zona de Quintana Roo propuesta como reserva de la Biosfera. CIQRO-SEDUE, México. Clemente, M.A., Lange, D., Del-Claro, K., Prezoto, F., Campos, N.R. & Barbosa, N.C. (2012). Flower-visiting social wasps and plants Interaction: network pattern and environmental complexity. Psyche, 2012, id.478431. Clemente, M.A., Lange, D., Dáttilo, W., Del-Claro. K., Prezoto, F. (2013). Social wasp-flower visiting guild interactions in less structurally complex habitats are more susceptible to local extinction. Sociobiology, 60: 337-344. Dáttilo. W. (2012). Different tolerances of symbiotic and nonsymbiotic ant-plant networks to species extinctions. Network Biology, 2: 127-138. Dáttilo, W., Guimarães, P.R. & Izzo, T.J. (2013). Spatial structure of ant-plant mutualistic networks. Oikos, 122: 1643- 1648. Deguines, N., Jono, C., Baude, M., Henry, M., Julliard, R. & Fontaine, C. (2014). Large-scale trade-off between agricultural intensification and crop pollination services. Frontiers in Ecology and the Environment, 12: 212-217. Deloya, C., Parrra-Tabla, V. & Delfín-González, H. (2007). Fauna de coleópteros Scarabaeidae Laparosticti y Trogidae (Coleoptera: Scarabaeoidea) asociados al bosque mesófilo de montaña, cafetales bajo sombra y comunidades derivadas en el centro de Veracruz, México. Neotropical Entomology, 36: 5-21. Díaz-Castelazo, C., Rico-Gray, V., Oliveira, P.S. & Cuautle, M. (2004). Extrafloral nectary-mediated ant-plant interactions in the coastal vegetation of Veracruz, Mexico: richness, occurence, seasonality, and ant foraging patterns. Ecoscience, 11: 472-481. Donatti, C.I., Guimarães, P.R., Galetti, M., Pizo, M.A., Marquitti, F.M.D. & Dirzo, R. (2011). Analysis of a hyper- diverse seed dispersal network: modularity and underlying mechanisms. Ecology Letters, 14: 773-781. Dormann, C.F., Fründ, J., Blütgen, N. & Gruber, B. (2009). Indices, graphs and null models: analyzing bipartite ecological networks. Open Ecology Journal, 2: 7-24. Dupont, Y.L., Padrón, B., Olesen, J.M. & Petanidou, T. (2009). Spatio-temporal variation in the structure of pollination networks. Oikos, 118:1261-1269. Dyer, L.A., Walla, T.R., Greeney, H.F., Stireman, J.O. & Hazen, R.F. (2010). Diversity of interactions: a metric for studies of biodiversity. Biotropica, 42: 281-289. Fonseca, C.R. & Ganade, G. (1996). Asymmetries, compartments and null interactions in an Amazonian ant- plant community. Journal of Animal Ecology, 65: 339-34. García, E. (1981). Modificaciones al Sistema Climático de Köppen, Universidad Nacional Autónoma de México (UNAM), México DF. Guimarães, P.R. & Guimarães, P. (2006). Improving the analyses of nestedness for large sets of matrices. Environmental. Modelling Software, 21: 1512-1513. Guimerà R, Sales-Pardo M, Amaral LAN (2004) Modularity from fluctuations in random graphs and complex networks. Physical Review E, 70:25-101. Guimerà, R. & Amaral, L.A.N. (2005). Functional cartography of complex metabolic networks. Nature, 433: 895-900. Hagen, M., Kissiling, W.D., Rasmussen, C., de Aguiar, M.A.M., Brown, L., Carstensen, D.W., Alves-dos-Santos, I., Dupont, Y.L., Edwards, F.K., Genini, J., Guimarães, P.R., Jenkins, G.B., Jordano, P., Kaiser-Bunbury, C.N., Ledger, M., Maia. K.P., Marquitti, F.M.D., McLaughlin, O., Morellato, L.P.C., O’Gorman, E.J., Trojelsgaard, K., Tylianakis, J.M., Vidal, M.M., Woodward, G. & Olesen, J.M. (2012). Biodiversity, species interactions and ecological networks in a fragmented world. Advances in Ecological Research, 46: 89-210. Hansen, D.M., Olesen, J.M. & Jones, C.G. (2002). Trees, birds and bees in Mauritius: exploitative competition between introduced honeybees and endemic nectarivorous birds. Journal of Biogeography, 29: 721-734. Hernández-Yáñez, H., Lara-Rodríguez, N., Díaz-Castelazo, C., Dáttilo, W. & Rico-Gray, V. (2013). Understanding the complex structure of a plant-floral visitor network from different perspectives in coastal Veracruz, Mexico. Sociobiology, 60: 329- 336. INEGI (Instituto Nacional de Estadística, Geografía e Informática) (1988) Síntesis geográfica, nomenclátor y anexo cartográfico del estado de Veracruz, México. INEG, México DF. Joner, F., Specht, G., Muller, S.C. & Pillar, V.D. (2011). Sociobiology 62(1): 56-64 (March, 2015) 63 Functional redundancy in a clipping experiment on grassland plant communities. Oikos, 120: 1420-142. Jordano, P. (1987). Patterns of mutualistic interactions in pollination and seed dispersal: connectance, dependence, and coevolution. American Naturalist, 129: 657-677. Jordano, P., Bascompte, J. & Olesen, J.M. (2003). Invariant properties in coevolutionary networks of plant-animal interactions. Ecology Letters, 6: 69-81. Jordano, P., Vázquez, D. & Bascompte, J. (2009). Redes complejas de interacciones mutualistas planta-animal. In: Medel R, Aizen M, Zamora R (eds), Ecología y Evolución de las Interacciones Planta-Animal: conceptos y aplicaciones. Editorial Universitaria, Santiago, Chile. Lange, D., Dáttilo, W. & Del-Claro, K. (2013). Influence of extrafloral nectary phenology on ant–plant mutualistic networks in a neotropical savanna. Ecological Entomology, 38: 463-469. Lewinsohn, T.M., Prado, P.I., Jordano, P., Bascompte, J. & Olesen, J.M. (2006). Structure in plant-animal interaction assemblages. Oikos, 113: 174-184. Loreau, M., Naeem, S., Inchausti, P., Bengtsson, J., Grime, J.P., Hector, A., Hooper, D.U., Huston, M.A., Raffaelli, D., Schmid, B., Tilman, D. & Wardle, D.A. (2001). Biodiversity and ecosystem functioning: current knowledge and future challenges. Science, 294: 804-808. Marquitti, F.M.D., Guimarães, P.R., Pires, M.M. & Bittencourt, L.F. (2013). MODULAR: Software for the autonomous computation of modularity in large network sets. Ecography, 37: 221-224. Méndez-Larios, I., Ortiz, E. & Villaseñor, J.L. (2004). Las Magnoliophyta endémicas de la porción xerofítica de la provincia florística del Valle de Tehuacán-Cuicatlán, México. Anales del Instituto de Biologia UNAM, 75: 87-104. Méndez-Larios, I., Villaseñor, J.L., Lira, R., Morrone, J.J., Dávila, P. & Ortiz, E. (2005). Toward the identification of a core zone in the Tehuacán–Cuicatlán biosphere reserve (México), based on parsimony analysis of endemicity of flowering plant species. Interciencia, 30: 267-274. Mello, M.A.R., Marquitti, F.M.D., Guimarães, P.R., Kalko, E.K.V., Jordano, P. & de Aguiar, M.A.M. (2011). The modularity of seed dispersal: differences in structure and robustness between bat–and bird–fruit networks. Oecologia, 167: 131-140. McPeek, M.A. (1996). Trade-offs, food-web structure, and the coexistence of habitat specialists and generalists. American Naturalist, 148:124-138. Michener, C.D. (2000). The Bees of the World. John Hopkins University Press, Baltimore. Mouchet, M.A., Villeger, S., Mason, N.W.H. & Mouillot, D. (2010). Functional diversity measures: an overview of their redundancy and their ability to discriminate community assembly rules. Functional Ecology, 24: 867-876. Olesen, J.M. & Jordano, P. (2002). Geographic patterns in plant- pollinator mutualistic networks. Ecology, 83: 2416-2424. Olesen, J.M., Eskildsen, I. & Venkatasamy, S. (2002). Invasion of pollination networks on oceanic islands: importance of invader complexes and endemic super generalists. Diversity and Distributions, 8: 181-192. Olesen, J.M., Bascompte, J., DuPont, Y.L. & Jordano, P. (2007). The modularity of pollination networks. Proceedings of the National Academy of Sciences USA, 104: 19891-19896. Olesen, J.M., Bascompte, J., Elberling, H. & Jordano, P. (2008). Temporal dynamics in a pollination network. Ecology, 89: 1573- 1582. Poole, R.W. (1974). An Introduction to Quantitative Ecology. McGraw.Hill, Ltd., Kogakusha. Ramos-Jiliberto, R., Domínguez, D., Espinoza, C., López, G., Valdovinos, F.S., Bustamante, R.O. & Medel, R. (2010). Topological change of Andean plant-pollinator networks along an altitudinal gradient. Ecological Complexity, 7: 86-90. Rico-Gray, V. (1993). Use of plant-derived food resources by ants in the dry tropical lowlands of coastal Veracruz, Mexico. Biotropica, 25: 301-315. Roubik, D.W. (1992). Ecology and Natural History of Tropical Bees. Cambridge University Press, Cambridge Rzedowski, J. (1978). Vegetación de México. Editorial Limusa, México. Sánchez-Galván, I.R., Díaz-Castelazo, C. & Rico-Gray, V. (2012). Effect of hurricane Karl on a plant-ant network occurring in coastal Veracruz, Mexico. Journal of Tropical Ecology, 28: 603-609. Santos, G.M.M., Aguiar, C.M.L., Genini, J., Martins, C.F., Zanella, F.C.V. & Mello, M.A.R. (2012). Invasive Africanized honeybees change the structure of native pollination networks in Brazil. Biological Invasions, 14: 2369-2378. Steffan-Dewenter, I., Münzenberg, U. & Tscharntke, T. (2001). Pollination, seed set and seed predation on a landscape scale. Proceedings of the Royal Society B, 268: 1685-1690. Stephens, C.R., Giménez, G.J., González, C., Ibarra-Cerdeña, C.N., Sánchez-Cordero, V. & González-Salazar, C. (2009). Using biotic interaction networks for prediction in biodiversity and emerging diseases. PLoS One 4: e5725. Thompson, J.N. (2005). The geographic mosaic of coevolution. University of Chicago Press, Chicago. Tylianakis, J.M., Tscharntke, T. & Lewis, O.T. (2007). Habitat modification alters the structure of tropical host-parasitoid food webs. Nature, 445: 202-205. Valiente-Banuet, A., Flores-Hernández, N., Verdú, M. & Dávila, P. (1998). The chaparral vegetation in Mexico under a non-mediterranean climate: the convergence and madrean- VA Ramírez-Flores et al. - Bee-plant Networks 64 tethyan hypotheses reconsidered. American Journal of Botany, 85: 1398-1408. Vázquez, D. & Aizen, M. (2004). Asymmetric specialization: a pervasive feature of plant-pollinator interactions. Ecology, 85: 1251- 1257. Vázquez, D., Blüthgen, N., Cagnolo, L. &, Chacoff, N.P. (2009) Uniting pattern and process in plant–animal mutualistic networks: a review. Annals of Botany, 103: 1445-1457. Vergara, C.H. & Ayala, R. (2002). Diversity, phenology and biogeography of the bees (Hymenoptera, Apoidea) of Zapotitlán de las Salinas, Puebla. Journal of the Kansas Entomological Society, 75: 16-30 Warren, M.S., Hill, J.K., Thomas, J.A., Asher, J., Fox, R., Huntley, B., Roy, D.B., Telfer, M.G., Jeffcoate, S., Harding, P., Jeffcoate, G., Willis, S.G., Greatorex-Davies, J.N., Moss, D. & Thomas, C.D. (2001). Rapid responses of British butterflies to opposing forces of climate and habitat change. Nature, 414:65-69. Zavala, J.A. (1982). Estudios ecológicos en el valle semiárido de Zapotitlán, Puebla y clasificación numérica de la vegetación, basada en atributos binarios de presencia o ausencia de las especies. Biótica, 7: 99-120. Appendix Species (plants, bees) in the core of the networks per sampled site La Mancha, Veracruz Plants Bees Bauhinia divaricata Apis melifera Bidens pilosa Lasioglossum sp1 Randia laetevirens Trigona nigra Turnera diffusa Ceratina sp1 Lantana camara Euglossa viridissima Waltheria indica Palafoxia lindenii Hyptis suaveolens Tecoma stans Piscidia piscipula Carrillo Puerto, Quintana Roo Plants Bees Gymnopodium floribundum Centris sp. Cassia fistula Augochlora sp. Croton campechianus Trigona fulviventris Parthenium hysterophorus Totutla, Veracruz Plants Bees Bidens pilosa Augochloropsis sp1 Elephantopus molis Ceratina sp2 Smallanthus maculatus Bombus medius Sida rhombifolia Lasioglossum sp1 Marsypianthes chamaedrys Paratetrapedia apicalis Cyperaceae sp1 Apis melifera Ceratina sp1 Trigona fulviventris Lasioglossum sp8 Augochlora aurifera Ceratina nautlana Melissodes tepaneca Lasioglossum costale Trigona corvina Trigona nigerrima Tehuacan-Cuicatlan valley, Puebla Plants Bees Salvia polystachya Plebeia frontalis Ipomoea pauciflora Agapostemon cockerelli Sanvitalia fruticosa Melissodes gilensis crenata Echinocactus platyacanthus Augochlorella neglectula Hechtia podantha Lasioglossum sp16 Gymnosperma glutinosum Centris atripes Tecoma stans Lasioglossum sp6 Wigandia urens Halictus ligatus Flaveria pringlei Paraugochloropsis vesta Prosopis laevigata Ashmeadiella zapotitlana Viguiera dentata Macrotera azteca Verbesina neotenoriensis Megachile zapoteca Beaucarnea gracilis Paraugochloropsis metallica Macroptilium atropurpureum Anthidium maculifrons Nissolia microptera Anthophora californica Exomalopsis pueblana Agapostemon sp1 Lasioglossum sp13 Centris transversa Ceratina sp3 Ceratina nautlana Lasioglossum sp18 Megachile otomita Melissodes tepaneca