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

OPEN

# Plant health status effects on arbuscular mycorrhizal fungi associated with *Lavandula angustifolia* and *Lavandula intermedia* infected by Phytoplasma in France

Marie-Noëlle Binet<fnref n="1" />, Camille Marchal<fnref n="1" />, Justine Lipuma², Roberto A. Geremia<fnref n="1" />, Olivier Bagarri³, Bert Candaele⁴, Delphine Fraty⁴, Benjamin David², Sophie Perigon<fnref n="1" />, Viviane Barbreau⁵ & Bello Mouhamadou<fnref n="1" />✉

We investigated root communities of arbuscular mycorrhizal fungi (AMF) in relation to lavender (*Lavandula angustifolia*) and lavandin (*Lavandula intermedia*) health status from organic and conventional fields affected by Phytoplasma infection. The intensity of root mycorrhizal colonization was significantly different between diseased and healthy plants and was higher in the latter regardless of agricultural practice. This difference was more pronounced in lavender. The root AMF diversity was influenced by the plant health status solely in lavender and only under the conventional practice resulting in an increase in the AMF abundance and richness. The plant health status did not influence the distribution of root AMF communities in lavandin unlike its strong impact in lavender in both agricultural practices. Finally, among the most abundant molecular operational taxonomic units (MOTUs), four different MOTUs for each plant species were significantly abundant in the roots of healthy lavender and lavandin in either agricultural practice. Our study demonstrated that the plant health status influences root colonization and can influence the diversity and distribution of root AMF communities. Its effects vary according to plant species, can be modified by agricultural practices and allow plants to establish symbiosis with specific AMF species.

Arbuscular mycorrhizal fungi (AMF, Glomeromycota) are obligate symbiotic fungi which colonize roots of the majority of land plants, including the most important food crops<fnref n="1" />. They provide essential ecosystem services by improving plant nutrient availability, plant diversity and soil structure formation². They also enhance plant resistance to abiotic stresses³ as well as to biotic stresses⁴,⁵. Many studies have shown the important role of AMF against a wide range of phytopathogens resulting in a reduction or suppression of plant disease symptoms⁶. These remarkable properties make the AMF key players in cropping systems to increase crop yields and plant defenses within the framework of sustainable agriculture. Studying the diversity and the drivers shaping AMF communities in roots of agricultural crops could help developing strategies involving these fungi in relation to plant health⁷.

Although scarcely studied, plant health status seems to play an important role on AMF community structure and diversity. For instance, clear relationship between root-associated fungal communities' composition and plant health status in *Pisum sativum* has been reported⁷. They found that AMF, described as able to strengthen the defense of plants against phytopathogenic agents or improving their fitness, were more abundant and diverse

<footnote n="1">¹Laboratoire d'Écologie Alpine, CNRS UMR 5553, Université Grenoble Alpes, CS 40700, 38058 Grenoble Cedex 09, France. ²MYCOPHYTO, 50 Avenue de la Plaine, 06250 Mougins, France. ³Université Européenne des Senteurs et des Saveurs, Couvent des Cordeliers, 04300 Forcalquier, France. ⁴Centre Régionalisé Interprofessionnel d'Expérimentation en Plantes Aromatiques et Médicinales, Lieu-dit Les Quintrands, Route de Volx, 04100 Manosque, France. ⁵Collège Henri Wallon, 17 rue Henri Wallon, 38400 Saint Martin d'Hères, France. ✉email: bello.mouhamadou@univ-grenoble-alpes.fr</footnote>
in healthy pea plants. Likewise, an important shift in diversity and composition of root AMF communities when comparing healthy to diseased plants in Salix bioenergy plantations was shown⁸. In both studies, the effects of plant health status were investigated on sites characterized by the same mode of agricultural practices. To our knowledge there is no information on the importance of the plant health status on root AMF communities associated to plant crops under different agricultural systems such as conventional or organic practices. Both practices are often increasingly developed to improve crop yields and are also reported to affect AMF communities⁹. For example, organic farming systems via the crop rotation, the cover cropping or a lack of synthetic pesticides and inorganic fertilizers promote root colonization and abundance of spores in root systems maintaining alongside a greater diversity of AMF communities in the soil¹⁰. Conversely, in conventional farming the negative effects of high input cropping systems or conventional flooding on AMF communities have been reported due probably to the disturbance of the hyphal networks they provoke¹¹,¹². However, since conflicting results have also been described, showing very weak effects or even the lack of effects of agricultural practices on the community composition of AMF in roots¹³,¹⁴, we suggested that possible interactive effects with other biotic or abiotic factors such as the host plant health status could exist, thereby impacting AMF communities.

The purpose of this article was to study the responses of AMF to plant health status on two aromatic plants, *L. angustifolia* and *L. intermedia* cultivated under different agricultural practices in south-eastern France. In these lavender agricultural systems, the yellow decline of lavender caused by the stolbur phytoplasma '*Candidatus Phytoplasma solani*' reduced seriously plant productivity¹⁵. Otherwise, different agricultural practices such as conventional or organic farming were used for both *L. angustifolia* and *L. intermedia* so as to find a strategy to reduce the impact of the disease. In both practices, diseased and healthy plants coexisted. The latter are likely the result of either the absence of contact with phytoplasma or biotic interactions such as the establishment of a mycorrhizal symbiosis with specific AMF communities able to suppress plant damage as reported elsewhere¹⁶⁻¹⁹. By studying the root AMF communities associated to the healthy and diseased lavender and lavandin, we tested the hypothesis that the AMF colonization, the diversity and the composition of root AMF communities could be influenced by the plant health status. We also sought to determine whether these plant health status effects varied according to the agricultural practice and to the plant species identity.

## Results

**AMF root colonization.** In lavender, whatever the agricultural practice, the mycorrhizal root colonization was significantly higher in the healthy plant compared to the diseased plant. It was 27% and 15% in the healthy and diseased plants, respectively, in organic practice and it was 30% and 5% in the healthy and diseased plants, respectively, in conventional practice (Fig. 1a,b). Variance partitioning indicated that the plant health status explained about 30% of the variability of the mycorrhizal root colonization in organic or conventional practices (Table S1). Similar patterns were observed in lavandin roots with a higher root colonization in the healthy plant compared to the diseased plant. In organic practice, it was 42% and 30% in the healthy and diseased plants, respectively and in conventional practice, it was 12% and 7% in the healthy plant and the diseased plant respectively. (Fig. 1c,d). However, the contribution of the plant health status in the variability in root colonization was lower compared to that observed in lavender roots and accounted approximately for 15% in each of the agricultural practices (Table S1).

**AMF diversity.** First, rarefaction analysis performed on the data indicated that the taxon accumulation curves reached saturation (Fig. S1). The AMF diversity was studied by focusing on the abundance and richness of fungal communities. Regarding the abundance of AMFs (corresponding to the abundance of number of sequences of each OTU), the results were different depending on the agricultural practice in lavender. We found an effect of the plant health status only in conventional practice which was reflected by a higher abundance in the healthy plants compared to the diseased plant (Fig. 2a,b). In lavandin, however, whatever the agricultural practice, the plant health status had no effect on the abundance of AMF (Fig. 2c,d).

Regarding the AMF richness, we observed a plant health status effect only in lavender in conventional practice. A greater AMF richness was observed in healthy plants compared to diseased plants (Fig. 3a–d).

**AMF distribution.** In lavender roots, multiple factor analysis of the AM fungal communities, based on the abundance of AMF OTUs, identified two distinct AMF groups determined by the plant health status (Fig. 4a,b). In both agricultural practices, the first axis, which clearly separated the fungal communities according to the plant health status, accounted for 34.8% and 42% of the variation in organic and conventional practices, respectively. The second axis, which did not represent more than 16% of the variation, less clearly separated the AMF communities according to the plant health status. Permanova analysis indicated that the plant health status explained significant portion of the variance in fungal distribution (Table S2), 16% and 19% in organic and conventional practices, respectively. The other factors (abiotic or biotic) which influence AMF communities were represented by “residuals” and accounted for 84% and 81% in organic and conventional practices, respectively. By contrast to lavender, the plant health status had no significant effect on the structure of AMF communities in the lavandin roots whatever the agricultural practice (Fig. 4c,d and Table S2).

**AMF community composition.** The global composition of AMF communities is displayed in Fig. S2. Globally, Glomeraceae dominated lavender and lavandin roots followed by unclassified Glomeromycetes. Cladoideoglomeraceae, Acaulosporaceae and Gigasporaceae were present in small proportion and their distribution was variable depending on the plant species, the plant health status or the type of agricultural practice (Fig. S2).

The distribution of the most abundant AMF OTU i.e. present in at least 2/3 of the samples was also studied. First, their sequences, deposited in Genbank (MN510817-MN510849), were analyzed by phylogeny (Fig. 5)
<figure><img src="image_2.png" /></figure>

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

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

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Arbuscular mycorrhizal fungal colonization in healthy and diseased plants of *L. angustifolia* (a, b) and *L. intermedia* (c, d) under different agricultural practices. Different letters indicate significant differences (Test Wilcoxon; *P* value < 0.05) according to the health status.</figcaption></figure>

using the reference sequences available in the Genbank databases to confirm the taxonomic assignments and validate the specific delimitation of each MOTU. The most abundant AMF MOTUs found for each plant species were related to the plant health status in each agricultural practice. For example, in the lavender roots (Fig. 6a), 21 MOTUs were found with a variable distribution according to the plant health status in either agricultural practice. Among them, the most abundant MOTUs were *Rhizophagus intraradices* and *Glomus* sp2 in diseased lavender and *Glomus* sp5 and *Glomus* sp7 in healthy lavender under conventional practice on the one hand and *Glomus* sp2, *Glomus* sp4, *Glomus* sp6 and *Glomus* sp8 in diseased lavender and *Glomus macrocarpum* and *Funneliformis mosseae* in healthy lavender under organic practice on the other hand.

Likewise, in lavandin roots (Fig. 6b), 23 abundant MOTUs were found with variable patterns according to the plant health status in either agricultural practice. Among them, one MOTU (*Rhizophagus irregularis*) and 3 MOTUs (*Glomus aggregatum*, *Glomus* sp2 and *Rhizophagus intraradices*) were statistically more abundant in diseased and healthy lavandins respectively under conventional practice whereas 2 MOTUs (*Glomus aggregatum* and *Rhizophagus* sp5) and one MOTU (*Rhizophagus* sp1) were statistically more abundant in diseased and healthy lavandins respectively in organic practice.
<figure><img src="image_6.png" /><figcaption>**Figure 2.** Abundance of AMF community in roots of healthy and diseased plants of *L. angustifolia* (a, b) and *L. intermedia* (c, d) under different agricultural practices. Different letters indicate significant differences according to the plant health state (Test Wilcoxon; *P* value < 0.05).</figcaption></figure>

**Discussion**

Till now, studies that have investigated the plant health status effects on root AMF communities have been carried out on sites characterized by one agricultural management type⁷,⁸. These studies do not take into account the importance of different agricultural practices which are increasingly used to enhance crop yield and which are moreover described to influence structure and diversity of root AMF communities⁹⁻¹³. To bridge this gap, we investigated the specific effects of the plant health status of two Mediterranean crop plants, lavender and lavandin on the AMF colonization intensity and the structure of root AMF communities from organic and conventional fields affected by Phytoplasma infection. We assessed whether these effects differed according to the agricultural practices and were specific to each plant species.

We found an influence of the plant health status of both plant species on the intensity of root mycorrhizal colonization resulting in the establishment of a more effective symbiosis in the healthy plant compared to the
<figure><img src="image_7.png" /><figcaption>**Figure 3.** Richness of AMF community in roots of healthy and diseased plants of *L. angustifolia* (a, b) and *L. intermedia* (c, d) under different agricultural practices. Different letters indicate significant differences according to the plant health state (Test Wilcoxon; *P* value < 0.05).</figcaption></figure>

diseased plant as demonstrated elsewhere²⁰,²¹. Several factors could explain the implementation of symbiosis between AMF and host plants such as the root exudates²² or the ability of the host plant to supply C to AMF¹. We suggested that lavender and lavandin produce root exudates whose quality and quantity, specific to each species, could vary according to the plant health status thus impacting AM symbiosis. It is also possible that diseased lavender and lavandin could reduce the C supply necessary for AMF development due to the modification in carbohydrate synthesis and transport as shown in phytoplasma-infected plants²³,²⁴. Plant health status effect was evidenced in both types of agricultural practices, however the fact that it was significantly greater in lavender than in lavandin (30% vs 15%, Table S1) could suggest the influence of the host plant identity on root AMF colonization²⁵⁻²⁷.

The mycorrhizal root colonization may be the reflection of the diversity of root AMF in host plants. Interestingly, we found that the plant health status impacted significantly the root AMF diversity by increasing the AMF abundance and richness in healthy lavender. However, this result, consistent with prior studies⁷,⁸, seemed
<figure><img src="image_8.png" /><figcaption>**Figure 4.** Multiple factor analysis of root AMF communities associated to *L. angustifolia* and *L. intermedia* roots. (a, b) Corresponded to *L. angustifolia* in organic and conventional practices, respectively; (c, d) corresponded to *L. intermedia* in organic and conventional practices, respectively.</figcaption></figure>

to depend on the type of agricultural practices as it was seen only under the conventional practice whereas it was lacking in organic practice. We suggested that the plant health status and agricultural practices produce interactive effects on the diversity of root AMF communities in lavender. In lavandin however, consistent with the weaker effect of the plant health status on the root colonization, no effect of the plant health status was observed whatever the mode of agricultural practice. These contrasting results between lavender and lavandin may be related to the specificity of each species resulting probably from differences in root exudation patterns.

The survey of the distribution of root AMF communities partially confirms the patterns observed on fungal diversity, in particular the absence of plant health status effect on AMF communities associated with lavandin but its strong impact on the AMF associated with lavender as reported⁸. There are certainly other drivers, not measured in this study, which affect AMF communities including probably the agricultural practice, but in lavander, we did not detect the difference of the plant health status effect according to the agricultural practice as
<figure><img src="image_9.png" /><figcaption>**Figure 5.** Phylogenetic analyses performed with the neighbour-joining method based on aligned sequences of 28S ribosomal RNA gene. Analyses were performed as described⁴⁰. The AMF species from which the sequences (Genbank accession number in brackets) have been recovered from GenBank and used as reference sequences are indicated in red. The Genbank accession numbers from this work are underlined. Bootstrap values exceeding 50% are shown on the branches.</figcaption></figure>

observed in the case of the diversity study. However, this result is not unusual as the sensitivity of the distribution or the diversity of AMF toward different drivers is not necessarily the same as reported elsewhere²⁸.
(a)

<table><thead><tr><th></th><th>Conventional diseased</th><th>Conventional healthy</th><th>Organic diseased</th><th>Organic healthy</th></tr></thead><tbody><tr><td>Rhizophagus intraradices</td><td>●</td><td>●</td><td>●</td><td>●</td></tr><tr><td>Rhizophagus sp2</td><td></td><td></td><td>●</td><td>●</td></tr><tr><td>Rhizophagus sp1</td><td>●</td><td>●</td><td>●</td><td>●</td></tr><tr><td>Rhizophagus sp3</td><td>●</td><td>●</td><td>●</td><td>●</td></tr><tr><td>Glomus sp1</td><td></td><td></td><td></td><td>-</td></tr><tr><td>Glomus sp2</td><td>●</td><td>●</td><td>●</td><td>-</td></tr><tr><td>Glomus sp3</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Glomus sp4</td><td></td><td></td><td>●</td><td>-</td></tr><tr><td>Glomus sp5</td><td></td><td>●</td><td></td><td>-</td></tr><tr><td>Glomus sp6</td><td></td><td></td><td>●</td><td>-</td></tr><tr><td>Glomus sp7</td><td></td><td>●</td><td></td><td>-</td></tr><tr><td>Glomus sp8</td><td></td><td></td><td>●</td><td>-</td></tr><tr><td>Glomus macrocarpum</td><td></td><td></td><td>-</td><td>●</td></tr><tr><td>Glomus sp9</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Glomus sp10</td><td></td><td>-</td><td>-</td><td>●</td></tr><tr><td>Rhizophagus sp6</td><td>●</td><td>●</td><td></td><td>-</td></tr><tr><td>Rhizophagus sp4</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Septoglomus sp2</td><td>●</td><td>●</td><td></td><td>-</td></tr><tr><td>Rhizophagus sp8</td><td>●</td><td>-</td><td>-</td><td>-</td></tr><tr><td>Rhizophagus sp5</td><td></td><td></td><td>●</td><td>-</td></tr><tr><td>Funneliformis mosseae</td><td></td><td></td><td>-</td><td>●</td></tr></tbody></table>

(b)

<table><caption>**Figure 6.** Composition of the most abundant AMF MOTUS in healthy and diseased plants of *L. angustifolia* (a) and *L. intermedia* (b) under different agricultural practices. The proportion of each MOTU (represented by circle) was compared between diseased and healthy conditions using a Chi-squared test (for conventional and organic agriculture, respectively).</caption><thead><tr><th></th><th>Conventional diseased</th><th>Conventional healthy</th><th>Organic diseased</th><th>Organic healthy</th></tr></thead><tbody><tr><td>Septoglomus sp1</td><td>●</td><td>●</td><td>-</td><td>-</td></tr><tr><td>Septoglomus jasnowskae</td><td></td><td></td><td>●</td><td>-</td></tr><tr><td>Glomus aggregatum</td><td></td><td>●</td><td>●</td><td>-</td></tr><tr><td>Septoglomus viscosum</td><td></td><td></td><td>●</td><td>●</td></tr><tr><td>Rhizophagus sp13</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Rhizophagus intraradices</td><td></td><td>●</td><td>●</td><td>-</td></tr><tr><td>Rhizophagus sp1</td><td></td><td></td><td>-</td><td>●</td></tr><tr><td>Rhizophagus sp3</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Glomus sp2</td><td></td><td>●</td><td></td><td>-</td></tr><tr><td>Rhizophagus sp11</td><td>●</td><td>●</td><td>-</td><td>●</td></tr><tr><td>Glomus sp7</td><td></td><td></td><td>●</td><td>●</td></tr><tr><td>Glomus macrocarpum</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Glomus sp10</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Rhizophagus sp6</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Claroideoglomus sp1</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Claroideoglomus sp2</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Claroideoglomus sp3</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Claroideoglomus sp4</td><td></td><td></td><td>-</td><td>-</td></tr><tr><td>Rhizophagus sp4</td><td>●</td><td>●</td><td>●</td><td>●</td></tr><tr><td>Rhizophagus sp7</td><td>●</td><td>●</td><td>●</td><td>●</td></tr><tr><td>Rhizophagus irregularis</td><td>●</td><td>●</td><td>●</td><td>●</td></tr><tr><td>Rhizophagus sp5</td><td></td><td></td><td>●</td><td>●</td></tr><tr><td>Funneliformis mosseae</td><td></td><td></td><td>-</td><td>-</td></tr></tbody></table>

The plant health status effects were also analyzed for the most abundant MOTUs present in at least 2/3 of the samples. Only MOTUs belonging to Glomeraceae and Claroideoglomeraceae families met these criteria. The plant health status effects were reflected more by differences in the MOTU abundance rather than the selection
<table><caption>**Table 1.** Description of management practices in organic and conventional fields for each plant species, *L. angustifolia* and *L. intermedia* and soil characteristics. The herbicides and insecticides correspond to propyzamid, cycloxydim, fluroxypyr or amidosulfuron and chlorpyriphos-methyl for pest protection. The reduced tillage was performed on 5–10 cm depth. For each plant species, means with the same letter within columns are not significantly different (t-test; *P* < 0.05). Organic matter, total Carbone and pH were measured from sampled soils for each site by Aurea, France according to standard methods.</caption><thead><tr><th>Plant species</th><th>Management</th><th>GPS</th><th>Fertilizer applied</th><th>Pesticide use</th><th>Inter-row management</th><th>Organic matter (%)</th><th>Total carbone (%)</th><th>pH</th></tr></thead><tbody><tr><td rowspan="2">L. intermedia</td><td>Organic for 20 years</td><td>43° 54' 5.15" N<br>6° 11' 3.34" E</td><td>None</td><td>None</td><td>No tillage Plant<br>cover (Fabaceae and<br>Poaceae)</td><td>3.20 ± 0.10<sup>a</sup></td><td>1.84 ± 0.01<sup>a</sup></td><td>8.25 ± 0.05</td></tr><tr><td>Conventional > 20 years</td><td>43° 52' 23.2" N</td><td>30:30:40 kg NPK<br>ha<sup>-1</sup></td><td>Herbicides and<br>insecticides</td><td>Reduced tillage</td><td>2.70 ± 0.00<sup>a</sup></td><td>1.57 ± 0.01<sup>a</sup></td><td>8.25 ± 0.05</td></tr><tr><td rowspan="2">L. angustifolia</td><td>Organic for 4 years</td><td>44° 7' 46.45" N<br>5° 29' 25.98" E</td><td>None</td><td>None</td><td>Reduced tillage</td><td>5 ± 0.20<sup>a</sup></td><td>2.90 ± 0.09<sup>a</sup></td><td>8.2 ± 0.00<sup>a</sup></td></tr><tr><td>Conventional > 20 years</td><td>44° 8' 8.12" N<br>5° 28' 1.38" E</td><td>30:30:40 kg NPK<br>ha<sup>-1</sup></td><td>Herbicides and<br>insecticides</td><td>Reduced tillage</td><td>4.2 ± 0.30<sup>a</sup></td><td>2.46 ± 0.18<sup>a</sup></td><td>8.25 ± 0.05</td></tr></tbody></table>

of specific MOTUs since some MOTUs considered as absent (Fig. 6a,b) are nevertheless present but in less than 1/3 of samples. Second, among the affected MOTUs, four different MOTUs for each species were significantly abundant in the roots of healthy lavender and lavandin in either agricultural practice. More interestingly, among selected MOTUs dominating in healthy plants, some such as *Funneliformis mosseae* or *Glomus macrocarpum* are described as being able to alleviate disease-related symptoms or plant stresses caused by abiotic conditions⁴,²⁹. Plant health status may contribute to the selection of specific MOTUs in the roots of healthy plants, the amplification and use of which as an inoculum could probably help to strengthen or improve the resistance of plants against plant pathogens.

## Conclusions

Taken together, our results show the importance of plant health status on root AMF communities. This factor seems to promote in varying degrees the development of AMF symbiosis in plant roots and influence the diversity and composition of AM fungal communities. It may be affected by agricultural practices and depends on the host plant species identity. It is not excluded to think that the abundant MOTUs in the roots of healthy plants characterized here could contribute to the protection of these plants against Phytoplasma by immunizing the plant to a pathogen as well as other abiotic stresses.

## Materials and methods

**Site description.** The study was conducted in four fields corresponding to two different agricultural practices of two species, the lavender (*L. angustifolia*, cultivar Rapido, seed derived populations) and the sterile hybrid lavandin (*L. intermedia*, cultivar Grosso, vegetatively propagated by cuttings) in the South-East of France. The coverage phytoplasma infections observed in these studied fields reached about 20–30% of the total plants. The description of management practices in organic and conventional fields for each plant species and soil characteristics are shown in Table 1. In both organic and conventional fields affected by Phytoplasma infection, two plant health status coexisted:—the asymptomatic plants defined as “healthy plants” showing no visible disease symptoms and—the symptomatic plants corresponded to “diseased plants” exhibiting typical symptoms of phytoplasma disease with yellowing and reduction of inflorescences on the whole plant¹⁵. These symptoms were examined by professional experts from CRIEPPAM (Manosque, France).

**Experimental approach.** In area of 20 m × 20 m in each management field, in July 2017, we randomly collected 12 roots samples from the underground part of asymptomatic or symptomatic plants aged 2 years. All the resulting 96 roots samples (2 plant health status: healthy and diseased × 2 management types: organic and conventional × 2 plant species: lavender and lavandin × 12 replicates) were thoroughly rinsed with tap water and were used for both the determination of level of mycorrhizal colonization and DNA extraction.

**Level of mycorrhizal root colonization.** Roots were stained as described³⁰. The AMF colonization intensity was estimated under light microscopy³¹ and based on the MYCOCALC program (https://www.dijon.inra.fr/mychintec/Mycocalc-prg/download.html).

**DNA extraction, PCR amplification and illumina-based sequencing.** For each plant species, the 12 root samples were pooled together two by two giving six samples for DNA extraction. The roots have been dried and ground with liquid nitrogen and 100 mg were used to extract DNA with the Fast prep DNA Spin kit (MPBIO, France) according to the manufacturer's instructions. The AMF-specific universal primers FLR3 and FLR4 was used to amplify the partial large subunit region³²,³³. The primers were extended by sample specific tags of 8 nt length to allow parallel sequencing of multiple samples. The PCR reactions were carried out as reported³⁴ by adjusting the annealing temperature (58 °C) and the elongation time (30 s) according to the FLR3 and FLR4
primers. The amplicons were pooled for the sequencing by using equivalent molarities amongst amplicons. The library construction and sequencing (Illumina MiSeq 150-bp pair-end) were carried out by Fasteris (Geneve, Switzerland).

**Bioinformatic analysis.** Reads assembly and primary filtering were done using the OBItools package³⁵. Short (<100 nt) or rare (occurrency <2) reads were removed, resulting in 138,593 unique sequences for 3,917,531 reads in total. Highly similar unique sequences were clustered into Molecular Operational Taxonomic Units (MOTUs) by computing pairwise similarities with Sumatra (OBItools package, https://metabarcoding.org/sumatra) and forming 98% similarity clusters³⁶ with the Markov Cluster algorithm (MCL) classification process³⁷. The ecotag function of OBItools was used to taxonomically assign the 3260 obtained MOTUs to an AMF database. It was generated from the r136 EMBL fungi database with an in silico PCR (function “ecoPCR”, OBItools package) using the FLR3 and FLR4 primers.

The R software 3.5.3 was then used to aggregate unique sequences belonging to the same MOTUs. From these sequences, we removed sequences of rare MOTUs (found in only one sample), sequences of the samples that have a Bray–Curtis scatter distance which exceeds 95% of their confidence interval and sequences of MOTUs not assigned to AMF clades. The dataset was then normalized by sample using Hellinger standardization method («descostand» function, package vegan). The final contingency table displays the relative abundance (normalized read counts) of 190 MOTUs and 180 among the lavender and lavandin samples, respectively.

**Statistical analysis.** Normality of the data was tested by Shapiro-Wilk test. Accordingly, the root colonization intensity (M %) and AMF diversity were analyzed by Levene's test (for the homogeneity of variance) and Student's t-tests or Wilcoxon test. The distribution of AMF communities was evaluated by Multiple Factor Analysis (MFA). The significance of the effect of plant health status on the root colonization intensity, diversity and structure of AMF communities were assessed by Permanova followed by Permutest. For each plant species, the dominant MOTUs (found in at least 2/3 of samples) were selected. Their composition was compared between diseased and healthy conditions—for conventional and organic agriculture respectively—using the Chi-squared test³⁸, function «chisq.test») and Adjusted standardized residuals. All statistical analyses were performed using R package vegan³⁹.

Received: 11 July 2020; Accepted: 5 November 2020

Published online: 20 November 2020

## References

1. Smith, S. E. & Read, D. J. *Mycorrhizal Symbiosis* 3rd edn. (Academic Press, London, 2008).

2. Gianinazzi, S. et al. Agroecology: the key role of arbuscularmycorrhizas in ecosystem services. *Mycorrhiza* **20**, 519–530 (2010).

3. Lenoir, I., Fontaine, J. & Sahraoui, A. L. H. Arbuscularmycorrhizal fungal responses to abiotic stresses: a review. *Phytochem* **123**, 4–15 (2016).

4. Song, Y., Chen, D., Lu, K., Sun, Z. & Zeng, R. Enhanced tomato disease resistance primed by arbuscularmycorrhizal fungus. *Front. Plant Sci.* **6**, 786 (2015).

5. Van Geel, M. et al. Abiotic rather than biotic filtering shapes the arbuscularmycorrhizal fungal communities of European semi-natural grasslands. *New Phytol.* **220**, 1262–1272 (2018).

6. Varma, A., Prasad, R. & Tuteja, N. *Mycorrhiza—Nutrient Uptake* (Biocontrol, Ecorestoration Fourth Edition, Springer, 2017).

7. Yu, L., Nicolaisen, J., Larsen, J. & Ravnskov, S. Molecular characterization of root-associated fungal communities in relation to health status of *Pisum sativum* using barcoded pyrosequencing. *Plant Soil* **357**, 395–405 (2012).

8. Corredor, A. H., Van Rees, K. & Vujanovic, V. Host genotype and health status influence on the composition of the arbuscularmycorrhizal fungi in Salix bioenergy plantations. *For. Ecol. Manag.* **314**, 112–119 (2014).

9. Martinez, N. & Johnson, N. C. Agricultural management influences propagule densities and functioning of arbuscularmycorrhizas in low- and high-input agroecosystems in arid environments. *Appl. Soil Ecol.* **46**, 300–306 (2010).

10. Hontoria, C., García-González, I., Quemada, M., Roldánd, A. & Alguacil, M. M. The cover crop determines the AMF community composition in soil and in roots of maize after a ten-year continuous crop rotation. *Sci. Total Environ.* **660**, 913–922 (2019).

11. Lumini, E., Vallino, M., Alguacil, M. M., Romani, M. & Bianciotto, V. Different farming and water regimes in Italian rice fields affect arbuscularmycorrhizal fungal soil communities. *Ecol. Appl.* **21**, 1696–1707 (2011).

12. Manoharan, L., Rosenstock, N. P., Williams, A. & Hedlund, K. Agricultural management practices influence AMF diversity and community composition with cascading effects on plant productivity. *Appl. Soil Ecol.* **115**, 53–59 (2017).

13. Dai, M., Bainard, L. D., Hamel, C., Gan, Y. & Lynch, D. Impact of land use on arbuscularmycorrhizal fungal communities in rural Canada. *Appl. Environ. Microbiol.* **79**, 6719–6729 (2013).

14. Aghili, F. et al. Wheat plants invest more in mycorrhizae and receive more benefits from them under adverse than favorable soil conditions. *Appl. Soil Ecol.* **84**, 93–111 (2014).

15. Gaudin, J., Semetey, O., Foissac, X. & Eveillard, S. Phytoplasmatiter in diseased lavender is not correlated to lavender tolerance to stolburphytoplasma. *Bull. Insectol.* **64** (Supplement), S179–S180 (2011).

16. Kamińska, M., Klamkowski, K., Berniak, H. & Treder, W. Effect of arbuscularmycorrhizal fungi inoculation on aster yellows phytoplasma-infected tobacco plants. *Sci. Hortic.* **125**, 500–503 (2010).

17. Battle, A. et al. Tolerance increase to *Candidatus phytoplasma prunorum* in mycorrhizal plums fruit trees. *Bull. Insectol.* **64**, 125–126 (2011).

18. D'ameli, R. et al. Increased plant tolerance against chrysanthemum yellows phytoplasma (*Candidatus Phytoplasma asteris*) following double inoculation with Glomusmosseae BEG12 and *Pseudomonas putida* SIPf1Rif. *Plant. Pathol.* **60**, 1014–1022 (2011).

19. Fiorilli, V. et al. Omics approaches revealed how arbuscularmycorrhizal symbiosis enhances yield and resistance to leaf pathogen in wheat. *Sci. Rep.* **8**, 9625 (2018).

20. Bødker, L., Kjøller, R., Kristensen, K. & Rosendahl, S. Interactions between indigenous arbuscularmycorrhizal fungi and *Aphano-mycyes euteiches* in field-grown pea. *Mycorrhiza* **12**, 7–12 (2002).
21. Al-Askar, A. A. & Rashad, Y. M. Arbuscularmycorrhizal fungi: a biocontrol agent against common bean Fusarium root rot disease. *Plant Pathol. J.* **9**, 31–38 (2010).

22. Hugoni, M., Luis, P., Guyonnet, J. & Haichar, F. Z. Plant host habitat and root exudates shape fungal diversity. *Mycorrhiza* **28**, 451–463 (2018).

23. Bertaccini, A. & Duduk, B. Phytoplasma and phytoplasma diseases: A review of recent research. *Phytopathol. Mediterr.* **48**, 355–378 (2009).

24. Stierlin, E., Nicolè, F., Costes, T., Fernandez, X. & Michel, T. Metabolomic study of volatile compounds emitted by lavender grown under open-field conditions: a potential approach to investigate the yellow decline disease. *Metabolomics* **16**, 31 (2020).

25. Lopez-Garcia, A. et al. Plant traits determine the phylogenetic structure of arbuscularmycorrhizal fungal communities. *Mol. Ecol.* **26**, 6948–6959 (2017).

26. Alguacil, M. M., Díaz, G., Torres, P., Rodríguez-Caballero, G. & Roldan, A. Host identity and functional traits determine the community composition of the arbuscularmycorrhizal fungi in facultative epiphytic plant species. *Fungal Ecol.* **39**, 307–315 (2019).

27. Neuenkamp, L. et al. The role of plant mycorrhizal type and status in modulating the relationship between plant and arbuscularmycorrhizal fungal communities. *New Phytol.* **220**, 1236–1247 (2018).

28. Alguacil, M. M., Torrecillas, E., García-Orenes, F. C. & Roldán, A. Changes in the composition and diversity of AMF communities mediated by management practices in a Mediterranean soil are related with increases in soil biological activity. *Soil Biol. Biochem.* **76**, 34–44 (2014).

29. Giri, B. & Mukerji, K. G. Mycorrhizal inoculant alleviates salt stress in *Sesbania aegyptiaca* and *Sesbania grandiflora* under field conditions: evidence for reduced sodium and improved magnesium uptake. *Mycorrhiza* **14**, 307–312 (2004).

30. Phillips, J. M. & Hayman, D. S. Improved procedure for clearing roots and staining parasitic and vesicular-arbuscularmycorrhizal fungi for rapid assessment of infection. *Trans. Br. Mycol. Soc.* **55**, 158–163 (1970).

31. Trouvelot, A., Kough, J. L. & Gianinazzi-Pearson, V. Mesure du taux de mycorhization VA d'un système radiculaire. Recherche de méthodes ayant une signification fonctionnelle. In *Physiological and genetical aspects of mycorrhizae* (eds Gianinazzi-Pearson, V. & Gianinazzi, S.) 217–221 (INRA Press, Paris, 1986).

32. Gollotte, A., van Tuinen, D. & Atkinson, D. Diversity of arbuscularmycorrhizalfungicolonisingroots of the grassspecies *Agrostis capillaris* and *Lolium perenne* in a fieldexperiment. *Mycorrhiza* **14**, 111–117 (2004).

33. Binet, M. N. et al. Responses of above- and below-ground fungal symbionts to cessation of mowing in subalpine grassland. *Fungal Ecol.* **25**, 14–21 (2017).

34. Mouhamadou, B. et al. Effects of two grass species on the composition of soil fungal communities. *Biol. Fertil. Soils* **49**, 1131–1139 (2013).

35. Boyer, F. et al. Obitools: a unix- inspired software package for DNA metabarcoding. *Mol. Ecol. Resour.* **16**, 176–182 (2016).

36. Lentendu, G. et al. Assessment of soil fungal diversity in different alpine tundra habitats by means of pyrosequencing. *Fungal Div.* **49**, 113–123 (2011).

37. van Dongen, S. Graph clustering by flow simulation. Ph.D. thesis, University of Utrecht (2000).

38. Thompson, L. A. S-PLUS (and R) manual to accompany Agresti's Categorical Data Analysis (2002), 2nd ed (2009).

39. Oksanen, J., Kindt, R., Legendre, P., O'Hara, B. & Gavin, L. vegan: Community Ecology Package. *R package* version 1.15–4 (2009).

40. Mouhamadou, B. et al. Molecular screening of xerophilic *Aspergillus* strains producing mycophenolic acid. *Fungal Biol.* **121**, 103–111 (2017).

## Acknowledgments

This study was conducted as part of FEADER project Mycolav and was supported by European Union and Région Sud-Provence-Alpes-Côte d'Azur. We thank Eric Chaisse, Marie Fontaine and Guillaume Vial for help and assistance in the fields. We address our special thanks to Bruno, Charlotte, Louis, Nael, Ilian and Lila.

## Author contributions

M.-N.B. and B.M. were responsible of this project. They wrote the whole manuscript and reviewed it. C.M. fulfilled the bioinformatic and statistical analyses. R.G. dealt with the bioinformatic analyses. O.B., B.C., D.F., B.D., J.L. and S.P participated in the experiments. V.B. participated in the preparation of the manuscript and in the reviewing of the article.

## Competing interests

The authors declare no competing interests.

## Additional information

Supplementary information is available for this paper at https://doi.org/10.1038/s41598-020-77240-6.

Correspondence and requests for materials should be addressed to B.M.

Reprints and permissions information is available at www.nature.com/reprints.

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

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

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons licence, and indicate if changes were made. The images or other third party material in this article are included in the article's Creative Commons licence, unless indicated otherwise in a credit line to the material. If material is not included in the article's Creative Commons licence and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this licence, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2020