Abstract

In order to reduce the amounts of pesticides used, and thereby their associated risks, new generations of less environmentally dangerous molecules with lower weight are currently being used in the mixtures sprayed on crops. Few studies have been made, however, to analyse their impact on the soil, and more particularly on the microorganisms living in the soil which maintain the essential functions of this ecosystem. By taking a microcosmic approach, we were able to assess the impact of the maize herbicides “cocktail” Mesotrione and S-metolachlor on global soil microbial activity, biomass, and structures, by using the formulated compounds, respectively, Callisto and Dual Gold (both registered brands of Syngenta). Our results highlighted a synergetic effect in “cocktail” microcosms resulting in an increase in the Mesotrione herbicide dissipation time and in an impact on the microbial community at onefold field rate equally to more than a single herbicide used at tenfold field rate.

1. Introduction

For several decades herbicides have been used in spray form to reduce unwanted vegetation in crops. However, the main indirect casualty of these pesticides, following their spraying, is the soil, and, more particularly, the microorganisms it contains which are essential to maintain the functions of the soil ecosystem [1, 2]. So nowadays, European laws are converging to reduce the amounts of pesticides (Directive Cadre sur l’Eau 2000/60/CE, http://www.eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=CELEX:32000L0060:fr:html. Directive Cadre sur les pesticides 91/414/CEE, http://www.eur-lex.europa.eu/LexUriServ/LexUriServ.do?uri=CELEX:31991L0414:fr:html. Plan Ecophyto 2018, http://www.agriculture.gouv.fr/ECOPHYTO-2018/) used, and, in doing so, reduce their associated risks. The lack of controls in pesticides targets has caused the revaluation of many first-generation molecules and their replacement with new generations of molecules with lower weight, which, therefore, are less dangerous to the environment [3]. To ensure the efficiency of these new pesticides, the trend is to multiply the molecules used for the same crop in combinations called “cocktails.” Unfortunately, comparing to data on single pesticides (e.g., [1]), there are few studies concerned with the impact of these “cocktails” on the environment, particularly regarding their indirect impact on soil microbial communities. Some studies [4, 5], however, showed that when two herbicides were used simultaneously, their persistence in the soil was generally prolonged, and the environmental contamination potential increased. Among others [5, 6], Tejada showed an increased toxic effect of “cocktails” on biological microbial activity. In order to better understand the impact of pesticides, it is essential that the potential synergetic effects of combined herbicide mixtures on soil microbial communities are investigated further.

In this context, to replace the herbicide Atrazine, which was banned in 2003 [7], the triketone herbicide Mesotrione (2-[4-(methylsulfonyl)-2-nitrobenzoyl]-1,3-cyclohexanedione) is often used in combination with the chloroacetanilide herbicide S-metolachlor (2-chloro-N-(2-ethyl-6-methylphenyl)-N-[(1S)-2-methoxy-1-methylethyl] acetamide) on maize crops.

Mesotrione is a chemical derivative of a natural phytotoxin produced by Callistemon citrinus [3] which inhibits a critical enzyme, p-hydroxy-phenylpyruvate dioxygenase (HPPD), in carotenoid biosynthesis [8]. It was recently shown that it could have some impact on soil microbial communities when used alone [9] and that some of these metabolites could be more toxic than the parent molecule [10].

S-metolachlor is classified as an inhibitor of very long chain fatty acid (VLCFA) formation (HRAC (2007), http://www.hracglobal.com/). It interferes with normal cell development and inhibits both cell division and cell enlargement [8]. The choice of S-metolachlor among the other herbicides sprayed together with Mesotrione, was made mainly because of its known difficult biodegradation [8] and because of its targets, which differ from those of Mesotrione.

Although we may know the impact of such pesticides when they are applied alone [1], it is necessary to study the “cocktails” used on fields to evaluate the real impact of pesticides on soil microbial communities. No studies were found in the literature which made reference to the effects of the combined application of Mesotrione+S-metolachlor on soil microorganisms. The objective of this paper, therefore, was to study the effect of the combination of these two molecules on soil microbial communities and ecological processes through the use of two herbicides commonly used for maize crops, being the formulated compounds Callisto and Dual Gold.

We assessed the impacts of the herbicide “cocktail” on global soil microbial activity, biomass, and structures. We also focused on the effect of the “cocktail” on potentially sensitive communities, such as phototrophic microorganisms and those involved in the nitrogen cycle. Our work was carried out using these herbicides at onefold field-rate, being the normal conditions of spraying, and at a tenfold field-rate, to replicate the higher pressure of pollution encountered sometimes in the first millimetres of soil after spraying. To conduct this study we privileged a microcosm approach as a reliable substitute for field research into the effects of pesticides on soil ecosystems [11].

2. Materials and Methods

2.1. Experimental Design

Topsoil (0–5 cm) was collected in January 2010 from a fallow chernozem soil (determined using standard protocols of Arras Soil Analysis Laboratory, INRA, France; and according to FAO classification) located in the Limagne plain near Clermont-Ferrand (France). The soil was then sieved (mesh size 2 mm) and preincubated under experimental conditions for two weeks, in large polyvinyl chloride (PVC) tanks.

The effects of Mesotrione and S-metolachlor herbicides on soil microbial communities when applied alone or in a “cocktail,” were studied using two different doses: onefold field rate: 1 ×FR and tenfold field rate: 10 ×FR of formulated Mesotrione (0,45 mg·kg−1 s o i l d r y w e i g h t ( d w ) of active ingredient (a.i.) for 1 ×FR and 4,5 mg·kg−1 soildw (a.i.) for 10 ×FR of Callisto, 100 g·L−1 Mesotrione, Syngenta Crop Protection) and formulated S-metolachlor (5,4 mg·kg−1 soildw (a.i.) for 1 ×FR and 54 mg·kg−1 soildw (a.i.) for 10 ×FR of Dual Gold, 915 g·L−1 S-metolachlor, Syngenta Crop Protection). The formulated herbicides were diluted in distilled water, and soil pH was measured after herbicide application to ensure a nonoccurrence of pH effects.

Each treatment (C1 and C10 for Callisto applied at 1 ×FR and 10 ×FR, respectively; DG1 and DG10 for Dual Gold applied at 1 ×FR and 10 ×FR, resp., M1 (C1 + DG1) and M10 (C10 + DG10) for the Callisto/Dual Gold Mix applied at 1 ×FR and 10 ×FR, resp.) was carried out on different soil subsamples in separate PVC tanks, by spraying and then mixing in the herbicides with a small rake. The control soil (T) was sprayed with distilled water only. Each soil sub-sample was divided into 20 microcosms (190 g of soil in PVC cylinders, height = 8 cm, diameter = 10 cm), leading to a total of 140 microcosms. Soil water content was adjusted to 35% of maximum water holding capacity (WHC), making it comparable to field moisture, and was monitored throughout the experiment by weighing the soil of each microcosm every week. The microcosms were covered with porous plastic film allowing air exchange while limiting soil dehydration and incubated in an experiment room with a day/night cycle (14 : 10) at 20°C ± 2°C for 93 days. Five randomly chosen microcosms from each treatment set were sacrificed on days 2, 8, 42, and 93. For each microcosm, microbial biomass, microbial respiration, N - N O 3 and N - N H 4 + (except at day 93), and photosynthetic pigments were obtained from fresh soil aliquots. The remaining soils were stored at −20°C for DNA analysis, herbicide quantification, enzymatic activity analysis, and carbohydrate and protein measurements.

2.2. Herbicide Dissipation

At the beginning of the experiment, soil samples were checked for the absence of Mesotrione, S-metolachlor, and other pesticides (multiresidues117 scan, LDA 26 Laboratory). In order to check the amount of herbicide initially spread on the soil and track its dissipation, Mesotrione and S-metolachlor concentrations were monitored for each treatment for three out of five replicate microcosms on every sacrificing day. Mesotrione was extracted and concentrations were quantified, as described by Crouzet et al. [9]. S-metolachlor was extracted by adding 10 mL of MeOH/H2O (4/1) onto 5 g of fresh soil, and stirring for 24 h using an orbital shaker at 50 rpm and at room temperature. Herbicide concentration of soil liquid extracts was quantified by high-performance liquid chromatography (HPLC) on an Agilent Series 1100 chromatograph (Courtaboeuf, France) equipped with a reverse-phase column (C18 Zorbax Eclipse column, 150 mm × 4.6 mm, 5 μm) at 22°C and a diode array detector set at 𝜆 = 2 5 4  nm for Mesotrione and 𝜆 = 2 2 0  nm for S-Metolachlor. The mobile phase was composed of acetonitrile (A) and aqueous H3PO4 (0.01% v/v, pH = 2.9) (B) at a flow rate of 1 mL·min−1. Gradient: 0–5 min: 5% A; 5–30 min: 5–70% A (linear); 30–32 min 70–95% A; 32–33 min: 95–5% A: 33–36 min: 5% A.

2.3. Soil Chemical Analysis

Physicochemical characteristics of Limagne soil have already been described in a previous study [9] and checked in this study for any evolutions since the last analyses.

Evolutions of carbohydrates and proteins in soil samples were followed, according to Dubois et al. [12] and Lowry et al. [13], respectively, to estimate the availability of the principal organic compounds which ensure microbial development. Results were read spectrophotometrically on a SFM 25 (Kontron Instruments).

Measurements of ammonium ( N H 4 + ) and nitrate ( N O 3 ) concentrations were performed within hours following the sampling. For N H 4 + , 2 g of soil were mixed with 10 mL of 1 M KCl and stirred for 1 h at 250 rpm at room temperature. For N O 3 , 10 g of soil were mixed with 50 mL of deionized water at 40°C and stirred for 1 h at 250 rpm at room temperature. After centrifugation at 4000 ×g during 10 min, supernatants were filtered on a cellulose filter (Whatman 42). N H 4 + measurements were performed spectrophotometrically on a SFM 25 (Kontron Instruments), directly after filtration, following the Berthelot reaction method [14]. N O 3 measurements were performed by ionic chromatography [15]. Samples were run in a Waters 501 HPLC pump (Millipore) at 0.5 mL·min−1 and conductivity was followed by a Waters 431 detecting device (Millipore). Data was recorded and analysed by TotalChrom Navigator software (PerkinElmer).

2.4. Estimated Soil Microbial Activity and Biomasses

The estimated activity of soil microbial communities was quantified using the substrate-induced respiration (SIR) method. SIR was measured according to Crouzet et al. [9] with 6 mg glucose g−1 soil as performed by Lin and Brookes [16], using a gas chromatograph (Agilent 3000 μGC, Agilent Technology). SIR results were not converted into biomass values [17] and used as a proxy of the activity level of the heterotrophic active microbial communities.

The total soil microbial biomass was quantified using the fumigation/extraction method described by Vance et al. [18] with slight modifications. Briefly, 10 g of fresh soil was fumigated in a bell jar vacuum system containing 250 mL of CHCl3 stabilized with 0.6% EtOH at 37°C during 24 h. The total carbon (TC) extraction from fumigated and nonfumigated soil samples was performed by adding 40 mL of 30 mM K2SO4 and stirring for 40 min at 250 rpm at room temperature. After centrifugation at 4000 ×g during 10 min, supernatants were filtered on 1.2 μm glass microfiber filters (GF/C Whatman). Analyses of TC and inorganic carbon (IC) were performed on a total organic carbon analyzer TOC-5050A (Shimadzu). Each sample was run a minimum of 6 times to get correct values ( 𝑃 value < 0.05). Indirect microbial biomass was obtained by using Lin and Brookes general formula: Bc = 2,22 × Ec (with Bc, the microbial C. and Ec the K2SO4-extracted C. in fumigated soil less the K2SO4-extracted C. in nonfumigated soil) [16].

The photosynthetic microbial biomass was quantified by photosynthetic pigment extraction and measurement (chlorophyll a, active chlorophyll 𝑎 , and carotenoids). Pigments were extracted within hours following the sampling, according to Tsujimura et al. [19]. The concentrations of pigments were measured spectrophotometrically from absorbance readings (480, 630, 645, 663, and 750 nm) on a CE 2021-2000 series (CECIL Instruments) using the equation described by SCOR-UNESCO Working Group 17 [20].

2.5. Soil Microbial Structures
2.5.1. Bacterial Genetic Structure by PCR-DGGE

Bacterial community structure from each microcosm was investigated on genomic DNA extracted and purified from 1 g of soil using the Ultra Clean Soil DNA isolation kit extraction (MoBio Laboratories) according to the manufacturer’s instructions for maximum yields. The quality of DNA extractions was checked on NanoDrop 1000 (Thermo Fisher Scientific). The amplification was run with the primers 968f-GC and 1401r (Table 1) targeting the V6–V8 regions of bacterial 16S rDNA [21]. The PCR reaction (50 μL) contained 1 ×PCR buffer (Bioline), 2.5 mM MgCl2 (Bioline), 200 μM of each dNTP (Promega), 0.3 μM of each primer, 500 ng·μL−1 of bovine serum albumin (BSA, Sigma), 1.5 units of BIOTAQ DNA polymerase (Bioline), and 40 ng of soil DNA extract. The samples were amplified in MJ Research thermocycler (PTC-200). PCR products were checked by 0.8% (wt/v) agarose gel electrophoresis, stained with ethidium bromide, and quantified with the DNA quantitation Kit (Sigma) using a Kontron spectrofluorometer (Kontron Instruments) at wavelength of 360 nm (excitation) and 460 nm (emission). 500 ng of each bacterial PCR product was loaded into 7.5% polyacrylamide (w/v) denaturing gels with linear gradients of 50–60% (100% denaturant contains 7 M urea and 40% formamide) and run in 1X TAE buffer at a temperature of 60°C at 80 V for 21 h on a TTGEK-2401-220 (CBS Scientific).

2.5.2. Fungal Genetic Structure by PCR DGGE

Fungal community structure was investigated on genomic DNA extracted and purified from a pool of 7.5 g of soil from the five randomly chosen microcosms for each treatment with PowerMax Soil DNA isolation kit extraction (MoBio Laboratories) according to the manufacturer’s instructions for maximum yields. Quality of DNA extractions was checked on NanoDrop 1000 (Thermo Fisher Scientific). The amplification was run with the primers ITS5-GC and ITS2 (Table 1) targeting the ITS and 5.8S rDNA region [22]. The PCR reaction (50 μL) contained 1 ×PCR buffer (Qiagen), 2.5 mM MgCl2 (Qiagen), 200 μM of each dNTP (Promega), 0,3 μM of each primer, 500 ng·μL−1 of bovine serum albumin (BSA, Sigma), 1.25 units of HotStart Taq DNA polymerase (Qiagen), and 40 ng of soil DNA extract. The samples were amplified in UnoCycler thermocycler (VWR). PCR products were checked as for bacteria. 500 ng of each fungal PCR product was loaded into 7.5% polyacrylamide (w/v) denaturing gels with linear gradients of 25–50% (100% denaturant contains 7 M urea and 40% formamide) and run in 1X TAE buffer at a temperature of 59°C at 70 V for 16 h on the same apparatus as for the bacterial PCR products.

After electrophoresis, the DGGE gels were stained in 1X TAE buffer containing 1/20000 dilution of Gel Star (Lonza, Rockland, ME-USA) and digitized using a BioSpectrumAC Imaging System (UVP).

2.6. Statistical Analyses

All soil microbial measurements were calculated on the basis of oven-dried soil weight and expressed as the mean values (±SEM) of five replicates per treatment and per date, except for bacterial DGGE results, where only three replicates were analysed, and for fungal DGGE results. Statistical procedures were carried out using PAST software [23]. Normality distribution and homogeneity of variances were checked for all data. One-way ANOVA was used to test treatment effects on the microbial activities and biomasses, for a given date after herbicide application. Differences between individual means were then compared using Tukey’s honest significant difference (HSD) post hoc test. Significance was set at 𝑃 0 . 0 5 . Gel analyses were performed using Gel ComparII software (Applied Maths, Kortrijk, Belgium). According to Fromin et al. [24], total band intensity for each lane was normalized among lanes and each DNA band was considered in the analysis if its relative intensity accounted for more than 0.5% of total lane intensity. Data obtained was compared to control results for each sampling date, using the Student’s t-test for two samples. Based on band position and band relative intensity analysis, the diversity of the microbial communities was analysed by the richness ( 𝑅 ) estimated as the number of bands, the Shannon-Weaver index (H) according to the following equation: 𝐻 = P i (ln Pi) [25], where Pi is the ratio between specific band intensity and total intensity of all bands in a lane sample.

3. Results

3.1. Herbicide Dissipation

Soil herbicide concentrations just after contamination were very similar to the expected values for all treatments (data not shown). Herbicide dissipation was strongly dependent on application rates. At the same field rate, soil dissipation kinetics of S-metolachlor were similar between microcosms containing Dual Gold applied alone or in a “cocktail” (Figure 1), with a DT50 of 28 days at 1 ×FR and 85 days at 10 ×FR. Cumulative losses at day 93 represented around 89% and 53% of the initial S-metolachlor input for 1 ×FR- and 10 ×FR-treated soils, respectively. At day 42, microcosms containing Callisto showed Mesotrione concentration at 1 ×FR to be under the detection limit of the device used for the analyses, in both application methods. For the Callisto at 10 ×FR, Mesotrione concentrations at days 42 and 93 were statistically lower in microcosms where herbicides had been applied alone. Cumulative losses at day 93 represented around 97.5% and 80% of the initial Mesotrione 10 ×FR input for single application and “cocktail” treated soils, respectively. The “cocktail” increased the DT50 of Mesotrione (DT50 of 40 days in “cocktail” against 28 days in alone application).

3.2. Effects on Edaphic Parameters

Among the initial physicochemical soil characteristics, organic matter (OM) concentration was at 29.9 g·kg−1 with a pH value of 6.5 (supplementary data 1).

N H 4 + concentrations increased during the first 42 days of experimentation (Figure 2). Two days after the treatment, N H 4 + concentrations were statistically higher in C10 microcosms, compared to T, M1, and M10. The results evolved at day 8 and 42, with statistically higher concentrations of N H 4 + in DG10, M1, and M10 microcosms, compared to T and C1.

Whatever the treatment or the sampling date, no effect was observed on N O 3 (average values ranged from 62.3 to 107.8 μg·g−1dw) and proteins concentrations (average values ranged from 0.67 to 1.33 mg·g−1dw) in microcosms (data not shown).

Carbohydrate concentrations stayed quite similar, whatever the treatment and the sampling date (ranging from 2.18 to 2.62 mg·g−1dw), except at day 8 where statistical increases of glucose equivalent in DG10, M1, and M10 microcosms were observed, compared to control microcosms (Figure 3).

3.3. Effects on Microbial Biomasses and Activity

A statistical decrease of microbial carbon was observed between day 2 (values ranged from 385 to 449 μg·g−1dw) and day 93 (values ranged from 319 to 384 μg·g−1dw) (supplementary data 2). At each sampling date, the total soil microbial carbon value showed no statistical differences between treatments.

The general pattern observed for the SIR experiment was a global increase in CO2 of 0.015% per hour for all microcosms. No difference was observed between control and treatments, whichever day was considered (data not shown).

The photosynthetic microbial biomass, assessed by measuring photosynthetic pigments, showed no statistical difference between treatments until the 42nd day of experimentation, where total chlorophyll, active chlorophyll 𝑎 , and carotenoid concentrations displayed similar patterns of responses for the different treatments. At day 93 (Figure 4), whereas chlorophyll 𝑎 concentrations in T and C1 microcosms were statistically higher than other treatments, carotenoid concentrations in T microcosms were only statistically higher than the DG10 and M10 microcosms. Active chlorophyll 𝑎 was monitored to give an estimation of active biomass of autotrophic communities and showed statistically less active chlorophyll 𝑎 in C10, M1, and M10 microcosms compared to T. Moreover, in C1 microcosms, active chlorophyll 𝑎 was always higher than other treatments whereas carotenoid was different to DG1, DG10, M1, and M10 microcosms.

3.4. Effects on Microbial Community Structures

The diversity of the results obtained from DGGE gels showed a large band number, with 24 to 53 for bacterial gels and from 34 to 48 for fungal gels. Community structures showed a complex diversity with major bands commonly found in all microcosms and a lot of minor bands. Some communities, therefore, seemed to be stimulated by the presence of herbicides in microcosms with one or both herbicides (supplementary data 3).

At day 2, the bacterial diversity results (Table 2) showed that richness in C10 microcosms was much higher than control (T). At the same day, however, the Shannon-Weaver index ( 𝐻 ) showed a decrease of overall diversity in C1 microcosms and an increase in C10, DG10, and M1 microcosms. The only differences observed at days 8 and 42 were a decrease of DG1 band numbers. At day 93, C10 and M1 microcosms presented poor richness compared to T microcosms, and in M1 and M10 microcosms the diversity indicator ( 𝐻 ) was statistically lower.

A global decrease of richness in fungal communities was observed between day 2 and day 93 ( 𝑃 0 . 0 5 ), but there was no change in the 𝐻 index.

4. Discussion

This study focused on analysing the unknown effects on soil microbial diversity and activities of a “cocktail” of herbicides used for maize crops. As far as we know, it is the first study using the Mesotrione+S-metolachlor mixture, despite the fact that these molecules are recommended in many pesticide applications to control undesirable vegetation. The effects recorded may have resulted from the combined impact of herbicides, adjuvants, and degradation products.

4.1. Herbicides Fate in Soil

Despite the fact that nonwater extractable residues have not been quantified in this study, the known Koc values of both herbicides led us to consider that there was little possibility of such formations occurring [26, 27] (Agritox, http://www.dive.afssa.fr/agritox/index.php/). Furthermore, as Mesotrione and S-metolachlor have little sensitivity to photodegradation [28] (Agritox, http://www.dive.afssa.fr/agritox/index.php/), we can conclude that biodegradation could explain their decreasing concentration in soil microcosms (Figure 1). Moreover, the 35% of maximum WHC represented optimal conditions for microbial growth and allowed a good distribution of herbicides in the soil for biodegradation. Mesotrione decrease was faster than S-metolachlor decrease, which complies with Crouzet et al. [9] DT50 results for Mesotrione, and Ma et al. [29] for S-metolachlor. Mesotrione dissipation was longer when associated with S-metolachlor and its adjuvants at 10 ×FR. This observation has already been noted by other authors who have demonstrated that the combined use of herbicides in a “cocktail” [4, 5], or the association of herbicides and fungicides [30], could increase their persistence in soil. We can subsequently conclude that the S-metolachlor and/or additive compounds (adjuvants) contained in Dual Gold, which were certainly different from those in Callisto, could act on Mesotrione availability by modifying solubility or sorption parameters [31] and/or impact microorganisms potentially capable of Mesotrione degradation, resulting in herbicide persistence in the soil. Additional compounds, however, may be responsible for the effects observed in the “cocktail” microcosms: for example, alcohol ethoxylates (AEOs) and alkylamine ethoxylates (ANEOs) pesticide adjuvants have demonstrated their toxicity on small eukaryotes and bacteria [32] even at low concentrations.

4.2. How do Microbial Communities Respond to Herbicide Applications?

Herbicides were originally designed for their effects on plants. They are, however, potentially capable of directly or indirectly altering soil microbial communities which include organisms sensitive to the parent molecule and/or degradation products [1, 33]. In this study, the pollution pressure, in terms of active ingredients, was found to remain until the end of the experiment for all the treated microcosms, except for C1 (Figure 1).

When applied alone and at the recommended field rate, Callisto (C1) showed no statistical differences compared to the control microcosms, except for a small decrease in bacterial diversity assessed using the Shannon-Weaver index (Table 2). These results are close to those observed by Crouzet et al. [9], who concluded that this herbicide when applied at the recommended field rate, had a nonsignificant impact on soil microbial communities. The soil used in our microcosms was different to that used by Crouzet et al. in terms of percentage of sand and silt and pH values (supplementary data 1). The increase of the porosity (37% of sand against 23% in 2007) could have modified the retention properties of molecules in the soil by increasing the Mesotrione access to microorganisms, thereby ensuring faster degradation and short-term impacts. The results obtained for Dual Gold treatment (DG1) with the same application dose, showed a higher impact on microbial communities, with a decrease of bacterial richness at days 8 and 42, and a decrease of chlorophyll 𝑎 concentrations at day 93. We can conclude, therefore, that these herbicides, when applied at the recommended field rate, acted differently on microbial communities and produced nondesired effects. Dual Gold impacted microbial communities longer than Callisto, however, at the end of the experimentation, C1 and DG1 (except for chlorophyll 𝑎 concentration) were similar to T microcosms.

According to the SIR results, microbial communities were equally active in control and treated microcosms throughout the whole experiment. The SIR results showed the same pattern as the one observed in the Zabaloy et al. study [34] for 2,4-D and metsulfuron-methyl at 10 ×FR, but differed from others [9, 35], where pesticide supplementation increased the respiration activity. Currently, no general pattern of microbial responses can be inferred regarding pesticide identity, doses applied, exposure time, soil type, or other environmental factors [3638]. Moreover, the total microbial biomass and protein concentrations showed no statistical differences between treatments during the experimentation at any given sampling date.

The absence of effects in SIR, microbial carbon, and protein measurements could be explained by the hypothesis submitted by Wardle and Parkinson [39], whereby a microbial community under the effects of herbicides was presumably in a considerable state of flux, with susceptible microbes being killed and others, thereby, having a readily available source of carbon.

With reference to the results concerning carbohydrates (Figure 3), herbicides in “cocktail” microcosms (M1 and M10) at day 8 could produce such an effect by killing microorganisms, such as bacteria and fungi, which would result in a decrease of bacterial richness and diversity (Table 2) and an increase in sugar concentrations. Moreover, increases of labile substrates and organic matter could be used by the resistant population (supplementary data 3) and influence the nitrogen cycle, which would explain the increase in ammonium found in “cocktail” microcosms at days 8 and 42. Similar results were observed with C10 treatments as early as day 2, where a short-term effect on bacterial communities had been noticed, with increases in some populations (Table 2) and ammonium concentrations (Figure 2). This could not be considered, however, as a general effect of pesticides, as underlined by Lo in his review [1], where pesticides produced both increases and decreases in nitrogen activities.

Another explanation for this increase in ammonium concentration could be an inhibition of the nitrification process resulting in a reduced use of ammonium by microbial communities. Notwithstanding the above, herbicide treatments did not modify the amounts of nitrates. This leads us to consider the possibility of small effects made by ammonium shifts, not recorded by our analyses, given the high values of nitrate measured, equivalent to 102 to 103 superior than ammonium concentrations. Additional data would be necessary, however, to link this increase in ammonium to a specific community. Although it was not one of the objectives of this study, it would be prudent to use specific primers to target those microbial communities involved in the nitrogen cycle to provide data to confirm this hypothesis [40, 41].

At the end of the experimentation, bacterial structure results (Table 2) showed a decrease in microbial diversity for “cocktail” microcosms. This is not linked to microbial abundance results and might be masked by the global effects on microbial communities [39]. It underlined, however, a long-term synergetic effect of the application of “cocktails.”

Another impact is clearly recorded (Figure 4) on the same day, being the colonization of the surface of “cocktails” and C10 microcosms by phototrophic communities. The impact of Callisto on photosynthetic pigments, particularly on the HPPD enzyme, in some phototrophic microbial communities, inhibited the colonization of such communities throughout the experimentation period. The “cocktail” affects total and active chlorophyll 𝑎 in ideal growth conditions [19], primarily because of the duality of the targets covered by the two herbicides for phototrophic communities.

These results underlined the fact that a simple “cocktail” at 1 ×FR can impact the microbial community as much as a single molecule used at 10 ×FR. Furthermore, “cocktail” results are closer to those obtained for 10 ×FR Dual Gold applications (DG10) than to 10 ×FR Callisto applications (C10).

5. Conclusion

The results obtained in this study clearly showed that when applied alone at the recommended field rate doses, the herbicides Callisto and Dual Gold exert only few consistent minor effects on soil microbial communities. The study confirmed, however, that a “cocktail” of such herbicides applied at 1 ×FR, has an effect on microbial communities also confirmed at 10 ×FR.

The results of this study show the deleterious effect of even low quantities of herbicide “cocktails” on soil microbial communities and must be used with caution despite the supposedly low dangerousness of these new molecules for the environment. Knowing that pesticide “cocktails” of herbicides, fungicides, insecticides, and many other molecules are the current trend in agricultural practices, it is essential that we investigate more precisely their impact with more sensitive methodologies and approaches, and focus on specific microbial communities which ensure key functional steps in biogeochemical cycles. Functional communities involved in N-cycling provide good microbial ecology models to study the role of size and structure of microbial communities on corresponding process rates and ecosystem functions [40, 42]. Moreover, phototrophic communities could be used as bioindicator of pesticide “cocktails” stresses. Any further studies should not underestimate the number of pesticides which can be applied in combinations or “cocktails,” and must focus on specific parameters and specific microbial communities to obtain an overall understanding of their potential impact.

Conflict of Interests

The authors declare that there is no conflict of interests.

Acknowledgments

The authors would like to thank D. Tottey, English native translator, for the English review of this paper, B. Debote from Crouël INRA (Institut National de la Recherche Agroalimentaire) for providing soil samples, G. Maillot from LPMM (Laboratoire Photochimie Moléculaire et Macromoléculaire) for technical assistance on carbon measurements, J. Lavedrine from UBP (Université Blaise Pascal) Chemistry Department for technical assistance on nitrate measurements, R. Nouaille from LMGE for technical assistance on the SIR measurements, and S. Fontaine from Crouël INRA for technical assistance on fumigation/extraction experiments.

Supplementary Materials

The supplementary material provide information on (1) the physicochemical properties of the soil used in the study, (2) the carbon values from the soil microbial biomass and (3) a set of fingerprints data used to obtain the table 2 results.

Supplementary Data 1. Physicochemical properties of the Limagne basin soil (chernozem soil according to FAO classification), determined using standard protocols of the Arras Soil Analysis Laboratory (INRA, France).

Supplementary Data 2. Carbon from soil microbial biomass (μg.g-1dw), according to Lin and Brookes formula [16], depending on soil treatments: control (T), Callisto® 1×FR (C1) and 10×FR (C10), Dual Gold® 1×FR (DG1) and 10×FR (DG10), and Callisto®/Dual Gold® Mix 1×FR (M1) and 10×FR (M10) and sampling dates.

Supplementary Data 3. DGGE gels of (a) bacterial communities (n=3) at 2 days after treatment: control (T), Dual Gold® 10×FR (DG10), and Callisto®/Dual Gold® Mix 1×FR (M1) and 10×FR (M10), and (b) fungal communities at 42 and 93 days after treatment: control (T), Callisto® 1×FR (C1) and 10×FR (C10), Dual Gold® 1×FR (DG1) and 10×FR (DG10), and Callisto®/Dual Gold® Mix 1×FR (M1) and 10×FR (M10).

  1. Supplementary Data 1
  2. Supplementary Data 2
  3. Supplementary Data 3