Linking Coleopteran Diversity With Agricultural Management of Maize-Based Agroecosystems in Oaxaca, Mexico

Biodiversity is known to be influenced by agricultural practices in many ways. However, it is necessary to understand how this relation takes place in particular agroecosystems, sociocultural contexts and for specific biological groups, especially in highly biodiverse places. Also, in order to systematically study and track how biodiversity responds or changes with agricultural practices, it is necessary to find groups that can be used as practical indicators. We conduct a study of beetle (Coleoptera) diversity in maize-based agricultural plots with heterogeneous management practices in the Central Valleys of Oaxaca, Mexico, a region with outstanding biodiversity and a long agricultural history. We use a mixture of local knowledge and multivariate statistics to group the plots into two broad and contrasting management categories (traditional vs. industrialized). Then, we present an analysis of Coleopteran diversity for each category, showing higher levels across different diversity indexes for the traditional plots. Specifically, Coleopteran guilds associated with natural pest control and soil conservation are more common in traditional plots than in industrialized ones, while herbivorous beetles are more abundant in the second. Also, our results let us postulate the Curculionidae family as an indicator of both management type and overall Coleopteran diversity in the agricultural lands of the study site. We discuss our results in terms of the agricultural matrix quality and its role in strategies that favor the coexistence of culturally meaningful agricultural systems and local biodiversity.


INTRODUCTION
Agriculture is a diverse and dynamic process that involves a large variety of production practices. These practices are not assembled at random, but occur as sets of land management activities, "management types, " that are interdependent, adapted to each other and that function as a system with specific goals (Andow and Hidaka, 1989;Vandermeer, 2011). Depending on management type, agriculture can affect in different ways the maintenance of biodiversity at the local and regional level. Indeed, it has been suggested that some types of agricultural management can coexist and even favor biodiversity conservation if it allows for the temporary or permanent establishment of local species in agroecosystems (Perfecto et al., 2009;Fahrig et al., 2011;Kremen and Merenlender, 2018).
However, the relationship between specific agricultural management types with biodiversity and ecosystem services, like pest control and pollination, is not straightforward. On the one hand, agricultural landscapes can rarely be associated with a single and clear management type. In the tropics and subtropics, for example, one regularly finds patchy landscapes with mixed management types that incorporate some industrial characteristics but keep many features of traditional agriculture (Perfecto et al., 2009;Álvarez et al., 2014;Ramos, 2020;Urrutia et al., 2020). It is thus a challenge to organize such heterogeneity into typologies that allow us to test the role of management, as well as of specific agricultural practices, on the maintenance of biodiversity and its associated ecosystem services (Álvarez et al., 2014). On the other hand, the response of biodiversity to agricultural management depends on numerous factors, like the species under study, landscape structure, type of crop, climatological conditions, etc. (Lindenmayer et al., 2012;Wilson et al., 2017;Winter et al., 2018;Qiu et al., 2019). In order to overcome a simplistic view that conceives agriculture as an absolute antagonist to biodiversity, it is thus necessary to assess the effect of agricultural management in different contexts and case studies, and to identify, favor or adapt strategies that allow for its coexistence with biodiversity.
In ecological terms, agricultural managements with high planned diversity and low input dependence are usually associated with heterogeneous landscapes that enable diverse ecosystem services (Perfecto, 2003;Fahrig et al., 2011;Kremen and Merenlender, 2018;Urrutia et al., 2020). Under certain conditions, these types of agroecosystems have been shown to maintain or even increment the biodiversity at a landscape or plot level, exhibiting a high diversity of organisms that are, in turn, good indicators of overall biodiversity and ecosystem functioning (Perfecto, 2003;Gallé et al., 2019;Happe et al., 2019;Qiu, 2019). In turn, it has been suggested that the reduction of plant species in some agroecosystems leads to a decrease in the number and type of habitats suitable for associated biodiversity, including arthropods (Tylianakis et al., 2007;Letourneau et al., 2011;Fabian et al., 2013;Isbell et al., 2017;Liere et al., 2017). Such reduction in plant diversity, and the concomitant reduction in potential habitats, is often driven by industrialized managements in which plants other than the main crop are not grown and are usually not tolerated and landscape structure is simplified (Ekroos et al., 2010;Cizek et al., 2012;Karp et al., 2018).
In particular, the loss of arthropod diversity has been reported to result in the loss of key ecosystemic services, like natural pest control, pollination, or soil conservation (Bianchi et al., 2006;Cizek et al., 2012;Poveda et al., 2012;Wan et al., 2019;Flores-Gutierrez et al., 2020). Within arthropods, in this work we will focus on Coleoptera (beetles), since their diversity and abundance has been found to correlate with other vertebrate and invertebrate taxa (Pearson and Cassola, 1992;Holland, 2002;Ohsawa, 2010), as well as with with key environmental variables (Lassau et al., 2005;Campanelli and Canali, 2012) and agricultural management practices (Kromp, 1999;Holland and Luff, 2000;Burgio et al., 2015;Flores-Gutierrez et al., 2020). Beetles are cosmopolitan, they define a great variety of niches, exhibit diverse life strategies and they take part in several ecological processes making them important matter and energy flux regulators (Hunt et al., 2007;Bouchard et al., 2011). Moreover, beetles have been postulated as bioindicators of agroecosystem functioning (Pizzolotto et al., 2018), as they include both potential pests and natural pest enemies, as well as pollinators and detritivorous species, among others (Hunt et al., 2007;Brooks et al., 2012).
In spite of the key role of beetles within agroecosystems, it remains to be tested if their overall diversity and abundance, as well as that of guilds potentially involved in agroecosystem services or disservices, exhibit differences between different types of management in widely distributed maize-based agroecosystems in the tropics. Addressing this issue may also lead to the postulation of specific groups of beetles as indicators of management, overall coleopteran diversity and agroecosystem functioning in this type of agroecosystem. Maize-based agroecosystems like the Mesoamerinca milpa, which includes maize, beans and other plants, are fundamental in terms of: (i) food production and nutrition (e.g., Bellon et al., 2018;GIAHS-FAO, 2020), (ii) the vast cultural diversity and knowledge associated to their reproduction, as illustrated for instance in their uncountable gastronomic expressions (CEMDA, 2017;Mora Van Cauwelaert, 2017;Lyver et al., 2019;GIAHS-FAO, 2020), and (iii) their potential role in conservation in the ample and highly biodiverse areas in which they have been historically practiced (Myers et al., 2000;Oviedo, 2002;Bellon et al., 2018).
In this work, we sampled the beetle community in plots belonging to a highly heterogeneous, peasant-driven agricultural landscape (Urrutia et al., 2020). We worked in the Central Valleys of Oaxaca, Mexico, a region recognized for its extraordinary biological and cultural diversity, and that is in many ways representative of small-scale agriculture landscapes in southern Mexico (Mora Van Cauwelaert, 2017;Ramos, 2020;Urrutia et al., 2020). In the area where this study was conducted, a relatively diverse scarab beetle community has been previously described, and was found to be as rich in croplands as in contiguous forests (Ramírez-Ponce et al., 2019). Nevertheless, the relationship between agricultural management and the diversity of beetle taxons and guilds, as well as their potential ecosystemic services, had not been assessed prior to our study. Given the evidence summarized above, we expected to find a higher diversity of Coleoptera in plots with traditional practices involving relatively high plant diversity (and habitats) and little use of external inputs, such as pesticides. Within this diversity, we expected to find a better representation of guilds associated to ecosystemic services in more diversified and traditionally managed agroecosystems. We also expected to find specific taxons that could function as practical biological indicators of the management type and the overall coleoptera diversity.

Study Area
This study was conducted at Villa de Zaachila (Zaachila hereafter), in the Central Valleys of Oaxaca, Mexico. It is a semiurban population located 17 km southeast of the state's capital. The history of landscape management of Zaachila begins with the Zapotec peoples, about 3,500 years ago. Today, agricultural plots are mostly managed by small landholders for family or local consumption. Many of these plots represent a stategiven usufruct land, the so-called ejidos, in which variable degrees of collective management occurs (INEGI, 2007;Mora Van Cauwelaert, 2017). The most represented land use types are agriculture, which covers 48% of its total area, secondary forest with 23% and urban zones with 19% (Urrutia et al., 2020). According to government data, the municipality is composed of 1,669 ha which are distributed among 1,521 small landholders (heads of family), which is to say every family has an average of 1 ha of land. This small-scale agricultural scheme combines very diverse practices in a rather fine-scale mosaic (INEGI, 2007; Figure 1). In order to organize this heterogeneous management context into types or categories that were amenable to further analyses, we used a previously reported method to generate management typologies (Álvarez et al., 2014; see Plot Categorization section below).

Plot selection
Using the methodology developed by Álvarez et al. (2014), we looked for 16 plots belonging to contrasting types of management. The methodology seeks to group heterogeneous plots based upon management information provided by local key informants in a reproducible way. It starts with an initial hypothesis, in which plots that are thought to belong to qualitatively different management categories are selected by key informants. Then, thorough information about their management is obtained through interviews and, finally, a statistical multivariate test is performed in order to confirm or refute the hypothesized categorization. In the following we elaborate these steps.
During September 2016 we worked with a group of small landholders, key informants in the community who provided us information about the study site (Marshall, 1996a,b). Together, we set out to find plots belonging to two qualitatively different management types. First, a traditional type, characterized by a low use of commercial inputs and the planting of diverse, local seeds. Even though the exact combination of agricultural practices in each plot actually varied, we postulated that by sharing those two features they would constitute an identifiable category. Second, we looked for the industrialized type, which would be formed by plots with a relatively high use of machinery (for the oaxacan Central Valleys context), agrochemicals and few types of hybrid seeds. With this in mind, and according to their knowledge of local practices, the local small landholders led us to eight pairs of plots composed of what they judged was one plot under traditional management and one industrialized plot. In order to make the sample as representative as possible and to reduce biases stemming from possible landscape heterogeneities, we required plots in each pair to be as close to one another as possible (when possible even side to side) and for them to cover the North, South, East and West parts of the locality (Figure 1). This assignment of plots to management type represented our initial hypothesis and was later tested with a statistical analysis performed over a set of variables that were obtained from detailed interviews to the owners of each plot, as we will explain in the subsequent sections.

Plot Typology Validation
Still following the methodology by Álvarez et al. (2014), we conducted semi-structured interviews with the owners of each plot in order to collect information on their agricultural management. We asked for information about the planted crops in the current and previous years, use of industrial inputs along the year (fertilizer, pesticides, herbicides, and machinery), rotation techniques, attitude toward non-planted herbs and trees, irrigation use, purpose of the crops (such as family use or market sales), relative economic importance for the family, and so forth. The complete list of variables, definitions and data type can be read in in the Supplementary Material (Supplementary Table 2).
Interviews lasted between 1 and 2 h and they were recorded in audio and summarized in written questionnaires. Some of the variables were also systematically verified by sight in the plots, such as the presence of trees or herbs and the number of crops planted at the current time, but most were only obtainable through the interviews since they involved actions taken in the previous months or other non-observable aspects. Prior to this process, an informative letter was given to each of the participants and all of them gave their informed consent to participate in the study. Because many plots were not managed by single individuals but by families, we interviewed more than one person by plot when it was possible, mainly couples who were together in charge of the plots and households. In total, we interviewed 21 people, with ages ranging from 39 to 78 years old (average age was 60, standard deviation 10.9 years), composed of 16 men and 5 women.
Out of all the information obtained, we chose those management-related variables for which we had clear and unequivocal responses from all interviewees. Two variables (number of crops and number of crop varieties) were recorded as absolute richness values and the rest were recorded as presence/absence values. We are aware of the loss of information that this process implied but we hold it necessary for having a comparable set of data for all plots and to proceed as suggested by Álvarez et al. (2014). This gave a result of 22 variables, which were poured into a quantitative table that can be found in the Supplementary Material (Supplementary Table 3). We used this table to test the initial categorization of plots proposed jointly with key informants (traditional and industrialized). We created a correlogram of Spearman coefficients in order to look for relationships among management practices. Then, due to the type of variables in our study, we used a Factor Analysis of Mixed Data (FAMD) for organizing the plots according to their management practices. This allowed us to confirm or refute the categorization hypothesized by our key informants in order to form the groups that were later used as treatments to compare beetle diversity. All tests were done using the FactoMineR and corrplot packages from the coding language R 1.1.383 (R Core Team, 2014). All scripts used for this investigation are publicly available at: https://github.com/laparcela/ColeopteraZaachila.

Diversity of Coleoptera
Plot extension was around 50 × 200 m, and we established five quadrants measuring 1.5 × 1.5 m in each one. In order to randomly locate the quadrants within plots, we first established a random point in each of the long sides of the plot and drew an imaginary line connecting them. At this line's center, we established the first sampling quadrant. Then we established another two quadrants halfway between this center and each of the line's extremes. Afterwards, a second imaginary line was drawn perpendicularly to the first line, crossing its center. The fourth and fifth quadrants were established at the extremes of this second line. Following this procedure, we obtained five quadrants within each plot, three at the plot interior and two at the plot borders (A diagram can be found in the Supplementary Figure 1). Sampling was made this way because it was done as part of a broader study and this design fitted all sampling purposes (work in preparation). Beetle collection was made before any other sampling in order to avoid disturbing individuals.
During September 2016, toward the end of the rainy season and when crops were almost ready for harvesting, all plots were sampled through the use of sweep nets in the five quadrants. Two plots were sampled each day, starting at 7:00 a.m. and until ∼11:00 a.m. This was kept constant in order to avoid extra noise from a day's normal temperature variation. Due to time limitations, we only had one sampling date per site, which represents a shortcoming for our study. Likewise, it is important to note that the use of sweep nets biases the sampling toward beetles that are active in the above-ground part of the agroecosystem, mainly on stems, leaves, flowers, fruits or sometimes flying in the air. Collected specimens were fixed with 70% ethanol and specimens from the five quadrants were pooled together in a single sample for each plot and then taken to a laboratory for identification. Individuals were identified to family level using identification keys (White, 1983;Tripplehorn and Johnson, 2005) and to morph level afterwards. Diversity was quantified as abundance and richness at both family and morph levels. The complete dataset for this investigation is publicly available at: https://knb.ecoinformatics.org/view/urn %3Auuid%3A0b5a89ef-93f8-493c-8bff-e4843b9a45b9.
Diversity in the industrialized and the traditional plots was compared using one-way ANOVAs, Renyi diversity analysis [using the vegan and BiodiversityR packages from the coding language R 1.1.383 (R Core Team, 2014)] and direct quantifications at several taxonomic levels. Following the ANOVA tests, which were performed using each plot as a separate sample, we pooled all plots belonging to the industrialized group in one, and all plots belonging to the traditional group in another sample, and compared their overall richness. A Renyi diversity analysis was performed on these compound samples. Renyi's analysis relates in a graphical way many widespread ecological indexes, which are located along the horizontal axis while the vertical axis shows the corresponding value for Renyi's index (H-alpha) 1 . Afterwards, we grouped the families into four trophic guilds in order to assess functional differences between types of management. The guilds comprised herbivores, predators, detritivores, and polyphagues, the latter including cases when individuals are able to exploit a wide variety of feeding resources and cases when juveniles and adults belong to different categories. The complete table of families belonging to each guild can be found in Supplementary Table 4 (based on the public dataset mentioned above). We also performed a Canonical Correlation Analysis (CCA) in order to test the relation between a chosen set of management variables and a set of beetle families. In order to reduce the number of variables and meet the test's requirements, we chose those management practices whose contributions to our PCA's first two dimensions were above a reference value that corresponds to the expected value if the contributions were uniform. Variance-inflation factors were calculated to eliminate variables causing multicollinearity. Additionally, we focused only on the most abundant families inside each trophic guild. This procedure left us with 8 management variables: number of crops, irrigation, quelites, trees, compost, and industrial fertilizer, pesticide and herbicide. These were tested for correlation with 9 beetle families: Cleridae, Coccinellidae, Carabidae, Chrysomellidae, Curculionidae, Nitidulidae, Anobiidae, Cantharidae, and Phalacridae.
Finally, we looked for a beetle family that could be a useful indicator of the agricultural management type. We required it to be non-exclusive of a particular type of management, but to show a clear and different response to the management categories. In order to test its potential as an indicator, we analyzed the relation between its abundance and (i) plot management and (ii) overall beetle diversity ( Figure 8B). For this, we adjusted a second degree polynomial to its abundance in each plot (we took the natural logarithm of this measure in order to comply with the analysis assumptions) against its position along the management gradient (the coordinate each plot occupied in the horizontal axis of Figure 3A). We also analyzed the relation between such family's abundance and the total richness of beetles in each plot.

Plot Categorization
The correlogram showed relations among sets of management practices, as expected from the theoretical postulation of management types (Vandermeer, 2011;Álvarez et al., 2014). We found a positive correlation (Figure 2, groups of blue circles) between practices commonly associated to a traditional management, such as managed border, use of manure, animal traction, number of crops, quelites (non-planned herbs that  are frequently useful to small landholders' families in different ways), number of varieties, presence of trees and use of organic fertilizers. Likewise, there is a strong positive correlation between practices associated with an industrialized management, such as use or industrial fertilizers, herbicides and pesticides, irrigation and monoculture. We also found negative correlations (Figure 2, groups of red circles) between irrigation and fallow time, and subsistence agriculture and monoculture.
The FAMD we conducted explained 46.2% of the total variance in its two first dimensions (27.4 and 18.8%, respectively). Figure 3B corresponds to the variables factor map, where we show that the most correlated variables with the first dimension were the number of varieties within crops (0.896 correlation), number of crops (0.872 correlation), and the monoculture planting scheme (−0.722 correlation); while the second dimension was most correlated with the presence of a fallow period (0.9 correlation), irrigation (−0.84 correlation) and use of industrial pesticides (−0.62 correlation) (see Supplementary Table 2 for a more detailed definition of each variable). Moreover, all the variables with a positive correlation to the first dimension were practices associated to a traditional management (Altieri et al., 1997): use of local seeds, crop rotation, fallow period, animal traction, presence of quelites, managed border, presence of trees, use of green manure, compost regular manure, and organic inputs. On the other hand, the variables with a negative correlation were all associated with an industrial model of agriculture (Perfecto et al., 2009): a monoculture scheme, irrigation and use of industrial inputs (fertilizers, herbicides, and pesticides).
The first dimension of the FAMD thus depicts a management gradient that goes from the most industrialized on its negative side to the most traditional on its positive side ( Figure 3A). Based upon this result, we separated the plots falling to the right and left side of the first dimension's origin, thus discretizing the management gradient into two separate categories (industrial and traditional). By doing this, we found that 15 out of 16 plots fell in the expected side from the management category initially hypothesized by key informants. Thus, we sampled and compared coleopteran diversity among these two groups, leaving aside the one plot that did not match the original hypothesis. Note that the vertical axis was not used to further divide the plots in more groups, nevertheless it is worth saying that the variables with higher coefficients in this dimension were fallow period with a positive sign, and irrigation and pesticide use with a negative sign.

Beetle Diversity
In Figure 4 we present an overview of the sampled diversity of families in the order Coleoptera. In total, we found 1,168 individuals belonging to 25 families, which we then identified as 80 different morphs. The most abundant families were Chrysomelidae, Cantharidae and Curculionidae, with the first two being well-represented in both management types and the last mainly in the traditional plots. We found one exclusive family to the industrialized management: Trogossitidae. Likewise, there were four families exclusively found in the traditional plots: Bostrichidae, Buprestidae, Languriidae, and Tenebrionidae.
Overall, coleopteran abundance did not show significant differences among plots from the industrialized and the traditional groups, but both family and morph richnesses were significantly higher in the traditional plots (p = 0.0265 and p = 0.0203, respectively) ( Figure 5). Figure 6A depicts the results of pooling together all plots belonging to the same management category. We present morph richness inside each family in these two compound samples. Renyi diversity was higher for the traditional category, as it was significantly more diverse across all the Renyi's spectra of indexes (confidence intervals at 95%) ( Figure 6B).
Functional groups varied strongly between management types. Traditional plots had 30% less abundance of herbivorous families, 27% more predators, 78% more detritivores, and 7% less polyphages. Herbivores were dominated by Chrysomelidae and Curculionidae; predators were mainly Cleridae, Coccinellidae and Carabidae; detritivores were dominated by Nitidulidae and Anobiidae, and polyphagues were mainly Cantharidae and Phalacridae.
The CCA showed that management variables accounted for 73% of the variance in beetle community composition, with the strongest predictor variables being: industrial herbicide, followed by industrial pesticide, quelites, industrial fertilizer, and number of crops (p < 0.05). Figure 7 shows that the most correlated pairs are Curculionidae with number of crops, and Cantharidae with the use of industrial fertilizer and herbicide, while families closer to the center of the graph show weaker correlations. The most abundant predators, Cleridae, correlate the most with the presence of quelites; while the most abundant herbivores, Chrysomellidae, correlate the most with industrial herbicide and fertilizer. Nitidulidae, the most abundant detritivores, correlate better with both quelites and number of crops. The full model was statistically significant (p < 0.001).

Curculionidae As An Indicator of Management Type and Diversity
Curculionidae was the only family that complied with the requirements for working as an indicator: it was non-exclusive of any management type and showed a strong and different response to each one, it was also abundant and relatively easy to identify (Figure 8A), so we investigated its potential as an indicator of management and overall beetle diversity. The adjustment of a second degree polynomial between Cucrculionid abundance in each plot and its position along the management gradient gave a positive relation (R = 0.50428, p < 0.05), suggesting that plots with a more traditional management tend to have a higher abundance of Curculionidae ( Figure 8B). Additionally, in Figure 8C we show the relation between curculionid abundance and the total richness of beetles in each plot. We first took richness at a family level but found no significant results, nevertheless, the same comparison at a morph  level did show a significant positive relation between the two variables (R = 0.70715, p < 0.05).

DISCUSSION
In the first part of our study, we used the method by Álvarez et al. (2014) to define two broad management categories, traditional and industrialized, based upon several agricultural practices. Through the multivariate statistical approach we followed, we identified two variables to be most important for the organization of plots into these categories: total number of crops and total number of varieties grown. From an ecological point of view, crop diversity is involved in many processes both below or above ground. It is also associated with structural diversity (Del Río et al., 2003) and therefore, to the set of habitats available for the biota (feeding, refuge, reproduction spots, etc.) inhabiting or passing through the plots (Tylianakis et al., 2007;Letourneau et al., 2011;Fabian et al., 2013;Isbell et al., 2017;Liere et al., 2017). In a practical sense, this result is important given the complexity of defining different types of management in heterogeneous landscapes. Indeed, both variables can be measured with relative simplicity by academics, technicians or small landholders, and they seem to work as "umbrella" variables for many other practices associated with agricultural management, making them useful guidelines for future work that seeks to study different types of management in the region.
Another remarkable aspect is the difference in the grouping patterns of the industrialized and the traditional plots. As seen in Figure 4, the industrialized plots fall much closer to each other than the traditional plots. This means that even though both categories have a relatively high internal variation, the industrialized plots are far more similar to each other than the traditional ones, reflecting the tendency to homogenization that is characteristic of industrialized agriculture, in contrast with the heterogeneity and context-specificity of traditional agriculture (Gliessman, 2015). This homogenization tendency is clear, even for our study area, where the level and complexity of the machinery and inputs used in the industrialized plots is not as high as in other regions of Mexico, like El Bajio and the North of Mexico. The high variability found inside the traditional group of plots is surprising given the small scale of this work (Figure 1), which further illustrates the heterogeneity that can be found even in this type of landscape (Urrutia et al., 2020). This shows that the traditional management is not static or "archaic, " as there is significant variation in the practices that small landholders adopt through time and space.
As for coleopteran diversity, the traditional plots showed higher values at both family and morph level, regardless of the diversity measure used, and in spite of the fact that many of the industrialized plots were close to them (which could have otherwise brought diversity to industrialized plots). This indicates that coleopterans are significantly sensible to agricultural management even at small scales, a fact worth noting given that local communities are influenced by their surroundings and that their diversity depends on the species pool at a broad scale and on many landscape variables (Duelli et al., 1999;Gabriel et al., 2010). Moreover, management variability inside traditional and industrialized categories was also large, as we have discussed above, which makes significant intercategory differences more striking. In addition to having a higher diversity of beetles in general, the traditional plots had more exclusive species, which points to their importance as reservoirs of rare species.
In terms of trophic guilds, the traditional plots had less herbivores and more predators and detritivores, which results in a more functionally rich community with possible positive repercussions for agriculture. Namely, the bigger number of predators could be controlling the herbivorous populations and preventing potential pests, while detritivorous species could be sustaining healthier nutrient cycling processes. This may play a critical role for production, since it is estimated that 18-40% of crop losses can be attributed to herbivore arthropods (Oerke, 2006;FAO, 2017;Sharma et al., 2017), and predatory arthropods have proven important drivers of autonomous pest control mechanisms (Zhang et al., 2007;Vandermeer et al., 2010;Letourneau et al., 2011;Rusch et al., 2016). Our results are consistent with other studies that have found more natural enemies in diversified agroecosystems than in those with an industrialized management (Tooker and Frank, 2012;Torcat-Fuentes et al., 2018;Wan et al., 2019;Flores-Gutierrez et al., 2020). The reasons for this may include that natural enemies are more dependent on secondary sources of food and refuges or nesting sites, which are available in diverse environments (Bianchi et al., 2006;Chaplin-Kramer et al., 2011;Karp et al., 2018). Additionally, natural enemies appear to be more vulnerable to pesticides than herbivores, on the one hand because they too are killed by them (Theiling and Croft, 1988), but additionally, because their food sources (namely, other arthropods) are periodically depleted by pesticides while herbivores' food sources remain unaltered (Perfecto et al., 2009). The combination of predators and detritivorous beetles may be crucial in providing the ecosystemic services of pest regulation and soil fertility (Power, 2010;Tscharntke et al., 2012), which can in turn translate into a decrease in the need for agricultural inputs (Nabhan and Buchmann, 1997;Naylor and Erlich, 1997).
Moreover, the bigger amount of herbivores that we found in industrialized plots further suggests that this combination of practices is not associated to the alleviation of pest related problems. This has been documented before, with herbivore communities being sometimes unresponsive to management industrialization (Flores-Gutierrez et al., 2020) and other times being actually increased by it (Knight and Norton, 1989;Krauss et al., 2011;Meehan et al., 2011). Indeed, our CCA showed that the most abundant herbivores, the Chrysomellidae, were positively correlated to the use of industrial fertilizers and herbicides and negatively correlated to the number of crops, characteristics that together describe an industrialized management. On the other hand, the most abundant predators and detritivores, Cleridae and Nitidullidae, correlated positively with the presence of quelites and number of crops, while relating negatively to industrial inputs and irrigation, again strongly suggesting that a traditional management benefits this guild. Finally, little can be said of the polyphague guild in our analysis since it comprises families with too wide a variety of eating strategies; identification to genre or species levels would be required in order to refine the categorization of these taxa.
Regarding the family Curculionidae, we found that the abundance of this family increases as management tends to be more traditional and that it also reflects the diversity of coleopterans in general. Nevertheless, this should be taken with caution because it was the case for most, but not all the sampled plots (Figure 3). This said, curculionids have an easily recognizable morphology (mainly because of the anterior projection in their head) and in fact it is safe to say that all small landholders in Mexico are well-familiarized with them, so they can indeed be widely used as indicators of overall  Frontiers in Sustainable Food Systems | www.frontiersin.org coleopteran diversity or agricultural management type. The family is abundant in most of the world, and more markedly in the tropics. They mainly feed on vegetal tissue, though they can also eat lichens, algae and fungi (Morrone, 2014). According to Zimmerman (1994), all angiosperms are probably consumed by at least one curculionid species. However, they are not just herbivores, for various plants depend on curculionids for pollination (Morrone, 2014). This family has been used as an indicator of coleopteran diversity before (see Ohsawa, 2010), but we have no knowledge of it being used in Mexico or in maize-based agroecosystems. Because of its cosmopolitan nature, it would probably be found to be a good estimator of coleopteran diversity in a wide arrangement of environments. On the other hand, as an indicator of agricultural management, Silva et al. (2002) also found this family to be more diverse under a management that would fit into our traditional category.
In the agricultural context, curculionids are often thought of as a pest, because they feed on a wide range of crops (maize, beans, avocado, cotton, rice, etc.) (Morrone, 2014). However, it is worth noting that none of the small landholders interviewed in this study regarded them as an important problem (see Silva Aparicio et al., 2003 for a similar case). In general, they mentioned that curculionids bit some of the stored grains, but the amount was not considered significant and these grains were generally used to feed fowl. Also, many of them said to avoid damage by storing the grains with local herbs and flowers which repel curculionids, while only a few of the more industrialized or market-oriented producers declared the use of industrial pesticides as part of their storage practices. In any case, none of the interviewed small landholders reported significant damage by curculionids in the field.
While our work provides significant insights into the relationship between agricultural management and beetle diversity, it has limitations that should be kept in mind. We only conducted one sampling (during the rainy season), which is restricted considering the high inter-annual, inter-geographic and inter-crop variability found in coleopterans (Finn et al., 1999;Andresen, 2003). For example, in another agroecosystem, Flores-Gutierrez et al. (2020) reported that arthropod abundance of herbivores was more than 100-times greater during the rainy season than during the dry season, so our results must be taken cautiously. Also, the use of sweep nets biases samples toward foliage-dwelling insects, so it would be advisable to use other sampling techniques in further studies, as they all have shortcomings and can complement each other (Doxon et al., 2011). Our sample size was relatively small (16 plots), although the collaboration of key informants and the recollection of a large amount of beetles seemed to compensate for this, as shown by the clear differences between management types and beetle diversity (see Blanco andCastro, 2007 andÁlvarez et al., 2014 for discussions on sampling for qualitative information). Given that this was a field study, we had little control over environmental variables, which made statistical replicates impossible and limited the analysis possibilities. As a counterpart, those results that were statistically significant even under the great amount of environmental noise are probably a reflection of strong tendencies.
A previous work in Zaachila, Oaxaca, has shown that patches with different land use types, such as agricultural patches, grasslands and forests, all harbor a high beetle diversity but that the species arrangement is markedly different in each of them (Ramírez-Ponce et al., 2019). In this study we were able to further differentiate agricultural land uses by organizing the large heterogeneity of management practices into two broad types. From these types, the one characterized by the use of local landraces and little external inputs harbored a larger coleopteran diversity than its more industrialized counterpart. In another work, we showed that Zaachila's landscape is characterized by a high spatial heterogeneity given by small and multiple patches, a feature that is probably shared by other areas in the tropics with similar historic and physical characteristics, and that differentiates them from the more commonly studied landscapes of the temperate regions of the world (Urrutia et al., 2020). We also showed that agricultural patches, especially rainfed-agriculture patches, dominate this rural landscape and that, given their atomized pattern and large connectivity, they play a central role in the potential migration and recolonization among the adjacent forest patches (Urrutia et al., 2020). Together, these two studies suggest that a traditional management of the agricultural patches of Zaachila contributes to a high quality matrix for biodiversity conservation. We believe that the results of this work are relevant for the conservation of biodiversity in other contexts where small landholders are still the principal actors of agriculture and where traditional practices still survive. These and other works support the key role of traditional campesino agriculture in the conservation of biodiversity, including agrobiodiversity and its cultural expressions in landscapes driven by small landholders (Perfecto et al., 2009;Mora Van Cauwelaert, 2017;Bellon et al., 2018), a matter that is crucial since most of the world's biodiversity lies in these regions.

DATA AVAILABILITY STATEMENT
The datasets presented in this study can be found in online repositories. The names of the repository/repositories and accession number(s) can be found in the article/Supplementary Material.

ETHICS STATEMENT
The studies involving human participants were reviewed and approved by Ethics Committee of the Ecology Institute at UNAM. Written informed consent for participation was not required for this study in accordance with the national legislation and the institutional requirements.

AUTHOR CONTRIBUTIONS
MB, LJ-B, CG and TL conceived the ideas and designed methodology. CG and TL collected the data. CG analyzed the data. CG and MB led the writing of the manuscript. All authors contributed to the article and approved the submitted version. LJ-B acknowledged financial support form UNAM-DGAPA-PAPIIT (IA202515) and was grateful for NOT being part of the Sistema Nacional de Investigadores.