# A ‘Landscape physiology’ approach for assessing bee health highlights the benefits of floral landscape enrichment and semi-natural habitats

Cédric Alaux, Fabrice Allier, Axel Decourtye, Jean Francois Odoux, Thierry Tamic, Mélanie Chabirand, Estelle Delestra, Florent Decugis, Yves Le Conte,

Mickaël Henry

► To cite this version:

Cédric Alaux, Fabrice Allier, Axel Decourtye, Jean Francois Odoux, Thierry Tamic, et al. A ‘Landscape physiology’ approach for assessing bee health highlights the benefits of floral landscape enrichment and semi-natural habitats. Scientific Reports, 2017, 7, 10 p. 10.1038/srep40568. hal-01606556

HAL Id: hal-01606556

https://hal.science/hal-01606556

Submitted on 26 May 2020

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
SCIENTIFIC REPORTS

<figure><img src="image_2.png" /></figure>

OPEN

# A ‘Landscape physiology’ approach for assessing bee health highlights the benefits of floral landscape enrichment and semi-natural habitats

Received: 28 April 2016

Accepted: 08 December 2016

Published: 13 January 2017

Cédric Alaux<fnref n="1" />,², Fabrice Allier²,³, Axel Decourtye²,³,⁴, Jean-François Odoux⁵, Thierry Tamic⁵, Mélanie Chabirand⁵, Estelle Delestra⁶, Florent Decugis<fnref n="1" />, Yves Le Conte<fnref n="1" />,² & Mickaël Henry<fnref n="1" />,²

Understanding how anthropogenic landscape alteration affects populations of ecologically- and economically-important insect pollinators has never been more pressing. In this context, the assessment of landscape quality typically relies on spatial distribution studies, but, whether habitat-restoration techniques actually improve the health of targeted pollinator populations remains obscure. This gap could be filled by a comprehensive understanding of how gradients of landscape quality influence pollinator physiology. We therefore used this approach for honey bees (Apis mellifera) to test whether landscape patterns can shape bee health. We focused on the pre-wintering period since abnormally high winter colony losses have often been observed. By exposing colonies to different landscapes, enriched in melliferous catch crops and surrounded by semi-natural habitats, we found that bee physiology (i.e. fat body mass and level of vitellogenin) was significantly improved by the presence of flowering catch crops. Catch crop presence was associated with a significant increase in pollen diet diversity. The influence of semi-natural habitats on bee health was even stronger. Vitellogenin level was in turn significantly linked to higher overwintering survival. Therefore, our experimental study, combining landscape ecology and bee physiology, offers an exciting proof-of-concept for directly identifying stressful or suitable landscapes and promoting efficient pollinator conservation.

Anthropogenic effects on landscape (habitat loss, fragmentation and degradation) expose most insect pollinators to new and enduring environmental challenges and are primary drivers of their decline<fnref n="1" />⁻⁵. This represents a major conservation issue because insect pollination is vitally important to the maintenance of biodiversity and crop production⁶,⁷. Therefore, there is an urgent need to understand how landscape alteration affects those populations, and to promote landscape restoration, notably via agri-environment schemes (incentives for farmer to benefit the environment)⁸⁻¹³.

Traditionally, studies have focused on the relation between species distribution (e.g. presence/absence, abundance) and landscape patterns<fnref n="1" />,³,⁶,¹⁴⁻¹⁹. However, while informative, the assessment of disturbances is limited because the health state of the population is not considered and the deleterious effects of landscape alteration can only be detected once the population has started to decline²⁰,²¹. A more powerful approach would be to characterize the specific mechanisms underlying the population response by combining physiological and ecological knowledge²¹. Indeed, the persistence of a population can be inferred by the health conditions of individuals within the population, and their physiological responses to environmental changes can provide an early indication of a stressful landscape²⁰,²². But most importantly it provides a cause-and-effect relationship between

<footnote n="1">¹INRA, UR406 Abeilles et Environnement, Domaine Saint-Paul, CS 40509, 84914 Avignon, France. ²UMT PrADE, CS 40509, 84914 Avignon, France. ³ITSAP-Institut de l'Abeille, Domaine Saint-Paul, CS 40509, 84914 Avignon, France. ⁴ACTA, CS 40509, 84914 Avignon, France. ⁵INRA, UE1255 Entomologie, 17700 Surgères, France. ⁶ADAPIC, Cité de l'Agriculture, 45921 Orléans, France. Correspondence and requests for materials should be addressed to C.A. (email: cedric.alaux@inra.fr)</footnote>
landscape quality and population response, which has the potential to directly contribute to decision-making and support conservation policy²⁰⁻²².

Floral resource availability in different landscape contexts has been linked to colony growth or productivity and variation in nutritional variables in both bumble bees²³⁻²⁵ and honey bees²⁶⁻³⁰. However, knowledge on the connection between landscape quality, notably landscape enrichment with floral resources, and bee health, is clearly limited. We therefore used a 'Landscape physiology' approach, integrating physiological data with landscape ecology²⁰, to test i) the connection between bee health and landscape quality, and ii) whether agri-environment schemes can provide benefits to bee health. For that purpose, we exposed honey bee colonies (*Apis mellifera*) to different agricultural landscapes, either enriched or not by melliferous catch crops (environmentally friendly practices to promote bee forage) and surrounded by semi-natural habitats. We then assessed the link between the landscape quality (catch crop and semi-natural habitats), bee physiology and the consequential colony survival.

The study was performed during the pre-wintering period because severe winter mortality recently observed in honey bee colonies³¹,³² suggests that preparation for overwintering is especially challenging. Indeed, sufficient energetic reserves must be stored at the individual and colony level for a successful overwintering³³. Bee health was assessed by determining fat body mass and the gene expression level of vitellogenin²⁸,³⁴. Both are physiological features of winter bees that arise during the autumn in temperate regions as an adaptation for surviving throughout the winter period. Indeed, winter bees have greater nutrient storage in the fat body and tolerance to oxidative stress than summer bees, due to the storage protein vitellogenin³⁵,³⁶. This ubiquitous protein, produced in the fat body, acts as an antioxidant and promotes the longevity of bees³⁷. Its level is high in young bees but exhibits a negligible decline over time in winter bees as compared to summer bees³⁸, likely explaining why winter bees are long-lived (several months) as compared to summer bees (4–6 weeks). Since fat body growth and vitellogenin production are both triggered by pollen intake³⁹,⁴⁰, we hypothesized a connection between their levels and the landscape-wide floral resource availability. In addition, we assessed the infestation levels of *Varroa destructor* as this parasitic mite is known to have detrimental effects on overwintering survival³⁵.

## Results and Discussion

### Ecophysiological basis of bee health.

Before wintering, colonies were set up either inside (n = 184 colonies split into 10 apiaries) or outside a melliferous catch crop area (n = 166 colonies split into 8 apiaries), following a control-vs-treatment experimental design (supplementary Figs S1 and S2). The paired control-vs-treatment experimental set-up was designed to avoid concomitant variations of semi-natural habitat land cover while varying catch crop treatment. Catch crop and semi-natural land cover values were therefore maintained uncorrelated (n = 18, Pearson's correlation r = -0.17, P = 0.49).

We first explored the dataset for possible confounding effects due to year (winters 2012–13 and 2013–14) or colony allocation between treatments. No inter-annual variation in overwintering survival was detected when considering year effect alone (generalized linear mixed models (GLMM), df = 347, z = 1.19, P = 0.23) nor in combination with the fully parameterized survival model including physiological, brood and landscape covariates (df = 165, z = -0.426, P = 0.67). We could therefore consider the apiaries from different years as independent replicates within each beekeeping set-up. Furthermore, the random allocation of colonies led to apiaries with brood initial state independent from their landscape context, either considering colonies as independent entities (simple linear model, df = 345, catch crop land cover: t = 0.88, P = 0.38; semi-natural habitat land cover: t = 1.17, P = 0.24), or specifying a random grouping structure to properly account for the non-independencies within beekeeping set-ups (GLMM, df = 343, catch crop land cover: t = -1.03, P = 0.30; semi-natural habitat land cover: t = 1.04, P = 0.30).

We then performed a path analysis to disentangle the direct and indirect dependencies of bee physiological traits and overwintering survival (n = 350 colonies) on landscape quality (gradients of catch crop from 0 to 0.315 km² and semi-natural habitats from 0.04 to 2.652 km²), brood area and *Varroa* infestation level. Path analysis helps reconstruct the most plausible chain of causal links in multivariate datasets by assessing conditional independences among indirectly linked variables⁴¹,⁴². We identified the simplest path that did not deviate from conditional independence expectations while including only significant links (path analysis d-separation test, Fisher C = 30.08, df = 28, P = 0.36). This path model, that best captured the ecophysiological causal links behind overwintering survival, involved all studied variables (Fig. 1 and Table 1).

According to the path analysis, the initial colony level of brood had a direct and positive influence on *Varroa* infestation level and fat body mass of bees (GLMM with Gaussian distribution, df = 334, t = 3.16, p = 0.0017 and df = 169, t = 3.13, p = 0.0021, respectively; Fig. 1, Table 1 and supplementary Fig. S3). This could be easily explained by the fact that i) *Varroa* mites reproduce in brood cells, and ii) brood requires feeding by nurse bees, who take up nutrients from the fat body for secreting brood food via hypopharyngeal glands⁴³. In addition, brood production was strongly determined by the initial colony level of brood (df = 341, t = 12.67, p < 0.001; supplementary Fig. S4) but did not affect bee physiological traits.

Regarding bee health, our data showed that overwintering survival was positively influenced by vitellogenin level (GLMM with binomial family distribution, df = 170, z = 2.41, P = 0.016; Figs 1 and 2a), as previously evidenced in different environmental contexts²⁸,³⁴. This suggests that this phospholipoglycoprotein involved in survival traits such as oxidative stress resilience, cellular immunity, and longevity³⁷ can be employed as a predictive biomarker for monitoring honey bee populations. As expected, higher vitellogenin level was itself linked to higher fat body contents (df = 169, t = 2.94, P = 0.003; Figs 1 and 3) and overwintering was negatively influenced by *Varroa* infestation level during the pre-winter period (df = 170, z = -2.78, P = 0.005; Figs 1 and 2b). However, we did not find a link between vitellogenin and *Varroa* levels (t = -0.67, P = 0.50), contrary to a previous study, which showed that colonies that were not treated against *Varroa* (high mite infestation rate) had lower levels of vitellogenin in the fall compared to treated colonies³⁴. One probable explanation is that our colonies were all
<figure><img src="image_3.png" /><figcaption>**Figure 1. Path model revealing the ecophysiological basis of honey bee colony overwintering survival.**</figcaption></figure>

Significance level is indicated next to each link. All links stand for positive effects, except *Varroa* infestation level that negatively affects overwinter survival. When at least two links reach the same box, their thickness is proportional to their effect coefficient (standardized on their respective range for direct comparison). Total explained variance (r²) is indicated in the box for each response variable in the causal chain. Landscape influence on bee health is highlighted in grey.

treated against the mite by beekeepers and thus exhibited rather low *Varroa* infestation rates and negligible effects on vitellogenin levels⁴⁴.

Brood production, which was promoted by the presence of catch crops but not semi-natural habitats (t = 2.75, p = 0.0063 and t = 1.83, p = 0.068, respectively, Table 1 and supplementary Fig. S4), did not contribute to a better overwintering survival (df = 170, z = 2.32, p = 0.98). However, we found that landscape quality, whether melliferous catch crop or semi-natural habitat land covers, further facilitated the above-mentioned physiological causal chain, with a 1.6 and 1.9 times greater effect on vitellogenin and fat body, respectively, of semi-natural habitats as compared to catch crop (catch crop: fat body: df = 169, t = 2.84, P = 0.005; vitellogenin: t = 3.49, P < 0.001; semi-natural habitats: fat body: df = 169, t = 2.95, P = 0.003; vitellogenin: t = 3.10, P = 0.002; see Figs 1 and 4, and standardized estimates in Table 1 for effect size comparison). We therefore showed that bee health is affected by landscape patterns. More specifically, landscape enrichment with catch crops and semi-natural habitats promoted the physiological development of winter bees, which then enhanced the probability of overwintering survival. A similar result has been found in the United States when comparing the winter survival of colonies exposed to high or low potential forage²⁸⁻³⁰.

We could reconstruct a consistent ecophysiological causal chain behind overwintering survival, but there was still substantial variability inherent in the studied system (e.g. variability in bee physiological traits and brood production). In addition, many factors (e.g. pathogens, genotype and weather) can cause large variability in bee health and therefore induce winter losses³⁵. This likely explained why the studied landscape variables did not directly influence final overwintering survival (catch crop: df = 170, z = -1.61, P = 0.91; semi-natural habitats: df = 170, z = 1.53, P = 0.89; Table 1), even if they demonstrably improved the physiological state of bees during the pre-wintering period. Thus, the access to nutritional resources may be important for winter preparation but does not necessarily prevent the detrimental effects of other environmental and sanitary factors before or during the winter (e.g. *Varroa* effects, as shown above and by Dolezal *et al.*³⁰).

This study provides a first tentative reconstruction of the ecophysiological basis of overwintering survival. Although the number of different landscape contexts actually covered in our study is currently moderate (18 contexts, including 13 with detailed physiological data), we could recover here the most important correlate of overwintering mortality (*Varroa* mites). We have also shown that the landscape quality is liable to influence the physiological state of honey bees. This candidate model should however be reevaluated with a greater amount of apiaries, enlarged range of biogeographical contexts and foraging distances to confirm that physiological state may be used as a surrogate of future colony survival or collapse risk. This should be coupled with a scale dependency analysis to model and assess the optimal spatial management grain for effective landscape restoration.

**Effects of landscape enrichment on pollen diets.** We performed palynological identifications and nutritional evaluation as an *a posteriori* analysis for confirming the use of melliferous catch crop by bees, and as an additional source of information to help understand the observed links between landscape and colony state of health. Landscape enrichment with the melliferous catch crop treatment (Table 2) improved the physiological state of bees during the pre-winter period, likely due to higher pollen diet diversity rather than better nutritional abundance and quality (Table 3). In the absence of pollen resource landscape enhancement, pollen diet was
<table><caption>**Table 1. Statistical details of selected or missing paths coefficients from the path model analysis.** The estimate standardized to data range and corresponding P-value are shown for each path. The path model includes survival, physiological, *Varroa* infestation, brood and landscape variables.</caption><thead><tr><th>Response</th><th>Predictor</th><th>Estimate</th><th>SE</th><th>DF</th><th>Statistics</th><th>P</th></tr></thead><tbody><tr><td colspan="7">Selected path coefficients (estimates standardized to data range)</td></tr><tr><td>Overwintering survival</td><td>Varroa infestation</td><td>-5.552</td><td>1.996</td><td>170</td><td>z = -2.782</td><td>0.0054</td></tr><tr><td>Overwintering survival</td><td>Vitellogenin</td><td>3.046</td><td>1.265</td><td>170</td><td>z = 2.408</td><td>0.016</td></tr><tr><td>Varroa infestation</td><td>Initial brood area</td><td>0.052</td><td>0.017</td><td>334</td><td>t = 3.161</td><td>0.0017</td></tr><tr><td>Vitellogenin</td><td>Melliferous catch crops</td><td>0.114</td><td>0.033</td><td>169</td><td>t = 3.486</td><td>0.0006</td></tr><tr><td>Vitellogenin</td><td>Semi-natural habitats</td><td>0.178</td><td>0.057</td><td>169</td><td>t = 3.105</td><td>0.0021</td></tr><tr><td>Vitellogenin</td><td>Fat body</td><td>0.201</td><td>0.068</td><td>169</td><td>t = 2.944</td><td>0.0036</td></tr><tr><td>Fat body</td><td>Initial brood area</td><td>0.162</td><td>0.052</td><td>169</td><td>t = 3.131</td><td>0.0021</td></tr><tr><td>Fat body</td><td>Semi-natural habitats</td><td>0.186</td><td>0.063</td><td>169</td><td>t = 2.948</td><td>0.0037</td></tr><tr><td>Fat body</td><td>Melliferous catch crops</td><td>0.099</td><td>0.035</td><td>169</td><td>t = 2.836</td><td>0.0051</td></tr><tr><td>Final brood area</td><td>Initial brood area</td><td>0.554</td><td>0.044</td><td>341</td><td>t = 12.671</td><td>&lt;0.001</td></tr><tr><td>Final brood area</td><td>Melliferous catch crops</td><td>0.109</td><td>0.04</td><td>341</td><td>t = 2.749</td><td>0.0063</td></tr><tr><td colspan="7">Missing path coefficients (estimates standardized to data range)</td></tr><tr><td>Final brood area</td><td>Semi-natural habitats</td><td>0.1404</td><td>0.0766</td><td>340</td><td>t = 1.8322</td><td>0.0678</td></tr><tr><td>Varroa infestation</td><td>Semi-natural habitats</td><td>-0.1227</td><td>0.2221</td><td>333</td><td>t = -0.5526</td><td>0.5809</td></tr><tr><td>Overwintering survival</td><td>Semi-natural habitats</td><td>0.7922</td><td>0.5186</td><td>170</td><td>z = 1.5276</td><td>0.8994</td></tr><tr><td>Varroa infestation</td><td>Melliferous catch crops</td><td>1.7131</td><td>1.2636</td><td>333</td><td>t = 1.3557</td><td>0.1761</td></tr><tr><td>Overwintering survival</td><td>Melliferous catch crops</td><td>-3.4473</td><td>2.1434</td><td>170</td><td>z = -1.6083</td><td>0.9143</td></tr><tr><td>Vitellogenin</td><td>Initial brood area</td><td>0.1431</td><td>0.0966</td><td>168</td><td>t = 1.4806</td><td>0.1406</td></tr><tr><td>Overwintering survival</td><td>Initial brood area</td><td>0.0367</td><td>0.2319</td><td>170</td><td>z = 0.1582</td><td>0.3832</td></tr><tr><td>Fat body</td><td>Final brood area</td><td>0.5433</td><td>0.552</td><td>168</td><td>t = 0.9843</td><td>0.3263</td></tr><tr><td>Vitellogenin</td><td>Final brood area</td><td>-0.1932</td><td>0.1111</td><td>167</td><td>t = -1.7386</td><td>0.0839</td></tr><tr><td>Varroa infestation</td><td>Final brood area</td><td>-0.2167</td><td>0.1587</td><td>330</td><td>t = -1.3656</td><td>0.173</td></tr><tr><td>Overwintering survival</td><td>Final brood area</td><td>0.6236</td><td>0.269</td><td>170</td><td>z = 2.3183</td><td>0.9837</td></tr><tr><td>Varroa infestation</td><td>Fat body</td><td>-0.0224</td><td>0.0323</td><td>167</td><td>t = -0.6945</td><td>0.4883</td></tr><tr><td>Overwintering survival</td><td>Fat body</td><td>-0.0548</td><td>0.0432</td><td>170</td><td>z = -1.2679</td><td>0.838</td></tr><tr><td>Varroa infestation</td><td>Vitellogenin</td><td>-0.1077</td><td>0.1606</td><td>166</td><td>t = -0.6709</td><td>0.5032</td></tr></tbody></table>

largely dominated by climbing ivy (*Hedera helix*) (71% of pollen diet volume on average, and >90% in half of the samples) (Tables 2 and 3). When available in the foraging range, melliferous catch crop accounted for, on average, 52% of a colony total pollen diet composition, thereby offsetting significantly the broad dominance of climbing ivy pollen (−39% relative change) and increasing significantly pollen species diversity (+87% in species diversity, +77% in species evenness, Table 3). However, the daily pollen intake did not increase significantly and the overall pollen diet energy and protein contents remained unchanged (Table 3). Although, we could not exclude the possibility that the diets provided by the melliferous catch crop treatment were of higher quality regarding other nutrients (e.g. lipids, amino acids, vitamins).

Monofloral and diversified diet, have long been suspected to differentially affect bee health. While some studies revealed, via experimental manipulation of diet composition, that pollen diet diversity increases bee immunocompetence³⁹ and reduces disease susceptibility⁴⁰,⁴⁵, it still remains unclear whether the benefits of resource diversity apply at the field level, which is the key for supporting decision-making toward the management of natural resources. Although dedicated experiments will be needed to elucidate this question, our results suggest that pollen diversity might provide benefits to bee health in natural conditions and at the landscape level. Those results were further supported by the even stronger effect of semi-natural habitats, which generally offer a great diversity of floral resources⁹,⁴⁶,⁴⁷.

The underlying mechanisms linking pollen diversity to bee health have yet to be determined. However, De Groot showed that honey bees require a set of essential amino acids in specific proportions for normal growth and development⁴⁸. In addition, bumble bee larvae become heavier when fed with polyfloral pollen diets compared to larvae fed with monofloral diets, even with higher protein content⁴⁹. Therefore, it is possible that an increase in environmental plant species diversity optimizes the occurrence, diversity and/or proportion of specific pollen nutrients (proteins, amino acids, lipids, starches, sterols, vitamins and minerals) that are required for the development of certain physiological traits, such as fat body and vitellogenin production.

## Conclusions

The ubiquity of habitat degradation requires identifying potential landscape patterns that may act as stressors for bees and providing recommendations for habitat restoration¹⁰,⁵⁰. Within the scope of this experimental design, we found that bee health is better influenced by semi-natural habitats than by landscape enrichment with catch crops. This suggests that, when considering habitat restoration, artificial bee pastures may be designed as a complementary management measure intended to support semi-natural habitat protection and restoration.
<figure><img src="image_4.png" /><figcaption>**Figure 2.** Representation of the path model links showing the influence of *Varroa* infestation and vitellogenin levels on overwintering survival. Survival probability was influenced positively by vitellogenin level (a) and negatively by *Varroa* infestation level (b). The continuous lines show model predictions. For the binary response variable (colony overwintering survival), data are represented as mean ± SE after being pooled into groups of consistent sizes. Tick marks show the position of raw data along the horizontal axis.</figcaption></figure>

<figure><img src="image_5.png" /><figcaption>**Figure 3.** Relationship between fat body content and vitellogenin levels.</figcaption></figure>

Our study also indicates that applying ecophysiological approaches to honey bee and native bee conservation might be complementary to the more conventional distribution-based studies, which have shown that agri-environmental schemes promoting the conservation of semi-natural habitats or the development of flower-rich field margins in farming areas are favourable to bee abundance and diversity⁵¹⁻⁵⁵ (but see ref. 56).
<figure><img src="image_6.png" /><figcaption>**Figure 4.** Representation of the path model links showing the influence of landscape quality variables on bee physiological traits. Melliferous catch crop and semi-natural habitats positively influenced the fat body (a) and vitellogenin levels (b). Fat body and vitellogenin values were averaged per apiary and plotted as a function of landscape metrics. Trends are depicted by the regression planes. A slight horizontal jitter was applied to separate overlying apiary data with equal flowering catch crop treatments.</figcaption></figure>

Adding a physiological dimension to the environmental variables will benefit the assessment of population health and sustainability<sup>20,21</sup>. In conclusion, this work highlights landscape ecophysiology as a promising field of research for better understanding the influence of the environment on pollinator health and setting the stage for more effective pollinator conservation.

## Methods

**Experimental set-up and colony monitoring.** We monitored the overwintering survival of 350 honey bee colonies (163 and 187 colonies over winters 2012–13 and 2013–14, respectively) in an intensive farming system from central western France, Centre French region (see supplementary Table S1, Figs S1 and S2). Colonies belonged to three volunteer professional beekeepers from different parts of the region located 80 to 110 km apart. In each beekeeping set-up and each year, a 1.5 km radius area (approximate honey bee foraging range in autumn<sup>57</sup>) was experimentally enriched with 5.0 to 31.5 ha of melliferous catch crops during the pre-wintering period (mid-September to mid-October). Colonies were set up to prepare for winter either inside the melliferous areas (n = 184 colonies split into 10 apiaries) or 8–10 km farther away with exclusively non-melliferous catch crop fields within the foraging range (n = 166 colonies split into 8 apiaries), following a control-vs-treatment experimental design with 2 to 4 simultaneous apiary monitoring (supplementary Table S1 and Figs S1 and S2). Details on the species composition of catch crops (*Avena sativa*, *Brassica juncea*, *Helianthus annuus*, *Phacelia tanacetifolia*, *Sinapis alba*, *Trifolium alexandrinum*, *Vicia benghalensis*, *Vicia sativa*) are shown in the supplementary Table S2. The allocation of colonies among apiaries was random, but all colonies within a given paired (control-vs.-treatment) experimental set-up originated from the same professional beekeeper, and therefore had shared the same management history and honey production (rapeseed and acacia) during the season. Colonies were different between years.

The monitoring started during the week preceding the expected onset of catch crop flowering (mid to end of September, depending on the study year). All colonies received a standardized *Varroa* mite treatment (APIVAR®,
<table><caption>**Table 2.** Average species and morphotype composition of the pollen volume collected by honeybees in melliferous and non-melliferous catch crop treatments. In the "Other pollen species" section, pollen originating from crop plants are indicated by an asterisk; Other pollens are likely originating from semi-natural habitats.</caption><thead><tr><th>Pollen species or morphotypes</th><th>Melliferous catch crop treatments</th><th>Absence of melliferous catch crop treatments</th></tr></thead><tbody><tr><td colspan="3">Pollen species or morphotypes potentially originating from the melliferous catch crops</td></tr><tr><td>Helianthus annuus (Asteraceae)</td><td>0.1%</td><td></td></tr><tr><td>Phacelia tanacetifolia (Hydrophyllaceae)</td><td>3.2%</td><td></td></tr><tr><td>Trifolium alexandrinum (Fabaceae)</td><td>16.9%</td><td></td></tr><tr><td>type Sinapis/Brassica sp. (Brassicaceae)</td><td>32.0%</td><td>2.3%</td></tr><tr><td colspan="3">Other pollen species</td></tr><tr><td>Ammi majus (Apiaceae)</td><td>&lt;0.1%</td><td></td></tr><tr><td>Castanea sativa (Fagaceae)</td><td>&lt;0.1%</td><td>&lt;0.1%</td></tr><tr><td>Datura stramonium (Solanaceae)</td><td></td><td>&lt;0.1%</td></tr><tr><td>Fagopyrum esculentum (Polygonaceae)*</td><td></td><td>&lt;0.1%</td></tr><tr><td>Guizotia abyssinica (Asteraceae)*</td><td>&lt;0.1%</td><td>0.1%</td></tr><tr><td>Hedera helix (Araliaceae)</td><td>43.7%</td><td>73.3%</td></tr><tr><td>Mercurialis annua (Euphorbiaceae)</td><td></td><td>2.0%</td></tr><tr><td>Papaver rhoeas (Papaveraceae)</td><td>&lt;0.1%</td><td></td></tr><tr><td>Plantago lanceolata (Plantaginaceae)</td><td></td><td>0.1%</td></tr><tr><td>Raphanus sativus (Brassicaceae)*</td><td>3.7%</td><td>18.9%</td></tr><tr><td>Reseda lutea (Resedaceae)</td><td>0.1%</td><td>0.1%</td></tr><tr><td>type Bellis perennis (Asteraceae)</td><td></td><td>&lt;0.1%</td></tr><tr><td>type Brassica napus (Brassicaceae)*</td><td>0.1%</td><td></td></tr><tr><td>type Calendula sp. (Asteraceae)</td><td></td><td>0.1%</td></tr><tr><td>type Cichorium sp. (Asteraceae)</td><td>0.1%</td><td>0.6%</td></tr><tr><td>type Crepis sp. (Asteraceae)</td><td></td><td>0.1%</td></tr><tr><td>type Pinus sp. (Pinaceae)</td><td>&lt;0.1%</td><td>&lt;0.1%</td></tr><tr><td>type Raphanus raphanistrum (Brassicaceae)</td><td>0.6%</td><td>1.9%</td></tr><tr><td>type Rhamnus sp. (Rhamnaceae)</td><td>&lt;0.1%</td><td></td></tr><tr><td>type Rubus sp. (Rosaceae)</td><td>0.1%</td><td></td></tr><tr><td>type Solanum sp. (Solanaceae)*</td><td>0.1%</td><td>0.6%</td></tr><tr><td>type Trifolium pratense (Fabaceae)</td><td>&lt;0.1%</td><td></td></tr><tr><td>type Veronica sp. (Scrophulariaceae)</td><td></td><td>0.1%</td></tr><tr><td>type Viola sp. (Violaceae)</td><td></td><td>0.7%</td></tr></tbody></table>

Véto-pharma) and their initial state of brood development was documented by tallying the amount of beehive frame sides (out of 20) covered with brood. Four weeks later, at the end of flowering (mid to end of October), brood area was reassessed and adult honey bees were sampled to determine their physiological state and *Varroa* infestation levels (see methods below). Finally, as cold days arrived, colonies were gathered back into their respective beekeeping domain to standardize overwintering conditions. Overwintering was considered successful for colonies that remained operational for beekeeping activity in early spring, therefore excluding collapsed and weak or orphan colonies. Colony survival was related with physiological state, brood area, *Varroa* infestation and the surrounding landscape characteristics.

**Landscape quality.** Beside catch crop land cover, landscape quality was documented by quantifying the extent of permanent semi-natural habitat (woodlots and hedgerows) within the 1.5 km foraging range around the apiaries. Geographical information on semi-natural habitats was obtained from the French national remote sensing database on vegetation layers (Institut Géographique National) and processed with the Quantum GIS mapping software version 2.2.

**Physiological traits and Varroa infestation assessment.** To determine the phoretic *Varroa* mite infestation rates of colonies, around 200 bees were collected and washed with soapy water (TEEPOL) to dislodge the mites for counting⁵⁸. Infestation rate was reported as the number of mites per 100 adult bees (n = 339 colonies from 18 apiaries).

Additionally, at the end of the flowering period, around 100 adult bees per colony (n = 175 colonies from 13 apiaries) were collected on brood frames, placed in dry ice and stored at -80 °C. Fat body quantification was performed on a pool of 30 abdomens using the ether extraction method described in⁵⁹. The expression level of vitellogenin was determined by quantitative RT-PCR. For each colony, three pools of 10 abdomens were each homogenized in 1 ml of Qiazol reagent (Qiagen) with a TissueLyser (Qiagen) (4 × 30 s at 30 Hz). The homogenates were incubated for 5 min at room temperature and after centrifugation (12,000 g for 30 s at 4 °C) the 3
<table><caption>**Table 3.** Comparison of mean (±SD) pollen diet characteristics between catch crop treatments. MCC: Melliferous catch crop; Species richness S stands for the number of species or morphospecies; Species diversity H' approaches 0 as differences increase among pollen species (or morphospecies) relative abundances, and gets higher as relative abundances become more equal; Species evenness J' standardizes H' back to the range [0, 1] for easier interpretation, with a value of 1 meaning perfect equidistribution of abundances among species.</caption><thead><tr><td>Pollen diet properties</td><td>MCC treatments</td><td>Absence of MCC treatments</td><td>Kruskal-Wallis statistics</td></tr></thead><tbody><tr><td>Daily pollen trap content (g)</td><td>23.82 ± 33.13</td><td>13.19 ± 16.82</td><td>χ² = 0.68, n = 25, P = 0.41</td></tr><tr><td>Pollen diet volume (%) from MCC species</td><td>51.99 ± 27.91</td><td>2.35 ± 5.58</td><td>χ² = 16.74, n = 25, P < 0.001</td></tr><tr><td>Species richness (S)</td><td>5.2 ± 2.1</td><td>4.4 ± 1.8</td><td>χ² = 0.83, n = 25, P = 0.36</td></tr><tr><td>Species diversity (Shannon H')</td><td>1.18 ± 0.53</td><td>0.63 ± 0.48</td><td>χ² = 5.47, n = 25, P = 0.019</td></tr><tr><td>Species evenness (Pielou J')</td><td>0.55 ± 0.24</td><td>0.31 ± 0.23</td><td>χ² = 4.97, n = 25, P = 0.026</td></tr><tr><td>Energy content (kcal/kg of dry matter)</td><td>5085.6 ± 266.7</td><td>5024.4 ± 161.9</td><td>χ² = 0.14, n = 17, P = 0.70</td></tr><tr><td>Protein content (% of dry matter)</td><td>27.78 ± 2.00</td><td>27.85 ± 2.58</td><td>χ² = 0.15, n = 17, P = 0.69</td></tr><tr><td>Hedera helix pollen diet volume (%)</td><td>43.33 ± 27.91</td><td>71.42 ± 33.87</td><td>χ² = 4.05, n = 25, P < 0.044</td></tr></tbody></table>

supernatants were pooled (167 µl each, giving a 501 µl supernatant). RNA extraction was then carried out as indicated in the RNeasy Plus Universal kit (Qiagen). cDNA synthesis and analysis of vitellogenin expression level was performed as in⁴⁰. Cycle threshold values of vitellogenin were normalized to the geometric mean of the housekeeping genes actin and eIF3-S8 using the comparative quantification method (delta Ct method). We used published sequences of primers for vitellogenin⁶⁰, actin⁶¹ and eIF3-S8⁶².

**Palynological validation and nutritional characteristics.** To further ascertain the use of melliferous catch crops by foraging honey bees and to document the possible contribution of pre-wintering nutrition to bee health and overwintering survival, we sampled pollen from one to five (depending on apiary size) randomly chosen colonies per apiary (n = 27 colonies from 10 apiaries). On the second week of catch crop flowering, standard pollen traps were placed at the hive entrances to collect pollen loads during three to six consecutive days (depending on microclimatic conditions), in order to cover three sampling days suitable for foragers (temperature >15 °C, low wind and no rain) and to reach at least 15 g pollen samples. Samples were stored at -20 °C. Two 4-g subsamples per colony were used for the pollen species identification as described in Requier *et al.*⁴⁶. Briefly, pollen samples were diluted in water and mounted onto microscope slides, which were stained with Fuschin, and examined at 400x magnification. Pollen was identified to genus and, when possible, to species. At least 300 pollen grains were counted and identified on each slide. The identification process was duplicated and averaged for improving accuracy. The relative contribution of each pollen species to the total volume of collected pollen in a sample was determined based on pollen grain size⁴⁶. Protein and energy content were determined on 5 g of homogenized pollen subsamples the year of sampling as in Requier *et al.*⁴⁶.

**Statistical analyses.** Data were analyzed using generalized linear mixed models (GLMMs) with either a binomial family error distribution for binary data (overwintering success) or a Gaussian family for quantitative data, beforehand log₂-corrected to recover normality whenever necessary (brood area, vitellogenin level). By specifying a random grouping structure in the dataset, GLMMs allow proper accounting for the non-independency of colonies from the same apiary and the beekeeping set-up.

We then produced a tentative path model linking survival with the studied variables: physiological traits, *Varroa* infestation, brood area and landscape. In particular, we expected that *i*)* overwintering survival would increase with vitellogenin level, itself increasing with fat body content and that *ii*)* each step of this causal chain would be potentially influenced either negatively by *Varroa* infestations, or positively by landscape quality and initial or final (post-flowering) brood area. We also considered the possible concomitant effects of initial brood area and landscape quality on final brood area before winter. Once computed, the tentative path model was refined by dropping non-significant links and by sequentially adding any link that was initially ignored until the path model was judged statistically supported by the data. New links were added by order of increasing P-value, and deviation from expected conditional independence assessed using the d-separation test¹¹ specially suited for generalized mixed models⁴². All quantitative explanatory variables were standardized to a range [0, 1], so that coefficient estimates can be readily compared to determine the most influential explanatory variables in the candidate path models.

We finally performed a *posteriori* comparisons of pollen diet characteristics (pollen intake (g.day⁻¹), composition (species richness S, Shannon diversity index *H*' and Pielou's evenness index *J*') and nutritional properties (energy and protein content)) among catch crop treatments using Kruskal-Wallis tests. All analyses were performed with the R software version 3.1.1⁶³.

**References**

1. Biesmeijer, J. C. et al. Parallel declines in pollinators and insect-pollinated plants in Britain and the Netherlands. *Science* **313**, 351–354 (2006).

2. Goulson, D., Nicholls, E., Botias, C. & Rotheray, E. L. Bee declines driven by combined stress from parasites, pesticides, and lack of flowers. *Science* **347**, 1255957 (2015).

3. Ollerton, J., Erenler, H., Edwards, M. & Crockett, R. Extinctions of aculeate pollinators in Britain and the role of large-scale agricultural changes. *Science* **346**, 1360–1362 (2014).

4. Winfree, R., Bartomeus, I. & Cariveau, D. P. Native Pollinators in Anthropogenic Habitats. *Ann Rev Ecol Evol Syst* **42**, 1–22 (2011).
5. Potts, S. G. et al. Global pollinator declines: trends, impacts and drivers. *Trends Ecol Evol* **25**, 345–353 (2010).

6. Kremen, C. et al. Pollination and other ecosystem services produced by mobile organisms: a conceptual framework for the effects of land-use change. *Ecol Lett* **10**, 299–314 (2007).

7. Klein, A. M. et al. Importance of pollinators in changing landscapes for world crops. *Proc R Soc Lond B* **274**, 303–313 (2007).

8. Vaudo, A. D., Tooker, J. F., Grozinger, C. M. & Patch, H. M. Bee nutrition and floral resource restoration. *Curr Opin Insect Sci* **10**, 133–141 (2015).

9. Wratten, S. D., Gillespie, M., Decourtye, A., Mader, E. & Desneux, N. Pollinator habitat enhancement: Benefits to other ecosystem services. *Agr Ecosyst Environ* **159**, 112–122 (2012).

10. Decourtye, A., Mader, E. & Desneux, N. Landscape enhancement of floral resources for honey bees in agro-ecosystems. *Apidologie* **41**, 264–277 (2010).

11. Batary, P., Dicks, L. V., Kleijn, D. & Sutherland, W. J. The role of agri-environment schemes in conservation and environmental management. *Conserv Biol* **29**, 1006–1016 (2015).

12. Henry, M. et al. Spatial autocorrelation in honeybee foraging activity reveals optimal focus scale for predicting agro-environmental scheme efficiency. *Ecol Model* **225**, 103–114 (2012).

13. Decourtye, A. et al. Why Enhancement of Floral Resources in Agro-Ecosystems Benefit Honeybees and Beekeepers? In *Ecosystems Biodiversity* (eds Grillo, O. & Venora, G.) Ch. 16, 371–388 (InTech, 2011).

14. Scheper, J. et al. Local and landscape-level floral resources explain effects of wildflower strips on wild bees across four European countries. *J Appl Ecol* **52**, 1165–1175 (2015).

15. Kleijn, D. et al. Mixed biodiversity benefits of agri-environment schemes in five European countries. *Ecol Lett* **9**, 243–254 (2006).

16. Scheper, J. et al. Environmental factors driving the effectiveness of European agri-environmental measures in mitigating pollinator loss - a meta-analysis. *Ecol Lett* **16**, 912–920 (2013).

17. Holzschuh, A., Steffan-Dewenter, I., Kleijn, D. & Tscharntke, T. Diversity of flower-visiting bees in cereal fields: effects of farming system, landscape composition and regional context. *J Appl Ecol* **44**, 41–49 (2007).

18. Kleijn, D. et al. Ecological effectiveness of agri-environment schemes in different agricultural landscapes in the Netherlands. *Conserv Biol* **18**, 775–786 (2004).

19. Williams, N. M., Minckley, R. L. & Silveira, F. A. Variation in native bee faunas and its implications for detecting community changes. *Conserv Ecol* **5**, 1–21 (2001).

20. Ellis, R. D., McWhorter, T. J. & Maron, M. Integrating landscape ecology and conservation physiology. *Landscape Ecol* **27**, 1–12 (2012).

21. Cooke, S. J. et al. What is conservation physiology? Perspectives on an increasingly integrated and essential science. *Conserv Physiol* **1**, 1–23 (2013).

22. Madliger, C. L. & Love, O. P. The Power of Physiology in Changing Landscapes: Considerations for the Continued Integration of Conservation and Physiology. *Integr Comp Biol* **5**, 545–553 (2015).

23. Crone, E. E. & Williams, N. M. Bumble bee colony dynamics: quantifying the importance of land use and floral resources for colony growth and queen production. *Ecol Lett* **19**, 460–468 (2016).

24. Goulson, D., Hughes, W. O. H., Derwent, L. C. & Stout, J. C. Colony growth of the bumblebee, *Bombus terrestris*, in improved and conventional agricultural and suburban habitats. *Oecologia* **130**, 267–273, (2002).

25. Williams, N. M., Regetz, J. & Kremen, C. Landscape-scale resources promote colony growth but not reproductive performance of bumble bees. *Ecology* **93**, 1049–1058 (2012).

26. Odoux, J. F. et al. ECOBEE: a tool for long-term honey bee colony monitoring at the landscape scale in West European intensive agroecosystems. *J Apicul Res* **53**, 57–66 (2014).

27. Smith, G. W. et al. Bee Abundance and Nutritional Status in Relation to Grassland Management Practices in an Agricultural Landscape. *Environ Entomol* **45**, 338–347 (2016).

28. Smart, M., Pettis, J., Rice, N., Browning, Z. & Spivak, M. Linking Measures of Colony and Individual Honey Bee Health to Survival among Apiaries Exposed to Varying Agricultural Land Use. *Plos One* **11**, e0152685 (2016).

29. Smart, M. D., Pettis, J. S., Euliss, N. & Spivak, M. S. Land use in the Northern Great Plains region of the US influences the survival and productivity of honey bee colonies. *Agr Ecosyst Environ* **230**, 139–149 (2016).

30. Dolezal, A. G., Carrillo-Tripp, J., Miller, W. A., Bonning, B. C. & Toth, A. L. Intensively Cultivated Landscape and Varroa Mite Infestation Are Associated with Reduced Honey Bee Nutritional State. *Plos One* **11**, e0153531 (2016).

31. Steinhauer, N. A. et al. A national survey of managed honey bee 2012–2013 annual colony losses in the USA: results from the Bee Informed Partnership. *J Apicul Res* **53**, 1–18 (2014).

32. Laurent, M., Hendrikx, P., Ribière-Chabert, M. & Chauzat, M.-P. *A pan-European epidemiological study on honeybee colony losses 2012–2014*. p.44 (ANSES, 2016).

33. Hahn, D. A. & Denlinger, D. L. Meeting the energetic demands of insect diapause: nutrient storage and utilization. *J Insect Physiol* **53**, 760–773 (2007).

34. Dainat, B., Evans, J. D., Chen, Y. P., Gauthier, L. & Neumann, P. Predictive markers of honey bee colony collapse. *Plos One* **7**, e32151 (2012).

35. Döke, M. A., Frazier, M. & Grozinger, C. M. Overwintering Honey Bees: Biology and Management. *Curr Opin Insect Sci* **10**, 185–193 (2015).

36. Amdam, G. V., Fennern, E. & Havukainen, H. Vitellogenin in Honey Bee Behavior and Lifespan In *Honeybee Neurobiology and Behavior* (eds Galizia, C.G., Eisenhardt, D. & Giurfa, M.) 17–29 (Springer, 2011).

37. Seehuus, S. C., Norberg, K., Gimsa, U., Krekling, T. & Amdam, G. V. Reproductive protein protects functionally sterile honey bee workers from oxidative stress. *Proc Natl Acad Sci USA* **103**, 962–967 (2006).

38. Fluri, P., Luscher, M., Wille, H. & Gerig, L. Changes in the weight of the pharyngeal gland and haemolymph titres of juvenile hormone, protein and vitellogenin in worker honey bees. *J Insect Physiol* **28**, 61–68 (1982).

39. Alaux, C., Ducloz, F., Crauser, D. & Le Conte, Y. Diet effects on honeybee immunocompetence. *Biol Lett* **6**, 562–565 (2010).

40. Di Pasquale, G. et al. Influence of pollen nutrition on honey bee health: Do pollen quality and diversity matter? *PLoS One* **8**, e72016 (2013).

41. Shipley, B. The AIC model selection method applied to path analytic models compared using a d-separation test. *Ecology* **94**, 560–564 (2013).

42. Shipley, B. Confirmatory path analysis in a generalized multilevel context. *Ecology* **90**, 363–368 (2009).

43. Crailsheim, K. The protein balance of the honey bee worker. *Apidologie* **21**, 417–429 (1990).

44. Amdam, G. V., Hartfelder, K., Norberg, K., Hagen, A. & Omholt, S. W. Altered physiology in worker honey bees (Hymenoptera: Apidae) infested with the mite *Varroa destructor* (Acari: Varroidae): A factor in colony loss during overwintering? *J Econ Entomol* **97**, 741–747 (2004).

45. Foley, K., Fazio, G., Jensen, A. B. & Hughes, W. O. Nutritional limitation and resistance to opportunistic Aspergillus parasites in honey bee larvae. *J Invertebr Pathol* **111**, 68–73 (2012).

46. Requier, F. et al. Honey bee diet in intensive farmland habitats reveals an unexpectedly high flower richness and a major role of weeds. *Ecol Appl* **25**, 881–890 (2015).

47. Odoux, J. F. et al. Territorial biodiversity and consequences on physico-chemical characteristics of pollen collected by honey bee colonies. *Apidologie* **43**, 561–575 (2012).
48. de Groot, A. P. Protein and amino acid requirements of the honey bee (Apis mellifera L.). *Physiol Comp Oecol* **3**, 197-285 (1953).

49. Tasei, J. N. & Aupinel, P. Nutritive value of 15 single pollens and pollen mixes tested on larvae produced by bumblebee workers (Bombus terrestris, Hymenoptera: Apidae). *Apidologie* **39**, 397-409 (2008).

50. Kleijn, D. & Sutherland, W. J. How effective are European agri-environment schemes in conserving and promoting biodiversity? *J Appl Ecol* **40**, 947-969 (2003).

51. Blaauw, B. R. & Isaacs, R. Flower plantings increase wild bee abundance and the pollination services provided to a pollination-dependent crop. *J Appl Ecol* **51**, 890-898 (2014).

52. Carvell, C., Meek, W. R., Pywell, R. F., Goulson, D. & Nowakowski, M. Comparing the efficacy of agri-environment schemes to enhance bumble bee abundance and diversity on arable field margins. *J Appl Ecol* **44**, 29-40 (2007).

53. Kennedy, C. M. et al. A global quantitative synthesis of local and landscape effects on wild bee pollinators in agroecosystems. *Ecol Lett* **16**, 584-599 (2013).

54. Pywell, R. F. et al. Effectiveness of new agri-environment schemes in providing foraging resources for bumblebees in intensively farmed landscapes. *Biol Conserv* **129**, 192-206 (2006).

55. Williams, N. M. et al. Native wildflower plantings support wild bee abundance and diversity in agricultural landscapes across the United States. *Ecol Appl* **25**, 2119-2131 (2015).

56. Wood, T. J., Holland, J. M. & Goulson, D. Pollinator-friendly management does not increase the diversity of farmland bees and wasps. *Biol Conserv* **187**, 120-126 (2015).

57. Couvillon, M. J., Schurch, R. & Ratnieks, F. L. W. Waggle Dance Distances as Integrative Indicators of Seasonal Foraging Challenges. *PLoS One* **9**, e93495 (2014).

58. Fries, I., Aarhus, A., Hansen, H. & Korpela, S. Comparisons of diagnostic methods for detection of Varroa jacobsoni in honey bee (Apis mellifera) colonies at low infestation levels. *Exp Appl Acarol* **10**, 279-287 (1991).

59. Wilson-Rich, N., Dres, S. T. & Starks, P. T. The ontogeny of immunity: development of innate immune strength in the honey bee (Apis mellifera). *J Insect Physiol* **54**, 1392-1399 (2008).

60. Fischer, P. & Grozinger, C. M. Pheromonal regulation of starvation resistance in honey bee workers (Apis mellifera). *Naturwissenschaften* **95**, 723-729 (2008).

61. Amdam, G. V., Norberg, K., Fondrk, M. K. & Page, R. E. Jr. Reproductive ground plan may mediate colony-level selection effects on individual foraging behavior in honey bees. *Proc Natl Acad Sci USA* **101**, 11350-11355 (2004).

62. Ament, S. A. et al. Mechanisms of stable lipid loss in a social insect. *J Exp Biol* **214**, 3808-3821 (2011).

63. Team, R. D. C. R: A language and environment for statistical computing. (URL http://www.R-project.org/, Vienna, Austria, 2014).

## Acknowledgements

Special thanks go to Nathalie Moreau for nutritional analyses, Cyril Vidau for field assistance, Marion Salignon and Lucie Salem for help in physiological analysis, Sarah Rivière for support on landscape analysis, two anonymous referees for comments that improved the manuscript and Cynthia McDonnell for the English editing. We are also grateful to the volunteer farmers and beekeepers involved in the study. This work was funded by the French Ministry of Agriculture program (CASDAR, INTERAPI program n° 1176) and a grant from INRA SPE department (CLIMBEE).

## Author Contributions

C.A., F.A. and M.H. conceived the study, F.A., A.D., J.-F.O., T.T., M.C., E.D. and F.D. conducted the experiments, C.A. and M.H. analyzed the data, Y.L.C. contributed to reagents, C.A. and M.H. wrote the manuscript.

## Additional Information

**Supplementary information** accompanies this paper at http://www.nature.com/srep

**Competing financial interests:** The authors declare no competing financial interests.

**How to cite this article:** Alaux, C. et al. A 'Landscape physiology' approach for assessing bee health highlights the benefits of floral landscape enrichment and semi-natural habitats. *Sci. Rep.* **7**, 40568; doi: 10.1038/srep40568 (2017).

**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_7.png" /></figure>

This work is licensed under a Creative Commons Attribution 4.0 International License. The images or other third party material in this article are included in the article's Creative Commons license, unless indicated otherwise in the credit line; if the material is not included under the Creative Commons license, users will need to obtain permission from the license holder to reproduce the material. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/

© The Author(s) 2017