1. Introduction
The extensive application of herbicides is a result of constant advancements in the studies of application technology and the expansion of cultivation areas, mainly for food, fiber, cellulose, and bioenergy production [
1,
2]. Several researchers have dedicated themselves to studying the impact of herbicides on the environment, obtaining differentiated responses regarding toxicity to non-target organism groups and the real effect on the sustainability of productive systems [
3,
4]. Reports of adverse effects on organisms and the functioning of soil ecosystems are observed [
5], as well as persistence in soil and areas with accumulations of water, modulating interactions among non-target organisms, which is also considered negative [
6].
The progress observed in recent decades regarding the development of safer molecules and the reduction of product losses to the environment through innovation in herbicide application techniques is undeniable [
7].
Approaches to the impact of herbicide use have strongly considered the negative impact of molecules on soil microorganisms [
7]. However, a very small portion of research discusses with authority the concentrations used, as well as the final effect on soil microbial dynamics [
8,
9].
Considering tropical agriculture during grain cultivation, we observe the use of two groups of herbicides, one applied pre-emergence and the other post-emergence of plants. These herbicide groups complement each other in weed control, as while one acts on the seed bank and prolongs the effect over time, the other acts directly on treated plants, controlling present vegetation and ensuring the formation of organic material through plant residue [
10,
11]. The mixture of metolachlor and glyphosate herbicides has these working effects, and for this reason, it was chosen for this research.
Metolachlor(2-chloro-N-(2-ethyl-6-methylphenyl)-N-[(1S)-2-methoxy-1-methyl] acetamide) is a herbicide belonging to the chloroacetamide class. In the United States, metolachlor is among the three herbicides with the highest volume of use [
12]. It has prolonged persistence in the soil, commonly used in crops such as
Eucalyptus sp.,
Zea mays L.,
Glycine max (L),
Solanum tuberosum L.,
Beta vulgaris L.,
Helianthus annuus L.,
Solanum lycopersicum L., and various others, as well as on lawns.
In several agriculturally important countries, metolachlor has been particularly effective against
Amaranthus sp.
Amaranthus palmeri was well controlled in corn crops in the United States with metolachlor [
13]. In Brazil, metolachlor has been one of the main herbicides applied pre-emergence, mainly for grass and some broadleaf weed control [
14,
15]. In Brazil,
Amaranthus palmeri is the main pest listed by the Ministry of Agriculture [
16], and metolachlor is one of the control options. The applied doses of metolachlor average 960 g/ha for weed control in grain legume crops [
17]. A study indicates that the safe application of metolachlor at a rate of 480 g/ha can be applied approximately 20 days before crop sowing to control
Lolium rigidum in wheat selectively [
18]. Thus, the dosage of metolachlor may vary depending on factors such as the specific crop, targeted weeds, soil clay content, organic matter, and local conditions and regulations [
19]. The herbicide’s persistence, which refers to its control power for a longer period, is highly dependent on climate and soil conditions [
20]. Under laboratory conditions, metolachlor has higher persistence [
21,
22]. Thus, research should consider field-observed results as being more promising.
An innovative approach advocated in this research is the study of the impact of metolachlor on soil microbial dynamics. Research on the impact of herbicides on microorganisms is constantly conducted, but innovations in integrated weed management have shifted the perspective of the hypotheses tested. Phytoremediation green technology is a good example of research where the potential for microbial degradation is associated with remedial species tolerant to herbicide residues [
23]. The adoption of phytoremediators in modern agriculture has advanced significantly, with a strong trend towards increasing use due to the need for better utilization of tropical agricultural areas, with incentives for efficient crop rotation [
15,
24]. In this sense, it is necessary to know species adapted to summer or winter to protect fields between main crops, covering the soil and competing with weeds while reducing residues from herbicides previously applied pre-emergence [
25,
26]. Therefore, phytoremediation is an innovation desired by farmers who wish to maintain the positive effect of residual herbicides while eliminating the impact on sensitive crops, known as the carryover effect [
24,
27].
A highly debated issue regarding the future use of herbicides is the impact of these molecules on agricultural sustainability. There is a contradiction between implementing phytoremediation technologies to allow for the rational use of pre-emergence herbicides and the expectation of greater preservation of microbial diversity, as herbicides are considered harmful to microorganisms. To assist in studies of this complexity, genetic studies have been proposed to comprehensively evaluate the real impact of herbicides on soil health, especially in intensive agriculture systems such as tropical ones [
24,
28].
Studies on metagenomics have been used to identify taxonomic and functional groups of microorganisms, genes, and metabolic pathways associated with xenobiotic manipulation, which typically involve an assessment of the microbiota and the identification of markers indicating the presence of individuals, genes, and enzymes involved in the manipulation process [
23,
29]. Herbicide application in soil can result in various effects on different components of the soil microbial community, with both beneficial and adverse consequences. The nature of these effects, whether direct or indirect, determines the overall impact on the microbial community [
30]. The consequences involve herbicides being detrimental to microorganisms, leading to a decrease in microbial biomass, subsequently affecting soil heterotrophic respiration, organic matter decomposition, and the activity of microbes involved in the nutrient cycle [
31,
32]. On the other hand, herbicide addition can have a direct positive effect on soil microbes, providing a resource that supports their growth, for example, by inhibiting some less important groups in soil resilience [
33]. Additionally, the impact of herbicides on plants can indirectly influence microorganisms [
34].
We hypothesize that the phytoremediation process, considered a green technology for tropical agriculture, may be key in the relationship between herbicide residues and soil health as measured by functional genomics of the plant rhizosphere-associated microbiota. Thus, the objective of this study was to analyze the impact of using metolachlor as a case study of herbicides and phytoremediating plants on soil biodiversity and resilience.
2. Materials and Methods
The experiment was conducted in the field at the experimental area of the Faculty of Agricultural Sciences of the Federal University of the Valleys of Jequitinhonha and Mucuri—UFVJM, located at the JK campus, Diamantina, MG. The geographic coordinates of the area are latitude 18°10′ S and longitude 43°30′ W, with an altitude of 1388 m. The region has a mesothermal climate, classified as a Subtropical high-altitude—Cwb according to the Köppen classification, characterized by mild and humid summers and cold and dry winters. Before the experiment, the area was desiccated with glyphosate (1440 g/ha) and cultivated with corn and beans without the application of pre-emergence herbicides.
The soil where the experiment was installed is classified as Orthic Quartzarenic Neosol, with a sandy loam texture. Its chemical characteristics include a pH (H
2O) of 5.20, 2.64 cmolc dm
−3 of H + Al
+, 1.19 cmolc dm
−3 of Ca
2+, 0.15 cmolc dm
−3 of Mg
2+, 8.10 mg dm
−3 of P, 36.0 mg dm
−3 of K, and 1.60 dag kg
−1 of organic matter [
35]. The experiment was conducted in a randomized complete block design with three replications.
Ten treatments were planned to combine doses of metolachlor (in the ready mix with glyphosate, called Sequence
®) and two green manure species, selected according to their potential for remediating areas with herbicides:
A. sativa (cultivate: URS-Taura/BRSEE-Ds
®) [
36] and
M. sativa (cultivate: URS/BRSEE-Ds
®) [
37]. The list of treatments is presented in
Table 1.
After soil preparation and fertilization, commercial recommendations for green manures were followed [
38,
39]. The commercial mixture Sequence
® was applied, composed of metolachlor (353.8 g/L) and glyphosate (265.7 g/L). For this research, we did not consider the effect of glyphosate due to its strong soil adsorption [
40]. So, the treatments were conditioned to the residual effect of metolachlor. We opted for metolachlor doses, equivalent to half (1.5 L/ha) and the full average (3.0 L/ha) recommended dose for major agricultural and forestry crops. For application, a backpack sprayer equipped with a transversal bar with four nozzles at a 110° angle and a flow rate of 0.15 gallons per minute, equivalent to 150.0 L ha
−1 of solution, was used. Three days after application, the predetermined plots were sown with the species
A. sativa and
M. sativa. The planting density was 7.0 g and 2.0 g seeds/m
2, respectively, for
A. sativa and
M. sativa [
36,
37].
After 60 days of herbicide application, three plants were sampled per plot, and the roots were removed to collect the rhizosphere soil. These samples were collected and kept at −80 °C for metagenomic analysis. Samples were taken only from plots growing
A. sativa and
M. sativa, with and without herbicide, and from native vegetation soil with no previous herbicide contamination (No-0). Metagenomic DNA was extracted and purified. Sequencing was performed on the Illumina Hiseq 2000 platform (Illumina, San Diego, CA, USA) with amplification of the 16S rRNA and ITS genes. Data pre-filtering was performed using the Readfq program v8
https://github.com/cjfields/readfq (accessed on 22 January 2024) to remove sequences containing low-quality bases and those that did not reach a proportion greater than 10.0%. Subsequently, high-quality sequences from each sample were assembled (MEGAHIT), generating “Scaffolds” with fragments larger than or equal to 500 bp. The sequences were compared to the genomic database of bacteria, fungi, and archaea using the Blastx tool. Genetic sequences were processed and compared with the database of fungi and bacteria to identify and quantify the species present in the samples. The relative abundance values (%) of microorganisms were determined by the proportion of the microorganism among all the groups to which it belongs. The symbol (*) represents the microbial group below the top 10 occurrences of microorganisms at the family and genus levels.
The domain was estimated using the Gini-Simpson indices [
41], and the richness was quantified by the number of different species present in the sample, according to Chapo1. The species distribution quantified the diversity based on the number of species present in the samples, and for this, we used the Shannon Index [
42]. Additionally, to indicate the distribution of the different species in the sample, we estimated using the Pielou index, determining the uniformity [
43]. The diversity of soil microorganisms was estimated as a percentage based on dominance, richness, species distribution, and uniformity between species indices. To determine correlations between rhizosphere microbiota, phytoremediating plants, and herbicide doses, principal component analysis, Pearson correlation network in cords, and neural networks were performed using R Software v4.4.1.
Principal component analysis (ACP) of microbial biodiversity was performed based on the microbial phyla associated with phytoremediating plants to determine changes in the studied samples. The separation of groups and subgroups was determined by the significance cutoff point (above 0.9 or below −0.9). The distribution of microbial phyla for the plants A. sativa and M. sativa was determined with the Pearson correlation network (above 0.6 or below −0.6) and grouped by the string procedure established according to the groups formed by the ACP. The distribution of microbial families associated with A. sativa and M. sativa plants in soil phytoremediation areas containing metolachlor residues was estimated with Pearson’s correlation with a significance cutoff point (above 0.6 or below −0.6), represented graphically by the neural network.
After 80 days of herbicide application, soil samples were collected to determine herbicide residue in the soil by chromatography. Samples were collected from plots with herbicide application, i.e., plots with remediation species and uncultivated areas that received herbicide application and were weeded throughout the period. For the analytical determination of residues, extraction was performed by the modified QuEChERS method [
44], and analysis was carried out by ultra-high-performance liquid chromatography coupled with tandem mass spectrometry (UHPLC-MS/MS). The extracts were cleaned using dispersive solid-phase extraction (d-SPE). For the d-SPE phase, magnesium sulfate and the sorbents octadecylsilane (C18) and primary secondary amine (PSA) were used. The extracts were then shaken, centrifuged, and filtered. Before analysis by UHPLC-MS/MS, the samples were diluted five times in ultrapure water. The samples were analyzed on the Waters (USA) UHPLC-MS/MS system, equipped with a liquid chromatograph; triple quadrupole MS detector, model Xevo TQ; electrofoil ionization source/interface; nitrogen peak generator; solvent control system (binary pump system) for high-pressure gradient operation; analytical column Acquity UPLC
® BEH C18 (50 × 2.1 mm, 1.7 μm) from Waters (USA); and data acquisition system with MassLynx software v4.1 (Waters, Milford, MA, USA). Selected reaction monitoring was used for quantification and identification of substances to be analyzed. The mobile phase used was (A) water:methanol (98:2,
v/
v) and (B) methanol, both containing 5 mmol L
−1 of ammonium format and 0.1% formic acid (
v/
v), with a flow rate of 0.225 mL min
−1 and an injection volume of 10 µL [
45]. The manipulation estimate was determined by the average value of each area for a qualitative analysis. The values presented in the figures are given as means ± standard errors.
3. Results
The remediation species Avena sativa and Medicago sativa exhibited satisfactory growth in the areas, regardless of metolachlor residues, confirming their tolerance to this herbicide.
Soil remediation containing metolachlor was observed in areas managed with green manures
A. sativa (AM1 and AM2) and
M. sativa (MM1 and MM2) for 80 days, where the residual concentration in the soil depends on the remediation species (
Figure 1A,B). Remediation efficiency was particularly increased in areas AM2 and MM2, where the herbicide application dose was 1061.40 g L of metolachlor (
Figure 1A,B). Residual concentrations of metolachlor were 45.5% (NoM1) and 75.0% (NoM2), at 0.012 mg kg
−1 and 0.022 mg kg
−1, respectively, dissipating 25.0% (AM1) and 54.5% (AM2) of the molecule in the presence of
A. sativa (
Figure 1A). In the presence of
M. sativa, there was a dissipation of 7.7% (MM1) and 36.4% (MM2) of the soil molecule (
Figure 1B). The dynamics of soil microbiota in association with phytoremediating plants were diverse and selective, with the potential for herbicide molecule biodegradation at high concentrations in areas AM2 and MM2.
In the control soil sample (from Cerrado) where there was a concentration of 1061.40 g/ha of metolachlor and persistence of 0.012 (AM1), 0.014 (AM2), 0.009 (MM1), and 0.0010 (MM2) mg/kg
−1 of the molecule, there was high microbial biodiversity (
Table 2).
Soil samples A0 and AM1 showed lower species distribution uniformity according to the Pielou index, reaching −17.44% compared to areas with medium and high herbicide concentrations (
Table 2). For
A. sativa, the richness of the associated microbial species was higher in AM1 and AM2, with 116–121 species. The distribution of microbial species in area A0 was 19.51% lower than in AM2 and No-0. In
M. sativa, the highest richness was observed in M0 and MM1, ranging from 5 to 13 compared to MM2 and No-0.
M. sativa exhibited higher values of evenness, richness, distribution, and dominance compared to
A. sativa (
Table 2). In the MM1 areas, there was a 10.6% greater richness compared to the No-0 area. In the MM2 treatment, richness decreased by 4.1% compared to the MM1 area. The distribution of species in the native Cerrado area was 4.8% lower than in the rhizospheric environment of
M. sativa, regardless of metolachlor presence (
Table 2). For
M. sativa, the response to dominance was the same, regardless of the treatment. Overall, the values for the global ecological index indicate the following decreasing sequence for biodiversity: MM1 (80.0%) > M0 (79.1%) > No-0 (62.83%) > A0 (61.9%) > MM2 (61.06%) > AM2 (60.17%) > AM1 (59.29%).
The principal component analysis of soil phytoremediation with
A. sativa and
M. sativa highlighted the relationships between the abundance value of microbial phyla concerning herbicide doses and uncultivated Cerrado soil, forming four phylogenetic groups (
Figure 2A,B). Regarding
A. sativa, native Cerrado soil samples (No-0) formed cluster I, with three subgroups. Subgroup I is composed of Acidobacteriota, Chloroflexota, Chytridiomycota, Dormibacterota, Eremiobacterota, and Mucoromycota (
Figure 2A). Subgroup II is formed by Verrucomicrobiota, and subgroup III by Ascomycota (
Figure 2A). Cluster II, related to the environment with higher metolachlor dose, is formed by Planctomycetota, Proteobacteria, and Mortierellomycota. Cluster III, related to the lower metolachlor dose, is formed by Basidiomycota, and cluster IV,
A. sativa without herbicide, is formed by Rozellomycota, Gemmatimonadota, and Actinobacteria. For
M. sativa, cluster II, also related to the higher herbicide dose, was formed by Actinobacteria, Proteobacteria, and Rozellomycota. Cluster III,
M. sativa without herbicide, was formed solely by Chytridiomycota. Cluster IV, regarding the lower herbicide dose, was formed by Basidiomycota, Gemmatimonadota, Planctomycetota, and Mortierellomycota (
Figure 2B). The percentage of the interaction of microbial phyla in the rhizospheric soil of
A. sativa was 28.8%, highly positive (>0.6) among the 15.0 observed phyla (
Figure 3A). Null interactions (0.0), moderately correlated (between 0.3 and 0.5), and highly negative correlations (<0.6) were, respectively, 0.5%, 56.6%, and 14.2%, highlighting the dynamism among organisms (
Figure 3A). Positive correlations between microbial phyla with
M. sativa were 1.5 times lower compared to
A. sativa (
Figure 3B). In rhizospheric soils of
M. sativa, microorganisms of the phyla Mortierellomycota x Dormibacterota, Gemmatimonadota x Actinobacteria, and Ascomycota x Rozellomycota show null correlations (
Figure 3B).
A total of 94 families of fungi and soil bacteria were found to be associated with
A. sativa and
M. sativa in the presence and absence of the herbicide molecule:
Acetobacteraceae,
Acidimicrobiales,
Acidobacteriaceae,
Acidobacteriales,
Alphaproteobacteria,
Aspergillaceae,
Azospirillaceae,
Baltobacteraceae,
Bionectriaceae,
Bryobacteraceae,
Bulleribasidiaceae,
Burkholderiaceae,
Capnodiales,
Caulobacteraceae,
Chaetomiaceae,
Chaetosphaeriaceae,
Chaetothyriaceae,
Chaetothyriales,
Chrysozymaceae,
Chthoniobacteraceae,
Chthoniobacterales,
Cladosporiaceae,
Coniochaetaceae,
Coniothyriaceae,
Cucurbitariaceae,
Cunninghamellaceae,
Cyanobacteriia,
Cyphellophoraceae,
Deinococcaceae,
Dermatophilaceae,
Didymellaceae,
Didymosphaeriaceae,
Dormibacteraceae,
Dormibacteria,
Enterobacteriaceae,
Filobasidiaceae,
Gaiellaceae,
Gemmataceae,
Gemmatimonadaceae,
Geodermatophilaceae,
Helotiaceae,
Herpotrichiellaceae,
Hydnodontaceae,
Hypocreaceae,
Isosphaeraceae,
Jatrophihabitantaceae,
Ktedonobacteraceae,
Labraceae,
Lasiosphaeriaceae,
Limnocylindrales,
Lycoperdaceae,
Magnaporthaceae,
Microbacteriaceae,
Micrococcaceae,
Micromonosporaceae,
Mortierellaceae,
Mucoraceae,
Mycobacteriaceae,
Mycosphaerellaceae,
Nectriaceae,
Niessliaceae,
Nocardioidaceae,
Ophiocordycipitaceae,
Phaeosphaeriaceae,
Plectosphaerellaceae,
Pleosporaceae,
Pleosporales,
Powellomycetaceae,
Rhizobiaceae,
Rhynchogastremataceae,
Sclerotiniaceae,
Solirubrobacteraceae,
Solirubrobacterales,
Sordariales,
Sordariomycetes,
Sphingomonadaceae,
Spizellomycetales,
Sporidiobolaceae,
Stachybotryaceae,
Stellaceae,
Steroidobacteraceae,
Streptomycetaceae,
Streptosporangiaceae,
Sympoventuriaceae,
Teichosporaceae,
Tepidisphaeraceae,
Thermoleophilia,
Thyridariaceae,
Trichocomaceae,
Trichomeriaceae,
Trimorphomycetaceae,
Tumebacillaceae,
Vicinamibacterales, and
Xanthobacteraceae (
Figure 4A,B).
The families of microorganisms responded satisfactorily to the biotic and abiotic pressure imposed by soil agricultural disturbance, even with the application of metolachlor and the cultivation of green manures (
Figure 4).
The uniformity of bacterial species distribution in the rhizosphere of
A. sativa was different from the other treatments for this species (
Table 3). In areas cultivated with
A. sativa containing metolachlor residues, all ecological indices showed reductions—29.0%, 32.0%, 33.0%, and 24.0%, respectively, for uniformity, richness, distribution, and dominance (
Table 3). This effect was not observed in soil samples from the area with
M. sativa cultivation (
Table 3).
The abundance of bacteria varied among treatments, with the most prevalent families presented in
Table 4. Families with abundance greater than 10.0% were
Burkholderiaceae (10.63% for A0),
Micrococcaceae (11.01% for A0, 10.2% for AM1, 14.27% for AM2, and 36.6% for M0, 17.5% for MM1), and
Xanthobacteraceae (12.56% for AM2, and 17.74% for MM2). The observed
Bryobacteraceae family was less abundant in areas with
M. sativa. Similarly,
Ktedonobacteraceae was not observed in areas with
A. sativa (
Table 4).
Despite the low percentage of undetermined families (typically new ones or those with DNA fragments insufficient for cataloging), the highest percentage of abundance for predominant bacterial genera was considered undetermined, ranging from 31.0% (for M0) to 47.12% (for A0). Also, a higher percentage of uncultivated bacterial genera were observed, composed of species with ecological importance but difficult to cultivate in vitro (
Table 4). The most predominant bacterial genera were
Bradyrhizobium sp. (4.78%),
Massilia sp. (2.06%) in AM2;
Sphingomicrobium sp. (2.69%),
Mycobacterium sp. (2.29%),
Arthrobacter sp. (1.2%),
Pseudarthrobacter sp. (0.8%) in AM1; and
Pseudarthrobacter sp. (3.47%),
Streptomyces sp. (2.59%),
Sphingomicrobium sp. (2.01%),
Baekduia sp. (0.89%) in A0.
The uniformity of fungal species distribution also varied to a greater degree in the rhizosphere of
A. sativa compared to
M. sativa (
Table 5). The values were lower for soil samples of
A. sativa without metolachlor and under the effect of the lower herbicide dose. A similar effect was observed for the distribution variable (
Table 5). Regarding
M. sativa, the values for uniformity, richness, distribution, and dominance were relatively consistent among treatments; however, they were higher than those observed for native Cerrado soil (
Table 5).
Among the abundance percentages of fungal families, the highest average values (above 10.0%) were observed for
Aspergillaceae (10.02%),
Didymellaceae (10.05%) in
A. sativa soils, and
Chaetomiaceae (10.27%) and
Nectriaceae (17.71%) in
M. sativa soils (
Table 6). Among these averages, the highest percentages observed were for
Aspergillaceae in AM2 (12.18%),
Didymellaceae in AM2 (13.19%),
Chaetomiaceae in MM2 (13.69%), and
Nectriaceae in M0 (28.88%) (
Table 6). In the assessment of each treatment, other high values were observed for
Mortierellaceae in AM1 (19.57%),
Aspergillaceae in M0 (10.49%),
Cladosporiaceae in M0 (10.06%%) and MM1 (14.61%),
Didymellaceae in M0 (13.84%), and, with the greatest discrepancy,
Rhynchogastremataceae in MM1 (20.0%) (
Table 6). Under native Cerrado soil, the highest percentage of abundance was observed for
Aspergillaceae (21.14%) and
Nectriaceae (18.60%). The percentages of the abundance of undetermined fungal families were higher for AM1 (17.35%), AM2 (18.13%), and Cerrado native soil (11.06%) (
Table 6).
No significant percentage of the families
Filobasidiaceae and
Rhynchogastremataceae was observed in soil samples of
A. sativa. On the other hand,
Mortierellaceae did not show significant abundance in soils with
M. sativa (
Table 6).
The genera
Fusarium sp. and
Penicillium sp. were observed with higher abundance in the rhizosphere of
A. sativa; however, their presence was significantly reduced in the presence of metolachlor (
Table 6). The genera
Cladosporium sp. and
Papiliotrema sp., also abundant in
A. sativa, showed higher activity and abundance under the effect of the lower herbicide dose; however, their occurrence decreased with the increase in dose. In
M. sativa, the occurrence of
Papiliotrema sp. and
Saitozyma sp. was not significant. On the other hand,
Fusarium had an abundance percentage of 11.94 in the rhizospheric soil of
M. sativa but decreased in the presence of the herbicide (
Table 6).
4. Discussion
Despite being tolerant to herbicides, the tested species showed divergent potential for remediation. The more effective soil exploration by grasses results in a higher rate of nutrient and molecule absorption [
43], which could explain the higher efficiency of
A. sativa. This species has an optimized and fasciculate root length in the soil, allowing for the exploration of a larger volume of particles with molecule sorption. Its tolerance to metolachlor and various other herbicides has already been reported [
46].
Due to its taproot system,
M. sativa has secondary roots that can intercept the molecule to be remediated in lower concentrations. However, the presence of exudates acting in the remediation process should be considered. Compounds such as alginate oligosaccharides are beneficial for alfalfa root growth [
47]. According to Ref. [
48], many herbicides have been formulated with different carriers, including alginate, which allows for greater sorption.
In both situations, plants can select communities of microorganisms to compose the rhizospheric system, and their interaction can have a direct impact on the extent of herbicide degradation [
49]. The interaction between microorganisms and plants remediating environments with herbicide residues is underexplored, but it can be synergistic and positive [
50].
The herbicide should be strictly applied according to the regimen that defines its optimal dosage to mitigate any potential disturbance in soil homeostasis. A dose of 1.200 g of metolachlor can contribute to the increase in soil dehydrogenase activity promoted by the microbiota [
51]. In overdoses ranging from 15 to 400 g per kg of soil, the dynamics and ecological state of the soil are modified [
52].
The high uniformity of species distribution in environments AM2, No-0, and AM1 demonstrates that the herbicide is not capable of interfering with nutritional consumption. Under controlled conditions, it was found that part of the microbial energy consumption can come from the content of the metolachlor molecule, which has 15 carbon atoms and one nitrogen atom [
53]. On the other hand, microbial specificity or restriction is more affected by the rhizospheric effect of
A. sativa than by the presence of the herbicide. The restrictive association, due to microbial specificity, has already been observed for
A. sativa [
54], and the herbicide seems to provide greater abundance to the rhizosphere.
The presence
of A. sativa modulates the richness of the microbial population. This result leads to the hypothesis that the presence of the herbicide in the soil stimulated the formation of a more diverse community (
Table 2). Practically speaking, our findings indicate soils that are more diverse and balanced, with a dominance between 0.92 and 0.97 according to the Gini-Simpson index. In other words, an equitable and resilient rhizospheric microbial community signals soil with ecological health in all treatments (
Table 2).
Bacterial and fungal richness increases during the early stages of
M. sativa cultivation [
55]. The presence of the herbicide appears to interfere with rhizosphere dynamics, promoting a greater number of species. Although this increase is less apparent at higher concentrations, it still surpasses the richness observed in the absence of the molecule.
The reduction in species distribution in native Cerrado soil indicates that opting for legume cultivation will promote a highly resilient, stable, and balanced microbial population, evidencing the conditioning of healthy soil with a wide variety of ecological niches, even with the likely interaction with herbicide residues, after 80 days of application (
Table 2).
Thus, regarding the integration of the ecological indices analyzed, we can infer that the presence of the legume (M. sativa) provides greater microbial biodiversity, surpassing the environment cultivated with A. sativa and the native Cerrado environment. However, under the effect of the dose of 1061.4 g/ha of metolachlor, this biodiversity equals the other treatments, indicating that recommendations for M. sativa cultivation may be interesting in areas where metolachlor is used if the practiced doses are low.
In our research, the principal component analysis highlights groups with greater similarity to the imposed treatment [
56] in the cultivation of
A. sativa and
M. sativa with the presence of metolachlor. Overall, the cultivation of these species groups of microbial phyla isolates a special grouping for native Cerrado soil. It has been observed that native tropical soils exhibit greater microbial diversity compared to systems that promote disturbance, such as those of intense agriculture [
57]. Our findings have already proven that cultivation with green manures such as
A. sativa and, especially,
M. sativa, can mitigate negative effects, maintaining good biodiversity of microbial species, even with low doses of metolachlor. For phylogenetic diversity, the formation of groups that demonstrate the effect of agricultural activity is clear, separating, to a greater extent, the native Cerrado environment. The variation in the abundance of rare species “depends on rare microbes” and is linked to changes in anthropized areas. For this reason, they are observed in environments with low agricultural activity [
58]. However, they are important because they drive key processes in geochemical cycles and are significant actors in the degradation of xenobiotics.
The behavior of common microorganisms remains consistent regardless of the cultivated species, with groups observed for the species without herbicide or with the lowest dose and a new group for samples under higher doses of metolachlor. Thus, the differentiated effect for areas with the highest dose of the herbicide is confirmed.
We verified through principal component analysis that phyla such as Proteobacteria are present in the environment with higher doses of metolachlor for both plant species. However, the phyla Planctomycetota, Mortierellomycota, and Basidiomycota are examples of groups that tolerate the herbicide at higher doses in the rhizospheric environment of A. sativa but not for M. sativa, indicating a remedial action resulting from the possible plant–microorganism interaction. Similarly, the phyla Rozellomycota and Actinobacteriota were significant in the rhizosphere of M. sativa under the effect of the highest herbicide dose; however, in A. sativa, this was only under the effect of the highest dose or in the absence of the herbicide. Our findings are relevant for selecting groups of microorganisms to produce bioinputs used in the decontamination of areas, also because the recommendation of the microbial group is highly dependent on the remediating plant species.
The distribution of phyla occurs depending on the phytoremediating plant. The herbicide possibly moderates it, as the phylum Rozellomycota was highly correlated with Ascomycota in the presence of
M. sativa, indicating synergy between the organisms (
Figure 3A). This synergy may explain the remediation potential of
M. sativa. Soils with abundant fungi have been reported as promising for the rapid degradation of metolachlor [
59].
Despite agricultural activities causing disturbances in the soil environment, our findings demonstrate that the cultivation of green manures improves soil organic matter and increases microbial richness, especially when considering distinct families such as grasses and legumes [
60]. These microbial communities, associated with phytoremediators, play a fundamental role in the degradation of organic compounds and contribute to the resilience of the ecosystem, even in agricultural systems with pre-emergent herbicide applications such as metolachlor.
The pressure exerted by agricultural cultivation and herbicide application has altered the interaction between organisms, where rare groups in lower abundance are essential in soil suppressiveness and homeostasis [
57]. The negative impact of metolachlor on microorganisms has been reported for in vitro experiments and at doses above those practiced in the field [
61,
62]. Our findings indicate that, at appropriate doses, the presence of metolachlor does not amplify the effects of green manure cultivation practices, and at the highest concentration, even though they alter some microbial responses, the interactions and microbiological quality are maintained.
The bacterial sensitivity can be partly explained by specific functional characteristics, a recurring situation when bacterial cells are active in natural communities [
63]. The removal of metolachlor by remediation promotes bacterial alteration in the soil [
64], a situation identified in areas with
A. sativa, and partially explains the observed values. On the other hand, the decrease in indices in areas with this plant species without herbicide indicates a likely restriction of this microbial group associated with the plant. The defense of
A. sativa against certain bacterial groups was confirmed by other studies with gene expression patterns [
65].
Regarding
M. sativa, it is likely that the positive interaction between the plant and the microbial community outweighs any potential effect of herbicide presence. Studies indicate that in the rhizosphere, between rhizobacteria and
M. sativa plants, there is an active interaction at the level of extracellular amino acid exchange [
66]. According to these authors, plant roots absorb the amino acids excreted by bacteria, and these bacteria absorb the amino acids excreted by
M. sativa roots.
There is complexity in working with uncultivable and undetermined individuals due to the difficulty in measuring their functional potential and biotechnological applications. However, in our findings, there is a correlation between uncultivable and undetermined organisms in areas containing metolachlor (AM2 and MM2). The occurrence of possible bacterial genera involved in the biodegradation of the molecule was noticeable, as confirmed for
Streptomyces lydicus [
67],
Pseudomonas sp. [
68].,
Nitrososphaera,
Nitrosospira,
Sulfuricella,
Gemmatimonas,
Candidatus sp.,
Burkholderia,
Bradyrhizobium,
Accumulibacter,
Entotheonella,
Cupriavidus,
Azoarcus,
Thiobacillus, and
Sideroxydans [
69].
The richness and diversity ratio of fungi increased in the presence of metolachlor, associated with green manures. The effects of the interaction between green manures and microorganisms, combined with climatic effects, are more relevant predictors in the final composition of the fungal community [
70].
Symbiotic fungal species that colonize the roots of green manures are common, with reports of herbicide degradation such as atrazine. The interaction between the roots of green manures and fungi such as
Aspergillus niger,
Phanerochaete sp., and
Pleurotus sp. enabled the addition of sugars that aided these microorganisms in breaking down the atrazine molecule [
71].
Regarding
A. sativa, the decrease in the abundance of unidentified microorganisms is associated with the selection pressure imposed by this species, reducing occurrence, and in a synergistic effect with the herbicide prevailing, however, in some families such as
Mortierellaceae. The increase in the abundance of
Mortierellaceae has been reported as a result of metolachlor bioremediation [
72]. It helps explain the potential of
A. sativa as a phytoremediator of this herbicide.
Some fungal specimens undergo broader impacts on the composition and structure of the community, both biotically through the suppressiveness of control agents such as
Papiliotrema sp.,
Humicola sp.,
Trichoderma sp., and
Epicoccum and abiotically through the presence of the herbicide molecule [
73]. The genera
Humicola sp.,
Trichoderma sp.,
Pleurophragmium sp., and
Saitozyma sp. were the most abundant in soils of
A. sativa under the effect of the highest dose of metolachlor, indicating possible plant–microorganism interaction in the phytoremediation process. These genera are involved in enzyme production and nutrient cycling directly in phytopathogen control [
74] and have been observed as agents in herbicide decomposition [
75]. Although we associate higher abundance values with degradation potential, we do not dismiss the fact that species with low percentages, such as
Trichoderma, are widely known for metolachlor degradation [
62]. In our study,
Trichoderma was highlighted in all treatments; however, metolachlor application significantly altered the abundance of this fungus, increasing it in
A. sativa roots and decreasing it in
M. sativa roots. Thus, this variation serves as a good example that the association between fungi and roots of remediating species should be further explored so that agricultural practices using bioinputs can be more accurately directed.