# Rapid responses of winter aphid-parasitoid communities to climate warming

Kévin Tougeron, Maxime Damien, Cécile Le Lann, Jacques Brodeur, Joan van Baaren

► To cite this version:

Kévin Tougeron, Maxime Damien, Cécile Le Lann, Jacques Brodeur, Joan van Baaren. Rapid responses of winter aphid-parasitoid communities to climate warming. Frontiers in Ecology and Evolution, 2018, 6 (OCT), pp.173. 10.3389/fevo.2018.00173. hal-01937630

HAL Id: hal-01937630

https://univ-rennes.hal.science/hal-01937630

Submitted on 18 Jul 2019

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
<figure><img src="image_2.png" /></figure>

# Rapid Responses of Winter Aphid-Parasitoid Communities to Climate Warming

Kévin Tougeron<fnref n="1" />,²*, Maxime Damien<fnref n="1" />,³, Cécile Le Lann<fnref n="1" />, Jacques Brodeur² and Joan van Baaren<fnref n="1" />

Consequences of inter-annual environmental fluctuations, including those associated with climate change, can have a knock-on effect from individual to community scale. In particular, changes in species seasonal phenology can modify the structure and composition of communities, with potential consequences on their functioning and the provision of ecosystem services. In mild climate areas, aphids can be present in cereal fields throughout the winter, which allows aphid parasitoids to remain active. Using a 9-year dataset of aphid-parasitoid winter trophic webs in cereal fields of Western France, we report that the community structure and composition that prevailed before 2011 have recently shifted toward a more diversified community, with the presence of two new braconid parasitoid species (*Aphidius ervi* and *Aphidius avenae*), a few hyperparasitoid species and one aphid species (*Metopolophium dirhodum*). Modifications in minimal winter temperatures and frequency of frost events across the years partially explain observed community changes, although no clear climatic trend can be emphasized. Strong bottom-up effects from the relative abundance of aphid species also determine the relative abundance of parasitoid species each winter. Strong compartmentalization in parasitoid preference for host is reported. We suggest the recent modifications in parasitoid community composition to be linked to shifts in diapause expression (reduction or arrest of the use of winter diapause) and to host availability throughout the year. We highlight the implications for natural biological control in cereal fields. Perspectives are proposed to predict the composition of future host-parasitoid communities in the climate change context.

**Keywords:** overwintering strategies, diapause, species diversity, competition, biological control

## INTRODUCTION

Climate change impacts the geographic distribution, diversity, and abundance of organisms (Walther et al., 2002; Parmesan, 2006). In particular, climate warming can strongly influence their seasonal phenology, migration pattern, number of generations per year, and overwintering strategy (Roy and Sparks, 2000; Altermatt, 2010; Bale and Hayward, 2010). In temperate areas temperatures are increasing faster in winter than in summer, leading to overall milder, shorter, and later winter periods (IPCC, 2014). Plastic and adaptive responses of organisms to new thermal environments

**OPEN ACCESS**

**Edited by:**
Paul-André Calatayud,
Institut de Recherche Pour le
Développement (IRD), France

**Reviewed by:**
Pablo Almaraz,
Instituto de Ciencias Marinas de
Andalucla (ICMAN), Spain
Régis Babin,
Centre de Coopération Internationale
en Recherche Agronomique Pour le
Développement (CIRAD), France

*Correspondence:
Kévin Tougeron
tougeron.kevin@gmail.com

**Specialty section:**
This article was submitted to
Population and Evolutionary
Dynamics,
a section of the journal
Frontiers in Ecology and Evolution

**Received:** 20 July 2018
**Accepted:** 10 October 2018
**Published:** 29 October 2018

**Citation:**
Tougeron K, Damien M, Le Lann C, Brodeur J and van Baaren J (2018) Rapid Responses of Winter Aphid-Parasitoid Communities to Climate Warming. *Front. Ecol. Evol.* 6:173. doi: 10.3389/fevo.2018.00173

<footnote n="1">¹ Univ Rennes, CNRS, ECOBIO (Écosystèmes, Biodiversité, Évolution) - UMR 6553, Rennes, France, ² Institut de Recherche en Biologie Végétale, Département de Sciences Biologiques, Université de Montréal, Montreal, QC, Canada, ³ UMR 1355-7254 INRA-CNRS, Institut Sophia Agrobiotech, Université Nice Sophia Antipolis, Sophia Antipolis, France</footnote>
could modify species interactions such as competition, predation and parasitism and impact the structure and stability of communities (Hughes, 2000).

In the context of the global diversity crisis, studies increasingly focus on how trophic networks respond to global changes (Parmesan, 2006; Chaianunporn and Hovestadt, 2015). Indeed, species interactions within communities support the majority of ecosystem services and must be considered as study systems *per se* (Montoya et al., 2003). In some cases, food web structure, and composition are quite fragile and are likely to rapidly change in the context of climate warming; understanding how and why these food webs vary in space and time is a central objective in community ecology (Facey et al., 2014). New species appear while others disappear from food webs and changes in species interactions between trophic levels occur (Tylianakis et al., 2008; Chaianunporn and Hovestadt, 2015).

Parasites are omnipresent in almost every food web (Dobson et al., 2008) and their interactions with hosts greatly contribute to ecosystem functioning (Lafferty et al., 2008). Their ecology is tightly associated with their hosts' and likely to be influenced by climate change, threatening the provisioning of ecosystem services such as natural biological pest control by insect parasitoids (Hance et al., 2007; Jeffs and Lewis, 2013). The impacts of climate change or inter-annual variations in climatic conditions on host-parasitoid communities remain little explored compared to other types of food webs (e.g., plant-herbivore networks; Singer and Parmesan, 2010) and there have been few attempts at predicting their future structure and composition under different scenarios of climate change (Jeffs and Lewis, 2013).

In regions characterized by mild winter temperatures, the absence of lethal frosts allows aphids and their parasitoids to remain active and reproduce throughout winter. In cereal crops of Western France, aphid-parasitoid communities in winter were historically (over the past three decades) composed of the two parasitoid species *Aphidius rhopalosiphi* De Stefani-Perez and *Aphidius matricariae* Haliday and the two aphid species *Rhopalosiphum padi* (L.) and *Sitobion avenae* (Fabricius). From late spring to fall, additional species were present, including the parasitoids *Aphidius ervi* Haliday and *Aphidius avenae* Haliday and the aphid *Metopolophium dirhodum* (Walker) (Rabasse et al., 1983; Krespi, 1990; Krespi et al., 1997). These seasonal variations in aphid and parasitoid species occurrence seem to be consistent across Western Europe in cereal crops (Lumbierres et al., 2007; Honek et al., 2018) and likely reflect thermal niche separation (Le Lann et al., 2011; Andrade et al., 2016). The parasitoid *A. avenae* shows less cold resistance and more heat resistance than *A. rhopalosiphi* (Le Lann et al., 2011), while the aphid *R. padi* prefers cooler conditions and is more cold resistant than *S. avenae* (Jarošík et al., 2003; van Baaren et al., 2010; Alford et al., 2016).

Andrade et al. (2016) reported that in host-parasitoid communities of Western France parasitoid species usually not encountered during winter are now being active throughout the season and exploiting anholocyclic aphids (i.e., aphids that have parthenogenetic reproduction all-year-long). In the present study, we have adopted a community-wide approach to analyze the effects of long-term (inter-annual) variations in

climatic conditions on aphid and parasitoid species occurrence and relative abundance. Using a 9-year dataset, we explored community assembly rules by first describing temporal changes in winter aphid-parasitoid associations, and then linking these changes to modifications in winter climatic conditions and to other abiotic and biotic factors such as shifts in species interactions. We expect warmer winters to be associated with occurrence and higher abundance of *A. avenae* and *A. ervi* for parasitoids and *S. avenae* and *M. dirhodum* for aphids. We expect colder winters to be associated with the presence of *A. rhopalosiphi* and *A. matricariae* for parasitoids and *R. padi* for aphids.

# MATERIALS AND METHODS

## Data Collection

Data consist in aphid-parasitoid pairs of species gathered from different studies conducted in the Long Term Ecological Research (LTER) ZA Armorique, France (48°29' N–1°35'W), each winter from 2009/10 to 2017/18, at variable dates from late-November to mid-March of each year (excepted in 2009/10 when sampling was conducted only in January and February). Data from winter 2009/10 to winter 2012/13 were obtained from Andrade et al. (2016) and Eoche-Bosy et al. (2016), data of 2013/14 from Tougeron et al. (2016), data of 2014/15 from Tougeron et al. (2017), data of 2015/16 from Damien et al. (2017) and data of 2016/17 and 2017/18 from unpublished field results. In winter 2010/2011, no parasitoids nor aphids were found in the fields due to frost conditions during 15 consecutive days in November (Andrade et al., 2016) so this winter was excluded from the dataset to minimize unbalanced analyses on community data.

Mean, mean minimum and mean maximum daily temperature data per sampling year in the LTER were obtained from Météo France (2018). Additionally, we calculated the number of frost events (i.e., occurrence of at least 3 consecutive days with negative mean temperatures, which could be lethal for most species) and mean duration of frost events (days). Highly correlated variables >70% were not used for our analyses. We thus only used the mean minimal temperature (correlated with mean temperature and mean maximum temperatures) and the mean duration of frost events (correlated with the number of frost events) for statistical analyses. Mean temperature was used in graphic representations as a generic measure for interannual climatic variations. Raw data can be found in **Supplementary Table S1**.

## Sampling and Quantitative Food-Webs

In each of the studies from which data was collected, sampling was performed following the protocol of Andrade et al. (2016), with differences in the location of sampled fields due to crop rotations. In brief, sampling was conducted every 10 days in six to fifteen cereal fields each year; mainly winter wheat, but also barley and triticale. Sampling effort was consistent across the years. All aphids and aphid mummies (i.e., exoskeleton of dead aphid containing a developing parasitoid) were randomly collected during a 1-h period over an approximate area of 1,000 m². Aphid density was very low in winter; around one aphid/m².
Accordingly, parasitism rate during all winters was high (60-90%), underlying rarity of hosts and high competition levels among parasitoids. As the aphid-parasitoid network is stable over the winter sampling period (Andrade et al., 2016; Damien et al., 2017), data were pooled for the entire winter season. Live aphids were brought back to the laboratory and reared on winter wheat until mummification or death. All mummies were maintained in gelatin capsules at ambient temperature (17-20°C) until parasitoid emergence. Emerging adult primary parasitoids and aphid hosts were then identified to the species based on morphological characters (Hullé et al., 2006). Hyperparasitoids (secondary parasitoids) were identified to the genus level. Parasitoids and hyperparasitoids that emerged more than 25 days after sampling, representing each year <25% of the total number of sampled mummies, were excluded from the food-web analysis to avoid accounting for diapausing individuals when characterizing winter-active communities. Important differences in sample sizes are due to variations in aphid densities and different climatic conditions among years.

To examine trophic interactions between host and parasitoid species, quantitative food webs using the relative abundance (%) of each species were constructed for each winter following the methodology of Memmott et al. (1994). Hyperparasitoids were assigned to the same trophic level than primary parasitoids because it was impossible to assess in which parasitoid host species they developed.

## Analyses

Food webs were compared among years using several quantitative and qualitative metrics calculated using the bipartite (Dormann et al., 2009) and the codyn R packages (Hallett et al., 2016): Connectance—the overall complexity of the food web (realized proportion of potential links); Web Asymmetry—the balance between numbers of parasitoid and aphid species (negative values indicate more species in higher than in lower trophic-level); H2—the level of specialization within a network, from 0 (no specialization) to 1 (perfect specialization); Generality—the weighted mean number of aphid species exploited by each parasitoid species; Vulnerability—the weighted mean number of parasitoid species attacking a given aphid species.

Then, we performed a non-metric multidimensional scaling (NMDS) analysis to group years by climatic similarities based on a distance matrix; accordingly, years were grouped by four based on their distance to the more extreme years on the NMDS. Following these analyses, years were characterized as either mild or cold winters, and these categories were used to group years on the following Principal Component Analysis (PCA) representation. A first PCA was performed to separate each year of sampling based on selected climatic variables; mean minimal temperatures, and mean duration of frost events. Another PCA was performed to separate each year of sampling based on aphid and parasitoid species relative abundances. Finally, a Canonical Correspondence Analysis (CCA) was performed to assess relationships between aphid and parasitoid species and the climatic variables matrixes. ANOVA-like permutation tests for CCA (vegan) were used to assess the significance of

constraints. Only primary parasitoids from the genus *Aphidius* were considered for analyses on climatic variations.

As we wanted to account for species co-occurrences in our analyses, species-by-species models were not appropriate. Instead, we used a community approach by analyzing separately the effects of the selected climatic variables across years on parasitoid species (matrix containing the relative abundance of the four parasitoid species) and aphid species (matrix containing the relative abundance of the three aphid species). To do so, Bray-Curtis dissimilarity indexes in species relative abundances of each sampling year were calculated (separately for parasitoids and aphids) and fitted to linear models as response variables using the Adonis-Permanova function from the R package *vegan* (Oksanen et al., 2015) calculating permutation test with pseudo-F ratios. For aphids, the minimal temperature and the mean duration of frost events were used as explanatory variables. For parasitoids, we also included the relative abundance of each of the three aphid species as explanatory factors in the linear models. All analyses were performed using R software (R Core Team, 2017).

## RESULTS

Changes in species richness and relative abundances were observed from winter 2009/10 to winter 2017/18 in the aphid-parasitoid food webs, with important inter-annual variations (Figure 1A). The year 2009/10 was similar to the past three decades, as described in introduction, with *A. rhopalosiphii* and *A. matricariae* being the only two parasitoid species active in winter and exploiting *S. avenae* and *R. padi*. In addition to these two aphid species, the aphid *M. dirhodum* was present every winter starting 2011/12 at 12%, and it represented up to 68% of the aphid species relative abundance in winter 2015/16 (Figure 1A). There was high variability in aphid proportions between years, for each species.

*A. avenae* was observed for the first time in the winter 2011/12 with a relative abundance of 52%. *A. ervi* was observed in the network in 2013/14 with a relative abundance of 5%. Both species have since remained present in the network, although at variable relative abundances. *A. rhopalosiphii* was present every winter while *A. matricariae* occurrence and relative abundance were highly variable over the years. *Ephedrus plagiator* (Nees) and *Diaeretiella rapae* (M'intosh), two generalist species (Hullé et al., 2006), were anecdotally reported in winter 2011/12 and 2012/13. Hyperparasitoids from the genera *Alloxysta*, *Asaphes* and *Phaenoghyphis* were also detected in two out of eight winters.

Mean temperature varied across years and ranged from 4.2°C in 2009/10, to 7.9°C in 2015/16, and was graphically concurrent with variations in aphid and parasitoid relative abundances (Figure 1B).

Food-web metrics are summarized in Table 1. The food-web connectance was ≥0.6 every year; all links the first year and almost all potential links the following years between each parasitoid and aphid species were observed. In most years, the food web was asymmetric, with more parasitoid species than aphid species. The degree of specialization within the food-web (H2 index) tended to decrease over the years. Accordingly, the
<figure><img src="image_3.png" /></figure>

**FIGURE 1** | (A) Quantitative food webs of parasitoid and aphid community composition in winters 2009/10 to 2017/18 (there was no insects in 2010/11). Upper and lower bars represent the relative abundance (%) of parasitoid, including hyperparasitoids, and aphid species, respectively. The thickness of the arrow between bars is proportional to the interaction strength between a pair of species. Total number of individuals (N) used to construct each food web is shown for each year.

(B) Rankplots showing the relative abundances (%) of aphids (left panel) and *Aphidius* parasitoids (right panel) each winter from 2009/10 to 2017/18. Mean winter temperature (°C) is shown in black.
<table><caption>**TABLE 1** | Number of species in each trophic level (primary parasitoids and aphids) and food-web metrics for each sampled winter.</caption><thead><tr><th>Year</th><th>2009/10</th><th>2011/12</th><th>2012/13</th><th>2013/14</th><th>2014/15</th><th>2015/16</th><th>2016/17</th><th>2017/18</th></tr></thead><tbody><tr><td>Number of primary parasitoid species</td><td>2</td><td>3</td><td>5</td><td>4</td><td>4</td><td>4</td><td>3</td><td>4</td></tr><tr><td>Number of aphid species</td><td>2</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td><td>3</td></tr><tr><td>Connectance</td><td>1.00</td><td>0.89</td><td>0.60</td><td>0.83</td><td>0.92</td><td>0.92</td><td>0.89</td><td>0.83</td></tr><tr><td>Web asymmetry</td><td>0.00</td><td>0.00</td><td>-0.25</td><td>-0.14</td><td>-0.14</td><td>-0.14</td><td>0.00</td><td>-0.14</td></tr><tr><td>H2</td><td>0.53</td><td>0.71</td><td>0.41</td><td>0.36</td><td>0.24</td><td>0.22</td><td>0.34</td><td>0.24</td></tr><tr><td>Generality</td><td>1.63</td><td>1.43</td><td>2.18</td><td>2.41</td><td>2.67</td><td>2.26</td><td>2.25</td><td>2.50</td></tr><tr><td>Vulnerability</td><td>1.77</td><td>2.00</td><td>1.96</td><td>2.47</td><td>2.41</td><td>2.73</td><td>2.26</td><td>2.64</td></tr></tbody></table>

Connectance, the overall complexity of the food web (realized proportion of potential links); Web Asymmetry, the balance between numbers of parasitoid and aphid species (negative values indicate more species in higher than in lower trophic-level); H2, the level of specialization within a network, from 0 (no specialization) to 1 (perfect specialization); Generality, the weighted mean number of aphid species exploited by each parasitoid species; Vulnerability, the weighted mean number of parasitoid species attacking a given aphid species.

generality and vulnerability indexes tend to increase over time; each parasitoid species attacked more aphids and each aphid species was exploited by more parasitoid species over the years, in general (Table 1).

NMDS analysis showed that years 2009/10, 2012/13, 2013/14, and 2017/18 can be grouped together as “cold winters” whereas years 2011/12, 2014/15, 2015/16, and 2016/17 can be described as “mild winters.” This clustering was supported by 94.7% inertia on the PC1 of the PCA (Figure 2A). Graphically, coldest winters (i.e., decreasing minimal temperatures and increasing duration of frost events) were overall associated with the co-occurrence and higher relative abundances of both *A. rhopalosiphi* and *A. matricariae* parasitoids and *R. padi* aphids, while warmest winters were associated with co-occurrence and higher abundance of both *A. avenae* and *A. ervi* parasitoids and *S. avenae* aphids. *M. dirhodum* was highly abundant in the warmest winter (2015/16) (Figures 1, 2B). However, the aphid-parasitoid community PCA was only supported by 43% inertia on the PC1, indicating that species partition on this figure is only partially explained by the sampling year (Figure 2B).

We found a marginally non-significant influence of the selected climatic data (mean minimal temperature and duration of frost events) on global aphid and parasitoid relative abundances across years (CCA ANOVA-like permutation test, $F = 1.67$, df = 2, $p = 0.06$). In details, aphid abundances were significantly affected by changes in mean minimal temperatures (Permanova, $F = 2.03$, df = 1, $R^2 = 0.22$, $p = 0.03$) and in mean duration of frost events across the years ($F = 2.57$, df = 1, $R^2 = 0.27$, $p = 0.04$), in the way described in the precedent paragraph. Parasitoid abundances were marginally significantly affected by changes in mean minimal temperatures ($F = 2.67$, df = 1, $R^2 = 0.20$, $p = 0.05$) and significantly affected by changes in mean duration of frosts across the years ($F = 14.72$, df = 1, $R^2 = 0.15$, $p = 0.01$), as described above. Changes in parasitoid abundances across the years were not affected by the abundances of the aphid *S. avenae* ($F = 1.7$, df = 1, $R^2 = 0.16$, $p = 0.16$) but was significantly affected by the abundances of *M. dirhodum* ($F = 35.5$, df = 1, $R^2 = 0.37$, $p=0.01$) and *R. padi* ($F = 17.2$, df = 1, $R^2 = 0.18$, $p = 0.009$).

Preferential associations between aphid hosts and primary parasitoid species occurred as shown by the arrows on Figure 2B (pooled data across years). The aphid *S. avenae* was mostly

associated with the parasitoids *A. avenae* and *A. ervi* while the aphid *R. padi* was mostly associated with the parasitoids *A. rhopalosiphi* and *A. matricariae*. *M. dirhodum* was not preferentially associated with any parasitoid in the food web (Figure 2B).

## DISCUSSION

Our results illustrate how climatic changes during winter have rapidly, over nine years, translated into modifications in species composition within an aphid-parasitoid community. By contrast with the community typically described over the past decades (Rabasse et al., 1983; Krespi et al., 1997), two parasitoid species *A. avenae* and to a minor extent *A. ervi*, and one aphid species *M. dirhodum* are now active in Western France throughout winter together with other species of the community. The winter trophic network composition in cereal fields is getting similar to what is usually described in spring in this area. The winter food web has become more diversified in aphid and parasitoid species and, while the connectance (realized links) remains stable over time, the degree of specialization tends to decrease, suggesting that parasitoids exploit aphids in function of their relative abundance, as reported in spring (Andrade et al., 2016). This may be due to increasing aphid densities in winter on cereal crops, leading to lower competition pressure among parasitoids.

Changes in occurrence do not arise from recent modifications in distribution range of the species, since *A. avenae*, *A. ervi*, and *M. dirhodum* have been commonly observed for several decades in spring at the same location (Krespi, 1990; Andrade et al., 2016). Our results suggest a recent shift in overwintering strategy in *A. avenae* and *A. ervi* parasitoid populations with some individuals remaining active throughout the winter rather than entering diapause. This hypothesis is supported by results from a laboratory experiment showing that diapause incidence in both parasitoid species was low (<15%), even when parasitoids were reared under fall-like temperature conditions that usually induce high levels of diapause (Tougeron et al., 2017). Variations in species composition in the food web over the years may arise from differences in thermal niches; the most cold-resistant species usually remained active during winter (e.g., *A. rhopalosiphi*, *A. matricariae*, and *R. padi*) whereas less
<figure><img src="image_4.png" /></figure>

cold-resistant species (e.g., *A. avenae*, *A. ervi*, *M. dirhodum*, and hyperparasitoids) were probably in diapause and mostly active from spring to fall (Krespi, 1990; Le Lann et al., 2011; Alford et al., 2016; Andrade et al., 2016; Tougeron et al., 2017). Overwintering temperature may now be warm enough to allow niche overlapping of all these species during winter. No clear trends in temperature change seem to appear over 9 years of study so our data can only lead to partial conclusion on a climate-change effect on aphids and parasitoids. However, we recently showed that from 1976 to 2015, the daily average temperature increased by an average of 1.13 °C in winter. Furthermore, days with cold spell events during winter decreased in frequency since the 1970s (Tougeron et al., 2017). Thus, it is important to state that the interannual variation in climatic conditions observed over 9 years may not be larger than the long-term trend in the study area.

It has been shown that fine-scale intra-seasonal temperature variations (i.e., temperature experienced by the insect during its development) played an important role in shaping local aphid-parasitoid communities in Western France (Andrade et al., 2015, 2016). For example higher developmental temperatures were associated with increasing abundance in *A. avenae*, *S. avenae* and decreasing abundance in *R. padi* (Andrade et al., 2016). Winter 2016/17 was on average warmer than other winters but important cold spells occurred in December and January, which may have conducted to higher abundances of *R. padi* and *A. rhopalosiphi*, the more cold tolerant species in the food-web (Le Lann et al., 2011; Alford et al., 2016), and quasi-disappearance of

*M. dirhodum* from the system, through environmental filtering. Such thermal extremes and microclimatic variations may reduce or eliminate any advantages of global warming for some species (Ma et al., 2015; Sgrò et al., 2016) and may impede evaluation and prediction of climate change effects on community dynamics (Bailey and van de Pol, 2016; Blonder et al., 2017).

We have shown that both mean minimal temperatures and mean duration of frost events over the winter are predictors of winter aphid abundances and of their variation in occurrence among years. Honek et al. (2018) also demonstrated that temperature in winter was an important predictor of maximum abundances of cereal aphids during the weeks following sampling. However, change in mean minimal temperature and duration of frost events only slightly contributed to the trend observed in parasitoid relative abundance changes over the years. Stochastic effects or other environmental variables than temperature such as host-parasitoid interactions may better explain inter-annual variations in species abundances and occurrences during winter. For instance, we showed high level of host-parasitoid compartmentalization within the food web; the variation in relative abundance of some species was highly correlated with abundance of other species, suggesting bottom-up effects on parasitoid abundance. The importance of host species in shaping parasitoid response to climate warming must therefore be accounted (Barton and Ives, 2014). Parasitoids and their hosts may also be influenced by microclimatic refuges in the landscape (Tougeron et al., 2016; Alford et al., 2017), by surrounding plant covers (Gagic et al.,
2012; Damien et al., 2017) or by plant quality (Honek et al., 2018).

Modifications of the parasitoid guild could also be due to shifts in competition for hosts following the addition of new species. Indeed, female parasitoids show seasonal variations in foraging behavior (Roitberg et al., 1992) and can adapt their foraging strategies to competition or host-patch quality (Outreman et al., 2005; Le Lann et al., 2008; Barrette et al., 2010; Moiroux et al., 2015). In winter, it has been demonstrated that female parasitoids adopt generalist strategies due to shortage of optimal hosts, leading to high competition, whereas spring parasitoids usually display specialist strategies by selecting optimal host species (Eoche-Bosy et al., 2016). The recent addition of *A. avenae* and *A. ervi* in the overwintering food web, which are good competitors at exploiting *S. avenae*, may have reduced the abundance of *A. matricariae* and *A. rhopalosiphii* (Le Lann et al., 2012; Andrade et al., 2016; Eoche-Bosy et al., 2016).

Climate warming challenges coexistence and interaction between ecologically related species, as well as community stability and ecosystem functioning (van der Putten et al., 2004; Tylianakis et al., 2008). In cereal fields, overwintering reproduction of aphid parasitoids plays an important role in suppressing early cereal aphid populations in spring (Langer et al., 1997; Plantegenest et al., 2001; Honek et al., 2018). Increasing number of parasitoid species during winter due to climate warming could enhance aphid natural biological control through increasing already high parasitism rate in winter, even if consequences of niche overlapping between parasitoid species via addition of new species in the food-web are difficult to predict. The presence of non-diapausing hyperparasitoids, reported for the first time in winter in Western France in 2012/13 (Tougeron et al., 2017), may reduce the efficiency of biological control in the fields, although hyperparasitoids can sometimes stabilize primary parasitoid populations (Tougeron and Tena, in press). In Spain, characterized by relatively warm winters, hyperparasitism remains high throughout the year which disrupts biological control by primary parasitoids in orchards (Gómez-Marco et al., 2015). With an expected temperature increase from 0.5 to 2°C in the next decades (Karl and Trenberth, 2003), occurrences of new species in food-webs such as hyperparasitoids can be more common, either through shifts in geographic distribution (e.g., biological invasions) or shifts in phenology (e.g., reduction of diapause expression) (Tougeron and Tena, in press).

By drawing from the effects of interannual climatic condition variations over nine years, and based upon studies conducted in this area decades ago, we can observe that changes in host-parasitoid winter communities are not random, which ultimately helps at determining potential actions of long-term climate-change. Based on the data currently available on this host-parasitoid system, we observed the homogenization of winter and spring aphid-parasitoid communities, most likely due to change in phenology after increasing temperatures during winter and decreasing duration of frost events. High variations between years may underline a transition period between two episodes of stable communities, although the food-web could also remain unstable due to variations in climate at fine temporal scale. Data over a longer time period will help at establishing stronger conclusions on climate-change effects on insect communities in this area. Predictive analyses on the community structures should now integrate local changes in overwintering strategies of one or more species to identify the potential effect of climate change on ecosystem services provided by parasitoids.

## AUTHOR CONTRIBUTIONS

KT and MD collected and analyzed the data. KT wrote a first version of the manuscript. KT, MD, CL, JB, and JvB made substantial contributions to the manuscript.

## ACKNOWLEDGMENTS

We thank all collaborators at Ecobio lab who helped with field work, S. Chollet and J.-S. Pierre for their assistance in analyzing the community data. KT was funded by the French Region Bretagne and the Canada Research Chair in biological control awarded to JB. MD was funded by the FLEUR project and INRA. This study was supported by the LTER France Zone Atelier Armorique.

## SUPPLEMENTARY MATERIAL

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fevo.2018.00173/full#supplementary-material

**Supplementary Table S1** | Data on species relative abundance and climatic variations across the years.

## REFERENCES

Alford, L., Androdias, A., Franco, T., Pierre, J.-S., Burel, F., and van Baaren, J. (2016). Potential host manipulation by the aphid parasitoid *Aphidius avenae* to enhance cold tolerance. *PLoS ONE* 11:e0168693. doi: 10.1371/journal.pone.0168693

Alford, L., Tougeron, K., Pierre, J.-S., Burel, F., and Van Baaren, J. (2017). The effect of landscape complexity and microclimate on the thermal tolerance of a pest insect. *Insect Sci.* 25, 905–915. doi: 10.1111/1744-7917.12460

Altermatt, F. (2010). Climatic warming increases voltinism in European butterflies and moths. *Proc. R. Soc B* 277, 1281–1287. doi: 10.1098/rspb.2009.1910

Andrade, T. O., Krespi, L., Bonnardot, V., van Baaren, J., and Outreman, Y. (2016). Impact of change in winter strategy of one parasitoid species on

the diversity and function of a guild of parasitoids. *Oecologia* 180, 877–888. doi: 10.1007/s00442-015-3502-4

Andrade, T. O., Outreman, Y., Krespi, L., Plantegenest, M., Vialatte, A., Gauffre, B., et al. (2015). Spatiotemporal variations in aphid-parasitoid relative abundance patterns and food webs in agricultural ecosystems. *Ecosphere* 6, 1–14. doi: 10.1890/ES15-00010.1

Bailey, L. D., and van de Pol, M. (2016). Tackling extremes: challenges for ecological and evolutionary research on extreme climatic events. *J. Anim. Ecol.* 85, 85–96. doi: 10.1111/1365-2656.12451

Bale, J. S., and Hayward, S. A. L. (2010). Insect overwintering in a changing climate. *J. Exp. Biol.* 213, 980–994. doi: 10.1242/jeb.037911

Barrette, M., Boivin, G., Brodeur, J., and Giraldeau, L.-A. (2010). Travel time affects optimal diets in depleting patches. *Behav. Ecol. Sociobiol.* 64, 593–598. doi: 10.1007/s00265-009-0876-5
Barton, B. T., and Ives, A. R. (2014). Direct and indirect effects of warming on aphids, their predators, and ant mutualists. *Ecology* 95, 1479–1484. doi: 10.1890/13-1977.1

Blonder, B., Moulton, D. E., Blois, J., Enquist, B. J., Graae, B. J., Macias-Fauria, M., et al. (2017). Predictability in community dynamics. *Ecol. Lett.* 20, 293–306. doi: 10.1111/ele.12736

Chaianunporn, T., and Hovestadt, T. (2015). Evolutionary responses to climate change in parasitic systems. *Glob. Chang. Biol.* 21, 2905–2916. doi: 10.1111/gcb.12944

Damien, M., Le Lann, C., Desneux, N., Alford, L., Al Hassan, D., Georges, R., et al. (2017). Flowering crops in winter increases pest control but not trophic link diversity. *Agric. Ecosyst. Environ.* 247, 418–425. doi: 10.1016/j.agee.2017.07.015

Dobson, A., Lafferty, K. D., Kuris, A. M., Hechinger, R. F., and Jetz, W. (2008). Homage to linnaeus: how many parasites? How many hosts? *Proc. Natl. Acad. Sci. USA.* 105, 11482–11489. doi: 10.1073/pnas.0803232105

Dormann, C. F., Fründ, J., Blüthgen, N., and Gruber, B. (2009). Indices, graphs and null models: analyzing bipartite ecological networks. *Open Ecol. J.* 2, 7–24. doi: 10.2174/1874213000902010007

Eoche-Bosy, D., Outreman, Y., Oliveira Andrade, T., Krespi, L., and Baaren, J. (2016). Seasonal variations of host resources influence foraging strategy in parasitoids. *Entomol. Exp. Appl.* 161, 11–19. doi: 10.1111/eea.12494

Facey, S. L., Ellsworth, D. S., Staley, J. T., Wright, D. J., and Johnson, S. N. (2014). Upsetting the order: how climate and atmospheric change affects herbivore–enemy interactions. *Curr. Opin. Insect Sci.* 5, 66–74. doi: 10.1016/j.cois.2014.09.015

Gagic, V., Hänke, S., Thies, C., Scherber, C., Tomanović, Ž., and Tscharntke, T. (2012). Agricultural intensification and cereal aphid–parasitoid–hyperparasitoid food webs: network complexity, temporal variability and parasitism rates. *Oecologia* 170, 1099–1109. doi: 10.1007/s00442-012-2366-0

Gómez-Marco, F., Urbaneja, A., Jaques, J. A., Rugman-Jones, P. F., Stouthamer, R., and Tena, A. (2015). Untangling the aphid-parasitoid food web in citrus: can hyperparasitoids disrupt biological control? *Biol. Control* 81, 111–121. doi: 10.1016/j.biocontrol.2014.11.015

Hallett, L. M., Jones, S. K., MacDonald, A. A. A., Jones, M. B., Flynn, D. F. B., Ripplinger, J., et al. (2016). Codyn: an R package of community dynamics metrics. *Methods Ecol. Evol.* 7, 1146–1151. doi: 10.1111/2041-210X.12569

Hance, T., van Baaren, J., Vernon, P., and Boivin, G. (2007). Impact of extreme temperatures on parasitoids in a climate change perspective. *Annu. Rev. Entomol.* 52, 107–126. doi: 10.1146/annurev.ento.52.110405.091333

Honek, A., Martinkova, Z., Saska, P., and Dixon, A. F. G. (2018). Aphids (Homoptera: Aphididae) on winter wheat: predicting maximum abundance of *Metopolophium dirhodum*. *J. Econ. Entomol.* 111, 1751–1759. doi: 10.1093/jee/toy157

Hughes, G. E. (2000). Biological consequences of global warming: is the signal already apparent? *Trends Ecol. Evol.* 15, 56–61. doi: 10.1016/S0169-5347(99)01764-4

Hullé, M., Turpeau, E., and Chaubet, B. (2006). *Encyclop'aphid, a key for Aphids and Their Parasitoids*. INRA (French Institute for Agronomical Research).

IPCC (2014). "Climate change 2013 the physical science basis," in *Final Draft Underlying Scientific-Technical Assessment: Working Group I Contribution to the IPCC Fifth Assessment Report*, eds T.F. Stocker, D. Qin, G.-K. Plattner, M. Tignor, S.K. Allen, J. Boschung, A. Nauels, Y. Xia, V. Bex, and P. M. Midgley (Cambridge, UK; New York, NY: Cambridge University Press), 1535.

Jarošík, V., Honěk, A., and Tichopad, A. (2003). Comparison of field population growths of three cereal aphid species on winter wheat. *Plant Protect Sci.* 39, 61–64. doi: 10.17221/3827-PPS

Jeffs, C. T., and Lewis, O. T. (2013). Effects of climate warming on host-parasitoid interactions: effects of climate warming. *Ecol. Entomol.* 38, 209–218. doi: 10.1111/een.12026

Karl, T. R., and Trenberth, K. E. (2003). Modern global climate change. *Science* 302, 1719–1723. doi: 10.1126/science.1090228

Krespi, L. (1990). *Etude de la Biocénose Parasitaire des Pucerons des Céréales Dans le Bassin de Rennes : cas Particulier d'Aphidius Uzbekistanicus Luz*. Ph.D. dissertation, Université de Rennes 1, France.

Krespi, L., Dedryver, C.-A., Creach, V., Rabasse, J.-M., Le Ralec, A., and Nenon, J.-P. (1997). Variability in the development of cereal aphid parasitoids and hyperparasitoids in oceanic regions as a response to climate and abundance of hosts. *Popul. Ecol.* 26, 545–551.

Lafferty, K. D., Allesina, S., Arim, M., Briggs, C. J., De Leo, G., Dobson, A. P., et al. (2008). Parasites in food webs: the ultimate missing links. *Ecol. Lett.* 11, 533–546. doi: 10.1111/j.1461-0248.2008.01174.x

Langer, A., Stilmant, D., Verbois, D., and Hance, T. (1997). Seasonal activity and distribution of cereal aphid parasitoids in Belgium. *Entomophaga* 42, 185–191. doi: 10.1007/BF02769896

Le Lann, C., Outreman, Y., Van Alphen, J. J. M., Krespi, L., Pierre, J.-S., and Van Baaren, J. (2008). Do past experience and competitive ability influence foraging strategies of parasitoids under interspecific competition? *Ecol. Entomol.* 33, 691–700. doi: 10.1111/j.1365-2311.2008.01017.x

Le Lann, C., Roux, O., Serain, N., Van Alphen, J. J. M., Vernon, P., and Van Baaren, J. (2011). Thermal tolerance of sympatric hymenopteran parasitoid species: does it match seasonal activity? *Physiol. Entomol.* 36, 21–28. doi: 10.1111/j.1365-3032.2010.00758.x

Le Lann, C., Visser, B., van Baaren, J., van Alphen, J. J. M., and Ellers, J. (2012). Comparing resource exploitation and allocation of two closely related aphid parasitoids sharing the same host. *Evol. Ecol.* 26, 79–94. doi: 10.1007/s10682-011-9498-2

Lumbierres, B., Starý, P., and Pons, X. (2007). Seasonal parasitism of cereal aphids in a Mediterranean arable crop system. *J. Pest Sci.* 80, 125–130. doi: 10.1007/s10340-006-0159-0

Ma, G., Hoffmann, A. A., and Ma, C.-S. (2015). Daily temperature extremes play an important role in predicting thermal effects. *J. Exp. Biol.* 218, 2289–2296. doi: 10.1242/jeb.122127

Memmott, J., Godfray, H., and Gauld, I. D. (1994). The structure of a tropical host-parasitoid community. *J. Anim. Ecol.* 63, 521–540. doi: 10.2307/5219

Moiroux, J., Boivin, G., and Brodeur, J. (2015). Temperature influences host instar selection in an aphid parasitoid: support for the relative fitness rule. *Biol. J. Linnean Soc.* 115, 792–801. doi: 10.1111/bij.12545

Montoya, J. M., Rodríguez, M. A., and Hawkins, B. A. (2003). Food web complexity and higher-level ecosystem services. *Ecol. Lett.* 6, 587–593. doi: 10.1046/j.1461-0248.2003.00469.x

Oksanen, J., Blanchet, F. G., Kindt, R., Legendre, P., Minchin, P. R., O'Hara, R. B., et al. (2015). vegan: Community Ecology Package. R package version 2.3-0. Available online at: http://CRAN.R-project.org/package=vegan

Outreman, Y., Le Ralec, A., Wajnberg, E., and Pierre, J.-S. (2005). Effects of within- and among-patch experiences on the patch-leaving decision rules in an insect parasitoid. *Behav. Ecol. Sociobiol.* 58, 208–217. doi: 10.1007/s00265-004-0895-1

Parmesan, C. (2006). Ecological and evolutionary responses to recent climate change. *Annu. Rev. Ecol. Evol. Syst.* 37, 637–669. doi: 10.1146/annurev.ecolsys.37.091305.110100

Plantegenest, M., Pierre, J. S., Dedryver, C. A., and Kindlmann, P. (2001). Assessment of the relative impact of different natural enemies on population dynamics of the grain aphid *Sitobion avenae* in the field. *Ecol. Entomol.* 26, 404–410. doi: 10.1046/j.1365-2311.2001.00330.x,

R Core Team (2017). *R: A language and Environment For Statistical Computing*. R Foundation for Statistical Computing, Vienna.

Rabasse, J.-M., Dedryver, C.-A., Gelle, A., Lafont, J.-P., and Molinari, I. (1983). Biologie des pucerons des céréales dans l'Ouest de la France. III. Action des hyménoptères parasites sur les populations de *Sitobion avenae* F., *Metopolophium dirhodum* Wlk. et Rhopalosiphum padi L. Agronomie 3, 779–790. doi: 10.1051/agro:19830809

Roitberg, B. D., Mangel, M., Lalonde, R. G., Roitberg, C. A., van Alphen, J. J. M., and Vet, L. (1992). Seasonal dynamic shifts in patch exploitation by parasitic wasps. *Behav. Ecol.* 3, 156–165. doi: 10.1093/beheco/3.2.156

Roy, D. B., and Sparks, T. H. (2000). Phenology of British butterflies and climate change. *Glob. Chang. Biol.* 6, 407–416. doi: 10.1046/j.1365-2486.2000.00322.x

Sgrò, C. M., Terblanche, J. S., and Hoffmann, A. A. (2016). What can plasticity contribute to insect responses to climate change? *Annu. Rev. Entomol.* 61, 433–451. doi: 10.1146/annurev-ento-010715-023859

Singer, M. C., and Parmesan, C. (2010). Phenological asynchrony between herbivorous insects and their hosts: signal of climate change or pre-existing adaptive strategy? *Philos. Trans. R. Soc. B* 365, 3161–3176. doi: 10.1098/rstb.2010.0144
Tougeron, K., Le Lann, C., Brodeur, J., and van Baaren, J. (2017). Are aphid parasitoids from mild winter climates losing their winter diapause? *Oecologia* 183, 619–629. doi: 10.1007/s00442-016-3770-7

Tougeron, K., and Tena, A. (in press). Hyperparasitoids as new targets in biological control in a global change context. *Biological Control*. doi: 10.1016/j.biocontrol.2018.09.003

Tougeron, K., van Baaren, J., Burel, F., and Alford, L. (2016). Comparing thermal tolerance across contrasting landscapes: first steps towards understanding how landscape management could modify ectotherm thermal tolerance. *Insect Conserv. Divers.* 9, 171–180. doi: 10.1111/icad.12153

Tylianakis, J. M., Didham, R. K., Bascompte, J., and Wardle, D. A. (2008). Global change and species interactions in terrestrial ecosystems. *Ecol. Lett.* 11, 1351–1363. doi: 10.1111/j.1461-0248.2008.01250.x

van Baaren, J., Le Lann, C., and and, J. M., van Alphen, J. (2010). "Consequences of Climate Change for Aphid-Based Multi-trophic Systems," in *Aphid Biodiversity Under Environmental Change*, eds P. Kindlmann, A. F. G. Dixon, and J. P. Michaud (Dordrecht: Springer Netherlands), 55–68.

van der Putten, W. H., de Ruiter, P. C., Martijn Bezemer, T., Harvey, J. A., Wassen, M., and Wolters, V. (2004). Trophic interactions in a changing world. *Basic Appl. Ecol.* 5, 487–494. doi: 10.1016/j.baae.2004.09.003

Walther, G.-R., Post, E., Convey, P., Menzel, A., Parmesan, C., Beebee, T. J., et al. (2002). Ecological responses to recent climate change. *Nature* 416, 389–395. doi: 10.1038/416389a

**Conflict of Interest Statement:** The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Copyright © 2018 Tougeron, Damien, Le Lann, Brodeur and van Baaren. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.