A Metagenomic Study Highlights Phylogenetic Proximity of Quorum-Quenching and Xenobiotic-Degrading Amidases of the AS-Family

Abstract

Quorum-sensing (QS) signals of the N-acylhomoserine lactone (NAHL) class are cleaved by quorum-quenching enzymes, collectively named NAHLases. Here, functional metagenomics allowed the discovery of a novel bacterial NAHLase in a rhizosphere that was treated with γ-caprolactone. As revealed by rrs-DGGE and rrs-pyrosequencing, this treatment increased the percentage of the NAHL-degrading bacteria and strongly biased the structure of the bacterial community, among which Azospirillum dominated. Among the 29 760 fosmids of the metagenomic library, a single one was detected that expressed the qsdB gene conferring NAHL-degradation upon E. coli and decreased QS-regulated virulence in Pectobacterium. Phylogenetic analysis of the 34 orfs of the fosmid suggested that it would belong to an unknown Proteobacterium - probably a γ-proteobacterium. qPCR quantification of the NAHLase-encoding genes attM, qsdA, and qsdB revealed their higher abundance in the γ-caprolactone-treated rhizosphere as compared to an untreated control. The purified QsdB enzyme exhibited amidase activity. QsdB is the first amidase signature (AS) family member exhibiting NAHLase-activity. Point mutations in the AS-family catalytic triad K-S-S abolished the NAHLase activity of QsdB. This study extends the diversity of NAHLases and highlights a common phylogenic origin of AS-family enzymes involved in the degradation of natural compounds, such as NAHLs, and xenobiotics, such as nylon and linuron.

About the Authors:

Mélanie Tannières

Affiliation: Institut des Sciences du Végétal, UPR2355, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France

Amélie Beury-Cirou

Affiliations Institut des Sciences du Végétal, UPR2355, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France, Comité Nord Plants de Pommes de Terre (CNPPT), Semences, Innovation, Protection, Recherche et Environnement (SIPRE), Achicourt, France

Armelle Vigouroux

Affiliation: Laboratoire d’Enzymologie et Biochimie Structurales (LEBS), UPR3082, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France

Samuel Mondy

Affiliation: Institut des Sciences du Végétal, UPR2355, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France

Franck Pellissier

Affiliation: Institut de Chimie des Substances Naturelles, UPR2301, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France

Yves Dessaux

Affiliation: Institut des Sciences du Végétal, UPR2355, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France

Denis Faure

* E-mail: [email protected]

Affiliation: Institut des Sciences du Végétal, UPR2355, Centre National de la Recherche Scientifique, Gif-sur-Yvette, France

Introduction

N-acylhomoserine lactones (NAHLs) are diffusible signals used by many Proteobacteria to correlate gene expression to cell density via a regulatory pathway called quorum-sensing (QS) [1]–[2]. Aside from the enzymes implicated in biosynthesis of NAHLs, others are able to cleave or modify NAHLs, hence to disrupt QS-signaling [3]. They have been identified in bacteria and eukaryotes, and are collectively called NAHLases or quorum-quenching enzymes [4]. According to their structures and enzymatic activities, NAHLases are classified in different families: the lactonases of the Zn-hydrolase family, such as AiiA, AttM/BlcC, AhlD, and QlcA [5]–[9]; the lactonases of the alpha/beta hydrolase fold family such as AidH and AiiM [10]–[11]; the lactonases of the phosphotriesterase family with a single NAHLase member, QsdA [12]; the amidases of the amidohydrolase cluster of the beta-lactam acylases, such as AiiD, PvdQ, AhlM, AiiC, and QuiP [13]–[17]; and the short-chain reductase family with one NAHLase member, the reductase BpiB09 [18]–[19]. The biological role(s) of NAHLases generally remains unclear, except in some firmicutes where they may contribute to the degradation of toxic NAHL-derivatives [20], or in the plant pathogen Agrobacterium tumefaciens in which they slightly modulate QS-functions [21]–[22]. In spite of this lack of information on the biological roles of NAHLases, these enzymes have been used successfully to quench QS-regulated functions. For instance, virulence in the plant pathogen Pectobacterium could be reduced or abolished either via the production of transgenic plants expressing NAHLase-encoding gene [23] or by the selection of bacterial isolates or populations exhibiting NAHLase activities [24]–[27].

Functional metagenomics, which includes screening of environmental-DNA libraries for enzymatic activities or metabolite synthesis, emerged as a powerful approach to explore and exploit the natural biodiversity of microorganisms [28]–[30]. In the investigated metagenomic libraries, the gene(s) of interest may be present at a low frequency. To circumvent this limitation, two main different strategies were implemented, such as (i) the development of cloning vectors that allow a positive selection or an easiest screening and (ii) the enrichment of a microbial community with a biological function of interest under appropriate conditions [31]–[32]. With respect to QS, functional metagenomics extended knowledge of the diversity of both the NAHL-producing and NAHL-degrading enzymes [9], [18]–[19], [33]–[34]. This approach also revealed biosynthesis of novel mimics which activate QS-response [35]. To our knowledge, a biased-environment, enriched in NAHL-degrading bacteria, was never used for discovering novel NAHLases by functional metagenomics.

In this work, γ-caprolactone (GCL), a biodegradable compound structurally related to NAHLs [26], was introduced in the rhizosphere of Solanum tuberosum to increase the percentage of NAHL-degrading bacteria in the plant environment. The GCL-induced bias on bacterial diversity was verified by a combination of rrs-DGGE and rrs-pyrosequencing. Then, genes involved in NAHL-degradation were searched among a library of 29,760 clones using a functional screening. A qsdB-encoded amidase was identified and characterized. This work is a first report of a NAHL-amidase that belongs to the amidase signature (AS) family, members of which are usually involved in the degradation of various xenobiotic compounds, such as nylon and linuron herbicide.

Materials and Methods

Plant Culture

One hundred plants of S. tuberosum var. Allians, which were recovered from cultures performed under axenic conditions, were placed into holes (3 cm space to each other) of batch lids. Each of the non-sterile batches (40×60×8 cm) contained 13 L of the nutritive solution Hydrobloom (Cellmax, UK) with nitrogen at 0.80 g/L, phosphore at 0.56 g/L and potassium at 1.48 g/L as major components. The solution was diluted from a concentrated stock solution (x250) with non-sterile water from the public water system. Planted batches were placed in the greenhouse (Comité Nord Plants de Pomme de Terre, Bretteville-du-Grand-Caux) under natural light at 10–15°C (night) and 25–30°C (day). Treatments (0.4 g/L) with γ-caprolactone (GCL; CAS # 695-06-7; Sigma-Aldrich) were performed at 1 and 28 days. This GCL concentration was selected because of its capacity to stimulate the growth of HAHL-degrading bacteria and its low toxicity to potato plant. A second batch was used as an untreated control. The two batches (untreated, GCL-treated) were simultaneously analyzed at 42-days. GCL concentration in batch and plant tissues was determined by HPLC-MS as previously described [36].

Identification of Bacterial Isolates for NAHL-production and NAHL-degradation

At 42 days, one gram of roots (fresh weight) was suspended in 10 ml of sterile 10 mM MgSO4, diluted and spread onto TSA medium (AES) for counting members of the cultured bacterial community. Three samples were analyzed from each batch. From each of the six samples, thirty TSA-isolates were grown in 96-microwell plates and individually assayed for production of NAHLs with the A. tumefaciens biosensor NT1(pZNLR4) [37]. All the isolates were tested for their capacity to disrupt perception of the QS-signal hexanoylhomoserine lactone (C6HSL) using the biosensor Chromobacterium violaceum CVO26 [38].

rrs-DGGE and rrs-pyrosequencing

PCR-amplicons of the rrs 5′-region between positions 341 and 534 (according to E. coli sequence) were submitted to DGGE analysis (Microbial Insights, Rockford, TN, USA) as previously described [39]. Sequence identification of selected bands was performed using the sequence match facility of the Ribosomal Database Project (http://rdp.cme.msu.edu/) according to NCBI taxonomy (http://www.ncbi.nlm.nih.gov/Taxonomy/).

Experimental procedures for rrs-amplification using 8F- and 357R-derived primers, 454-pyrosequencing (Roche/454 FLX Pyrosequencer at GATC biotech, Konstanz, Germany), and data analysis were previously described by Cirou et al. [36]. The obtained sequences were selected by their length (>150 nt) and their quality score (90% of nucleotides with a quality score>25) using the Greengenes website (http://greengenes.lbl.gov/cgi-bin/nph-index.cgi), and analyzed using the RDP's pyrosequencing pipeline at the Ribosomal Database Project Website (http://pyro.cme.msu.edu/). The reads were aligned with the Infernal Aligner software and complete linkage clustering method were performed to define OTU. Singleton OTUs were kept for the calculation of ecological indexes. CHAO1 and Shannon's Index were calculated at 0.03 and 0.05 distance.

Quantification of qsdA, qsdB and attM Genes by qPCR

At 42 days, 10 L of nutritive solution from the untreated and GCL-treated batches were filtrated through a filter paper to eliminate large plant debris, then centrifuged at 17,600 g. Total genomic DNA was extracted using the DNesay Blood and Tissue Kit (Qiagen) according to the manufacturer’s protocol and quantified with a spectrophotometer (NanoDrop ND1000, Labtech, France). Relative abundance of the qsdA, attM and qsdB genes was estimated by quantitative PCR (qPCR) in 96-well plates using a LightCycler 480 (Roche). Primers were designed for conserved region of the genes qsdA (5′-ACGAGCATGTCTTCGTTCTG and 5′-GGATCGACGA TCGTGCTGAT), attM (5′-TGACATCGGCCGGATCGAAA and 5′-ACGGCGGC AACGCGATTGAA), and qsdB (5′-GAGTGCCCAGGAACTTCACG and 5′-CCTTGAT CAGGAAGGGCACG). Primer couples generated fragments of ca. 150 bp in length. Calibration curves of qsdA, attM, and qsdB were defined with genomic DNA of Rhodococcus erythropolis R138, Agrobacterium tumefaciens C58 and E. coli DH5α harboring the p90H6 fosmid, respectively. Composition of the PCR mix for each sample was as follows: 5 µL of SYBER Green I Master Mix (Roche), reverse primer (0.5 µM), forward primer (0.5 µM) and 1 µL of DNA sample at 8 ng/µl. Calculation of the relative abundance of the investigated genes attM, qsdA, and qsdB in batches was normalized according level of attM gene in untreated batch and took into account the variation of DNA concentration in batches.

Construction and Screening of the Metagenomic Library

The metagenomic library of the microbial community collected at 42 days from the GCL-treated batch was obtained as follows. DNA was partially digested with Sau3A. The resulting fragments (40–50 kbp) were cloned into the fosmid EpiFos harboring a chloramphenicol-resistance gene, and transferred into E.coli DH10B (LibraGen, Toulouse, France). A total of 29,760 clones constituted the metagenomic library which was distributed into 310 microwell plates. Each clone was individually cultivated in 96-well plates and tested for inactivation of the C6HSL signal (12.5 µM), according the previously described procedures [9].

Sequencing and Phylogenetic Analyses of the NAHL-degrading Fosmid p90H6

Complete sequence of the p90H6 fosmid (46 kbp) harboring a 39 kbp insert was achieved by Sanger sequencing. To this end, a library of 2.0–2.5 kbp fragments was constructed, cloned and sequenced by GATC Biotech (Mulhouse, France). Assembling (7x in average coverage) was performed using DNASTAR Lasergene, and appropriate primers were designed for gap-closure. DNA map was generated by Geneious and orf determination by Bioedit softwares. For each of the putative proteins encoded by the p90H6 insert, 50 homologous sequences were retrieved from GenBank (http://www.ncbi.nlm.nih.gov) after identification by BLASTp. Multiple sequence alignments were carried out with ClustalW. The evolutionary history was inferred using the Maximum Likelihood method based on the JTT matrix-based model [40]. The percentage of replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) was calculated [41]. Initial tree(s) for the heuristic search were obtained automatically as follows. When the number of common sites was <100 or less than one fourth of the total number of sites, the maximum parsimony method was used; otherwise BIONJ method with MCL distance matrix was used. Evolutionary analyses were conducted inMEGA5 [42].

NAHL-degradation Assays by Entire Cells and Cell-free Extract

The pMTHindIII and pMTXhoI plasmids were pME600 derivatives harbouring a 4.5 kbp and 2.3 kbp qsdB-insert, respectively. Overnight cultures of E.coli strain DH5α harboring the p90H6, pME6000 empty vector or the two plasmids pMTHindIII and pMTXhoI were adjusted to OD600 nm 2.0. The resulting suspensions were incubated at 30°C with C6HSL at 25 µM, OC8HSL and C8HSL at 400 nM in LBm medium. Assays were performed up to 48 h, and residual concentration of NAHL was determined using the NAHL biosensors A. tumefaciens biosensor NT1(pZNLR4) and C. violaceum CVO26.

NAHL-degradation assay was also performed using cell-free extract of E. coli harboring the pME6000 and pMTXhoI plasmids. Detection of residual C6HSL was done by HPLC-MS. Overnight cultures of the E.coli strains were centrifuged, washed twice in KPBS buffer (10 mM Na2HPO4 1.76 mM KH2PO4 pH 6.8), and sonicated in the presence of a protease inhibitor cocktail (Sigma P8465). Cell debris were removed by centrifugation (20 min, 4°C, 14,000 g), and the resulting crude cell extract was filtered through polyethersulfone column (10 kDa Vivaspin500). In NAHL-degradation assay, the reaction mixture contained cell free extract at a final concentration of 10 mg/mL and C6HSL at 50 µM in a final volume of 1 ml of KPBS buffer. Assays were incubated at 30°C up to 24 h, and filtered through polyethersulfone column (10 kDa Vivaspin500). Thirty to 50 µL were analyzed by HPLC (Waters Alliance 2690) conjugated with LC-MS/MS (Waters ZQ Mass Spectrometer with single quadrupole system and electrospray ionization). To detect C6HSL, a Gemini C18 5 µm (2.0×150 mm) column was eluted with an isocratic mixture (80/20) of H2O/HCOOH (0.1%) and CH3CN/HCOOH (0.1%). Quantification was performed according a calibration curve generated with pure C6HSL.

Soft Rot Tuber Assays

Cultures of Pectobacterium atrosepticum CFBP6276 and its derivatives harboring the pME6000 and pMTXhoI plasmids were washed twice in NaCl 0.8%. Fourteen tubers of S. tuberosum var. Allians were inoculated with 107 CFU of the plant pathogen, and incubated at 25°C. Five days after infection, the tubers were cut in the middle, and observed. The aggressiveness of the strains was categorized into four classes of diameters of the maceration zone: less than 1 mm, between 1 and 2 mm, between 2 and 5 mm, more than 5 mm. The Kruskal–Wallis test (α = 0.05) allowed the statistical analysis of the maceration categories.

Construction and Purification of the His-tag Derivatives of QsdB, LibA and NylA

The qsdB, libA (G3K3F8 UniProt KB) and nylA (P13398 UniProt KB) genes were synthezised by Eurofins (Germany). Codon usage was optimized for heterologous expression in E. coli. Appropriate restriction sites and six terminal His-codons were incorporated for cloning the genes in the pET11d expression vector. Mutations in the catalytic triad K-S-S of QsdB were constructed using the pET11d-qsdB vector as a template and the QuikChange® II XL site-directed mutagenesis kit (Stratagene). The synthetic forward primers 5′-TTCCTGATCGCAGACCTGGTT (the altered nucleotides underlined), 5′-TCGGGTGGTGCCTCTGGAGGG, and 5′-GGCGGTGGTGCTATCC GCATT, and their appropriate reverse primers were designed for the mutations K70A, S147A, and S171A, respectively. All constructed alleles were verified by DNA-sequencing (GATC Biotech, Mulhouse, France).

QsdB and LibA expression plasmids were transformed into E.coli Rosetta pLysS (Novagen), and NylA expression plasmid in E.coli C41 (Lucigen). The cells were grown in 2YT media at 37°C to OD600 nm 0.6; the temperature was reduced to 28°C, and expression was induced with 0,5 mM IPTG during 5 h. The cells were harvested and stored at −20°C. The pellet was suspended in buffer A (50 mMTris–Cl buffer, pH 8, 20 mM imidazole, 10% glycerol and 500 mM NaCl) and disrupted by sonication. The lysate was cleared by centrifugation at 20,000 g for 30 min at 4°C. The supernatant was filtered through a 0.22 µm membrane (Millipore Stericup) and loaded on a 5 ml Ni-NTA agarose column (GE Healthcare). After washing with 10% of buffer B (50 mMTris–Cl pH 8, 300 mM imidazole, 10% glycerol and 500 mM NaCl), elution was performed with 100% of buffer B. Fractions containing the proteins were pooled, and a buffer exchange was performed on a Superdex 200 10/300 (GE Healthcare) or by dialysis (SlideA Lyzer, Pierce) during 4 hours in Tris 50 mMTris–Cl pH 8, 150 mM NaCl and 10% glycerol. Fractions were analyzed by SDS-PAGE, and those containing the proteins were pooled, concentrated using Vivaspin-10 kDa (Sartorius) and stored at −20°C.

HPLC-MS Identification of Homoserine Lactone Released by the Amidase QsdB

Purified QsdB were incubated at 0.1 mg/mL with 12.5 µM of C6HSL for 24 hours, then 200 µL of the assay was filtered through polyethersulfone column (10 kDa Vivaspin500). The eluate was diluted 1/5 in acetonitrile (CH3CN), and finally, 5 µL were in the HLC column for detection of homoserine lactone, as an amidase product of C6HSL. To detect HSL molecule, a SeQuant Zic-pHilic 5 µM (150 mm×4.6 mm) column was eluted by isocratic mixture (80/20) of H2O/HCOOH (0.1%) and CH3CN/HCOOH (0.1%). Five µL of standard HSL in buffer elution 1% and acetonitrile/water at a ratio 80/20 at different concentrations were also injected through the column to generate a calibration curve.

Results

High Proportion of NAHL-degrading Bacteria in GCL-treated Environments

GCL (0.4 g/L) was introduced twice (1-day and 28-day) in the rhizosphere of Solanum tuberosum var. Allians grown under hydroponic conditions (Fig. 1A). At 42-day, cultured bacteria (CFU/g of fresh root) were significantly more abundant in the GCL-treated batch than in the untreated one (Fig. 1B). Among cultured bacteria, the proportion of C6HSL-degrading bacteria reached 3% in control batch, and 36% in the treated batch (Fig. 1C). In contrast, no significant modification of the percentage of NAHL-producing bacteria was observed in the different batches (Fig. 1D). Partial rrs-sequencing showed that all isolated NAHL-degrading bacteria belong to the Rhodococcus genus. Using high performance liquid chromatography-mass spectrometry (HPLC-MS), GCL concentrations were measured in nutritive solutions and plant tissues of untreated and GCL-treated cultures at 42-day. While GCL was introduced at 0.8 g/L (2×0.4 g/L), it was detected at 3.2 10−4 g/L in the nutritive solution and at 4.9 10−4 mg/g of fresh weight in plant tissues recovered from treated batch at 42-day. GCL was undetectable in the untreated samples. However, over the same period, GCL concentration remained stable in a sterile nutritive solution.

[Figure omitted. See PDF.]

Figure 1. NAHL-degrading community in GCL-treated plant cultures.

(A) Potato plants were cultivated under untreated and GCL-treated conditions during 42 days; vertical arrows indicate the two applications of GCL at 0.4g/L. At 42-day, total cultured bacteria (B), and percentage of NAHL-degraders (C) and NAHL-producers (D) were enumerated and calculated under both the GCL-treated and untreated conditions. The mean of three replicates are shown. Statistically different values (Mann & Whitney, α = 0.05) are noted by different letters.

https://doi.org/10.1371/journal.pone.0065473.g001

Drastic Reshaping of Bacterial Diversity Upon GCL-treatment

GCL-mediated modifications of bacterial diversity in plant environment were explored by rrs-DGGE and rrs-pyrosequencing. After 42 days, total DNA was extracted from the bacterial cells collected from the nutritive solution and used as a template for rrs-amplifications. DGGE profiles revealed a drastic remodeling of bacterial populations after GCL treatments (Fig. 2A). Under GCL-treated conditions, major bands identified by sequencing were generated from members of the Azospirillum genus, while in untreated control they originated from members of the genera Hydrogenophaga and Acetivibrio.

[Figure omitted. See PDF.]

Figure 2. GCL-induced changes in the bacterial rrs-diversity.

At 42-day, rrs-diversity was evaluated by DGGE (A) and pyrosequencing (B, C). In A, for each of the numbered bands, one of the closest rrs sequences was indicated and characterized by GenBank number, the name and taxonomical position of the bacteria of origin (Firm, Firmicutes; α, α-Proteobacteria; β, β-Proteobacteria; γ, γ-Proteobacteria), and the similarity index (SI) calculated at Ribosomal Database Project (http://rdp.cme.msu.edu/). The rrs-pyrosequencing data were analyzed at the class level (B) and, within the class of alphaproteobacteria, at the genus level (C).

https://doi.org/10.1371/journal.pone.0065473.g002

The rrs-pyrosequencing (NCBI accession number PRJNA196890) gave a deeper view of bacterial community and confirmed the high abundance of Azospirillum-related bacteria in the GCL-treated samples. All 86,960 usable sequences (36,847 and 50,113 associated to the untreated and GCL-treated samples, respectively) were analyzed by hierarchical clustering to define clusters at 97% and 95% sequence identity. These served as Operational Taxonomy Units (OTUs) for calculation of Chao1 and Shannon indexes. In the untreated samples, cluster number reached 834 (610 with 95% cut-off), whereas it was only 271 (188 with 95% cut-off) in GCL-treated samples. Using the 97% cut-off, the Chao1 and Shannon values were higher in the untreated sample (1204 and 4.03, respectively) than in treated samples (397 and 1.94). Similar data were obtained with the 95% cut-off: the Chao1 and Shannon values were 824 and 4.00 under untreated condition and 309 and 1.29 under GCL-treatment. These data indicate that GCL-treatment strongly stimulated the growth of a limited number of OTUs, and hence reduced the diversity of the resident community. Diversity analysis of the GCL-induced bias was refined by identifying and comparing rrs-sequences at the class and genus levels using tools available from the Ribosomal Database Project Website. While α-proteobacterial rrs-sequences dominated in GCL-treated samples (more than 95%), those identified in the untreated samples belonged mostly to β-proteobacteria (26%), unclassified bacteria (21%) and γ-proteobacteria (18%) (Fig. 2B). As mentioned above, Azospirillum appeared as the dominant genus that represented up to 92% of the identified α-proteobacteria in treated samples (Fig. 2C). In contrast, the α-proteobacterial diversity of untreated batch samples was larger, with more than 10 genera of which emerged Sphingobium (30%) and Rhizobium (15%).

Identification of a Metagenomic Fosmid Encoding NAHLase Activity

Bacterial DNA was extracted from the GCL-treated batch at 42-day, and used to construct a metagenomic library of 29 760 fosmids, each of which exhibited a DNA insert of 40 to 50 kbp. Hence, the metagenomic library represented up 1.5 Gbp of environmental DNA. All metagenomic clones were individually tested for their capacity to disrupt QS using hexanoyl-homoserine lactone (C6HSL) as a prototypic signal and Chromobacterium violaceum CVO26 as a bacterial biosensor. A single fosmid, p90H6, conferring NAHL-degradation ability was identified and sequenced.

The fosmid p90H6 contains an insert of 39,632 bp encoding 34 complete ORFs (NCBI accession number JQ292794). Protein sequence comparison of all the deduced ORFs was performed using BLAST searches upon available sequences in GenBank database. This analysis revealed that out of 34 ORFs, 31 showed highest similarity with known or predicted proteins from bacteria belonging to the phylum of Proteobacteria. Furthermore, a phylogenetic analysis of 50 homologous proteins was carried out for each of the 34 ORFs allowing the identification of the closest orthologs which are showed in the Table 1. Phylogenetical approach confirmed the relationship between DNA insert of p90H6 and the class of Proteobacteria, and especially γ-proteobacteria. Two remarkable gene clusters were identified. The first one from orf12 to orf18 could be involved in metabolism of acetophenone-related compounds; the second one from orf25 to orf34 in the catabolism of the plant hormone indole-3-acetic acid. Another remarkable feature was the similarity of orf1, which is located immediately downstream the lacZ-promoter of the fosmid, with amidase or putative amidase genes. This orf, studied in this work, was named qsdB for quorum-sensing degradation.

[Figure omitted. See PDF.]

Table 1. ORFs of the p90H6 DNA-insert.

https://doi.org/10.1371/journal.pone.0065473.t001

Relative Abundance of qsdB in GCL-treated Bacterial Community

In parallel to the enzymatic characterization of QsdB, the relative abundance of genes encoding NAHL-degrading enzymes was compared by qPCR in the DNA extracted from the GCL-treated and untreated rhizosphere. In addition to qsdB, we quantified qsdA from Rhodococcus erythropolis [12] the proportion of which increased after GCL-treatment [36], and attM which is present in Agrobacterium [6] and related α-proteobacteria. The lowest level of attM in the untreated condition was arbitrary chosen as a reference level. The relative abundance of qsdA and more notably that of qsdB and attM were higher in GCL-treated batch than in untreated batch (Fig. 3). However, in the GCL-treated condition, the relative abundance of qsdB was higher than that of attM and qsdA.

[Figure omitted. See PDF.]

Figure 3. Relative abundance of NAHLase-encoding genes in GCL-treated and untreated plant cultures.

Relative abundance of the qsdB (A), qsdA (B) and attM (C) genes in GCL-treated and untreated batches at 42-day was measured by qPCR. The attM intensity under the untreated condition (C) was used as a normalized reference (arbitrary value = 1) for calculation of the relative abundance of all genes.

https://doi.org/10.1371/journal.pone.0065473.g003

Characterization of qsdB and Encoded NAHL-amidase Activity

The GC content of the metagenomic DNA-sequence reached a mean value of 67.2%, and appeared to be roughly stable throughout the whole insert, though some GC% differences could be seen in between some ORFs (Fig. 4A). Two pME6000 plasmid derivatives, pMTHindIII and pMTXhoI harboring qsdB on 4.5 kbp and 2.3 kbp inserts, respectively, were constructed (Fig. 4A). In addition to qsdB, the pMTXhoI smallest insert contains an additional orf encoding a putative inner-membrane translocator (117 amino acids) with no expected hydrolytic properties. E. coli cells harbouring p90H6, pMTHindIII or pMTXhoI were able to inactivate over 95% of the introduced NAHLs, which were C6HSL, octanoylhomoserine lactone (C8HSL), and 3-oxo-octanoylhomoserine lactone (OC8HSL), in less than 24 hours (Fig. 4B). A HPLC-MS analysis confirmed the disappearance of the QS-signal C6HSL in the presence of cell-free extract from the E. coli strain harboring the pMTXhoI plasmid (Fig. 4C).

[Figure omitted. See PDF.]

Figure 4. Physical map of the fosmid p90H6 and sub-cloning of qsdB coding for NAHLase activity.

In A, GC content (%) and orientation (+/−) of the 34 putative orfs of the fosmid p90H6, and physical map of the pME6000-derivatives pMTHindIII and pMTXhoI barboring qsdB (orf1). In B, residual level of C6HSL, OC8HSL or C8HSL measured in the presence of E. coli strain DH5α harboring the empty vector pME6000, the fosmid p90H6, and the constructed pMTHindIII and pMTXhoI (symbols of the three formers are superimposed in the graphs).Three replicates were done. In C, C6HSL analysis by HPLC/MS: mass, retention time, and quantification of C6HSL before (t = 0) and 24 hours after incubation in the presence of cell-free extracts of E. coli strains DH5α(pME6000) and DH5α (pMTXhoI).

https://doi.org/10.1371/journal.pone.0065473.g004

To evaluate the quorum-quenching capacity of the pMTXhoI plasmid, this later was introduced into the plant pathogen P. atrosepticum CFBP6276, in which expression of virulence factors is controlled by QS with OC8HSL as a major signal [43]. The constructed P. atrosepticum (pMTXhoI) exhibited a ten fold decrease of the NAHL level (Fig. 5A) and a reduced level of symptoms on potato tubers (Fig. 5B) as compared to the control strain P. atrosepticum (pME6000). Finally, the QsdB enzyme was purified using a His-tag procedure (Fig. 6A), verified for its capacity to inactivate QS-signals (Fig. 6B), and tested for amidase activity using HPLC-MS identification and quantification of homoserine lactone as a reaction product (Fig. 6C). After 24-h incubation in the presence of QsdB, 100% of the introduced C6HSL (13 µM) were converted into the amidase product homoserine lactone, the concentration of which reached 14±1 µM.

[Figure omitted. See PDF.]

Figure 5. QsdB-mediated quorum-quenching in P. carotovorum.

In A, relative abundance of NAHLs (OC8HSL as the main QS-signal) in wild-type P. carotovorum CFBP6276 (Pca6276), and P.carotovorum harboring the empty vector pME6000 and the pMTXhoI plasmid expressing QsdB. In B, symptoms on potato tubers (n = 14) inoculated with Pca6276, Pca6276 (pME6000), Pca6276 (pMTXhoI) or NaCl 0.8% as a negative control were classified according to their importance (see material and methods) and compared using Kruskal-Wallis test (P<0.05). Different letters indicate statistical differences between the compared plant pathogens.

https://doi.org/10.1371/journal.pone.0065473.g005

[Figure omitted. See PDF.]

Figure 6. Purification and amidase activity of QsdB.

In A, SDS-PAGE analysis of the S-200 fraction containing the purified QsdB. In B, degradation kinetics of OC8HSL and C6HSL in the presence of QsdB at 0.1 mg/ml; the protein-free buffer was used as a control. Three replicates were done. In C, Amidase activity of QsdB toward C6HSL releases hexanoic acid and homoserine lactone which was characterized and quantified by HLPC/MS as described in the material and method section.

https://doi.org/10.1371/journal.pone.0065473.g006

QsdB Belongs to the Amidase Signature (AS) Family and Requires the Catalytic Triad K-S-S for NAHLase Activity

BLAST analysis highlighted the AS family enzyme OctHD (47% identity with QsdB) among the 50 closest relatives of QsdB. A refined phylogenetic analysis of QsdB confirmed its position within the AS family (Fig. 7). QsdB belongs to a sub-cluster which encompasses the NylA proteins from Ruegeria, Arthrobacter and Pseudomonas, the ω-laurolactam and ω-octalactam hydrolases from Rhodococcus, the linuron-degrading enzyme LibA of Variovorax, and two predicted amidases in Burkholderia cenocepacia J2315 and Burkholderia sp. 383. To evaluate whether the NAHLase activity is conserved within the AS family, the proteins LibA and NylA were purified using the same His-tag approach as described for QsdB. The purified LibA and NylA did not inactivate NAHLs suggesting a probable specialization of these enzymes towards their xenobiotic substrates.

[Figure omitted. See PDF.]

Figure 7. QsdB within the amidase signature (AS) family.

Phylogenetic relationship between twenty AS-family members including the metagenomic NAHL-amidase QsdB was inferred using the Maximum Likelihood method based on the JTT matrix-based model. Percentage of the replicate trees in which the associated taxa clustered together in the bootstrap test (500 replicates) is shown next to the branches. The tree is drawn to scale, with branch lengths measured in the number of substitutions per site.

https://doi.org/10.1371/journal.pone.0065473.g007

Three residues, collectively named the catalytic triad K-S-S, have been described as being involved in the amidase activity of several AS family members [44]–[46]. A sequence alignment revealed that these amino-acid residues were also conserved in QsdB at positions Lys70, Ser147, and Ser171 (Fig. 8A). To evaluate their implication in NAHLase activity, three mutants of QsdB affected in each of the K-S-S residues were constructed, purified and tested against NAHLs (Fig. 8B). All constructed protein mutants lost their NAHLase activity as compared to the wild type QsdB, a result that demonstrates the involvement of the AS-family triad K-S-S in the degradation of the QS signals.

[Figure omitted. See PDF.]

Figure 8. The AS-family catalytic triad K-S-S is required for NAHLase activity of QsdB.

In A, the amino acids K70, S147 and S171 of the AS-family catalytic triad are underlined in the QsdB sequence. In B, concentration of residual OC8HSL after 24 hr-incubations in the presence of the wild type protein QsdB and its constructed derivatives K70A, S147A and S171A, all at 0.1 mg/ml; the protein-free reaction buffer was used as a control.

https://doi.org/10.1371/journal.pone.0065473.g008

Discussion

The novel metagenomic NAHL-amidase QsdB identified in this study is distantly related to the previously identified NAHLases and belongs to the AS family that comprises most than 200 proteins from diverse origins including Bacteria, Archaea and Eukarya [44]. The amidase activity of AS enzymes requires the K-S-S triad [44]–[45], which is also involved in the NAHLase activity of QsdB. Despite their phylogenetic proximity and common ability to hydrolyse amide bonds, individual representatives of the AS family show very distinctive substrate-specificities [45]. It has been suggested that this is due to the involvement of residues outside the catalytic triad K-S-S in substrate binding [44]. Remarkably, QsdB belongs to an AS sub-cluster, which encompasses the bacterial amidases NylA, OctHD and LibA cleaving xenobiotic compounds such as nylon and the herbicide linuron [47]–[49]. To our knowledge, QsdB is the first enzyme of the NylA-OctHD-QsdB-LibA sub-cluster to be involved in the degradation of natural compounds, i.e. NAHLs. One can hypothesized that QsdB-related NAHLases would be ancestors from which have evolved the enzymes cleaving structurally-related xenobiotics, such as linuron and nylon. The pre-existing skeleton of AS-family NAHLase would therefore explain, in part, why nylon- and linuron-degrading enzymes have rapidly emerged in microbial communities even though such compounds have been only recently introduced in the soil as a consequence of industrial and agricultural activities [50]–[51]. Such evolved genes encoding for substrate-specialized enzymes could be thereafter propagated by horizontal transfer [50].

QsdB-encoding gene belongs to a large environmental-DNA insert (39 kpb) of which most of the ORFs (31/34) were phylogenetically associated with Proteobacteria. In addition to qsdB, this environmental DNA harbored two remarkable regions, from orf12 to orf18 and from orf25 to orf34, which could be involved in the metabolism of acetophenone-related compounds and indole-3-acetic acid. Nine of the 10 genes which are known to be involved in the catabolism of indole-3-acetic acid [52] are present in the metagenomic fosmid. Their syntheny designated γ-proteobacteria as the closer bacteria. The presence of genes encoding for auxin catabolism strongly suggests that the trapped environmental DNA could belong to a plant-associated bacterium. This feature is in line with the rhizospheric origin of the analyzed community. Noticeably, this metagenomic DNA did not show robust similarities with the known genomes of Azospirillum (α-proteobacteria) which was the most abundant genus in the analyzed environment, neither with Rhodococcus (Actinobacteria) which was the dominant genus among the isolated NAHL-degrading bacteria. Previous studies showed that, among bacterial isolates recovered from GCL-treated rhizosphere, the NAHL-degrading Rhodoccocus grow on GCL as a sole carbon source [53], while Azospirillum assimilate GCL, but can not inactivate NAHLs [26]–[27]. The qPCR analyses in metagenomic DNA confirmed that qsdB was most abundant than the genes qsdA which belongs to Rhodococcus and attM which is associated to Agrobacterium and related α-proteobacteria. However, the frequency of trapped NAHLase-expressing fosmid among the constructed library is quite similar to that observed in previous studies [9], [33]. The apparent low frequency of positive NAHL-degrading clones may be explained by the combined constraints of large size of the inserts (50 kbp) and necessity to drive the gene expression from the promoter located at the cloning site of the vector in an E. coli host (see reviews about functional metagenomics in [32], [54]).

Finally, this study increases our knowledge on the dynamics of bacterial community upon hydroponic culture of potato plants in green-house. While the previous 454-pyrosequencing analyses of rrs-amplicons focused on the early stages (0 to 4 weeks) of plant growth culture [36], this work focused on a later stage at 7 weeks. In the two analyzed conditions (with or without GCL-treatment), the α-proteobacterium Azospirillum dominated the bacterial 7-week community, while Acidovorax and Flavobacterium were the most abundant genera in the earliest communities [36]. This work highlights the remarkable capacity of the Azospirillum genus to take advantage of the simultaneous constraints that are plant and aquatic environments, as suggested by comparative genomics [55]. This analysis should be repeated in an independent experiment for robustly deciphering the community dynamics.

In a conclusion, this work highlights the power of metagenomics to access unknown functional diversity, and reveals a common phylogenetical origin between AS-family amidases involved in the degradation of cell-to-cell signals and xenobiotics.

Acknowledgments

The authors thank A. Raffoux (ISV, CNRS) for a helpful contribution in screening of the metagenomic library, and O. Thoison (ICSN, CNRS) for analysis of the NAHL by-products.

Author Contributions

Conceived and designed the experiments: MT AB-C DF. Performed the experiments: MT AB-C AV FP. Analyzed the data: MT AB-C SM YD DF. Wrote the paper: MT AB-C DF.

Citation: Tannières M, Beury-Cirou A, Vigouroux A, Mondy S, Pellissier F, Dessaux Y, et al. (2013) A Metagenomic Study Highlights Phylogenetic Proximity of Quorum-Quenching and Xenobiotic-Degrading Amidases of the AS-Family. PLoS ONE 8(6): e65473. https://doi.org/10.1371/journal.pone.0065473

References

1. Fuqua WC, Winans SC, Greenberg EP (1994) Quorum sensing in bacteria: the LuxR/LuxI family of cell density-responsive transcriptional regulators. J Bacteriol: 269–275.

2. Boyer M, Wisniewski-Dyé F (2009) Cell-cell signalling in bacteria: not simply a matter of quorum. FEMS Microbiol Ecol 70: 1–19.

3. Uroz S, Dessaux Y, Oger P (2009) Quorum sensing and quorum quenching: the yin and yang of bacterial communication. Chembiochem 26: 205–216.

4. Zhang LH (2003) Quorum quenching and proactive host defense. Trends Plant Sci 8: 238–244.

5. Dong YH, Xu JL, Li XZ, Zhang LH (2000) AiiA, an enzyme that inactivates the acylhomoserine lactone quorum-sensing signal and attenuates the virulence of Erwinia carotovora. Proc Natl Acad Sci USA 97: 3526–3531.

6. Zhang HB, Wang LH, Zhang LH (2002) Genetic control of quorum sensing signal turnover in Agrobacterium tumefaciens. Proc Natl Acad Sci U S A 99: 4638–4643.

7. Carlier A, Uroz S, Smadja B, Fray R, Latour X, et al. (2003) The Ti plasmid of Agrobacterium tumefaciens harbors an attM paralogous gene, aiiB, also encoding N-Acyl homoserine lactonase activity. Appl Environ Microbiol 69: 4989–4993.

8. Park SY, Lee SJ, Oh TK, Oh JW, Koo BT, et al. (2003) AhID, an N-acylhomoseine lactonase in Arthrobacter sp., and predicted homologues in other bacteria. Micobiology 149: 1541–1550.

9. Riaz K, Elmerich C, Moreira D, Raffoux A, Dessaux Y, et al. (2008) A metagenomic analysis of soil bacteria extends the diversity of quorum-quenching lactonases. Environ Microbiol 10: 560–570.

10. Mei GY, Yan XX, Turak A, Luo ZQ, Zhang LQ (2010) AidH, an alpha/beta-hydrolase fold family member from an Ochrobactrum sp. strain, is a novel N-acylhomoserine lactonase. Appl Environ Microbiol 76: 4933–4942.

11. Wang WZ, Morohoshi T, Ikenoya M, Someya N, Ikeda T (2010) AiiM, a novel class of N-acylhomoserine lactonase from the leaf-associated bacterium Microbacterium testaceum. Appl Environ Microbiol. 76: 2524–2530.

12. Uroz S, Oger PM, Chapelle E, Adeline MT, Faure D, et al. (2008) A Rhodococcus qsdA-encoded enzymes defines a novel class of large-spectrum quorum-quenching lactonases. Appl Environ Microbiol 74: 1357–1366.

13. Lin YH, Xu JL, Hu J, Wang LH, Ong SL, et al. (2003) Acyl-homoserine lactone acylase from Ralstonia strain XJ12B represents a novel and potent class of quorum-quenching enzymes. Mol Microbiol 47: 849–860.

14. Huang JJ, Petersen A, Whiteley M, Leadbetter JR (2006) Identification of QuiP, the product of gene PA1032, as the second acyl-homoserine lactone acylase of Pseudomonas aeruginosa PAO1. Appl Environ Microbiol 72: 1190–1197.

15. Park SY, Kang HO, Jang HS, Lee JK, Koo BT, et al. (2005) Identification of extracellular N-acylhomoserine lactone acylase from a Streptomyces sp and its application to quorum quenching. Appl Environ Microbiol 71: 2632–2641.

16. Huang JJ, Han J, Zhang L, Leadbetter JR (2003) Utilization of acyl-homoserine lactone quorum signals for growth by a soil pseudomonad and Pseudomonas aeruginosa PAO1. Appl. Environ Microbiol 69: 5941–5949.

17. Romero M, Diggle SP, Heeb S, Camara M, Otero A (2008) Quorum quenching activity in Anabaena sp. PCC 7120: identification of AiiC, a novel AHL-acylase. FEMS Microbiol Lett 280: 73–80.

18. Bijtenhoorn P, Mayerhofer H, Müller-Dieckmann J, Utpatel C, Schipper C, et al. (2011a) A Novel Metagenomic Short-Chain Dehydrogenase/Reductase Attenuates Pseudomonas aeruginosa Biofilm Formation and Virulence on Caenorhabditis elegans. PLoS One 6: 262–278.

19. Bijtenhoorn P, Schipper C, Hornung C, Quitschau M, Grond S, et al. (2011b) BpiB05, a novel metagenome-derived hydrolase acting on N-acylhomoserine lactones. J Biotechnol 155: 86–94.

20. Kaufmann GF, Sartorio R, Lee SH, Rogers CJ, Meijler MM, et al. (2005) Revisiting quorum sensing: Discovery of additional chemical and biological functions for 3-oxo-N-acylhomoserine lactones. Proc Natl Acad Sci USA 102: 309–314.

21. Haudecoeur E, Tannières M, Cirou A, Raffoux A, Dessaux Y, et al. (2009a) Different regulation and roles of lactonases AiiB and AttM in Agrobacterium tumefaciens C58. Mol. Plant-Microbe Interact 22: 529–537.

22. Haudecoeur E, Planamente S, Cirou A, Tannières M, Shelp BJ, et al. (2009b) Proline antagonizes GABA-induced quenching of quorum-sensing in Agrobacterium tumefaciens. Proc Natl Acad Sci USA 106: 14587–14592.

23. Dong YH, Wang LH, Xu JL, Zhang HB, Zhang XF, et al. (2001) Quenching quorum-sensing-dependent bacterial infection by an N-acyl homoserine lactonase. Nature 411: 813–817.

24. Uroz S, D'Angelo-Picard C, Carlier A, Elasri M, Sicot C, et al. (2003) Novel bacteria degrading N-acylhomoserine lactones and their use as quenchers of quorum-sensing-regulated functions of plant-pathogenic bacteria. Microbiology 149: 1981–1989.

25. Dong YH, Zhang XF, Xu JL, Zhang LH (2004) Insecticidal Bacillus thuringiensis silences Erwinia carotovora virulence by a new form of microbial antagonism, signal interference. Appl Environ Microbiol 70: 954–60.

26. Cirou A, Diallo S, Kurt C, Latour X, Faure D (2007) Growth promotion of quorum-quenching bacteria in the rhizosphere of Solanum tuberosum. Environ Microbiol 9: 1511–1522.

27. Cirou A, Raffoux A, Diallo S, Latour X, Dessaux Y, et al. (2011) Gamma-caprolactone stimulates growth of quorum-quenching Rhodococcus populations in a large-scale hydroponic system for culturing Solanum tuberosum. Res Microbiol 162: 945–950.

28. Kennedy J, O’Leary ND, Kiran GS, Morrissey JP, O’Gara F, et al. (2011) Functional metagenomic strategies for the discovery of novel enzymes and biosurfactants with biotechnological applications from marine ecosystems. J Appl Microbiol 111: 787–799.

29. Monier JM, Demanèche S, Delmont TO, Mathieu A, Vogel TM, et al. (2011) Metagenomic exploration of antibiotic resistance in soil. Curr Opin Microbiol. 14: 229–235.

30. Piel J (2011) Approaches to capturing and designing biologically active small molecules produced by uncultured microbes. Annu Rev Microbiol 65: 431–453.

31. Gabor EM, de Vries EJ, Janssen DB (2004) Construction, characterization, and use of small-insert gene banks of DNA isolated from soil and enrichment cultures for the recovery of novel amidases. Environ Microbiol 6: 948–958.

32. Taupp M, Mewis K, Hallam SJ (2011) The art and design of functional metagenomic screens. Curr Opin Biotechnol. 22: 465–472.

33. Schipper C, Hornung C, Bijtenhoorn P, Quitschau M, Grond S, et al. (2009) Metagenome-derived clones encoding two novel lactonase family proteins involved in biofilm inhibition in Pseudomonas aeruginosa. Appl Environ Microbiol 75: 224–233.

34. Hao Y, Winans SC, Glick BR, Charles TC (2010) Identification and characterization of new LuxR/LuxI-type quorum sensing systems from metagenomic libraries. Environ Microbiol 12: 105–117.

35. Guan C, Ju J, Borlee BR, Williamson LL, Shen B, et al. (2007) Signal mimics derived from a metagenomic analysis of the gypsy moth gut microbiota. Appl Environ Microbiol 73: 3669–3676.

36. Cirou A, Mondy S, An S, Charrier A, Sarrazin A, et al. (2012) Efficient biostimulation of native and introduced quorum-quenching Rhodococcus erythropolis populations is revealed by a combination of analytical chemistry, microbiology, and pyrosequencing. Appl Environ Microbiol 78: 481–92.

37. Cha C, Gao P, Chen YC, Shaw PD, Farrand SK (1998) Production of acyl-homoserine lactone quorum sensing signals by gram-negative plant-associated bacteria. Mol Plant Microbe Interact 11: 1119–1129.

38. McClean KH, Winson MK, Fish L, Taylor A, Chhabra SR, et al. (1997) Quorum sensing and Chromobacterium violaceum: exploitation of violacein production and inhibition for the detection of N-acylhomoserine lactones. Microbiology 143: 3703–3711.

39. D’Angelo-Picard C, Faure D, Carlier A, Uroz S, Raffoux A, et al. (2004) Dynamics of bacterial populations in the rhizosphere of tobacco plants producing – or not – the quorum sensing signals hexanoyl- and 3-oxohexanoyl-homoserine lactone. FEMS Microbiol Ecol 51: 19–29.

40. Jones DT, Taylor WR, Thornton JM (1992) The rapid generation of mutation data matrices from protein sequences. Computer Applications in the Biosciences 8: 275–282.

41. Felsenstein J (1985) Confidence limits on phylogenies: an approach using the bootstrap. Evolution 39: 783–791.

42. Tamura K, Peterson D, Peterson N, Stecher G, Nei M, et al. (2011) MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol 28 2731–2739.

43. Smadja B, Latour X, Faure D, Chevalier S, Dessaux Y, et al. (2004) Involvement of N-acylhomoserine lactones throughout plant infection by Erwinia carotovora subsp. atroseptica (Pectobacterium atrosepticum). Mol Plant Microbe Interact 17: 1269–1278.

44. Labahn J, Neumann S, Büldt G, Kula MR, Granzin J (2002) An alternative mechanism for amidase signature enzymes. J Mol Biol 322: 1053–1064.

45. Valiña AL, Mazumder-Shivakumar D, Bruice TC (2004) Probing the Ser-Ser-Lys catalytic triad mechanism of peptide amidase: computational studies of the ground state, transition state, and intermediate. Biochemistry 43: 15657–72.

46. Negoro S, Shibata N, Tanaka Y, Yasuhira K, Shibata H, et al. (2012) Three-dimensional structure of nylon hydrolase and mechanism of nylon-6 hydrolysis. J Biol Chem. 287: 5079–5090.

47. Kanagawa K, Negoro S, Takada N, Okada H (1989) Plasmid dependence of Pseudomonas sp. strain NK87 enzymes that degraded 6-aminohexanoate-cyclic dimer. J Bacteriol 171: 3181–3186.

48. Fukuta Y, Koizumi S, Komeda H, Asano Y (2010) A new aryl acylamidase from Rhodococcus sp. strain Oct1 acting on o-lactams: its characterization and gene expression in Escherichia coli. Enzyme and Microb Technol 46: 237–245.

49. Bers K, Leroy B, Breuglemans P, Albers P, Lavigne R, et al. (2011) A novel linuron hydrolase identified by genomic-proteomic analysis of phenylurea herbicide mineralization by Variovorax sp. strain SRS16. Appl Environ Microbiol 77: 8754–8764.

50. Okada H, Negoro S, Kimura H, Nakamura S (1983) Evolutionary adaptation of plasmid-encoded enzymes for degrading nylon oligomers. Nature 306: 203–206.

51. Tsuchiya K, Fukuyama S, Kanzaki N, Kanagawa K, Negoro S, et al. (1989) High homology between 6-aminohexanoate-cyclic-dimer hydrolases of Flavobacterium and Pseudomonas strains. J Bacteriol 171: 3187–3191.

52. Leveau JHJ, Gerards S (2008) Discovery of a bacterial gene cluster for catabolism of the plant hormone indole 3-acetic acid. FEMS Microbiol Ecol 65: 238–250.

53. Barbey C, Crépin A, Cirou A, Budin-Verneuil A, Orange N, et al. (2012) Catabolic pathway of gamma-caprolactone in the biocontrol agent Rhodococcus erythropolis. J Proteome Res 11: 206–216.

54. Faure D, Tannières M, Mondy S, Dessaux Y (2011) Recent contributions of metagenomics to studies on quorum-sensing and plant-pathogen interactions. In Metagenomics: current innovations and future trends. D. Marco (ed) Caister Academic Press, Norfolk (UK), 253–263.

55. Wisniewski-Dyé F, Borziak K, Khalsa-Moyers G, Alexandre G, Sukharnikov LO, et al. (2011) Azospirillum genomes reveal transition of bacteria from aquatic to terrestrial environments. PLoS Genet 7(12): e1002430.