An expert-assisted citizen science program involving agricultural high schools provides national patterns on bee species assemblages

Ecology studies often require large datasets. The benefits of citizen science for collecting such datasets include the extension of spatial and temporal scales, and cost reduction. In classical citizen science, citizens collect data and send them directly to scientists. This may not be possible for the many biological groups for which specimen identification is difficult and requires high-level expertise. Here we report the results of an expert-assisted citizen science program where teachers from 20 French agricultural high schools collected bees, which were identified to species level by a panel of expert bee taxonomists. Overall the dataset included 70 collections (year × sampling site combinations) that resulted in 4574 specimens belonging to 195 species. We analysed this dataset using data freely available at a national scale on agriculture intensity and landscape composition. We found that species richness increased with increasing proportion of herbaceous semi-natural elements; species dominance decreased with increasing crop diversity; the proportion of above ground nesting species and specimens increased as the intensity of agricultural practices decreased. Comparing the results obtained with identification to species level and those obtained with higher taxa or parataxonomic approaches, we found that the loss of taxonomic resolution resulted in the non-significance of some results on the effects of environmental variables on bee assemblage-level attributes. Our study suggests that identification to species level is of great importance to detect the effects of global change on bees and that an expert-assisted citizen science paradigm could provide relevant results to guide conservation measures at a national scale.


Introduction
Data on species distribution provide baseline information in biogeography and population trend studies (Kerr et al. 2007;Cardoso et al. 2011).In the context of global change, these data have gained importance in the assessment of the effects of anthropogenic disturbances on biodiversity (Butchart et al. 2010).Such studies often require large datasets that cannot be gathered by individual research teams (Devictor et al. 2010;Dickinson et al. 2010).For centuries, amateur naturalists have collected data on species occurrence (Miller-Rushing et al. 2012;Pocock et al. 2015).I n recent years, Abstract Ecology studies often require large datasets.The benefits of citizen science for collecting such datasets include the extension of spatial and temporal scales, and cost reduction.I n classical citizen science, citizens collect data and send them directly to scientists.This may not be possible for the many biological groups for which specimen identification is difficult and requires highlevel expertise.Here we report the results of an expertassisted citizen science program where teachers from 20 French agricultural high schools collected bees, which were identified to species level by a panel of expert bee taxonomists.Overall the dataset included 70 collections (year × sampling site combinations) that resulted in 4574 specimens belonging to 195 species.We analysed this citizens collect data and send them directly to scientists.In the expert-assi sted ci ti zen sci ence paradi gm, citizens provide material to experts that analyse it (e.g., identification of specimens to species level) and these experts send the data to scientists (Fig. 1).We here report the first results of an ongoing citizen science program where a specific group of participants, namely teachers from agricultural high schools in France, conducted standardized surveys of bees, which were then identified to species level by expert bee taxonomists.This approach combines the benefits of data collection by a national network of volunteers and those of data reliability to species level provided by the contribution of experts.We first study the effects of agriculture intensity and landscape composition on bee species assemblages.Then we compare the results obtained with the species level (i.e. with the help of experts) to those obtained with higher taxa or parataxonomic approaches (i.e.without the help of experts and participants being able to provide data directly to scientists).We hypothesise that (1) bee taxonomic richness will decrease, and inversely dominance (the proportion of the most abundant taxon) will increase, with increasing agriculture intensity and decreasing amount of semi-natural elements; (2) increasing agriculture intensity and decreasing amount of semi-natural elements will influence the functional composition of bee assemblages, with parasitic, oligolectic, and above ground nesting species being negatively affected; (3) the loss of the species level resolution will result in the nonsignificance of some results on the effects of environmental variables on bee assemblage-level attributes.

Study sites
From 2009 onwards, volunteer biology and ecology teachers from French agricultural high schools were asked to collect wild bees following a standardized protocol.Contributions came from most continental French areas with land cover dominated by agricultural land-use (Fig. 2).The 20 contributing schools were located in rural areas and comprised a farm for educational purposes.They encompassed the diversity of situations currently occurring in French agriculture, in terms of production types (farms devoted to annual or perennial crops, to livestock, or mixed farming systems), agriculture intensity, and relative surface of semi-natural elements in the landscape (see "Environmental variables").

Bee sampling and identification
Netting, pan traps and trap nests are common methods used to sample bees (Westphal et al. 2008).I n this study, citizen science programs have flourished in ecology, with many benefits such as an extension of spatial and temporal scales (including data collection on private areas such as gardens), time and cost reduction, and reconnection of people with nature leading to an increasing public awareness about environmental issues (Bonney et al. 2009;Devictor et al. 2010;Dickinson et al. 2010;Birkin and Goulson 2015;van der Wal et al. 2015).
Birds and butterflies dominate as study groups in citizen science (Schmeller et al. 2009;Devictor et al. 2010;Dickinson et al. 2010;Roy et al. 2015).I n these groups, comprehensive identification tools exist and identification to species level is relatively easy.For many other biological groups, identification to species level is virtually impossible without collecting specimens and identifying them using both reference collections and relevant literature.I n those cases, the alternatives to identification to species level by citizen scientists are the identification at higher taxonomic (e.g., genus, family) or at parataxonomic levels (sensu Krell 2004) sometimes based on photographs (Deguines et al. 2012;Casanovas et al. 2014) or the help of experts to reach the species level (Gardiner et al. 2012 for the Coccinellidae).
Bees (Hymenoptera: Apoidea) constitute one of these species-rich and diverse groups (926 species recorded in continental France, Kuhlmann et al. 2015) for which identification to species level is inherently challenging.Moreover, limited and sparse identification tools are available (Brown and Paxton 2009;Patiny et al. 2009).The identification to species level is therefore dependent on the help of a panel of expert bee taxonomists, i.e. persons recognised as an authority in their field, typically academics or museum scientists but sometimes also non-academics with a deep interest in the natural history of a particular taxon.Bees have become the focus of much interest over the last decade because of concerns about species decline and the expected consequences on wild plant and crop pollination (Steffan-Dewenter et al. 2005;Biesmeijer et al. 2006;González-Varo et al. 2013;Ollerton et al. 2014;Goulson et al. 2015;Garibaldi et al. 2016).Although long-term recording allowed the assessment of species trends in some countries (see e.g.Biesmeijer et al. 2006;Ollerton et al. 2014;Thomas et al. 2015 for British studies based on bee records collected by the UK BWARS-Bees Wasps & Ants Recording Society), data are severely lacking regarding bee species distribution (Brown and Paxton 2009;Patiny et al. 2009;Nieto et al. 2014) and community composition (Winfree et al. 2011) to understand the effects of global change on bee populations and thereby be able to design meaningful conservation strategies (Cardoso et al. 2011).
Our overall objective was to test an expert-assisted citizen science paradigm to study the effects of global change on bee species assemblages in agricultural landscapes over a large area.I n the class i cal c i t i zen sc i ence parad i gm, 1 3 (Institut National de la Recherche Agronomique), Avignon, France, the scientific coordinator of the program.Participants chose the number of sampling sites (ranging from one to seven per school, mean ± S.E.2.25 ± 0.35) and their location in the school.Participants were asked to make, under suitable weather conditions (minimum of 15 °C, low wind, no rain), one 24-h capture session each month from March to October, in order to obtain an overall assessment of the bee assemblage composition over the whole flying season (Banaszak et al. 2014).In practice, some participants could not carry out the required number of sessions or others made more than one session in a month.In order to derive meaningful assemblage summary statistics, we discarded the four annual datasets made up of <5 sampling months.This resulted in 70 sampling site × year combinations, 6 in 2009, 24 in 2010, and 40 in 2011 (Online Resource 1).For the sake of simplicity, we further defined a collection as the captures made at a given sampling site during a given year.The number of sampling dates (that is 24-h sessions) for a collection ranged from 5 to 12. Participants were asked to keep the same number and the same location for sampling sites across years, but in practice these parameters changed somewhat as some locations became inaccessible, or enthusiastic participants wanted to increase their sampling effort after their first year in the project.We totalled 45 different sampling sites over the 20 schools, which were sampled we chose pan trapping as a standardized sampling method because: (1) this method is particularly well suited to minimize collector biases (Westphal et al. 2008).These biases could have been considerable with netting in our case given the high number of participants, most of which had no prior entomological experience; (2) pan traps can provide reliable data to assess the overall species richness of a study site even though they are known to be poor in capturing species of genera like Bombus and Colletes (Roulston et al. 2007;Westphal et al. 2008;Wilson et al. 2008, but see;Wood et al. 2015); (3) pan traps are cheap and simple to use (LeBuhn et al. 2013).Traps were made of 500-ml plastic bowls that were sprayed inside with a UV-reflecting paint and were mounted on a wooden pole.One disadvantage of pan traps is that their effectiveness is affected by the local environmental context, in particular the local floral resource availability (Westphal et al. 2008;Wilson et al. 2008;Baum and Wallen 2011).To cope with this bias, traps were placed near flowers when there were flowers at the site (flowers were not necessarily present all year long), at a height slightly above that of the average vegetation, and in the sun inasmuch as was possible so as to be clearly visible.A sampli ng si te consisted of a set of three pan traps (a blue, a white, and a yellow) and sampling was initiated by filling each bowl with 400 ml of water with a drop of detergent.I dentical pan traps were supplied to all schools by I NRA In the classical citizen science paradigm, citizens collect data that are directly sent to scientists.Classical citizen science requires various skill levels but the identification of biological material at a useful level remains doable with common tools (e.g., books, identification keys).In the expertassisted citizen science paradigm, citizens provide the material (e.g., specimens, photographs) to experts that analyse it (e.g., identification of specimens to species level) and these experts send the data to scientists.The involvement of experts is required when identification for the biological group under study is difficult, requires high level of knowledge and experience, and specific tools (e.g., reference collection).Grey elements illustrate the expert-assisted targeted citizen science program described in this paper 1 3 Kuhlmann et al. (2015).Honey bees (Apis mellifera) were caught, but this species was not considered as its abundance is largely determined by beekeeping.Therefore we use "bees" synonymously with "wild bees" in the following.

Assemblage-level attributes
For each collection, we calculated two sets of attributes.The first set required the identification of bee specimens to species level, i.e. with the help of expert bee taxonomists (expert-assisted citizen science paradigm).The attributes of the second set did not require the identification to species level and could be obtained by teachers alone, provided that they received the course on bee identification to genus level (classical citizen science paradigm).over 1, 2, or the 3 years (2009, 2010, and 2011).These 45 sampling sites were located in 25 different municipalities (local administrative units).
At the onset of the program, all the participants followed a 5-day course on bee biology and systematics provided by scientists and bee experts, and which included training on techniques to prepare specimens recovered from pan traps and identification to genus level.Afterwards, all bees caught during the program were pinned, labelled and pre-identified to genus by the teachers, centralized and double-checked at INRA Avignon, and then sent to a panel of expert bee taxonomists to be identified to species level.Specimens were deposited in the collection of I NRA Avignon and, for the most common species, in a reference collection in each participating school.Taxonomy followed the nomenclature of additional subdivision reduces the taxonomic resolution that is lost when using genera instead of species (van Rijn et al. 2015).The inter-tegular distance (ITD, the distance between the two insertion points of the wings) is often used as a proxy for body size in bees (Cane 1987).ITD values for our dataset were obtained from Fortel et al. (2014).Each genus was subdivided into three size classes: ITD <2.6, 2.6-4.6, and >4.6 mm.We defined the parataxonomic richness as the number of parataxonomic units caught over 1 year at the sampling site; -Parataxonomic richness, with seven subdivisions according to body size: This is the same approach as above but instead of three, seven subdivisions were considered following van Rijn et al. ( 2015): ITD <1, 1-1.5, 1.5-2, 2-2.5, 2.5-3, 3-3.5, and >3.5 mm; -Genus dominance: This is the proportion of specimens for the most abundant genus; -Parataxonomic dominance, with three subdivisions according to body size: This is the proportion of the most abundant parataxonomic unit when considering three ITD subdivisions; -Parataxonomic dominance, with seven subdivisions according to body size: This is the proportion of the most abundant parataxonomic unit when considering seven ITD subdivisions.

Environmental variables
We related our ecological dataset to environmental data freely available at the national scale.We used information derived from the High Nature Value (HNV) indicator regarding agriculture intensity and the Corine Land Cover (CLC) database regarding the landscape composition.

Agriculture intensity at the municipality level
We assigned each of the 25 municipalities to four variables regarding agriculture intensity: "Crop diversity" (a proxy for the crop rotation system), "Extensive farming practices" (an estimation of pesticide, mineral fertilizers, and irrigation use), "Landscape elements" (an estimation of the relative area of semi-natural elements), and an overall index of agriculture intensity.The three first variables were scored from 0 (low crop diversity, high input level, and low availability of semi-natural elements, respectively) to 10.The overall intensity index was obtained in summing the three scores (with low index values corresponding to high agriculture intensity).These four variables were derived from the French HNV dataset.European Union countries have been required to identify HNV farmlands, i.e. areas that include semi-natural elements, low-intensity farming and diverse, First set of attributes (expert-assisted citizen science paradigm) -Species richness: This is the number of species caught over 1 year at the sampling site; -Species dominance: This is the proportion of the most abundant species (also known as Berger-Parker index); -Ecological traits: Each species was described according to three ecological traits that have been shown to be important to determine the response of bees to environmental disturbances (e.g., Moretti et al. 2009;Williams et al. 2010;Winfree et al. 2011;Sheffield et al. 2013;Rader et al. 2014; Kremen and M'Gonigle 2015): reproducti ve strategy, trophi c speci ali zati on, and nesting behaviour.Information was compiled from the literature (Westrich 1989;Amiet et al. 1999Amiet et al. , 2001Amiet et al. , 2004Amiet et al. , 2007Amiet et al. , 2010;;Michener 2007).Occasionally, ecological traits can be inferred from the genus (e.g., all Nomada species are parasitic), but the species level is usually necessary to determine ecological traits.For examples, very similar (cryptic) species can have opposite behaviour (for example very similar Bombus species are either parasitic or non-parasitic) or a single genus may comprise both oligolectic and polylectic species (e.g., the genus Andrena).
For the reproductive strategy, we separated nonparasitic from parasitic species and the subsequent classes were made only for non-parasitic species.For trophic specialisation regarding pollen use, we distinguished between oligolectic species collecting pollen on plant species from a single family and polylectic species collecting pollen on several plant families.For nesting behaviour, species nesting below ground were separated from those nesting above ground in cavities (e.g., plant stems, wood holes, snail shells).Then, we calculated the proportion of species and the proportion of specimens for modalities that have been shown to be negatively affected by environmental disturbances, that is "parasitic", "oligolectic", and "above ground nesting" modalities (Biesmeijer et al. 2006;Williams et al. 2010;Rader et al. 2014).

Second set of attributes (classical citizen science paradigm)
-Genus richness: This is the number of genera caught over 1 year at the sampling site; -Parataxonomic richness, with three subdivisions according to body size: We used this attribute to test whether the "classical citizen science" approach could be improved by further subdividing genera using body size.Doing this, we obtained parataxonomic units, i.e. based on external morphology (sensu Krell 2004).This inspected against fitted values to ensure residual normality and homoscedasticity requirements were fulfilled.Models were fitted using the maximum log-likelihood criteria, and the significance of environmental variables was assessed a posteriori using likelihood ratio deletion tests.When necessary, unequal sampling effort (number of sampling dates in an annual survey) was accounted for by specifying sampling effort as an offset term.This was required for analyses of species, genus and parataxonomic richness, but deemed unnecessary for dominance attributes and ecological trait proportions.Analyses were conducted using R software version 3.1.1(R Core Team 2014), using the package lme4 (Bates et al. 2014).

Overall taxonomic and functional composition of the bee dataset
Participants collected a total of 4574 specimens representing 195 species (Online Resources 1, 2).The family Halictidae largely dominated the captures (73.3 % of specimens and 31.3 % of species), followed by Andrenidae (16.2 % of specimens and 23.6 % of species) and Apidae (6.8 % of specimens and 22.6 % of species).Colletidae, Megachilidae and Melittidae all represented <2 % of the specimens and 6.1, 15.4, and 1.0 % of the species, respectively.The genus Lasioglossum was the most abundant, representing 55.4 % of the specimens.The genus Andrena represented 15.6 % of the specimens, but it was the most species rich with 43 species (vs.36 species of Lasioglossum).Lasioglossum malachurum was a superabundant species, representing 20.0 % of all specimens.I t was detected in 88.6 % of collections and dominated in 45.7 %.The second most abundant species, L. morio, represented 5.7 % of specimens, was detected in 67.1 % of collections and dominated in 11.4 %.These two species were followed by four species of Halictidae (Lasioglossum glabri usculum, L. pauxi llum, Hali ctus scabi osae, H. tumulorum), each making between 4 and 5 % of the specimens.The following species, i.e. the most abundant nonhalictid species, was Andrena flavipes that made-up 3.1 % of specimens and was present in 52.9 % of collections.
Regarding the ecological traits, 30 parasitic and 165 nonparasitic species were caught.These non-parasitic species included 47 above ground versus 116 below ground nesting species (information was missing for two species), and 36 oligolectic versus 120 polylectic species (with diet information lacking for nine species).All seven most abundant species (species that accounted for more than 3 % of the total abundance) were non-parasitic, below ground nesting, polylectic, and social bees (with the exception of H. scabiosae and A. flavipes that are solitary).
small-scale mosaics of land-use types.The details of how to implement this indicator were at the discretion of each state (CEC 2006).In France, indicators aggregating statistics of agricultural holdings were calculated at the municipality level (Pointereau et al. 2010;Doxa et al. 2012;Deguines et al. 2014).The national mean for the overall intensity index was 12.20 (S.E.±1.8 × 10 −4 ) and its values ranged from 1 to 30 when considering all 36,027 French municipalities.In our dataset, the average value was nearly the same as the national average (mean ± S.E.12.42 ± 0.77) and the values ranged from 4.92 to 25.12.

Proportion of semi-natural elements at the landscape scale
To examine the effect of landscape composition, we obtained land cover data from the CLC 2006 database (Bossard et al. 2006).CLC is a geo-referenced database including the main habitats for the European Union countries in contiguous polygons classified according to 44 different landcover categories in the finest classification.We quantified the landscape composition in 100 and 500 m radius windows centred on sampling sites by using the geographical information system package ArcGI S 10.1 (ESRI 2012).We chose these window sizes as relevant scales for flight and foraging distances in bees (Gathmann and Tscharntke 2002).For each sampling site, we calculated the proportion of semi-natural herbaceous elements.In our case, this variable included the CLC classes "Pastures", "Complex cultivation patterns", "Land principally occupied by agriculture, with significant areas of natural vegetation", "Moors and heathland", and "Transitional woodland-shrub".The proportion of herbaceous semi-natural elements ranged from 0 to 100 % (mean ± S.E.37.7 ± 5.3) in 100 m radius windows and from 0 to 99.05 % (mean ± S.E.39.9 ± 3.3) in 500 m radius windows.We also considered a variable including both herbaceous and woody elements (forests), but models returned non-significant results (not shown) and so they were not considered further.

Statistical analyses
The assessment of assemblage-level responses to environmental variables was performed within a generalized linear mixed model (GLMM) framework in order to cope with the non-independence of surveys carried out by the same contributing school at different sampling sites and repeated over several years.We therefore parameterized GLMMs with assemblage-level attributes as response variables, environmental variables as fixed effect candidate covariates and school, municipality and sampling site identities as a suite of hierarchically nested random grouping variables.Normality requirements of continuous data were tested using the Shapiro-Wilk statistics.Model residuals were further 1 3 pesticide use decreased.Similarly, the proportion of the most abundant parataxonomic unit, when considering three subdivisions according to body size, decreased with decreasing agriculture intensity and with increasing crop diversity (Table 2).

Discussion
We present the results of a monitoring program conducted at the national scale and which involved a targeted group of citizens (teachers from agricultural high schools), researchers, and expert bee taxonomists (Fig. 1).Our dataset provided general patterns regarding functional and taxonomic composition of bee assemblages, as well as responses of bees to agriculture intensity and landscape composition.Comparing the results obtained through expert-assisted and classical citizen science paradigms, i.e. the results with identification to species level and those obtained with higher taxa or parataxonomic approaches, we found that the loss of taxonomic resolution resulted in the non-significance of some results on the effects of environmental variables on bee assemblage-level attributes.

Taxonomic and functional composition of assemblages
Over 3 years of sampling in 20 agricultural high schools and 4574 individuals caught, this program provided information on the distribution of 195 species, representing 21.1 % of the 926 species recorded in continental France (Kuhlmann et al. 2015).The dataset did not include threatened species at the European level (Nieto et al. 2014), but some species were uncommon at the national scale such as Andrena apicata, A. hattorfiana, and A. ventricosa (David Genoud, pers. comm.).Although the 45 sampling sites encompassed a great diversity of climate, agricultural, landscape and local scale contexts, we found some recurrent features across sites and years.The dominance of Halictidae, Lasioglossum sp. in particular, appeared as a general trend.This feature has been observed in various places across the world [e.g., Marini et al. (2012) in Italy; Morandin and Kremen (2013) in USA; Fortel et al. (2014) in France; Rader et al. (2014) in New-Zealand; Saunders and Luck (2014) in Australia; Pisanty and Mandelik (2015) in Israel; Le Féon et al. (2016) in Argentina].These species are especially well caught by pan traps but their high abundance is also observed when bees are sampled by netting (e.g., Rollin et al. 2015).The dominance of Lasioglossum species in bee assemblages may be due to their eusociality (in some species) and their polylecty.Beyond these results at the family and genus levels, our results provide insights on the status of species at the national scale.In particular, Lasioglossum malachurum was virtually omnipresent.Such information on the identity of Regarding richness attributes that could be obtained by the teachers without the help of expert bee taxonomists (classical citizen science paradigm), the dataset comprised 25 genera, 31 parataxonomic units when considering three subdivisions according to body size, and 49 parataxonomic units when considering seven subdivisions.
The dominant species made-up between 11.5 and 71.1 % of the total abundance (mean ± S.E.30.8 ± 1.5), and it was always a non-parasitic and below ground nesting species.In most cases, it was social (54 collections) and polylectic (65 collections).Oligolectic species were dominant in some rare cases (Dasypoda hirtipes and Tetralonia malvae both in two collections).L. malachurum dominated in 32 collections, and L. morio, H. scabiosae and H. tumulorum dominated in 8, 5 and 5 collections, respectively.

Comparison between the expert-assisted and the classical citizen science paradigms when relating bee assemblage-level attributes and environmental variables
Species richness increased with the increasing proportion of herbaceous semi-natural elements in 100 m radius windows (Table 1; Fig. 3a).Species dominance decreased with increasing crop diversity (Fig. 3b).The proportion of above ground nesting species and specimens increased with decreasing agriculture intensity (significant only for the proportion of specimens), with increasing crop diversity, and when the intensity of fertilizer and pesticide use decreased (Fig. 3c).There was no significant relationship between attributes regarding either parasitic status or trophic specialization with environmental variables (Table 1).
Considering the attributes obtained without the help of expert bee taxonomists (classical citizen science paradigm), there was no significant relationship between genus and parataxonomic richness and environmental variables (Table 2; Fig. 4).The genus dominance decreased with decreasing agriculture intensity, with increasing crop diversity (Fig. 3d) and when the intensity of fertilizer and Our study suggests that identification to species level is of great importance for the detection of the effect of environmental variables on bee richness.Furthermore, as many ecological traits are species-specific rather than genus-or family-specific, comprehensive functional approaches cannot be reached without species-level identification (Williams et al. 2010).

Relationships with agriculture intensity and landscape composition
When ecological data collection occurs at a local scale, information on landscape or agricultural practices may be obtained simultaneously, for example by site mapping from aerial photographs associated to ground-checking and farmer surveys, respectively.For studies at the regional or national scales, this may not be possible and one must seek for environmental data already available (see also Deguines et al. 2012Deguines et al. , 2014)).
I n France, information on agriculture intensity derived from HNV is given at the municipality level (i.e. 15 km 2 on average), which may not be the most relevant scale when considering bee assemblages.Similarly, CLC data are subject to several drawbacks, notably poor spatial definition (the minimum area for an element to be represented is 25 hectares), annual crops pooled together without information on cultivated species (e.g., no information on mass-flowering crops, that can play a major role on bee assemblages, Diekötter et al. 2010;Le Féon et al. 2013).
Despite these major drawbacks, analyses of our empirical dataset gave significant results that confirmed at a national scale those of previous studies.The decrease in species richness associated with decreasing amount of semi-natural elements has been mentioned in several studies (see Winfree et al. 2011 for a review).The particular sensitivity of above ground nesting species to environmental disturbances has also been reported previously (see Williams et al. 2010 for common versus uncommon species could be trivial regarding well-known groups such as butterflies but remain useful given the lack of knowledge on the distribution and status of bees in France. From a functional point of view, non-parasitic, polylectic, and below ground nesting species dominated.However parasitic, oligolectic and above ground nesting species each represented a substantial proportion of species and specimens suggesting that farmland on agricultural schools may harbour populations of species with more specific ecological requirements that are also more sensitive to environmental disturbances (Rader et al. 2014).

Comparison between expert-assisted and classical citizen science paradigms
When comparing the results obtained with identification to species level to the results obtained with higher taxa (genus) or parataxonomic (genus associated with body size subdivisions) approaches, we found that the effect of crop diversity on dominance remained significant when replacing the species level by alternatives that did not require the help of the expert bee taxonomists.On the other hand, the effect of landscape composition on bee richness was not detected when losing the species level resolution.
Previous citizen science programs on bees have provided interesting results on the broad effects of land-use changes (Deguines et al. 2012), on the role of hedgerows for bees (Kremen et al. 2011), and on the nesting ecology and the status of bumblebee species in the UK (Lye et al. 2012).These programs were based on photographic collections (Deguines et al. 2012) and on field observations by trained citizens (Kremen et al. 2011) or members of the public (Lye et al. 2012) and provided information at higher taxonomic levels (genus, family) or at morphospecies or species group levels. 1 3 This aspect has been little investigated in bees (but see Sheffield et al. 2013;Marini et al. 2014).Thus our results that dominance at the species, genus and parataxonomic levels decreased with crop diversity provides new insights on the influence of agricultural practices on bee assemblages.I ndeed, it provides support to the new CAP policy that stresses the importance of crop diversity in order for farmers to be eligible for EU subsidy (European Union 2013).

Conservation implications
Getting biodiversity data at large spatial and temporal scales is the first stage of any conservation plan (Cardoso a review).Our result may be due to the direct vulnerability of these species to pesticide use (Vaughan et al. 2014) or to the hidden effects associated with agriculture intensification such as the loss of micro-habitats suitable for nest establishment (Kremen and M'Gonigle 2015).
The presence of a few very abundant and many rare species in biotic assemblages is an universal law in ecology (McGill et al. 2007) and, indeed, it has been often noticed in bees (e.g., Williams et al. 2001).However, in various groups, the degree of dominance of the most abundant taxa, as well as their identity, respond to environmental changes, sometimes more rapidly than does species richness (Caruso et al. 2007;Hillebrand et al. 2008;Tolkkinen et al. 2013).programs have been proven to be very useful tools for such purposes (Devictor et al. 2010;Dickinson et al. 2010).Regarding bees, we showed that citizen science programs combining the collection of biological material by a targeted network of citizens and the identification to species level et al. 2011; Goulson et al. 2015).The national scale often has the advantage of representing both a large scale from an ecological point of view and a policy-relevant scale for the future implementation of conservation measures (Woodcock et al. 2014;Budge et al. 2015).Citizen science ).Models contained an assemblage-level attribute as the response variable, an environmental variable as fixed effects, and school, municipality and sampling site identities as a suite of hierarchically nested random effects.In each case, the first four environmental variables were derived from the High Nature Value data and calculated at the municipality level.The proportion of herbaceous semi-natural elements was derived from Corine Land Cover data ( 2006) and calculated at two spatial scales (100 and 500 m radius centred on sampling sites) 1 3 (often crucial for most conservation issues) by expert bee taxonomists could provide relevant results to guide conservation measures at a national scale.Moreover, through the involvement of agricultural high schools, our study focused on agricultural landscapes where bee conservation has been proven to be essential because of the large areas they occupy and the benefits of bee pollination on entomophilous crops (Deguines et al. 2012).

Fig. 1
Fig.1Conceptual framework (adapted fromDevictor et al. 2010) of the expert-assisted citizen science paradigm (in black).In the classical citizen science paradigm, citizens collect data that are directly sent to scientists.Classical citizen science requires various skill levels but the identification of biological material at a useful level remains doable with common tools (e.g., books, identification keys).In the expertassisted citizen science paradigm, citizens provide the material (e.g.,

Fig. 2
Fig. 2 Location of the 20 agricultural high schools involved in the study in continental France.The number of collections (sampling site × year combinations) is given for each school (n = 70 collections in total)

Fig. 3
Fig.3Relationship between (a) bee species richness and the proportion of herbaceous semi-natural elements in 100 m radius windows; b the proportion of the most abundant species (species dominance) and the index of crop diversity at the municipality level; c the proportion of above ground nesting bee species and the index of extensive farming practices; d the proportion of the most abundant genus (genus dominance) and the index of crop diversity at the municipality level.

Fig. 4
Fig. 4 Number of units for species richness and for alternative (para) taxonomic attributes (black bars), and P-value for the generalized linear mixed model including the assemblage-level attribute and the proportion of herbaceous semi-natural elements in 100 m radius windows as environmental variable (see Tables 1, 2 for details).The white bar denotes a significant result (P ≤ 0.05) whilst the grey bars denote nonsignificant results.'Parataxonomic richness (7)' and 'Parataxonomic richness (3)' are richness measures based on either 7 or 3 subdivisions of the inter-tegular distance (which is used as a surrogate for body size) = 70 collections and df = 1 in all cases.Significant P-values are indicated in bold (*P ≤ 0.05; **P ≤ 0.01).Models contained an assemblage-level attribute as the response variable, an environmental variable as fixed effects, and school, municipality and sampling site identities as a suite of hierarchically nested random effects.In each case, the first four environmental variables were derived from the High Nature Value (HNV) data and calculated at the municipality level.The proportion of herbaceous semi-natural elements was derived from Corine Land Cover data (2006) and calculated at two spatial scales (100 and 500 m radius centred on sampling sites)

Table 2
Results of the generalized linear mixed models for the assemblage-level attributes based on the identification of bee specimens at genus level or at parataxonomic levels (i.e. that did not require the help of expert bee taxonomists) = 70 collections and df = 1 in all cases.Significant P-values are indicated in bold (*P ≤ 0.05; **P ≤ 0.01