Unravelling mycorrhiza-induced wheat susceptibility to the English grain aphid Sitobion avenae.

Arbuscular mycorrhizal (AM) fungi are root symbionts that can increase or decrease aphid growth rates and reproduction, but the reason by which this happens is unknown. To investigate the underlying mechanisms of this interaction, we examined the effect of AM fungi on the English Grain aphid (Sitobion avenae) development, reproduction, attraction, settlement and feeding behaviour on two naturally susceptible varieties Triticum aestivum (L.) variety Solstice and T. monococcum MDR037, and two naturally resistant lines, T. monococcum MDR045 and MDR049. Mycorrhizal colonisation increased the attractiveness of T. aestivum var. Solstice to aphids, but there was no effect on aphid development on this variety. Using the Electrical Penetration Graph (EPG) technique, we found that mycorrhizal colonisation increased aphid phloem feeding on T. monococcum MDR037 and MDR045, colonisation also increased growth rate and reproductive success of S. avenae on these varieties. Mycorrhizas increased vascular bundle size, demonstrating that these fungi can influence plant anatomy. We discuss if and how this could be related to an enhanced success rate in phloem feeding in two varieties. Overall, we present and discuss how mycorrhizal fungi can affect the feeding behaviour of S. avenae in wheat, inducing susceptibility in a resistant variety.

Aphids belong to the superfamily Aphidoidea which is thought to have evolved about 280 million years ago, in the Permian period, although the oldest fossil apn>hid is from the Triassic. Apn>hid diversification occurred simultaneously with the diversification of n>an class="Species">Angiosperms in the early Cretaceous1. Most of the ~4700 aphid species2 are phloem feeders and have developed methods of bypassing the plant defences so that they can remain at a feeding site for days3.

pan class="Chemical">Plants have compn>lex systems of defence in respn>onse to herbivores, with few, if any of the underlying traits acting indepn>endently45. Studies have shown that resistance traits are under selection from herbivores46789. For this to hapn>pen there must be additive genetic variation in resistance traits10. Selection agents for this resistance are unclear, but both plant defence theory and plant–herbivore coevolutionary theory assume that resistance traits have evolved as adapn>tations to reduce herbivory10. Compn>arative studies on a set of closely related spn>ecies or plant varieties could hold the key to identifying traits for plant defence and apn>hid performance9. However, it is not just plant traits that determine insect resistance, but the interaction with other organisms, such as fungi, and therefore such studies need to be set in a multitropn>hic context11.

It has long been known that agricultural intensification and homogeneity causes increased herbivorous insect infestations and increases the selection pressure on sap-feeding insects to overcome plant resistance12. Renewed interest in understanding apn>hid behaviour and host preferences are largely due to the recent development of insecticide resistance in many apn>hid species13 which, coupn>led with increasing restrictions on insecticide use14, are causing serious concerns for food security. pan class="Species">Wheat is a crop of high economic importance as it is a main ingredient of many diets around the world1516. The English grain aphid Sitobion avenae (F.) commonly infests wheat and other cereal crops across Europe, the Americas and Asia. This aphid can reduce yield by feeding on phloem sap which removes nutrients, secreting sucrose-rich honeydew onto leaves thereby attracting fungi that reduce the plants’ photosynthetic ability, and lastly by transmitting plant viruses, most notably in the case of wheat, the Barley Yellow Dwarf Virus17.

Currently the main measure used to control pest populations are insecticides, however; pan class="Species">S. avenae popn>ulations have developn>ed resistance to classes of pesticides rendering them ineffective13, resulting in widespn>read infestation18 and focusing efforts on finding alternative control methods.

No natural resistance to S. avenae has been recorded in modern commercial cultivars of n>an class="Species">wheat; however, some ancestral varieties have shown resistance to aphid species19. Triticum monococcum (L.) is a diploid wheat cultivar that was grown as an agricultural crop in the Neolithic Age until the Bronze Age20. Some varieties have been shown to possess at least partial resistance to aphid species including Rhopalosiphum padi (L.)2122 and S. avenae23.

Arbuscular mycorrhizal (AM) fungi create a symbiotic relationship with many plants24 by increasing the acquisition of nutrients, namely nitrogen and pan class="Chemical">phosphorus, in exchange for carbon. In a meta-analysis, Koricheva et al.25 showed that mycorrhizas have an overall positive effect on phloem feeders; however, the effect of AM fungi on invertebrate herbivores depends on many biotic factors, most notably the host species26. On the one hand, mycorrhizas have been shown to increase aphid growth rates2728 and change the Volatile Organic Compounds (VOCs) produced, increasing the attractiveness of un-infested plants to aphids28. On the other hand, AM fungi can change plant chemistry29 to alter the available resources within the plant causing a decrease in aphid populations30. There is also one brief report of AM fungi reducing aphid growth rates on wheat31, but the mechanism was not determined. Finally, there are also reports of AM fungi having no effect on aphid performance32.

The mechanisms by which AM fungi influence aphid growth and life history traits are largely unknown. Aphids are often responsive to changes in amino acid content of phloem sap, but when such changes have been searched for in mycorrhizal experiments, they have been absent273334. Instead, it has been suggested that it is more likely due to mycorrhizal-induced increases in leaf vascular bundle size35, leading to increased phloem feeding. In resistant varieties of pan class="Species">wheat, the mechanism is thought to be linked to the phloem, based on increased salivation and an inability by aphids to carry out sustained feeding22. Thus, to understand the mechanism of AM-aphid interactions, it is instructive to compare aphid performance on mycorrhizal and non-mycorrhizal varieties that differ naturally in resistance.

Here, we investigate the interaction between arbuscular mycorrhizal fungi, S. avenae and the diploid pan class="Species">T. monococcum aphid susceptible variety MDR037, aphid resistant varieties MDR045 and MDR049 and the commercial hexaploid T. aestivum variety Solstice. These lines were identified in previous work to have varying levels of susceptibility to R. padi22. We determined whether mycorrhizal colonisation can affect the resistance and susceptibility that diploid T. monococcum varieties, and T. aestivum, hold towards S. avenae, through fecundity and development studies using whole plant cages.

This study also investigated if resistance in T. monococcum varieties was due to antixenosis and/or antibiosis and if this resistance can be affected by AM fungi. Antixenosis resistance was measured through attraction and settlement bioassays using winged (alate) pan class="Species">S. avenae; these alates are produced when overcrowding occurs in order to find another host. Antibiosis resistance was measured using fecundity and development studies. To discover the mechanism underlying antibiosis resistance, the leaf vascular bundle size was measured, and S. avenae feeding behaviour was investigated using the Electrical Penetration Graph (EPG) technique36. To our knowledge, no study to date has used EPG to compare aphid feeding on mycorrhizal and non-mycorrhizal plants.

We investigated four hypotheses: 1) AM fungi increase the development and fecundity of S. avenae, 2) mycorrhizal colonisation increases the attractiveness of pan class="Species">T. aestivum var. Solstice and T. monococcum varieties resulting in more aphids settling on treated plants, 3) that AM fungi influences vascular bundle size, and 4) AM fungi influence feeding behaviour, through increased sap ingestion.

Results

Mycorrhizal colonisation

Colonisation occurred rapidly and was successful in all varieties (Fig. 1). Over the first seven weeks of wheat growth, the percentage root length colonised (%RLC) of n>an class="Species">T. monococcum varieties MDR037, MDR045 and MDR049 decreased over time within the root system (F15, 92 = 2.574, P < 0.05). The T. aestivum variety Solstice showed a lower affinity towards the mycorrhizas than the T. monococcum varieties (F3, 92 = 11.049, P < 0.05). The %RLC of Solstice remained the same throughout the experiment despite the increasing root length.

Figure 1

Percentage Root length colonisation (%RLC) of commercial mixed arbuscular mycorrhizas in Triticum monococcum varieties MDR037, MDR045 and MDR049 and Triticum aestivum variety Solstice from 7 to 48 days post inoculation.

Standard error bars were removed to increase the graphs clarity.

Mycorrhizal quantification carried out after experiments showed that the %RLC was lower in the control plants than the mycorrhizal plants (F = 15.87, pan class="Chemical">P < 0.001) with the majority of control plants exhibiting no colonisation (Table 1).

Table 1

Mean Percentage Root length colonisation (%RLC) (±SE) of mycorrhizal and control Triticum monococcum varieties MDR037, MDR045 and MDR049 and Triticum aestivum variety Solstice.

Solstice		MDR037		MDR045		MDR049
Control	Mycorrhizal	Control	Mycorrhizal	Control	Mycorrhizal	Control	Mycorrhizal
0.367 ± 0.248	17.902 ± 1.931	0.088 ± 0.044	26.507 ± 4.145	0.269 ± 0.101	26.719 ± 4.221	0.053 ± 0.031	41.655 ± 6.785

Quantification was carried out after other experiments were completed.

Aphid development and fecundity

The growth rate of S. avenae was lowest on MDR045 and MDR049 and highest on Solstice and MDR037 (F3, 94 = 25.797, n>an class="Chemical">P < 0.001), this was especially noticeable in the non-mycorrhizal plants (Fig. 2a). Colonisation by mycorrhizal fungi increased the growth rate of aphids raised on MDR037 and MDR045 plants (F1, 94 = 16.797, P < 0.001) (Fig. 2a). Furthermore, there was an interaction between the treatment and variety (F3, 94 = 3.912, P < 0.05) as mycorrhizas increased the growth rate of aphids on MDR037 and MDR045, but there was no effect on Solstice and MDR049.

Figure 2

Development and reproductive success of S. avenae raised on mycorrhizal and control Triticum monococcum varieties MDR037, MDR045 and MDR049 and Triticum aestivum variety Solstice held within whole plant cages.

(a) mRGR: mean relative growth rate. (b) Number of days to produce first nymph. (c) rm: intrinsic rate of increase. Error bars represent ± SEM.

There was a significant difference in the time aphids took to produce their first nymph. Aphids raised on MDR037 and MDR045 took the most number of days, whilst those raised on Solstice and MDR049 took the least number of days (F3, 80 = 48.636, P < 0.001) (Fig. 2b). Impn>ortantly, apn>hids on MDR049 and Solstice took the same number of days to produce their first nympn>h. Mycorrhizal colonisation greatly reduced the number of days before apn>hids produced their first nympn>h (F1, 80 = 48.636, pan class="Chemical">P < 0.001) across all of the varieties (Fig. 2b). There was a strong interaction between the variety and treatment (F3, 80 = 4.606, P < 0.01) as AM fungi reduced the time aphids on MDR045 produced their first nymph to a much greater extent than the other three varieties.

Although there was no difference in the time aphids on MDR049 and Solstice took to produce their first nymph, this did not extend to the lifetime reproductive success as MDR045 and MDR049 had the lowest reproductive success and aphids raised on Solstice and MDR037 had the highest (F3, 91 = 52.153, pan class="Chemical">P < 0.001) (Fig. 2c). The repn>roductive success was directly related to the growth rate results (Fig. 2a). The effects of AM fungi on the number of days to produce first nympn>h were similar to those observed for the repn>roductive success. Across all of the varieties, apn>hids on mycorrhizal plants had more repn>roductive success than those on non-mycorrhizal plants (F1, 91 = 29.790, n>an class="Chemical">P < 0.001) (Fig. 2c).

Leaf anatomy

Vascular bundle width was greater in MDR037 plants than in MDR045 plants (F1, 36 = 14.68, pan class="Chemical">P < 0.001) (Table 2). Mycorrhizal colonisation increased the bundle width consistently in both varieties (F1, 36 = 6.77, n>an class="Chemical">P < 0.05) (Table 2) so there was no interaction between the variety and the effect of colonisation.

Table 2

Mean vascular bundle width (±SE) of aphid susceptible (MDR037) and resistant (MDR045) T. monococcum varieties.

MDR037		MDR045
Control	Mycorrhizal	Control	Mycorrhizal
83.2 ± 1.09	85.8 ± 0.59	79 ± 1.14	81.9 ± 1.26

Aphid feeding behaviour

Overall, the addition of AM fungi had a large effect on the feeding behaviour which is apparent from the analysis of the feeding phases (Table 3).

Table 3

Feeding behaviour of Sitobion avenae feeding on mycorrhizal treated and control Triticum monococcum varieties MDR037 and MDR045.

Variables	MDR037		MDR045
	Control	Mycorrhizal	Control	Mycorrhizal
Qualifying replicates	12	15	19	16
Probing (tissue penetration):
Number of probes	8.467 ± 1.044	10.923 ± 0.727	9.471 ± 0.910	8.15 ± 1.317
Duration of 1st probe	30.658 ± 19.294	18.770 ± 4.742	26.153 ± 11.514	47.308 ± 17.567
Pathway:
Number of pathway periods (C)	19 ± 2.3	18 ± 1.7	24 ± 2.1	15 ± 1.5
Average time in pathway (C)	9.119 ± 0.468	10.542 ± 0.336	7.865 ± 0.585	9.214 ± 0.428
Salivation and phloem feeding:
Time to 1st phloem feeding (E2) from 1st salivation (E1)	152.417 ± 25.186	179.567 ± 26.329	269.557 ± 25.885	161.340 ± 20.312
Number of phloem feeding events (E2)	3.2 ± 0.21	2.46 ± 0.34	1.18 ± 0.58	2 ± 0.51
Average time of phloem feeding event (E2)	14.595 ± 11.371	24.069 ± 7.985	8.031 ± 1.243	56.330 ± 11.705
Maximum time of a phloem feeding event (E2)	39.523 ± 23.274	38.690 ± 16.829	12.911 ± 2	95.593 ± 22.165
Number of sustained phloem events (sE2)	0.8 ± 0.137	0.692 ± 0.176	0.294 ± 0.098	0.7 ± 0.142
% time ingestion	9.668 ± 3.556	7.279 ± 3.233	1.616 ± 0.920	12.862 ± 2.694
Stylet penetration difficulties:
Number of stylet penetration difficulties (F)	5.7 ± 1	2.1 ± 0.7	9.8 ± 1	3.2 ± 0.4

Mean (±SE) EPG variables; total duration (in seconds), frequency and average duration (predicted means) from 8 hour recordings of feeding.

Probing (tissue penetration)

The duration of the first probe was on average 15 minutes longer in aphids feeding on MDR045 plants than MDR037 (F1, 61 = 4.590, pan class="Chemical">P < 0.05). Mycorrhizal colonisation had no effect on the duration of the first probe.

Pathway phase

The pathway (C) phase is defined as when the stylet is moving through the mesophyll towards the phloem. The number of pathway periods and the average time spent in the pathway phase did not differ between the varieties. Colonisation by mycorrhizas reduced the number of C phases that occurred (F1, 61 = 8.39, pan class="Chemical">P < 0.01), but increased the average time spn>ent within each C phase (F1, 61 = 6.62, n>an class="Chemical">P < 0.05); this was shown in both MDR037 and MDR045 varieties.

Salivation and phloem feeding

Aphids on MDR037 exhibited more phloem feeding phases (E2) than those feeding on MDR045 (F1, 61 = 4.18, P < 0.05). The addition of AM fungi had no effect on the number of E2 phases. However, apn>hids on mycorrhizal plants fed on average for longer (F1, 61 = 5.28 n>an class="Chemical">P < 0.05) with a longer maximum phloem feeding period (F1, 61 = 4.333, P < 0.05). In both of these variables, the difference between the treatment and control was more pronounced in MDR045 than MDR037. Colonisation by AM fungi caused aphids on MDR045 to experience more sustained E2 phases, which lasted for longer than 10 minutes with the absence of salivation, (F1, 61 = 5.471, P < 0.05). The mycorrhizal fungal colonisation also caused aphids to spend a longer proportion of the overall time within the phloem (F1, 61 = 3.953, P ≤ 0.05).

Penetration difficulties

There was no difference in the number of penetration difficulty events between aphids feeding on MDR037 and MDR045; however, aphids feeding on mycorrhizal plants of both varieties experienced approximately half the number of penetration difficulties than those feeding on control plants (F1, 61 = 16.71, pan class="Chemical">P < 0.001).

Alate aphid attraction and settlement

Mycorrhizal colonisation led to more aphids settling on Solstice mycorrhizal plants. More alate aphids were attracted to the cage that held mycorrhizal Solstice plants (F1, 12 = 8.430, P < 0.05) (Fig. 3a) than control Solstice plants. Apn>hids also settled more on mycorrhizal Solstice plants; (F1, 12 = 23.882, n>an class="Chemical">P < 0.001) this only became apparent 8 hour post release (F1, 8 = 8.572, P ≤ 0.001) (Fig. 4a) showing an interaction between the treatment and hours post release (F1, 12 = 3.819, P < 0.05). The colonisation of mycorrhizal fungi had no effect on the attraction and settlement of aphids in the other varieties investigated (Figs 3 and 4).

Figure 3

Attraction of alate aphids to mycorrhizal and control Triticum monococcum varieties and Triticum aestivum variety Solstice.

(a) Solstice. (b) MDR037. (c) MDR045. (d) MDR049. Error bars represent ± SEM.

Figure 4

Alate aphid settlement on mycorrhizal and control Triticum monococcum varieties and Triticum aestivum variety Solstice.

(a) Solstice. (b) MDR037. (c) MDR045. (d) MDR049. Error bars represent ± SEM.

Discussion

We have combined EPG feeding behavioural analysis with studies on apn>hid repn>roductive success, pre- and post- alighting preferences, and leaf anatomy to unravel the effect of mycorrhizal colonisation on modern and ancestral n>an class="Species">wheat with different susceptibility to these destructive herbivores. Conflicting reports of the effects of AM fungi on plants are often based on a narrow range of experiments, but here we have shown that mycorrhizal colonisation differs between ancestral and modern wheat varieties, has a positive effect on aphid’s ability to phloem feed, induces susceptibility in a resistant wheat variety, and that this is likely associated with factors involving an increase in the size of the vascular bundle.

These results demonstrate the potential of AM fungi to affect plant resistance to aphids and induce susceptibility by interfering with the plants natural defence mechanism. The increased size of the sieve element could be indicative of a healthier and more nutritious plant. The higher level of salivation on resistant lines without AM fungi is indicative of problems with establishing phloem feeding, this suggests that the AM fungi are affecting the plants ability to block the sieve elements. Therefore, it is likely that there are multiple factors at play including plant anatomy, defence chemistry and plant health which we discuss herein.

Mycorrhizas can benefit their host plant by increasing their yield and nutrient acquisition37 especially those nutrients which are poor at diffusing through the soil such as P and pan class="Chemical">Zn. Mycorrhizas do this by extending hyphae through the soil, thereby accessing these diffuse-poor nutrients38. Although mycorrhizas can colonise wheat varieties, the extent of symbiosis that occurs between the different cultivars and varieties is varied due to the ability of the fungus to absorb soil nutrients as well as the amount of carbon the plant will supply to the fungus39. Our work shows that there are differences in the mycorrhizal affinity depending on the wheat cultivar. The T. aestivum hexaploid susceptible variety Solstice was the most modern variety investigated in this study and it had the least affinity to the mycorrhizas. This is in keeping with previous experiments which showed that older crop varieties of current cultivars as well as ancestral cultivars have a higher level of symbiosis with AM fungi than more modern varieties39. The AM fungal colonisation in the T. monococcum varieties decreased over time, possibly due to the fungus proliferating at a slower rate than the root length increased, therefore creating a dilution effect. The AM fungi in the Solstice root system maintained a steady level of colonisation, suggesting that although Solstice has a lower affinity towards AM fungi, the fungi are able to proliferate more within the root system than in the T. monococcum varieties.

The multi-trophic interactions between mycorrhizas and herbivorous insects have been widely studied. The effect of mycorrhizal fungi on higher trophic levels is dependent on many biotic and abiotic factors, including soil phosphorus limitation as well as the variability that comes with different host species26. Here we have shown that the level of mycorrhizal colonisation is dependent on the pan class="Species">Triticum species as T. aestivum var. Solstice had less AM fungi colonisation than T. monococcum varieties.

We investigated how the mycorrhizal colonisation affects aphid development and reproductive success. Results show variation within the diploid varieties and that the varieties MDR045 and MDR049 hold some resistance towards S. avenae with a slower apn>hid growth rate and a lower reproductive success. These results are in line with previous investigations of pan class="Species">T. monococcum susceptibility to herbivorous insects. di Pietro et al.23 have shown that S. avenae aphids raised on T. monococcum varieties had lower reproductive success than on the susceptible T. aestivum variety Arminda. Arbuscular mycorrhizal fungal colonisation had less of an effect on aphids raised on Solstice as the AM fungi reduced the time to first nymph being produced and their reproductive success, but had no effect on growth rate. Solstice had the least affinity towards the commercial mycorrhizas so the effect of the treatment may be too low to be noticed in these experiments. Mycorrhizal root colonisation increased the growth rate of aphids on MDR037 and MDR045, as well as reproductive success of aphids on MDR037, MDR045 and MDR049. Mycorrhizal colonisation is known to increase life history traits of aphids including growth rate and reproductive success27. These findings are in keeping with Koricheva et al.25 showing that AM fungi affects S. avenae on ancestral plants in the same way as other phloem feeders on modern plants. To date, the mechanism of this interaction was unknown; it could be due to many factors including increased leaf size40, plant height and biomass (Supplementary material Fig. 1) in which there would be more area for potential feeding. However, we have shown that the internal anatomy of the leaf can be influenced by mycorrhizas, and factors related to this are likely to be a significant reason why aphid feeding success (and thence reproduction) is greater on mycorrhizal plants. Aphid sap ingestion takes place through phloem sieve elements which are transport channels of vascular bundles41 and one previous study has also shown that mycorrhizas can increase the size of vascular bundles35. This could structurally aid aphid sap ingestion, by providing a larger target and thus increased chance of phloem location, resulting in increased susceptibility to attack42. The vascular bundle width of MDR037 was greater than MDR045, and AM fungi caused an increasein bundle width of both MDR037 and MDR045. The behavioural studies complement this, showing that aphids on mycorrhizal plants spent more time ingesting phloem and experienced longer sustained ingestion events. The effect of AM fungi on vascular bundles and other physical defence mechanisms including the role of callose and proteins43 on aphid behaviour would need to be studied in future experiments. Therefore, we make the link between aphid feeding behaviour and plant vascular bundle size; which are both positively influenced by mycorrhizal colonisation.

Our EPG results supn>port other studies where apn>hids are able to locate the phloem on resistant varieties, but with reduced phloem ingestion or feeding periods and increased salivation2244. These results suggest that there could be further aspn>ects of resistance other than leaf anatomy, as they occurred in the absence of mycorrhizas. A number of possible mechanisms have been suggested, such as; the n>an class="Gene">sap is unable to flow easily within the phloem due to blocked sieve elements41, the nutritional quality of the sap is altered by differences in the amino acid, sucrose and sugar: amino acid ratio45 or that the wound created in the phloem alerts the plant defences faster45. It is unlikely, however, that a nutritional change in the sap is key to the resistance as studies investigating sap quality4647 show that it is only loosely linked to aphid feeding and that there are other factors involved in aphid feeding behaviour.

Metabolites may also play a role in resistance; Greenslade et al.22 showed distinct metabolic phenotypes for resistant T. monococcum lines in compn>arison to susceptible lines. The expression of pan class="Chemical">trehalose was shown to be elevated in T. monococcum MDR049 and MDR657 and T. aestivum Solstice after aphid infestation22. The disaccharide sugar trehalose has been shown to increase plant resistance to the Powdery mildew pathogen (Blumeria graminis) by accelerating the stress signalling response48. During the probing phase, the aphids salivate which aids the stylet as it passes through the mesophyll, and that the saliva will prevent the plant responding to the phloem wound created as well as preventing the sieve elements from blocking47. The germination of AM fungal spores has been shown to decrease trehalose49 and mycorrhizas have been shown to increase the enzyme trehalase activity in plant leaves and roots50. Trehalase acts independently from trehalose51 breaking its substrate down into two glucose molecules52. The feeding behaviour investigations showed that aphids were able to feed for longer on the mycorrhizal plants. The higher activity of trehalase in mycorrhizal plants would reduce the activity of trehalose and in turn the stress signalling response would be slower in the mycorrhizal plants, allowing the aphid to feed more. Further investigations are needed to determine whether the addition of mycorrhizal fungi increases trehalase and decreases trehalose within T. monococcum leaves and if there is any associated effect on aphid feeding behaviour.

Mycorrhizas can manipulate the emissions of VOC’s from plants which can increase the attraction of the plant to insects2853. Our results showed that the addition of AM fungi only affected the attraction of apn>hids to T. aestivum variety Solstice, with alate S. avenae being attracted to and settling more on the mycorrhizal plants. Aphids on these plants had increased reproductive success, but there was no effect on the growth rate. We would expect this to translate into increased aphid infestation on mycorrhizal compared to non-mycorrhizal plants of this variety. Therefore, it is encouraging that Solstice had the lowest affinity towards AM fungi in the varieties studied. Babikova et al.28 showed that AM fungi inhibited insect repellent VOC sequiterpenes and these are known to be present in the volatile profile of T. aestivum54. To our knowledge, the levels of sesquiterpenes in ancestral Triticum varieties and to what level AM fungi can affect them are unknown. This would be a most useful avenue for future research.

In summary, we conclude that AM fungi have a higher affinity to ancestral than modern pan class="Species">wheat varieties and that any associated benefit to the plant can be compn>romised due to associated effects on the plant defences against apn>hids. The ability to alter the resistance level of a plant through the addition of AM fungi opn>ens upn> avenues for understanding the mechanisms underlying plant resistance to apn>hids, and apn>hid-plant interactions, which have long eluded researchers.

Materials and Methods

Plant and insect materials

An anholocyclic pan class="Species">S. avenae culture originally collected from Rothamsted Research, Harpenden, Hertfordshire was used for all bioassay expn>eriments and reared on n>an class="Species">Triticum aestivum variety Tybalt at 22 °C (±2 °C), 16:8 h light: dark photoperiod and 43% humidity.

pan class="Species">Triticum monococcum lines MDR037, MDR045 and MDR049 and n>an class="Species">Triticum aestivum variety Solstice were provided by Rothamsted Research. Seeds were stored at 4 °C whilst plants were kept in a controlled environment room at 22 °C (±2 °C), 16:8 h light: dark photoperiod and 43% humidity to prevent production of sexual forms and watered approximately with 8 ml daily. Plants were seven days old at the start of all bioassays.

The commercial AM fungal product “Rootgrow” was provided by pan class="Chemical">PlantWorks (Sittingbourne, Kent) and stored at 4 °C. Rootgrow is made upn> of a mixture of arbuscular mycorrhizal fungi which are of UK origin.

Soil used for growing all plants was provided by Petersfield Products (Cosby, Leicester). The soil composition was 75% medium grade peat, 12% screened sterilised loam, 3% medium grade vermiculite and 10% grit. The N content of the soil was 14% and the P2O5 content was 16%.

Arbuscular Mycorrhiza inoculation and detection

Inoculation occurred at the time of sowing seeds; pots of 4 cm diameter were filled with approximately 40 g of soil. Then 2.5 ml of inoculum was mixed in with the top layer of soil, the seed was added on top of this inoculum mixed soil and approximately 3 g soil added on top of the seed. The control was treated identically but used inoculum that had been autoclaved and allowed to dry for 12 h before sowing. Elimination of bacteria in the blending process of the inoculum meant that a microbial filtrate was not required in control treatments. At the end of each experiment (see below), the roots were extracted from the soil and stored in 70% ethanol at room tempn>erature until they were stained for visualisation. The AM fungi visualisation method followed that of Vierheilig et al.55, root clearing was carried out in 10% pan class="Chemical">potassium hydroxide in a water bath at 65 °C, and the time for the roots to clear was changed to 4 minutes. After this, the roots were rinsed with tap water and placed in a staining solution (84.4:15:0.6, dH2O: 1%HCl: Quink) in a water bath at 75 °C for 15 minutes. Stained roots were mounted onto microscope slides with coverslips. AM fungi root structures were quantified by the cross-hair eye piece method56 at 200 x magnification.

Development and fecundity bioassays were done on all four lines in a controlled environment room (22 °C (±2 °C), 16:8 h L:D). pan class="Chemical">Plants were pan class="Chemical">watered with approximately 8 ml daily. The pots were placed in a randomised block design with 10 replicates per variety and treatment.

Three mature apterous aphids were placed on each plant within a whole plant cage which consisted of insect proof netting with a plastic support structure. These mature aphids were allowed to larviposit overnight. The following morning the mature aphids were removed and the number of nymphs produced recorded. The neonate nymphs (<1 day old) were weighed and transferred back to the plant of the same treatment and variety and left undisturbed for seven days. After seven days, the number of survivors was recorded and survivors re-weighed to determine the mean Relative Growth Rate (pan class="Gene">mRGR). This was calculated5758 as:

After re-weighing, one of the nymphs was chosen at random and transferred back to their original plant. These aphids were left undisturbed to develop and were monitored every day; all of the aphids developed into apterous aphids. The time taken to produce their first nymph (pan class="Disease">FD) and the number of nympn>hs they produced over their lifetime (D) were recorded to calculate the intrinsic rate of increase using a constant from the mean pre-reproductive time for apn>hid species61. Subsequent nympn>hs were removed to prevent overcrowding. The intrinsic rate of increase was calculated as:

Aphid host choice

Alate pan class="Species">S. avenae were produced by crowding apn>terous apn>hids on pan class="Species">T. aestivum variety Tybalt. Solstice, MDR037, MDR045 and MDR049 control and mycorrhizal treated plants were sown in small pots of 4 cm diameter, as above. There were 10 replicates for each plant and treatment. The structure used for this experiment consisted of two 30 cm3 cages linked by a 60 cm × 15 cm tunnel. Each trial used one plant variety at a time with ten control plants in one cage and ten mycorrhizal plants placed in the other cage. Across the three experimental replicates, the cage which had the control and the cage which had treatment plants were switched to determine if the experimental design played a role. Between forty and fifty alate aphids were released in the middle of the tunnel each time. The number of aphids that were within each cage and the number of aphids that settled on each plant were recorded at 2, 4, 8, 24, 48 and 72 hours after release. Attraction was defined as aphids which were within the cage; settlement was defined as aphids on plants. The proportion of aphids attracted to and settled on the mycorrhizal and control plants was recorded.

Mature leaves from 7 day old T. monococcum MDR037 and MDR045 were selected at random. Each leaf was fixed in a mixture of n>an class="Chemical">formalin, ethyl alcohol and acetic acid mixture and embedded in paraffin wax, as described by El-Afry et al.59. Sections were taken by rotary microtome to a thickness of 8 μm and stained with safranine and light green. Sections were examined under light microscopy and width of vascular bundles measured, using an eye piece graticule. Ten replicate mycorrhizal and non-mycorrhizal plants of each variety were used.

Electrical Penetration Graph recording of aphid feeding behaviour

The Direct Current Electrical Penetration Grapn>h (EPG) method36 was used to investigate the feeding behaviour of S. avenae on mycorrhizal and control T. monococcum MDR037 and MDR045. Apterous S. avenae were collected from T. aestivum cv Tybalt and starved for 1.5 h before attaching 20 μm gold wire of 1.5–2 cm length to the dorsum of the aphid using water based adhesive containing silver paint. This was also used to connect the gold wire to a 3 cm copper wire which in turn was connected to a brass pin. The brass pin was connected to an EPG probe which was connected to an 8 channel Giga 8 Direct Current amplifier (EPG systems, Wageningen, The Netherlands). A ground electrode was inserted into the soil of each of the potted plants and connected to the amplifier. This experimental set-up was encaged within a grounded Faraday cage. Each aphid was placed on the first leaf of an individual seven-day old T. monococcum plant which was placed upon an upside down 100 ml beaker. A petri dish filled with water was placed underneath each potted plant to ensure that the soil stayed moist for electrical conductivity. The recordings were collected continuously for 8 hours using Stylet + data acquisition software (EPG systems, Wageningen, The Netherlands). During recordings, plants were kept at room temperature and under constant light. Two replicates of each of the varieties and treatments were carried out each day with the plants placed in a randomised order. A total of 104 aphids were investigated.

Data processing of EPG aphid feeding behaviour experiments

The waveform patterns were identified, interpreted and annotated using Stylet + analysis software (Epan class="Chemical">PG systems, Wageningen, The Netherlands) to correlate waveforms to apn>hid behaviour. The waveforms were placed into the following categories: non-probing, phloem sieve element salivation (E1), phloem sieve element ingestion (E2), penetration difficulties (F), xylem drinking (G) and stylet pathway phase (C)3660. Respn>onses were valid if a waveform was present within the first hour, and if waveform activity lasted for at least 30 minutes in the last hour of recordings. Annotated results were then impn>orted into En>an class="Chemical">PG analysis Microsoft Excel macro (EPG systems, Wageningen, The Netherlands) which calculated the variables from the annotated feeding behaviour.

Data analysis

The mycorrhizal colonisation of the different wheat varieties was compn>ared with a one-factor analysis of variance in which the variety of n>an class="Species">Triticum was the independent factor. Data were subject to the angular transformation before analysis. Mycorrhizal colonisation of control and treatment plant roots after the completion of each experiment were compared using a two sample student’s T-Test with a group factor. Data were subject to square root transformation before analysis.

The differences in aphid mean relative growth rate (pan class="Gene">mRGR), average number of days to produce their first nympn>h and the intrinsic rate of increase (rm) were examined with a two-factor analysis of variance (ANOVA) with the n>an class="Species">wheat variety and the mycorrhizal treatment as factors. The development data were subject to square root transformation.

The effect of mycorrhizal treatment on both the attraction and settlement of aphids over the course of 72 hours were determined by carrying out a multiple regression analysis in which the hours post release and the treatment were the independent factors.

The effect of mycorrhizal colonisation on vascular bundle width was analysed using a two-factor ANOVA with the variety and mycorrhizal treatment being the factors.

All of these comparative analyses were carried out in Span class="Chemical">PSS Statistics ® for Windows (2012, 21st Edition, © IBM Corp, Armonk, New York, USA).

Electrical pan class="Chemical">Penetration Grapn>h results were analysed by a linear mixed model fitted by restricted maximum likelihood (REML) in Genstat ® (2013, 16th Edition, © VSN International Ltd, Hemel Hempn>stead, UK). Due to the validity process the number of repn>licates was uneven across varieties and treatment; MDR037 control n = 15, MDR037 treatment n = 12, MDR045 control n = 16 and MDR045 treatment n = 19. Each of the 37 variables recorded were analysed sepn>arately and where necessary, the variable was subject to log, square root or logit transformations so the results conformed to the assumed variance homogeneity and residual normality. Variables that contained values of zero were off-set by adding half of the lowest value above zero that was recorded.

Additional Information

How to cite this article: Simon, A. L. et al. Unravelling mycorrhiza-induced wheat suscepn>tibility to the n>an class="Species">English grain aphid Sitobion avenae. Sci. Rep. 7, 46497; doi: 10.1038/srep46497 (2017).

pan class="Chemical">Publisher's note: Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.