# Glomerales Dominate Arbuscular Mycorrhizal Fungal Communities Associated with Spontaneous Plants in Phosphate-Rich Soils of Former Rock Phosphate Mining Sites

Amandine Ducousso-Détrez, Robin Raveau, Joel Fontaine, Mohamed Hijri, Anissa Lounès-Hadj Sahraoui

► To cite this version:

Amandine Ducousso-Détrez, Robin Raveau, Joel Fontaine, Mohamed Hijri, Anissa Lounès-Hadj Sahraoui. Glomerales Dominate Arbuscular Mycorrhizal Fungal Communities Associated with Spontaneous Plants in Phosphate-Rich Soils of Former Rock Phosphate Mining Sites. Microorganisms, 2022, 10 (12), pp.2406. 10.3390/microorganisms10122406. [https://hal.inrae.fr/hal-03932239](https://hal.inrae.fr/hal-03932239)

HAL Id: hal-03932239

https://hal.inrae.fr/hal-03932239

Submitted on 10 Jan 2023

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
<figure><img src="image_2.png" /></figure>

Article

# Glomerales Dominate Arbuscular Mycorrhizal Fungal Communities Associated with Spontaneous Plants in Phosphate-Rich Soils of Former Rock Phosphate Mining Sites

Amandine Ducousso-Détrez <fnref n="1" />,<fnref n="2" />, Robin Raveau <fnref n="3" /> ID, Joël Fontaine <fnref n="1" /> ID, Mohamed Hijri <fnref n="2" />,<fnref n="4" /> ID
and Anissa Lounès-Hadj Sahraoui <fnref n="1" />,* ID

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

check for
updates

**Citation:** Ducousso-Détrez, A.; Raveau, R.; Fontaine, J.; Hijri, M.; Lounès-Hadj Sahraoui, A. Glomerales Dominate Arbuscular Mycorrhizal Fungal Communities Associated with Spontaneous Plants in Phosphate-Rich Soils of Former Rock Phosphate Mining Sites. *Microorganisms* **2022**, *10*, 2406. https://doi.org/10.3390/microorganisms10122406

Academic Editors: Annamaria Bevivino and Maria Maddalena Del Gallo

Received: 2 October 2022
Accepted: 30 November 2022
Published: 5 December 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

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

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

**Abstract:** Arbuscular mycorrhizal fungi (AMF) are key drivers of soil functioning. They interact with multiple soil parameters, notably, phosphorus (P). In this work, AMF communities of native plants grown spontaneously on former mining sites either enriched (P sites) or not enriched with P (nP sites) by mining cuttings of rock phosphate (RP) were studied. No significant differences were observed in the root mycorrhizal rates of the plants when comparing P and nP sites. The assessment of AMF diversity and community structure using Illumina MiSeq metabarcoding and targeting 18S rDNA in roots and rhizospheric soils showed a total of 318 Amplicon Sequence Variants (ASVs) of Glomeromycota phylum. No significant difference in the diversity was found between P and nP sites. Glomeraceae species were largely dominant, formed a fungal core of 26 ASVs, and were persistent and abundant in all sites. In the P soils, eight ASVs were identified by indicator species analysis. A trend towards an increase in Diversisporaceae and Claroideoglomeraceae and a reduction in Paraglomeraceae and Glomeraceae were noticed. These results provide new insights into AMF ecology in former RP mining sites; they document that P concentration is a driver of AMF community structures in soils enriched in RP long term but also suggest an influence of land disturbance, ecosystem self-restoration, and AMF life history strategies as drivers of AMF community profiles.

**Keywords:** arbuscular mycorrhizal fungi; phosphorus; rock phosphate; diversity; taxonomic composition

## 1. Introduction

Many efforts have been deployed to sustain agricultural production systems while reducing their environmental footprint. Ecosystemic contributions of microbial resources to improve and maintain soil fertility are among the strategies that have been utilized. In particular, arbuscular mycorrhizal fungi (AMF) have gained a large interest as biological input in agroecosystems [1–3]. Indeed, AMF are ubiquitous fungi that form symbiotic associations with most plant species, including important crop plants, such as vegetables and cereals [4,5]. They are obligate biotrophs, receiving their carbon sources in the form of photosynthetic sugars and fatty acids from plant hosts [6,7]. In exchange, elements such as phosphorus (P), nitrogen, sulfur, or microelements, such as copper and zinc, may also be transferred to the host via arbuscules, favoring plant nutrient supply [1,8,9]. Mycorrhizal symbiosis can also contribute to plant protection against biotic or abiotic stresses [10–13]. In addition, evidence also indicate that AMF can be critical components to improve soil structure, water retention capacity, and soil fertility and reduce soil nutrient

<footnote n="1">¹ Unité de Chimie Environnementale et Interactions sur le Vivant (UCEIV), Université du Littoral Côte d'Opale, UR 4492, SFR Condorcet FR CNRS 3417, CEDEX, 62228 Calais, France</footnote>

<footnote n="2">² Institut de Recherche en Biologie Végétale (IRBV), 3 AgroBioSciences, Université de Montréal, Montréal, QC H1X 2B2, Canada</footnote>

<footnote n="3">³ INRAE, UMR SAVE, Bordeaux Science Agro, ISVV, 33882 Villenave d'Ornon, France</footnote>

<footnote n="4">⁴ African Genome Center, Mohammed VI Polytechnic University (UM6P), Ben Guerir 43150, Morocco</footnote>

<footnote>* Correspondence: anissa.lounes@univ-littoral.fr</footnote>
leaching [1,14–17]. Concomitantly, the multiple benefits that AMF confer to their hosts and agroecosystems’ functioning have raised opportunities for their application as commercial biofertilizers. In particular, AMF increase the uptake of relatively immobile ions (such as phosphate ions) for their host plant due to an extensive extraradical hyphal network that extends beyond the zone of direct uptake by roots. In this way, they contribute to greater soil prospection and enable mycorrhizal plants to access more soil P resources. Therefore, AMF have great importance to productivity and plant growth in most ecosystems [1,18]. However, despite the great potential of AMF and the evidence of their relevance for agriculture [19], a deeper understanding of AMF ecology is of primary importance to improving engineering and the use of AMF-based bioinoculants in ecosystems.

In such context, significant improvements have been made in the detailed analysis of AMF communities among complex habitats using AMF gene markers and high-throughput sequencing technology for taxonomic identification, phylogenetic reconstruction, or quantification of AMF [20–23]. Thus, due to such improvements, much compelling evidence has reported that different abiotic and biotic factors interfere with AMF communities at either a global or local scale [24–27]. In particular, the host plant can be considered a strong factor in shaping AMF communities [28–30], but these are also known to be sensitive to interactions with other soil organisms [31,32]. Edaphic properties [24,29,33] and climatic parameters [10,25,34,35] are also recognized as strong drivers of AMF communities. In addition, some studies have suggested a high diversity of AMF communities in natural habitats [36,37] while a lower diversity tends to be observed in agricultural ecosystems [38]. Notably, AMF communities’ characteristics may vary with land use (usage intensity, cultivation duration, pasture, etc.) and with agronomic practices for soil management (chemical inputs, plowing, tillage, crop rotation, land conversion, etc.) [39,40]. In particular, the responses of AMF to soil P were largely investigated [41–43]. So, research efforts on the understanding of P as a key driver of AMF colonization, diversity, and community composition have gained interest, in particular, for identifying guides in the development of AMF inoculants and their use with a reduction in the use of chemical fertilizers [44,45]. In particular, how AMF communities interact with rock phosphate (RP) has attracted considerable attention because RP is a natural, low industrial-cost phosphate input that has more ecological and sustainable acceptability than chemical fertilizers. However, P as a key player in shaping AMF communities still remains a controversial subject in the literature. Thus, the impact of P concentrations on AMF colonization as well as on the diversity and composition of AMF communities is regularly discussed [28,39].

Thus, the current work aims at characterizing the diversity and composition of AMF communities across a former mining area of RP where two types of sites occurred: sites enriched or not enriched with phosphate deposits originating from the mining activity. The AMF community profiles in the soil and root samples were determined using MiSeq amplicon sequencing targeting 18S rDNA that generates Amplicon Sequence Variants (ASVs). The diversity and composition of AMF were analyzed at several levels of taxonomic resolution, and the potential role of the P gradient and soil restoration process in AMF community composition was discussed.

## 2. Materials and Methods

### 2.1. *Study Area and Experimental Design*

This work was carried out in the region of the “Phosphatières du Quercy” (southern France; N 44.351827, E 1.691021). This area is characterized by numerous paleokarst whose fillings are enriched in phosphorite (i.e., a phosphatic ore containing variable proportions of tricalcium phosphate). There, in the past, mining extraction was carried out in the form of open-pit mines (Figure S1a), and some excavated spoils were abandoned near the extraction points. Thus, each selected site corresponds to a former mining point (i.e., a site without enrichment with RP ores, referred hereafter as “nP site”) associated with its spoil area (i.e., a site with ore deposits due to mining exploitation process, referred to as “P site”). Today,
in this former mining area, the dominant landscapes are dry grasslands that are punctually wooded (Figure S1).

The P concentrations are known to be more or less high in each karst phosphatic fill [46]. Three nP were therefore chosen, expecting significantly high P concentrations. The choice of locations was also constrained (i) by the possibility of locating potential P and nP sites in close proximity to each other from the ground surface and (ii) by the necessary presence of the same set of species constituting the plant sampling. Thus, four native herbaceous mycotrophic species were sampled: *Ranunculus bulbosus*, *Bromus sterilis*, *Taraxacum officinale*, and *Dactylis glomerata*, and three plants per species were collected. From these 3 plants, the soil closely attached to roots (i.e., rhizospheric soil) was manually removed by gentle agitation and pooled into a single soil fraction. Then, the roots were washed with sterile water to get rid of remaining soil particles and then pooled into one root fraction per species (Figure 1). Soil and roots were then, respectively, kept at 4 °C and −20 °C until further DNA extraction. This was repeated for each plant species in each site, leading to 24 root samples (6 sites × 4 plant species) and 24 rhizospheric soil samples. Detailed soil properties for each site are available in supplementary data (Table S1).

<figure><img src="image_5.png" /><figcaption>**Figure 1.** Sampling sites and experimental design.</figcaption></figure>

## 2.2. Estimation of Arbuscular Mycorrhizal Root Colonization

Immediately after plant harvesting, from each plant species and site, fractions of fresh roots were stained using nonvital Trypan blue (0.5 g·L⁻¹) according to the method described by [47]. From each stained root batches (24 batches), fragments of 1 cm were sampled, and 3 intersections per root fragment were examined through microscopic observations
(Nikon Eclipse E600, ×100 magnification). The rate of root colonization by AMF structures (vesicles and arbuscules) was quantified as described by [48].

## 2.3. Soil Physicochemical Properties

In each site, the soil physicochemical properties (pH, soil texture, chemical composition) were measured by the CIRAD-US Analyse laboratory (Montpellier, France) using inductively coupled plasma spectrometry, atomic emission spectrometry, and X-ray fluorometry. The 6 soils were characterized by near-neutral pH, ranging between 6.9 and 7.3. Regarding granulometry, all topsoil (0–30 cm) resulted in high levels of clay with additionally higher percentages of coarse sands in P soils. In nP soils, total P concentrations ranged from 1057.1 to 1496.3 mg·kg⁻¹ while in P soils, they ranged from 2880.0 to 13,927.9 mg·kg⁻¹. Likewise, available Olsen P concentrations in P soils were between 46.1 mg·kg⁻¹ to 339.5 mg·kg⁻¹ while nP soils were characterized by values ranging from 5.04 to 12.82 mg·kg⁻¹. Thus, by comparing matched soils from the same localization, P and Olsen P ratios ranging from 2.07 to 13.05 and from 9.14 to 41.96 Pi, respectively, were obtained (Table S1), which allowed to qualify the P sites as high-P sites, and the nP sites as low-P sites.

## 2.4. DNA Extraction, PCR, and Sequencing

### 2.4.1. DNA Extraction

Genomic DNA extraction from the rhizospheric soil samples was directly performed from 250 mg of soil using a Nucleospin Soil® kit (Macherey-Nagel, Düren, Germany) according to the manufacturer's instructions.

Genomic DNA extraction from root samples was operated from roots initially frozen in liquid nitrogen and then reduced to fine powder. Then, the extraction from 250 mg of ground roots was performed using Cetyltrimethylammonium bromide (CTAB, 1.4 M NaCl₂ 100 mM Tris-HCl, pH 8.0; 20 mM EDTA, pH 8.0; 2% CTAB), Polyvinylpyrrolidone (PVP 1% *w*/*v*), β-Mercaptoethanol (5% *v*/*v*), and activated charcoal (0.5% *w*/*v*) extraction (30 min; 55 °C). A centrifugation step was then carried out (10 min; 16,000 × *g*), and lysate extraction was performed with chloroform: isoamyl alcohol (24:1). DNA precipitation was obtained in isopropanol (1 h incubation; 25 °C) followed by centrifugation (10 min; 700 × *g*). After three washings (with ice-cold ethanol (70%) followed by centrifugation—10 min at 900 × *g*) and air-drying at room temperature (approx. 90 min; 20 °C), the DNA pellet was dissolved in 50 µL of TE buffer (10 mM Tris-HCl, pH 8.0; 1.0 mM EDTA, pH 8.0) [13,49,50].

All extractions were performed in triplicate. DNA quality of root and soil DNA extracts was assessed using 1% (*w*/*v*) agarose gel, and measures of the 260/280 nm and 260/230 nm ratios were performed with a SpectraMax® iD3 device (Molecular Devices LLC, Sunnyvale, CA, USA). After measures of DNA concentrations, root and soil DNA extracts were normalized to 25 ng·L⁻¹ for further analyses. The extracted DNA was stored at −20 °C until further use.

### 2.4.2. PCR Targeting the 18S rRNA Gene of AMF

The DNA of AMF was specifically targeted using Nested-PCR (Surecycler 8800, Agilent Technologies, Les Ulis, France) with a set of two primer pairs as previously published by [50]. Thus, the first-round PCR was performed using the AMF-discriminating primer pair AML1 (3'-ATCAACTTTCGATGGTAGGATAGA-5') and AML2 (3'-GAACCCAAACAC TTTGGTTTCC-5') which generate amplicons of the small 18S subunit of the rRNA gene of about 800 bp in length [51,52]. The PCR conditions were as follows: initial denaturation at 94 °C for 3 min followed by 35 cycles at 94 °C for 1 min, 45 °C for 1 min, and 72 °C for 1 min, and a final elongation step at 72 °C for 5 min [13]. PCR reactions were performed in a reaction volume of 25 µL, and reagents were as follows: 5 µL of Q5® High-Fidelity DNA Polymerase (New England Biolabs France, Evry, France), 0.8 µL of each primer (0.4 µM), 1 µL of dNTPs (0.2 mM), 1 µL of DMSO, 1 µL of BSA (100 µg·mL⁻¹), and 1 ng of DNA template [13,51].
From the 800 bp length amplicons obtained after the first-round PCR, which are not compatible with Illumina sequencing technology, nested PCRs were performed using an in-house set of internal primers nu-SSU-0595-5-F (ACACTGACGACATGGTTCTACA CG-GTAATTCCAGCTCCAATAG) and nu-SSU-0948-3-R (TACGGTAGCAGAGACTTGGTCT TTGATTAATGAAAACATCCTTGGC). This set of primers flanked the V4 region of 18S and hence downsized the length of the amplicons to about 400 bp, which is compatible with MiSeq 300 paired-end sequencing [53]. The two primers were complemented with CS1 and CS2 barcoded adapters [13,53]. The reaction mixture was the same as previously described, but 1 ng of DNA from the first-round PCR as template was used. The reaction conditions for the second PCR were as follows: initial denaturation at 94 °C for 3 min followed by thirty-five cycles at 94 °C (1 min), 58 °C (1 min), and 72 °C (1 min), and a final elongation step at 72 °C (5 min).

### 2.4.3. Sequencing

From genomic DNA of root and soil samples collected for each plant species harvested in each site, PCR reactions were performed in triplicates for each sample, and these PCR products were then pooled together. Sequencing was processed by Genome Quebec Innovation Centre (Montreal, QC, Canada) by using an Illumina MiSeq platform generating 2 × 300 bp paired-end reads.

### 2.4.4. Nucleotide Sequence Accession Number

The 18S rRNA gene sequences of the raw data set have been deposited in the NCBI Sequence Read Archive database under the project accession number PRJNA877825.

## 2.5. *Bioinformatic and Statistical Analyses*

Bioinformatic process and statistical analyses were operated in the R 4.0.2 software (R Core Team, 2019) environment. The DADA2 pipeline (v. 1.16) [54], an open-source program implemented in R package (https://benjjneb.github.io/dada2/tutorial.html accessed on 1 October 2022), was used to process the sequencing data. Briefly, sequence reads were filtered and trimmed using optimized parameter settings as recommended. Sequence reads were dereplicated, denoised, and merged using DADA2 default parameters. Then, sequences were aligned and categorized to infer amplicon sequence variants (ASVs) grouping amplicon sequencing data by using 100% of sequence identity [55]. Sequences presenting only once (singletons) in the whole data set were eliminated. Validity of sequencing depth was controlled from rarefaction curves, computed using the “rarecurve” function from the Vegan package in R.

The taxonomic assignment of ASVs was performed following a previously established two-step approach [13,53]. In the first step, the Silva v132 database formatted for DADA2 was used to assign ASVs from kingdom to genus (minimum bootstrap 80) with the assignTaxonomy () command and the blast parameters as constitutively formatted for DADA2 [56]. From the ASVs assigned as AMF with Silva database, affiliation of each ASV was assessed at the genus level by performing a BLAST analysis against the NCBI database and MARJAAM, a web-based database containing referenced Glomeromycota DNA sequence data. The ASVs identified as non-Glomeromycota at the phylum level at the end of these two stages of assignment were excluded from further analyses.

In a second step, aiming at a refined taxonomic identification of each ASV, a phylogenetic tree was subsequently constructed as proposed by Stefani et al., 2020 [53]: multiple alignment including ASVs identified as Glomeromycota and multiple consensus sequences was computed using the web portal Kalign [57] (http://msa.sbc.su.se/ accessed on 1 October 2022). Then, a maximum-likelihood tree was obtained using RAxML v8.2.10 [58], through the CIPRES web portal [59]. Visualization of the output was finally obtained with the FigTree v1.4.4 program.

With the aim to identify the influence of former RP inputs on AMF communities, we studied AMF communities across the six sampling sites as well as inside root or
rhizospheric soil using proxy richness, alpha diversity, distribution, and abundance of taxa. Then, data relative to P sites were compared to those of nP sites, inside each location, or after pooling data from the 3 P sites on the one hand and the 3 nP sites on the other hand. The Chao1 richness estimator and alpha diversity indices (Shannon and Simpson) [60–62] were computed from the plot_richness() function using the phyloseq R package.

The normality and homoscedasticity of the data were assessed from the Shapiro-Wilk and Bartlett tests, respectively. If both conditions were verified, ANOVA analysis complemented with a posthoc Dunn test was carried out. Otherwise, a Kruskal-Wallis nonparametric test ("kruskal.test" function in R) was used. The significance of the statistical analyses was considered for $\alpha = 0.05$.

Venn diagrams were constructed with the Vegan package to visualize the numbers of ASVs that are either specific to one site or shared between paired sites inside the mining area; significance of the difference in ASVs numbers was then conducted using chi2 test. Indicator species analysis was performed using the multipatt function from indicspecies package [63] to compute the indicator value index (IndVal), i.e., a measure of specificity (based on abundance values) and fidelity (computed from presence data) of each ASV to a targeted clustering group of samples related to a targeted ecological condition [60]. This index, therefore, led to identifying ASVs (hereafter referred to as "indicator ASVs") that can be considered closely related to the ecological condition of their group [64]. The significance of IndVal of each species was assessed by a random permutation procedure from 9999 permutations with a significance at the $\alpha = 0.05$ level (IndVal's p_values with the padjust() function).

## 3. Results

### 3.1. Mycorrhizal Colonization of Plant Roots

The AMF root colonization was estimated in each sampling site from 135 microscopic observations of different individuals per plant species. Root colonization rates ranged from 13 to 63% across the different sampling sites (L1-P vs. -nP, L2-P vs. -nP, and L3-P vs. -nP) and plant species (*Ranunculus bulbosus*, *Bromus sterilis*, *Taraxacum officinale*, or *Dactylis glomerata*). However, when considering all individuals of the four plant species as a pool (to reduce a putative effect only linked to plant species), no significant difference in mycorrhizal colonization between plants growing in P versus nP sites was observed, whether the estimate was made for the entire mining area or within each individual localization. Thus, considering the entire mining area, the mean percentages reached 42 and 43% in plant roots sampled in P versus nP sites, respectively (Table 1).

<table><caption>**Table 1.** Mycorrhizal rates of the four plant species (*Ranunculus bulbosus*, *Bromus sterilis*, *Taraxacum officinale*, and *Dactylis glomerata*) grown on the six sites (L1-P and -nP, L2-P and -nP, and L3-P and -nP). Data based on 135 root fragments from three plants per species in each sampling site.</caption><thead><tr><th rowspan="2">Sites</th><th rowspan="2">Total P<br>Concentration<br>(mg kg<sup>-1</sup>)</th><th colspan="5">Mycorrhizal Rates (%)</th><th rowspan="2">Average P<br>versus nP</th></tr><tr><th><i>Ranunculus<br>bulbosus</i></th><th><i>Taraxacum<br>officinale</i></th><th><i>Dactylis<br>glomerata</i></th><th><i>Bromus<br>sterilis</i></th><th>Average<br>Per Site</th></tr></thead><tbody><tr><td>L1-P</td><td>2860</td><td>38</td><td>40</td><td>33</td><td>51</td><td>41</td><td></td></tr><tr><td>L2-P</td><td>13,928</td><td>20</td><td>27</td><td>42</td><td>18</td><td>27</td><td>42</td></tr><tr><td>L3-P</td><td>10,739</td><td>53</td><td>53</td><td>62</td><td>63</td><td>58</td><td></td></tr><tr><td>L1-nP</td><td>1380</td><td>53</td><td>58</td><td>40</td><td>56</td><td>52</td><td></td></tr><tr><td>L2-nP</td><td>1067</td><td>47</td><td>13</td><td>5</td><td>16</td><td>32</td><td>43</td></tr><tr><td>L3-nP</td><td>1496</td><td>58</td><td>44</td><td>33</td><td>42</td><td>44</td><td></td></tr></tbody></table>

### 3.2. Analysis of AMF Diversity

After Illumina Miseq sequencing, the global dataset resulted in a total of 4,795,981 raw MiSeq reads (yielded across 24 soil samples and 24 root samples, Table S2). After quality
filtering and chimera removal, a total of 3,527,278 sequences were retained as the AMF 18S rDNA dataset (Table S3). The rarefaction curves obtained after quality filtering reached a plateau, indicating we had a good representation of the microbial community, as most of the abundant species were represented and the sequencing depth effort was adequate to progress to further AMF community analysis (Figure S2).

From the AMF 18S rDNA dataset and using the Silva v132 database formatted for DADA2, 392 ASVs were inferred and identified as Glomeromycota. Among them, 74 ASVs were not assigned at the phylum level after inferring the phylogenetic tree (Figure S3). These ASVs were therefore excluded from further analysis. Thus, a total of 318 ASVs were kept for further comparisons of the AMF communities within P and nP soils.

### 3.3. Richness and Alpha-Diversity

The diversity (Shannon and Simpson indexes) and richness (Chao1 estimator) were used to compare AMF communities across P versus nP sites as well as between root or rhizosphere biotopes. A significant difference was observed between the two habitats, with higher diversity in roots compared to rhizosphere soil biotopes. In contrast, no significant difference was recorded comparing the P and nP sites either for the entire mining area (Figure 2) or within each localization (data not shown). No impact of P versus nP status was noticed in either the root or rhizosphere biotopes (data not shown).

<figure><img src="image_6.png" /><figcaption>**Figure 2.** Richness (Chao 1 estimator) and diversity (Shannon and Simpson) indexes comparing AMF communities of (A) P versus nP soils, root and soil data being pooled; (B) root versus rhizospheric soil samples, P and nP site data being pooled. Within each graph, different letters denote significant differences using the Kruskal-Wallis nonparametric test (*p* < 0.05).</figcaption></figure>

### 3.4. Taxonomic Assignment of ASVs

The 318 ASVs identified as Glomeromycota were assigned to Glomerales (81%), Diversisporales (12%), Paraglomerales (3%), and Archaeosporales (2%) (Figure 3A) while no taxonomic assignment at order level was obtained for 2% of them. Among them, five families were identified: Glomeraceae, Claroideoglomeraceae, Diversisporaceae, Paraglomeraceae, and Archeosporaceae. About 20.7% of the ASVs could be affiliated to a genus, mainly identified as *Rhizophagus*, *Funneliformis*, and *Glomus* (18, 15, and 14 ASVs,
respectively), but also, to a lesser extent, as *Archaeospora*, *Paraglomus*, *Septoglomus*, and *Claroideoglomus* (Table S4). With assignment up to the species rank, 39 ASVs (12.2%) could be named in line with referenced strains in the published databases. They were representative of eight species and six genera, namely *Archaeospora trappei*, *Claroideoglomus lamellosum*, *Funneliformis mosseae*, *Glomus indicum*, *G. iranicum*, *Rhizophagus irregularis*, *R. vesiculiferus*, and *Septoglomus africanum* (Figure 3B). Among them, the prevalent ones in terms of both read count and ASV numbers were *F. mosseae* (10 ASVs totaling together 12,591 reads, i.e., 2.7% of the readings), *G. indicum* (9 ASVs with 2698 reads, i.e., 0.6%), and *R. vesiculiferus* (6 ASVs with 8033 reads, i.e., 1.75%). Among ASVs assigned up to the species level, ASV 36, referred to as *F. mosseae*, was dominant in terms of relative abundance (2.6%) (Table S5).

<figure><img src="image_7.png" /><figcaption>**Figure 3.** Taxonomic distribution of the 318 ASVs identified as Glomeromycota inside the entire mining area (A) at the order level, (B) at the species level, and (C) taxonomic distribution of the 318 ASVs at the family level inside each sample. “NA” category: ASV that did not obtain taxonomic assignment at the studied level.</figcaption></figure>
Regarding variations across the samples (Figure 3C), Glomeraceae were largely dominant in all samples (i.e., across root and soil samples for each site) with Diversisporaceae and Paraglomeraceae mainly observed in soil.

### 3.5. ASVs Shared by the Six Mining Sites

When having a look at the different experimental conditions, a core of 26 ASVs (about 8.2% of the total ASV number) shared by the six sampling sites was identified (Table S6). These ASVs, ubiquitously present, correspond to about 64% of the total read number retained after the bioinformatic processing, taxonomic, and phylogenetic assignment. Among them, three were dominant (ASV61, ASV62, and ASV208), contributing to 57.1% of the read number of the 26 ASVs. With the exception of two ASVs assigned to Diversisporaceae, all core ASVs belonged to Glomeraceae. Two of them were identified up to the species level with a taxonomic identity inferred to *R. vesiculiferus*.

### 3.6. Dissimilarity in the AMF Community Composition and Identification of Indicator Species Comparing P and nP Sites

By comparing P versus nP sites, ASVs that were specific (i.e., occurring exclusively across the P versus nP sites) and ASVs that were shared (i.e., observed at both P and nP sites) were observed, as shown in Figure 4. Thus, considering the complete mining area, we identified 100 ASVs that are shared by the two profiles, representing 31.4% of the total ASV number but 86.5% of the total read number. In contrast, 138 ASVs (43.4%), accumulating 8.6% of reads, were identified as nP-site-specific (hereafter referred to as “nP-specific ASVs”). Concomitantly, 80 ASVs (25.2%) were P specific, displaying 4.85% of the reads (Figure 4).

Considering each localization, the same trend was observed, with P sites showing a lower number of specific ASVs compared to non-P sites. The most contrasting difference was observed in location two, where the number of nP-specific ASVs was almost twice the amount of P-specific ASVs (52 ASVs, or 41.6%, versus 28 ASVs, or 22.4%, respectively) (Figure 4). Variations in ASVs abundance among P vs. nP sites were also observed, as depicted in Figure 5, focusing on the 50 most abundant ASVs. In addition, it highlights the higher abundance of ASV 61, 62, and 208 in all sites.

Moreover, indicator species analysis was performed to search ASVs indicative of the following sample groups: the P versus nP group, the root versus soil group, as an individual group, or in combination (Table S7). Thus, in particular, seven nP-indicator ASVs and one P-indicator ASV were identified. All of them were assigned to Glomeraceae.

<figure><img src="image_8.png" /><figcaption>**Figure 4.** *Cont.*</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>**Figure 4.** Venn diagrams showing the overlap of the AMF communities across P and nP sites, considering the entire mining area (A) or each localization individually (L1 localization, L2, and L3). Shared and specific ASV numbers and their relative abundance (in percent) are given inside the circles; the read numbers are given out of the circles.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 5.** Heatmap representation showing the abundance of the 50 most abundant ASVs across all the sampling sites. The 50 most abundant ASVs include the 26 core ASVs shared by the six sampling sites *: Core ASVs shared by all sampling sites.</figcaption></figure>
Regarding the taxonomic classification, no Archaeosporaceae member was observed in P soils. Paraglomeraceae and Glomeraceae tend to be reduced in P soils in terms of ASV percentages (Figure 6). In contrast, Diversisporaceae and Claroideoglomeraceae families tended to have increased with the P content. Lower percentages of assigned ASV were recorded at genus and species levels in P soils (Figure 6).

<figure><img src="image_11.png" /><figcaption>**Figure 6.** Distribution at the family, genus, and species level of the shared and specific ASVs identified in P and nP sites across the total mining area, NA: no-assigned ASVs.</figcaption></figure>

In addition, comparisons between the three P sites were made. Different ASVs between the three P sites were identified with a higher number of site-specific ASVs and few shared ASVs, as shown in the Venn diagrams (Figure S4A). The same trend occurred comparing the nP sites (Figure S4B).

**4. Discussion**

In this study, we investigated AMF root colonization and taxonomic profiles of AMF communities in a former mining area comparing P and nP sites. These sites were characterized by total P concentrations ranging from 1067.1 to 13,927.9 mg·kg⁻¹ while the values conventionally reported in the literature generally range from 40 to 3000 mg·kg⁻¹ [65–67]. Additionally, the soluble P (P Olsen) concentrations ranged from 5 to 339.5 mg·kg⁻¹, which is largely above those displayed, for instance, in soils supplemented with mineral or organic fertilizers [67]. Such high concentrations in soluble P content are scarce and unexpected. Indeed, in the phosphatieres, soils are enriched with apatite, a form of highly insoluble soil P, and aluminous phosphate parageneses, clays, and sands [46]. As the soluble orthophosphate ions act as chemical ligands for soil compounds, P is generally considered to be highly unavailable in such soils and measured soil pH values (neutral). It is therefore appropriate to question the putative origin of such soluble P concentrations. Our work on the bacterial and fungal communities of the same sites revealed the presence of numerous taxa known to host phosphate-solubilizing microorganisms, and phosphate-solubilizing bacteria have also been successfully isolated (data not shown). These microorganisms may be factors of P release from its unavailable form. However, this would imply that there is no feedback to regulate microbial P solubilizing metabolic activity in the presence of high
phosphate substrate concentrations. However, to the best of our knowledge, there is no data in the literature on this point. In any case, such variations in P content provided the opportunity to examine the potential role of P content as a driver of AMF communities.

## 4.1. High P Contents in the Mining Area do Not Impact AMF Root Colonization

In this study, no difference in AMF colonization rates was observed when comparing P to nP soils, which therefore excluded a negative impact of high P levels on AMF colonization in the reclaimed soils studied.

Yet, a decrease in mycorrhizal colonization when P availability is high has been described by numerous authors [68–75]. Additionally, decreases in the abundance and richness of AMF communities were also observed in soil and roots after P fertilization [76–79]. However, some authors underlined that the decreased AMF colonization upon P fertilization was not systematic [80–84]. Moreover, a moderate amount of P fertilizer could improve AMF diversity while higher amounts of P reduce it [85,86], suggesting that the AMF response to soil P could be related to the P concentration in the soil [87]. On the contrary, Higo et al. [39] observed that the P fertilizer level did not impact either the AMF root colonization, diversity, or community structure.

Such data, therefore, highlight that the AMF response to P is a complex process and underline that the P level is not a reliable predictor of AMF occurrence, root colonization, or community patterns [88,89]. In particular, the AMF response to P may be strongly impacted by sampling time [90], agricultural management [39], or host plant species [28,42,83]. For instance, Tang et al., 2016 [83] showed that for faba bean, an excess of P fertilization did not significantly affect the mycorrhizal rates when compared to a low P and no P treatment. Conversely, the durum wheat root colonization was significantly lower at a high P level, and the arbuscular rate drastically decreased at both low and high P fertilization compared with the unfertilized treatment. In this line, we question the role of the site mining history possibly acting as a stronger driver of AMF communities compared to the P level.

## 4.2. AMF Community Characteristics across Sites: Consequences of Site Mining History Rather Than P Concentrations

In this report, we observed a large taxonomic similarity between AMF communities at the family level whether the soils were enriched in P or not. Glomeraceae was largely prevalent throughout the entire mining ecosystem whereas Archeosporaceae, Diversisporaceae, Paraglomeraceae, and Claroideoglomeraceae were under-represented. Paraglomeraceae and Gigasporaceae were not detected. In addition, we highlighted 26 ASVs shared by all sampling sites, three of which were very abundant.

Several hypotheses may be put forward to explain the low richness or diversity in terms of taxa at the family level: either the absence or under-representation of some AMF families and in contrast, the over-representation of Glomeraceae.

Firstly, Glomeraceae is classically recognized as the most widespread family in global, natural, and managed ecosystems [91,92]. In addition, anthropogenic activities are known to have considerable influence on AMF communities. Notably, intensive agricultural practices, such as soil disturbance by tillage [93] or P fertilization, can adversely affect AMF colonization and diversity [75,88,94,95]. Notably, some *Glomus* spp. are commonly found in soils subjected to fertilization and disturbance whereas others, especially *Scutellospora* spp., are indicative of minimally disturbed soils [96]. In the same way, some results interestingly suggest that soil AMF communities may differ according to land use history and time after land use conversion, suggesting that their diversity and composition may be a legacy of ancestral adaptations to historical habitat and exhibit niche conservatism [33,96]. It has also been suggested, that intensively managed agroecosystems are greatly simplified compared to natural ecosystems. For instance, comparisons between AMF communities in agrosystems managed in different ways showed that families such as Gigasporaceae or Acaulosporaceae are more frequently reduced or eliminated than Glomeraceae [97,98].
In theory, the AMF's life-history strategy (LHS) has been proposed to explain the dynamics, variations, and prevalence of some taxa among AMF communities [99,100]. This concept identifies K-strategists (i.e., AMF competitors that evolved traits to enhance survival in stable environments where competition is high) and r-strategists, or ruderals, that evolved traits to survive under stressful conditions and in disturbed environments, investing their energy mainly in the production of numerous offspring. Notably, Glomeraceae may show an opportunistic behavior similar to r-strategists [101–103]. Thus, taxa resilience and the prevalence of Glomerales in disturbed environments could be explained through abilities such as their higher sporulation, ability to produce high numbers of new hyphae from injured hyphae after mechanical disruptions with healing mechanisms, and anastomosis processes [103–105]. These properties are probably more suited to perturbed environments, and may enable them to spread and reestablish more quickly [106,107]. In addition, it has been shown that representatives of the *Glomus* genus promptly respond to carbon resource shortages with direct investment in the production of spores for reproduction, as is expected for r-strategists [108]. Conversely, Gigasporaceae may resemble the LHS of K-strategists; they mainly spread or colonize through spores from an intact mycelium within a long life cycle [109], which does not facilitate recovery and resilience from disturbance.

With regard to these theoretical assumptions, homogeneity across AMF communities as reported in our study might be explained partially by the landscape history, i.e., firstly by the drastic soil perturbation due to mining activities and secondly, by the restoration and revegetalization of soil during more than one century. Thus, we hypothesize that the soil recolonization by AMF would have taken place with the extensive development of AMF species functionally adapted to rapid recolonization and the stress of high P concentrations. Here, Glomeraceae could be relevant candidates while other AMF families could have been drastically or progressively suppressed, becoming low or nondetectable today under the influence of persistent high P concentrations in the soil for more than a century. In particular, no representative of Gigasporaceae has been identified in the sampled soils. Yet, it has been shown that Gigasporaceae provide functional complementarity to other AMF families in P uptake [110]. In our conditions, due to high soluble P concentrations in the sampling soils, we hypothesize that additional benefits might not have been supplied to the host by such taxa, resulting in their absence. Indeed, the assimilation of P by AMF-colonized plants reflects the sum of P uptake directly by plant cells and indirectly via the AMF pathway [69]. This, however, becomes subsidiary to plants when the P availability is high in the soil matrix and when direct plant uptake is sufficient to provide enough Pi without the help of the indirect mycorrhizal pathway. Consequently, AMF may no longer be essential for P uptake while their other functions (contribution to the nutrient or water supply as well as stress resistance) would remain major, hence, the persistent presence of AMF in P-rich soils. Therefore, we also hypothesize that, besides the predominance of Glomeraceae, the decline or loss of some AMF taxa, and the decrease in overall diversity and afferent mycorrhizal functioning, Glomeraceae recruited by plants may display a large functional diversity and complementation to gain the maximal benefits from the AMF communities in favor of the soil and host plant and for maintaining ecosystem functions.

In line, the occurrence of a representative of *Funneliformis* described in the literature as a ruderal and stress indicator taxon [100] and as a generalist AMF [111] may be compatible with the disturbed land story of sampling sites.

Following our results, the meaning and role filled by the set of persistent ASVs capable of coping with radically contrasting and altered environments can also be questioned. One could hypothetically suggest that this microbial core could correspond to taxa with basic functions (such as nutrient transfer, mainly C and N) and functional features (not yet identified) that confer some level of robustness, resistance, resilience, or functional redundancy in the face of high P content, environmental disturbances, and progressive soil restoration over time [112–114]. However, data about the taxonomic affiliation of the core members are insufficient to go further about their functionalities. Moreover, data in the
literature allowing functional inference to a given AMF species are still limited since the taxonomy and identification of AMF species are still to be refined.

### 4.3. A slight Impact of Soil P on AMF Assemblages Is Resolved at ASV Level

In this study, we identified shared and specific ASVs characterizing the P and nP sites; on the other hand, a lower richness and abundance of AMF communities were observed for the P profile compared to the nP profile. This was particularly pronounced in one site (the S3 site) which was furthermore characterized by the highest concentrations of total and soluble P.

Moreover, through an indicator species analysis, we statistically identified ASVs that are indicative of P and nP sites based on the indicator value calculation as a measure of specificity and fidelity of each ASV to the targeted clustering groups, i.e., P vs. nP site group [63]. Notably, both P- and nP-indicator ASVs were identified but with a lower number of P-indicator ASVs. As a whole, these results clearly differentiate the P and nP sites and therefore highlighted that the AMF community is responsive to P as an edaphic factor in driving AMF communities when analysis is performed with a stringent level of taxonomic resolution.

## 5. Conclusions

In this study, we identified a large similarity of AMF community compositions at the family rank as well as the persistence of a set of ASVs along the drastic P gradient inside the mining area sampled. This leads to the conclusion that AMF community profiles were robust across the contrasting soil P levels. Nevertheless, the shared and specific ASVs from the P vs. nP sites were also identified as well as indicator ASVs of the P vs. nP sites which, therefore, differentiated the P vs. nP sites. As a whole, these results highlighted the AMF community assemblages' responsivity to P but also showed that P is currently not a strong or prevalent ecological pressure exerted on the AMF community composition across the sampling area.

Conversely, our results suggest that contemporary AMF community assemblages may prevalently reflect life-history strategies and resilience of AMF communities, in response to mining history (i.e., land disturbance then ecosystem restoration), favoring Glomerales that behave as ruderals, able to rapidly respond to stressful conditions through their ability to disperse, colonize, and persist in soils and roots over time during reclamation of land-mined soils.

In conclusion, this work opens new opportunities to advance our in-depth understanding of the analysis of P's impact on AMF communities. It also improves the knowledge of AMF diversity in natural P-rich soil, thus providing novel perspectives for the design of AMF inoculants to be used in combination with RP fertilization and the development of more sustainable agroecosystems.

**Supplementary Materials:** The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/microorganisms10122406/s1. Figure S1: A: The remaining material from past open-pit mines, B, C, D: Landscape in localization 1, 2, and 3 respectively. Figure S2: Rarefaction curves obtained from the 18S rRNA gene dataset A: Curves obtained from the 3,527,278 sequences yielded after quality filtering and chimera removal (24 soil samples and 24 root samples). B: Rarefaction curves obtained after a random resampling step, computed to compensate for the uneven sequencing efforts of different samples and to keep an adequate sequencing depth to capture most of the AMF diversity from soil and root samples. Thus, the curves were obtained for a standardized 6000 sequences per sample. Figure S3: taxonomic assignment and RaxML phylogenetic tree of mycorrhizal AVSs. Different colors differentiate each family assigned. Red to pink: Diversisporales; yellow to orange: Paraglomerales; dark orange: Archaeosporales; green: nonassigned; blue and purple: Glomerales; and grey: others. Figure S4: Numbers of specific and shared ASVs comparing (A) the P sites and (B) the nP sites of the three localizations: Low membership characterized by a high number of ASVs specific to one site and few overlapping ASVs. Table S1: Total P and available Pi contents in soils sampled across three localizations (L1, L2, L3); each harboring "P" site covered in the past by PR
deposits or “nP” sites without deposit. Table S2: Distribution of number of reads in each treatment and plant species. Table S3: Taxonomic distribution of the ASVs identified as AMF by the Silva database and those identified de novo by the phylogenetics analysis. Table S4: Taxonomic distribution patterns of the 318 ASVs identified as Glomeromycota at the genus (A) and species (B) taxonomic rank. Table S5: ASVs with assignment up to species taxonomic level (39 ASVs). Table S6: Shared ASVs (26) by the six sampling sites. Table S7: ASVs identified after indicator species analysis.

**Author Contributions:** A.D.-D. conceived and designed the sampling study then operated the data acquisition and analysis. R.R. contributed to carrying out the phylogeny data analysis. J.F., M.H., and A.L.-H.S. contributed to supervision. A.D.-D. drafted the manuscript, and all authors were involved in revision of the final version. All authors have read and agreed to the published version of the manuscript.

**Funding:** This work was carried out within the framework of the ALIBIOTECH and BiHauts Eco de France projects which are funded by the European Union, the French State, and the French Region of Hauts-de-France as well as the TRIPLET project financed by A2U. This research was also funded by the Natural Sciences and Engineering Research Council of Canada (NSERC) to MH, Grant Number RGPIN-2018-04178. The authors wish to thank the "Université de Montréal" and the "Région des Hauts de France" for providing financial support to the ADD PhD thesis.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** The original contributions presented in the study are publicly available.

**Acknowledgments:** The authors wish to thank the "Préfecture du Lot" and the "Réserve Naturelle Nationale d'intérêt géologique du Lot" for the sampling authorization.

**Conflicts of Interest:** The authors declare no conflict of interest.

## References

1. Chen, M.; Arato, M.; Borghi, L.; Nouri, E.; Reinhardt, D. Beneficial Services of Arbuscular Mycorrhizal Fungi—From Ecology to Application. *Front. Plant Sci.* **2018**, *9*, 1270. [CrossRef] [PubMed]

2. van der Heijden, M.G.A.; Martin, F.M.; Selosse, M.A.; Sanders, I.R. Mycorrhizal ecology and evolution: The past, the present, and the future. *New Phytol.* **2015**, *205*, 1406–1423. [CrossRef] [PubMed]

3. Trabelsi, D.; Mhamdi, R. Microbial Inoculants and Their Impact on Soil Microbial Communities: A Review. *BioMed Res. Int.* **2013**, *2013*, 863240. [CrossRef] [PubMed]

4. Wang, B.; Qiu, Y.L. Phylogenetic distribution and evolution of mycorrhizas in land plants. *Mycorrhiza* **2006**, *16*, 299–363. [CrossRef]

5. Smith, S.E.; Read, D.J. *Mycorrhizal Symbiosis*; Academic Press: Cambridge, MA, USA, 2010.

6. Smith, S.E.; Read, D. *Mycorrhizal Symbiosis*, 3rd ed.; Academic Press: Cambridge, MA, USA, 2008.

7. Luginbuehl, L.H.; Menard, G.N.; Kurup, S.; Van Erp, H.; Radhakrishnan, G.V.; Breakspear, A.; Oldroyd, G.E.; Eastmond, P.J.J.S. Fatty acids in arbuscular mycorrhizal fungi are synthesized by the host plant. *Science* **2017**, *356*, 1175–1178. [CrossRef]

8. Liu, G.; Pfeifer, J.; de Brito Francisco, R.; Emonet, A.; Stirnemann, M.; Gübeli, C.; Hutter, O.; Sasse, J.; Mattheyer, C.; Stelzer, E.J.N.P. Changes in the allocation of endogenous strigolactone improve plant biomass production on phosphate-poor soils. *New Phytol.* **2017**, *217*, 784–798. [CrossRef]

9. Wang, D.; Lv, S.; Jiang, P.; Li, Y. Roles, Regulation, and Agricultural Application of Plant Phosphate Transporters. *Front. Plant Sci.* **2017**, *8*, 817. [CrossRef]

10. Weber, S.E.; Diez, J.M.; Andrews, L.V.; Goulden, M.L.; Aronson, E.L.; Allen, M.F.J.F.E. Responses of arbuscular mycorrhizal fungi to multiple coinciding global change drivers. *Fungal Ecol.* **2019**, *40*, 62–71. [CrossRef]

11. Diagne, N.; Ngom, M.; Djighaly, P.I.; Fall, D.; Hocher, V.; Svistoonoff, S. Roles of Arbuscular Mycorrhizal Fungi on Plant Growth and Performance: Importance in Biotic and Abiotic Stressed Regulation. *Diversity* **2020**, *12*, 370. [CrossRef]

12. Ferrol, N.; Tamayo, E.; Vargas, P. The heavy metal paradox in arbuscular mycorrhizas: From mechanisms to biotechnological applications. *J. Exp. Bot.* **2016**, *67*, 6253–6265. [CrossRef]

13. Raveau, R.; Lounès-Hadj Sahraoui, A.; Hijri, M.; Fontaine, J. Clary Sage Cultivation and Mycorrhizal Inoculation Influence the Rhizosphere Fungal Community of an Aged Trace-Element Polluted Soil. *Microorganisms* **2021**, *9*, 1333. [CrossRef] [PubMed]

14. Cavagnaro, T.R.; Bender, S.F.; Asghari, H.R.; van der Heijden, M.G. The role of arbuscular mycorrhizas in reducing soil nutrient loss. *Trends Plant Sci.* **2015**, *20*, 283–290. [CrossRef] [PubMed]

15. Sellal, Z.; Ouazzani, A.; Najoua, M.; Ouarraqi, M.; Karima, S.; Jamila, D.; Rachid, B.; El Modafar, C.; Allal, D. Effect of an Endomycorrhizal Inoculum on the Growth of Argan Tree Plants. *Int. J. Environ. Agric. Biotechnol.* **2017**, *2*, 928–939. [CrossRef]

16. Moradi, M.; Naji, H.R.; Imani, F.; Behbahani, S.M.; Ahmadi, M.T.J.J.o.A.E. Arbuscular mycorrhizal fungi changes by afforestation in sand dunes. *J. Arid. Environ.* **2017**, *140*, 14–19. [CrossRef]
17. Powell, J.R.; Rillig, M.C. Biodiversity of arbuscular mycorrhizal fungi and ecosystem function. *New Phytol.* **2018**, *220*, 1059–1075.
[CrossRef]

18. Bargaz, A.; Lyamlouli, K.; Chtouki, M.; Zeroual, Y.; Dhiba, D. Soil Microbial Resources for Improving Fertilizers Efficiency in an Integrated Plant Nutrient Management System. *Front. Microbiol.* **2018**, *9*, 1606. [CrossRef]

19. Hijri, M. Analysis of a large dataset of mycorrhiza inoculation field trials on potato shows highly significant increases in yield. *Mycorrhiza* **2016**, *26*, 209–214. [CrossRef]

20. Schüßler, A.; Walker, C. *The Glomeromycota*; The Royal Botanic Garden Kew: London, UK, 2010.

21. Redecker, D.; Schüßler, A.; Stockinger, H.; Stürmer, S.L.; Morton, J.B.; Walker, C. An evidence-based consensus for the classification of arbuscular mycorrhizal fungi (Glomeromycota). *Mycorrhiza* **2013**, *23*, 515–531. [CrossRef]

22. Öpik, M.; Davison, J. Uniting species-and community-oriented approaches to understand arbuscular mycorrhizal fungal diversity. *Fungal Ecol.* **2016**, *24*, 106–113. [CrossRef]

23. Opik, M.; Vanatoa, A.; Vanatoa, E.; Moora, M.; Davison, J.; Kalwij, J.M.; Reier, U.; Zobel, M. The online database MaarJAM reveals global and ecosystemic distribution patterns in arbuscular mycorrhizal fungi (Glomeromycota). *New Phytol.* **2010**, *188*, 223–241.
[CrossRef]

24. Roy, J.; Mazel, F.; Sosa-Hernández, M.A.; Dueñas, J.F.; Hempel, S.; Zinger, L.; Rillig, M.C.J.N.P. The relative importance of ecological drivers of arbuscular mycorrhizal fungal distribution varies with taxon phylogenetic resolution. *New Phytol.* **2019**, *224*, 936–948. [CrossRef] [PubMed]

25. Davison, J.; Moora, M.; Semchenko, M.; Adenan, S.B.; Ahmed, T.; Akhmetzhanova, A.A.; Alatalo, J.M.; Al-Quraishy, S.; Andriyanova, E.; Anslan, S.J.N.P. Temperature and pH define the realised niche space of arbuscular mycorrhizal fungi. *New Phytol.* **2021**, *231*, 763–776. [CrossRef] [PubMed]

26. Hazard, C.; Gosling, P.; Van Der Gast, C.J.; Mitchell, D.T.; Doohan, F.M.; Bending, G.D.J.T.I.j. The role of local environment and geographical distance in determining community composition of arbuscular mycorrhizal fungi at the landscape scale. *ISME J.* **2012**, *7*, 498–508. [CrossRef] [PubMed]

27. Větrovský, T.; Kohout, P.; Kopecký, M.; Machac, A.; Man, M.; Bahnmann, B.D.; Brabcová, V.; Choi, J.; Meszárošová, L.; Human, Z.R.; et al. A meta-analysis of global fungal distribution reveals climate-driven patterns. *Nat. Commun.* **2019**, *10*, 5142. [CrossRef]

28. Gosling, P.; Mead, A.; Proctor, M.; Hammond, J.P.; Bending, G.D. Contrasting arbuscular mycorrhizal communities colonizing different host plants show a similar response to a soil phosphorus concentration gradient. *New Phytol.* **2013**, *198*, 546–556. [CrossRef]

29. Zhao, H.; Li, X.; Zhang, Z.; Zhao, Y.; Yang, J.; Zhu, Y.J.P. Species diversity and drivers of arbuscular mycorrhizal fungal communities in a semi-arid mountain in China. *PeerJ* **2017**, *5*, e4155. [CrossRef]

30. Torrecillas, E.; Alguacil, M.M.; Roldán, A. Host Preferences of Arbuscular Mycorrhizal Fungi Colonizing Annual Herbaceous Plant Species in Semiarid Mediterranean Prairies. *Appl. Environ. Microbiol.* **2012**, *78*, 6180–6186. [CrossRef]

31. Paudel, S.; Longcore, T.; MacDonald, B.; McCormick, M.K.; Szlavecz, K.; Wilson, G.W.; Loss, S.R. Belowground interactions with aboveground consequences: Invasive earthworms and arbuscular mycorrhizal fungi. *Ecology* **2016**, *97*, 605–614. [CrossRef]

32. Tedersoo, L.; Bahram, M.; Zobel, M. How mycorrhizal associations drive plant population and community biology. *Science* **2020**, *367*, eaba1223. [CrossRef]

33. Vasar, M.; Davison, J.; Sepp, S.-K.; Oja, J.; Al-Quraishy, S.; Bueno, C.G.; Cantero, J.J.; Fabiano, E.C.; Decocq, G.; Fraser, L.; et al. Global taxonomic and phylogenetic assembly of AM fungi. *Mycorrhiza* **2022**, *32*, 135–144. [CrossRef]

34. Rasmussen, P.U.; Hugerth, L.W.; Blanchet, F.G.; Andersson, A.F.; Lindahl, B.D.; Tack, A.J.M. Multiscale patterns and drivers of arbuscular mycorrhizal fungal communities in the roots and root-associated soil of a wild perennial herb. *New Phytol.* **2018**, *220*, 1248–1261. [CrossRef] [PubMed]

35. Berruti, A.; Bianciotto, V.; Lumini, E. Seasonal variation in winter wheat field soil arbuscular mycorrhizal fungus communities after non-mycorrhizal crop cultivation. *Mycorrhiza* **2018**, *28*, 535–548. [CrossRef]

36. Öpik, M.; Moora, M.; Zobel, M.; Saks, Ü.; Wheatley, R.; Wright, F.; Daniell, T. High diversity of arbuscular mycorrhizal fungi in a boreal herb-rich coniferous forest. *New Phytol.* **2008**, *179*, 867–876. [CrossRef] [PubMed]

37. Bonfim, J.A.; Vasconcellos, R.L.F.; Gumiere, T.; de Lourdes Colombo Mescolotti, D.; Oehl, F.; Nogueira Cardoso, E.J.B. Diversity of Arbuscular Mycorrhizal Fungi in a Brazilian Atlantic Forest Toposequence. *Microb. Ecol.* **2016**, *71*, 164–177. [CrossRef] [PubMed]

38. Bainard, L.D.; Dai, M.; Gomez, E.F.; Torres-Arias, Y.; Bainard, J.D.; Sheng, M.; Eilers, W.; Hamel, C. Arbuscular mycorrhizal fungal communities are influenced by agricultural land use and not soil type among the Chernozem great groups of the Canadian Prairies. *Plant Soil* **2015**, *387*, 351–362. [CrossRef]

39. Higo, M.; Azuma, M.; Kamiyoshihara, Y.; Kanda, A.; Tatewaki, Y.; Isobe, K. Impact of Phosphorus Fertilization on Tomato Growth and Arbuscular Mycorrhizal Fungal Communities. *Microorganisms* **2020**, *8*, 178. [CrossRef]

40. Cofré, N.M.; Ferrari, A.E.; Becerra, A.; Domínguez, L.; Wall, L.G.; Urcelay, C. Effects of cropping systems under no-till agriculture on arbuscular mycorrhizal fungi in Argentinean Pampas. *Soil Use Manag.* **2017**, *33*, 364–378. [CrossRef]

41. Bodenhausen, N.; Somerville, V.; Desirò, A.; Walser, J.-C.; Borghi, L.; van der Heijden, M.G.A.; Schlaeppi, K. Petunia- and Arabidopsis-Specific Root Microbiota Responses to Phosphate Supplementation. *Phytobiomes J.* **2019**, *3*, 112–124. [CrossRef]

42. Gomes, E.A.; Lana, U.G.P.; Quensen, J.F.; de Sousa, S.M.; Oliveira, C.A.; Guo, J.; Guimarães, L.J.M.; Tiedje, J.M. Root-Associated Microbiome of Maize Genotypes with Contrasting Phosphorus Use Efficiency. *Phytobiomes J.* **2018**, *2*, 129–137. [CrossRef]
43. Robbins, C.; Thiergart, T.; Hacquard, S.; Garrido-Oter, R.; Gans, W.; Peiter, E.; Schulze-Lefert, P.; Spaepen, S. Root-Associated Bacterial and Fungal Community Profiles of Arabidopsis thaliana Are Robust Across Contrasting Soil P Levels. *Phytobiomes* **2018**, *2*, 24–34. [CrossRef]

44. Rillig, M.C.; Sosa-Hernández, M.A.; Roy, J.; Aguilar-Trigueros, C.A.; Vályi, K.; Lehmann, A. Towards an Integrated Mycorrhizal Technology: Harnessing Mycorrhiza for Sustainable Intensification in Agriculture. *Front. Plant Sci.* **2016**, *7*, 1625. [CrossRef] [PubMed]

45. Ferrol, N.; Azcón-Aguilar, C.; Pérez-Tienda, J.J.P.S. Arbuscular mycorrhizas as key players in sustainable plant phosphorus acquisition: An overview on the mechanisms involved. *Plant Sci.* **2019**, *280*, 441–447. [CrossRef] [PubMed]

46. Billaud, Y. Les Paragenèses Phosphatées du Paléokarst Des Phosphorites du Quercy. Y. Billaud. 1982. Ph.D. Thesis, University Claude-Bernard (Lyon 1) n° 1162, Karstologia, Villeurbanne, France, 1983. Volume 1, no 1. p. 77.

47. Phillips, J.M.; Hayman, D.S. Improved procedures for clearing roots and staining parasitic and vesicular-arbuscular mycorrhizal fungi for rapid assessment of infection. *Trans. Br. Mycol. Soc.* **1970**, *55*, 158-IN18. [CrossRef]

48. Trouvelot, A.; Kough, J.L.; Gianinazzi-Pearson, V. Estimation of VA mycorrhizal infection levels. Research for methods having a functional significance. In Proceedings of the European Symposium on Mycorrhizae, Dijon, France, 1–5 July 1985.

49. Abu-Romman, S. Comparison of methods for isolating high quality DNA from sage (Salvia officinalis). *J. Med. Plants Res.* **2011**, *5*, 938–941.

50. Aleksic, J.M.; Stojanovic, D.; Banovic, B.; Jancic, R. A simple and efficient DNA isolation method for Salvia officinalis. *Biochem. Genet.* **2012**, *50*, 881–892. [CrossRef]

51. Renaut, S.; Daoud, R.; Masse, J.; Vialle, A.; Hijri, M. Inoculation with Rhizophagus irregularis Does Not Alter Arbuscular Mycorrhizal Fungal Community Structure within the Roots of Corn, Wheat, and Soybean Crops. *Microorganisms* **2020**, *8*, 83. [CrossRef]

52. Lee, J.; Lee, S.; Young, J.P. Improved PCR primers for the detection and identification of arbuscular mycorrhizal fungi. *FEMS Microbiol. Ecol.* **2008**, *65*, 339–349. [CrossRef]

53. Stefani, F.; Bencherif, K.; Sabourin, S.; Hadj-Sahraoui, A.L.; Banchini, C.; Seguin, S.; Dalpe, Y. Taxonomic assignment of arbuscular mycorrhizal fungi in an 18S metagenomic dataset: A case study with saltcedar (Tamarix aphylla). *Mycorrhiza* **2020**, *30*, 243–255. [CrossRef]

54. Callahan, B.J.; McMurdie, P.J.; Rosen, M.J.; Han, A.W.; Johnson, A.J.; Holmes, S.P. DADA2: High-resolution sample inference from Illumina amplicon data. *Nat. Methods* **2016**, *13*, 581–583. [CrossRef]

55. Callahan, B.; McMurdie, P.; Holmes, S. Exact sequence variants should replace operational taxonomic units in marker-gene data analysis. *ISME J.* **2017**, *11*, 2639–2643. [CrossRef]

56. Callahan, B.J.S. Silva taxonomic training data formatted for DADA2. *Silva Version* **2018**, *132*. [CrossRef]

57. Lassmann, T.; Frings, O.; Sonnhammer, E.L.L. Kalign2: High-performance multiple alignment of protein and nucleotide sequences allowing external features. *Nucleic Acids Res.* **2009**, *37*, 858–865. [CrossRef] [PubMed]

58. Stamatakis, A. RAxML version 8: A tool for phylogenetic analysis and post-analysis of large phylogenies. *Bioinformatics* **2014**, *30*, 1312–1313. [CrossRef] [PubMed]

59. Miller, M.A.; Pfeiffer, W.; Schwartz, T. The CIPRES science gateway: A community resource for phylogenetic analyses. In Proceedings of the 2011 TeraGrid Conference: Extreme Digital Discovery, Salt Lake City, UT, USA, 18–21 July 2011; p. 41.

60. Kim, B.R.; Shin, J.; Guevarra, R.; Lee, J.H.; Kim, D.W.; Seol, K.H.; Lee, J.H.; Kim, H.B.; Isaacson, R. Deciphering Diversity Indices for a Better Understanding of Microbial Communities. *J. Microbiol. Biotechnol.* **2017**, *27*, 2089–2093. [CrossRef] [PubMed]

61. Shannon, C.E. A Mathematical Theory of Communication. *Bell Syst. Tech. J.* **1948**, XXVII, 379–423. [CrossRef]

62. Simpson, E.H. Measurement of diversity. *Nature* **1949**, *163*, 688. [CrossRef]

63. De Cáceres, M.; Legendre, P. Associations between species and groups of sites: Indices and statistical inference. *Ecology* **2009**, *90*, 3566–3574. [CrossRef]

64. De Cáceres, M.; Legendre, P.; Moretti, M.J.O. Improving indicator species analysis by combining groups of sites. *Oikos* **2010**, *119*, 1674–1684. [CrossRef]

65. Frossard, E.; Condron, L.M.; Oberson, A.; Sinaj, S.; Fardeau, J.C. Processes governing phosphorus availability in temperate soils. *J. Environ. Qual.* **2000**, *29*, 15–23. [CrossRef]

66. Francioli, D.; Schulz, E.; Lentendu, G.; Wubet, T.; Buscot, F.; Reitz, T. Mineral vs. Organic Amendments: Microbial Community Structure, Activity and Abundance of Agriculturally Relevant Microbes Are Driven by Long-Term Fertilization Strategies. *Front. Microbiol.* **2016**, *7*, 1446. [CrossRef]

67. Niederberger, J.; Kohler, M.; Bauhus, J. Distribution of phosphorus fractions with different plant availability in German forest soils and their relationship with common soil properties and foliar P contents. *Soil* **2019**, *5*, 189–204. [CrossRef]

68. Nouri, E.; Breuillin-Sessoms, F.; Feller, U.; Reinhardt, D. Phosphorus and nitrogen regulate arbuscular mycorrhizal symbiosis in Petunia hybrida. *PLoS ONE* **2014**, *9*, e90841. [CrossRef] [PubMed]

69. Smith, S.E.; Smith, F.A. Roles of arbuscular mycorrhizas in plant nutrition and growth: New paradigms from cellular to ecosystem scales. *Annu. Rev. Plant Biol.* **2011**, *62*, 227–250. [CrossRef] [PubMed]

70. Olsson, P.A.; Baath, E.; Jakobsen, I. Phosphorus effects on the mycelium and storage structures of an arbuscular mycorrhizal fungus as studied in the soil and roots by analysis of Fatty Acid signatures. *Appl. Environ. Microbiol.* **1997**, *63*, 3531–3538. [CrossRef] [PubMed]
71. Breuillin, F.; Schramm, J.; Hajirezaei, M.; Ahkami, A.; Favre, P.; Druege, U.; Hause, B.; Bucher, M.; Kretzschmar, T.; Bossolini, E.; et al. Phosphate systemically inhibits development of arbuscular mycorrhiza in Petunia hybrida and represses genes involved in mycorrhizal functioning. *Plant J.* **2010**, *64*, 1002–1017. [CrossRef] [PubMed]

72. Kobae, Y.; Ohtomo, R. An improved method for bright-field imaging of arbuscular mycorrhizal fungi in plant roots. *Soil Sci. Plant Nutr.* **2016**, *62*, 27–30. [CrossRef]

73. Balzergue, C.; Puech-Pagès, V.; Bécard, G.; Rochange, S.F. The regulation of arbuscular mycorrhizal symbiosis by phosphate in pea involves early and systemic signalling events. *J. Exp. Bot.* **2011**, *62*, 1049–1060. [CrossRef]

74. Peine, M.; Vitow, N.; Grafe, M.; Baum, C.; Zicker, T.; Eichler-Löbermann, B.; Schulz, S.; Schloter, M.; Leinweber, P. Effect of triple superphosphate and biowaste compost on mycorrhizal colonization and enzymatic P mobilization under maize in a long-term field experiment. *J. Plant Nutr. Soil Sci.* **2019**, *182*, 167–174. [CrossRef]

75. Williams, A.; Manoharan, L.; Rosenstock, N.P.; Olsson, P.A.; Hedlund, K. Long-term agricultural fertilization alters arbuscular mycorrhizal fungal community composition and barley (Hordeum vulgare) mycorrhizal carbon and phosphorus exchange. *New Phytol.* **2016**, *213*, 874–885. [CrossRef]

76. Liu, L.; Gundersen, P.; Zhang, T.; Mo, J. Effects of phosphorus addition on soil microbial biomass and community composition in three forest types in tropical China. *Soil Biol. Biochem.* **2012**, *44*, 31–38. [CrossRef]

77. Camenzind, T.; Hempel, S.; Homeier, J.; Horn, S.; Velescu, A.; Wilcke, W.; Rillig, M.C. Nitrogen and phosphorus additions impact arbuscular mycorrhizal abundance and molecular diversity in a tropical montane forest. *Glob. Chang. Biol.* **2014**, *20*, 3646–3659. [CrossRef] [PubMed]

78. Chen, Y.-L.; Zhang, X.; Ye, J.-S.; Han, H.-Y.; Wan, S.-Q.; Chen, B.-D. Six-year fertilization modifies the biodiversity of arbuscular mycorrhizal fungi in a temperate steppe in Inner Mongolia. *Soil Biol. Biochem.* **2014**, *69*, 371–381. [CrossRef]

79. Lin, X.; Feng, Y.; Zhang, H.; Chen, R.; Wang, J.; Zhang, J.; Chu, H. Long-Term Balanced Fertilization Decreases Arbuscular Mycorrhizal Fungal Diversity in an Arable Soil in North China Revealed by 454 Pyrosequencing. *Environ. Sci. Technol.* **2012**, *46*, 5764–5771. [CrossRef] [PubMed]

80. Wakelin, S.; Mander, C.; Gerard, E.; Jansa, J.; Erb, A.; Young, S.; Condron, L.; O’Callaghan, M. Response of soil microbial communities to contrasted histories of phosphorus fertilisation in pastures. *Appl. Soil Ecol.* **2012**, *61*, 40–48. [CrossRef]

81. Treseder, K.K. The extent of mycorrhizal colonization of roots and its influence on plant growth and phosphorus content. *Plant Soil* **2013**, *371*, 1–13. [CrossRef]

82. Fiorilli, V.; Vallino, M.; Biselli, C.; Faccio, A.; Bagnaresi, P.; Bonfante, P. Host and non-host roots in rice: Cellular and molecular approaches reveal differential responses to arbuscular mycorrhizal fungi. *Front. Plant Sci.* **2015**, *6*, 636. [CrossRef]

83. Tang, X.; Placella, S.A.; Daydé, F.; Bernard, L.; Robin, A.; Journet, E.-P.; Justes, E.; Hinsinger, P. Phosphorus availability and microbial community in the rhizosphere of intercropped cereal and legume along a P-fertilizer gradient. *Plant Soil* **2016**, *407*, 119–134. [CrossRef]

84. Yu, P.; Wang, C.; Baldauf, J.A.; Tai, H.; Gutjahr, C.; Hochholdinger, F.; Li, C. Root type and soil phosphate determine the taxonomic landscape of colonizing fungi and the transcriptome of field-grown maize roots. *New Phytol.* **2018**, *217*, 1240–1253. [CrossRef]

85. Higo, M.; Takahashi, Y.; Gunji, K.; Isobe, K. How are arbuscular mycorrhizal associations related to maize growth performance during short-term cover crop rotation? *J. Sci. Food Agric.* **2017**, *98*, 1388–1396. [CrossRef]

86. Kahiluoto, H.; Ketoja, E.; Vestberg, M.; Saarela, I. Promotion of AM utilization through reduced P fertilization 2. Field studies. *Plant Soil* **2001**, *231*, 65–79. [CrossRef]

87. Liu, M.; Liu, X.; Cheng, B.-S.; Ma, X.-L.; Lyu, X.-T.; Zhao, X.-F.; Ju, Y.-L.; Min, Z.; Zhang, Z.-W.; Fang, Y.-I. Selection and evaluation of phosphate-solubilizing bacteria from grapevine rhizospheres for use as biofertilizers. *Span. J. Agric. Res.* **2016**, *14*, e1106. [CrossRef]

88. van der Heijden, M.G.A.; Dombrowski, N.; Schlaeppi, K. Continuum of root-fungal symbioses for plant nutrition. *Proc. Natl. Acad. Sci. USA* **2017**, *114*, 11574–11576. [CrossRef] [PubMed]

89. Venegas, R.A.P.; Lee, S.-J.; Thuita, M.; Mlay, D.P.; Masso, C.; Vanlauwe, B.; Rodriguez, A.; Sanders, I.R. The phosphate inhibition paradigm: Host and fungal genotypes determine arbuscular mycorrhizal fungal colonization and responsiveness to inoculation in cassava with increasing phosphorus supply. *Front. Plant Sci.* **2021**, *12*, 693037. [CrossRef] [PubMed]

90. Higo, M.; Isobe, K.; Kondo, T.; Yamaguchi, M.; Takeyama, S.; Drijber, R.A.; Torigoe, Y. Temporal variation of the molecular diversity of arbuscular mycorrhizal communities in three different winter cover crop rotational systems. *Biol. Fertil. Soils* **2015**, *51*, 21–32. [CrossRef]

91. Oehl, F.; Laczko, E.; Bogenrieder, A.; Stahr, K.; Bösch, R.; van der Heijden, M.; Sieverding, E. Soil type and land use intensity determine the composition of arbuscular mycorrhizal fungal communities. *Soil Biol. Biochem.* **2010**, *42*, 724–738. [CrossRef]

92. Brearley, F.Q.; Elliott, D.R.; Iribar, A.; Sen, R. Arbuscular mycorrhizal community structure on co-existing tropical legume trees in French Guiana. *Plant Soil* **2016**, *403*, 253–265. [CrossRef]

93. Säle, V.; Aguilera, P.; Laczko, E.; Mäder, P.; Berner, A.; Zihlmann, U.; van der Heijden, M.G.A.; Oehl, F. Impact of conservation tillage and organic farming on the diversity of arbuscular mycorrhizal fungi. *Soil Biol. Biochem.* **2015**, *84*, 38–52. [CrossRef]

94. Bhadalung, N.N.; Suwanarit, A.; Dell, B.; Nopamornbodi, O.; Thamchaipenet, A.; Rungchuang, J. Effects of long-term NP-fertilization on abundance and diversity of arbuscular mycorrhizal fungi under a maize cropping system. *Plant Soil* **2005**, *270*, 371–382. [CrossRef]
95. García de León, D.; Davison, J.; Moora, M.; Öpik, M.; Feng, H.; Hiiesalu, I.; Jairus, T.; Koorem, K.; Liu, Y.; Phosri, C.; et al. Anthropogenic disturbance equalizes diversity levels in arbuscular mycorrhizal fungal communities. *Glob. Chang. Biol.* **2018**, *24*, 2649–2659. [CrossRef]

96. Faggioli, V.S.; Cabello, M.N.; Grilli, G.; Vasar, M.; Covacevich, F.; Öpik, M. Root colonizing and soil borne communities of arbuscular mycorrhizal fungi differ among soybean fields with contrasting historical land use. *Agric. Ecosyst. Environ.* **2019**, *269*, 174–182. [CrossRef]

97. Gosling, P.; Jones, J.; Bending, G.D.J.M. Evidence for functional redundancy in arbuscular mycorrhizal fungi and implications for agroecosystem management. *Mycorrhiza* **2015**, *26*, 77–83. [CrossRef] [PubMed]

98. Higo, M.; Isobe, K.; Miyazawa, Y.; Matsuda, Y.; Drijber, R.A.; Torigoe, Y. Molecular diversity and distribution of indigenous arbuscular mycorrhizal communities colonizing roots of two different winter cover crops in response to their root proliferation. *J. Microbiol.* **2016**, *54*, 86–97. [CrossRef] [PubMed]

99. Hart, M.M.; Reader, R.J.; Klironomos, J.N.J.M. Life-history strategies of arbuscular mycorrhizal fungi in relation to their successional dynamics. *Mycologia* **2001**, *93*, 1186–1194. [CrossRef]

100. Chagnon, P.L.; Bradley, R.L.; Maherali, H.; Klironomos, J.N. A trait-based framework to understand life history of mycorrhizal fungi. *Trends Plant Sci.* **2013**, *18*, 484–491. [CrossRef]

101. de Souza, F.A.; Dalpé, Y.; Declerck, S.; de la Providencia, I.E.; Séjalon-Delmas, N. Life history strategies in Gigasporaceae: Insight from monoxenic culture. In *In Vitro Culture of Mycorrhizas*; Springer: Berlin/Heidelberg, Germany, 2005; pp. 73–91.

102. De La Providencia, I.E.; De Souza, F.A.; Fernández, F.; Delmas, N.S.; Declerck, S. Arbuscular mycorrhizal fungi reveal distinct patterns of anastomosis formation and hyphal healing mechanisms between different phylogenetic groups. *New Phytol.* **2004**, *165*, 261–271. [CrossRef]

103. De La Providencia, I.E.; Fernández, F.; Declerck, S. Hyphal healing mechanism in the arbuscular mycorrhizal fungi Scutellospora reticulata and Glomus clarum differs in response to severe physical stress. *FEMS Microbiol. Lett.* **2007**, *268*, 120–125. [CrossRef]

104. Giovannetti, M.; Azzolini, D.; Citernesi, A.S. Anastomosis Formation and Nuclear and Protoplasmic Exchange in Arbuscular Mycorrhizal Fungi. *Appl. Environ. Microbiol.* **1999**, *65*, 5571–5575. [CrossRef]

105. Simard, S.W.; Beiler, K.J.; Bingham, M.A.; Deslippe, J.R.; Philip, L.J.; Teste, F.P. Mycorrhizal networks: Mechanisms, ecology and modelling. *Fungal Biol. Rev.* **2012**, *26*, 39–60. [CrossRef]

106. Allison, S.D.; Martiny, J.B.H. Resistance, resilience, and redundancy in microbial communities. *Proc. Natl. Acad. Sci. USA* **2008**, *105*, 11512–11519. [CrossRef]

107. Oehl, F.; Sieverding, E.; Ineichen, K.; Mäder, P.; Boller, T.; Wiemken, A. Impact of Land Use Intensity on the Species Diversity of Arbuscular Mycorrhizal Fungi in Agroecosystems of Central Europe. *Appl. Environ. Microbiol.* **2003**, *69*, 2816–2824. [CrossRef]

108. IJdo, M.; Schtickzelle, N.; Cranenbrouck, S.; Declerck, S.J.F.m.e. Do arbuscular mycorrhizal fungi with contrasting life-history strategies differ in their responses to repeated defoliation? *FEMS Microbiol. Ecol.* **2010**, *72*, 114–122. [CrossRef] [PubMed]

109. Schalamuk, S.; Cabello, M. Arbuscular mycorrhizal fungal propagules from tillage and no-tillage systems: Possible effects on Glomeromycota diversity. *Mycologia* **2010**, *102*, 261–268. [CrossRef] [PubMed]

110. Tian, H.; Drijber, R.A.; Li, X.; Miller, D.N.; Wienhold, B.J.J.M. Arbuscular mycorrhizal fungi differ in their ability to regulate the expression of phosphate transporters in maize (*Zea mays* L.). *Mycorrhiza* **2013**, *23*, 507–514. [CrossRef] [PubMed]

111. Öpik, M.; Moora, M.; Liira, J.; Zobel, M. Composition of root-colonizing arbuscular mycorrhizal fungal communities in different ecosystems around the globe. *J. Ecol.* **2006**, *94*, 778–790. [CrossRef]

112. Jiao, S.; Chen, W.; Wei, G.J.A. Resilience and assemblage of soil microbiome in response to chemical contamination combined with plant growth. *Appl. Environ. Microbiol.* **2019**, *85*, e02523-18. [CrossRef]

113. Griffiths, B.S.; Philippot, L.J. Insights into the resistance and resilience of the soil microbial community. *FEMS Microbiol. Rev.* **2013**, *37*, 112–129. [CrossRef]

114. Hodgson, D.; McDonald, J.L.; Hosken, D.J.J. What do you mean, 'resilient'? *Trends Ecol. Evol.* **2015**, *30*, 503–506. [CrossRef]