Coriander (Coriandrum sativum) Cultivation Combined with Arbuscular Mycorrhizal Fungi Inoculation and Steel Slag Application Influences Trace Elements-Polluted Soil Bacterial Functioning

The cultivation of aromatic plants for the extraction of essential oils has been presented as an innovative and economically viable alternative for the remediation of areas polluted with trace elements (TE). Therefore, this study focuses on the contribution of the cultivation of coriander and the use of arbuscular mycorrhizal fungi (AMF) in combination with mineral amendments (steel slag) on the bacterial function of the rhizosphere, an aspect that is currently poorly understood and studied. The introduction of soil amendments, such as steel slag or mycorrhizal inoculum, had no significant effect on coriander growth. However, steel slag changed the structure of the bacterial community in the rhizosphere without affecting microbial function. In fact, Actinobacteria were significantly less abundant under slag-amended conditions, while the relative proportion of Gemmatimonadota increased. On the other hand, the planting of coriander affects the bacterial community structure and significantly increased the bacterial functional richness of the amended soil. Overall, these results show that planting coriander most affected the structure and functioning of bacterial communities in the TE-polluted soils and reversed the effects of mineral amendments on rhizosphere bacterial communities and their activities. This study highlights the potential of coriander, especially in combination with steel slag, for phytomanagement of TE-polluted soils, by improving soil quality and health.


Introduction
Soils provide fundamental ecosystem functions and play an important role in water and nutrient cycling, food production and the provision of renewable materials. Their ability to store carbon is essential for coping with climate change, and soil biodiversity is critical for soil fertility and biodiversity in general. Unfortunately, increasing pollution has made large areas of soil unsuitable for cultivation, a hazard to wildlife and humans, and uninhabitable for many living microorganisms. It is widely recognized that trace elements (TE), such as lead (Pb), cadmium (Cd), and zinc (Zn), can inhibit biological processes in ecosystems and affect the biodiversity, abundance, activity and function of soil microorganisms. These pollutants originate from human activities and, in particular, from smelters that have contributed to the production of significant amounts of dusts, resulting in severe contamination of soils with TE [1]. They represent an increased toxicity for microorganisms living in the soil and are also responsible for soil acidification, small size of soil aggregates and poor soil structure due to the lack of organic matter (OM), humic acids and, by extension, nutrients [2][3][4]. Therefore, it is crucial to reduce TE-toxicity in polluted soils to restore soil ecosystem services. Phytomanagement technologies, such as phytostabilisation, are suitable options to manage agricultural soils polluted with TE ecologically and economically [5][6][7]. Therefore, a crucial point for successful phytomanagement is the selection of TE-tolerant plant species that grow rapidly and produce enough biomass to allow recovery [6,8]. Previous studies have shown that the cultivation of aromatic plants, such as coriander (Coriandrum sativum), should be a promising option for phytomanagement of TE-polluted soils [6,9]. Indeed, this aromatic belonging to the Apiaceae family is cultivated for the production of essential oils (EO), which are an important way of valorizing plant biomass. The EO produced on TE-polluted soils have antibacterial, antifungal and antioxidant activities, and can be used for various applications [10]. Although it has been demonstrated that coriander can grow in TE-polluted soils, this plant species accumulates Zn and Cd in its aboveground parts, which makes it impossible to recycle residues from the distillation of coriander (shoots or seeds), in the methanization process, or as compost [6,9,11].
For example, in the context of circular economy, the incorporation into the soil of mineral basic amendments could improve soil parameters, such as pH, nutrient mobilization or soil structure and microbial composition [12][13][14]. In addition, mineral basic amendments, such as slags, can contribute to mineral nutrition and thus improve crop yields, but also play an important role in stabilizing soil aggregates and resistance to plant diseases [12,13,15]. Some studies have shown the role of slag in TE immobilization, especially Cd, Zn and Pb [13,16,17]. Furthermore, steel slag has been described to increase microbial biomass, soil enzyme activity and the number of copiotrophic bacterial taxa, and is believed to provide beneficial ecosystem services to plants [14].
Therefore, an adequate supply of mineral, organic or biological amendments could be interesting, first to reduce TE-mobility and bioavailability, and also to improve the growth of coriander and soil function [11].
Biological amendments, such as arbuscular mycorrhizal fungi (AMF) inoculum, have also been shown to be useful for plant protection. Indeed, AMF enable host plants to grow vigorously under stressful conditions, especially in TE-polluted soils [18]. They contribute significantly to the uptake of soil nutrients, especially phosphorus [19], but also nitrogen, calcium, sulfur, potassium and zinc, thereby increasing plant biomass and improving plant tolerance to soil pollution [20,21]. In addition, AMF produce glomalin, an insoluble glycoprotein that has soil binding properties, but also melanin and chitin, which reduce metal bioavailability by metal immobilization into hyphae [18,22,23]. Soil microbial communities are indeed crucial for soil structure, OM degradation and C cycling [3,24]. They are often used to assess the degree of soil remediation [8,25]. Many recent studies have described a synergistic effect of AMF inoculation and biological and mineral amendments (olive residues, biochar, phosphate sludge from marble waste, compost, manure) on plant growth in TE-polluted soils [26][27][28][29][30]. So far, there are few scientific works dealing with the combined application of AMF and mineral amendments (steel slag), and their effects on microbial communities and their activities are generally poorly studied [16].
Therefore, under in situ conditions, this study focuses on the potential benefits of adding amendments (mineral amendment alone or in combination with AMF inoculum) to improve coriander growth and microbial function of historically TE-polluted soils. To this end, we investigated plant growth parameters, soil bacterial biomass (phospholipid fatty acid, PLFA), the functional activity (community-level physiological profile) and structure (16S rRNA metabarcoding) of rhizosphere bacterial communities.

Coriander Growth Parameters
Although the plants had the same height after 11 weeks of cultivation regardless of the modality, we found that plants inoculated with AMF had a lighter aboveground dry weight than those non-inoculated (Table 1). In addition, the type and dose of the amendment did not alter the dry weight of coriander. Table 1. Growth parameters and mycorrhizal rate of coriander. Data are presented as means ± standard error. The means were obtained from five replicates (n = 5) for height and dry weight of plants, and three replicates (n = 3) for mycorrhization rate. Means followed by the same letter are not significantly different, by one-way ANOVA comparison and post-hoc Fisher test (α = 0.05).

Mycorrhizal Rate
The mycorrhizal rate of coriander roots ranged from 31.7 to 51%, and was significantly higher in plots amended with 1.5 t ha −1 oven slag (51%) than in plots amended with 1.5 t ha −1 ladle slag (31.7%). No difference was observed between NI and I conditions, indicating the presence of native AMF. The dose of amendment did not affect the mycorrhization rate.

Soil Microbial Biomass
Total bacterial biomass was not significantly different between the two unamended conditions. Analysis of variance (ANOVA) of the Gram-positive bacteria biomass showed an increase with the lower dose of amendment (1.5 t ha −1 ) and with a combined effect of dose (1.5 t ha −1 ) and vegetation (Table 2). Moreover, the biomass of Gram-positive bacteria was significantly higher (8.14 µg g −1 of dry weight (DW) of soil) in soil amended with 1.5 t ha −1 oven slag and under vegetated conditions than under the same non-vegetated conditions (2.54 µg g −1 of DW of soil) ( Table 2). A significant effect of vegetation was observed on Gram-negative bacterial biomass. In fact, this increase was significantly higher in vegetated soils and NI condition supplemented with 1.5 t ha −1 ladle slag (5.30 µg g −1 soil) than under unvegetated modality (2.31 µg g −1 soil, Table 2). Table 2. Influence of amendments (ladle slag and oven slag used at 1.5 and 45 t ha −1 ), inoculation and vegetation on Gram-positive, Gram-negative and total bacterial biomass (µg g −1 of dry weight (DW) soil), and Gram-positive/Gram-negative ratio in trace element-polluted soil. Data are presented as means ± standard error. The means were obtained from three replicates (n = 3). Means followed by the same letter are not significantly different, by one-way ANOVA comparison and post-hoc Fisher test (α = 0.05). The total bacterial biomass followed the same trend as the Gram-positive bacterial biomass. Indeed, a significant increase in total bacterial biomass was observed in the amended soil with 1.5 t ha −1 oven slag and under vegetated conditions (11.92 µg g −1 of DW soil) than in the same non-vegetated plots (5.05 µg g −1 of DW soil). Examination of the ratio of the biomass of Gram-positive to Gram-negative bacteria showed no significant difference between diverse modalities ( Table 2).

Bacterial α-Diversity of Soil Biotopes
To further investigate the effects of slags, vegetation and AMF inoculation on rhizospheric bacterial communities, the bacterial OTUs (Operational Taxonomic Units) were taxonomically assigned. The bacterial OTUs in the soil were assigned to 993 genera belonging to 33 phyla. Actinobacteria was the most represented phylum in each modality and accounted for between 41 and 67% of the total bacterial abundance ( Figure 1). Actinobacteria associated with Proteobacteria, Gemmatimonadota, Myxococcota and Acidobacteriota accounted for over 85% of bacterial OTUs. In general, we found that the abundance of the different bacterial phyla, and thus the structure of the bacterial community, was affected by the conditions tested. One of the most important changes concerned the high abundance of Actinobacteria, which was even more pronounced under the vegetated conditions, than under the unvegetated conditions, except for the not inoculated and unamended condition (Table S2). In contrast, a decrease of the Gemmatimonadota and Acidobacteriota abundances was observed under all vegetated conditions ( Figure 1).
In addition to these changes in the abundance of the bacterial phyla, we also noted changes within phyla ( Figure 2  Bacterial α-diversity was calculated independently for each modality. According to the Chao1 index, soil bacterial species richness remains stable across modalities (average 325.3). Similarly, both bacterial richness and diversity (Shannon index) were not affected by the addition of slag, mycorrhizal inoculation or vegetation (Table 3). Table 3. Influence of amendments (ladle slag and oven slag used at 1.5 and 45 t ha −1 ), inoculation and vegetation on bacterial richness and diversity indexes. Data are mean ± SD (n = 3) for each condition. Means followed by the same lowercase letter are not significantly different, by one-way ANOVA comparison and post-hoc Fisher test (α = 0.05). Actinobacteria associated with Proteobacteria, Gemmatimonadota, Myxococcota and Acidobacteriota accounted for over 85% of bacterial OTUs. In general, we found that the abundance of the different bacterial phyla, and thus the structure of the bacterial community, was affected by the conditions tested. One of the most important changes concerned the high abundance of Actinobacteria, which was even more pronounced under the vegetated conditions, than under the unvegetated conditions, except for the not inoculated and unamended condition (Table S2). In contrast, a decrease of the Gemmatimonadota and Acidobacteriota abundances was observed under all vegetated conditions (Figure 1).

State
In addition to these changes in the abundance of the bacterial phyla, we also noted changes within phyla ( Figure 2). For example, within the phylum Actinobacteria, we found that the genera Pseudarthrobacter, Arthrobacter and Streptomyces were enriched under vegetated conditions, whereas the genus Gemmatimonadota was favored under nonvegetated conditions.

Metabolic Potential of Soil Microbial Communities
To understand the impact of taxonomic changes on the functionality of these microbial communities, a study was conducted on their ability to degrade different carbon sources. Overall, the metabolic potential of the soil microbial communities of each sample, expressed as average well color development (AWCD), was significantly higher under vegetated conditions than under the unvegetated conditions. The lowest value was found under the not inoculated, unvegetated and oven slag amended condition at 45 t ha −1 . In addition, there was no significant effect of inoculation, type of amendment or dose on AWCD level (Table 4).  Figure 2. Abundances of the 25 most represented bacterial genus in soil. NI, non-inoculated; I, inoculated.

Metabolic Potential of Soil Microbial Communities
To understand the impact of taxonomic changes on the functionality of these microbial communities, a study was conducted on their ability to degrade different carbon sources. Overall, the metabolic potential of the soil microbial communities of each sample, expressed as average well color development (AWCD), was significantly higher under vegetated conditions than under the unvegetated conditions. The lowest value was found under the not inoculated, unvegetated and oven slag amended condition at 45 t ha −1 . In addition, there was no significant effect of inoculation, type of amendment or dose on AWCD level (Table 4).
A multivariate ANOVA showed that functional richness was significantly higher in vegetated conditions than in the unvegetated conditions. Neither the inoculation, nor the type or dosage of amendment had a significant influence on the values of functional richness (Table 4).

Discussion
The cultivation of aromatic plants for the EO extraction has been presented as an innovative and economically viable alternative for the phytoremediation of TE-polluted areas [6]. However, the effects of growing aromatic plants, such as coriander, in combination with mineral (steel slag) and microbial (AMF) amendments introduced into the soil on the structure and potential function of bacterial communities in the TE-polluted rhizosphere, have never been studied. Thus, our results showed that coriander was able to establish and develop well on heavily TE-polluted soil, confirming previous work [6]. The addition of steel slag had no overall effect on plant growth. This result could be explained by the low nutrient requirement of coriander, but also by the OM content of this agricultural soil [6,31]. Consequently, it is likely that the slag did not provide any supplemental nutrient gain to the plant. Indeed, the beneficial effect of slag is controversial. Some studies have shown that these mineral amendments could contribute to mineral supply and consequently improve crop yields, while others showed no change or diverse responses, depending on the experimental conditions [13,15,[32][33][34].
We have shown that coriander was able to form a mycorrhizal symbiosis with the native AMF present in our TE-polluted soil. Therefore, no significant difference in root mycorrhizal colonization rates was found between NI and I conditions, while AMF inoculation had no significant effect on shoot growth. In contrast to our results, some studies have shown that the introduction of AMF inoculum into TE-polluted and unpolluted soils can lead to a significant improvement in overall mycorrhizal rates, indicating different responses in experiments depending on the plant species, but also on the type of mycorrhizal inoculum [6,35,36]. We cannot exclude that the natural presence of native AMF in the TE-polluted soils masked the contribution of the exogenous mycorrhizal inoculum on coriander growth. In contrast, an increase in mycorrhizal colonization rates was observed in the presence of oven slag. This increase could be due to a reduction in TE-bioavailability or changes in pH or soil texture, as suggested by Hu et al. [16].
Although no benefits were found for the growth of coriander after mineral amendment and inoculation with AMF, several studies reported the crucial role played by these mineral and microbial amendments, as well as vegetation, in restoring soil functions that might be disturbed by TE-pollution [16,26]. In our results, the high bacterial diversity and richness indexes were comparable to some indexes obtained on unpolluted soils [37,38]. Moreover, despite the soil TE pollution, the functionality of the AWCD-based microbial communities indicated that the studied soils had a high metabolic capacity, which is typical of functional soils [39][40][41][42]. However, numerous studies have also shown a reduction in the abundance, activity, and diversity of soil microbial communities following TE exposure [42][43][44]. For example, Fontaine et al. [45] demonstrated a relatively low ratio of fungi-to-bacteria in this agricultural polluted soil, indicating an extensively managed soil in combination with tillage, high fertilization and a low C:N ratio favored bacteria. In this study, more Grampositive than Gram-negative bacteria were observed under all conditions. Adaptation to an environment rich in TE leads microorganisms to exhibit activities of biosorption, bioprecipitation, extracellular sequestration, transport mechanisms and/or chelation [46]. Gram-positive bacteria are known for their diverse physiological and metabolic properties that allow them to thrive in a wide range of environments. Compared to Gram-negative bacteria, they appear to accumulate higher TE-concentrations on their cell walls [47]. Therefore, these bacteria are able to adapt easier to adverse conditions, such as pollution, and should even be used for TE-removal, which could explain their frequent presence in this soil [44,48]. Specifically, Actinobacteria, Proteobacteria and Gemmatimonadota, were the most abundant bacterial phyla in our TE-polluted soils. The Actinobacteria phylum has been highlighted for its TE-tolerance and even investigated for potential applications in remediation [49]. It plays a role in the degradation of polysaccharides, e.g., cellulose, xylan and chitin, OM turnover and in nutrient cycling [50]. Proteobacteria, the second most abundant phylum, is known to contain a wide range of functional bacteria, including fast-growing bacteria and many plant promoters, as well as growth-promoting bacteria [51]. Furthermore, Proteobacteria taxa have been described to be more abundant in heavily TEpolluted soils [44,52]. The latter phylum, Gemmatimonadota, is widespread in ecosystems around the world, making it an important member of soil bacterial communities [53]. Therefore, the long-term response of the bacterial community to TE results from the ability of these microorganisms to develop TE-tolerance without affecting overall community structure, as proposed by Brandt et al. [54].
The addition of amendments, such as steel slags or mycorrhizal inoculum, altered the structure of the rhizosphere bacterial community without affecting its metabolic potential, which was assessed by the AWCD parameter. In fact, the low dose of the mineral amendment increased bacterial biomass and caused a shift in the bacterial community structure. The Actinobacteria phylum was significantly less abundant under the slag-amended conditions, while the relative proportion of Gemmatimonadota increased. These results are in agreement with those of Gremion et al. [40] on a TE-polluted soil, but also with the functional redundancy of these communities [55]. Mineral granulometry has been shown to act as a physical support for the fixation of microorganisms (i.e., bacteria, fungi, etc.) and plants, as well as a nutrient reserve, and that the abundance of Acidobacteria decreases with increasing pH [56]. In this respect, and due to the alkaline nature of slags, the pH and mineralosphere can be important factors influencing microbial communities [56,57].
Cultivation of coriander also resulted in a change in the bacterial biomass and community structure. A significant increase in Gram-negative bacterial biomass under vegetated conditions and in the proportion of Actinobacteria phylum was detected. It is well known that plants secrete specific root exudates to shape the microbial community structures in the rhizosphere. These compounds, including sugars, amino acids and organic acids, have various effects on bacterial community composition and microbiome assembly [58]. In addition, aromatic plants are also known to be able to secrete secondary metabolites in their root exudates [59][60][61]. Due to their anti-microbial properties, some of these metabolites probably influence bacterial richness and diversity by targeting specific communities [62]. Some compounds may therefore affect microbial communities in the rhizosphere [60]. This is especially true for mono-and sesquiterpenes found in aromatic plants [60].
We noted an increase in the genus Streptomyces under vegetated conditions, a genus that has been described for its abilities as plant growth-promoting bacteria, but also for TE biosorption, making it particularly interesting in a phytostabilization context [63,64]. Moreover, the relative abundance of some genera (Sphingomonas, Gemmatimonas, Kribbella, . . . ) was very low, but they could play an important role in bio-geochemical cycles [65,66]. In fact, rare taxa are more modulated in abundance than dominant taxa [66], and therefore could be important factors in the resilience and adaptability of soil microorganisms [66]. Interestingly, some bacterial phyla have shown the ability to limit the accumulation of TE in plants. Indeed, some Actinobacteria and Proteobacteria can alter the bioavailability of heavy metals through redox, precipitation or sorption and desorption reactions [49]. These phyla, such as Firmicutes and Bacteroidetes, are known to possess TE-tolerance genes encoding proteins involved in efflux and sequestration of TE ions, such as AT-Pases, D-protein or B-protein (copper tolerance and transport, respectively) and oxidative stress attenuation proteins, such as superoxide dismutase, alkyl hyperoxide reductase or mycothiol reductase [67,68].
Therefore, coriander in TE-polluted soil activated the resilience of bacterial communities. This resilience of microbial communities and functions in TE-polluted soils is usually related to a shift from sensitive to tolerant species, genetic modifications leading to TEresistance, a transfer of genes encoding resistance or tolerance to TE, or a decrease in TE-bioavailability [69]. The recruitment of bacterial species belonging to the genera Pseudarthrobacter, Arthrobacter and Streptomyces by root exudates of coriander to the detriment of Gemmatimonadota bacteria, was associated with an increase in soil functionality that corresponds to such a shift. In this case, species belonging to Gemmatimonadota could be considered as indicator species for resistance to TE. Similarly, species belonging to Pseudarthrobacter, Arthrobacter and Streptomyces could be more sensitive to TE, but can characterize the resistance of TE-polluted soils when suitable substrates are available, such as those provided by root exudates in phytomanagement.

In Situ Experimental Design
Study Area. The study took place on an agricultural soil, historically polluted by TE (50 • 25 55.5 N, 3 • 02 25.5 E; elevation 23 m) in the "Hauts de France" region (northern France). The experimental site was located 600 m away from a former Pb and Zn smelter, Metaleurop Nord, whose activity generated significant amounts of TE-rich dust, resulting in TE-pollution of the topsoil (0-30 cm) [70,71]. The topsoil is characterized as silt loam, with a slightly alkaline water pH (7.9). It contains high TE total concentrations of Cd, Pb and Zn (7, 394 and 443 ppm, respectively), of which the NH 4 NO 3 extractable fractions account for 0.114, 0.078, and 0.930 ppm, respectively [6]. These concentrations were approximately 17, 11 and 6-fold higher, respectively, than the regional background values for agricultural soils [72].

Experimental Design
The in situ experimental design was divided into 20 equal plots of 9 m 2 , to investigate the effects of one biological (AMF) and two mineral amendments on soil function, with or without planting an annual herbaceous aromatic plant species, coriander (Coriandrum sativum L.). The mineral amendments, ladle and oven slags, were by-products of steel production, often used in agriculture due to their buffering capacity [73]. Their main constituent was CaO and Fe 2 O 3 , respectively. The biological amendment was the commercial mycorrhizal inoculum SYMBIVIT ® (INOCULUMplus, Bretenière, France), which contained six AMF species: Rhizophagus irregularis, Funneliformis mosseae, Glomus microaggregatum, Glomus claroideum, Claroideoglomus etunicatum and Funneliformis geosporum. The surface of each test area was 9 m 2 (3 m × 3 m), and each test area was randomly arranged and separated from each other by a 1 m wide corridor. A total of 20 conditions were designed (Table S1). Before sowing, the soil of the test plots was amended or not with 2 different concentrations of 1.5 and 45 t ha −1 mineral amendment (ladle and oven slags) and/or the biological amendment at a rate of 10 kg ha −1 . These plots were then sown (7.5 kg ha −1 ) or not (unvegetated condition) with coriander (seeds were provided by Iteipmai, France) in March 2019. The unvegetated plots were kept in this condition by regular manual weeding.

Sample Collection
At the beginning of the in situ experiment, 5 soil samples representing the original condition were taken with a soil auger in the 0-20 cm soil horizon, and stored at −20 • C until analysis. Three weeks after sowing, germination rate was assessed by counting the number of coriander plants in 3 random squares (1 m 2 ) for each vegetated plot. After 11 weeks of cultivation, 5 new soil samples were collected in each plot from the 0-20 cm soil horizon and stored at −20 • C until analysis; the aboveground biomass and roots of 5 plants were also harvested from the vegetated area. The height of the aerial parts was measured first, before the plants were dried and weighed, while the fresh roots were washed and stored at 4 • C until stained for AMF root colonization determination. The roots were cleared in 10% (w/v) KOH and stained with 0.05% (w/v) trypan blue (Alfa Aesar) [74] to determine root colonization [75]. Briefly, root fragments were placed on microscopic slides and mycorrhizal structures (arbuscules and vesicles) were counted for each fragment observed. Three counts per fragment were performed. The percentage of total colonization was equal to the number of intersections with mycorrhizal roots × 100/number of total intersections. This percentage of colonization was expressed as the percentage of arbuscules, vesicles and total colonization.

Fatty Acid Analysis
After removing plant debris from the soil samples, 3 g of freeze-dried soil (3 samples per test area) were used for the analysis of fatty acid content. Lipid extraction was performed according to Frostegård et al., with a mixture of chloroform:methanol:citrate buffer-0.15M, pH 4.0-(1:2:0.8 v/v/v), under agitation for 2 h [76]. The lipid material was fractionated on Solid Phase Extraction (SPE) columns containing silica (6 mL volume, 500 mg sorbent, Interchim, Montluçon, France) by successive elutions with chloroform, acetone and methanol (1:2:1, v/v/v) into neutral lipids (NLFA), glycolipids and polar lipids containing phospholipid (PLFA) [76]. The NLFA and PLFA were trans-esterified with 0.2 M KOH in methanol to obtain free fatty acid methyl esters. The resulting fatty acid methyl esters were analyzed using a gas phase chromatography-mass spectrometer (GC-MS) Shimadzu QP-2010 Ultra (Shimadzu, Japan) equipped with a single quadrupole mass detector and simultaneously coupled with a flame ionization detector (FID). Samples were analyzed in split mode (80:1 ratio) on a ZB-1MS fast capillary column (10 m length × 0.1 mm inner diameter × 0.1 µm phase thickness, 100% dimethylpolysiloxane, Zebron, Phenomenex, Torrance Calif, CA, USA) using helium as the carrier gas at a constant linear velocity (40 cm s −1 ). The injector temperature was 280 • C, and the detector temperatures were 330 • C and 280 • C for FID and for the ion source, respectively. The temperature program started with an initial temperature of 175 • C and increased by 25 • C every minute to reach a final temperature of 275 • C, which was maintained for 0.5 min. The ionization mode was electronic impact at 70 eV and the mass range between 50.0 and 400.0 u was scanned. The single impact monitoring (SIM) mode was used simultaneously. Quantification of fatty acids was performed using nonadecanoic acid methyl ester (C19:0, Sigma-Aldrich) as an internal standard. Fatty acids were identified by comparing their relative retention time with that of commercial fatty acid methyl ester standards (C4-C24:1, Sigma-Aldrich), and by comparing them with spectra either obtained from commercial standards and/or published in the literature (NIST Standard Reference Database). Phospholipid fatty acid (PLFA) i15:0, a15:0, i16:0, i17:0, a17:0 and PLFA cy17:0, C18:1ω7 and cy19:0 were used as indicators of Gram-positive and negative bacterial biomasses, respectively [77].

DNA Extraction
Genomic DNA was extracted in triplicate directly from 250 mg fresh soil (n = 60) using the NucleoSpin Soil ® kit (Macherey-Nagel, Düren, Germany) according to the manufacturer's instructions. The quality of the extracted DNA was checked using 1% (w/v) agarose gels. Quantification of the extracted DNA was performed using a SpectraMax ® iD3 spectrophotometer (Molecular Devices LLC, Sunnyvale, CA, United States). The extracted DNA was stored at −20 • C until use.
Polymerase Chain Reaction and sequencing of bacterial 16S rRNA gene Hypervariable regions V3-V4 of the 16S rRNA gene were amplified from 1 ng of genomic DNA using a PCR thermal cycler (Agilent Surecycler 8800) with the forward primer CS1_341_F (ACACTGACGACATGGTTCTACACCTACGGGNGGCWGCAG) and the reverse primer CS2_805_R (TACGGTAGCAGAGACTTGGTCTCTGACTACCAGGGTATCTAATC) [78]. Two independent PCR reactions were performed per DNA sample.

Illumina MiSeq Sequencing
All PCR products were pooled (n = 60) and sent to the Genome Quebec Innovation Centre (Montreal, QC, Canada) for sequencing, using Illumina MiSeq platform, producing paired-end 2 bp × 300 bp.

Data Analysis-Bioinformatic Processing
Sequences were loaded into the Galaxy instance (v.2.3.0) of the Genotool bioinformatics platform (http://sigenae-workbench.toulouse.inra.fr, accessed on 27 December 2022) to be processed in the pipeline FROGS (Find Rapidly OTU with Galaxy Solution, [79]). Sequences with ambiguous bases (N) or without the specific primers were removed. The Cutadapt software [80] was used to search and cut primer sequences with less than 10% difference. The number of sequences identified was 2,660,014 and sequence clustering was performed using the SWARM algorithm (v2.1.5; [81]. A first denoising step was performed to build very fine clusters with minimal differences (d = 1), and a second one with an aggregation distance of 3. The resulting representative sequences for each cluster, or OTUs (Operational Taxonomic Units), were subjected to chimera detection using the VSEARCH algorithm [82] and eliminated, as well as those with an abundance <0.001% of the total abundance. In the end, 48,251 OTUs were retained, corresponding to 2,264,807 sequences. The taxonomic classification of each OTU against the Silva database (v138.1) was done using RDPClassifier [83]. The abundance of the different samples ranged from 65,070 to 16,243. As a result, a rarefaction of OTUs at species level was performed with 16,243 as the rarefaction value ("rrarefy" function of the vegan package). Following this, diversity indices were determined. Comparisons of means were made ( Figures Shannon, and Chao1). As the Shapiro-Wilk group tests for the different indices indicated significant p-values (samples not following a normal distribution), comparisons via Wilcoxon-Mann-Whitney and Kruskal-Wallis tests were performed.

Nucleotide Sequence Accession Numbers
The 16S rRNA gene sequences of the whole dataset have been deposited in NCBI Sequence Read Archive (SRA) database and can be found under accession number PRJNA918629.

Soil Bacterial Metabolic Profiles
The metabolic potential of soil communities was determined from community-level physiological profiles (CLPPs) using Biolog EcoPlates (Biolog Inc.; [55]). For each soil sample, an EcoPlate containing 31 different carbon sources (and one blank with no carbon source) in triplicate was inoculated, and incubated for 196 h at 25 • C in an OmniLog ® system (Biolog Inc.). The rate of carbon source utilization was recorded by the reduction of tetrazolium, a redox color indicator that changes from colorless to purple and is detected at a wavelength of 590 nm. Data were recorded every 15 min of incubation and saved in OmniLog ® units generated by Biolog Data Analysis software (v1.7). Values for each well were calculated by subtracting the blank values from each well of the plate. A single absorbance time point at 50 h was used for the comparisons, as recommended in [84]. The metabolic potential of the soil microbial communities in each sample, expressed as the average well color development (AWCD), was calculated at the determined time point by dividing the sum of the optical density data by 31 (number of substrates). The total number of wells in a replicate with an absorbance above 25 OmniLog ® units was counted to determine the functional richness of the soil microbial community.

Statistical Analysis
Statistical analyses were performed using XLSTAT software (v2021.1.1; Addinsoft, Paris, France). Normality of the data was first subjected to the Shapiro-Wilk test. Depending on the result of this test, parametric (ANOVA) or non-parametric (Kruskal-Wallis) tests were used, followed by Fisher's post-hoc test (Least Significant Difference, LSD). Data expressed as percentages (AMF colonization rate) were converted to arcsine values (ASIN function in Microsoft ® Excel v16.50) before statistical analysis. Bacterial community analysis was performed using the R environment (v4.2.2; http://www.r-project.org/, accessed on 27 December 2022). Permutational multivariate analysis of variance (PERMANOVA), based on the Bray-Curtis dissimilarity, were performed with 1000 permutations using the "adonis" function of the vegan [85] R package. Bacterial richness (Chao1) and diversity (Shannon) indexes were calculated using the Phyloseq (v1.36.0; [86]) R package.

Conclusions
This study investigated the effect of coriander cultivation and introduction of AMF in combination with mineral amendments (steel slags) on the bacterial function of the rhizosphere, an aspect that is poorly understood and still poorly studied. The introduction of amendments, such as steel slags or mycorrhizal inoculum, into the soil had no significant effect on coriander growth. However, these amendments altered the structure of bacterial communities in the rhizosphere, without affecting their metabolic potential. In fact, Actinobacteria were significantly less abundant under the slag-amended conditions (55% on average), while the relative proportion of Gemmatimonadota increased (10% on average). In contrast, soil vegetation by coriander increased Gram-negative bacterial biomass (mean 2.96 to 3.87 µg g −1 of DW soil) and the proportion of Actinobacteria phylum (62%). Thus, the vegetation improved soil quality and health by increasing the diversity of bacterial metabolic functions and diversity. Our results showed that coriander activated the resilience of bacterial communities in TE-polluted soils. Consequently, our results suggest that steel slag, which contributes to the circular economy through its recycling, in combination with soil vegetation that includes aromatic plants such as coriander, potentially provides long-term ecological benefits to TE-polluted soils and should be considered as a reliable tool for future TE-polluted soil remediation. Further studies are needed to confirm the contribution of steel slag to assisted TE-phytostabilization and its persistence in soil.