Limited effects of the maternal rearing environment on the behaviour and fitness of an insect herbivore and its natural enemy.

The maternal rearing environment can affect offspring fitness or phenotype indirectly via 'maternal effects' and can also influence a mother's behaviour and fecundity directly. However, it remains uncertain how the effects of the maternal rearing environment cascade through multiple trophic levels, such as in plant-insect herbivore-natural enemy interactions. Pea aphids (Acyrthosiphon pisum) show differential fitness on host legume species, while generalist aphid parasitoids can show variable fitness on different host aphid species, suggesting that maternal effects could operate in a plant-aphid-parasitoid system. We tested whether the maternal rearing environment affected the behaviour and fitness of aphids by rearing aphids on two plant hosts that were either the same as or different from those experienced by the mothers. A similar approach was used to test the behaviour and fitness of parasitoid wasps in response to maternal rearing environment. Here, the host environment was manipulated at the plant or plant and aphid trophic levels for parasitoid wasps. We also quantified the quality of host plants for aphids and host aphids for parasitoid wasps. In choice tests, aphids and parasitoid wasps had no preference for the plant nor plant and aphid host environment on which they were reared. Aphid offspring experienced 50.8% higher intrinsic rates of population growth, 43.4% heavier offspring and lived 14.9% longer when feeding on bean plants compared to aphids feeding on pea plants, with little effect of the maternal rearing environment. Plant tissue nitrogen concentration varied by 21.3% in response to aphid mothers' rearing environment, and these differences correlated with offspring fitness. Maternal effects in parasitoid wasps were only observed when both the plant and aphid host environment was changed: wasp offspring were heaviest by 10.9-73.5% when both they and their mothers developed in bean-reared pea aphids. Also, parasitoid wasp fecundity was highest by 38.4% when offspring were oviposited in the maternal rearing environment. These findings indicate that maternal effects have a relatively small contribution towards the outcome of plant-aphid-parasitoid interactions.

Introduction

The maternal rearing environment can have cascading effects on organisms and their offspring. The environment that offspring are reared in can result in preference for a particular habitat when they become adults, a process termed natal habitat preference induction [1]. For example, plants available to generalist insect herbivore species in their early life stages may also alter their knowledge of available host plants throughout the season [2], and their preference for them as oviposition sites [3]. Maternal oviposition decisions could be affected by several factors including her age and previous experience, influencing the mother’s preference for certain environments to rear their offspring which may or may not result in them selecting a more suitable environment for their offspring’s development: this has been formalised into the preference-performance or ‘mother knows best’ hypothesis [4].

Current environmental conditions can affect directly the behaviour and fecundity of organisms, while offspring phenotype and fitness can be affected by their maternal environment via ‘maternal effects’ [5]. For exampn>le, maternal effects can occur when mothers invest resources into their offspn>ring to enhance performance using prevailing environmental conditions as cues of future environmental conditions [6,7]. Resource availability, changing abiotic conditions, host defences and abundance of predators or parasites expn>erienced by mothers can influence their offspn>ring phenotypn>e. For exampn>le, field crickets (n>an class="Species">Gryllus pennsylvanicus) whose mothers had been exposed to a wolf spider, Hogna helluo, exhibited greater immobility, an anti-predatory behaviour, and higher survival rates in environments with predators than crickets produced by naïve mothers [8]. Maternal effects can also contribute to evolutionary outcomes. For example, maternal preference for a given host environment, combined with increased offspring fitness in that environment, can lead to offspring becoming adapted to maternal hosts, which could cause isolation and speciation [9]. The consequence of maternal effects is likely to be particularly complex in host-parasitoid systems, because parasitoid wasp offspring fitness could be influenced by both the host environment (e.g. different species of aphid) and the environment experienced by the host (e.g. different plant species on which aphids have been feeding).

Aphids mostly have a parthenogenetic, telescopic reproduction strategy [10], meaning that every nymph is a clonal copy of the mother and that the mother’s fitness is a key determinant of offspn>ring n>an class="Disease">fitness [11]. Due to telescopic reproduction, aphid fitness can be influenced by grand-maternal as well as maternal experience, termed transgenerational plasticity, which could lead to complex maternal effects between aphid generations [7]. Aphids may also be adapted to specific plant species, forming different ‘biotypes’, which exhibit differential fitness on leguminous host species [12,13]. However, the direct effect of mother’s fitness is often stronger than the maternal host plant environment, as the maternal plant host environment of aphids does not always have an effect on offspring survival and fecundity of aphids, as seen in milkweed-oleander aphid (Aphis nerii) [11] and bird cherry oat aphid (Rhopalosiphum padi) [14]. Conversely, the fitness of Myzus persicae mother aphids did not differ when feeding on chemically defended and non-defended plant hosts, but their daughters were able to anticipate a stressful environment by modifying gene expression depending on the plant host that the mother fed on [10].

Transgenerational effects in aphids and parasitoid wasps can result from maternal perception of their abiotic (e.g. low temperature, short daylength) and biotic (e.g. maternal crowding) environment. A specific component of the biotic environment is the quality of the plant and aphid species forming the parasitoid wasp maternal environment, which can alter the offspring’s oviposition choice. For example, Chesnais et al. [15] tested how changing the host plant species of black bean aphids (n>an class="Species">Aphis fabae) influenced oviposition by a parasitoid wasp (Aphidius matricariae). They showed that mother parasitoid wasps were more attracted to the plant host environment that produced offspring with the lowest fitness, but oviposition frequency was highest on the plant species that resulted in the fittest offspring [15]. Although maternal effects were not tested specifically in the latter study, they have potential to influence the regulation of aphid populations [16,17]. Parasitoid conditioning to previously-experienced host types can influence subsequent preference or willingness to oviposit [18,19], indicating that acquired and innate preferences could influence parasitoid decisions with downstream consequences for the fitness of their offspring. The mechanism underlying these effects on parasitoid wasp fitness is unknown but could be due to direct plant effects on wasp behaviour or indirect effects on aphid quality for parasitism. One way that plants could influence aphid quality for parasitism is through the provision of nutrients. Manipulation of plant quality via nitrogen fertilisation had limited effect on fitness of the parasitoid Diaeretiella rapae attacking two aphid species (Myzus persicae and Brevicoryne brassicae) despite significant effects on aphid fitness and aphid nutritional quality[20]. However, the effect of changing host identity at plant and aphid trophic levels on fitness of parasitoid mothers and their offspring still remains to be explored.

In this study, we investigate whether the maternal environment affects the behaviour and pan class="Disease">fitness of pan class="Species">pea aphid (pan class="Species">Acyrthosiphon pisum) and the generalist parasitoid wasp Aphidius ervi. Although information is limited in the literature, both generalist and specialist insects are likely to experience fitness costs when transferred between different maternal hosts, with higher costs detected for specialist aphids [21-23]. We tested attractiveness of maternal and alternative environments to adult aphids and parasitoid wasps, and quantified effects on their fecundity. Maternal effects were tested by first evaluating the impact of changing pea aphid host environment between faba bean (Vicia faba) and pea (Pisum sativum) on pea aphid offspring fitness. Second, we assessed whether changing the maternal environment of parasitoid wasps at a single trophic level, using pea aphids reared on either bean or pea, affected fitness of parasitoid wasps and their offspring. Third, we tested parasitoid wasp fitness in response to manipulation of the parasitoid wasp host environment at two trophic levels using pea aphids reared on bean or potato aphid (Macrosiphum euphorbiae) reared on tomato (Solanum lycopersicum). We hypothesised that the maternal environment would have cascading effects on the fitness of aphid and parasitoid wasp offspring. We predicted that, first, adult pea aphids and adult parasitoid wasps (Generation 0, G0) would prefer the plant or plant+aphid host environment that formed their own host environment (Prediction 1). Second, we predicted that pea aphid and parasitoid wasp offspring (Generation 1, G1) would have higher fitness on the plant or plant+aphid host environment that formed the G0 host environment (Prediction 2). Finally, we predicted that wasp mothers (Generation 0, G0) would have higher fecundity on the plant or plant+aphid host environment that formed her own developmental host environment (Prediction 3).

Methods

Experimental design

To test the predictions, choice tests and performance assays were performed on pea aphids and parasitoid waspn>s. We tested the first prediction that adult G0 (Generation 0) pan class="Species">pea aphids and parasitoid wasps prefer the host environment they had been reared in. Adult G0 insects were given a choice comprising either the host environment they had been reared in or an alternative host environment (Fig 1). Three comparisons were undertaken: i) pan class="Species">pea aphid preferences for different plant hosts (bean vs. pea); and parasitoid wasp preferences for ii) pea aphids on different plant hosts (bean vs. pea; termed ‘plant comparison’) or iii) different plant+aphid host combinations (pea aphid on bean vs. potato aphid on tomato; termed ‘plant-aphid comparison’).

Fig 1

The experimental design of choice tests and performance assays.

The host environment that the G0 insects experienced and the host environments they were offered in choice tests and that G1 insects were reared on in performance assays. In choice tests, adult G0 insects were transferred to an olfactometer connected to two host environments (see Fig 2). In performance assays, adult G0 insects were transferred to the same or alternative host environment immediately before nymph deposition or oviposition of G1. By transferring reproductive G0 adults, G1 nymphs had no prior exposure to their maternal environment.

Fig 2

The olfactometer set up for choice tests.

(A) A diagram of the experimental set-up for the choice experiments and (B) Annotated photograph of an olfactometer displaying size measurements.

pan class="Chemical">Performance assays were conducted for pan class="Species">pea aphids and parasitoid wasps to test the predictions that the mother G0 host environment affects offspring G1 performance (pan class="Chemical">Prediction 2) and her fecundity (Prediction 3). This assay compared G1 performance in the G0 host environment with G1 performance in an alternative host environment (Fig 1). The experiment comprised a cross-over design of four plant or plant+aphid host combinations. Three comparisons were undertaken: i) offspring pea aphid fitness was examined in response to two different maternal plant hosts (bean vs. pea) and parasitoid wasp offspring fitness was examined in response to maternal experience of ii) pea aphids on two different plant hosts (bean vs. pea; plant comparison) or iii) two different plant+aphid host combinations (pea aphid on bean vs. potato aphid on tomato; plant-aphid comparison). G0 pea aphids and parasitoid wasps were raised to adulthood on their designated host environment and transferred to either the same host environment or an alternative host environment and allowed to deposit G1 pea aphid nymphs or oviposit G1 parasitoid wasps. Transfer of reproductive G0 adults ensured that G1 offspring had no prior exposure to their maternal environment. For aphids, G0 adults were removed within 24–48 h of G1 nymph deposition to minimise the G0 mother’s effect on G1 offspring. We also measured the dry weight of G2 nymphs produced from each G1 adult. G0 mother parasitoid wasps were removed from their oviposition environment when they had oviposited in 30 aphids (see details below).

pan class="Chemical">Pn>lant and aphid host quality was assessed by quantifying tissue dry weight and n>an class="Chemical">nitrogen concentrations although these data could not be collected for parasitoid wasp performance assays as this would have compromised collection of mummification data.

Plant and insect rearing

Bean (pan class="Species">Vicia faba cv. Sutton Dwarf) and pea (n>an class="Species">Pisum sativum cv. Douce Provence) seeds were planted individually into 3L pots filled with compost (Everris, UK, 2016; base fertiliser: 156 mg/L total N, 78 mg/L pan class="Chemical">P and 259 mg/L K, controlled release fertiliser: 480 mg/L total N, 117 mg/L P and 298 mg/L K and 36 mg/L Mg). Tomato (Solanum lycopersicum cv. Money Maker) seeds were germinated at 22°C then transferred individually into 1L pots filled with compost (as already described). Plants were grown in a greenhouse with 16:8 h (L:D) and 20°C:14°C and watered daily. Aphid clonal lines were confirmed to be free of facultative endosymbionts, such as Hamiltonella defensa [24], which can reduce their susceptibility to parasitism. Aphids, including G0, were reared on excised bean or pea cuttings (pea aphids, Acyrthosiphon pisum, line LL01) or excised tomato cuttings (potato aphid Macrosiphum euphorbiae, line AK15/01), which were replaced weekly, in ventilated plastic cups for at least four generations prior to use. This approach produced both winged and non-winged aphids in the same cup suitable for use in experiments. Winged aphid density was generally low, due to low numbers of aphids and use of high quality plant material, and these culture conditions produced a small number of winged aphids, which were sufficient in number for the olfactometer studies. Between 10–20 cups were used to culture aphids and aphids from different culture cups were mixed and randomly selected for experiments. Parasitoid wasps (Aphidius ervi) were obtained from Syngenta (Fargro, West Sussex, UK; batch number 31901) and all wasps used in experiments, including G0, were reared on pea aphids on bean or pea plants, or potato aphids on tomato plants, using aphids that were less than four days old, for at least one generation prior to experiments. When mummies had formed on plants, the leaves that the mummies had formed on were removed and placed into ventilated plastic boxes until used in experiments. Adult parasitoid wasps were fed a diluted honey solution (50% v/v) presented in cotton wool. Female parasitoid wasps, aged two to five days old and presumed mated, were used in experiments. Insect cultures were maintained at 16:8 h (L:D), 20°C:14°C and 70% humidity.

Insect choice experiments

Two-way choice tests were conducted with adult G0 winged aphids or parasitoid wasps (Fig 1). Experiments were performed under rearing conditions (see above) between 10:00 and 12:00 h using a two-armed olfactometer connected to the appropriate treatments. Three-week old plants were enclosed in polyethylene terephthalate (PET) bags and sealed at the plant base with inert plastic ties. Three days prior to parasitoid wasp choice tests, 30 aphid nymphs (2nd to 3rd instar) reared on the appropriate plant host were placed in mesh-covered clip cages (25 mm internal diameter) that were fixed onto leaves of experimental plants. Each arm of the olfactometer was attached to a small hole created in the corner of the PET bag enclosing experimental plants and sealed using polytetrafluoroethylene (PTFE) tape. Inert fabric, attached by PTFE tape, covered the ends of the olfactometer to prevent insect escape but allowing air flow (see Fig 2). After one hour, six winged adult pea aphids, which had been starved for a minimum of two hours, or six adult female parasitoid wasps, were placed in the centre of the olfactometer chamber and the position of each of the six insects was noted every two minutes for a 60-minute observation period. Insects were only used once. Olfactometers were cleaned for each choice test using a dilute solution of teepol, then rinsed with deionised water. Plants used in pea aphid choice tests were immediately harvested and dried at 60°C for two days to quantify total plant dry weight. A subsample of the dried leaves was randomly selected from each plant and ball milled to a fine powder. Aphids used in the parasitoid wasp choice tests were removed from plants, frozen at -20°C and freeze-dried to quantify dry weight. The tissue nitrogen concentration of milled plant material and freeze-dried aphids was determined by elemental analysis using a CE-440 Elemental Analyzer (Exeter Analytical Inc., North Chelmsford, Massachusetts, USA).

Insect performance assays

Pea aphids

Two G0 adult pan class="Species">pea aphids were transferred from plant cuttings to a 3-week old plant of the apn>propn>riate host environment (Fig 1) and secured by a clipn> cage. Expn>eriments were conducted in glasshouse conditions (as described above for plant rearing conditions). Following production of the first cohort of G1 nympn>hs (within 24–48 h of transferring the G0 adult), the G0 apn>hids and all but three G1 nympn>hs were removed. The developn>ment of three G1 nympn>hs was monitored to adulthood, after which apn>hids were removed from the plant so that a single apn>hid remained in each cage. n>an class="Chemical">Performance of the one remaining G1 aphid was monitored for a maximum of twelve days of adulthood in order to calculate intrinsic rate of population growth (Rm) [25] by:

We measured the number of days the G1 aphid took to deposit its first nymph (development time) and the number of nymphs (G2) the adult G1 aphid produced daily (fecundity). If G1 pea aphids died during the expn>eriment, the number of days the G1 apn>hid had lived was recorded as ‘survival’. Additionally, the dry weight of all G2 nympn>hs produced by each G1 adult was measured. n>an class="Chemical">Plants used in performance assays were harvested whole and dried at 60°C for two days to quantify shoot and root dry weight and leaf tissue nitrogen concentration as described above.

Parasitoid wasps

G0 wasps reared in each plant or plant+aphid host environment were transferred into a parasitism ‘arena’ compn>rising an excised leaf, abaxial surface upn>permost, immobilised in 1% pan class="Chemical">agarose (w/v) in a plastic petri dish (9 cm internal diameter). The leaf was infested with 30 aphid nymphs (2nd or 3rd instar); the plant or plant+aphid host environment presented in the arena reflected either the same G0 parasitoid wasp developmental host environment or an alternative host environment (Fig 2). Assays were performed under insect rearing conditions (described above). One mother G0 parasitoid wasp was introduced per arena and allowed to attack nymphs and oviposit until she had oviposited in 30 aphid nymphs; the time it took her to do this was recorded. Aphid nymphs that had been attacked were transferred onto plant cuttings of the same plant host environment used in the arena, contained in ventilated plastic cups. When an attacked nymph was removed, it was replaced with a naïve nymph to maintain a constant aphid density in the arena.

The number of mummies that developed in attacked aphids (after 14 d) was recorded as a measure of the G0 wasps’ fecundity. Whether or not a wasp emerged from each mummy was recorded as a measure of offspring G1 pan class="Disease">fitness, along with the sex and dry weight of each G1 waspn>.

Statistical analyses

Binomial general linear mixed models (GLMMs) were used to test for effects of host environment on pan class="Species">pea aphid and parasitoid waspn> behaviour (n>an class="Chemical">Prediction 1) in choice tests, using date of assay (tests were conducted over 3–7 d) and replicate as random factors. The proportion of insects (six in total, before excluding non-responders) spending time in each sector of the olfactometer chamber was calculated for each of the thirty time points for ten replicates per treatment, or eight for parasitoid wasps reared on potato aphids and tomato plants, and used in the analysis. pan class="Chemical">Paired t-tests were used to test for differences in tissue dry weight and nitrogen concentration between bean and pea plant pairs used in pea aphid choice tests and for aphids used in parasitoid wasp choice tests.

For pan class="Species">pea aphids, two-way analysis of variance (ANOVA) was used to test for maternal effects of host environment on G1 intrinsic rate of popn>ulation growth and G2 nympn>h weight (n>an class="Chemical">Prediction 2). Fixed effects in the models included the G0 and G1 host environment, and the interaction, and the position of the pea aphid assay in the greenhouse was included as a random effect. Twelve replicates were used. Two-way ANOVAs were also used to test the bean and pea host dry weight and tissue nitrogen concentration with the same fixed and random effects, as described above. The correlation was tested between G1 intrinsic rate of population growth and host plant tissue pan class="Chemical">nitrogen concentration. Pea aphid survival (measured as the number of days the G1 aphid lived) was analysed by fitting the survival data to a Cox proportional hazards regression model. The effect of the G0 and G1 host environment, and the interaction, was tested using analysis of deviance. During model simplification, analysis of deviance (using a Type 2 Wald chi-squared test) and Akaike’s Information Criterion (AIC) were used to check for model suitability.

Maternal effects on wasp survival (emerged vs. un-emerged wasps) and the sex (male vs. female) of the emerged G1 parasitoid wasps were tested using separate binomial GLMMs and logit links (due to using binary datasets) (Prediction 2). Maternal effects on G1 parasitoid waspn> dry weight were tested using a linear mixed model (LMM). All models included the G0 and G1 host environment, and the interaction, as fixed effects and the date on which the parasitoid waspn> assay was performed and the assay replicate number as a random effect. Ten repn>licates were used. In addition, the LMM used to analyse G1 parasitoid waspn> dry weight included G1 waspn> sex as a fixed factor. The ‘Anova’ function was used to obtain chi-squared, degrees of freedom and p values to establish significant differences between levels of each fixed effect. The effects of host environment on G0 fecundity (i.e. time taken for G0 parasitoid waspn>s to attack 30 apn>hid nympn>hs and number of waspn> mummies) were analysed using two-way ANOVAs, with G0 and G1 host environment, and the interaction, as fixed effects and the date on which the waspn> assay was performed as a random effect (n>an class="Chemical">Prediction 3). Model simplification was carried out as described above. Note that data collected for G1 parasitoid wasps which developed on bean-reared pea aphids, whose G0 parasitoid mothers also developed on bean-reared pea aphids, were used for statistical analysis in both the plant and plant-aphid comparisons.

The paired t-tests and ANOVAs were conducted using GenStat [26]. Data satisfied the requirements of parametric testing for equal homogeneity of variance and normal distribution. Survival analysis was conducted using RStudio version 3.2.4 ‘Very Secure Dishes’ [27] using the ‘coxph’ and ‘Surv’ functions in the ‘Survival’ package [28]. GLMMs and LMMs were carried out using ‘glmer’ and ‘lmer’ functions, respectively, using the ‘pan class="Chemical">lme4’ [29], ‘car’ [30] and ‘lmerTest’ [31] packages in RStudio [27]. Maternal effects were indicated when the interaction between the G0 and G1 host environment was significant at the 5% level.

Results

Do adult (G0) insects prefer their rearing host environment (Prediction 1)?

G0 pan class="Species">pea aphidsn> showed no preference for bean or pea plants irrespn>ective of the plant host environment they experienced (p>0.05; Fig 3A, S1 Table). Bean plants used in choice tests were significantly heavier than pea plants, but tissue n>an class="Chemical">nitrogen concentration was similar for the two species (for statistical outputs see S2 Table).

Fig 3

Insect host environment preferences.

The percentage (%) of time spent in each half of the olfactometer chamber, exhibited by (A) pea aphids for bean (filled bars) or pea (open bars) plants, (B) parasitoid wasps for pea aphids on bean (dark grey bars) or on pea plants (light grey bars) in the plant comparison, and (C) parasitoid wasps for pea aphids on bean (dark grey bars) or potato aphids on tomato plants (light grey hatched bars) in the plant-aphid comparison. Values are means (± SEM) of n = 10 for all choice tests, except for G0 wasps reared on potato aphids on tomato plants when n = 8. *signifies bars are significantly different. For statistical summary see S1 Table.

The percentage (%) of time spent in each half of the olfactometer chamber, exhibited by (A) pea aphids for bean (filled bars) or pea (open bars) plants, (B) parasitoid waspn>s for n>an class="Species">pea aphids on bean (dark grey bars) or on pea plants (light grey bars) in the plant comparison, and (C) parasitoid wasps for pea aphids on bean (dark grey bars) or potato aphids on tomato plants (light grey hatched bars) in the plant-aphid comparison. Values are means (± SEM) of n = 10 for all choice tests, except for G0 wasps reared on pan class="Species">potato aphids on tomato plants when n = 8. *signifies bars are significantly different. For statistical summary see S1 Table.

Parasitoid Wasps: Plant comparison

G0 wasps that develon class="Chemical">ped in pea aphids reared on bean or pea also showed no pn>reference for either plant typn>e (p>0.05; Fig 3B, S1 Table). Bean-reared pea aphids used in the choice tests were significantly heavier but had the same nitrogen concentration as pea-reared pea aphids (aphid dry weight: T18 = 2.96, p = 0.008; aphid nitrogen concentration: T18 = 0.56, p = 0.585; S1 Fig).

Parasitoid Wasps: Plant-aphid comparison

n class="Chemical">Pn>arasitoid wasps that had developed in bean-reared n>an class="Species">pea aphids showed a preference for tomato-reared potato aphids over bean-reared pea aphids (z = 2.445, p = 0.015; Fig 3C; S1 Table). No preference was observed for parasitoid wasps that had developed on tomato-reared potato aphids (p>0.05; Fig 3C, S1 Table). Pea aphids used in the choice tests were also heavier than potato aphids (T17 = 11.90, p<0.001; S1 Fig), but had significantly lower nitrogen concentration (6.21% for pea aphids vs. 7.21% for potato aphids, on average: T17 = 4.35, p<0.001; S1 Fig).

Is insect offspring performance driven by the maternal host environment (Prediction 2)?

There was little evidence for maternal effects on G1 pan class="Species">pea aphid pan class="Disease">fitness parameters. Values for G1 pan class="Species">pea aphid fitness were highest when reared on bean plants: intrinsic rate of G1 population growth was higher (G1: F1 = 90.62, p<0.001), G2 nymphs were heavier (G1: F1 = 48.48, p<0.001) and G1 aphids survived longer (G1: X21 = 21.95, p<0.001; Fig 4A–4C; S3 Table). However, when the G0 environment was bean, this also resulted in a higher intrinsic rate of G1 population growth compared to when the G0 environment was pea (G0: F1 = 15.47, p<0.001; Fig 4A; S3 Table), indicating a maternal effect. At the end of the performance assays, pea plants had a larger shoot dry weight (but not total weight) than bean plants, and bean plants had a higher leaf nitrogen concentration than pea plants (for statistical outputs S4 Table). There was a significant interaction between the pea aphid G0 and G1 host environment on plant tissue nitrogen concentration (G0*G1: F1 = 15.96, p<0.001). This was due to lower N concentrations in pea plants when infested with G1 pea aphids whose G0 mothers had been reared on pea plants (Fig 5). Indeed, the pattern observed for intrinsic rate of G1 population growth mirrored that of the host plant leaf tissue nitrogen concentration (Pearson’s correlation = 0.737, p<0.001).

Fig 4

Performance of pea aphid offspring (G1).

(A) G1 intrinsic rate of population increase (Rm), (B) G2 nymph dry weight and (C) G1 survival. For (A) and (B), values are means (± SEM). Bars that share the same letter are not significantly different from each other. For (A), (B) and (C) the number of G1 aphid assays was n = 11 for G0 reared on bean and G1 reared on bean, n = 10 for G0 reared on bean and G1 reared on pea and G0 reared on pea and G1 reared on bean, and n = 9 for G0 reared on pea and G1 reared on pea.

Fig 5

Leaf tissue nitrogen concentration of plants used in pea aphid performance assays.

Bean (shaded bars) and pea (hatched bars). Values are means (± SEM) of n = 12. Bars that share the same letter are not significantly different from each other.

Parasitoid wasps: Plant comparison

G1 parasitoid wasp survival did not vary with the mother G0 nor offspring G1 host environment (p>0.05; Fig 6A; for statistical outputs see S5 Table). The majority of G1 emerged parasitoid wasps were male (74 ±0.2%), and G1 host environment had a significant effect on the sex of successfully emerged G1 parasitoid wasps (G1: X21 = 4.817, p = 0.028; Fig 6B, S5 Table), as more males were produced from bean-reared aphids. G1 parasitoid wasp dry weight was not dependent on sex (p = 0.062), but the interaction between the G0 and G1 host environment was significant (G0*G1: X21 = 8.603, p = 0.003; Fig 6C, S5 Table), with the heaviest wasps experiencing a bean-reared pan class="Species">pea aphid G0 and G1 host environment and lightest waspn>s expn>eriencing a pea-reared n>an class="Species">pea aphid G1 host environment.

Fig 6

Performance of offspring (G1) parasitoid wasps in relation to G0 and G1 host environment.

Percentage of G1 wasps that emerged from their mummy and survived to adulthood in (A) the Plant Comparison (p>0.05) and (D) the Plant-Aphid Comparison; the sex of G1 wasps in (B) the Plant Comparison and (E) the Plant-Aphid Comparison; and the weight of G1 wasps in (C) the Plant Comparison and (F) the Plant-Aphid Comparison. Values are means (± SEM) of n = 10 for G0 wasps. nf and nm represent the number of females and males, respectively. For statistical summaries see S5 Table.

pan class="Chemical">Percentage of G1 wasps that emerged from their mummy and survived to adulthood in (A) the n>an class="Chemical">Plant Comparison (p>0.05) and (D) the Plant-Aphid Comparison; the sex of G1 wasps in (B) the pan class="Chemical">Plant Comparison and (E) the Plant-Aphid Comparison; and the weight of G1 wasps in (C) the Plant Comparison and (F) the Plant-Aphid Comparison. Values are means (± SEM) of n = 10 for G0 wasps. nf and nm represent the number of females and males, respectively. For statistical summaries see S5 Table.

Parasitoid wasps: Plant-aphid comparison

G1 parasitoid wasp survival to adulthood was only explained by the G0 host environment (G0: X21 = 5.179, p = 0.023, Fig 6D, S5 Table). In accordance with the plant comparison, the majority of parasitoid wasps were male (74±0.3%), but the interaction between the G0 and G1 plant+aphid developmental host environment was significant (G0*G1: X21 = 5.193, p = 0.023; Fig 6E, S5 Table): A smaller male sex bias was observed for wasps with a tomato-reared n>an class="Species">potato aphid G0 host environment and a bean-reared pea aphid G1 host environment compared with other host environment combinations. G1 parasitoid wasp weight was not explained by G1 parasitoid wasp sex (p = 0.480), but sex did interact with the G1 host environment: Females were heavier with a bean-reared pea aphid G1 host environment whilst males were heavier with a tomato-reared potato aphid G1 host environment. Also, the G0 and G1 host environment interaction was significant (G0*G1: X21 = 13.856, p<0.001; Fig 6F, S5 Table), with the heaviest wasps experiencing a bean-reared pea aphid G0 and G1 host environment and lightest wasps experiencing a bean-reared pea aphid G0 and tomato-reared potato aphid G1 host environment.

Is mother wasp fecundity driven by her host environment (Prediction 3)?

In the plant comparison, mother G0 parasitoid wasps attacked bean-reared pea aphid nympn>hs more rapn>idly than pea-reared pan class="Species">pea aphid nymphs (G1: F1 = 18.01, p<0.001; Fig 7A), but the number of mummies produced per G0 parasitoid (and the number of subsequent G1 parasitoids that emerged from mummies) were highest when the maternal and offspring rearing environments were bean-reared pea aphid (G0*G1: F1,39 = 6.92, p = 0.014; Fig 7B). In the plant-aphid comparison, mother G0 parasitoid wasps attacked bean-reared pea aphids significantly faster than they attacked tomato-reared potato aphids (G1: F1 = 7.58, p = 0.010; Fig 7C), but the greatest number of mummies (and the number of subsequent G1 parasitoids that emerged from mummies) was observed when the G0 and G1 plant+aphid developmental host environments were the same (G0*G1: F1,39 = 22.07, p<0.001; Fig 7D).

Fig 7

Performance of mother (G0) parasitoid wasps in relation to the G0 and G1 host environment.

Quantified as the time (mins) taken for G0 wasps to oviposit into 30 aphids in (A) the Plant Comparison and (C) the Plant-Aphid Comparison and the number of G1 mummies formed by G0 wasps successfully ovipositing into aphids in (B) the Plant Comparison and (D) the Plant-Aphid Comparison. Values are means (± SEM) of n = 10. Bars that share the same letter are not significantly different from each other.

Quantified as the time (mins) taken for G0 wasps to oviposit into 30 aphids in (A) the Plant Compn>arison and (C) the n>an class="Chemical">Plant-Aphid Comparison and the number of G1 mummies formed by G0 wasps successfully ovipositing into aphids in (B) the Plant Comparison and (D) the Plant-Aphid Comparison. Values are means (± SEM) of n = 10. Bars that share the same letter are not significantly different from each other.

Discussion

Here we report experiments that test whether the maternal environment affects the behaviour and fitness of n>an class="Species">pea aphids and parasitoid wasps, which are both ecologically and commercially important insects. By manipulating host environments at different trophic levels, we revealed contrasting effects of maternal host environment on pea aphids and their parasitoid natural enemy, which aligned partially with our predictions.

We predicted that, given a choice, n class="Chemical">pan class="Species">pea aphids and pn>arasitoid wasps would preferentially select the host environment that they experienced during development. However, in all but one case, aphids and waspn>s showed no preference for either maternal or alternative host environments, and this finding did not support our first prediction. Only parasitoid waspn>s that had developed on bean-reared n>an class="Species">pea aphids showed any preferences, which was for tomato-reared potato aphids. However, the fecundity of parasitoid mothers and fitness of their offspring (assessed by their weight) were unaffected in this particular host environment.

We had expected pan class="Species">pea aphids to prefer the host plant they were reared on, as previous studies found that n>an class="Species">pea aphid genotypes specialised to different legume plants preferred the legume they were adapted to rather than other leguminous species [12]. The olfactometer methodology we used relied on insects responding to the volatile organic compounds released by plants (see [32]), and prevented pea aphids from probing potentially favourable parts of the plant to assess host quality, including phloem sap quality [33], which ultimately affects aphid feeding. Aphids used in choice tests were reared on plant cuttings, rather than whole plants, which may have exposed the aphids to a variety of plant volatile cues and hence influenced aphid behaviour towards volatiles in the olfactometer tests. The leaf tissue nitrogen and pan class="Chemical">carbon concentrations and C:N ratio (data not shown) were similar between bean and pea plants used in choice tests, which may have contributed to the lack of pea aphid preference, although recent work has also found that leaf phosphorus concentrations do not drive aphid attractiveness to bean plants [34].

Like aphids, parasitoid wasps use a variety of methods to assess host quality in order to make the ultimate decision to oviposit using a variety of visual and physical cues. These include aphid size (instar), shape, colour and movement [35,36], probing aphid hosts and aphid defensive behaviours, like body raising, kicking and body rotation [37]. Waspn>s were unable to use these cues to test for host quality in our choice tests and were only able to use plant volatiles as cues. We found that waspn>s reared on bean-reared n>an class="Species">pea aphids preferred tomato-reared potato aphids, suggesting preference for volatiles released by tomato plants. Aphidius ervi has been shown to respond to similar volatiles, including methyl salicylate, released from both pea aphid-infested bean plants [38] and potato aphid-infested tomato plants[39], and we therefore need to undertake more detailed analyses of volatile compounds released from host plant leaves to provide a mechanistic understanding of our findings.

The pan class="Chemical">nitrogen concentration of pan class="Species">pea aphids reared on either bean or pea plants were similar, but pan class="Species">potato aphids had higher nitrogen concentration than pea aphids in choice tests. However, pea aphids were bigger than potato aphids, and hence did not offer a larger nitrogen resource overall. Plant nitrogen status has been shown to influence volatile emissions: in soybean plants (Glycine max), nitrogen-starved plants had a similar volatile composition to that of nitrogen-fertilised plants, but three compounds were released in different quantities; however, these differences were undetectable when plants were attacked by fall armyworms (Spodoptera frugiperda) and parasitoid wasps (Cotesia marginiventris) showed no preference between infested plants that had been nitrogen starved or fertilised [40]. The olfactory cues influencing parasitoid wasp choice in our study remain to be elucidated, particularly how they might relate to plant and aphid nutrient status.

The performance of aphid offspring, measured by nymph weight and survival, experiencing the same host environment as their mother’s host environment did not differ from those experiencing an alternative plant host environment. However, the intrinsic rate of population increase was explained by the maternal and offspring environment (although not by an interaction), lending a little support for our second prediction. Pea aphid offspn>ring had a higher intrinsic rate of popn>ulation increase, heavier nympn>hs and survived longer on a bean host environment. Maternal environment might be anticipn>ated to be a key determinant of offspn>ring n>an class="Disease">fitness in clonally-reproducing aphids. Previous research has shown, however, that pea aphids perform as well or better on bean plants regardless of their maternal hosts [13] and McLean et al. [7] also found no maternal effects on offspring fecundity of host-adapted pea aphid genotypes across multiple generations when host swapping between Lathyrus pratensis and bean plants.

Aphid infestation [41] or infestation by different aphid species [42,43] often leads to differences in plant tissue nutrient concentrations, including nitrogen and soluble amino acids. Here, we repn>ort that the prior plant host environment of an individual aphid can affect the nutrient concentration of their host plant and that this is positively correlated with apn>hid intrinsic rate of popn>ulation growth. Lower leaf tissue n>an class="Chemical">nitrogen concentrations and aphid intrinsic rate of population growth were observed in pea plants harbouring pea aphids whose mothers had fed on pea rather than bean plants. Aphids tend to be nitrogen limited [44] and these findings indicate that maternal plant quality can affect aphid physiology, feeding and fitness in a way that influences the aphid’s ability to utilise future plant host resources. Further work might reveal the mechanism of this effect, for example via induced resource sequestration [45], where resources are re-allocated to different plant structures. Indeed, measuring nitrogen content of phloem sap may shed further light on this phenomenon.

Our second prediction, that parasitoid wasp offspring fitness proxies would be highest when offspn>ring expn>erienced their mother’s environment, was not supported by offspn>ring waspn> emergence and waspn> sex ratio data when only maternal plant environment was manipn>ulated (i.e. plant compn>arison). However, in the plant-apn>hid compn>arison, parasitoid offspn>ring had a lower male bias with a n>an class="Species">tomato-reared potato aphid maternal host environment and bean-reared pea aphid offspring host environment compared to other host environment combinations. A. ervi has a haploid-diploid reproduction strategy [46], allowing mothers to choose whether to oviposit sons or daughters into aphid hosts. This choice is often decided by assessment of host quality, which could include the lipid content of aphid hosts [47]. In dense populations mothers choose to produce more males to outcompete mating rivals for females [46, 48]; mass rearing in the laboratory can simulate these conditions [48], which could explain why the majority of wasps observed in our study were male. Alternatively, mother wasps, although presumed mated based on observations in the wasp cultures, might not have been mated prior to experiments. However, mothers with a developmental host environment of tomato-reared potato aphids oviposited more daughters in bigger pea aphids that represent a more plentiful resource for offspring. Offspring wasp body mass is primarily determined by aphid body mass [49], but in our study the offspring host environment also interacted with the maternal host environment in both the plant and plant-aphid comparisons. Specifically, the biggest wasps were those that developed in pea aphids reared on bean in both the maternal and offspring generations. This finding demonstrates that bean reared pea aphids are an overall better host, but these beneficial effects are maximised when mother wasps also developed in bean-reared pea aphids: this supports our second prediction and evidences potential maternal effects. Body mass is an important fitness indicator for parasitoid wasps because it is often positively correlated with longevity, host and mate searching rate, fecundity and ability to parasitize [50]. The benefits of producing larger and fitter offspring could explain why G0 wasps oviposited more rapidly in the largest aphid hosts (bean-reared pea aphids) in both the plant and plant-aphid comparisons, but this preference was not observed in choice tests. However, smaller wasps were produced in a different plant+aphid host environment compared to only a different plant host environment: Offering wasps an aphid host that represented a resource of differing quality seemed to have a stronger impact on wasp fitness than changing the plant aspect of the host environment, although these factors are linked, making interpretation complex.

These findings can help inform biological control methods. For example, by choice of a favourable host environment combination, mass rearing of parasitoid wasps could be optimised to maximise parasitoid wasp mothers’ fecundity and improve the potential of their offspring to regulate aphid pests. However, the fitness of the offspn>ring is likely to be maximised using larger apn>hid hosts, which could result in mother-offspn>ring conflict if mothers preferentially ovipn>osit in smaller apn>hids, resulting in short-term negative maternal effects on offspn>ring n>an class="Disease">fitness that compromise long-term aphid biocontrol. Further work is needed to understand the impact of environmental heterogeneity and trophic complexity on maternal effects. For example, positive maternal effects on offspring fitness are likely to be compromised by fluctuations in the insect environment that lead to unpredictability in host quality and suitability. An interesting avenue for further work is to determine whether the positive maternal effects observed in this study could affect other trophic levels (e.g. hyperparasitoids) and lead to cascading maternal effects through the food web. The hypotheses tested in this study could be applied to additional trophic groups, although challenging to test empirically due to the experiment size doubling each time a treatment level is added.

Our study shows that the plant and aphid host environments of wasps can interact to affect wasp fitness, but the exact outcomes are compn>lex and inconsistent. Specifically, waspn>s reared on bean plants and pan class="Species">pea aphids appeared to have the highest pan class="Disease">fitness, possibly as pea aphids themselves had the highest fitness on bean plants. In order to better gather evidence for maternal effects in wasps, more plant and aphid host environment combinations need to be tested. This information could help facilitate the most appropriate rearing environments used in the biocontrol industry.

Quality of aphids used in G0 wasp choice tests.

Aphid dry weight (mg) for (A) the pan class="Chemical">Plant Comparison and (C) the pan class="Chemical">Plant-Aphid Comparison. Aphid pan class="Chemical">nitrogen concentration for (B) the Plant Comparison and (D) the Plant-Aphid Comparison. Values are means (± SEM) of n = 20 for pea aphids reared on bean plants and n = 19 pea aphids reared on pea plants in the Plant Comparison and n = 18 for pea aphids reared on bean plants and n = 18 for potato aphids reared on tomato plants for the Plant-Aphid Comparison. ** p<0.01, *** p<0.001.

(TIF)

Click here for additional data file.

Statistical summaries of the pea aphid and wasp choice tests shown in Fig 3.

Negative estimates and z values indicate preference for the ‘alternative’ compared to the ‘same’ host environment (Fig 1) presented in the choice tests. Significant preferences are highlighted in bold.

(DOCX)

Click here for additional data file.

Details on plants used in G0 pea aphid choice tests.

Weight (g) and leaf pan class="Chemical">nitrogenn> concentration (% dry mass) of 3-week old bean. Values are means (± SEM) of n = 12 plants. Significant differences are highlighted in bold.

(DOCX)

Click here for additional data file.

Statistical summaries of two-way ANOVAs for pea aphid performance.

G1 intrinsic rate of pon class="Chemical">pulation increase (Rm), G2 nymph dry weight and G1 survival for pan class="Species">pea aphids used in performance assays (Fig 4). Significant differences are highlighted in bold.

(DOCX)

Click here for additional data file.

Details on plants used in G1 pea aphid performance assays.

Weight (g) and leaf pan class="Chemical">nitrogenn> concentration (% dry mass) of six-week old bean and pea plants, after being infested with G1 n>an class="Species">pea aphids for three weeks. Significant differences are highlighted in bold.

(DOCX)

Click here for additional data file.

Summary statistics of (general) linear mixed models ((G)LMMs) for wasp performance.

GLMMS for G1 wasp survival and sex, and LMMs for G1 wasp weight were used to test for maternal effects for the pan class="Chemical">Plant and n>an class="Chemical">Plant-Aphid Comparison (Fig 6). Statistical outputs are provided from the most simplified models (see text for details). Significant results are highlighted in bold.

(DOCX)

Click here for additional data file.