Bees for Development: Brazilian Survey Reveals How to Optimize Stingless Beekeeping

Stingless bees are an important asset to assure plant biodiversity in many natural ecosystems, and fulfill the growing agricultural demand for pollination. However, across developing countries stingless beekeeping remains an essentially informal activity, technical knowledge is scarce, and management practices lack standardization. Here we profited from the large diversity of stingless beekeepers found in Brazil to assess the impact of particular management practices on productivity and economic revenues from the commercialization of stingless bee products. Our study represents the first large-scale effort aiming at optimizing stingless beekeeping for honey/colony production based on quantitative data. Survey data from 251 beekeepers scattered across 20 Brazilian States revealed the influence of specific management practices and other confounding factors over productivity and income indicators. Specifically, our results highlight the importance of teaching beekeepers how to inspect and feed their colonies, how to multiply them and keep track of genetic lineages, how to harvest and preserve the honey, how to use vinegar traps to control infestation by parasitic flies, and how to add value by labeling honey containers. Furthermore, beekeeping experience and the network of known beekeepers were found to be key factors influencing productivity and income. Our work provides clear guidelines to optimize stingless beekeeping and help transform the activity into a powerful tool for sustainable development.


Survey data
A draft survey questionnaire was designed in 2012 in collaboration with several researchers and successful beekeepers from Brazil. The survey included questions on the species kept, the number of colonies, management techniques, and sales of bee products, among others, as well as an informed consent to participate in this study that was signed by all participants. As our study did not involve indigenous people, and does not disclose any personal information from its participants, we did not seek approval by an ethics committee. This initial questionnaire was tested across Rio Grande do Norte, a small Brazilian State with a rich tradition in meliponiculture [35]. The analysis of the data retrieved from this initial examination allowed us to refine the survey, eliminating meaningless questions, adding new ones, and rewording to obtain more accurate answers (see S1 Questionnaire for final questionnaire). We used this questionnaire to personally interview several beekeepers from different Brazilian States. A pdf version of the final questionnaire was also distributed to research institutions and extension officers throughout Brazil, whom interviewed local beekeepers. In addition, an online version of the survey (created using Google Forms) was placed in the University of São Paulo´s Bee Lab website, and disseminated through the main Brazilian social networks related to meliponiculture (Yahoo and Facebook groups). The project was presented in local meliponiculture conferences and the survey promoted via an article published in the main beekeeping journal of Brazil [36]. Because many Brazilian stingless beekeepers live in remote rural areas lacking internet, telephone, and access to postal services, we were not able to use these communication channels as effectively as previous surveys undertaken in Australia [17,37].
Survey data was collected during one year, between January 2013 and January 2014. Only data from beekeepers having at least one stingless bee colony at the moment of the interview were included in this study. Although respondents might have provided false or inaccurate answers in some cases, we have no reason to expect systematic response biases that could compromise our analyses. To minimize such potential biases, we carefully curated the data eliminating duplicated entries, contacting some beekeepers again to confirm certain answers, cross validating replies, and checking for outliers in each response. Four outlier observations were deleted because they were unrealistic, based on the other information provided and our current knowledge in meliponiculture (see S1 Dataset Rscripts for final dataset).

Statistical analyses
A model selection approach was used to relate production and income indicators to management practices and other confounding variables. Such an approach has gained substantial support in the natural sciences during the last decade, and is particularly well suited to analyze complex datasets, when several different competing hypotheses can be put forward [38,39]. Because we expected production and income to be mainly determined by the ability of beekeepers to multiply colonies, sell colonies, and sell honey, we first constructed three logistic models to identify which predictors explained whether beekeepers multiply colonies, sell honey, or sell colonies. We then constructed nine additional models relating the same predictors to nine different continuous indicators of productivity and income. We considered thirtynine different predictors (categorical and continuous) which described the beekeepers' social background, property characteristics, beekeeping experience, management practices, and other confounding factors (see S1 Table for details). For some models we also included the total number of colonies kept as an additional predictor, to account for variation in the size of beekeeping operations. In each model, we first constructed model subsets containing groups of related predictors (education, property characteristics, or management practices), and included the main species as a random effect (assuming different species vary in optimal management practices). We initially considered groups of predictors separately, because the data contained many missing observations for particular combinations of predictors, and thus we could not construct one single full model. For each subset of predictors, a full model containing all or most predictors was compared to reduced models where each predictor was eliminated one by one. A predictor was excluded if the reduced model (lacking that particular predictor) was not significantly different from the full model (using likelihood ratio tests, alpha = 0.05). Predictors were eliminated sequentially, based on the magnitude of p-values, and each time a predictor was eliminated we constructed a new full model with the larger dataset (including the observations that were missing in the dropped predictor) and repeated the whole model selection process. All statistical comparisons between two models were conducted using the same dataset. To avoid overfitting the data, care was taken not to construct models containing less than 15 observations per predictor. Once we identified the key predictor variables in each subset, we constructed a complete-cases dataset with these predictors and ran a similar model selection protocol to obtain a final model. Interaction effects were included in these final models to test whether they improved the models without interactions. Finally, we used the Akaike Information Criterion (AIC) to compare the best models with and without random intercepts, and in some cases with a random slope (where we had a priori justification for considering a differential response to predictors across bee species). The best models were validated by plotting residual vs. fitted values, residual vs. predictors, by looking at the distribution of residuals, by checking for multicollinearity, and by comparing observed data to simulations from the models. All analyses and graphs were implemented with R, version 3.1.0 [40] (R-scripts for all models are available in S1 Dataset Rscripts).

Results
The final data set consisted of 251 observations from stingless beekeepers distributed across 20 Brazilian States (S2 Table). Men dominated the data set (only 14 of the interviewed beekeepers were woman), and nearly half of all respondents had jobs within the agriculture sector, were public servants or pensioners. The age of beekeepers ranged between 15 and 80, while the time keeping bees varied between less than a year and 54 years (see S3 Table for summary statistics). The total number of colonies kept by beekeepers ranged between 1 and 3500, with the majority keeping less than 100 colonies (Fig 1). While 66% of beekeepers multiplied at least one colony per year, only 25% sold colonies, 30% sold honey, and 14% sold both colonies and honey during the previous year. Although we registered a total of 19 species kept by the surveyed beekeepers, the most commonly kept species were Tetragonisca angustula, Melipona quadrifasciata, M. scutellaris, M. subnitida, and M. fasciculata (Fig 2, S4 Table).
Experience, technical skills related to bee management, and some property characteristics were found to be significant predictors of whether beekeepers multiply colonies, sell colonies, or sell honey ( Table 1). The best models for nine different continuous indicators of productivity and income are summarized in Table 2. Regression coefficients, p-values, and confidence intervals for all models are summarized in S5 Table, while S1 and S2 Figs show the relationships between response and predictor variables respectively. Years spent keeping bees and the number of known beekeepers were the most consistently included predictors in the best models. The number of colonies sold, the liters of honey sold, and the total number of hives kept were all positively associated with years spent keeping bees (Fig 3A). The number of known beekeepers was also positively associated with the total number of colonies kept and the number of multiplied colonies per year (Fig 3B). Beekeepers feeding their colonies with sugar syrup or honey had a larger number of colonies and multiplied more colonies per year (Fig 4A), whereas beekeepers employing established honey conservation methods sold more honey and had higher yearly earnings than those storing honey in a refrigerator or simply leaving it outside ( Fig 4B). Beekeepers who inspected their colonies more often and harvested honey with a  syringe lost fewer colonies per year (Fig 5A and 5B). The use of vinegar to control parasitic flies was associated with a larger number of colonies (Fig 5C), and selective breeding was found to  Beekeepers with more colonies of the main species, more years of experience keeping bees, and using an established honey conservation method, sell more honey per year. These effects hold across the main species kept, although the magnitude of the relationship between Years keeping bees and Liters of honey sold varies.
Earnings in R$ 3 75 LM Number of colonies, Labeling of honey containers, Honey conservation method and the interaction Number of colonies-Labeling of honey containers -Beekeepers with more colonies that label honey containers and use an established honey conservation method, have higher yearly earnings. The influence of the number of colonies on earnings is more pronounced among beekeepers that label honey containers. These effects do not vary across the main species kept.
Costs in R$ c 132 LM Number of colonies, Education level, Feeding frequency, Age, and the interaction Feeding frequency-Age -Yearly costs increase with the number of colonies, the level of education, the feeding frequency, and the age of beekeepers. Older beekeepers spend more money regardless how frequently they feed their colonies. These effects not vary across the main species kept.
The number of observations included in each model is provided (N) along with the model type, the model structure and its biological interpretation.
Regression coefficients, p-values, and confidence intervals for all models are summarized in S5 increase honey production per colony (Fig 5D). The total number of colonies kept was a key predictor of earnings. However, the relationship between earnings and the number of colonies was steeper among beekeepers labeling honey containers than among those that did not employ labels (Fig 6A). Similarly, we found that costs were explained by an interaction between the beekeeper's age and feeding frequency, whereby the relationship between age and costs was steeper among beekeepers that did not feed their colonies very often (Fig 6B).
The number of observations included in each model is provided (N) along with the model structure and its biological interpretation. All models are generalized linear mixed models (GLMM) with a Bernoulli distributed response variable (logistic regressions). Regression coefficients, p-values, and confidence intervals for all models are summarized in S5 Table. When asked to identify the main problem related to keeping stingless bees, more than half of all interviewed beekeepers pointed at the current legislation and the management skills (Fig 7). Only 13% of the interviewed beekeepers exploited alternative bee products other than honey and colonies, including propolis, pollen, boxes and wax (Fig 8). Finally, when asked to assess the status of wild stingless bee populations, 92% of beekeepers responded they believe there are now fewer wild bees than fifty years ago (this proportion also holds when considering only beekeepers older than 59 years).

Discussion
Survey data from 251 beekeepers scattered across Brazil revealed the influence of specific management practices and other confounding factors over productivity and income indicators. Specifically, beekeeping experience, networking with other beekeepers, and technical skills related to bee management were found to be important predictors of productivity and income, regardless of species-specific differences in management. Even though we focused on the general influence of management practices across species, we would like to stress out the importance of optimizing species-specific management practices. Although not many stingless bee species are actually managed [1], there are more than 400 recognized species in the Neotropics [41], each showing a different life history and distribution range [12,28]. Further efforts are needed to optimize husbandry of each managed species and facilitate the commercial use of other still unmanaged species. Our results could thus be used as a general guideline to optimize stingless beekeeping, bearing in mind that particular practices might be more or less effective across different species or groups of beekeepers. We also believe that extension programs aiming to optimize stingless beekeeping should consider local or indigenous traditions and beliefs, to enrich these efforts and preserve cultural diversity.
Although we collected data from beekeepers of 20 Brazilian States with varying socio-economic backgrounds, our data set should not be considered a representation of the population of stingless beekeepers in Brazil. There are about 1800 stingless beekeepers registered in the most popular social network (a Yahoo group), but the total number of stingless beekeepers in Brazil is estimated to be much higher (above 5000 according to the administrators of that group), as many live in remote rural areas without internet access and have never been registered by any agency. Nonetheless, we believe our data is suitable to model the relationships between management practices and production/income indicators. For instance, we did not find a strong dependency structure among predictors (S2 Fig.), so our models are appropriate to    assess the response of stingless bees to different management practices. Moreover, our study is the most extensive assessment of stingless beekeeping undertaken in a developing country [15,31,35,42], and the first large-scale effort aiming at optimizing stingless beekeeping for honey/colony production based on quantitative data.

Multiplication of colonies
Beekeepers can acquire new colonies either by collecting nests in the wild, by capturing new swarms using trap-nests [43], by buying them from other beekeepers, or by multiplying their own colonies. Colony multiplication, whereby one colony is split into two or more colonies, is considered the most sustainable approach to obtain new colonies, since it involves low economic costs (a new box) and does not require the removal of wild bees from their native habitat [1]. Colony multiplication is a key part of stingless bee management, since more colonies are needed to produce more honey or more colonies to sell. Our results show that beekeepers who received a course in meliponiculture are more likely to multiply colonies than those that did not, highlighting the importance of acquiring the necessary management skills to multiply colonies. This is because multiplying colonies of some species is not easy, requiring knowledge on the best moment to do it, and the ability to recognize and carefully extract the brood discs at the appropriate developmental stage. Moreover, newly multiplied colonies are usually weak and vulnerable to infestation by pests [1]. Our results also indicate that beekeepers keeping their bees close to native vegetation and feeding them with sugar or honey syrup, are also more likely to multiply colonies, which suggests an influence of the available feeding resources on colony multiplication success (see below).

Supplementary feeding and honey conservation method
The importance of supplementary feeding has been highlighted in many technical manuals [11,29,44], and it is regarded as an essential part of bee management. Supplementary feeding with sugar syrup or honey has been shown to influence colony growth, performance and survival, both in stingless bees and in honeybees [45][46][47]. Such an influence is corroborated by our data, since beekeepers feeding their colonies with sugar syrup had a larger number of colonies and multiplied more colonies per year (Fig 4A). Our data, however, does not allow separating the relative effect of feeding on colony growth, performance, or survival, because the observed pattern is the outcome of the interplay between them. We hypothesize that in places with a marked seasonality, where flower resources are scarce or completely absent during a given time period, supplementary feeding is particularly important to assure colony survival. The same applies to urbanized areas that have lost most native vegetation. For instance, we found that beekeepers who kept their bees close to native vegetation had more colonies, which could imply that the availability of natural resources influences colony growth and multiplication success. Future studies in arid regions of Northeastern Brazil or severely deforested areas, could further test this hypothesis controlling for other confounding factors not accounted for here. Additional research is also needed to assess the influence of different supplementary feeding diets on colony performance.
Honey conservation method was found to be a key predictor of the amount of honey sold and the beekeeper´s yearly earnings (Fig 4B). It is well known that stingless bee honey ferments much faster than honeybee honey due to its higher water content [48]. For this reason, different conservation methods have been developed to increase post-harvest stability and shelf life of stingless bee honey, including freezing, pasteurization, maturation, and dehumidification or dehydration [49]. Our findings highlight the importance of honey conservation, as beekeepers employing any of these established honey conservation methods sold more honey and earned more money than those that did not. Interestingly, our results also show that beekeepers leaving honey in the refrigerator (which is not considered an established conservation method), sold even less honey and earned less money than beekeepers that left honey out of the refrigerator. This finding could be explained by the beekeepers or their families being more tempted to eat the honey when found in the refrigerator, thus reducing the stock and increasing its exposure to potential contaminants.
Inspection frequency, honey harvest method, and the use of vinegar to control parasitic flies Phorid flies (Diptera, Phoridae) constitute one of the most devastating pests of stingless bee colonies [11,15,50]. Attracted by the odors emitted by stored pollen, the flies enter colonies and lay hundreds of eggs. These later become voracious larvae that deplete the colony´s food stores, causing a considerable damage and often the total collapse of the colony. To control phorid fly infestations many beekeepers use vinegar traps [1], which consist of small vinegar containers provided with lids pierced with small holes (allowing flies but not bees to enter). Vinegar attracts the flies, which enter the containers and drown. Our results show that beekeepers employing such traps have a larger number of colonies than those that do not (Fig 5C), suggesting a relationship between the use of vinegar traps and colony growth and survival.
Experimental work directly measuring colony growth and survival in colonies with and without vinegar traps is needed to support this pattern.
Although the use of vinegar traps was not found to be among the factors explaining colony losses (Table 2), the predictors explaining colony losses (inspection frequency and honey harvest method) are indirectly related to phorid fly infestation. Beekeepers that check their colonies frequently are able to detect infestations and control them before they become a major problem. In addition, they can also strengthen weak colonies with workers or brood from stronger ones, provide supplementary feeding if needed, or detect when a colony is queen-less. Frequent colony inspections thus seem important to minimize colony losses (Fig 5A). Honey harvest method was also found to influence colony losses. Beekeepers that collect honey by "flipping the box", not only kill a larger portion of the colony´s worker population, but also make the harvested colonies more attractive to phorid flies. This is because this technique involves the rupture and exposure of honey and pollen pots before flipping the box to let the honey drain into a collection container. The harvested colonies are thus left with exposed pollen, which makes them particularly attractive to phorid flies. This can explain why beekeepers employing a syringe or an electric pump to harvest honey, which cause a minor rupture of honey pots and leave pollen deposits untouched, have much lower colony losses (Fig 5B).

Labeling of honey containers and feeding frequency
Economic gains from selling honey or colonies were determined by the total number of colonies kept. This is not surprising, given that more colonies imply more honey and more colonies are available to sell. However, the practice of labeling honey containers was found to affect this relationship in an unexpected way. For instance, the relationship between earnings and the number of colonies was much stronger among beekeepers labeling honey containers than among those that do not employ labels (Fig 6A). This result highlights the importance of labeling honey containers (or investing more effort in branding the product), as a strategy to increase the yearly earning of beekeepers selling honey.
On the other hand, costs associated to keeping bees were found to be influenced by feeding frequency. Because the acquisition of sugar or honey for supplementary feeding implies input costs, our results could be explained by the fact that beekeepers feeding their bees more frequently need to buy more sugar or honey. Interestingly, we also found that this pattern was affected by the age of beekeepers. In general older beekeepers invested more into their bees, regardless how many colonies they kept and how frequently they feed them. Hence, older beekeepers seem to be more involved with managing their colonies, investing more money into things other than supplementary feeding, such as buying more colonies, building better boxes or improving the condition of their bee business.

Selective breeding
Selective bee breeding has had a tremendous impact on apiculture, resulting in healthier and more productive bees [51,52]. In stingless bees, no formal breeding programs have been implemented yet, and there is surprisingly little data on queen rearing and other breeding techniques [34]. Here we show that beekeepers undertaking some selective breeding have colonies that produce more honey (Fig 5D). It should be pointed out, however, that we did not assess whether beekeepers preferentially multiplied a selection of their colonies based on desirable traits (such as honey production, colony health, worker population, etc.). We only recorded if beekeepers were able to keep track of the genetic origin of their colonies by avoiding to mix brood from different colonies during colony multiplications. Mixing brood is a common practice, since the easiest way to strengthen a weak colony is offering it fresh brood from another (usually stronger) colony. Such mixing can result in the replacement of the original queen by a new queen, born from the newly introduced brood disc, thus making it impossible for the beekeeper to track maternal lineages. Our results should thus be interpreted as an association between the ability to maintain and track such maternal lineages and colony honey production [but see [53]].

Beekeeping experience and network of known beekeepers
Only a minority of the surveyed beekeepers sold honey (30%) or colonies (25%), a fact that illustrates that stingless beekeeping is still mainly a non-commercial practice. Motivating beekeepers to sell their products is essential to transform meliponiculture into a profitable activity. We found that the years of beekeeping experience and the number of known beekeepers are key determinants of whether beekeepers sell colonies or honey (Table 1). Years keeping bees was found to be an important predictor of the total number of hives kept, as well as the number of colonies sold and the liters of honey sold (Fig 3A). Beekeepers with more experience are likely to have acquired a broader set of technical skills throughout their lives. For instance, many beekeepers make and test different boxes, and try different sugar syrup recipes, feeding techniques, pest control methods, and many other management practices. More experienced beekeepers thus have had more time to test different methods, and establish the most effective set of management practices [15]. Likewise, more experience beekeepers have had more time to build a larger network of known beekeepers (S2 Fig.). Such network is important, as we show that the number of known beekeepers was associated to the total number of colonies kept and to the number of multiplied colonies per year (Fig 3B). The network of known beekeepers thus seems to be a source of technical knowledge on how to multiply colonies, as well as a source of colonies. Indeed, many beekeepers exchange knowledge and colonies, and some leading beekeepers are well known as promoters of meliponiculture, offering free hives and training to starters and amateurs across vast regions [11,54].
Interestingly, years of experience but not the network of known beekeepers was found to predict the number of colonies sold, the liters of honey sold, and whether beekeepers sold honey or not, suggesting that technical skills are the limiting factor determining sales. For instance, beekeepers that participated in a meliponiculture course were more likely to sell honey (in any amount). This explanation is further supported by the fact that demand for stingless bee products often exceeds supply [3,16], so sales seem to be limited by the beekeeper´s ability to produce enough honey or colonies. However, beekeepers cultivating crops on their property were found to be more likely to sell honey, indicating a possible role of the exposure to other agro-markets. The lack of an association between the number of known beekeepers and the probability of selling honey, can also be explained by the fact that contacting other beekeepers is more difficult in rural areas than in urban centers (honey sales were indeed found to be more common in rural areas). On the other hand, beekeepers with a larger network of known beekeepers were found to be more likely to sell colonies, a result that suggests that beekeepers are trading colonies among themselves. Indeed, many beekeepers confirmed that they had sold their colonies to other fellow beekeepers (S3 Table). Because the current legislation restricts the commercialization of colonies (see below), beekeepers might be more tempted to sell colonies to trusted beekeepers.

Main problems and opportunities
When asked to identify the main problem related to keeping stingless bees, more than half of all interviewees pointed at the current legislation and the lack of management skills (Fig 7). This is not surprising since the current Brazilian legislation severely restricts the commercial use of native bees and their products, imposing compulsory registration for stingless beekeepers with more than 50 colonies or those who wish to commercialize any bee products, and outlawing the unauthorized transportation of colonies [55]. Registrations and applications for transport permits are cumbersome and require beekeepers from remote rural areas to travel to urban centers. Currently, there is a movement to change this legislation, and some states have taken the lead in implementing regional regulations [56]. However, a change at the federal level is urgently needed to ease these restrictions, and acknowledge stingless beekeeping as a conservation and development tool. On the other hand, the need for more efficient management skills calls for more research devoted to optimize management practices, as well as more extension work to make the knowledge resulting from such research available to beekeepers [1,31].
We foresee an enormous potential in the stingless beekeeping industry. Once management practices are optimized and the activity is decriminalized, we expect that the market for stingless bee products will experience a significant expansion. Although stingless bee honey is frequently commercialized as medicine in developing countries like Brazil, Mexico and Bolivia, there seems to be a growing number of consumers that regard the product as a luxury item and pay much more for it than for conventional honeybee honey [16]. Moreover, the market for organic or special honey, pollen, wax and propolis is growing across Europe and the US [57]. Likewise, there seems to be a high demand for stingless bee colonies to supply other beekeepers, hobbyists and research/education institutions [17,58]. Although some of the interviewed beekeepers (13%) already exploit alternative bee products such as propolis, pollen, boxes, and wax (Fig 8), only one beekeeper offered greenhouse pollination services. In contrast, many stingless beekeepers in Australia already provide commercial pollination services [17]. We expect the number of stingless beekeepers offering pollination services to grow in the near future, as demand for pollination services increases [21].
The assessment of the status of wild stingless bee populations made by the interviewed beekeepers (92% believe there are now less wild bees than 50 years ago), supports recent reports on global bee declines [18,59,60]. We must point out that this was no quantitative assessment, as it was based entirely on the beekeeper´s opinion. However, we believe this opinion is informative, given that most beekeepers have a close relationship with their bees, constantly asses the natural resources they use, and frequently collect wild colonies. Many beekeepers map the location of wild colonies, and some of them keep track of a remarkable number of them. Our findings thus imply that many wild stingless bee populations might have declined during the last decades, although the factors driving these declines remain unknown. Like honeybee beekeepers, who have helped maintain honeybee populations in places were wild honeybees have all but disappeared [19], stingless beekeepers could help maintain or increase stingless bee populations in many tropical regions. Incentivizing colony multiplication by stingless beekeepers could be a strategy to increase bee populations and minimize the need to extract colonies from the wild.

Conclusions
Our findings have important implications for the stingless beekeeping industry. We identified particular management practices, which could help many beekeepers produce and sell more honey and colonies and earn more money, thus making the activity more profitable. Specifically, our results highlight the importance of teaching beekeepers to inspect and feed their colonies, how to best multiply them and keep track of genetic lineages, how to best harvest and preserve the honey, how to use vinegar traps to control phorid fly infestations, and how to add value by branding their honey. Future experimental work could further test the patterns documented here. Our work also stresses out the importance of beekeeping experience, suggesting that future extension programs could profit from engaging more experienced beekeepers in the design and implementation of technical courses. Importantly, our results emphasize the key role of the network of known beekeepers as a source of colonies and technical knowledge. We believe that the creation or consolidation of local and state associations of stingless beekeepers could enhance such networks, facilitating communication between beekeepers, and unifying efforts to push forward the commercialization of stingless bee products. Our work underlines the need for more research devoted to optimize management practices, as well as more extension work to transfer the generated knowledge to beekeepers. Such efforts could help transform stingless beekeeping into a powerful tool to achieve sustainable development, helping low-income communities improve their living conditions, contributing to the conservation of wild bees and plants, and assuring crop pollination services.   Table. Summary statistics for all variables. N tot represents the total number of beekeepers providing an answer for each question (not all beekeepers responded to all questions), while N reply represents the number of beekeepers in each reply category (only for categorical variables). (PDF) S4 Table. List of the main species kept by 246 Brazilian stingless beekeepers, their common names in Brazil, and the number and proportion of beekeepers rearing each one. Although beekeepers also reared other species, they were not included in this list since they were not registered as the main species kept. (PDF) S5 Table. Regression estimates, standard errors (SE), confidence intervals (CI), z or t-values, and p-values for all best models. Estimates are provided for each level of the categorical predictors. Note that all predictors included in these models significantly increased the modeĺ s likelihood at p < 0.05 (likelihood ratio tests). (PDF)