Effects of conventional and organic management on plant and insect communities in a traditional elephant garlic crop

Agricultural management has a great influence on biodiversity and its services in agroecosystems. In Europe, a relevant proportion of biodiversity is dependent on low-input agriculture. To assess the effects of agricultural management on biodiversity, in this study we surveyed the communities of arable plants, diurnal flying insects, and pollinators in three conventional and in two organic fields of a traditional Elephant garlic (Allium ampeloprasum L.) crop of the Valdichiana area, in Tuscany (central Italy). The sampling was carried out twice during the season: in spring, during crop growing, and in summer, after crop harvesting. We assessed the effects of the different agricultural management on the richness and composition (species occurrence and abundance) of the three communities using univariate and multivariate analyses. Concerning our specific case study, only plant species richness was significantly higher in organic fields (15.7 ± 2.7 species per plot), compared to conventional ones (5.4 ± 2.3 species per plot). Regarding community composition, only pollinators showed a marginally significant difference between conventional and organic fields. Conversely, the effect of specific fields significantly explained differences in composition of all the investigated groups (plants, total insects, and pollinators). The results suggest that, in our case study, the emerged differences in diversity of the investigated communities were mainly attributable to environmental and management factors related to single fields, more than to organic or conventional farming. Such evidence could be partly due to the very local scale of the study, to the heterogeneity of the surveyed fields, and to the reduced number of surveyed fields. Further investigation is therefore needed.


Introduction
Biodiversity is essential to maintain the functionality of ecosystems, and agricultural areas can have a high biodiversity value (Byrnes et al., 2014;Sutcliffe et al., 2015). In agroecosystems, the services provided by living beings are many and diverse, including for instance recycling of nutrients, regulation of microclimate, soil protection, pest control, and detoxification of noxious chemicals. In the last decades, agricultural ecosystems increased their dependence on external energy inputs due to a reduction in their self-regulation ability, which is largely dependent on the occurring biodiversity (Altieri, 1999). This is increasing the awareness that agricultural management must be set towards sustainable practices (MacLaren et al., 2020;Maraux et al., 2013). Thus, in the European Union, a great attention is being paid to reducing the loss of biodiversity in food production systems (European Commission, 2020aCommission, , 2020b. Arable land is one of the main land use types in agricultural landscapes, and arable plants significantly contribute to biodiversity in these habitats (Marshall et al., 2003). Such species are adapted to the regular disturbance of tillage and other agricultural practices (Holzner, 1978(Holzner, , 1982. Arable vegetation evolved under centuries of interactions between man and nature, acquiring a very characteristic species composition (Oppermann et al., 2012). Balanced and biodiverse arable plant communities are currently acknowledged to provide ecosystem services such as support for pest control, improvement of soil quality, and mitigation of yield losses (Adeux et al., 2019;Blaix et al., 2018;Storkey & Neve, 2018). Agricultural intensification caused a remarkable decline in arable plant diversity in the last decades, especially due to chemical weeding and fertilization, shifts from crop rotation to monoculture, and improved seed cleaning (Fanfarillo, Kasperski, et al., 2019;Meyer et al., 2013;Richner et al., 2015;Storkey et al., 2012).
Plant and insect communities are highly interdependent (Moreira et al., 2016;Nelsen et al., 2018). Insects have strong trophic relationships with vascular plants, that conversely may depend on them for pollination and seed dispersal (Bàrberi et al., 2010;Bronstein et al., 2006;Sakai, 2002), or even for the passive defense against herbivores (Grasso et al., 2015). In arable landscapes, arable plants are particularly important for insects, especially in agroecosystems dominated by wind-pollinated crops, like cereals (Bretagnolle & Gaba, 2015;Kovács-Hostyánszki et al., 2016). They may also support beneficial insects as pest antagonists and pollinators, or act as an alternative trophic resource for phytophagous insects otherwise feeding on crops (Barbercheck & Wallace, 2021;Norris & Kogan, 2000). Like other groups of organisms, insects are experiencing a worldwide decline in agricultural landscapes due to the increasing application of pesticides and fertilizers, to the spread of monocultures, and to the elimination of fragments of natural and semi-natural vegetation (Raven & Wagner, 2021).
Agricultural management can have contrasting effects on biotic communities in different crop types, for instance in annual and perennial crops (Brugisser et al., 2010). At the same time, different taxonomic groups can respond differently to management practices (Lüscher et al., 2014). Organic farming, which tends to be considered more sustainable than conventional farming, showed as well contrasting effects on biodiversity, in relation to different groups of organisms and in different agricultural landscapes (Bengtsson et al., 2005;Brittain et al., 2010;Caprio et al., 2015;Gabriel et al., 2013;Stein-Bachinger et al., 2021). Moreover, the same taxonomic group can be influenced in different ways by organic management regarding species richness, abundance, and composition (Brugisser et al., 2010;Hyvönen & Salonen, 2002;Masoni et al., 2017). Such evidence highlights the importance of assessing the effects of agricultural practices on biotic communities including more than one taxonomic group and exploring different aspects of biodiversity (Brugisser et al., 2010;Batáry et al., 2012;Gabriel et al., 2013;Puig-Monserrat et al., 2017;Raven & Wagner, 2021).
Currently, there is lack of knowledge about the impact of agricultural management on biodiversity in horticultural crops. Some available information is limited to the effects of practices on arable plant species and communities, mostly in oilseed rape, beet, and pea crops (Andreasen & Stryhn, 2008;Fried et al., 2008;Lososová et al., 2004). To fill this gap, in this work we investigated plant and diurnal flying insect communities in three conventional and two organic Valdichiana Elephant garlic (Allium ampeloprasum L.) fields in Tuscany, central Italy, both during crop growing and after harvesting. Our aims were to investigate the effects of i) organic vs conventional managements; ii) local-scale factors summarized through the field identity, and iii) the sampling season, on plant and insect communities in terms of species richness and composition (species occurrence and abundance). Based on past evidence mentioned above, we hypothesize that the lower intensity of agricultural practices in organic management does not necessarily increase species richness or favor more sensitive species in the studied biotic assemblages. Moreover, we speculate that the well-known arable plant species turnover along the season (Fried et al., 2008;Lososová et al., 2004) can be highlighted in our data, with a consequent shift also in insect species composition due to plant-insect interdependence.

Study area
The survey was carried out in five fields in a traditional Elephant garlic crop (Allium ampeloprasum L.) of Tuscany (central Italy), locally known as "Aglione Della Valdichiana". This traditional Elephant garlic was historically cultivated in the Valdichiana area, between eastern Tuscany and western Umbria, and recently experienced a revival. The Valdichiana is a wide valley that underwent major transformations across the centuries, especially regarding drainages for reclaiming cultivated land. Today, it is characterized by intensive cereal, fruit, and horticultural croplands, and the few remnants of naturalness are represented by residual wetlands (Lastrucci et al., 2010).
The five fields belong to the four farmers that joined the "Vero Aglione Della Valdichiana" project, which aims at a thorough geographic characterization of this traditional crop, e.g., through chemical fingerprinting, and includes the survey and characterization of biodiversity in its fields in relation to agricultural practices (VAV Working Group, 2021;Loppi et al., 2021;Vannini et al., 2021). The location of the surveyed fields and the study area are shown in Fig. 1.
The fields are located between 250 and 500 m a.s.l., in very different landscape contexts that span from intensive agricultural areas dominated by arable land to forested areas, through peri-urban complex mosaics. The bioclimate of the study area is partly Mediterranean and partly Temperate sub-Mediterranean, depending on the duration and intensity of summer drought. Mean annual precipitations range between 730 and 800 mm. Mean annual temperatures are around 13-14 °C (Azzari, 2006;Pesaresi et al., 2017). Soils are mainly alluvial, loamy-clayey, subalkaline to subacid (Hengl et al., 2017). Field size varied between 700 m 2 and 9,000 m 2 .

Collection of management data
Depending on the farm, the "Valdichiana" Elephant garlic is cultivated with conventional or organic techniques. Overall, it is a winter annual crop. Tillage, milling, and bulb sowing are carried out in October, and harvesting takes place in June. Fertilization and mechanical or manual hoeing are usually performed during crop growing, sometimes accompanied by pre-emergence chemical weeding. In the past, the crop was grown with traditional low-intensity methods, which are still widely used nowadays.
All the involved farmers compiled a register where they annotated the date and type of agricultural practices they carried out during the cultivation and, if any, after harvesting, including details on the used products (fertilizers, herbicides- Table 1). Full information on agricultural practices and the used products is available in Supplementary File 1.

Sampling design and method
The sampling was carried out at the end of May 2020 and repeated in the same field at the end of July 2020, regardless of the agricultural succession (fallow, mulching sheets, or summer-annual crop). The resampling after harvesting aimed at catching the influence of the different succession types on biodiversity, given the importance of stubble fields and post-cultural fallows as habitats for arable plants and animals and the seasonal variation in arthropods activity in arable fields (Feber et al., 2015;Kuussaari et al., 2011;Pinke & Pál, 2009). We randomly placed 2 × 2 m plots all over the surface of the fields. The randomized design was built separately once per season, so that summer plots were not the same that were surveyed in spring. The number of surveyed plots was higher in bigger fields, and varied between 3 and 9 per season (field a = 9 plots; fields b,c,d = 3 plots; field e = 4 plots). We surveyed a total number of 44 plots, 22 in May (hereafter "spring plots/data") and 22 in July (hereafter "summer plots/data"). In each plot, we sampled vascular plants and diurnal flying insects (hereafter "total insects"). Each survey was carried out at least 15 days after the last agricultural operation.
Regarding vascular plants, all the species occurring in the plot were recorded and attributed a cover value, according to the Braun-Blanquet seven-degree scale (Braun-Blanquet, 1964). Plant community data were stored in the "European Weed Vegetation Database" of the European Vegetation Archive (Küzmič et al., 2020). To avoid the influence of disturbance on insects communities, these were sampled before vascular plants, during sunny and windless days. Targeting flying insects, sampling was carried out by sweeping the vegetation ten times all over each plot's surface, using an entomological net (McCravy, 2018). Net sampling was chosen since being more efficient than other methods of catching flying insects, especially pollinators (Popic et al., 2013). After net sampling, we integrated data collection by means of sight hunting and catching, both on the ground and on vegetation, for 15 min per plot. All the entomological samplings were made by the same person. Aphids were not sampled. The number of individuals per each insect taxon collected or observed in the plot was counted. The collected insects were stored in ethanol. Despite the reduced sampling effort, this method allowed to quickly obtain comparable observations on the entomological communities from all the surveyed plots. All the collected samples were identified in the lab at the lowest possible taxonomic level based on Pignatti et al. (2017Pignatti et al. ( -2019 for plants, and several sources for insects (Supplementary File 1). The nomenclature of vascular plants was updated according to Bartolucci et al. (2018) and Galasso et al. (2018). The nomenclature of insects follows de Jong et al. (2014). Given the lack of a unique reference, pollinator insects were distinguished through an experience-based approach, i.e., an expert entomologist classified insects into pollinators and non-pollinators based on his knowledge and experience.

Calculations and statistical analyses
Differences in plant, total insect, and pollinator species richness were assessed by means of Generalized Linear Mixed Models (GLMMs) (Bolker et al., 2009) using a Poisson distribution as link function after checking that the data were not overdispersed (function "dispersion test" in the package "AER" of R-project) (Kleiber & Zeileis, 2008;R Core Team, 2020). The GLMMs were run through the function "glmer" in the package "lme4" (Bates et al., 2015). Management type and season were used as fixed factors, while the field was included as a random factor. This model type was chosen after comparing it, through likelihood ratio test, with a more complex one in which the factor "field" was nested in management type. However, compared models did not differ, allowing us to choose the simplest model for all the taxonomic groups.
Permutational Analysis of Variance (PERMANOVA) was performed in the software PRIMER v.6 (Clarke & Gorley, 2006) to test whether management type and the season significantly affected species composition of vascular plants, total insects, and pollinators. One PERMANOVA was performed separately for each of the three groups. PERMANOVA correctly calculates an appropriate pseudo-F statistic for each term in the model for multivariate or univariate datasets. Moreover, the permutation approach is free from many of the assumptions of parametric statistics  (3), manuring (2), manual hoeing (3),mulching sheets (1) C 3000 m 2 Complex peri-urban agricultural mosaic 250 m a.s.l. (Anderson et al., 2008). The effects of management type (2-level fixed effect: organic and conventional), field (5-level random effect nested within the management type factor), and season (2-levels fixed effect: spring and summer) were analyzed, as well as the interactions between these factors. The analyses were based on Bray-Curtis dissimilarities of ln(x + 1)-transformed abundance data (percentage cover for plants and number of individuals for insects). The following settings were used for all the tests: • 999 permutations of residuals under a reduced model and α = 0.05, which yield the best power and maintain the most accurate type-I error for multi-factorial designs (Anderson & Ter Braak, 2003); • type III Sum-of Squares, which is the default option in PERMANOVA and the most conservative one (Anderson et al., 2008).
The differences in community structure (presence and abundance of species) among the management types were visualized by non-metric multidimensional scaling (NMDS) using Bray-Curtis dissimilarity in the program PRIMER v.6. We plotted in the NMDS graphic species with correlations ≥0.4 (Spearman rank correlation).

Description of the surveyed fields
Among the surveyed fields, there was a high variability in the mean cumulative cover of vegetation (sum of the percentage covers of all the plant species occurring in each plot, calculated to provide a measure of the degree of vegetation stratification). Species richness per field varied a lot as well, especially regarding plant communities. Table 2 shows descriptive statistics summarizing the main features of the surveyed fields.

Total plant and insect taxonomic diversity
We detected a total amount of 93 plant species, 97 insect taxa, and 138 insect individuals. Plants were all identified to the species level. Regarding insects, 59 were identified to the species level, 24 to the genus level, 13 to the family level, and one could not be identified at any taxonomic rank. The 38 insect taxa that could not be identified at the species level could anyway be recognized as different from one another so that they could be treated as species in the analyses. Plant richness varied between 1 and 21 species per plot. Insect richness varied between 0 and 14 species per plot. The five most frequent plant species were Echinochloa crus-galli (40% of the plots), Solanum nigrum (34%), Convolvulus arvensis (29%), Erigeron sumatrensis (28%), and Juncus bufonius (26%). The five most frequent insect species were Apis mellifera (13.6% of the plots), Altica sp. (9.1%), Sphaerophoria scripta (9.1%), Adelphocoris lineolatus (6.8%), and Adelphocoris sp. (6.8%).

Effects of agricultural management on plant and insect communities
The results of the GLMMs showed that management type (organic = 15.7 ± 2.7 species per plot, conventional = 5.4 ± 2.3 species per plot) (Fig. 2) and season (spring = 9.5 ± 4.9 species per plot, summer = 6.9 ± 6 species per plot) significantly affected only vascular plant species richness (Table 3). Regarding insect richness, the effect of management type was marginally significant (p = 0.06). The additive diversity partitioning showed significant contributions of α-diversity at the plot level to the total diversity of all the taxonomic groups, while the contribution of β-diversity at the plot level was significant for plants and total insects and marginally significant for pollinators. The contribution of β-diversity at the field level was never significant, while that of β-diversity for management types was significant only for plants (Fig. 3).
The results of PERMANOVA showed that management type was marginally significant only for pollinators (p = 0.07). Factors "Field" and "Season x Field" significantly affected the composition of the three analyzed groups. The effect of "Season" was never significant, as well as the interaction "Season x Management type" (Table 4).
The surveyed communities grouped in the NMDS ordination space only according to the sampling field, especially regarding plant communities. Nevertheless, some plant species like Echinochloa crus-galli and Convolvulus arvensis correlated with conventional fields. Regarding total insects and pollinators, we highlighted a correlation between several species (e.g., Adelphocoris lineolatus, Apis mellifera, Coccinella septempunctata, Sphaerophoria scripta) and the organic field "c", the one with the highest mean vegetation cover (Fig. 4).

Discussion
The observed differences between the surveyed communities were mostly attributable to the field identity. This is probably due to our study being limited by the very local scale and by the low number of studied fields, which are common issues when studying traditional, non-commercial crops (Latini et al. 2020). Likely, such a low number did not allow to distinguish the effect of management type from that of the environment, or of specific agricultural practices, on the variability of the surveyed communities. However, part of the achieved results is consistent with previous knowledge of the subject.
Our results showed inconsistent effects of conventional and organic management on plant, total insect, and pollinator communities, both in terms of species richness and composition. This is in accordance with previous evidence highlighting that the responses of biotic communities to organic and conventional management are different depending on the taxonomic group, with plants being the most affected organisms (Fuller et al., 2005;Tuck et al., 2014;Rundlöf et al., 2016). We found that the field was the only factor to significantly explain differences in both richness and composition for all the surveyed communities, suggesting that the variability between fields under the same management type is relevant and has a great influence on biotic communities. This strong identity of the fields can be related to several factors, representing the complex system of agricultural and  Table 3 Effects of management type and season on the species richness of plant, total insect, and pollinator communities according to the built GLMMs ↑ = richness increase; ↓ = richness decrease°0 .05 < p < 0.1, **p ≤ .01, ***p ≤ .001

Plants Total insects Pollinators
Factor Management (organic) ↑*** ↑°n.s Season (summer) ↓** n.s n.s Fig. 3 Additive diversity partitioning at the plot (α1 and β1), field (β2), and management type (β3) levels. Significance of differences is based on a random model with 99 permutations. Significance: ** p < 0.01, ° 0.05 < p < 0.1 environmental filters that shape biotic communities in arable land. For instance, landscape features were shown to highly affect the composition of plant and insect communities in agroecosystems, being even more important than farm management (Ponce et al., 2011). In the case of local, traditional crops, fields can be very small. This can lead to a high influence of surrounding vegetation on arable plant communities, which are colonized by plant species from adjacent habitats . E.g., we found in field e, which is located near a wood edge, plant species from such habitats like Quercus spp. and Xeranthemum cylindraceum. Time from conversion to organic management may also play a role for plant communities, with several years needed to observe shifts in species composition (Weibull & Östman, 2003). Plant communities were the only ones to show an increased species richness under organic practices, while the richness of total insects and pollinators was not affected by management type. This result is consistent with evidence from the literature. In fact, different taxonomic groups are known to respond differently to agricultural management, with the effects of organic farming being often dependent on single species or functional traits (Fuller et al., 2005). Most of the evidence shows that organic management increases plant species richness (Chamorro et al., 2016;Ponce et al., 2011;Ryan et al., 2010;Stein-Bachinger et al., 2021;Tyšer et al., 2021). However, also low-input, conventional farming is able to maintain a higher number of plant species in arable fields (Berbeć et al., 2020). Insects as well are highly affected by agricultural management, but they strongly respond at the same time to other types of factors. For instance, landscape structure can have relevant effects on insect communities, with simple landscapes reducing their species richness in organic fields (Winqvist et al., 2011(Winqvist et al., , 2012. For this reason, comparisons of entomological communities in organic and conventional farmlands should be done in similar landscape contexts (Masoni et al., 2017). Our fields were instead inserted in very different landscape types (intensive croplands, woodlands, peri-urban), which can be one reason for the marginally significant effect of Regarding species composition, only pollinator communities marginally responded to the different management type. The insignificant effect of management type on the composition of plant and total insect communities can be related to the high variability of agricultural practices carried out in different fields under the same management type, a variability that might not have been caught due to the low number of surveyed fields. Conversely, the field identity explained compositional differences in all the three groups of organisms, as well as the interaction between season and field, the latter indicating changes in species composition along the season for some fields, probably related to the different post-cultural management. This suggests once again that the variability of practices within management types is able to mask the effect of organic and conventional farming. In fact, in contrast with our findings, all the surveyed groups of organisms are known to be affected by organic and conventional management in terms of species composition (Andersson et al., 2013;Hyvönen et al., 2003).
The association of certain species to fields could be due to single practices, more than to the overall management type, as suggested by the NMDS results. For instance, fields with recurrent tillage practices and herbicide use resulted to host highly competitive plant species such as Echinochloa crus-galli, Amaranthus retroflexus, and Convolvulus arvensis (Heap, 2014), and to be negatively related to species that can be found also in semi-natural habitats like Cota tinctoria, Filago pyramidata, and Sherardia arvensis (Popic et al., 2013). Frequent mechanical disturbance also unfavoured ground-nesting bees like Adelphocoris lineolatus, Lasioglossum sp., and Halictus scabiosae, which are known to be related to fields with minimum tillage (Kratschmer et al., 2018). A high cover of arable vegetation was instead associated with phytophagous species and their predators, e.g., Adelphocoris spp., Coccinella septempunctata, Pieris rapae, and Sphaerophoria scripta, as already highlighted elsewhere (Sáenz-Romo et al., 2019).

Conclusions
In this work, we found out that the plant, insect, and pollinator communities of the investigated fields differ inconsistently in organic and conventional fields. Our initial hypothesis was partially confirmed since organic management resulted to significantly affect the richness and composition of plant communities only. Shifts in species richness and composition of the surveyed communities seemed to be rather explained by a complex set of other factors, both agronomic and environmental, which are related to the field identity. Several fields being under the same management type showed a variability of the carried-out practices, e.g., during the season, that led to a low affinity between them in terms of management. In order to be caught and interpreted, such variability should be described in a higher number of study sites. The limitations in extent of our research might have led to a partial identification of the differences, in relation to management type, between biotic communities in the investigated crop. Thus, more research is needed to verify possible further differences in richness and composition of the studied communities and their drivers.

Conflict interest
The authors have no competing interests to declare that are relevant to the content of this article.
Open Access This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http:// creat iveco mmons. org/ licen ses/ by/4. 0/.