Veronika Storck1†
Sofia Nikolaki2†
Chiara Perruchon3Camille Chabanis4Angela Sacchi5Giorgia Pertile5Céline Baguelin4Panagiotis A. Karas3
Aymé Spor1Marion Devers-Lamrani1Evangelia S. Papadopoulou3Olivier Sibourg4
Cedric Malandain4
Marco Trevisan6Federico Ferrari5
Dimitrios G. Karpouzas3
George Tsiamis2*
Fabrice Martin-Laurent1*- 1AgroSup Dijon, Institut National de la Recherche Agronomique, Université Bourgogne Franche-Comté, Agroécologie, Dijon, France
- 2Department of Environmental and Natural Resources Management, University of Patras, Agrinio, Greece
- 3Laboratory of Plant and Environmental Biotechnology, Department of Biochemistry and Biotechnology, University of Thessaly, Larissa, Greece
- 4Enoveo srl., Lyon, France
- 5Aeiforia srl, Spinoff Università Cattolica del Sacro Cuore, Fidenza, Italy
- 6Department of Agronomy and Environmental and Chemistry, Catholic University of the Sacred Heart, Piacenza, Italy
Pesticides are intentionally applied to agricultural fields for crop protection. They can harm non-target organisms such as soil microorganisms involved in important ecosystem functions with impacts at the global scale. Within the frame of the pesticide registration process, the ecotoxicological impact of pesticides on soil microorganisms is still based on carbon and nitrogen mineralization tests, despite the availability of more extensive approaches analyzing the abundance, activity or diversity of soil microorganisms. In this study, we used a high-density DNA microarray (PhyloChip) and 16S rDNA amplicon next-generation sequencing (NGS) to analyze the impact of the organophosphate insecticide chlorpyrifos (CHL), the phenyl-urea herbicide isoproturon (IPU), or the triazole fungicide tebuconazole (TCZ) on the diversity and composition of the soil bacterial community. To our knowledge, it is the first time that the combination of these approaches are applied to assess the impact of these three pesticides in a lab-to-field experimental design. The PhyloChip analysis revealed that although no significant changes in the composition of the bacterial community were observed in soil microcosms exposed to the pesticides, significant differences in detected operational taxonomic units (OTUs) were observed in the field experiment between pesticide treatments and control for all three tested pesticides after 70 days of exposure. NGS revealed that the bacterial diversity and composition varied over time. This trend was more marked in the microcosm than in the field study. Only slight but significant transient effects of CHL or TCZ were observed in the microcosm and the field study, respectively. IPU was not found to significantly modify the soil bacterial diversity or composition. Our results are in accordance with conclusions of the Environmental Food Safety Authority (EFSA), which concluded that these three pesticides may have a low risk toward soil microorganisms.
Introduction
Soil is a limited and largely non-renewable resource providing a unique and complex habitat for a wide range of microorganisms which support key soil ecosystem functions and contribute to complex processes with impacts at the global scale (Blum, 2006; Bondeau et al., 2007; Graham et al., 2016). For instance, soil microorganisms contribute to food supply, water quality, carbon cycling, and climate regulation (Haygarth and Ritz, 2009; Hillel, 2009). In this context, the preservation of soil quality by agriculture was identified to be one of the greatest objectives of the twenty-first century (Lal, 2008; Jones et al., 2009; Lichtfouse et al., 2009; Rattan, 2015; Steffan et al., 2017). To achieve this objective, the main pressures on soil organisms were weighted by the European soil biodiversity expert group and human intensive exploitation was ranked first (Jeffery et al., 2010; Gardi et al., 2013). Pesticides are among major contributors to the pressure exerted on soil organisms. They are intentionally applied in conventional agriculture to protect crops from various pests. They can persist in soil, and from there they can be dispersed to other environmental compartments (Looser et al., 2000; Barbash, 2003; Smalling et al., 2013), and can harm non-target organisms (Pimentel, 1995; Zhang et al., 2011) including soil microorganisms as summarized by Puglisi (2012). The preservation of the soil bacterial diversity is important as it largely contributes to crucial ecosystem functions. While a well-preserved high soil bacterial diversity may ensure high rates of ecosystem processes under changing and even extreme environmental conditions, a low diversity may lead to weak ecosystem process rates under extreme environmental conditions. Thus, it is important to evaluate a potential impact of pesticides on the soil bacterial diversity and composition (Imfeld and Vuilleumier, 2012).
The evaluation of the ecotoxicological effects of pesticides on soil microorganisms remains a difficult task, as research results and their complex interpretation often differ from one study to another. Hence a more standardized approach both in terms of methods and risk assessment schemes is required. Effects of the organophosphate insecticide chlorpyrifos (CHL), the phenyl-urea herbicide isoproturon (IPU), and the triazole fungicide tebuconazole (TCZ) are particularly interesting to study as these pesticides are on the market for 53, 43 and 30 years respectively and they have been identified as common contaminants of natural water resources (Hayes and Laws, 1991; Oerke et al., 1994; Börner et al., 2009). As happens with all pesticides currently on the market, the EU regulatory assessment of the soil microbial ecotoxicity of CHL, IPU, and TCZ solely relies on tests looking at their ecotoxicological impact on carbon and nitrogen mineralization in microcosm studies. These tests, presented in the relevant regulatory documents by the Environmental Food Safety Authority (EFSA), who is in charge for the authorization of an active substance of a pesticide for all EU countries, suggest no significant impact on carbon and nitrogen mineralization in soil for CHL (EC, 2005) and effects of <25% for IPU and TCZ, which are considered as an acceptable low risk to soil microorganisms by the regulatory authorities (EFSA, 2008, 2014, 2015a). Carbon and nitrogen mineralization tests are global indicators of microbial activity which can be criticized as being relatively insensitive, hence not providing a realistic assessment of the impact of pesticides not only on microbial functions but mostly on the soil microbial diversity that should be preserved.
To date, several studies have investigated the impact of CHL, IPU, and TCZ on soil microbial mass, diversity or activity. For all three pesticies, studies on their effects were mostly performed in soil microcosms (CHL: Pozo et al., 1995; Chu et al., 2008; Fang et al., 2009; Dutta et al., 2010; Srinivasulu and Rangaswamy, 2013; IPU: Tag-El-Din, 1982; Harden et al., 1993; Kuriyal and Pandey, 1994; Perrin-Ganier et al., 2001; TCZ: Cycon et al., 2006; Bending et al., 2007; Ferreira et al., 2009; Munoz-Leoz et al., 2011; Anuradha et al., 2016; Wang et al., 2016) and seldomly at field scale (CHL: Eisenhauer et al., 2009; Gupta et al., 2013; IPU: Schuster and Schröder, 1990). Overall, most of the studies reported small to moderate (sometimes temporary) toxic effects of the three pesticides to soil microbial biomass, diversity or activity (CHL: Pozo et al., 1995; Singh et al., 2002; Chu et al., 2008; Eisenhauer et al., 2009; Fang et al., 2009; Dutta et al., 2010; Gupta et al., 2013; Srinivasulu and Rangaswamy, 2013; IPU: Tag-El-Din, 1982; Schuster and Schröder, 1990; Harden et al., 1993; Kuriyal and Pandey, 1994; Perrin-Ganier et al., 2001; TCZ: Cycon et al., 2006; Bending et al., 2007; Ferreira et al., 2009; Munoz-Leoz et al., 2011; Anuradha et al., 2016; Wang et al., 2016; TCZ: Cycon et al., 2006; Bending et al., 2007; Ferreira et al., 2009; Munoz-Leoz et al., 2011; Anuradha et al., 2016; Wang et al., 2016). It is noteworthy that several studies showed pesticide exposure to cause short-term inhibitory effects within the first month after treatment, which were recovered in subsequent months in accordance with pesticide dissipation (CHL: Pozo et al., 1995; Fang et al., 2009; TCZ: Ferreira et al., 2009; Munoz-Leoz et al., 2011).
Almost all ecotoxicological studies about pesticide effects on soil microorganisms applied classical approaches to quantify impacts on the abundance (microbial mass), structure [phospholipid fatty acid (PLFA) analysis, DNA fingerprinting approaches such as T-RFLP or A-RISA] or activity (soil respiration, soil enzymatic activity) of soil microbial communities. The latest developments of next-generation DNA sequencing (NGS) allows the extensive description of the diversity of the soil microbial community. Recent studies used NGS approaches to assess the impact of pesticides on the diversity of the soil microbial community mostly at microcosm scale level (Feld et al., 2015; Newman et al., 2016; Romdhane et al., 2016; Papadopoulou et al., 2018), while equivalent field scale studies performed under realistic exposure conditions are still lacking.
In this study, we used NGS and high-density microarrays (PhyloChip) to assess the potential ecotoxicological effects of these three pesticides on the diversity and composition of the soil bacterial community. This was done following a lab to field experimental approach, consisting of laboratory tests with higher pesticide exposure regimes (0x, 1x, 2x, and 10x the recommended agronomical pesticide dose, tier 1), and a more realistic exposure scheme at field scale (0x, 1x, 2x, and 5x the recommended dose, tier 2). To our knowledge, this is the first time that advanced molecular tools (high-density DNA microarray and 16S rDNA amplicon sequencing) were used to enable a high-resolution analysis of pesticide impacts on even less abundant members of the soil bacterial community.
Materials and Methods
Microcosm and Field Study
For the soil microcosm (lab) study, soil samples from an agricultural field situated in North Italy (45°05′20.8″N 9°45′59.4″E, Google Maps), an area intensively cultivated with winter cereal, were used. The soil was characterized as loamy sand (4.2% clay, 13.5% silt, 82.2% sand) with an organic carbon content of 1.5%, a microbial biomass of 160.2 mg C/kg soil dwt and a pH of 7.5. The field had not been treated with the studied pesticides for more than 5 years, while treatment by pesticides of the same chemical group cannot be excluded although this was not documented upon personal communication with the farmer. Soil samples were collected in July 2013 from the top 20 cm of the soil profile following the procedures described by the International Standardization Organization (ISO) for collection and handling of soil samples (ISO 10381-6, 2009). The soil was manually homogenized, partially air-dried, sieved (2 mm) and stored at 4°C for ~1 week before starting the experiment. For each pesticide dose combination, three replicates were prepared that were treated with an appropriate volume of aqueous solutions of CHL, IPU or TCZ (prepared from their commercial formulations Carposan® (CHL), Quintil® (IPU), or Folicur® (TCZ), respectively) aiming pesticide doses of 1x, 2x, or 10x the recommended agronomical dose [2.0 (CHL), 1.9 (IPU), and 0.6 (TCZ) mg/kg soil dwt]. Three soil samples were treated with water instead of pesticide solutions to serve as untreated controls (0x). Altogether, 150 experimental units were prepared [(3 untreated controls + 3 pesticides × 3 concentrations × 3 replicates) × 5 time points]. The water content of all soil samples was adjusted to 40% of the water holding capacity (WHC = 46.8 wt%), which was kept constant throughout the study. Each of the 150 experimental units (microcosms) consisted of 150 g wet soil which were placed in aerated (perforated) plastic bags and incubated in the dark at 20°C. Immediately after pesticide application (0 days) and 7, 42, 56, and 100 days later, experimental units were sacrificed and stored at −80°C until processed for DNA extraction.
The field study was performed on the site from where the soil samples for the microcosm study were collected. A detailed description of the set-up of the field study is given in Papadopoulou et al., (2016). Briefly, a completely randomized block design was followed with three replicate plots of 60 m2 (4 × 15 m) per pesticide-dose rate combination. I.e., the field plots were randomly ordered in a checkerboard pattern. The field was seeded with a mixture of cereals (60% Hordeum vulgare, 25% Triticum, and 15% Triticosecale) on November 7th, 2013. Replicate plots were treated with 1x, 2x, or 5x the recommended agronomical dose of each pesticide on the 12th of November 2013, while three plots were treated with water instead to serve as untreated controls (0x). This way, 30 experimental units (field plots) were prepared (3 untreated controls + 3 pesticides × 3 concentrations × 3 replicates). Soil samples (nine sub-samples per plot bulked to a composite sample per plot and sampling day) were collected at 0, 14, 35, 70, and 105 days after treatment from the top 20 cm of the soil profile and were subsequently used for DNA extraction. Sampled soil was manually homogenized, sieved (2 mm) and kept at −80°C until DNA extraction.
DNA Extraction
DNA was extracted from 300 mg soil samples (dry weight equivalent) following the ISO 11063 procedure (2012), derived from the method described by Martin-Laurent et al. (2001). The procedure involved three principal steps: (i) microbial cell lysis by physical (bead beating) and chemical (sodium dodecyl sulfate) actions, (ii) deproteination by precipitation, and (iii) DNA precipitation, washing and purification. Soil DNA extracts were purified using PVPP-columns and then Sepharose 4B-columns according to Petric et al. (2011b).
Microarray Sample Preparation and Analysis
For the microcosm study, the time points 0 and 56 days after pesticide treatment were examined. For the field study, 0 and 70 days after treatment were studied. The selection of these time points was based on the dissipation and metabolic patterns of the pesticides (Papadopoulou et al., 2016) and our intention to allow ample time for recovery of the soil microbiota after potential short-term effects.
Near-full-length 16S rDNA amplification for PhyloChip analysis was carried out using universal 16S rDNA primers 27F (5′-AGAGTTTGATCM TGGCTCAG-3′)and 1492R (5′-TACGGHTACCTTGTTACGACTT-3′)(Lane, 1991). A reaction mixture of 20 μL was prepared containing the PCR buffer (Takara Mirus Bio Inc., WI), 1.5 mM MgCl2, 0.25 mM of each dNTP, the primers at 0.3 mM each, and 1 U Taq polymerase (Takara Mirus Bio Inc., WI). PCR reactions were performed using a PTC-200 thermocycler (MJ Research Inc., USA) with a denaturation step of 10 min at 94°C, followed by 35 cycles for library construction and 30 cycles for the PhyloChip analysis of 1 min denaturation at 94°C, 1 min primer annealing at 52°C for the library construction, eight annealing temperatures between 48 and 58°C for the PhyloChip analysis, and finally 90 s extension at 72°C. The PCR was completed by a final extension of 10 min at 72°C. The size of the PCR product was determined by agarose gel electrophoresis using appropriate size markers.
16S rDNA microarray sample preparation was performed as previously described (Tsiamis et al., 2008). Briefly, per each pool of 16S rDNA amplicons, 0.5 μg were DNase I-fragmented and biotin-labeled. An aliquot of 0.1 μg was hybridized to PhyloChip at 60 rpm at 48°C overnight. Washing, staining and scanning of the PhyloChip was done according to the manufacturer's protocol (Affymetrix, USA). PhyloTrac version 2.1 was used to remove operational taxonomic units (OTUs) present on the PhyloChip but not detected in any of the samples. An OTU was regarded as present in the sample when over 90% of its assigned probe pairs were positive (PosFrac > 0.90). Data were accessed based on the hybridization intensity score in order to calculate relative abundance and as binary presence/absence for cluster analysis. Primer6+ version 6.1.13 was used to generate principal coordinate analysis (PCoA) to examine sample groupings. Statistical significance tests between groupings were performed using Permanova as implemented within Primer6+.
Amplicon Sequencing Samples Preparation and Data Analysis
One hundred and fifty microcosm and 150 field study DNA samples [5 time points × (3 controls + 3 pesticides × 3 doses × 3 replicates)] were analyzed using an amplicon Illumina next-generation sequencing approach. In a first-step PCR, fusion primers U341F (5′-CCTACGGGRSGCAGCAG-3′) and 805R (5′-GACTACCAGGGTATCTAAT-3′) (Klindworth et al., 2012) were used to amplify the hypervariable V3-V4 region of the bacterial 16S rDNA (464 bp). Primer sequences contained corresponding overhang adapters (forward adapter: TCGTCGGCAGCGTCAGATGTGTATAAGAGACAG, reverse adapter: GTCTCGTGGGCTCGGAGATGTGTATAAGAGACAG) needed to add multiplexing index-sequences in a second-step PCR.
Concerning the samples from the microcosm study, each PCR reaction was carried out in a volume of 25 μL containing 1 μL of template DNA, 5 μL Kapa 5x High-Fidelity Hot Start Buffer, 0.3 μL of each primer (25 μM), 0.75 μL dNTPs (10 mM), High-Fidelity DNA Kapa polymerase (1 unit/μL) and 17.15 μL water. PCR reactions were performed using a PTC-200 thermocycler (MJ Research Inc., USA). Cycling conditions were 5 min at 95°C, 30 cycles at 98°C for 20 s, 60°C for 30 s, and 72°C for 45 s, followed by a final extension of 5 min at 72°C. The size of the PCR products was visually checked with agarose gel electrophoresis using appropriate size markers. The PCR products were precipitated by polyethyleneglycol (20% PEG, 2.5 M NaCl) (Hartley and Bowen, 2003) and resuspended in 15 μL water. The concentrations of the PCR products were determined with a NanodropQuawell UV-VIS spectrophotometer (Q5000). The PCR amplicons were then used as templates in a second-step PCR for addition of the multiplexing index-sequences to the overhang adapters. Amplicons were purified with AMPure XP beads and re-suspended in 30 μL water. Products were quantified using the Qubit dsDNA High-Sensitivity assay (Life Technologies), and a Bioanalyzer High-Sensitivity DNA chip (Agilent). Amplicons were pooled at equimolar concentrations and the library was quantified using a Qubit Fluorometer and dsDNA HS Assay kit (Life Technologies, CA, USA). Size-selection was done by Pippin-Prep (Sage Science), where a size range of 350–450 bp was extracted. Reaction negatives were included in the pool. Size-selected fragments were purified and quantified in the same manner as described above prior to sequencing. Sequencing was performed at IMGM SA on an Illumina MiSeq platform using a 300 bp paired-end read chemistry (Illumina).
Concerning the samples from the field study, first-step PCR amplifications were carried out in 50 μL reaction volumes containing 2.5 U of Accuzyme polymerase (Bioline), 0.2 μM of each primer, 1x AccuBuffer with 2 mM Mg2+, 0.5 mM of each dNTP, 0.4 μg/μL of bovine serum albumin (BSA, Sigma) (Kreader, 1996) and 2 ng of soil DNA. The thermal cycling conditions were 95°C for 5 min, followed by 35 cycles of 95°C for 15 s, 58°C for 15 s, and 72°C for 1 min, with a final extension of 72°C for 5 min. The PCR products were visualized on a 1.5% agarose gel to verify the correct size of amplicons. Amplicon sizes of randomly selected samples were further analyzed by Bioanalyzer DNA 1000 chip (Agilent Technologies). The amplicons were cleaned-up with a Nucleospin gel and PCR clean-up kit (Macherey-Nagel) following the manufacturer's instructions. Subsequently, amplicon concentrations were determined with NanoDrop 2000c (Thermo Scientific). For the addition of multiplexing index-sequences to the overhang adapters, a second-step PCR was performed using a 384 libraries Nextera XT index kit v2 (Illumina). Second-step PCR amplifications were carried out in 50 μL reaction volumes containing 1x Titanium Taq polymerase (Clontech), 1x Titanium Taq PCR buffer, 0.2 mM of each dNTP, 5 μL of each Nextera XT index primer using the TrutSeq Index Plate Fixture (Illumina). All cleaned-up amplicons from the first-step PCR were diluted to 10 ng/μL and added to the second-step PCR reaction mix at a final concentration of 1 ng/μL. Thermal cycling conditions were 96°C for 3 min, followed by 8 cycles of 95°C for 30 s, 55°C for 30 s and 72°C for 30 s, with a final extension of 72°C for 5 min. Amplicon sizes (464 bp) of randomly selected sampled were then analyzed by Bioanalyzer DNA 1000 chip (Agilent Technologies). The libraries were cleaned-up before quantification using AMPure XP beads (Beckman-Coulter Genomics) following the manufacturer's instructions. Amplicon concentrations were then fluorometrically measured with Qubit (Invitrogen) and converted to nM, assuming an average library size of 464 bp: Amplicon concentration in ng/μL × (660 g/mol × 464 bp) × 106 = amplicon concentration in nM. Libraries with amplicon concentrations > 40 nM were first diluted to 40 nM and then all samples were brought to 4 nM. For pooling of the libraries, 2 μL of each sample were brought together and denatured with 0.2 N NaOH and diluted to 4 pM in hybridization buffer HT1, following the Illumina manufacturer's instructions. The pooled amplicon library was then combined with a PhiX control library (Illumina) before the Illumina MiSeq run (2 × 300 bp) was performed.
Illumina Sequencing Data Analysis
The raw data sequences obtained from the Illumina next-generation sequencing of the samples from the microcosm and the field study were processed using similar pipelines. Raw sequences were submitted to the NCBI under the Bioproject numbers PRJNA385243 (microcosm study) and PRJNA381453 (field study). Reads were de-multiplexed and converted to FASTQ format. For the microcosm study, Cutadapt 1.2.1 (Martin, 2011) was used to trim Illumina adapter sequences from FASTQ files. Reads were trimmed if 3 bp or more of the 3' end of a read matched the adapter sequence. Sickle version 1.200 (Joshi and Fass, 2011) was used to trim reads based on quality: Any read with a quality score of <20, or <10 bp size after trimming were discarded. Paired-end reads were assembled, trimmed by length and further corrected using PandaSeq (Masella et al., 2012). For the field study, reads were quality-controlled and assembled using PEAR (Zhang and Sun, 2014). Unassembled reads and once-assembled reads outside the expected range were discarded. All subsequent analyses were conducted in QIIME 1.9.0 (Caporaso et al., 2010a). Sequences were clustered into OTUs using USEARCH (Edgar, 2010) by open-reference OTU picking. Chimeras were detected and omitted using the program UCHIME (Edgar et al., 2011) with the QIIME-compatible version of the SILVA aligned version of the Gold database (Quast et al., 2013) for the microcosm study and with the Greengenes version of the Gold database for the field study. Taxonomy was assigned to representative sequences using the SILVA 111 release database (Quast et al., 2013) for the microcosm study and with the Greengenes 13.8 release for the field study. Representative sequences were aligned using the SILVA 111 core reference alignment using PyNAST (Caporaso et al., 2010b). In total, 4,414,311 (microcosm study) and 5,159,912 (field study) sequences in 150 samples each were analyzed using QIIME.
Several α-diversity indices, as well as indices depicting the population structure, were calculated with the QIIME pipeline (Caporaso et al., 2010a) based on the rarefied OTU table at a depth of 5,000 sequences/sample (observed species, PD whole tree, chao1 and simpson reciprocal). The indices were compared using analyses of variances (ANOVAs) in a factorial design followed by the Tukey HSD test (p < 0.05). Potential bacterial composition differences at phylum level were visualized using a comparative bar chart. Phyla whose total relative abundances in all samples were <0.05% were excluded from the analysis. Principal coordinate analysis (PCoA) of OTU weighted unifrac distance matrices, based on Bray-Curtis dissimilarity matrices, was performed at OTU (species) level and plotted (Hamady et al., 2010). ANOSIM analyses were performed in a factorial design to identify potentially significant differences at the community level between treatments of the same sampling day. Significantly variating OTUs were identified using the R package “pamR.” ANOVA on OTUs were followed by the Bonferroni test grouping treatments per each OTU.
Results
Evaluation of the Impact of Pesticides on the Soil Bacterial Community With Phylochip
In the microcosm study, a soil bacterial community composed of bacteria belonging to 414 families (related to the 63 bacterial phyla present on the PhyloChip) was detected. From these, six phyla were the most dominant in the microcosm study (Proteobacteria, Firmicutes, Actinobacteria, Bacteroidetes, Chloroflexi, and Acidobacteria) constituting approximately 90% of the total bacterial abundance in each sample (Supplementary Figure 1). PCoA and Permanova analyses indicated that there were no significant differences at OTU level in response to pesticide treatment, pesticide dose and time, (p = 0.88, p = 0.88, and p = 0.15, respectively).
In the field study, a bacterial community composed of members belonging to 402 families were detected via PhyloChip analysis. The six most dominant phyla in the field study were Proteobacteria, Actinobacteria, Firmicutes, Cyanobacteria, Bacteroidetes, and Acidobacteria, which made up around 95% of the bacterial abundance in each sample (Supplementary Figure 2). Interestingly, PCoA and Permanova analyses of the field samples at OTU level indicated a significant difference in the presence/absence of OTUs in response to the pesticide treatment (p = 0.013) but not for the pesticide dose (p = 0.458) or time 70 d after treatment (Figure 1B), while no significant differences were detected at 0 d after treatment in response to pesticide treatment and pesticide dose (p = 0.201, and p = 0.124, respectively) (Figure 1A). Interestingly, by performing pairwise tests between the pesticide treatments and the controls, significant differences were observed between the control and all three pesticides. As observed for the microcosm experiment, neither the pesticides dose nor the time had an effect on the bacterial composition when all samples of 0 d and 70 d were compared (p = 0.257, and p = 0.045, respectively).
Figure 1. Field study—PcoA on Bray Curtis similarity matrices based on presence and absence of bacterial OTUs in controls (0x dose) and samples treated with CHL, IPU, or TCZ (1, 2x or 5x recommended agronomical doses) and collected at (A) 0 days or (B) 70 days after treatment. For each treatment values of each of the three replicates are shown. Variance is expressed in % on each axis.
Evaluation of the Impact of Pesticides on the Diversity and Composition of the Soil Bacterial Community by 16S rDNA Amplicon Next-Generation Sequencing
The Impact of Pesticides on the α-Diversity of the Soil Bacterial Community
Analysis of the α-diversity indices “observed species,” “PD whole tree,” “chao1,” and “Simpson reciprocal” indicated that the soil bacterial α-diversity was significantly affected by time in the microcosm study (Table 1), whereas time effects on the α-diversity of soil bacteria were less marked in the field study (Table 2). This observation was further confirmed by ANOSIMs of PCoAs grouping only the untreated control samples of all time points in the microcosm study (p = 0.001, Supplementary Figure 6, Table 3) and in the field study (p = 0.003, Supplementary Figure 10, Table 3), respectively.
Table 1. Microcosm study—Bacterial α-diversity indices “observed species, PD whole tree, chao1 and Simpson reciprocal” at various time points (0, 7, 42, 56, and 100 days after treatment) in soil exposed to each of the three studied pesticides (CHL, IPU, or TCZ) applied at different doses (0x, 1x, 2x, or 10x the recommended agronomical dose).
Table 2. Field study—Bacterial α-diversity indices “observed species, PD whole tree, chao1 and Simpson reciprocal” at various time points (0, 14, 35, 70, and 105 days after treatment) in soil exposed to each of the three studied pesticides (CHL, IPU, or TCZ) applied at different doses (0x, 1x, 2x, or 5x the recommended agronomical dose).
Table 3. Microcosm and field study—ANOSIM of the PCoA ordinations of OTU weighted unifrac distance matrices for untreated (control) and treated [CHL, IPU or TCZ at 1x, 2x, or 10x (microcosms) or 1x, 2x, or 5x (field) doses] soil samples at different time points [0, 7, 42, 56, and 100 days (microcosms) or 0, 14, 35, 70 and 105 days (field)].
In contrast, regardless of the applied pesticide or dose, no significant differences in any of the diversity indices between control (0x) and pesticide-treated (1x, 2x, 10x) soil microcosms were recorded at each time point at both experimental scales (microcosm and field study), indicating that none of the three pesticides significantly modified the α-diversity of the soil bacterial community (Tables 1, 2).
The Impact of Pesticides on the β-Diversity of the Soil Bacterial Community
Based on the amplicon NGS analysis, Proteobacteria, Acidobacteria, Chloroflexi, Bacteroidetes, Actinobacteria, and Gemmatimonadetes constituted the six most abundant bacterial phyla in the soil samples from the microcosm and the field experiment representing together up to 85% of the abundance of all bacterial phyla. Further analysis showed that neither the applied pesticide and dose nor the time significantly affected the abundance of the most dominant bacterial phyla in the soil (Figures 2, 3).
Figure 2. Microcosm study—The relative abundance (%) of bacterial phyla in the untreated (control) and pesticide-treated (CHL, IPU, or TCZ at 10x dose) soil samples at different time points (0, 7, 42, 56, and 100 days after treatment). 16S rDNA sequences that could not be assigned to a phylum were grouped as “unassigned”.
Figure 3. Field study—The relative abundance (%) of bacterial phyla in the untreated (control) and pesticide-treated (CHL, IPU, or TCZ at 5x dose) soil samples at different time points (0, 14, 35, 70, and 105 days after treatment). Bacterial phyla whose abundance sum of all samples was <0.05% were grouped as “others.” 16S rDNA sequences that could not be assigned to a phylum were grouped as “unassigned”.
A more extensive analysis of the 16S rDNA amplicon sequences was performed by ANOSIMs (Table 3) on PCoAs (Supplementary Figures 3–10) using OTU weighted unifrac distance matrices. This analysis verified the lack of highly significant pesticide-induced differences in the composition of the soil bacterial community. However, at certain time points in the microcosm and the field study, transient pesticide-driven effects were identified. In the microcosm study, a significantly modified (p < 0.01) bacterial community was observed at 56 days in the CHL-treated samples (all three doses grouped together) compared to the corresponding untreated controls (Table 3, Supplementary Figure 3). Similarly, in the field study, TCZ induced significant transient changes (p < 0.01) in the composition of the soil bacterial community at 35 days (Table 3, Supplementary Figure 9). No significant changes in the soil bacterial β-diversity were observed in response to IPU-treatment in both microcosm and field experiments.
Analysis of the detected OTUs by the R package “pamR” was employed to identify the OTUs driving the significant changes detected by ANOSIM in the β-diversity of bacteria in the CHL- or TCZ-treated samples collected at 35 or 56 days from the field or microcosm study, respectively. In the microcosm study, eight OTUs were responsible for the differences between CHL-treated samples and untreated samples at 56 days: unknown species of the (i) genus Aeromicrobium, (ii) class Acidobacteria, (iii) genus Flexibacter, (iv)order Chloroacidobacterium, (v) order Sphingobacteriales, and three unknown species of the phylum Chloroflexi (Supplementary Figure 11). In the field study, two OTUs were found to be responsible for differences between TCZ-treated and control samples at 35 days: unknown species of the (i) class Acidobacteria and (ii) the order Acidimicrobiales (Supplementary Figure 12).
Discussion
The study evaluated at both lab and field scale the ecotoxicological impact of three commonly used pesticides, CHL, IPU, or TCZ, on the soil bacterial diversity initially using PhyloChip, a high density phylogenetic microarray, and then further by Illumina next-generation sequencing of 16S rDNA amplicons. Both approaches revealed that Proteobacteria, Acidobacteria, Bacteroidetes, and Actinobacteria dominated the total bacterial community in both the microcosm and the field study. This is in accordance with various previous studies which also reported the dominance of these phyla in a range of soils (Janssen, 2006; Acosta-Martinez et al., 2008; Wessen et al., 2010; Petric et al., 2011a; Merlin et al., 2015; Romdhane et al., 2016). The abundance of these phyla remained relatively constant throughout the microcosm and the field study.
Time Effects More Distinctive in the Microcosm Study Than in the Field Study
Diversity indices in the control samples (not exposed to pesticides) significantly decreased with time in the microcosm study. This decrease trend in the bacterial diversity was not seen in the field study, where only a few time points of some diversity indices were affected by time without showing a clear pattern of decrease. This is in accordance with one of the limitations of long-term soil microcosm experiments causing changes in microbial abundance, composition, diversity and activity due to the lack of nutrient inputs (Edwards et al., 1997). In our study, the duration of the microcosm experiment was extended to 100 days in order to accommodate the long persistence of CHL, IPU and TCZ (10x dose rate) in the given soil with DT90s of 96, 85 and 316 days, respectively (Papadopoulou et al., 2016).
The changes in α-diversity observed in the soil samples collected from the untreated plots in the field experiment could not be assigned to shortage of nutrients but most likely to the climatic conditions which were characterized by long cold periods (mean daily temperature <4°C between 15 and 70 days) and precipitation events at days 3 to 11 which might be responsible for the higher α-diversity values observed at 14 days. PCoA and ANOSIM analyses of the data from amplicon sequencing data at OTU level confirmed the significant effect of time on bacterial diversity in both experimental scales, with a more distinctive time effect on the β-diversity in the microcosm study (p = 0.001) than in the field study (p = 0.003).
Pesticide Effects on the Soil Bacterial Diversity
Although no systematic changes in the diversity and composition of the soil bacterial community were observed in response to CHL exposure, amplicon sequencing analysis identified a slight but significant change in the β-diversity of soil bacteria at 56 days in the microcosm experiment. Interestingly, this coincided with the formation of the maximum levels of 3,5,6-trichloropyridinol (TCP), the main hydrolysis product of CHL, in the soil which received the 10x dose (Papadopoulou et al., 2016). In contrast, in the field experiment, where CHL rapidly dissipated and high TCP concentrations were not observed (probably due to TCP leaching toward groundwater out of the examined soil layer), both the α- and β-diversity was not significantly altered when examined by NGS (although PhyloChip analysis proposed a slight but significant CHL effect at 70 days in the field study). One could therefore hypothesize that (i) changes in the bacterial diversity recorded at 56 days in the microcosm experiment in response to CHL exposure may not only be due to CHL (that had partially dissipated by that time) but could be also assigned to the formation of TCP and that (ii) resilience resulting in the recovery of the bacterial diversity by the end of the study may be attributed to the rapid dissipation of TCP and the gradual dissipation of CHL. This is in accordance with a number of studies reporting the toxicity of TCP to (micro)organisms (Feng et al., 1997; Cáceres et al., 2007). Although the review report of the European Commission (EC, 2005) and other studies (Singh et al., 2002) concluded that CHL does not impose unacceptable effects on soil microorganisms, our study reports a slight but significant ecotoxicological effect of CHL on the soil bacterial β-diversity only in the microcosm experiment (when analyzed by NGS), probably because CHL and TCP cannot dissipate out of the tested soil layer by transport processes (leaching, run-off, volatilization), as it can be the case in field experiments where CHL was not found to modify the bacterial composition when analyzed by NGS. The effect of CHL was only transient and the soil bacterial β-diversity recovered. The biological importance (as compared to other environmental stresses such as extreme weather conditions) and the consequences of this transient modification of the diversity in response to CHL exposure on the functions supported by the bacterial community remain unknown. Our findings are in the line with previous studies which also observed an impact of CHL on the diversity of soil bacteria using molecular fingerprinting (Gilani et al., 2010; Gupta et al., 2013) or PLFAs (Pozo et al., 1995). However, most of the studies did not investigate the parallel formation of pesticide transformation products which might have a secondary effect on soil microorganisms, and also relied on methods which do not have the discrimination capacity of PhyloChip and mostly amplicon sequencing.
IPU exposure did not affect the soil bacterial diversity, as determined by amplicon sequencing, in both microcosms and field experiments. Similar results were obtained by the PhyloChip with the only exception of the slight but significant effect of IPU on the presence/absence of OTUs at 70 d in the field experiment. Overall, our findings are in agreement with the conclusion of EFSA (2015b) “that IPU shows no unacceptable effect on carbon and nitrogen mineralization” and with other studies reporting small to moderate (sometimes temporary) effects of IPU on soil microorganisms (mainly based on decreasing microbial activity and biomass) (Tag-El-Din, 1982; Schuster and Schröder, 1990; Harden et al., 1993; Kuriyal and Pandey, 1994).
In the field experiment, TCZ exposure induced slight but significant changes in the β-diversity of the bacterial community after 35 days (NGS outcome) and 70 days (PhyloChip outcome) of exposure in the field study. TCZ was rapidly dissipated in the field experiment (60% dissipation after 35 days) (Papadopoulou et al., 2016). In addition, an important number and variety of TCZ transformation products, including triazole dead-end transformation products possibly interacting with the hormone regulation network of non-target organisms (Shalini et al., 2011; Rieke et al., 2014), were detected but not quantified (due to the absence of appropriate reference standards) in the samples of the field experiment (treated with the 5x agronomical dose) (Storck et al., 2016). Therefore, one could propose that observed changes in the bacterial composition are not only due to exposure to TCZ but also to its transformation products. Nevertheless, although the late resilience of the bacterial composition was not addressed by PhyloChip, one could observe by NGS results that the bacterial β-diversity was resilient suggesting the recovery of the bacterial community. Although the evaluations of EFSA (2008, 2014) for TCZ concluded that the fungicide does not entail an unacceptable risk for soil microorganisms (based on carbon and nitrogen mineralization tests), our study reports slight but significant and transient effect of TCZ on the composition of the bacterial community in the field experiment. This is in agreement with previous studies reporting TCZ effects on soil microorganisms, evaluated by analysis of soil enzyme activities (Anuradha et al., 2016), microbial activity and biomass (Cycon et al., 2006; Bending et al., 2007; Munoz-Leoz et al., 2011; Wang et al., 2016), or DGGE fingerprinting (Ferreira et al., 2009).
Our study showed that small and transient changes in the diversity of soil bacteria can be detected for two (CHL and TCZ) of the three tested pesticides by 16S rDNA amplicon NGS, while the use of PhyloChip indicated significant changes in response to all three tested pesticides only for the field study. The main advantage of the PhyloChip lies on its great sensitivity for rare species (Nikolaki and Tsiamis, 2013). Indeed, one could estimate that 500 ng of 16S rDNA PCR product comprised more than 600 billion sequences enabling the detection of not only dominant but also of rare amplicons (Katsaveli et al., 2012; Nikolaki and Tsiamis, 2013). On the contrary, the main disadvantages of PhyloChip are that only known phylotypes can be detected and that signals obtained for dominant phylotypes can be saturated leading to underestimation of their abundance.
Contrariwise to Romdhane et al. (2016) reporting big shifts in the bacterial abundance and diversity exposed to leptospermone, a triketone herbicide of natural origin, only slight effects were observed in our study. This may be due to the fact that, in contrary to a natural triketone for which the target (4-hydroxyphenylpyruvate dioxygenase, HPPD) is present in non-target organisms including soil microorganisms, the targets of the three pesticides evaluated in our study are (to our best knowledge) not found in soil bacteria (IPU targets photosystem II, CHL blocks acetylcholine neurotransmission, and TCZ blocks sterol 14α-demethylase and impedes any additional modes of action on sterol biosynthesis). Thus, no direct effects of the three tested pesticides on soil bacteria were expected. Nonetheless, the transformation products of the three tested pesticides that are known for their toxicity [namely 4-isopropyl-aniline (4-IA) for IPU, TCP for CHL, and triazole transformation products for TCZ] may directly impact soil microorganisms. In addition, indirect effects of the three pesticides can be expected notably for TCZ because of close interactions between fungal and bacterial species in soil.
In any cases, the interpretation of the biological significance of the ecotoxicological effects of the tested pesticides on the composition of the soil bacterial community remains difficult in the absence of data on the normal operating range (NOR) for each studied bioindicator (such as bacterial diversity), making it challenging to define specific protection goals to protect bacterial diversity and supported soil ecosystem functions (Bell et al., 2005).
Conclusions
Our study provides a first comprehensive lab-to-field assessment of the ecotoxocity of the three pesticides CHL, IPU, and TCZ on the soil bacterial diversity using advanced molecular methods such as PhyloChip and 16S rDNA amplicon NGS. Although the EFSA concluded that all three pesticides do not induce unacceptable changes on the soil microbial activity, one could observe that they induced minor (but significant) and transient changes in the composition of the soil bacterial community. These modern methods seem to have a good potential for the assessment of the toxicity of pesticides on the soil microbial diversity. More research is needed to define the NOR of the soil bacterial diversity in order to biologically interpret the importance of changes observed in response to pesticide exposure.
Author Contributions
VS and SN contribute to the production of the results, their exploitation and writing. CP, CC, GP, CB, PK, and EP contributed to the setup of field and lab experiments. In addition, GP and CB were in charge of soil DNA extraction. AnS and FF were in charge of pesticide residues extraction and quantification. AyS, MD-L, CM and GT were in charge of the exploitation of NGS results. DK, OS, GT and FM-L coordinated and supervised all the co-authors within the framework of the EU project LTH. They contribute actively to the preparation of this manuscript. All co-authors agreed to publish this paper to which they contribute.
Conflict of Interest Statement
CB, OS, and CM were employed by ENOVEO srl. AS, GP, and FF were employed by AEIFORIA srl.
The remaining authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Acknowledgments
The Ph.D. studies of VS performed at the Doctoral School of Environment and Health (ED E2S) at the University of Burgundy were funded by the French Ministry of Education and Research (grant number 2013-31). This work was done within the framework of the FP7-PEOPLE-2012-IAPP (Industry-Academia Partnerships and Pathways) Marie-Curie project Love-to-Hate funded by the European Commission (grant number 324349).
Supplementary Material
The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmicb.2018.01412/full#supplementary-material
References
Acosta-Martinez, V., Dowd, S., Sun, Y., and Allen, V. (2008). Tag-encoded pyrosequencing analysis of bacterial diversity in a single soil type as affected by management and land use. Soil Biol. Biochem. 40, 2762–2770. doi: 10.1016/j.soilbio.2008.07.022
Anuradha, B., Rekhapadmini, A., and Rangaswamy, V. (2016). Influence of tebuconazole and copper hydroxide on phosphatase and urease activities in red sandy loam and black clay soils. 3 Biotech 6:78. doi: 10.1007/s13205-016-0367-0
Barbash, J. E. (2003). “The geochemistry of pesticides,” in Treatise on Geochemistry (Washington, DC: Elsevier), 9.
Bell, T., Newman, J. A., Silverman, B. W., Turner, S. L., and Lilley, A. K. (2005). The contribution of species richness and composition to bacterial services. Nature 436, 1157–1160. doi: 10.1038/nature03891
Bending, G. D., Rodríguez-Cruz, M. S., and Lincoln, S. D. (2007). Fungicide impacts on microbial communities in soils with contrasting management histories. Chemosphere 69, 82–88. doi: 10.1016/j.chemosphere.2007.04.042
Blum, W. E. H. (2006). Soil resources - the basis of human society and the environment. Bodenkultur 57, 197–202.
Bondeau, A., Smith, C. M., Zaehle, S., Schapfhoff, S., Lucht, W., Cramer, W., et al. (2007). Modelling the role of agriculture for the 20th century global terrestrial carbon balance. Glob. Change Biol. 13, 679–706. doi: 10.1111/j.1365-2486.2006.01305.x
Börner, H., Schlüter, K., and Aumann, J. (2009). Pflanzenkrankheiten und Pflanzenschutz. Berlin: Springer.
Cáceres, T., He, W., Naidu, R., and Megharaj, M. (2007). Toxicity of chlorpyrifos and TCP alone and in combination to Daphnia carinata: the influence of microbial degradation in natural water. Wat. Res. 41, 4497–4503. doi: 10.1016/j.watres.2007.06.025
Caporaso, J. G., Bittinger, K., Bushman, F. D., DeSantis, T. Z., Andersen, G. L., and Knight, R. (2010b). PyNAST: a flexible tool for aligning sequences to a template alignment. Bioinformatics 26, 266–267. doi: 10.1093/bioinformatics/btp636
Caporaso, J. G., Kuczynski, J., Stombaugh, J., Bittinger, K., Bushman, F. D., and Costello, E. K. (2010a). QIIME allows analysis of high-throughput community sequencing data. Nat. Methods 7, 335–336. doi: 10.1038/nmeth.f.303
Chu, X., Fang, H., Pan, X., Wang, X., Shan, M., Feng, B., et al. (2008). Degradation of chlorpyrifos alone and in combination with chlorothalonil and their effects on soil microbial populations. J. Environ. Sci. 20, 464–469. doi: 10.1016/S1001-0742(08)62080-X
Cycon, M., Piotrowska-Seget, Z., and Kaczynska, A. (2006). Microbiological characteristics of a sandy loam soil exposed to tebuconazole and λ-cyhalothrin under laboratory conditions. Ecotoxicology 15, 639–646. doi: 10.1007/s10646-006-0099-8
Dutta, M., Sardar, D., Pal, R., and Kole, R. K. (2010). Effect of chlorpyrifos on microbial biomass and activities in tropical clay loam soil. Environ. Monit. Assess. 160, 385–391. doi: 10.1007/s10661-008-0702-y
EC (2005). Review Report for the Active Substance Chlorpyrifos. Finalised in the Standing Committee on the Food Chain and Animal Health at its meeting on 3 June 2005 in view of the inclusion of chlorpyrifos in Annex I of the Directive 91/414/EEC. SANCO/3059/99 - rev. 1.5.
Edgar, R. C. (2010). Search and clustering orders of magnitude faster than BLAST. Bioinformatics 26, 2460–2461. doi: 10.1093/bioinformatics/btq461
Edgar, R. C., Haas, B. J., Clemente, J. C., Quince, C., and Knight, R. (2011). UCHIME improves sensitivity and speed of chimera detection. Bioinformatics 27, 2194–2200. doi: 10.1093/bioinformatics/btr381
Edwards, C. A., Subler, S., Parmelee, R., Knacker, T., and Pokarzhevskij, A. A. (1997). Use of soil microcosms in assessing the effect of pesticides on soil ecosystems. Environ. Behav. Crop Protect. Chem. 28, 435–451.
EFSA (2014). Conclusion on the peer review of the pesticide risk assessment of the active substance tebuconazole. EFSA J. 12:3485. doi: 10.2903/j.efsa.2014.3485
EFSA (2015a). Conclusion on the peer review of the pesticide risk assessment of the active substance isoproturon. EFSA J. 13:4206. doi: 10.2903/j.efsa.2015.4206
EFSA (2015b). Conclusion on the peer review of the pesticide risk assessment of the active substance glyphosate. EFSA J. 13:4302. doi: 10.2903/j.efsa.2015.4302
Eisenhauer, N., Klier, M., Partsch, S., Sabais, A. C. W., Scherber, C., Weisser, W. W., et al. (2009). No interactive effects of pesticides and plant diversity on soil microbial biomass and respiration. Appl. Soil Ecol. 42, 31–36. doi: 10.1016/j.apsoil.2009.01.005
Fang, H., Yu, Y., Chu, X., Wang, X., Yang, X., and Yu, J. (2009). Degradation of chlorpyrifos in laboratory soil and its impact on soil microbial functional diversity. J. Environ. Sci. 21, 380–386. doi: 10.1016/S1001-0742(08)62280-9
Feld, L., Hjelmsø, M. H., Nielsen, M. S., Jacobsen, A. D., Rønn, R., Ekelund, F., et al. (2015). Pesticide side effects in an agricultural soil ecosystem as measured by amoA expression quantification and bacterial diversity changes. PLoS ONE 10:e0126080. doi: 10.1371/journal.pone.0126080
Feng, Y., Racke, K. D., and Bollag, J. (1997). Isolation and characterization of a chlorinated-pyridinol-degrading bacterium Appl Environ. Microbiol. 63, 4096–4098.
Ferreira, E. P. B., Dusi, A. N., Costa, J. R., Xavier, G. R., and Rumjanek, N. G. (2009). Assessing insecticide and fungicide effects on the culturable soil bacterial community by analyses of variance of their DGGE fingerprinting data. Eur. J. Soil Biol. 45, 466–472. doi: 10.1016/j.ejsobi.2009.07.003
Gardi, C., Jeffery, S., and Saltelli, A. (2013). An estimate of potential threats levels to soil biodiversity in EU. Glob. Change Biol. 19, 1538–1548. doi: 10.1111/gcb.12159
Gilani, S. T. S., Ageen, M., Shah, H., and Raza, S. (2010). Chlorpyrifos degradation in soil and its effect on soil microorganisms. J. Anim. Plant Sci. 20, 99–102.
Graham, E. B., Knelman, J. E., Schindlbacher, A., Siciliano, S., Breulmann, M., Yannarell, A., et al. (2016). Microbes as engines of ecosystem function: when does community structure enhance predictions of ecosystem processes? Front. Microbiol. 7:214. doi: 10.3389/fmicb.2016.00214
Gupta, S., Gupta, R., and Sharma, S. (2013). Impact of chemical- and bio-pesticides on bacterial diversity in rhizosphere of Vignaradiata. Ecotoxicology 22, 1479–1489. doi: 10.1007/s10646-013-1134-1
Hamady, M., Lozupone, C., and Knight, R. (2010). Fast UniFrac: facilitating high-throughput phylogenetic analyses of microbial communities including analysis of pyrosequencing and PhyloChip data. ISME J. 4, 17–27. doi: 10.1038/ismej.2009.97
Harden, T., Joergensen, R. G., Meyer, B., and Wolters, V. (1993). Soil microbial biomass estimated by fumigation-extraction and substrate-induced respiration in two pesticide-treated soils. Soil Biol. Biochem. 25, 679–683. doi: 10.1016/0038-0717(93)90107-M
Hartley, J. L., and Bowen, H. (2003). PEG precipitation for selective removal of small DNA fragments. Focus 25:18.
Hayes, W. J., and Laws, E. R. (1991). “Handbook of pesticide toxicology,” in Classes of Pesticides, Vol. 3, ed R. Krieger (San Diego, CA: Academic Press Inc.), 1571.
Haygarth, P., and Ritz, K. (2009). The future of soils and land use in the UK: soil systems for the provision of land-based ecosystem services. Land Use Policy 26, 187–197. doi: 10.1016/j.landusepol.2009.09.016
Hillel, D. (2009). The mission of soil science in a changing world. J. Plant Nutr. Soil Sci. 172, 5–9. doi: 10.1002/jpln.200800333
Imfeld, G., and Vuilleumier, S. (2012). Measuring the effects of pesticides on bacterial communities in soil: A critical review. Eur. J. Soil Biol. 49, 22–30. doi: 10.1016/j.ejsobi.2011.11.010
ISO 10381-6 (2009). Soil Quality - Sampling - Part 6: Guidance on the Collection, Handling and Storage of Soil Under Aerobic Conditions for the Assessment of Microbiological Processes Biomass and Diversity in the Laboratory.
Janssen, P. H. (2006). Identifying the dominant soil bacterial taxa in libraries of 16S rRNA and 16S rRNA genes. Appl. Environ. Microbiol. 72, 1719–1728. doi: 10.1128/AEM.72.3.1719-1728.2006
Jeffery, S., Gardi, C., Jones, A., Montanarella, L., Marmo, L., Miko, L., et al. (2010). European Atlas of Soil Biodiversity. European Commission; Publications Office of the European Union.
Jones, A., Stolbovoy, V., Rusco, E., Gentile, A.-R., Gardi, C., Marechal, B., et al. (2009). Climate change in Europe. 2. Impact on soil. A review. Agronomy for Sustainable Development 29, 423–432. doi: 10.1051/agro:2008067
Joshi, N., and Fass, J. (2011). Sickle: A Sliding-Window, Adaptive, Quality-Based Trimming Tool for FastQ Files. Available online at http://github.com/najoshi/sickle (Accessed February 1, 2016).
Katsaveli, K., Vayenas, D., Tsiamis, G., and Bourtzis, K. (2012). Bacterialdiversity in Cr(VI) and Cr(III)-contaminated industrialwastewaters. Extremophiles 16, 285–296. doi: 10.1007/s00792-012-0429-0
Klindworth, A., Pruesse, E., Schweer, T., Peplies, J., Quast, C., Horn, M., et al. (2012). Evaluation of general 16S ribosomal RNA gene PCR primers for classical and next-generation sequencing-based diversity studies. Nucleic Acids Res. 41:e1. doi: 10.1093/nar/gks808
Kreader, C. A. (1996). Relief of amplification inhibition in PCR with bovine serum albumin or T4 gene 32 protein. Appl. Environ. Microbiol. 62, 1102–1106.
Kuriyal, S. K., and Pandey, B. N. (1994). Studies of the effects of isoproturon on soil microflora. J. Indian Bot. Soc. 73, 47–49.
Lal, R. (2008). Soils and sustainable agriculture. A review. Agron. Sustain. Dev. 28, 57–64. doi: 10.1051/agro:2007025
Lane, D. (1991). Nucleic Acid Techniques in Bacterial Systematics. Michigan: Wiley, Michigan University Press.
Lichtfouse, E., Navarette, M., Debaeke, P., Souchère, V., and Alberola, C. (2009). Sustainable Agriculture, 1st Edn. Amsterdam: Springer.
Looser, R., Froescheis, O., Cailliet, G. M., Jarman, W. M., and Ballschmiter, K. (2000). The deep-sea as a final global sink of semivolatile persistent organic pollutants? Part II: organochlorine pesticides in surface and deep-sea dwelling fish of the North and South Atlantic and the Monterey Bay Canyon (California). Chemosphere 40, 661–670. doi: 10.1016/S0045-6535(99)00462-2
Martin, M. (2011). Cutadapt removes adapter sequences from high-throughput sequencing reads. EMBnet J. 17, 10–12. doi: 10.14806/ej.17.1.200
Martin-Laurent, F., Philippot, L., Hallet, S., Chaussod, R., Germon, J. C., Soulas, G., et al. (2001). DNA extraction from soils: old bias for new microbial diversity analysis methods. Appl. Environ. Microbiol. 67, 2354–2359. doi: 10.1128/AEM.67.5.2354-2359.2001
Masella, A. P., Bartram, A. K., Truszkowski, J. M., Brown, D. G., and Neufeld, J. D. (2012). PANDAseq: paired-end assembler for Illumina sequences. BMC Bioinformatics 13:31. doi: 10.1186/1471-2105-13-31
Merlin, C., Devers, M., Béguet, J., Boggio, B., Rouard, N., and Martin-Laurent, F. (2015). Evaluation of the ecotoxicological impact of the organochlorine chlordecone on soil microbial community structure, abundance, and function. Environ. Sci. Pollut. Res. 23, 4185–4198. doi: 10.1007/s11356-015-4758-2
Munoz-Leoz, B., Ruiz-Romera, E., Antigëdad, I., and Garbisu, C. (2011). Tebuconazole application decreases soil microbial biomass and activity. Soil Biol. Biochem. 43, 2146–2183. doi: 10.1016/j.soilbio.2011.07.001
Newman, M. M., Hoilett, N., Lorenz, N., Dick, R. P., Liles, M. R., Ramsier, C., et al. (2016). Glyphosate effects on soil rhizosphere-associated bacterial communities. Sci. Total Environ. 543, 155–160. doi: 10.1016/j.scitotenv.2015.11.008
Nikolaki, S., and Tsiamis, G. (2013). Microbial diversity in the era of omic technologies. BioMed Res. Int. 2013:958719. doi: 10.1155/2013/958719
Oerke, E.-C., Dehne, H. W., Schönebeck, F., and Weber, A. (eds.). (1994). “Crop production and crop protection,” in Estimated Losses in Major Food and Cash Crops (Amsterdam: Elsevier Science). doi: 10.1016/C2009-0-00683-7
Papadopoulou, E. S., Genitsaris, S., Omirou, M., Perruchon, C., Stamatopoulou, A., Ioannides, I., et al. (2018). Bioaugmentation of thiabendazole-contaminated soils from a wastewater disposal site: factors driving the efficacy of this strategy and the diversity of the indigenous soil bacterial community. Environ. Pollut. 233, 16–25. doi: 10.1016/j.envpol.2017.10.021
Papadopoulou, E. S., Karas, P., Nikolaki, S., Storck, V., Ferrari, F., Trevisan, M., et al. (2016). Dissipation and adsorption of isoproturon, tebuconazole, chlorpyrifos and their main transformation products under laboratory and field conditions. Sci. Total Environ. 569–570, 86–96. doi: 10.1016/j.scitotenv.2016.06.133
Perrin-Ganier, C., Schiavon, F., Morel, J.-L., and Schiavon, M. (2001). Effect of sludge-amendment or nutrient addition on the biodegradation of the herbicide isoproturon in soil. Chemosphere 44, 887–892. doi: 10.1016/S0045-6535(00)00283-6
Petric, I., Bru, D., Udikovic-Kolic, N., Hrsak, D., Philippot, L., and Martin-Laurent, F. (2011a). Evidence for shifts in the structure and abundance of the microbial community in a long-term PCB-contaminated soil under bioremediation. J. Hazard. Mater. 195, 254–260. doi: 10.1016/j.jhazmat.2011.08.036
Petric, I., Philippot, L., Abbate, C., Bispo, A., Chesnot, T., Hallin, S., et al. (2011b). Inter-laboratory evaluation of the ISO standard 11063 “Soil quality - Method to directly extract DNA from soil samples”. J. Microbiol. Methods 84, 454–460. doi: 10.1016/j.mimet.2011.01.016
Pimentel, D. (1995). Amounts of pesticides reaching target pests: Environmental impacts and ethics. J. Agric. Environ. Ethics 8, 17–29. doi: 10.1007/BF02286399
Pozo, C., Martinez-Toledo, V., Salmeron, V., Rodelas, B., and Gonzalez-Lopez, J. (1995). Effect of chlorpyrifos on soil microbial activity. Environ. Toxicol. Chem. 14, 187–192. doi: 10.1002/etc.5620140201
Puglisi, E. (2012). Response of Microbial Organisms (Aquatic and Terrestrial) to Pesticides. EFSA Supporting Publications.
Quast, C., Pruesse, E., Yilmaz, P., Gerken, J., Schweer, T., and Yarza, P. (2013). The SILVA ribosomal RNA gene database project: improved data processing and web-based tools. Nucleic Acids Res. 41, 590–596. doi: 10.1093/nar/gks1219
Rattan, L. (2015). Restoring soil quality to mitigate soil degradation. Sustainability 7, 5875–5895. doi: 10.3390/su7055875
Rieke, S., Koehn, S., Hirsch-Ernst, K., Pfeil, R., Kneuer, C., and Marx-Stoelting, P. (2014). Combination effects of (tri)azole fungicides on hormone production and xenobiotic metabolism in a human placental cell line. Int. J. Environ. Res. Public Health 11, 9660–9679. doi: 10.3390/ijerph110909660
Romdhane, S., Devers-Lamrani, M., Barthelmebs, L., Calvayrac, C., Bertrand, C., Cooper, J.-F., et al. (2016). Ecotoxicological impact of the bioherbicide Leptospermone on the microbial community of two arable soils. Front. Microbiol. 7:775. doi: 10.3389/fmicb.2016.00775
Schuster, E., and Schröder, D. (1990). Side-effects of sequentially-applied pesticides on non-target soil microorganisms: Field experiments. Soil Biol. Biochem. 22, 367–373. doi: 10.1016/0038-0717(90)90115-G
Shalini, K., Kumar, N., Drabu, S., and Kumar Sharma, P. (2011). Advances in synthetic approach to and antifungal activity of triazoles. Beilstein J. Org. Chem. 7, 668–677. doi: 10.3762/bjoc.7.79
Singh, B. K., Walker, A., and Wright, D. (2002). Degradation of chlorpyrifos, fenamiphos, and chlorothalonil alone and in combination and their effects on soil microbial activity. Environ. Toxicol. Chem. 21, 2600–2605. doi: 10.1002/etc.5620211211
Smalling, K. L., Fellers, G. M., Kleeman, P. M., and Kuivila, K. M. (2013). Accumulation of pesticides in pacific chorus frogs (Pseudacrisregilla) from California's Sierra Nevada Mountains, USA. Environ. Toxicol. Chem. 32, 2026–2034. doi: 10.1002/etc.2308
Srinivasulu, M., and Rangaswamy, V. (2013). Influence of insecticides alone and in combination with fungicides on enzyme activities in soils. Int. J. Environ. Sci. Technol. 10, 341–350. doi: 10.1007/s13762-012-0133-8
Steffan, J. J., Brevik, E. C., Burgess, L. C., and Cerdà, A. (2017). The effect of soil on human health: an overview. Eur. J. Soil Sci. 69, 159–171. doi: 10.1111/ejss.12451
Storck, V., Lucini, L., Mamy, L., Ferrari, F., Papadopoulou, E. S., Nikolaki, S., et al. (2016). Identification and characterization of tebuconazole transformation products in soil by combining suspect screening and molecular typology. Environ. Pollut. 208, 537–545. doi: 10.1016/j.envpol.2015.10.027
Tag-El-Din, A. (1982). Herbicides and respiration flush of soil biomass. Medical Faculty Landbouww University Gent 47, 191–199.
Tsiamis, G., Katsaveli, K., Ntougias, S., Kyrpides, N., Andersen, G., Piceno, Y., et al. (2008). Prokaryotic community profiles at different operational stages of a Greek solar saltern. Res. Microbiol. 159, 609–627. doi: 10.1016/j.resmic.2008.09.007
Wang, C., Wang, F., Zhang, Q., and Liang, W. (2016). Individual and combined effects of tebuconazole and carbendazim on soil microbial activity. Eur. J. Soil Biol. 72, 6–13. doi: 10.1016/j.ejsobi.2015.12.005
Wessen, E., Hallin, S., and Philippot, L. (2010). Differential responses of bacterial and archaeal groups at high taxonomical ranks to soil management. Soil Biol. Biochem. 42, 1759–1765. doi: 10.1016/j.soilbio.2010.06.013
Zhang, W., Jiang, F., and Ou, J. (2011). Global pesticide consumption and pollution: with China as a focus. Proc. Int. Acad. Ecol. Environ. Sci. 1, 125–144. doi: 10.0000/issn-2220-8860-piaees-2011-v1-0012
Keywords: microbial ecotoxicology, pesticide, soil, microorganisms, next-generation sequencing, DNA microarray
Citation: Storck V, Nikolaki S, Perruchon C, Chabanis C, Sacchi A, Pertile G, Baguelin C, Karas PA, Spor A, Devers-Lamrani M, Papadopoulou ES, Sibourg O, Malandain C, Trevisan M, Ferrari F, Karpouzas DG, Tsiamis G and Martin-Laurent F (2018) Lab to Field Assessment of the Ecotoxicological Impact of Chlorpyrifos, Isoproturon, or Tebuconazole on the Diversity and Composition of the Soil Bacterial Community. Front. Microbiol. 9:1412. doi: 10.3389/fmicb.2018.01412
Received: 31 August 2017; Accepted: 08 June 2018;
Published: 29 June 2018.
Edited by:
Guillermina Hernandez-Raquet, Institut National de la Recherche Agronomique (INRA), FranceReviewed by:
Roberta Fulthorpe, University of Toronto Scarborough, CanadaBruno B. Castro, University of Minho, Portugal
Copyright © 2018 Storck, Nikolaki, Perruchon, Chabanis, Sacchi, Pertile, Baguelin, Karas, Spor, Devers-Lamrani, Papadopoulou, Sibourg, Malandain, Trevisan, Ferrari, Karpouzas, Tsiamis and Martin-Laurent. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: George Tsiamis, gtsiamis@upatras.gr
Fabrice Martin-Laurent, fabrice.martin@inra.fr
†First authors.