Introduction {#s1}
============

Bacteria of the plant microbiome are exposed to many plant-derived compounds that are likely to play important roles in plant-bacteria communication and facilitate bacterial growth and colonization of plant niches. For example, lipid-derived molecules of both microbes and plants have recently been shown to play a role as signals for plant-bacteria communication \[[@R1]\]. In addition, fatty acids are also the precursors for the synthesis of compounds involved in intracellular plant signalling including jasmonic and azelaic acids, which are part of the cascade that triggers systemic defences \[[@R1]\].

Azelaic acid ([Fig. 1a](#F1){ref-type="fig"}) is a straight chained nine-carbon (C~9~) saturated dicarboxylic acid that occurs naturally in whole-grain cereals, rye and barley \[[@R3]\]. Azelaic acid in plants is derived from lipids via lipid peroxidation and primes systemic defence in *Arabidopsis*, suggesting that this mobile molecule can move systemically, most likely through the plant vasculature, and contribute to long-distance signalling in plant defence against pathogens. Azelaic acid specifically promotes disease resistance by priming salicylic acid (SA) signalling as part of systemic plant immunity \[[@R4]\]. Recently, plant-associated *[Pseudomonas syringae](http://doi.org/10.1601/nm.10982)* pathovars have been reported to most likely produce azelaic acid, hinting a role in plant-bacteria communication \[[@R6]\]. In soil, azelaic acid along with other short-chain fatty acids are considered the products of microbial enzymatic oxidation of unsaturated lipids and aromatic hydrocarbons \[[@R7]\]. Several reports have shown that bacteria such as *[Acinetobacter](http://doi.org/10.1601/nm.2765)* spp., *[P. fluorescens](http://doi.org/10.1601/nm.2606)* and *[P. aeruginosa](http://doi.org/10.1601/nm.2553)* are capable of degrading dicarboxylic acids and utilizing them as their sole carbon source \[[@R8]\]. Established work reported the isolation and identification of a *[Pseudomonas](http://doi.org/10.1601/nm.2552)* strain that can very efficiently degrade azelaic acid and utilize it as the sole source of carbon and energy \[[@R10]\]. The bacterial isolate was classified as *Pseudomonas azelaica* DSM 9128 and later reclassified as *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* \[[@R12]\].

![Structure of dicarboxylic acids and *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 phylogenetic tree. (a) Structure of dicarboxylic acids used in this study. (b) Phylogenetic tree showing the phylogenetic position of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* strain DSM9128 based on evolutionary distance of concatenated COGs listed in Table 3. Included in the tree are *[Pseudomonas](http://doi.org/10.1601/nm.2552)* and *[Marinobacter](http://doi.org/10.1601/nm.2853)* species, which have a homologue of AzeR (see text).](mic-166-73-g001){#F1}

Initial experiments from Janota--Bassalik and Wright using *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 \[[@R10]\] support the hypothesis that the trait of degrading the azelaic acid is a consequence of *β*-oxidation steps analogous to those of fatty acid degradation. Due to its possible role in plant-bacteria communication, it was of interest to isolate a bacterial regulator that can respond to the presence of azelaic acid and determine its presence in plant-associated bacteria. In this study, we report genetic evidence of azelaic acid degradation and the identification of a transcriptional regulator that responds to azelaic acid in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128. Three mutants possessing transposon insertions in genes of acyl-CoA ligase, of acyl-CoA dehydrogenase and in an isocitrate lyase display a growth deficiency in the presence of azelaic acid. Studies on transcriptional regulation of these genes resulted in the identification of a IclR family repressor that we designated as AzeR that specifically responds to azelaic acid. AzeR is widespread in *[Pseudomonas](http://doi.org/10.1601/nm.2552)* and *[Marinobacter](http://doi.org/10.1601/nm.2853)* species.

Methods {#s2}
=======

Bacterial strains, constructs, media and growth conditions {#s2-1}
----------------------------------------------------------

*[P. nitroreducens](http://doi.org/10.1601/nm.2652)* and *[Escherichia coli](http://doi.org/10.1601/nm.3093)* strains DH5α \[[@R13]\], CC118 \[[@R14]\] and HB101 \[[@R15]\] were routinely grown at 28 and 37 °C respectively, in Luria--Bertani broth/agar medium \[[@R16]\] or in defined M9 minimal medium \[[@R15]\]. All plasmid constructs used in this study are listed in [Table 1](#T1){ref-type="table"}. When azelaic acid was used as the carbon source, it was added to the growth medium at a final concentration of 0.1--0.2 % (5--10 mM). The following antibiotic concentrations were used: ampicillin (Amp) 100 µg ml^−1^, kanamycin (Km) 100 µg ml^−1^, gentamicin (Gm) 50 µg ml^−1^, tetracycline (Tc) 40 µg ml^−1^, nalidixic (Nx) 25 µg ml^−1^ and nitrofurantoin (Nf) 50 µg ml^−1^ for *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* and Amp100 μg ml^−1^, Km50 μg ml^−1^, Gm 15 µg ml^−1^ and Tc 15 µg ml^−1^ for *[E. coli](http://doi.org/10.1601/nm.3093)* strains.

###### 

Plasmids used in this study

+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| Plasmids          | Characteristics                                                                                                   | References   |
+===================+===================================================================================================================+==============+
| pGem T-Easy       | Cloning vector; Amp^R^ *LacZ* reporter                                                                            | (Promega)    |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pBluescript II KS | Cloning vector; Amp^R^ *LacZ* reporter, sito polylinker, da pUC19                                                 | (Stratagene) |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pLAFR3            | Broad-host-range cosmid cloning vector,                                                                           | \[[@R37]\]   |
|                   |                                                                                                                   |              |
|                   | IncP1; Tet^r^                                                                                                     |              |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pMP220            | Promoter probe vector, IncQ; Tet^R^                                                                               | \[[@R38]\]   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pKNOCK- Km        | Conjugative suicide vector; Km^R^                                                                                 | \[[@R14]\]   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pCRS530           | pAB2001 derivative containing CAS-GNm with flanking region NotI-HindIII-EcoRI-SmaI; Ap^r^, Km^r^                  | \[[@R19]\]   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pRK2013           | helper triparental conjugation, Km^r^                                                                             | \[[@R18]\]   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pCos1             | pLAFR3 carrying *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genomic DNA                         | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pCos 5            | pLAFR3 carrying *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genomic DNA                         | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pCos G            | pLAFR3 carrying *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genomic DNA                         | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pCos M            | pLAFR3 carrying *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genomic DNA                         | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pKnReg            | Internal PCR fragment of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 *aze*R cloned in Pknock-Km | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pMPReg            | Promoter of transcriptional regulator gene cloned in pMP220 vector                                                | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pMPLig            | Promoter of Acyl-CoA ligase cloned in pMP220 vector                                                               | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pMPDeh            | Promoter of Acyl-CoA dehydrogenase cloned in pMP220 vector                                                        | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+
| pMPIsocit         | Promoter of Isocitrate lyase cloned in pMP220 vector                                                              | This study   |
+-------------------+-------------------------------------------------------------------------------------------------------------------+--------------+

Recombinant DNA techniques {#s2-2}
--------------------------

DNA manipulations, including digestion with restriction enzymes, agarose gel electrophoresis, purification of DNA fragments, ligation with T4 DNA ligase, transformation of *[E. coli](http://doi.org/10.1601/nm.3093)*, colony hybridization and radioactive labelling by random priming, were performed as previously described \[[@R15]\]. Southern hybridizations were performed using Hybond-N+membrane (Amersham Pharmacia Biotech). Plasmids were purified using EuroClone columns (EuroClone S.p.A., Italy). Total DNA from *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* was isolated with the sarkosyl-pronase lysis method Better *et al.* \[[@R17]\]. Triparental matings to mobilize DNA from *[E. coli](http://doi.org/10.1601/nm.3093)* to *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* were carried out with the helper strain *[E. coli](http://doi.org/10.1601/nm.3093)* (pRK2013) \[[@R18]\]. PCR amplifications were performed on *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* genomic DNA using GoTaq Flexi DNA Polymerase (Promega, Madison, WI, USA) and the automated sequencing was performed by GATC-Biotech services (Kostanz, Germany). The oligonucleotide primers used in this study are listed in [Table 2](#T2){ref-type="table"}.

###### 

Primers used in this study

  Primer name          Sequence of primer                     References
  -------------------- -------------------------------------- ------------
  Tn5ext               5′-TATTGCTGAAGAGCTTGGC-3′              This study
  [arb]{.smallcaps}1   5′-GGCCACGCGTCGACTAGTACNNNNNNNNNN-3′   \[[@R21]\]
  Tn5int               5′-GCATCGCCTTCTATCGCCTTC −3'           This study
  [arb]{.smallcaps}2   5′-GGCCACGCGTCGACTAGTAC-3′             \[[@R21]\]
  prMutS1              5′-TATCAAAACGACATCAAAGCCGTTGC-3′       This study
  prMutS2              5′-TTAGTGGAACTGGTTCATGGTGTT-3′         This study
  prMutD1              5′-ATGCCCGCTCAACCGACC-3′               This study
  prMutD2              5′-TCACGGGGTCCTTCT-3′                  This study
  prRegBam             5′-AGGATCCGGAGGGTTGAGCAGACAT-3′        This study
  prRegEco             5′-AGAATTCCGCGTGTGCTCGTTGCG-3′         This study
  prLigBam             5′-AGGATCCTTGGGGTTGCTGCTGTCT-3′        This study
  prLigEco             5′-AGAATTCGGAGGGTTGAGCAGACAT-3′        This study
  prIsocitBam          5′-AGGATCCAGAACGGTTTCGCCCG-3′          This study
  prIsocitEco          5′-AGAATTCTTTCAGAGCGGCAACGG-3′         This study
  KnRegXho             5′-ACTCGAGGAAATCTCGCGGATCACCG-3′       This study
  KnRegXba             5′-ATCTAGATCATCGACAGCCAGTCACGG-3′      This study
  prDehEco             5′-AGAATTCACCGGTGGGTTGTCGACT-3′        This study
  prDehBam             5′-AGGATCCGCCAGCAACTGGTGTT-3′          This study
  ΔReg mut             5′-ATCACCGACCAATCGCCTTCG-3′            This study
  prRegSensBam         5′-AGGATCCGAGCAGGCGATCCAG-3′           This study

Construction of genomic transposon mutant bank and isolation of mutants unable to use azelaic acid as the carbon source {#s2-3}
-----------------------------------------------------------------------------------------------------------------------

Plasmid pCRS530 \[[@R19]\] carrying the mTn5-GNm (containing *gus*A-P~tac~- *npt*II-T~trpA~) was harboured in *[E. coli](http://doi.org/10.1601/nm.3093)* S17-1 \[[@R20]\] and biparental conjugation was performed in order to conjugate this plasmid into *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128. The *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* strains harbouring the mTn5-GNm into the chromosome were selected on Luria--Bertani agar plates containing Km, Nf and Nx and incubated for 18 h at 28 °C. Several conjugations were performed in order to obtain approximately 30 000 independent transposon insertions. Among these, 5000 mutants were screened for inability to grow on azelaic acid as the unique source of carbon and energy as follows: purified colonies from the mTn5-GNm genomic mutant bank were individually hand-picked and grown at 28 °C in replicate on minimal M9 plates containing glucose or azelaic acid (0.1%). Six *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 transposon mutants were selected, which were unable or could very poorly grow on azelaic acid as the sole carbon source, but retaining the ability to grow on glucose. Their azelaic acid utilization inability was confirmed by growth-curve studies.

Mapping of transposon insertion sites and complementation of genomic mutants {#s2-4}
----------------------------------------------------------------------------

The DNA sequence-flanking transposon mutants in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 were determined using an arbitrary PCR procedure, as previously described \[[@R21]\]. In this method, the DNA flanking the transposon insertion site was enriched in two rounds of amplification using primers specific to the ends of the mTn5-GNm element and primers of random sequence that annealed to chromosomal sequences flanking the transposon (for primers see [Table 2](#T2){ref-type="table"}). For mutants 9128E, 9128I, 9128F and 9128M the arbitrary PCR did not work and in order to identify/clone the chromosomal DNA flanking the transposon, digestion of genomic DNA was performed using restriction enzymes that do not recognize the DNA sequence of the transposon. Sequences obtained were subjected to homology searches using a local NCBI [blast]{.smallcaps} and also with the draft genome sequence of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 using a local [blast]{.smallcaps} algorithm confirming the genomic location of the mutants.

The mutants selected were complemented for their defective growth as follows: a genomic bank (cosmid library) of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 was constructed from a partial *EcoR*I digested *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genomic DNA and ligated into *Eco*RI digested cosmid pLAFR3. The concatameric constructs obtained were then packaged in phage lambda and transfected into *[E. coli](http://doi.org/10.1601/nm.3093)* HB101 cells using Gigapack III XL-4 packaging kit as recommended by the supplier (Stratagene-Agilent, Santa Clara, CA, USA). The genomic cosmid bank (7500 clones) was then screened for the gene of interest by hybridization with labelled DNA probes via colony blot hybridization and then verified by Southern analysis. Triparental matings between the *[E. coli](http://doi.org/10.1601/nm.3093)* HB1010 harbouring the cosmids selected were then conjugated into *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* using *[E. coli](http://doi.org/10.1601/nm.3093)* (pRK2013) as the helper strain.

Construction of the *azeR* mutant in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* {#s2-5}
-----------------------------------------------------------------------------------------

The *azeR* gene was inactivated in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* by amplifying via PCR an internal fragment of this gene using primers KnRegXho, KnRegXba listed in [Table 2](#T2){ref-type="table"}, cloning it in pKNOCK-Km \[[@R14]\] and then using it as a suicide plasmid via conjugation into *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* and selecting for insertion into the chromosome via homologous recombination. Km-resistant colonies were selected and the mutation of the targeted gene was confirmed by PCR.

Construction of plasmid gene promoter transcriptional fusions and *β*-galactosidase assays {#s2-6}
------------------------------------------------------------------------------------------

Transcriptional gene promoter activity studies of four promoters were performed in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)*. All constructs were made in the promoter probe vector pMP220 (Tc^R^); this is a plasmid which carries a promoterless *lacZ* gene and harbours the gene for tetracycline resistance. The four gene promoter regions were amplified from *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genomic DNA by using the primers listed in [Table 2](#T2){ref-type="table"}. The PCR amplified fragments were first cloned in pGEM-T Easy vector (Promega, Madison, WI, USA), sequenced and then excised as *Eco*RI/*Bam*HI fragments and cloned into the *Bgl*II/*Eco*RI sites in pMP220 obtaining pMPLig, pMPDeh, pMPIsocit, pMPReg. The constructs obtained are listed in [Table 1](#T1){ref-type="table"}.

*β-*galactosidase activities were determined essentially as described by Miller \[[@R16]\], with the modifications of Stachel *et al*. \[[@R22]\]. Each experiment was performed in triplicate. Gene promoter activity was determined by monitoring *β*-galactosidase activity in overnight cultures.

Draft genome sequence of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128. {#s2-7}
---------------------------------------------------------------------------------------

The genome sequence of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* was determined using a 250 bp paired-end library with the Illumina MiSeq sequencing system (University of Exeter Sequencing Service, Exeter, UK), which generated a total of 1 611 521 pairs of reads. Reads were assembled using SPAdes 3.9.03 \[[@R23]\] and the assembled sequence annotated using the NCBI Prokaryotic Genome Annotation Pipeline (PGAP). The *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 genome was also uploaded to the Integrated Microbial Genomes and Metagenomes (IMG/M) database and was automatically annotated using annotation pipeline IMG Annotation Pipeline v.4.16.6 \[[@R24]\]. The IMG genome ID is 2818991220.

Phylogenomics and AzeR homology analysis {#s2-8}
----------------------------------------

The amino acid sequence of AzeR of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 was queried by a [blast]{.smallcaps} search against all genomes in the US Department of Energy IGM/M database. All hits with greater than 60 % amino acid sequence identity were considered matches, and the species associated with each protein were taken for further comparison (*n*=192). Overall, our methodology compared sequence variation of 30 COGs (Clusters of Orthologous Genes) that are considered universal single-copy genes \[[@R25]\] to calculate the evolutionary distance between the species. Of the 192 species matches to AzeR, poor-quality genomes that were missing universal COGs as well as those with more than one copy of COG genes were removed from the analysis, with a final list of 143 genomes continuing in the pipeline. The COGs from each species were aligned using Clustal Omega, with FASTA-formatted output. Each alignment was trimmed using TrimAl v1.3 \[[@R26]\] using the 'automated1' flag. The COGs from each species were concatenated, all in the same order. These 143 sequences were inputted into fasttree 2.1.10 \[[@R27]\] using standard settings. The outputted tree was visualized using Interactive Tree of Life \[[@R28]\]. *[E. coli](http://doi.org/10.1601/nm.3093)* K12 MG1655 served as an outgroup. The same analysis was performed for comparison of relatives of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)*, using the marked COGs indicated in [Table 3](#T3){ref-type="table"}.

###### 

COGs used for *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* phylogenomics and AzeR phylogenetic tree

  -----------
  COG0012\*
  COG0016
  COG0018
  COG0049\*
  COG0052\*
  COG0080\*
  COG0081
  COG0087\*
  COG0088\*
  COG0090\*
  COG0092\*
  COG0093\*
  COG0094\*
  COG0096\*
  COG0097\*
  COG0098\*
  COG0103\*
  COG0124
  COG0172
  COG0185\*
  COG0186\*
  COG0197\*
  COG0200\*
  COG0201\*
  COG0202\*
  COG0215
  COG0256\*
  COG0495\*
  COG0522\*
  COG0525\*
  -----------

Those with a star \*were used for Fig. 1. All COGs in this table were used for Fig. 6.

Nucleotide sequence accession numbers {#s2-9}
-------------------------------------

This whole-genome shotgun project of the *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* genome has been deposited at DDBJ/EMBL/GenBank under the accession number VASG00000000. Sequence reads are available in the Sequence Read Archive under accession number SRR9041312.

Results {#s3}
=======

Draft genome sequence of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 {#s3-1}
--------------------------------------------------------------------------------------

In order to begin to study the possible regulatory response in bacteria to azelaic acid, the genome of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 was sequenced, since this strain was isolated and reported as being able to efficiently degrade and utilize azelaic acid as the sole carbon and energy source \[[@R10]\]. We obtained a total of 1 611 521 pairs of reads, which we assembled into 15 contigs with a total length of 6 707 527 bp and an N~50~ length of 782 887. The G+C content was 65.5 %, similar to that of other sequenced *[Pseudomonas](http://doi.org/10.1601/nm.2552)* genomes. Automated annotation of the *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* draft genome sequence assigned a total of 6055 candidate protein-coding genes. A total of 6 rRNA and 60 tRNA genes were also identified. The phylogenetic analysis ([Fig. 1b](#F1){ref-type="fig"}) showed that the strain of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 is most closely related to other strains of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* spp. such as *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* strain TX1, which was isolated from a rice field\'s recycle drainage, and has the ability to metabolize industrial surfactants that pollute aquatic environments \[[@R29]\].

Isolation and identification of genomic transposon mutants affected in azelaic acid degradation in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 {#s3-2}
----------------------------------------------------------------------------------------------------------------------------------------------------------------

We confirmed that *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* strain DSM 9128 can utilize azelaic acid as the sole carbon source since it could grow well in M9 minimal medium supplemented with 0.2 % w/v of azelaic acid as the unique source of carbon and energy. In order to identify the regulatory mechanism(s) in response to azelaic acid a *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* transposon mutant bank was constructed and screened for inability to grow on azelaic acid. In total, 5000 transposon mutants of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 were screened on solid media for their ability to grow on azelaic acid. Six mutants, designated 9128D, 9128E, 9128F, 9128I, 9128M, 9128S, which could not grow on minimal media with azelaic acid but grew on glucose as the sole carbon source were isolated. Mutants 9128S, 9128E, 9128I, and 9128F displayed a complete lack of growth phenotype whereas 9128D and 9128M were characterized by very slow growth compared to the wild-type strain. Growth curves were determined for mutants 9128M, 9128D and 9128S, which represent the three different classes along with the wild-type strain and the results show significant differences in growth ([Fig. 2a](#F2){ref-type="fig"}). It was therefore concluded that the identified mutants most likely had transposon insertions in loci, which were involved in the utilization of azelaic acid as the unique source for carbon in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* 9128.

![Growth rates of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 and mutant derivatives. (a) Growth rates of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 and mutant derivatives 9128M, 9128D, 9128S and 9128AZER grown in M9 minimal media with azelaic acid 0.2 % as the unique source of carbon and energy. (b) Growth rates of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 and mutant derivative 9128D and 9128D (pCos5) grown in M9 minimal media with azelaic acid 0.2 % as the unique source of carbon and energy. (c) Growth rates of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 and mutant derivative 9128M and 9128M (pCos1) grown in M9 minimal media with azelaic acid 0.2 % as the unique source of carbon and energy. (d) Growth rates of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 and mutant derivative 9128S and 9128S (pCosM) grown in M9 minimal media with azelaic acid 0.2 % as the unique source of carbon and energy.](mic-166-73-g002){#F2}

Characterization of the six transposon genomic mutants of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 affected in azelaic acid utilization {#s3-3}
------------------------------------------------------------------------------------------------------------------------------------------------------------

In order to identify the genes in the transposon insertion mutants that were defective in azelaic acid utilization, the location of the transposon in the mutants was mapped. Mutant 9128D had the transposon inserted in a gene encoding for an acyl-CoA ligase for long-chain fatty acids whereas mutant 9128S had the insertion in the isocitrate lyase gene ([Fig. 3](#F3){ref-type="fig"}). Mutants 9128E, 9128F, 9128I and 9128S carried the transposon in an isocitrate lyase gene but in different positions within its nucleotide sequence thus representing independent insertion mutants and strengthening the association of this gene with azelaic acid utilization. The transposon in the mutant 9128M was located in a locus encoding an acyl-CoA dehydrogenase, which was genetically adjacent to the acyl-CoA ligase, which had the transposon insertion of mutant 9128D ([Fig. 3b](#F3){ref-type="fig"}). The acyl-CoA dehydrogenase is part of an operon of three genes and the acyl-CoA ligase is part of an operon with the enoyl-CoA hydratase. The two operons are separated by a transcriptional regulator that belongs to the IclR family ([Fig. 3](#F3){ref-type="fig"}, designated as AzeR).

![Map of genetic loci harbouring the Tn5 mutations of the identified mutants. (a) Genetic locus harbouring the Tn5 mutants 9128D and 9128M, which were affected in growth in azelaic acid as the unique source of carbon. The position of the Tn5 is shown by a red triangle, the positions of the ORFs are shown and the encoding gene product is described below the open arrows. The AzeR IclR-family transcriptional regulator is indicated by a black-filled arrow. (b) Genetic locus harbouring the Tn5 mutants 9128E, 9128I and 9128S, which were affected in growth in azelaic acid as the unique source of carbon. The position of the Tn5 is shown by a red triangle, the positions of the ORFs are shown and the encoding gene product is described below the open arrows.](mic-166-73-g003){#F3}

Next, we tried to complement the mutations. A pLAFR3 cosmid genomic library of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM9128 was constructed as described in Methods. The cosmid colonies harboured in *[E. coli](http://doi.org/10.1601/nm.3093)* were then hybridized with a labelled DNA probe of the *iclR* regulatory gene located in between the two operons and four cosmids were identified and designated pCos1, pCos5, pCosM and pCosG, which contained the genomic loci harbouring the genes of interest.

It was then determined that pCos5 complemented the defect of growth phenotype on azelaic acid displayed by the acyl-CoA ligase 9128D mutant ([Fig. 2b](#F2){ref-type="fig"}); pCos1 and pCos5 could significantly complement (*P*\<0.001) the acyl-CoA dehydrogenase 9128M ([Fig. 2c](#F2){ref-type="fig"}) mutant and pCosM complemented the isocitrate lyase mutant 9128S ([Fig. 2d](#F2){ref-type="fig"}). It was therefore concluded that the phenotypes observed by the transposon mutants were due to the loci in which the insertion occurred.

AzeR is an IclR family repressor that responds to azelaic acid {#s3-4}
--------------------------------------------------------------

The main aim of this work was to study the regulation of genes in response to azelaic acid. The gene promoters regulating the transcription of the operons harbouring the acyl-CoA ligase and acyl-CoA dehydrogenase and the gene promoter of the isocitrate lyase were cloned in a promoter probe vector, which harboured a promoterless *β*-galactosidase (*lacZ*) gene. *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM9128 harbouring the plasmids constructs were assayed for promoter activity and it was established that the presence of 0.1 % w/v of azelaic acid dramatically induced the transcription from the promoters controlling the expression of the operons containing the acyl-CoA dehydrogenase and of the acyl-CoA ligase ([Fig. 4a](#F4){ref-type="fig"} and b). If other very closely related dicarboxylic acids were used instead, like suberic, sebacic and adipic acid ([Fig. 1a](#F1){ref-type="fig"}), no activation of the two promoters was observed indicating that the response was specific for azelaic acid ([Fig. 5a](#F5){ref-type="fig"}). Intergenically located between the operons, which harbour two loci identified in this study, is a locus encoding a putative *iclR* family transcriptional regulator ([Fig. 3a](#F3){ref-type="fig"}); it was therefore postulated that this regulator, which we designated as AzeR, might play a role in their regulation. A genomic knock-out mutant of this regulator was constructed, which we designated as strain 9128AZER. The gene promoter assays were then performed in the wild-type strain and the *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* 9128AZER establishing that the two promoters displayed significantly high constitutive expression even in the absence of azelaic acid ([Fig. 4a](#F4){ref-type="fig"} and b). It was therefore concluded that AzeR behaved as a negative transcriptional regulator, which dramatically de-repressed transcription of the two operons in the presence of azelaic acid. We also determined whether AzeR was autoregulating its own expression; gene promoter studies showed that it was negatively autoregulating its own transcription in response to azelaic acid ([Fig. 4c](#F4){ref-type="fig"}). The isocitrate lyase promoter on the other hand did not respond to azelaic acid and it was active and constitutive in all the conditions tested ([Fig. 4d](#F4){ref-type="fig"}). With the aim to understand at which concentration of azelaic acid the AzeR de-repressed the operon promoters, a calibration curve in response to different concentrations of azelaic acid was determined. The results obtained revealed that the genes involved in the azelaic acid degradation pathway were strongly regulated by the azelaic acid, starting from a concentration higher than 10 µM ([Fig. 5b](#F5){ref-type="fig"}).

![Gene-promoter studies in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)*. (a) *β*-galactosidase assays in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128, 9128AZER and 9128AZER(pAzeR) harbouring the gene promoter regulating the transcription of the acyl-CoA ligase (inactivated by transposon 9128D) fused to a promoterless *lacZ* in plasmid construct pMPLig. All strains were grown in minimal M9 media contain either 0.2 % glucose (M9 GLU) or azelaic acid (M9 AZ) as the unique source of carbon and energy. (b) *β*-galactosidase assays in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128, 9128AZER and 9128AZER(pAzeR) harbouring the gene promoter regulating the transcription of the acyl-CoA dehydrogenase (inactivated by transposon 9128M) fused to a reportless *lacZ* in plasmid construct pMPDeh. All strains were grown in minimal M9 media containing either 0.2 % glucose (M9 GLU) or azelaic acid (M9 AZ) as the unique source of carbon and energy. (c) *β*-galactosidase assays in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128, 9128AZER and 9128AZER(pAzeR) harbouring the gene promoter regulating the transcription of the *azeR* regulator gene fused to a reportless *lacZ* in plasmid construct pMPazeR. All strains were grown in minimal M9 media containing either 0.2 % glucose (M9 GLU) or azelaic acid (M9 AZ) as the unique source of carbon and energy. (d) *β*-galactosidase assays in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128, 9128AZER and 9128AZER(pAzeR) harbouring the gene promoter regulating the transcription of the isocitrate lyase fused to a reportless *lacZ* in plasmid construct pMPIsocit. All strains were grown in minimal M9 media containing either 0.2 % glucose (M9 GLU) or azelaic acid (M9 AZ) as the unique source of carbon and energy. The data are means of three repeats, error bars indicate standard deviation and and *t*-test was performed (\**P*\<0.05, \*\**P*\<0.01, \*\*\**P*\<0.001), using Prism 7 (GraphPad Software).](mic-166-73-g004){#F4}

![Gene promoter studies in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* grown in the presence of different dicarboxylic acids. (a) *β*-galactosidase assays in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 harbouring the gene promoter regulating the transcription of the acyl-CoA ligase (fused to a reportless *lacZ* in plasmid construct pMPLig) grown in minimal M9 medium in the presence of different dicarboxylic acids at a concentration of 0.2 %; C4 is succinic acid, C4 subst. is oxaloacetic acid, C5 is glutaric acid, C6 is adipic acid, C8 is suberic acid, C9 is azelaic acid, C10 is sebacic acid and C12 is dodecanedioic acid. (b) *β*-galactosidase assays in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 harboring the gene promoter regulating the transcription of the acyl-CoA ligase (inactivated by transposon 9128D) fused to a reportless *lacZ* in a plasmid construct pMPLig in M9 minimal medium with azelaic acid as the unique source of carbon and energy in different concentrations.](mic-166-73-g005){#F5}

AzeR is a novel bacterial regulator {#s3-5}
-----------------------------------

AzeR belongs to the IclR family of bacterial transcriptional regulators, which consists of both activators and repressors \[[@R30]\]. These regulators are approximately 250 amino acids long and possess two domains: a DNA-binding helix-turn-helix domain at their *N*-terminus and an effector/inducer domain at their C-terminus. IclR-family regulators have been implicated in the regulation of many different loci, which are involved in degradation of aromatic compounds, in the glyoxylate cycle, in multidrug resistance and in virulence \[[@R30]\]. In order to establish the presence of orthologues and homologues of AzeR in bacteria, the amino acid sequence was compared to available sequences in the databases. [Fig. 6](#F6){ref-type="fig"} depicts the highly similar proteins (\>60 %) in bacteria showing that although AzeR-like proteins occur almost exclusively in *[Pseudomonas](http://doi.org/10.1601/nm.2552)* species, they are also found in *[Marinobacter](http://doi.org/10.1601/nm.2853)* species. To our knowledge, none of these proteins has yet been studied or their biological role reported making AzeR a new IclR-family regulator in bacteria. The *[Pseudomonas](http://doi.org/10.1601/nm.2552)* species that harbour an AzeR-like regulator are known to be able to degrade many compounds for their energy and carbon source and many are plant-associated. *[Marinobacter](http://doi.org/10.1601/nm.2853)* are gamma-Proteobacteria that live in seawater and have been isolated and studied for their ability to degrade hydrocarbons \[[@R31]\]. Moreover, azelaic acid is enriched in marine aerosols of high biological activity in the ocean vs. regions of low biological activity (Bikkina *et al*. \[[@R32]\]). *[Marinobacter](http://doi.org/10.1601/nm.2853)* in these aquatic regions may be responding to azelaic acid thorough its AzeR regulon. AzeR therefore is confined to a few proteobacterial genera that are likely to be able to degrade and utilize azelaic acid as the sole source of carbon and energy.

![AzeR homologies in other bacteria. Phylogenetic tree showing organisms with a strong (\>60 % amino acid sequence identity) AzeR homologue. The evolutionary distance between the organisms was determined using concatenated COGs listed in [Table 3](#T3){ref-type="table"}.](mic-166-73-g006){#F6}

Discussion {#s4}
==========

The bacterial strain *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 was used as a model for the study of genes necessary for the catabolism of azelaic acid. In bacteria, which utilize fatty acids as a carbon source, azelaic acid is relatively less preferred over other carbon sources \[[@R7]\]. A genetic screen using a genomic transposon library resulted in the identification of a set of mutants, which were unable, or significantly altered, in the ability to utilize azelaic acid as a unique source of carbon and energy. Mutants 9128D and 9128M (in an acyl-CoA ligase, an acyl-CoA dehydrogenase genes, respectively) that showed slow growth phenotype can possibly have mutated a functionally redundant gene for the degradation of azelaic acid. The total absence of growth on azelaic acid observed in mutants 9128S, 9128E, 9128I and 9128F, which are null mutants in an isocitrate lyase gene, indicated that the locus of the transposon insertion is absolutely necessary for this catabolic pathway.

The acyl-CoA ligase catalyses the pre-step reaction of *β*-oxidation by joining the acyl-CoA to the fatty acids that activates the breakdown of the molecule \[[@R33]\]. The acyl-CoA dehydrogenase enzyme on the other hand, catalyses the second step of *β*-oxidation of fatty acids and in particular transfers electrons from the substrate to the cofactor FAD in order to form a molecule of enoyl-CoA acting as a substrate for subsequent reactions \[[@R34]\]. Inactivation of these two genes in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 results in the bacterium being significantly less efficient in the catabolism of azelaic acid (C9 dicarboxylic acid), as well as suberic acid (C8) and sebacic acid (C10). From *in silico* analysis it was observed that genes encoding for the acyl-CoA ligase and for the acyl-CoA dehydrogenase are part of a cluster of genes encoding for other enzymes in the pathway of *β*-oxidation ([Fig. 3a](#F3){ref-type="fig"}). Specifically, these genes are organized into two operons: one is composed of three genes, an enoyl-CoA hydratase, an acyl-CoA dehydrogenase and a long-chain fatty acid acyl-CoA dehydrogenase (this last locus identified and studied in this work) and the second one is composed of two genes, the long-chain fatty acid acyl-CoA ligase identified in this study and an enoyl-CoA hydratase. These two operons are separated by a nucleotide sequence encoding for a transcriptional regulator belonging to the IclR family.

Previous studies made in *[E. coli](http://doi.org/10.1601/nm.3093)* and in *[Rhodopseudomonas palustris](http://doi.org/10.1601/nm.1494)* suggest the presence in the genomes of these bacteria of a number of enzymes involved in *β*-fatty acid oxidation; some are organized in an operon called *fad* operon in *[E. coli](http://doi.org/10.1601/nm.3093)* and the *pim* operon in *[R. palustris](http://doi.org/10.1601/nm.1494)* \[[@R8]\] both encoding for five enzymes involved in long and medium fatty acid degradation, reminiscent of the ones found in *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* DSM 9128 reported here. The organization of the loci and the sequence of individual genes do not show however a particularly high homology.

An *iclR* family regulator is located adjacent to the two operons involved in azelaic acid degradation (see above). We established that this regulator in *P. nitroreducens,* which we designated as AzeR, is involved in the stringent negative regulation of the two adjacent operons in response to azelaic acid. As far as we know, this is the first transcriptional regulator responding to azelaic acid, an important molecule involved in plant--microbe interactions. The fact that it responds very poorly to other closely related dicarboxylic acids ([Fig. 5](#F5){ref-type="fig"}) indicates that *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* has evolved the ability to specifically respond to the presence of azelaic acid and via AzeR de-repress the transcription of the two operons ([Fig. 5](#F5){ref-type="fig"}). AzeR also negatively auto-regulates its own transcription in response to azelaic acid ([Fig. 5](#F5){ref-type="fig"}); auto-regulation of transcription factors is a common trait in bacteria since it allows a more efficient response to the stimulus affecting the activity of the regulator. Interestingly AzeR homologes are only found in *[Pseudomonas](http://doi.org/10.1601/nm.2552)* (homology from 70--99 %) and *[Marinobacter](http://doi.org/10.1601/nm.2853)* sp. (homology from 60--70 %); it remains to be established if these species are able to utilize azelaic acid as the carbon source and if AzeR is the regulator which is responsible for the regulation of the genetic loci. The ability to respond to azelaic acid by these two genera indicates that they encounter it in the environment that they colonize. As AzeR is an IclR-family regulator, which specifically responds to azelaic acid and de-represses the transcriptional regulation, it is likely that it binds azelaic at its C-terminus and this results in AzeR being unable to bind to DNA through its helix-turn-helix N-terminus DNA-binding domain, thereby allowing transcription of the target genes \[[@R30]\]. Biochemical and structural studies are needed to confirm this likely model of regulation by AzeR.

The isocitrate lyase mutant is characterized by the total lack of growth of *[P. nitroreducens](http://doi.org/10.1601/nm.2652)* with azelaic acid as the unique source of carbon and energy. This is an enzyme of the glyoxylate cycle, which acts within the TCA cycle allowing it to skip the last steps of the TCA cycle \[[@R36]\]. The growth on fatty acids requires the induction of this glyoxalate bypass since it prevents the loss of carbon as CO~2~. The absolute requirement of isocitrate lyase indicates that azelaic degradation also occurs via the glyoxalate shunt. The isocitrate lyase is not regulated by AzeR or induced by the presence of azelaic acid and displayed constitutive expression during all the conditions tested in this study.

In summary, this study has revealed the presence of an IclR family repressor in bacteria, designated as AzeR, which responds to azelaic acid and de-represses transcription of two operons involved in the degradation of azelaic acid. AzeR is widespread only in *[Pseudomonas](http://doi.org/10.1601/nm.2552)* spp. and in *Marinebacter* indicating that these two genera have evolved the ability to sense the presence and respond to azelaic acid and regulate gene expression and this might be an important train in plant--microbe interaction.

This work received no specific grant from any funding agency.

C.B. and S.G.J. are beneficiary of an ICGEB fellowship. We thank ICGEB for supporting this work. The bacterial *P. nitroreducens* strain DSM9128 sequenced in this work was obtained from the DSMZ collection and was isolated by Janot-Bassalik and Wright as cited in this work.

The authors declare that there are no conflicts of interest.

Abbreviations: AZ, azelaic acid; GLU, glucose.

Edited by: D. Grainger and D. Lee

[^1]: These authors contributed equally to this work
