Vanillic acid and methoxyhydroquinone production from guaiacyl units and related aromatic compounds using Aspergillus niger cell factories

Background The aromatic compounds vanillin and vanillic acid are important fragrances used in the food, beverage, cosmetic and pharmaceutical industries. Currently, most aromatic compounds used in products are chemically synthesized, while only a small percentage is extracted from natural sources. The metabolism of vanillin and vanillic acid has been studied for decades in microorganisms and many studies have been conducted that showed that both can be produced from ferulic acid using bacteria. In contrast, the degradation of vanillin and vanillic acid by fungi is poorly studied and no genes involved in this metabolic pathway have been identified. In this study, we aimed to clarify this metabolic pathway in Aspergillus niger and identify the genes involved. Results Using whole-genome transcriptome data, four genes involved in vanillin and vanillic acid metabolism were identified. These include vanillin dehydrogenase (vdhA), vanillic acid hydroxylase (vhyA), and two genes encoding novel enzymes, which function as methoxyhydroquinone 1,2-dioxygenase (mhdA) and 4-oxo-monomethyl adipate esterase (omeA). Deletion of these genes in A. niger confirmed their role in aromatic metabolism and the enzymatic activities of these enzymes were verified. In addition, we demonstrated that mhdA and vhyA deletion mutants can be used as fungal cell factories for the accumulation of vanillic acid and methoxyhydroquinone from guaiacyl lignin units and related aromatic compounds. Conclusions This study provides new insights into the fungal aromatic metabolic pathways involved in the degradation of guaiacyl units and related aromatic compounds. The identification of the involved genes unlocks new potential for engineering aromatic compound-producing fungal cell factories. Supplementary Information The online version contains supplementary material available at 10.1186/s12934-021-01643-x.

Currently, less than 1% of the produced vanillin is derived from natural sources, while the majority of vanillin is obtained through chemical synthesis, mainly from guaiacol [1,7]. Vanillin obtained through chemical synthesis is considered "artificial" based on European regulations. Therefore, new strategies and methods are being developed to obtain "natural" vanillin through biosynthesis using microorganisms [9,10].
The metabolism of vanillin and its derivatives by microorganisms, especially bacteria, has been well studied [11]. Vanillin is toxic in low concentrations for many microorganisms and therefore the ability to degrade vanillin is essential for microorganisms that live in natural habitats or are used in industrial processes where vanillin is present in significant amounts [12,13]. Several vanillin metabolic pathways have been described in microorganisms [11]. In bacteria, vanillin is converted to vanillic acid and this reaction is catalyzed by vanillin dehydrogenase, followed by decarboxylation to guaiacol by vanillic acid decarboxylase [14]. In addition to oxidative routes, vanillic acid can also be demethylated to protocatechuic acid by vanillate-o-demethylase oxidoreductase or hydroxylated to methoxyhydroquinone by vanillate hydroxylase [15][16][17]. The latter appears to be uncommon for bacteria [17]. All these conversions have also been observed in filamentous fungi, but in contrast bacteria the conversion towards methoxyhydroquinone appears to be common in fungi [15].
Due to the extensive study of the vanillin metabolic pathway in bacteria, many vanillin-producing bacterial systems have been described [9,18]. However, only one strategy producing a considerable amount of vanillin using filamentous fungi has been reported [10,19,20]. In this method Aspergillus niger is used to convert ferulic acid to vanillic acid, which is further converted by Pycnoporus cinnabarinus to vanillin. However, this method is not efficient since both fungi can also convert vanillic acid to methoxyhydroquinone. In order to engineer efficient fungal cell factories, the metabolic pathway genes need to be identified. Vanillate hydroxylase from the fungi Phanerochaete chrysosporium and Sporotichum pulverulentum has been characterized, but the gene encoding vanillate hydroxylase remains to be identified [21][22][23].
It has been shown that A. niger converts coniferyl alcohol to ferulic acid, which is then converted to vanillic acid and further to methoxyhydroquinone [15,[24][25][26]. Recently, we observed that ferulic acid was converted through the CoA-dependent β-oxidative pathway to vanillic acid [26]. Deletion of the CoA-dependent β-oxidative genes did not result in abolished growth on ferulic acid. Therefore, it is possible that another pathway is present in A. niger in which vanillin is an intermediate [24,25]. In this study, we aimed to identify the genes involved in the vanillin and vanillic acid metabolic pathway of A. niger. Whole-genome transcriptome data of A. niger N402 upon transfer to coniferyl alcohol, ferulic acid, vanillic acid and veratric acid was used to identify genes involved in the vanillic acid metabolic pathway. With the obtained knowledge we created two fungal cell factories that can produce vanillic acid and methoxyhydroquinone from guaiacyl lignin units and related aromatic compounds.

Results
Identification of candidate genes encoding putative vanillin dehydrogenase (vdhA), vanillate hydroxylase (vhyA), methoxyhydroquinone 1,2-dioxygenase (mhdA) and 4-oxo-monomethyl adipate esterase (omeA) To identify candidate genes involved in the metabolism of ferulic acid and vanillic acid, whole-genome transcriptome data of A. niger N402 grown on coniferyl alcohol, vanillic acid and veratric acid were generated. Veratric acid was selected since it was observed to be converted through the same pathway as vanillic acid [11,27]. Whole-genome transcriptome data from A. niger N402 grown on ferulic acid and no carbon source was obtained from Lubbers et al. [25]. Transcriptome data from A. niger N402 grown in minimal media (MM) without carbon source was used as control.
In total, 69 genes were upregulated (FPKM ≥ 10, fold change (log2) ≥ 2, p-value ≤ 0.01) in all conditions compared to the no carbon source control (Fig. 1), but no genes annotated as aldehyde dehydrogenases were upregulated in all four tested conditions. As no vanillin dehydrogenase has been identified in fungi, we performed a BLASTP analysis using the amino acid sequence of vanillin dehydrogenase (Vdh) from Pseudomonas sp. (Uniprot O05619) [28] as a query against the A. niger genome, resulting in 29 hits. A cut-off E-value of e ≤ − 60 was used to remove insignificant hits resulting in three genes (NRRL3_10496, NRRL3_6772 and NRRL3_3887). Only NRRL3_3887 was upregulated by coniferyl alcohol, while the other two were not upregulated in the tested conditions. Therefore, NRRL3_3887 was selected as a putative vdhA (Fig. 2a).
Allthough no vanillate hydroxylase-coding genes have been identified, a vanillate hydroxylase was partially characterized in S. pulverulentum showing that this enzyme needs FAD and NADH or NADPH as cofactor [22]. Therefore, we selected FAD-binding domain containing genes upregulated in coniferyl alcohol, ferulic acid and vanillic acid. In our dataset, only one FAD-binding domain containing gene (NRRL3_9897), annotated as a putative salicylate hydroxylase, was upregulated in all conditions and was selected as a putative vhyA (Fig. 2b).
Within the top 10 upregulated genes, a gene annotated as homogentisate 1,2-dioxygenase (NRRL3_10111) was upregulated under all tested conditions (Fig. 2c) and was selected as a putative mhdA. In addition, we observed that a putative carboxylesterase gene (NRRL3_10110) neighboring mhdA, was also highly upregulated under the same conditions (Fig. 2d) and was selected as a putative 4-oxo-monomethyl adipate esterase (omeA).

Deletion of the putative vdhA, vhyA mhdA and omeA genes result in reduced growth on vanillin and its metabolites
To verify if the candidate genes encode enzymes involved in the vanillin and vanillic acid metabolic pathway, deletion mutants of the candidate genes were made and screened for phenotypes on 14 aromatic compounds. Deletion of vdhA resulted in abolished growth on vanillin indicating that this gene is involved in the conversion of vanillin (Fig. 3). In addition, reduced growth on p-hydroxybenzaldehyde and protocatechuic aldehyde was also observed. Deletion of vhyA, mhdA and omeA resulted in reduced growth on coniferyl alcohol, ferulic acid and vanillic acid and abolished growth on vanillin, indicating that vhyA, mhdA and omeA are involved in the metabolic pathway. In addition, deletion of mhdA resulted in a red-brownish coloration of the medium when grown on ferulic acid and vanillic acid, which indicates that a compound, presumably methoxyhydroquinone, accumulates. After 10 days of growth, ΔmhdA appears to have slightly reduced growth and spore formation on vanillyl alcohol and veratric acid (Fig. 3). In addition, a light redbrownish coloration in the media was observed when grown on vanillyl alcohol indicating that it is converted towards methoxyhydroquinone. Deletion of vhyA, mhdA or omeA did not result in a phenotype on benzoic acid, caffeic acid, cinnamic acid, p-coumaric acid, p-hydroxybenzoic acid, p-hydroxybenzaldehyde, protocatechuic aldehyde or protocatechuic acid, indicating that these genes are not involved in the conversion of these aromatic compounds. To confirm the enzymatic function of VdhA, VhyA and MhdA, E. coli strains that produce VdhA, VhyA or MhdA were created. All three enzymes were isolated, purified, and visualized on SDS-PAGE by Coomassie brilliant blue staining (Additional file 1: Fig. S1) and Western blotting using a monoclonal antibody raised against the histidine-tag (Fig. 4a). The detected size of VhyA-His and MhdA-His corresponds with the expected masses of 48.1 and 55.5 kDa, respectively. However, the expected mass of VdhA (52.2 kDa) did not correspond with the slightly higher mass observed (Fig. 4a). Assays were performed and analyzed using HPLC, demonstrating that all enzymes had activity on their corresponding substrate (Fig. 4b).
The mixture containing MhdA showed a clear reduction in methoxyhydroquinone content and the formation of an unknown compound was detected (Fig. 5a). This compound was identified as 4-hydroxy-6-methoxy-6-oxohexa-2,4-dienoic acid using Orbitrap LC-HRMS. Identification of the detected compounds was carried out with Orbitrap Fusion high-resolution MS with both protonated and deprotonated molecules (Fig. 5b). The product peak at 3.76 min retention time was consistent with the expected product 4-hydroxy-6-methoxy-6-oxohexa-2,4-dienoic acid which showed an m/z

Deletion mutants accumulate vanillic acid, methoxyhydroquinone and 4-oxo-monomethyl adipate
To evaluate if the ∆mhdA deletion mutant can accumulate methoxyhydroquinone, an accumulation . All enzymatic reactions were incubated for 1 h at 30 °C experiment in liquid media with 5 mmol/L substrate was performed. In addition, ∆vdhA, ∆vhyA and ∆omeA deletion mutants were also tested for the accumulation of (aromatic) compounds. After 24 h incubation, supernatants were collected and analyzed with HPLC ( Table 1, Additional file 1: Fig. S2). The starting substrates (coniferyl alcohol, ferulic acid, vanillic acid and vanillin) were not detected in the supernatant of the reference strain indicating that these compounds were fully consumed within 24 h. Veratric acid and veratryl alcohol appeared to be slowly consumed compared to the other tested substrates.
Deletion of vdhA did not result in the accumulation of vanillin when grown on coniferyl alcohol or ferulic acid, which indicates these compounds are not converted to vanillin. Compared to the reference, vanillin consumption by ∆vdhA was clearly reduced (Table 1). However, a small amount of vanillyl alcohol was detected, which indicates that vanillin can also be converted to vanillyl alcohol. Trace amounts of vanillin were detected when ∆vdhA was grown on vanillyl alcohol, indicating that vanillyl alcohol can be reduced to vanillin.
Deletion of vhyA resulted in the accumulation of vanillic acid when grown on coniferyl alcohol, ferulic acid, vanillin, vanillyl alcohol, veratryl alcohol, veratric aldehyde and veratric acid, whereas vanillic acid itself was not consumed. Similar results were observed for ∆mdhA, which resulted in the accumulation of methoxyhydroquinone on these compounds and vanillic acid (Table 1, Additional file 1: Fig. S2, Additional file 2: Table S1). This indicates that the tested compounds are converted to vanillic acid and then further to methoxyhydroquinone.
Deletion of omeA resulted in the accumulation of an unknown product when grown on the tested compounds (Table 1). To investigate the accumulated product of ∆omeA, the unknown compound was extracted with ethyl acetate from the supernatant of ∆omeA grown on vanillic acid. The 1 H-NMR of the accumulated compound corresponded well with 4-oxo-monomethyl adipate (6-methoxy-4,6-dioxohexanoic acid) [29].

Phylogenetic study of VdhA and MhdA
To evaluate the presence of VdhA and MhdA in other fungi, a phylogenetic analysis was performed for each enzyme using selected ascomycete and basidiomycete genomes. VhyA was not analyzed since it was already included in the phylogenetic tree made for protocatechuic acid hydroxylase (PhyA, NRRL3_4659) [30].
Many aldehyde dehydrogenases were obtained as BLAST hits, but only a few clustered with A. niger VdhA (Additional file 1: Fig. S3). Homologs were observed in the genomes of Aspergillus japonicus, Aspergillus nidulans, Aspergillus oryzae, Aspergillus fumigatus, Phaeomoniella chlamydospora and Talaromyces stipitatus. Most fungi contained two or three genes annotated as homogentisate 1,2-dioxygenases (Additional file 1: Fig. S4). Some exceptions were found in A. nidulans and Penicillium subrubescens, which both contained four homologs, whereas Trichoderma reesei contained only one homolog. MhdA homologs were found in several Eurotiomycetes (A. japonicus, A. nidulans, A. oryzae, A. fumigatus, P. subrubescens, P. chrysogenum and T. stipitatus), the Sordariomycete Fusarium graminearum, and the Dothideomycete Cochliobolus lunatus. All these fungi had a homolog clustering with the homogentisate

Discussion
In this study, we identified four genes involved in the degradation of guaiacyl units of lignin (G-units), such as ferulic acid, vanillin and vanillic acid. In addition, we observed the conversion of methoxyhydroquinone to 4-hydroxy-6-methoxy-6-oxohexa-2,4-dienoic acid, which was not described before. Next to that, we demonstrated that the deletion mutants of this study can be used as cell factories for the production of vanillic acid and methoxyhydroquinone.
Vanillate hydroxylase activity has been observed in many ascomycetes and basidiomycetes and is suggested to play a major role in the degradation of the lignin G-units [31][32][33][34][35]. The identification of vhyA is an important finding which can unlock new strategies to engineer efficient fungal cell factories. For example, vanillin can be produced from ferulic acid by fungi using a two-step bioconversion process. In this process, ferulic acid or ferulic acid derived from sugar beet pulp and rice bran oil, was converted to vanillic acid by A. niger and was further reduced to vanillin by P. cinnabarinus [10,19,20,36]. However, a majority of the vanillic acid was lost in this process, because it is converted to methoxyhydroquinone by both fungi. Blocking the hydroxylation of vanillic acid to methoxyhydroquinone by deleting vhyA in both species could improve the vanillin yield in this process. It has been shown that VhyA is closely related to PhyA [30] which is involved in hydroxylation of protocatechuic acid to hydroxyquinol. In addition, PhyA is also involved in the degradation of tannic acid and gallic acid [37]. However, it was suggested that gallic acid is converted to 2-carboxy-cis,cis-muconate by PhyA, which does not correspond with the function observed for these enzymes [38]. Therefore, it is more likely that PhyA converts gallic acid to 1,2,3,5-tetrahydroxybenzene. Despite this, it is possible that VhyA is also involved in the degradation of other aromatic compounds that were not tested in our growth assay.
Deletion of vdhA abolishes growth on vanillin and results in reduced growth on p-hydroxybenzaldehyde and protocatechuic aldehyde, indicating that vdhA is also involved in the benzoic acid metabolic pathway [30,39]. It has been shown that bacterial vanillin dehydrogenases were also able to convert p-hydroxybenzaldehyde and protocatechuic aldehyde [40][41][42]. Deletion of vdhA did not result in reduced growth on coniferyl alcohol or ferulic acid, which indicates that these aromatic compounds are not converted to vanillin, while deletion of vhyA results in reduced growth on ferulic acid and coniferyl alcohol. This shows that both ferulic acid and coniferyl alcohol are converted towards vanillic acid, as previously observed [24][25][26], and revealed that vanillin is not an intermediate of the ferulic acid metabolic pathway. Deletion of the β-oxidative pathway genes did not result in abolished growth on ferulic acid [26]. A possible explanation is that alternative CoA-dependent β-oxidative genes are present in A. niger. In addition, the transcriptome data showed that vdhA was not induced by ferulic acid, but by coniferyl alcohol. It is possible that VdhA plays a role in the conversion of coniferyl aldehyde to ferulic acid. However, deletion of vdhA did not result in reduced growth on coniferyl alcohol, suggesting that other enzymes are also involved in this conversion.
Deletion of mhdA results in reduced growth on coniferyl alcohol, ferulic acid, vanillin and vanillic acid. In addition, deletion of mhdA resulted in accumulation of methoxyhydroquinone when grown on coniferyl alcohol, ferulic acid, vanillyl alcohol, vanillin, vanillic acid, veratryl alcohol, veratric aldehyde and veratric acid. Methoxyhydroquinone can be used as a building block to create epoxy resins or the thermoplastic poly(arylene ether sulfone) [5,7,43]. The enzyme assay with MhdA verified that the ring of methoxyhydroquinone is cleaved to 4-hydroxy-6-methoxy-6-oxohexa-2,4-dienoic acid, a reaction that has not been described before. At this moment, there are no descriptions of 4-hydroxy-6-methoxy-6-oxohexa-2,4-dienoic acid or its tautomer in literature and therefore the applications of this compounds are unexplored.
Phylogenetic analysis of VhyA and MhdA revealed that these enzymes are conserved in A. oryzae, Aspergillus flavus, Aspergillus terreus, A. fumigatus and A. nidulans, which were all able to convert vanillic acid to methoxyhydroquinone [15]. No homologs of Trichoderma species clustered with VhyA or MhdA, which correlates with the observation that Trichoderma converts vanillic acid to vanillin and vanillyl alcohol and not to methoxyhydroquinone. In A. nidulans, a homogentisate 1,2-dioxygenase (HmgA, AN1897) has been characterized [44]. However, BLASTP with the amino acid sequence of HmgA revealed that it shares only 50.8% identity with MhdA while NRRL3_9969 is 87% similar to HmgA. Alternative pathways were suggested in several Aspergillus species, in which vanillic acid was converted to protocatechuic acid or guaiacol [45], but these were not observed in A. niger [15]. Our study supports this observation since protocatechuic acid or guaiacol were not detected during the accumulation experiment with ΔvhyA or ΔmhdA on ferulic acid, vanillic acid and vanillin. Next to that, we previously showed that deletion of protocatechuic acid 3,4-dioxygenase (prcA) and/or hydroxyquinol 1,2-dioxygenase (hqdA, NRRL3_2644) does not result in reduced growth on ferulic acid, vanillin or vanillic acid [39]. Interestingly, hqdA was upregulated by coniferyl alcohol, ferulic acid and vanillic acid (Additional file 2: Table S2). However, deletion of hqdA did not result in a phenotype on ferulic acid nor did HqdA show activity on methoxyhydroquinone [39]. This suggests that methoxyhydroquinone can be converted to hydroxyquinol, which was also proposed for Paecilomyces variotii and S. pulverulentum [35,46]. However, the deletion of mhdA resulted in severely reduced growth on vanillin and vanillic acid and indicates that the conversion of methoxyhydroquinone to hydroxyquinol plays a minor role. The presence of this pathway could explain why deletion of mhdA does not result in complete accumulation of methoxyhydroquinone when grown on coniferyl alcohol, ferulic acid or vanillic acid. Based on previous observations and our transcriptome data, phenotypic profile and enzymatic assays, we suggest an updated version of the vanillic acid metabolic pathway of A. niger (Fig. 6).
In the A. niger genome sequence, the gene encoding for OmeA neighbors with mhdA and NRRL3_10109, annotated as an MFS-type tetracycline resistance protein, which was also upregulated by coniferyl alcohol, ferulic acid, vanillic acid and veratric acid (Additional file 2: Table S2). However, the involvement of this gene in the degradation of aromatic compounds remains to be studied.
The knowledge about the conversion of ferulic acid, vanillic acid and vanillin is important for new biotransformation strategies in fungi, for which better understanding of the aromatic metabolic pathways and the involved genes and enzymes involved is crucial. The identification and characterization of the novel A. niger enzymes VdhA, VhyA, MhdA and OmeA contributes greatly to a better understanding of these fungal aromatic metabolic pathways and revealed a novel metabolic pathway for the degradation of methoxyhydroquinone. In addition, we demonstrated that this knowledge can be applied to create fungal cell factories that can produce aromatic compounds such as vanillic acid and methoxyhydroquinone.

Growth conditions
Aspergillus niger strains used in this study are shown in Table 2. The fungi were grown on complete medium (CM) [51] agar (1.5% w/v) plates at 30 °C for 4 days. Spores were harvested with 10 mL N-(2-acetamido)-2-aminoethanesulfonic acid buffer. Minimal medium (MM) [51] agar (1.5% w/v) plates were inoculated with 10 3 freshly isolated spores. MM plates for growth profile experiments were supplemented with aromatic compounds as sole carbon source. Due to the toxicity of the aromatic compounds different concentrations were used for the growth profile, i.e. 2 mmol/L for benzoic acid, ferulic acid and vanillin, while 5 mmol/L was used for the remaining aromatic compounds. All aromatic compounds and chemicals were purchased from Sigma Aldrich.
Transfer experiments were performed as described previously [39]. Equal portions of mycelia were transferred to flasks containing 50 mL MM and 0.02% (w/v) coniferyl alcohol, vanillic acid or veratric acid. The cultures were incubated on a rotary shaker for 2 h at 30 °C, 250 rpm. Mycelia were harvested, dried between tissue paper to remove excess liquid and frozen in liquid nitrogen. Lubbers et al. Microb Cell Fact (2021) 20

RNA extraction and RNA sequencing
RNA was extracted as described previously [39]. The quality and quantity of RNA were determined by gel electrophoresis and RNA6000 Nano Assay using the Agilent 2100 Bioanalyzer (Agilent Technologies, Santa Clara, CA, USA). RNA sequencing was conducted by BGI Tech Solutions (Tai Po, Hong Kong) using the Illumina HiseqTM 2000 platform (Illumina Inc., San Diego, CA, USA). Transcriptome data of A. niger N402 grown on ferulic acid and no carbon source was obtained from GEO Accession number GSE135001 [25] and GSE13499 [39], respectively. Transcriptome analysis was performed as previously described [52]. The transcriptome data was stored at the NCBI Gene expression omnibus, under the GEO Accession number GSE154865.

Protoplast-mediated transformation of A. niger
All deletion strains were made through protoplast-mediated transformation of A. niger N593 ΔkusA by homologous recombination using deletion cassettes containing 900-1000 bp upstream and downstream of the target gene fused to the orotidine 5ʹ-phosphate decarboxylase (pyrG) gene from A. oryzae RIB40. The used primers are listed in Additional file 2: Table S3. The protoplastmediated transformation and purification of the selected transformants was performed as described previously [52]. Transformants were verified with diagnostic PCR using primers listed in Additional file 2: Table S3.

Accumulation of aromatic compounds
Transfer experiments were conducted as indicated above. Equal portions of mycelia were transferred to 250 mL flasks containing 50 mL MM and 5 mmol/L coniferyl alcohol, ferulic acid, vanillyl alcohol, vanillin, vanillic acid, veratryl alcohol, veratric aldehyde or veratric acid. The cultures were incubated on a rotary shaker for 24 h at 30 °C, 250 rpm. Supernatant samples were collected after 24 h and diluted 50 times with acetonitrile before HPLC analysis using the setup described previously [53].  [39]. After purification, a final concentration of 0.5 mmol/L FAD was added to VhyA.

Identification of 4-hydroxy-6-methoxy-6-oxohexa-2,4-dienoic acid
The MhdA reaction mixture, containing 0.1 mol/L phosphate buffer, pH 6.0, 250 µmol/L methoxyhydroquinone, 50 µmol/L FeSO 4 and 5 µL purified MhdA with a total volume of 5 mL, was analyzed with HPLC-MS before and after extraction with ethyl acetate. The extraction was performed with 10 mL reaction mixture by adjusting the pH to pH 2 and the water-soluble analytes were salted out with brine to ethyl acetate. The extraction was repeated three times and finally the product mixture was collected by evaporation of the solvent under vacuum in a rotary evaporator. The sample (5 mg) was dissolved in 1% formic acid in 1:10 acetonitrile:H 2 O as 1 mg/mL for high resolution LC-HRMS analysis. The LC-HRMS analysis was performed on Thermo Scientific Orbitrap Fusion mass spectrometer (San Jose, USA) connected to Thermo Scientific Dionex Ultimate 3000 ultrahigh performance liquid chromatograph (Germering, Germany). The sample injection volume was 5 μL. LC separation was done using Phenomenex Luna Ⓡ Omega Polar C18 (1.6 µm, 100 × 2.1 mm) column at 40 °C using a non-linear gradient (curve 7) of two mobile phases: 0.1% formic acid in water (A) and 0.1% formic acid in acetonitrile (B). The gradient with flow rate of 0.3 mL/min was run from 1% B at 0 min to 95% B at 15 min. After this the B eluent was kept 95% for 1 min, returned to 5% during 1 min and equilibrated at 5% for 3 min.
The ionization was done using thermal heated electrospray HESI probe both in the positive and negative ion mode. The instrumental parameters were set as follows: positive ion spray voltage 3500 V, negative ion spray voltage 2500 V, source temperature 300 °C, ion transfer tube temperature 350 °C, sheath gas 40, auxiliary gas 15 and sweep gas 0. Mass measurement was done with mass range m/z 70-350 using RF lens at 60% and quadrupole isolation (m/z 70-350) at resolution of 120,000. Mass accuracy of the instrument using external calibration for both positive and negative ion mode was specified to be ≤ 3 ppm.

Identification of 6-methoxy-4,6-dioxohexanoic acid
For the identification of 6-methoxy-4,6-dioxohexanoic acid, the accumulated product of ∆omeA grown on vanillic acid was isolated through ethyl acetate extraction. The organic phase was dried, dissolved in acetonitrile and analyzed for the presence of 6-methoxy-4,6-dioxohexanoic acid by HPLC. For NMR analysis, the aqueous sample was washed three times with CDCl 3 . The aqueous layer was concentrated to dryness in vacuo (50 °C) and redissolved in D 2 O. The 1 H-NMR spectrum was recorded on an Agilent MRF 400 equipped with a OneNMR probe and Optima Tune system. The signals of the major species correspond to 6-methoxy-4,6-dioxohexanoic acid (bisdeuterated at C5). The recorded spectrum is in good agreement with spectral data reported in the literature for the proposed structure [29].

Phylogenetic analysis
The amino acid sequences of VhyA and MhdA were used for BLASTP analyses on selected ascomycete and basidiomycete genomes as described previously (Additional file 2: Table S4) [54]. To reduce the amount of insignificant hits, a cutoff E-value of e−40 was used. Several amino acid sequences were curated manually or with the gene prediction software, Augustus (Additional file 2: Table S5) [55]. The bacterial enzymes, Vdh (Uniprot, O05619) from Pseudomonas sp. strain HR199 and LigV (Uniprot, A2PZP3) from Sphingomonas paucimobilis were added manually to the multiple alignment. The Maximum Likelihood, Neighbor Joining and Minimum Evolution trees were constructed using MEGA 7 with 500 bootstraps and complete deletion of gaps.