Body size, richness, and abundance of Staphylinidae unaffected by landscape composition and cropping system in a push–pull maize system in Kenya

The change from predominantly natural to agricultural landscapes has been shown to reduce pollination and biological control services and also affect the functional traits of the insects that provide such services. Research shows that increasing agricultural area at the landscape scale can have variable effects on functional traits such as body size in predators. Moreover, local factors such as farm management strategies have also been shown to affect functional traits in insects. Studies of local or landscape effects on functional traits are often contradictory in temperate regions, and absent from tropical regions. In Kenya, local pest management technologies, e.g. push–pull systems have been shown to increase biological control, but data is lacking on how local management could affect body size of natural enemies. We investigated how the proportion of natural areas in the landscape and local management (push–pull) affected the body size of rove beetles (Staphylinidae) in Kenyan maize fields. Elytral and pronotal measurements were used to quantify body size, and a mixed-effects model was used to test the effect of landscape complexity and local management on the body size, richness, and abundance of rove beetles. We found no effect of landscape composition or local management on body size, richness, or abundance of rove beetles. We hypothesize that the lack of landscape effects could be caused by less plasticity in functional traits in tropical versus temperate zones. Our results demonstrate that changes in land use and local management do not always lead to changes in the diversity and functional traits of predator communities.


Background
As global need for food and farmland increases, the way humans utilize the land is changing. Ecosystems worldwide are becoming increasingly homogenous due to increases in permanent pastureland and large-scale agriculture (Benton et al. 2003), and pastureland and arable crops together currently occupy 39% of the global terrestrial area (Alexander et al. 2015). Land use change has been identified as a main driver of arthropod biodiversity loss worldwide (Sánchez-Bayo and Wyckhuys 2019) negatively affecting community assemblages (Tscharntke et al. 2005;Gámez-Virués et al. 2015), and the ecosystem services they provide (Bianchi et al. 2006;Dainese et al. 2019;Grab et al. 2019). One important ecosystem service in agricultural systems is biological control. Biological control has been quantified as providing benefits valued as $572/ha globally and at $5.5 billion in the U.S. (Naranjo et al. 2015). Its promotion can reduce pesticide application by 70% and increase yields by 5% (Gurr et al. 2016). In agroecosystems, predation by natural enemies such as beetles is important for keeping pest pressure below economic injury levels (Letourneau et al. 2009;Rusch et al. 2015). Land use change and agricultural intensification is also known to affect natural enemies and biological control (Chaplin-Kramer et al. 2011;Veres et al. 2013;Zhao et al. 2015;Rusch et al. 2016). In general, (semi)natural habitats have been hypothesized to provide more floral resources, habitat refuges, protection from pesticides and alternate prey to sustain natural enemy populations (Geiger et al. 2010;Haan et al. 2019). Therefore, it is not surprising to find that landscapes dominated by agriculture have decreased richness and abundance of natural enemies (Dainese et al. 2019;Chaplin-Kramer et al. 2011;Veres et al. 2013) and less effective biological control (Dainese et al. 2019;Zhao et al. 2015).
The composition of the landscape is important not only due to its direct effect on natural enemy diversity and abundance, but also due to its potential interaction with local management strategies (Tscharntke et al. 2005;Batáry et al. 2011). One such strategy is the push-pull system developed by the International Centre of Insect Physiology and Ecology (icipe) in Kenya to control stemborers (Cook et al. 2007). In Kenya, lepidopteran stemborer pests frequently attack small maize fields-one of the main staple crops-causing as much as 88% yield loss during the worst infestations (Midega et al. 2015). In response to the devastation caused by stemborers, an integrated cropping system known as push-pull was implemented in the late 1990's. It consists of a maize field intercropped with Desmodium which acts as a deterrent (push) and surrounded by a border of Napier grass or Brachiaria sp. which act as attractant plants (pull). Push-pull farming has been successful in increasing yield and reducing stemborer damage and has even remained effective in the presence of the recent invasion by fall armyworm (Midega et al. 2017). Push-pull fields also have increased plant diversity and higher percent ground cover which has been shown to benefit natural enemy populations (Rivers et al. 2016), and specifically increase the abundance of generalist predators (Kebede et al. 2018).
Studies and meta-analyses have shown a strong interaction between farm management and the surrounding landscape (Grab et al. 2019;Batáry et al. 2011;Poveda et al. 2019;Perez-Alvarez et al. 2019). A 2011 meta-analysis by Batary et al. demonstrated that local agricultural management schemes such as flower strips and beetle banks were more effective at increasing biodiversity in simplified landscapes. This is because a farm surrounded by a high percentage of natural or semi-natural area (> 20%) is likely already experiencing the ecological benefits of higher arthropod biodiversity, and so follows a law of diminishing returns in regard to human-implemented increases in plant diversity (Batáry et al. 2011). However, most of this research comes from temperate regions, and data on natural enemies in tropical regions, particularly Sub-Saharan Africa, are largely lacking (but see (Midega et al. 2014;Mailafiya 2015;Ndakidemi et al. 2016;Stein et al. 2018)). Given that most of this work has been performed in temperate regions, we wanted to investigate if this hypothesis of landscape interactions with local management practices holds true for the push-pull maize cropping system in Kenya.
Furthermore, it is not just the abundance and diversity of natural enemies that is affected by land use change. Functional traits such as body size have also been shown to be affected by a decrease in the proportion of natural area (Perez-Alvarez et al. 2021). Functional traits are physical traits of organisms that have high variability and low heritability and can be influenced by environmental factors such as resource availability. They are advantageous for thinking about how organisms respond to changes in the environment in a way that directly impacts ecosystem functioning (Poff et al. 2006). Body size is a useful functional trait due to its variability and elastic response to resource availability (Gámez-Virués et al. 2015;Rusch et al. 2015;Perez-Alvarez et al. 2021;Renauld et al. 2016) and its ease for measurement. Functional traits can vary at the interspecific and also at the intraspecific level. For example, Gaméz-Virués et al. (2015) found interspecific variation where arthropod body size decreases with increasing proportion of natural area, a result supported by a recent study by Perez-Alvarez et al. (2021) on Carabidae. Body size can also be used to predict pest control in generalist predators, as demonstrated by Rusch et al. (2015). In generalist predators such as Carabidae, an increase in the body size of the community is positively related to predation rates in agroecosystems due to the increased abundance of a few larger-bodied species and the decrease of many smaller species (Perez-Alvarez et al. 2021).
We decided to study the effect of land use change and local management on functional traits in rove beetles (Coleoptera: Staphylinidae), both at the community level and for the most single abundant genus, given that this is a highly diverse group containing many generalist predators (Irmler et al. 2018), and its relatively high abundance in tropical systems. In comparison with Carabidae (Perez-Alvarez et al. 2021; Kotze and O'Hara 2003), there is a distinct lack of literature pertaining to body size in Staphylinidae (but see (Rusch et al. 2014)). In addition, a recent meta-analysis found that out of 133 publications studying rove beetles, only 10% were done in the Oriental, Australian and Afrotropical regions combined, making this study unique in its biogeographical context (Méndez-Rojas et al. 2021 previous study has found no effect of landscape context on the diversity of staphylinids (Clough et al. 2007), but this was somewhat contradicted by recent findings by Méndez-Rojas et al. 2021, who found that richness decreased in agricultural landscapes, but not in urban areas. Based on previous research stating an effect of landscape composition on body size, we predict that landscapes with a low percentage of natural area will have generalist predator communities (Staphylinidae) with larger body sizes both at the community level and for individual species within the community. We also hypothesize that landscapes with a low percentage of natural area will have lower species richness and lower predator abundance, building on the findings of Perez-Alvarez et al. 2021. Based on the hypothesis that landscape can impact the efficacy of agri-environmental management methods (Batáry et al. 2011), the goal of our study was to investigate the single and combined effects of push-pull local management and landscape composition on body size, abundance, and richness of rove beetles in a Kenyan maize agroecosystem.

Site selection and landscape characteristics:
In 2018, 19 sites were selected in four districts in Kenya: Homa Bay, Bondo, Lambwe, and Rachuonyo ( Fig. 1). The land cover surrounding each farm at a 250 m, 500 m, 750 m, and 1 km radius was classified as agricultural, semi-natural, developed, or open water using QGIS 3.4.12 (QGIS Development Team 2019) and Google Earth imagery from 2017 or 2018 (Google Earth, 2018). Landscape complexity was quantified as the percentage of natural area out of the total land area. Sites varied along a gradient of landscape complexity from 2.9% natural area to 25.3% natural area (mean 10.4% ± 6.5%). All farms were separated a minimum of 2 km from all other sites used in the study, with the exception of two farms that were approximately 400 m apart from each other. Each site had one control and one push-pull field that were located less than 30 m to one another. All fields were predominantly maize, with occasional intercrops of groundnuts (Apios americana), beans (Lablab purpureus), cassava (Manihot esculenta) and greengrams (Vigna radiata). Control fields and push-pull fields were located on the same farm, managed by the same farmer, and at 17 out of the 19 farms they shared one border. Push-pull fields were intercropped with Desmodium sp. and surrounded by a border of either Napier grass (Pennisetum purpureum) (4/19) or signalgrass (Brachiaria cv Mulato-II) (15/19). The fields were planted at the beginning of the long rains season (February-April) by the farmers using their conventional planting practices for control fields and following the guidance of our collaborators at icipe for the push-pull fields.

Specimen collection and measurement
Specimens were collected twice. The first time while maize was in the tasseling stage between May 15 and June 7 2018 and the second time while ears were developing on the plants between June 12 and June 27 2018. Collection was done with pitfall traps and sweep net sampling. Five pitfall traps consisting of a standard SOLO brand 473 ml cup with approximately 100 ml of soapy water were placed near the four corners and in the center of each field and left open for 24 h during both sampling periods. Each sampling period also included a sweep net sample consisting of 25 sweeps taken between two rows of maize near the middle of the field for both the control and the push pull sites. Samples were stored in 70% ethanol. In the second sampling period only, a ten-minute transect walk was performed, in addition to the pitfall traps and the sweep netting, during which all arthropods that could be collected by hand from maize plants were sampled. After conducting preliminary abundance measurements, we found that rove beetles (Staphylinidae) were abundant in the system when compared to other putative predator taxa. Of 3519 putative predators, ants (Formicidae) accounted for 63% of specimens, and are outside of the scope of this paper. Of the remaining 1290 non-ant specimens, the most abundant groups included spiders (38%), parasitoid wasps (19%), Staphylinidae (10%), Coccinellidae (6%). and Dermaptera (4%). Spiders may predate on eggs occasionally but are generally not considered important egg predators and based on rearing of eggs collected from the field, egg parasitism was extremely rare (Luttermoser et al. (in prep), unpublished data). Due to a lack of relevant literature on measuring body size in Staphylinidae and histerids having similarly truncate elytra, our method was based on the work of Stoffolano & Geden (1987) on Histeridae (Stoffolano and Geden 1987). To measure body size, ten initial measurements were performed on 28 specimens belonging to the Paederus genus. The measurements were: pronotal width at the head, maximum pronotal width, right elytral width at anterior end, maximum right elytral length, and the diagonal length of the right elytron from inner left (anterior) to outer right (posterior) points, head width across eyes, maximum pronotal length, and the lengths of the hind femora and tibiae. We selected the traits with the highest variance with the expectation that these traits would be most likely to reflect any variation due to local management or landscape composition (Table 1).
From this analysis, elytra length, elytra diagonal, and maximum pronotal length exhibited the highest variance and were selected to measure on the remaining specimens. Diagonal elytral length was measured from the anterior left corner of the elytron to the posterior right corner in direct dorsal view. Maximum elytral length was measured from the anterior median to posterior median point. Maximum pronotal length was measured from the anterior medial to posterior medial point (Fig. 2). Measurements were taken using an Olympus SZX 10 stereo microscope and DP22 Microscope Digital Camera. Olympus cellSens Standard 1.16 software was used to measure beetles to the nearest 0.01 mm. Specimens were sorted to morphospecies based on observable characteristics such as number and shape of flagellomeres, shape of pronotum, size, and coloration. 13 morphospecies were identified to subfamily based on Ferro et al. 2015. Vouchers were deposited in the Cornell University Insect Collection, CUIC #000,002,452-#00,002,482.

Statistical analysis
We used a mixed-effects model with farm ID as a random effect to test for the effect of local management (control vs push-pull), landscape complexity and their interaction on body size, richness, and abundance. We analyzed effects across all morphospecies, and for the genus Paederus specifically. Paederus was abundant across our sites (29 of 125 total specimens, or 23% of all specimens) and can be easily identified to genus.
We conducted separate analyses for all morphospecies and Paederus specifically to test whether the variation in body size was due to community variation or intraspecific variation within the single most numerous genus. We then ran a linear mixed model to test effects of landscape at the 250 m, 500 m, and 1 km scale using both community-level and Paederus-specific variation data to determine how landscape was influencing traits in the community. Data were analyzed in R-Studio, ver. 3.5.2 (R Core Team 2019) using packages lme4 (Bates et al. 2015) and nlme (Pinheiro et al. 2021).

Results
We found a total of 125 specimens of 13 Staphylinidae morphospecies. 3 out of 19 farms had no staphylinids on either visit. The subfamilies identified were Paederinae, Staphylininae, Aleocharinae, Tachyporinae, Oxytelinae, Osoriinae, and Omaliinae. We found that the percentage of natural area did not significantly affect body size both within Paederus and at the family level (Table 2). Additionally, local management (push-pull vs control) did not have a significant effect on body size (Table 2). We also found no significant effect of the interaction between landscape and local management on body size within Paederus and at the family level (Table 2). Additionally, we found the percentage of natural area, local management and their interaction had   Agriculture and Bioscience (2022) 3:54 no effect on abundance or richness at the morphospecies level (Table 3).

Discussion
Initially we predicted that an increase in the percentage of natural area would increase the richness and abundance of staphylinids but decrease the average body size. We also predicted that compared to a control field, a push-pull field would have higher Staphylinidae richness and abundance, and that this effect of push-pull would be stronger in a landscape with a lower proportion of natural area. Contrary to our predictions, we found no effect of proportion of natural area on body size, richness, or abundance. Additionally, we also did not find any interaction between landscape composition and push-pull and its effects on richness, abundance, or body size. Although our sample size can be considered relatively small, we have seen significant results in other studies with smaller sample sizes (Renauld et al. 2015), showing that the sample size of 112 individuals across 16 sites could still be rendering significant effects if those effects were strong. A larger sample size may have rendered other results, but given that the significance was so low, we trust that staphylinids do not vary as significantly. Moreover, staphylinids were a significant proportion of our non-ant predators, and we felt they were a potentially important predator group in our system and deserved to be analyzed. Despite considerable evidence showing natural enemies respond to changes in landscape composition (Gámez-Virués et al. 2015;Chaplin-Kramer et al. 2011;Veres et al. 2013;Perez-Alvarez et al. 2021;Palmu et al. 2014), neither richness nor abundance of Staphylinidae changed with an increase in the percentage of natural area. Other studies have also found that staphylinids do not exhibit a consistent response to agricultural intensification (Rusch et al. 2014;Méndez-Rojas et al. 2021;Palmu et al. 2014), andTscharntke et al. (2016) reviewed the mechanisms of why natural habitats can fail to support natural enemies, showing that positive effects of natural habitats on natural enemies and biological control are not always to be expected (Tscharntke et al. 2016). Although there are many mechanisms that could be responsible for the lack of landscape and local effects on Staphylinidae richness and abundance, our results suggests that the Staphylinidae community does not seem to be affected by these particular local practices and by the land-use change seen in the area of study, as will be discussed below.
Our result for body size contradicts the work of Gámez-Virués et al. (2015), who suggested that landscape simplification results in larger body size across arthropod assemblages. Possible explanations could be that Staphylinidae are more robust, less plastic, and are able to withstand the negative effects of conventional farming and agricultural intensification. Although the existing literature has found that staphylinid assemblages tend to be negatively affected by anthropogenic disturbance (Caballero et al. 2009;Martins et al. 2013) this response is highly variable among species (Caballero et al. 2009;Work et al. 2013). It is also possible that our landscape gradient is not extreme enough to detect a measurable difference. Percent natural area in our study ranged from 2.9 to 25.3%. Still, studies in the same region on parasitoids have shown that even smaller gradients in landscape composition can have an effect on biological control (Midega et al. 2014), and differences in functional traits between predators and parasitoids can be one of the reasons those different responses are detected. It may be that staphylinids in our system do not experience extreme enough anthropogenic disturbance to have a pronounced effect on either body size or community composition, or it may also be that the composition of Staphylinidae does not show a unified response to simplification, but individual taxa may respond differently. It is surprising that we did not find higher staphylinid abundance in the push-pull fields, given that rove beetles are known to be active on the surface and benefit from diverse ground cover environments due to the increase in decaying plant matter (see Ref. Rivers et al. 2016). Additionally, certain specimens (Paederus) were hand collected from maize, while the rest were caught in pitfall traps. This suggests a possible spatial stratification in the staphylinid assemblage, which may also be contributing to the non-unified response to landscape simplification.
The local management of push-pull also did not affect diversity measures or body size. Although somewhat unexpected, it is not completely surprising since the push-pull and conventional farm treatments are normally adjacent to each other. This proximity of the two management practices potentially also mitigates the effects of landscape simplification. Previous research has shown beetle banks and flower strips to be highly effective at conserving natural enemies (MacLeod et al. 2004), and research by Collins et al. (2002) demonstrated that predation increases in proximity to beetle banks (Collins et al. 2002). The configuration of the push-pull and control fields as directly adjacent to or very near each other may result in a beetle bank effect where arthropod predators and prey from the conventional system can find refuge in the Napier or Brachiaria grass that borders the push-pull fields. This is encouraging information for Kenyan farmers because it suggests that farmers may be able to experience the pest control benefits of push-pull without implementing it across all fields, as suggested by long term analysis of the data (Luttermoser et al. in prep). This could potentially save farmers both time and money while still effectively controlling stemborer pests and providing the benefits of agri-environmental management not just to the managed, but also to the adjacent fields.
Another hypothesis for the lack of effect from land use change and local management on rove beetle communities, is that this study was conducted in the tropics. Most of the work investigating landscape effects on arthropod communities has been done in temperate regions (Rusch et al. 2015;Chaplin-Kramer et al. 2011). However, it is known that the diversity and number of interactions increases towards the tropics (Hargreaves et al. 2019). Maybe this higher diversity makes tropical systems more resilient to disturbance and also makes communities of arthropods more resilient. Although we do not have enough evidence in this paper or the current literature to test this hypothesis, we propose that a meta-analysis would be a useful way to test this hypothesis more broadly, as soon as enough data are available to test it. It would be informative to conduct bioassays with Paederus species and Spodoptera eggs to assess their biocontrol potential. We know that Paederus was a highly abundant genus in our samples, and were hand collected from maize plants. Moreover, both larvae and adults of this genus are known to be generalist predators (Frank and Paederus 1987), so there is a possibility for egg predation in the field. We also suggest further studies on functional traits and cropping systems in tropical regions on other morphologically and taxonomically diverse groups such as Diptera, Hymenoptera, and other Coleoptera. These groups are potentially providing important ecosystem services on farms, and more research should investigate how their functional traits may vary with both landscape complexity as well as local farm management.

Conclusions
Our results show that in tropical agroecosystems local and landscape effects do not affect rove beetle communities in terms of richness, abundance, and body size. Although this does not match predictions generated from temperate regions, we think that this study along with others from the tropics could be the clue to better understand latitudinal gradients of biodiversity and the impacts that land-use change has on them.
The insert shows Kenya and the location of the study sites is shown by the white box.