Open access journal: http://periodicos.uefs.br/ojs/index.php/sociobiology ISSN: 0361-6525 DOI: 10.13102/sociobiology.v68i1.5451Sociobiology 68(1): e-5451 (March, 2021) Introduction Honeybees visiting flowers is such a common picture that it is hard to remember it is just a part of the pollination story. First, not all pollinators are honeybees. Pollination stands as one of the most important classes of animal–plant interaction, as ~85% of angiosperms are pollinated by animals (Ollerton et al., 2011). Such interactions were key drivers of the evolutionary radiation of both angiosperms and insects across the world (Crepet et al., 1991; Grimaldi, 1999). Within Hymenoptera, one of the largest insect orders, honeybees and other bee species are frequently cited as major pollinators (Cardinal et al., 2012; Kevan & Baker, 1983; Potts et al., 2010), but other insects such as wasps are also important pollinators (Heithaus, 1979; Freitas & Sazima, 2006; Sühs et al., 2009; Abstract Basic research assessing environmental effects on entire pollinator communities are still uncommon, particularly for rare and commercially unattractive plant- pollinator partners. We investigated the community of flower visitors of Erythroxylum myrsinites to check for potential pollinators and to check the extent of weather influence of visitor behaviour, registered as the number of visitors attending flowers of E. myrsinites. We then calculated species’ dominance and constancy and assessed location of pollen attachment in each visitor’s body. We correlated weather variables with the composition and abundance of visits carried out by the entire community and by most constant and dominant species. The wasps Polybia sericea, P. ignobilis and P. fastidiosuscula showed the highest values of constancy, dominance and attached pollen. There was a community-level effect of atmospheric pressure, solar radiation and wind speed on the number of visits. Atmospheric pressure affected the number of visits of eudominant species P. fastidiosuscula and P. sericea, while solar radiation affected the number of visits of P. ignobilis. Our results demonstrate the influence of weather variables on flower visiting insects and suggest the importance of native wasps in pollen transport and potential pollination. Sociobiology An international journal on social insects Rafael Barbizan Sühs1, Alexandre Somavilla2, Eduardo Luís Hettwer Giehl1 Article History Edited by Solange Augusto, UFU, Brazil Received 26 May 2020 Initial acceptance 15 December 2020 Final acceptance 29 December 2020 Publication date 22 February 2021 Keywords Community ecology, insect-community, insect-plant interaction, Polybia, wasps. Corresponding author Rafael Barbizan Sühs https://orcid.org/0000-0001-5921-7181 Laboratory of Human Ecology and Ethnobotany, Universidade Federal de Santa Catarina, Campus Florianópolis Av. Prof. Henrique da Silva Fontes, 2754 Trindade, Santa Catarina, Brasil. E-Mail: rbsuhs@gmail.com Somavilla et al., 2010; Rader et al., 2016). Second, not all insects visiting flowers are pollinators. Among other resources, flower visitors seek for pollen or nectar (Faegri & van der Pijl, 1971; Proctor et al., 1996) and are only pollinators if visits are frequent and pollen is regularly transferred to another flower (stigma) of the same species (Rech et al., 2014) or within the same individual. Last, the study of unusual plant-pollinator systems and factors driving pollination success may still uncover gaps in the “story of pollination”. Pollination success can be driven by factors such as flower and pollinator morphology and plant phenology or pollinator behaviour. Over short temporal and spatial scales, successful reproduction of plant pollinated species can result from the activity and identity of pollinators, whereas the behaviour and activity of pollinators can be constrained by 1 - Universidade Federal de Santa Catarina (UFSC), Florianopolis, Santa Catarina, Brazil 2 - National Institute for Amazonian Research (INPA), Manaus, Amazonas, Brazil RESEARCH ARTICLE - WASPS Weather variables affecting the behaviour of Insect flower visitors and main pollinators of Erythroxylum myrsinites Martius (Erythroxylaceae) RB Sühs, A Somavilla, ELH Giehl – Drivers of visits to E. myrsinites flowers2 abiotic factors (Case & Barrett, 2004). Air temperature and humidity, atmospheric pressure, solar radiation and wind intensity may influence insect visits, insect flight activity or even whether the insects leave their nests at all (Abrol, 1988; Burrill & Dietz, 1981; Fijen & Kleijn, 2017; Kearns & Inouye, 1993; Lundie, 1925; Nielsen et al., 2017). However, studies usually focus on the effect of weather variables either on few or a single species (e.g. honeybees: Burrill and Dietz, 1981; Lundie, 1925). This leaves a gap in the understanding of how weather affects the behaviour of entire pollinator communities (Fijen & Kleijn, 2017), particularly for rare and economically unattractive plant-pollinator partners. Erythroxylum myrsinites Mart. (Erythroxylaceae) is a shrub native to South America, occurring in Argentina, Brazil, Paraguay and Uruguay (Amaral Jr., 1980). In Brazil, E. myrsinites is found in Southern and South-Eastern regions, in the Atlantic Forest vegetation domain (Loiola & Costa- Lima, 2015) and is considered endangered in parts of its range (São Paulo 2016). Even though pollination of E. myrsinites is still underexplored, bees, wasps, and flies are known to be the main flower visitors and pollinators of several species of Erythroxylum (Barros, 1998; Freitas & Sazima, 2006; Rosas & Domínguez, 2009). Here we gathered information on the flower visitors of E. myrsinites and assessed what drives the activity and number of visitors on its flowers. By assuming pollination systems are normally evolutionary conserved (Macior 1982; Johnson et al. 1998; but see Ollerton 1996) and based on previous studies for Erythroxylum spp. (Barros, 1998; da Silva et al., 2007; Freitas & Sazima, 2006), we hypothesized that bees, wasps, and flies are the main flower visitors of E. myrsinites. Regarding the drivers of visitors’ activity, we expected air temperature and relative humidity, atmospheric pressure, solar radiation and wind intensity (e.g. Nielsen et al. 2017) to affect the community of visitors by influencing the number of visiting insects. By checking these hypotheses, we tried both to improve basic knowledge of a poorly known plant species and its interacting insect community and to apply community- level modelling techniques to assess environmental drivers of the behaviour of flower visitors. Materials and methods Study system The study was conducted in municipality of Santa Cruz do Sul (29º41’S – 52º26’W), Southern Brazil. Vegetation type is seasonal deciduous forest, being inside the Atlantic Forest domains. Following Koeppen’s classification (Köppen, 1931), the regional climate is Cfa (subtropical humid). Erythroxylum myrsinites individuals were found growing in the edge of an isolated forest remnant of ca. 9 ha at 45 m above sea level. Study plant Erythroxylum myrsinites is a light-demanding shrub and grows on moist soils, reaching up to 3 m in height (Amaral Jr., 1980). Flowers are solitary, small, light cream coloured and slightly perfumed (Amaral Jr., 1980). It is a distylous plant with brevistylous and longistylous flowers, very similar in size, produced from August to April (Amaral Jr., 1980). This kind of floral dimorphism, in which anthers and stigmas of the two floral morphs are in different positions, may promote pollen flux between the two flower types and increase cross-pollination (Barros, 1998). The fruit is a small drupe (ca. 7 mm length and 3.5 mm wide), red, monospermic and probably bird-dispersed. Fruiting spans from September to May (Amaral Jr., 1980). Sampling methods and data collection Flower visitors of seven individuals of E. myrsinites were collected during their complete flowering period between September and October 2008 (a different flowering period reported in literature). Visitors were collected with entomological nets, between 9h00 and 17h00. This temporal sampling extent was determined after preliminary observations where the period with the highest visitor activity was identified. Two additional observations were done at night to check if there were any nocturnal visitors, and no flower visitors were observed. Each sampling event lasted 30 min. Total time of observations and collections spanned over 32 h and was distributed along 16 days. Insects were deposited in individual vials and stored in the entomological collection of University of Santa Cruz do Sul (CESC). Subsequently, all collected individuals were examined under a stereoscopic microscope (at 40x) for identification and to verify if there was pollen attached to their integument and in which body part pollen grains were located (corbicula, face, frontal legs, mesosome, metasome, propleuron and scopa). From a fraction of these individuals, pollen was removed and visualised in optical microscope (at 1000x) to confirm if pollen grains belonged to E. myrsinites (same procedure than in Sühs et al. 2009). Weather data (air temperature and relative humidity, solar radiation, wind intensity and atmospheric pressure) were collected by a weather station (Davis®, model Vantage Pro Plus), located ca. 500 m away from the study site. Data was collected every 30 min and were provided by University of Santa Cruz do Sul (UNISC). Data analysis To verify the first hypothesis of this study, we recorded the identity, dominance, and constancy of flower visitors. As a complementary information regarding flower visitor behaviour, we also verified the number of individuals with pollen attached and the location of pollen on visiting species. Dominance (D) was calculated as: D = (i / t) × 100, where i = total number of individuals of i species; t = total number of collected individuals. Dominance was used to classify species in the following scheme: D > 10% = Eudominant; D > 5 ≤ 10% = Dominant; D > 2 ≤ 5% = Subdominant; D > 1 ≤ 2% = Recessive; D ≤ 1% = Rare (Sühs et al., 2009). Constancy (C) Sociobiology 68(1): e-5451 (March, 2021) 3 was calculated as: C = (P / N) × 100, where P = number of collections (number of records containing the species) and N = total number of records (Thomazini & Thomazini, 2002). When C ≤ 25% the species was considered Accidental, when C > 25 ≤ 50%, Accessory; and when C > 50%, Constant. We then listed species as potential pollinators when they showed high values of dominance and constancy, and when they showed E. myrsinites pollen on their bodies. We used generalized linear models for multivariate responses (GLMmv) with negative binomial errors to test the effects of weather predictors (air temperature, atmospheric pressure, relative humidity, solar radiation and wind speed) on community composition and on species individual responses. Following preliminary exploratory analyses, we decided to include also quadratic terms for all weather variables. Correlation amongst explanatory variables was previously checked using variance inflation factor (VIF) and redundant variables removed when VIF > 4, so that all kept variables were uncorrelated. All weather predictor variables were standardized before the analysis to avoid collinearity. P-values were calculated through 999 PIT-trap resamplings. Even though we applied a high number of tests on univariate responses (each species against predictors), we reported unadjusted P-values, but interpret the results with caution. To reduce the risk of type-I errors, we inform and discuss only responses of most abundant species (eudominant). Initial model contained all predictors, including quadratic terms. Model simplification was performed through backward elimination. Model selection was based on the sum of computed Akaike information criterion (AIC - by selecting models with lowest sum-of-AIC values) and validation was determined by graphical analysis of residuals. Analyses were performed on the R interface (R Development Core Team 3.0.1., 2013) using package mvabund (Wang et al., 2013) for model building. Results Erythroxylum myrsinites flowering period (since the opening of buds) lasted 19 days. The number of eudominant and subdominant visitors followed a bimodal pattern, peaking from 11h00 to 12h00 and from 13h00 to 14h00 and being low from 9h00 to 10h00, 12h00 to 13h00, and 16h00 to17h00 (Fig 1). We collected 376 visiting insects, 367 belonging to Hymenoptera and nine to Diptera. Nineteen species were recorded, 16 belonged to Hymenoptera and three to Diptera. Within Hymenoptera, ten species belonged to Vespidae, four to Halictidae and two to Apidae. Within Diptera, Syrphidae had two species and Stratiomyidae had one species (Table 1). Three visiting species were classified as eudominant and constant, all three belonging to the genus Polybia (Vespidae). These three species were: P. sericea (number of collected individuals, N = 182; constancy, C = 87.5%); P. ignobilis (N = 57; C = 57.8%); and, P. fastidiosuscula (N = 56; C = 56.3%). Four species were classified as subdominant and accessories: Apis mellifera (N = 16, C = 21.9%), Mischocyttarus rotundicollis (N = 14, C = 18.8%), Augochloropsis sp. (N = 13, C = 15.6%) and Brachygastra lecheguana (N = 9, C = 14.1%). There were no species classified in the dominant and accidental classes. The remaining 12 species were classified as either recessive or rare regarding dominance (Table 1). Fig 1. Number of individuals (NI) of eudominant and subdominant species collected per 0.5h periods over time of sampled days in Erythroxylum myrsinites’ flowers, Southern Brazil. RB Sühs, A Somavilla, ELH Giehl – Drivers of visits to E. myrsinites flowers4 Taxon N Do (%) Do Status Co (%) Co Status Pollen Location O.M. (%) C/S Fa Fl Ms Mt Pp Hymenoptera Vespidae Agelaia multipicta (Haliday, 1836) 3 0.8 Rr 4.7 Ac NA 0 0 0 0 0 0 Brachygastra lecheguana (Latreille, 1824) 9 2.4 Sd 14.1 Ac NA 2 2 0 0 0 0 Mischocyttarus rotundicollis (Cameron, 1912) 14 3.7 Sd 18.8 Ac NA 3 4 1 0 2 0 Pachodynerus guadulpensis (Saussure, 1853) 2 0.5 Rr 3.1 Ac NA 0 0 1 0 1 0 Polistes cavapyta Saussure, 1853 1 0.3 Rr 1.6 Ac NA 0 0 0 0 0 0 Polybia fastidiosuscula Saussure, 1854 56 14.9 Ed 56.3 Co NA 4 3 9 0 10 0 Polybia ignobilis (Haliday, 1836) 57 15.2 Ed 57.8 Co NA 23 32 35 0 22 26.3 Polybia platycephala Richards, 1951 1 0.3 Rr 1.6 Ac NA 0 0 0 0 0 0 Polybia scutellaris (White, 1841) 1 0.3 Rr 1.6 Ac NA 0 0 0 0 0 0 Polybia sericea (Olivier, 1792) 182 48.4 Ed 87.5 Co NA 66 144 153 0 110 62.1 Apidae Apis mellifera L. 1759 16 4.3 Sd 21.9 Ac 10 6 5 0 1 5 62.5 Plebeia sp. 6 1.6 Rc 9.4 Ac 2 0 0 0 1 0 33.3 Halictidae Augochlorella urania (Smith, 1853) 1 0.3 Rr 1.6 Ac 0 0 0 0 0 0 0 Augochloropsis sp. 13 3.2 Sd 15.6 Ac 2 1 3 3 1 1 33.3 Dialictus sp. 2 0.5 Rr 3.1 Ac 1 0 1 0 0 0 50 Neocorynura aenigma (Gribodo, 1894) 3 0.8 Rr 4.7 Ac 1 0 0 0 1 0 66.6 Diptera Stratiomyidae Chorisops sp. 7 1.9 Rc 7.8 Ac NA 0 0 0 0 0 0 Syrphidae Taxomerus sp.1 1 0.3 Rr 1.6 Ac NA 0 0 0 0 0 0 Taxomerus sp.2 1 0.3 Rr 1.6 Ac NA 0 0 0 0 0 0 Table 1. Recorded species in Erythroxylum myrsinites’ flowers, Southern Brazil. N = number of individuals, Do = Dominance, Do Status = Dominance Status (Ed = eudominant, Sd = subdominant, Rc = recessive and Rr = rare). Co = Constancy, Co Status = Constancy Status (Co = constant and Ac = accessory). Pollen Location (number of records): C/S (corbicula or scopa); Fa (face); Fl (frontal legs); Ms (mesosome/thorax); Mt (metasoma/abdomen); Pp (propleuron). O.M. (individuals with pollen analyzed in optical microscope). NA (structure not available in taxon). We observed that 70% of the visiting insects had pollen attached to their integument. We analysed 229 wasps (Vespidae) in detail (totalling 60.9% of all collected wasps) and the most common body structures showing attached pollen were mesosome (62% of inspected individuals), forelegs (56.7%), propleuron (44.5%) and face (29.4%) (Table 1). Of the collected wasps, 26.6% had no pollen attached. The highest number of records of pollen in body structures was recorded for Polybia sericea, followed by P. ignobilis (Fig 2). For bees (Apidae and Halictidae), 24 individuals (totalling 58.5% of all collected bees) had pollen attached to their bodies. We found pollen of E. myrsinites attached to all visiting insects inspected in optical microscope (N = 147). Out of those insects, 77% belonged to the wasp species P. sericea (Fig 3) and 10.1% to P. ignobilis (Vespidae), and 12.9% belonged to bees (Apidae and Halictidae). At the community level, the number of visiting insects correlated with atmospheric pressure (second- order parameter: LR = 42.5, p = 0.001), solar radiation (first-order: LR = 37.9, p = 0.001) and wind speed (first- order: LR = 35.5, p = 0.001; all results from the selected most parsimonious model described in detail in Table 2). Graphical analysis showed that residuals of the model were adequate. More information on model selection can be found in Supplementary file (Table S1). Sociobiology 68(1): e-5451 (March, 2021) 5 Fig 2. Number of records (NR) of pollen in body structures in individuals of eudominant and subdominant insect species in Erythroxylum myrsinites’ flowers, Southern Brazil. Fig 3. Images of Polybia sericea visiting Erythroxylum myrsinites flowers, Southern Brazil. The two images (A and B) show the visiting pattern of the wasp as well as pollen grains attached to its body. At the species level, weather variables had an overall effect on the eudominant species Polybia fastidiosuscula (LR = 8.21, p = 0.032) and Polybia sericea (LR = 8.09, p = 0.039), both affected by atmospheric pressure (second order). P. gnobolis, although affected by solar radiation (LR = 3.57, P = 0.041), was not affected by the set of weather variables (LR = 5.15, p = 0.125; Table 3). Wind speed did not affect significantly any eudominant visitor species. Graphical information about weather variables over the collection period is shown in Figure 4. Table 2. Selected most parsimonious Generalized Linear Model for multivariate response data (GLMmv) built for abundance of visitors to Erythroxylum myrsinites flowers, Southern Brazil. LR = Likelihood Ratio. Significant p-values are in bold. Variable LR p-value Model AIC Model test Statistics Model p-value Intercept 182.18 0.001 1141.2 101.5 0.0009 Atmospheric pressure² 42.54 0.001 Solar Radiation 37.9 0.001 Wind Speed 35.46 0.001 RB Sühs, A Somavilla, ELH Giehl – Drivers of visits to E. myrsinites flowers6 Model Parameters Polybia fastidiosuscula Polybia ignobilis Polybia sericea LR p-value LR p-value LR p-value Intercept 1.34 0.235 0 0.995 59.06 0.001 Atmospheric Pressure² 7.189 0.008 1.543 0.206 4.015 0.049 Solar Radiation 1.242 0.268 3.568 0.041 3.487 0.080 Wind Speed 2.747 0.097 1.331 0.233 0.858 0.368 Univariate Test Statistics 8.211 0.032 5.152 0.125 8.086 0.039 Table 3. Univariate estimates of deviance and associated P-values of modelled variables for eudominant species recorded on Erythroxylum myrsinites flowers, Southern Brazil. LR = Likelihood Ratio. Significant P-values are in bold. Fig 4. Variation of weather variables during the flowering period of Erythroxylum myrsinites, Southern Brazil. Lines represent daily averages and shaded areas the standard deviations. Discussion We found wasps are the most important flower visitors of E. myrsinites. The wasp Polybia sericea (Vespidae) was the most abundant species with high levels of constancy and pollen of the studied plant species attached to their integument, and thus can be considered the visitor with the highest potential to pollinate E. myrsinites. Other two species of wasps of the Sociobiology 68(1): e-5451 (March, 2021) 7 same genus (P. ignobilis and P. fastidiosuscula) also had high abundance and constancy. These eudominant species had their visiting activity associated with atmospheric pressure, solar radiation, and wind speed. The genus Erythroxylum has a variety of potential pollinators, especially flies, bees and wasps. We recorded all those groups visiting flowers of E. myrsinites, a pattern also reported by da Silva et al. (2007), for E. cf. macrophyllum. Bees and wasps were considered the main flower visitors of E. campestre St. Hil., E. suberosum St. Hil. and E. tortuosum Mart. (Barros, 1998). Wasps and flies were recognized as the main pollinators of E. microphyllum (Freitas & Sazima, 2006) and bees were considered the main pollinators for E. havanense (Rosas & Domínguez, 2009). In addition, our results suggest the wasp Polybia sericea potentially pollinates E. myrsinites, an additional instance to the intrinsic relationship between other species of Erythroxylum and wasps (Barros, 1998). Wind intensity, air temperature and solar radiation are important weather variables influencing insect visits and may directly affect the community of flower visitors (Abrol, 1988; Kearns & Inouye, 1993; Kjøhl et al., 2011; Nielsen et al., 2017). Our results showed that atmospheric pressure, solar radiation, and wind speed affected the number of visits of the entire community of insects found on flowers of E. myrsinites. Except for wind speed, the same set of variables (either together or alone) affected eudominant species and, thus, can be important to regulate success of E. myrsinites pollination. Moreover, models were of quadratic form, suggesting the existence of optimal values (peaks) in the number of visits. Air temperature did not affect the community of visitors, which may be due to its low variation along the studied period. Although weather conditions and abiotic stress are important factors in determining the success of insect visits (Case & Barrett, 2004; Kearns & Inouye, 1993), few studies directly evaluated their effects on the community level. Our results showed that 87% of collected individuals were wasps (Vespidae) while 11% were bees (Apidae and Halictidae). The same pattern was found for the species Schinus terebinthifolius Raddi, in the same forest fragment (Sühs et al., 2009). Visitation patterns found in both studies might be linked to local microclimatic conditions and visitor availability instead of a visitor preference for the flowers. However, the abundance and the seasonal pattern of common visitor species differ. In this study (during the spring), Polybia sericea was the most abundant species, with 55% of Vespidae individuals, followed by P. ignobilis and P. fastidiosuscula, both with 17% of collected Vespidae individuals. Conversely, these species individually accounted for less than 15% of collected Vespidae individuals and had great part of individuals collected in winter (Sühs et al., 2009). Although these studies were carried out in different years, it seems that there is a preference of distinct wasp species for flowers of either E. myrsinites or S. terebinthifolius. There was a high production of fruits in all sampled individuals after the study, indicating that pollination event occurred. Although not tested, based on our results, we can infer that pollination probably was mainly crossed and done by insects because floral dimorphism in Erythroxylum species – considered a mechanism of self-incompatibility (Pailler et al., 1998) – promotes the pollen flux between the two types of flowers, enhancing cross-pollination (Barros, 1998), and a previously reported rejection of pollen tubes resulting from self-pollination and interbreeding for several species of Erythroxylum (Barros, 1998; Ganders, 1979). Nevertheless, we encourage future studies in evaluating the effectivity of different groups of flower visitors as pollinators, which besides depending on number and constancy of visits, also depends on the quantity of pollen grain still available for the pollination of the flowers transported on the body (e.g. Garibaldi et al., 2013). Although carried out over a short time, this study provides novel data on flower visitors and possible pollinators of a little known and threatened plant species – at least in part of its range. Wasps were the main flower visitors, pollen carriers and potential pollinators of Erythroxylum myrsinites. The wasp Polybia sericea was considered the main potential pollinator of the species because of high values of dominance and constancy, and a high number of individuals with pollen attached to several of its body parts. On the community-level and for eudominant species, we found weather variables such as air humidity, atmospheric pressure, wind speed and solar radiation to affect the community of visitors. This study reinforces the role of wasps in pollen transport and potential pollination. Acknowledgements We thank David Warton for the suggestions with statistical analyses and the anonymous referees for their contributions. The authors declare no conflict of interest. Author’s contributions RB Sühs, ELH Giehl, A Somavilla – Conceptualization RB Sühs, ELH Giehl, A Somavilla – Methodology RB Sühs, A Somavilla – Formal analysis RB Sühs, A Somavilla – Investigation RB Sühs, A Somavilla – Data Curation RB Sühs, A Somavilla – Visualization RB Sühs, ELH Giehl, A Somavilla – Writing and revising Supplementary Material DOI: 10.13102/sociobiology.v68i1.5451.s3401 http://periodicos.uefs.br/index.php/sociobiology/rt/suppFiles/5451/0 References Abrol, D.P. (1988). Environmental factors influencing pollination activity of Apis mellifera on Brassica campestris. Journal of the Indian Institute of Science, 68: 49-52. RB Sühs, A Somavilla, ELH Giehl – Drivers of visits to E. myrsinites flowers8 Amaral Jr., A. (1980). Erythroxylaceae. In R. Reitz (Ed.), Flora Ilustrada Catarinense. (pp. 1-64). Barros, M.G. (1998). Sistemas reprodutivos e polinização em espécies simpátricas de Erythroxylum P. Br.(Erythroxylaceae) do Brasil. Revista Brasileira de Botânica, 21: 1-11. doi: 10.1590/S0100-84041998000200008 Burrill, R.M. & Dietz, A. (1981). The response of honey bees to variations in solar radiation and temperature. Apidologie, 12: 319-328. Cardinal, S. & Danforth, B.N. (2012). Bees diversified in the age of eudicots. Proceedings of the Royal Society B, 280. doi: 10.1098/rspb.2012.2686 Case, A.L. & Barrett, S.C.H. (2004). Environmental stress and the evolution of dioecy: Wurmbea dioica (Colchicaceae) in Western Australia. Evolutionary Ecology, 18: 145-164. Crepet, W.L., Friis, E.M., Nixon, K.C., Lack, A.J. & Jarzembowski, E.A. (1991). Fossil evidence for the evolution of biotic pollination. Philosophical Transactions of the Royal Society of London B: Biological Sciences, 333(1267): 187-195. da Silva, F.J.T., Schwade, M.R.M. & Webber, A.C. (2007). Fenologia , biologia floral e polinização de Erythroxylum cf macrophyllum (Erythroxylaceae), na Amazônia Central. Revista Brasileira de Biologia, 5: 186-188. Faegri, K. & van der Pijl, L. (1971). Principles of Pollination Ecology (2nd ed.). Elsevier. doi: 10.1016/B978-0-08-023160- 0.50018-9 Fijen, T.P. M. & Kleijn, D. (2017). How to efficiently obtain accurate estimates of flower visitation rates by pollinators. Basic and Applied Ecology, 19: 11-18. doi: 10.1016/j.baae. 2017.01.004 Freitas, L. & Sazima, M. (2006). Pollination biology in a tropical high-altitude grassland in Brazil: interactions at the community level. Annals of the Missouri Botanical Garden, 93: 465-516. Ganders, F.R. (1979). Heterostyly in Erythroxylum coca (Erythroxylaceae). Botanical Journal of the Linnean Society, 78: 11-20. doi: 10.1111/j.1095-8339.1979.tb02182.x Garibaldi, L.A., Steffan-Dewenter, I., Winfree, R., Aizen, M.A., Bommarco, R., Cunningham, S.A., Kremen, C., Carvalheiro, L.G., Harder, L.D., Afik, O., Bartomeus, I., Benjamin, F., Boreux, V., Cariveau, D., Chacoff, N.P., Dudenhöffer, J.H., Freitas, B.M., Ghazoul, J., Greenleaf, S., … Klein, A.M. (2013). Wild pollinators enhance fruit set of crops regardless of honey bee abundance. Science, 340(6127): 1608-1611. doi: 10.1126/science.1230200 Grimaldi, D. (1999). The co-radiations of pollinating insects and angiosperms in the Cretaceous. Annals of the Missouri Botanical Garden: 86: 373-406. Heithaus, E.R. (1979). Community structure of neotropical flower visiting bees and wasps: diversity and phenology. Ecology, 60: 190-202. Johnson, S., Linder, H. & Steiner, K. (1998). Phylogeny and radiation of pollination systems in Disa (Orchidaceae). American Journal of Botany, 85: 402-411. Kearns, C.A. & Inouye, D.W. (1993). Techniques for pollination biologists. University press of Colorado. University Press of Colorado. Kevan, P.G. & Baker, H.G. (1983). Insects as Flower Visitors and Pollinators. Annual Review of Entomology, 28: 407-453. doi: 10.1146/annurev.en.28.010183.002203 Kjøhl, M., Nielsen, A., Stenseth, N.C. & others. (2011). Potential effects of climate change on crop pollination. Food and Agriculture Organization of the United Nations (FAO). Köppen, W.P. (1931). Grundis der Klimakunde. In Berlin, Germany. Walter de Gruyter & Co. Loiola, M.I.B., & Costa-Lima, J.L. (2015). Erythroxylaceae in Lista de Espécies da Flora do Brasil. Jardim Botânico do Rio de Janeiro. retrieved from: http://floradobrasil.jbrj.gov.br/ jabot/floradobrasil/FB111 Lundie, A.E. (1925). The flight activities of the honeybee. United States Department of Agriculture, 1328: 1-38. doi: 10.1007/s13398-014-0173-7.2 Macior, L.W. (1982). Plant community and pollinator dynamics in the evolution of pollination mechanisms in Pedicularis (Scrophulariaceae). In Armstrong, J.A.,Powel, J.M. and Richards, A.J. (eds), Pollination and Evolution, Sydney, Royal Botanical Gardens, 29-45. Nielsen, A., Reitan, T., Rinvoll, A.W. & Brysting, A.K. (2017). Effects of competition and climate on a crop pollinator community. Agriculture, Ecosystems and Environment, 246: 253-260. doi: 10.1016/j.agee.2017.06.006 Ollerton, J. (1996). Reconciling Ecological Processes with Phylogenetic Patterns: The Apparent Paradox of Plant-- Pollinator Systems. Journal of Ecology, 84: 767-769. Ollerton, J., Winfree, R. & Tarrant, S. (2011). How many flowering plants are pollinated by animals? Oikos, 120: 321-326. Pailler, T., Humeau, L. & Thompson, J. D. (1998). Distyly and heteromorphic incompatibility in oceanic island species of Erythroxylum (Erythroxylaceae). Plant Systematics and Evolution, 213: 187-198. doi: 10.1007/BF00985199 Secretaria de Estado do Meio Ambiente. DOE 07-06-2016, Seção I, p. 69-71, Resolução SMA No 57, (2016). Potts, S.G., Biesmeijer, J.C., Kremen, C., Neumann, P., Schweiger, O. & Kunin, W.E. (2010). Global pollinator declines : trends , impacts and drivers. Trends in Ecology and Evolution, 25: 345-353. doi: 10.1016/j.tree.2010.01.007 Proctor, M., Yeo, P. & Lack, A. (1996). The natural history of pollination. HarperCollins Publishers. R Development Core Team 3.0.1. (2013). A Language and Sociobiology 68(1): e-5451 (March, 2021) 9 Environment for Statistical Computing. In R Foundation for Statistical Computing. http://www.r-project.org Rader, R., Bartomeus, I., Garibaldi, L.A., Garratt, M.P.D., Howlett, B.G., Winfree, R., Cunningham, S.A., Mayfield, M.M., Arthur, A.D., Andersson, G.K.S. & others. (2016). Non-bee insects are important contributors to global crop pollination. Proceedings of the National Academy of Sciences, 113: 146-151. Rech, A.R., de Avila Jr, R.S. & Schlindwein, C. (2014). Síndromes de polinização: especialização e generalização. In Biologia da Polinização (pp. 171–181). Rosas, F. & Domínguez, C.A. (2009). Male sterility, fitness gain curves and the evolution of gender specialization from distyly in Erythroxylum havanense. Journal of Evolutionary Biology, 22: 50-59. doi: 10.1111/j.1420-9101.2008.01618.x Somavilla, A., Sühs, R. B. & Köhler, A. (2010). Entomofauna associated to the floration of Schinus terebinthifolius Raddi (Anacardiaceae) in the Rio Grande do Sul state, Brazil. Bioscience Journal, 26: 956-965. Sühs, R.B., Somavilla, A., Köhler, A. & Putzke, J. (2009). Pollen vector wasps (Hymenoptera, Vespidae) of Schinus terebinthifolius Raddi (Anacardiaceae), Santa Cruz do Sul, RS, Brazil. Brazilian Journal of Biosciences, 7: 138-143. Thomazini, M.J. & Thomazini, A.P. (2002). Diversidade de abelhas (Hymenoptera: Apoidea) em inflorescências de Piper hispidinervum (C. DC.). Neotropical Entomology, 31: 27-34. Wang, Y., Naumann, U., Wright, S.T., Warton, D.I., Wang, A.Y. & Naumann, U. (2013). Package “mvabund”: Statistical methods for analysing multivariate abundance data. Methods in Ecology and Evolution, 3: 471-474.