Skip to main content

ORIGINAL RESEARCH article

Front. Bioeng. Biotechnol., 31 July 2017
Sec. Bioprocess Engineering

Proteomic Insights on the Metabolism of Penicillium janczewskii during the Biotransformation of the Plant Terpenoid Labdanolic Acid

\r\nIsabel MartinsIsabel Martins1Adlia Varela,Adélia Varela1,2Luís M. T. FrijaLuís M. T. Frija3Mnica A. S. EstevoMónica A. S. Estevão3Sbastien PlanchonSébastien Planchon4Jenny RenautJenny Renaut4Carlos A. M. AfonsoCarlos A. M. Afonso3Cristina Silva Pereira*\r\nCristina Silva Pereira1*
  • 1Instituto de Tecnologia Química e Biológica António Xavier, Universidade Nova de Lisboa (ITQB NOVA), Oeiras, Portugal
  • 2Instituto Nacional Investigação Agrária e Veterinária, Oeiras, Portugal
  • 3Instituto de Investigação do Medicamento (iMed.ULisboa), Faculdade de Farmácia, Universidade de Lisboa, Lisboa, Portugal
  • 4Environmental Research and Innovation (ERIN) Department, Luxembourg Institute of Science and Technology, Belvaux, Luxembourg

Plant terpenoids compose a natural source of chemodiversity of exceptional value. Many of these compounds own biological/pharmacological activity, others are regarded as unique chemical skeletons for the synthesis of derivatives with improved properties. Functional chemical modification of terpenoids through biotransformation frequently relies on the use of Ascomycota strains, but information on major cellular responses is still largely lacking. Penicillium janczewskii mediates a stereo-selective hydroxylation of labdanolic acid (LA)—terpenoid found abundantly in Cistus ladanifer—producing 3β-hydroxy-labdanolic acid with yields >90%. Herein, combined analyses of mycelial and extracellular differential proteomes demonstrated that the plant terpenoid increased stress responses, especially against oxidative stress (e.g., accumulation of superoxide dismutase) and apparently altered mitochondria functioning. One putative cytochrome P450 monooxygenase differentially accumulated in the secretome and the terpenoid bioconversion was inhibited in vivo in the presence of a P450 inhibitor. The stereo-selective hydroxylation of the plant terpenoid is likely mediated by P450 enzymes, yet its unequivocal identity remains unclear. To the best of our knowledge, this is the first time that proteomics was used to investigate how a plant terpenoid impacts the metabolism of a filamentous fungus during its efficiently biotransformation. Our findings may encourage the development of new strategies for the valorization of plant natural resources through biotechnology.

Introduction

Terpenoids compose a dissimilar group of natural compounds (>25,000 compounds identified so far) that structurally are saturated and unsaturated cyclic and aliphatic hydrocarbons with variable degrees of oxygenation, including alcohols, aldehydes, ketones, and carboxylic acids (Tholl, 2015). Typically these molecules are distributed into subclasses according to the number of isoprene units in their structure. Terpenoids are often regarded as valuable natural resource with far-reaching applications, from flavors, fragrances, pigments, and agrochemicals to pharmaceutical drugs. In particular, among the subclass of the labdane-type diterpenes (20 carbon atoms) several show cytotoxicity (Yang et al., 2010), bactericidal (Singh et al., 2010; Ghosh and Rangan, 2014), antioxidant (Kapewangolo et al., 2015), antifungal (Hawas et al., 2015; Mendoza et al., 2015), and anti-inflammatory properties (Huang et al., 2015; Kapewangolo et al., 2015).

Due to its richness in phytochemicals, the genus Cistus has been proposed as a model for the biosynthesis of labdane-type diterpenes and also as a natural source of high value pharmacological products (Papaefthimiou et al., 2014). In particular, labdanolic acid (LA) (Figure 1), a diterpene extracted in large quantities from the Cistus ladanifer L. (“Rock-rose”) (Martins et al., 2014a), displays anti-inflammatory properties (Jayaprakasam et al., 2007) and is used as a precursor compound in the chemical synthesis of numerous valuable compounds, e.g., Ambrox® (the prototype compound of all ambergris odorants) (Bolster et al., 2001). This emphasizes well the importance to better understand how this subclass of plant terpenoids impacts on the metabolism of microbes that may be used for their modification. In opposition to the chemical modification of terpenoids, biotransformation renders products with a “natural” label and usually requires mild-reaction conditions avoiding also generation of toxic byproducts (Bicas et al., 2009). Therefore, biotransformation constitutes a reference method for producing new terpenoid derivatives with improved properties (Silva et al., 2013; Schrader and Bohlmann, 2015). The capacity of fungi, as well as bacteria, to transform terpenoids yielding stereo- and regio-selective products and/or compounds functionalized at particular positions has been well documented (e.g., Frija et al., 2011; Mutafova et al., 2016). In a previous study, we reported that three Ascomycota fungal strains (out of the eight tested) were able to transform LA (Frija et al., 2013). In particular, Penicillium janczewskii mediated a stereo-selective hydroxylation of LA yielding 3β-hydroxy-labdanolic acid as the major product (Figure 1). This terpenoid derivative would be, in general, difficult to attain by classical chemical methods and, so far, has never been identified in natural extracts (Frija et al., 2013). 3β-Hydroxy-labdanolic acid was used for the chemical synthesis of, e.g., 3β-hydroxy-labd-8(17)-en-15-oic acid—a moderate bactericidal compound, previously identified in the stems of Moldenhawera nutans (Frija et al., 2013).

FIGURE 1
www.frontiersin.org

Figure 1. Hydroxylation of labdanolic acid at the 3β position mediated by Penicillium janczewskii.

Most studies on fungal biotransformation of terpenoids have essentially identified/characterized the terpenoid derivatives being produced, whereas their general impact on the fungal metabolism and the biotransformation pathways remain largely unknown, hampering the biotechnological development of this field. This study aims to specifically analyze P. janczewskii metabolism during the biotransformation of LA. With this goal in mind, herein we performed a differential proteomic analysis of both mycelial and extracellular proteomes, complemented by in vivo inhibition assays of cytochrome P450 monooxygenase activity. LA apparently increased stress responses, especially against oxidative stress in P. janczewskii cultures and its stereo-selective hydroxylation likely involves a P450 enzyme.

Materials and Methods

Materials and Chemicals

The LA extract was obtained from twigs of Cistus ladanifer as previously described (i.e., crude acidic fraction of which ≈16% is LA) (Frija et al., 2013). Chemicals were purchased from Sigma Aldrich (USA), except electrophoresis reagents such as dithiothreitol (DTT), Triton X-100, 3-[(3-cholamidopropyl) dimethylammonio]-1-propanesulfonate (CHAPS) and IPG buffers that were from GE Healthcare (Sweden). All solvents used were of the highest analytical grade and water was obtained from a Milli-Q system (Millipore).

Fungal Strain and Cultivation Conditions

Penicillium janczewskii K. M. Zalessky (CBS#3498) conidia were harvested and maintained as frozen suspensions (Carvalho et al., 2009). Cultures were initiated with 105 conidia/mL and incubated without agitation in the dark at 27°C (nine replicates). The fungal cultures (20 mL) were grown in a mineral minimal medium (Carvalho et al., 2009) supplemented, per liter, with 10 g of glucose, 0.25 g of urea, and 2.5 g of the LA extract (hereafter defined solely as LA medium), which was added from a concentrated stock in ethanol prior to the media sterilization. Final concentration of LA in this medium was ca. 0.04 mg/mL. Controls without LA were prepared in similar conditions. To keep the experimental conditions similar to those used in our previous study (Frija et al., 2013), after 60 days of incubation both mycelia and filtrate were recovered, immediately frozen in liquid nitrogen and stored at −80°C until further analyses (Martins et al., 2014b). Small aliquots of the culture broths were used to evaluate qualitatively by thin layer chromatography (TLC, see below) the biotransformation of the plant LA. The presence of 3β-hydroxy-labdanolic acid was evaluated by mass spectrometry (see below). In addition, to evaluate catalase activity in the secretome of cultures grown in LA medium, the ability of extracellular proteins (concentrated using 10 kDa centrifugal devices and resuspended in 0.1 M phosphate buffer, pH 7) to decompose H2O2 in vitro (Aebi, 1984) was qualitatively verified by the release of oxygen bubbles.

In Vivo Inhibition of Cytochrome P450

Cultures were initiated from spores (105/mL) in the minimal media containing, per liter, 10 g of glucose, 0.25 g of urea, and 0.02 g of pure LA (Sigma) (LA standard medium) (5 mL, 27°C, dark, no agitation). Metyrapone (2.0 mM) was added to half of these cultures due to its capacity to inhibit in vivo cytochrome P450 monooxygenase activity (Strauber et al., 2003). Culture filtrates were collected at the 21st day of incubation, immediately extracted with diethyl ether and LA quantitatively analyzed by high performance liquid chromatography (HPLC) (see below). The identity of the compounds was validated by mass spectrometry (see below).

Thin Layer Chromatography

Culture filtrates were extracted with diethyl ether (1:1; six times), the ensuing extracts dried under soft nitrogen flow, and re-solubilized in diethyl ether (1 mL). The extracts (10 µL) were resolved in TLC silica gel 60 F254 plates (Merck) with diethyl ether/n-hexane (3:1) as the mobile phase. At the end of the TLC run, the plates were air dried, immersed in a solution of 10% of phosphomolybdic acid in ethanol, and revealed by heating at 80–100°C during 5 min.

High Performance Liquid Chromatography

Labdanolic acid was analyzed by HPLC, using an Alliance 2695 Waters chromatographer (Waters Corporation, Milford, MA, USA), connected to a Differential Refractometer (Bromma, Sweden) detector. Data acquisition was accomplished with the Empower 2 software (Waters). Chromatographic separation was undertaken using a Symmetry C18 column (4.6 mm × 250 mm), 5 µm particle size (Waters) set at 28°C. Elution was carried out isocratically with water/acetonitrile (10:90) at a flow rate of 1.0 mL/min, and the injection volume was 95 µL. The retention time of the LA was compared with that of a pure standard (Sigma) for identification, and the peak area was used for quantification. LA retention time was 5.7 min, and the established quantification limits were 0.125–1 mg/mL.

High Performance Liquid Chromatography–Electrospray Ionization–High Resolution Mass Spectrometry (HPLC-ESI-HRMS)

The diethyl ether extracts of the fungal cultures grown in the standard LA medium were analyzed by microLC–MS using a Triple TOF 6600 MS system (Sciex) equipped with the DuoSprayTM ion source. Chromatographic separation was carried out in the Eksigent ekspert nanoLC425 in microflow using an HALO C18 (50 mm × 0.5 mm, 2.7 µm particle size, 90 Å) column from Eksigent at the flow rate of 10 µL/min. The mobile phase consisted of a solution of 0.1% formic acid (solvent A) and a solution of acetonitrile containing 0.1% formic acid (solvent B), set as follows: 20% B in 2 min, followed by a linear gradient of 20–95% B in 12 min, 2 min of 95% B, 2 min to return to the initial conditions, and 5 min to re-equilibrate the column. MS was operated in positive ionization mode, with TOF MS scan with an m/z range 100–1,000 for 500 ms for a total cycle time of 0.5 s. MS data were processed using the PeakView software using the extracted-ion chromatogram for the compounds of interest. For compound identity a Δ(m/z) ≤ 5 ppm as well as a low noise/signal ratio, was considered.

Protein Extraction

Established methods were used to extract mycelial (Martins et al., 2013) and extracellular proteins (Martins et al., 2014b). Mycelial [100 mg powder mixed with 20 mg poly(vinyl)polypyrrolidone from Merck] and extracellular proteins (concentrated using 10 kDa centrifugal devices) were precipitated with cold acetone containing 10% w/v trichloroacetic acid and 60 mM DTT. The washed pellets were dissolved in similar buffers (7 M urea, 2 M thiourea, 4% w/v CHAPS, 60 mM DTT, and 1.0% v/v of matching IPG buffer) except that 1% (w/v) Triton X-100 was also added to the buffer used in the extracellular proteins, and finally clarified by centrifugation.

Two-dimensional Gel Electrophoresis (2DE)

Protein extracts were quantified using bovine serum albumin (BioRad) as a standard. Mycelial and extracellular protein samples (80 and 50 µg, respectively) were loaded in precast 13 cm non-linear IPG strips (GE Healthcare) (pH 3–10 NL and 3–5.6 NL, respectively). Protein isoelectric focusing, SDS-PAGE, staining, image acquisition, and analysis were done as previously reported (Martins et al., 2013, 2014b). The gels were stained with flamingo dye (BioRad, USA) for image acquisition but for spot excision colloidal coomassie blue (Fluka, Switzerland) was used instead and higher quantities of protein were loaded, namely, 180 and 100 µg for mycelial and extracellular protein fractions, respectively. For each condition, three biological replicates (each accounting three batch cultures) were done.

Data Acquisition and Image Analyses

Gels stained with the flamingo dye were scanned in a Fuji TLA-5100 scanner to generate 64 bit images. The images were analyzed in SameSpots v2.0 (non-linear dynamics) accordingly to manufacturer’s instructions and included gel alignment to a reference image. Spot detection was refined by manual edition whenever necessary, and protein spots with areas lower than 1,500 were excluded from the analysis. Differential protein spots in the mycelial proteome were identified using SameSpots software (ANOVA p-value < 0.05). Aiming at a stringent statistical analysis (Valledor et al., 2010), normalized spot volumes were used for calculating mean values, SD, and CV across the two secretome gel sets (fungal grown in the LA medium or the control medium) (Data Sheet S1 in Supplementary Material). Only those spots showing high consistency between replicate gels (CV below 15%) were used to generate the list of p-values for the individual spots (ANOVA, X-Stat). This allowed to pinpoint the spots that were differentially accumulated in the secretome (Data Sheet S1 in Supplementary Material).

Protein Identification

Differential protein spots on LA in both sub-proteomes were manually excised from gels and processed using the Ettan Digester robot of the Ettan Spot Handling Workstation (GE Healthcare) (Carvalho et al., 2013). Samples (0.7 µL) were then spotted on MALDI-TOF target plates (Applied Biosystems), before the deposit of 0.7 µL CHCA (7 mg/mL in ACN 50%, TFA 0.1%). Peptide mass determinations were carried out using the Applied Biosystems 5800 Proteomics Analyzer (Applied Biosystems) using established protocols (Carvalho et al., 2013; Martins et al., 2013, 2014b). Proteins were identified, with ProteinPilot, by searching against the NCBInr database (restricted to fungi taxonomy, 1,291,260 sequences) with Mascot v2.3 (Matrix Science). Homology identification was retained with probability set at 95%. All identifications were confirmed manually.

Protein Functional Classification

The identified proteins were classified into categories according to the Munich Information Centre for Protein Sequences Functional Catalog.1 For proteins with no assigned function, homology searches were performed using the BlastP program against all non-redundant protein sequences deposited in the NCBI database2 to attribute putative functions. The cellular location of the extracellular proteins was predicted using WoLF-PSORT (Horton et al., 2007).

Results

Similar to our previous study (Frija et al., 2013), herein we observed that P. janczewskii mediates a stereo-selective hydroxylation of LA yielding predominantly 3β-hydroxy-labdanolic acid (Figure 2) that displays an ion mass (m/z) of 363.2506 [LA]Na+ (data not shown). In these static cultures (20 mL), not all the LA was hydroxylated after the two months of incubation (Figure 2, lane III).

FIGURE 2
www.frontiersin.org

Figure 2. Thin layer chromatogram of Penicillium janczewskii cultures during the biotransformation of labdanolic acid (LA) to 3β-hydroxy-labdanolic acid (LA-OH): pure commercial LA (I); diethyl ether extract of P. janczewskii cultures in the LA medium at the incubation time (time 0) (II) and after 60 days of growth (III); and diethyl ether extract of the negative control (IV, i.e., P. janczewskii cultures after 60 days of cultivation in media without the LA extract). [Eluent: diethyl ether/n-hexane (3:1).]

In this study, P. janczewskii mycelial and extracellular sub-proteomes after growth in LA or control media were analyzed using a 2DE approach (Figure 3). In both sub-proteomes, the protein extraction yields were slightly higher in the LA medium compared to the control yet retrieving a comparable number of total protein spots, notwithstanding the long incubation time used (Table 1). Comparative proteomic analyses were undertaken to identify in each sub-proteome, the protein spots that increased during growth in the LA medium compared to the control medium (i.e., differential accumulation, ANOVA analysis, p < 0.05). The number of differentially accumulated proteins spots in the mycelial proteome and the secretome was respectively, 18 and 20, retrieving 17 and 7 unique protein species (Data Sheet S2 in Supplementary Material); herein depicted in Tables 2 and 3, where proteins are grouped by functional categories. Regardless that the genome of this fungus is still unsequenced, a large fraction of the differential protein spots in either sub-proteome matched protein species annotated in Penicillium species already sequenced (e.g., P. chrysosporium, P. digitatum, and P. marnafreii) (Data Sheet S2 in Supplementary Material). Most of the remaining protein identifications were matched against sequences from Aspergillus species. The proteomic analysis provided a snap-shot view of major alterations provoked by the supplementation of the growth media with LA compared to control conditions (i.e., absence of LA). In general, the LA extract increased stress responses in P. janczewskii, especially against oxidative stress [e.g., superoxide dismutase (SD) and NAD(P)H-quinone oxidoreductase, Table 2]. Catalases were also show to differentially accumulate (Table 3); consistent with the detection of catalase activity (i.e., hydrogen peroxide decomposition assay) in the extracellular protein extracts from cultures grown in LA medium. A detailed description of important cellular events likely associated with the identified differential protein species (Tables 2 and 3) is presented in the Section “Discussion.”

FIGURE 3
www.frontiersin.org

Figure 3. Two-dimensional gels obtained for the mycelia (A) and secretome (B) of Penicillium janczewskii cultures grown for 60 days in the presence of labdanolic acid. The spots highlighted showed a statistically significant increase (p-value < 0.05) during cultivation in media supplemented with the plant terpenoid (protein identifications are depicted in Tables 2 and 3).

TABLE 1
www.frontiersin.org

Table 1. Protein extraction yield and number of identified spots in the 2DE gels.

TABLE 2
www.frontiersin.org

Table 2. Penicillium janczewskii mycelial proteins identified in the 2DE differentially accumulated protein spots [fold change (FC)] in the labdanolic acid medium when compared to the control medium.

TABLE 3
www.frontiersin.org

Table 3. Penicillium janczewskii extracellular proteins identified in the 2DE differentially accumulated protein spots [fold change (FC)] in the labdanolic acid (LA) medium when compared to the control medium.

Among the set of identified protein species up-accumulated in the presence of LA, only a putative cytochrome P450 monooxygenase (P450s or CYPs) [spot 126, fold change (FC) = 1.17, Table 3] came forward as capable of mediating the stereo-selective hydroxylation of the terpenoid. To investigate the involvement of P450 enzymes in the biotransformation of LA, their activity was inhibited in vivo using pure LA (0.02 mg/mL) as a substitute of the plant extract (containing ≈0.04 mg/mL of LA). In the absence of P450 inhibitors, the bioconversion rate of pure LA reached ca. 40% at the 21st day of incubation, whereas in their presence the conversion was significantly inhibited, dropping to ca. 18% (Table 4, HPLC). To further verify that the stereo-selective hydroxylation of LA has occurred, the presence of 3β-hydroxy-labdanolic acid was verified by its ion mass using high resolution mass spectrometry. The MS results showed that the ion mass of LA could be detected in cultures grown in either the absence or presence of the P450 inhibitor (Figures 4B,C), whereas that of 3β-hydroxy-labdanolic acid could only be detected in cultures devoid of the P450 inhibitor (Figures 4D,E).

TABLE 4
www.frontiersin.org

Table 4. Biotransformation yields of pure labdanolic acid (LA) in the absence or presence of metyrapone (inhibitor of P450 activity) in Penicillium janczewskii cultures as determined by liquid chromatography.

FIGURE 4
www.frontiersin.org

Figure 4. Spectrometric analyses of the biotransformation of pure labdanolic acid (LA) in the absence or presence of metyrapone (inhibitor of P450 activity) in Penicillium janczewskii cultures. Total ion chromatograms (A), extracted-ion chromatogram (XIC) obtained for LA (B) and for 3β-hydroxy-labdanolic acid (D); and the isotopic pattern obtained for LA (C) and for 3β-hydroxy-labdanolic acid (E), illustrating the data collected during the HPLC-ESI-HRMS analyses.

Discussion

LA Stereo-Specific Hydroxylation by P. janczewskii

Terpenoids chemical modifications through biotransformation, in particular stereo-selective hydroxylation of non-activated carbons in their structures, provide derivatives which would be usually difficult to achieve by chemical methods (Frija et al., 2013; Schrader and Bohlmann, 2015; Kemper et al., 2017). Fungi can biotransform terpenoids in particular product(s) and at specific rate(s) through multiple reactions strongly influenced by the cultivation conditions (Frija et al., 2011; Ghasemi et al., 2014; Mutafova et al., 2016). However, production yields of each distinct product from a stereo-selective hydroxylation of the terpenoid are usually very low (2–20%). Herein P. janczewskii could efficiently hydroxylate LA; regardless that not all LA was converted to LA-OH (Figure 2, lane III). When pure LA (0.02 mg/mL) was used as a substitute of the plant extract (containing ≈0.04 mg/mL of LA), the bioconversion rate reached ca. 40% at the 21st day of incubation (Table 4; Figure 4).

In most studies reported so far, the terpenoid containing medium is inoculated with mycelia (previously grown in a standard growth medium) and agitated batch or fed-batch cultivation conditions are used (Frija et al., 2011; Schrader and Bohlmann, 2015; Mutafova et al., 2016). Terpenoids are known to play major roles in plant defense mechanisms (Gershenzon and Dudareva, 2007; Mutafova et al., 2016), probably explaining why the LA plant extract inhibited the germination of P. janczewskii conidia in the agitated cultures. In our previous study, even at static conditions, germination of conidia was only observed in three out of the eight fungal strains tested (Frija et al., 2011). Several studies support the catabolic potential of P. janczewskii, e.g., for producing secondary metabolites (Madi and Katan, 1998), exo-inulinases (Pessoni et al., 2007), xylanases (Terrasan et al., 2010, 2016, 2017), and galactosidases (Zhang et al., 2011), as well as its ecological relevance in modulating plant–fungi interactions, either antagonizing (Madi and Katan, 1998) or stimulating (Kwasna, 2004) important plant pathogenic fungi.

Comparative proteomic analyses (Figure 3) of P. janczewskii sub-proteomes allowed the identification of 24 unique protein species that differentially accumulated during growth on the LA medium (Tables 2 and 3; Data Sheet S2 in Supplementary Material). A small number of distinct proteins spots (Tables 2 and 3) retrieved identical functional annotations that might result from posttranslational modifications of the same gene product (i.e., proteolysis, glycosylation, and phosphorylation) (Mann and Jensen, 2003) or from the presence of sequence-related isoforms encoded by distinct paralogs (Vödisch et al., 2009). Since P. janczewskii fully annotated genome sequence is not yet available on the basis of MS data alone, one cannot discriminate between these two possibilities.

Stress Responses to LA Exposure in P. janczewskii Cultures

The biological effect of labdane-type terpenoids (e.g., manoyl oxides) is usually associated with their lipophilicity, and the presence of specific end groups that may interact with cellular membranes altering their permeability and ultimately leading to severe damage and lipid peroxidation (Matsingou and Demetzos, 2007). Terpenoids usually lead to accumulation of intracellular reactive oxygen species (ROS) (Liu et al., 2013) and increase oxidative stress response both at transcript (e.g., Grosmannia clavigera expose to pine extract) (Hesse-Orce et al., 2010) and protein levels (e.g., Saccharomyces cerevisiae incubation with d-limonene) (Liu et al., 2013). ROS can severely damage the cells, causing oxidative damage of membranes, DNA and proteins, hence their levels are tightly regulated in fungi through a set of antioxidant responses both enzymatic (e.g., SOD, glutathione peroxidase, and catalase) and non-enzymatic (e.g., ubiquitinol and glutathione) (Aguirre et al., 2006; Marschall and Tudzynski, 2016). Several proteins species associated with antioxidant response, namely, SOD and catalase showed increased levels in the LA medium compared to control conditions (Tables 2 and 3). In particular, SODs (spots 3905 and 3938, FC = 3.9 and 4.0, respectively) that catalyze the dismutation of the superoxide radical, producing oxygen and/or hydrogen peroxide, constitute the first line of defense against oxidative stress in fungi (Belozerskaya and Gessler, 2007). The produced hydrogen peroxide can be subsequently decomposed to water and molecular oxygen by catalases. Two catalases showed a minor increase in P. janczewskii secretome on LA (spots 41 and 73, Table 3). This finding is consistent with their frequent detection in fungal secretomes undergoing oxidative stress conditions (Adav et al., 2012; Martins et al., 2014b). Supplementation of the growth media with catalase has been shown to reduce the ROS level and alleviated cell growth inhibition during S. cerevisiae cultivation in media containing d-limonene (Liu et al., 2013).

Growth media supplementation with sclareol lead to increased oxygen consumption on Botrytis cinerea cultures, probably due to altered mitochondrial function (Mendoza et al., 2015). In addition, α-pinene was shown to affect strongly the energy metabolism of mitochondria isolated from maize (Abrahim et al., 2003; Mendoza et al., 2015). This terpenoid provoked the uncoupling of oxidative phosphorylation and the inhibition of electron transfer, likely due to unspecific damage in the inner mitochondrial membrane (Abrahim et al., 2003; Mendoza et al., 2015). Similar findings were reported for β-pinene which strongly inhibited respiration in yeast cells and in their isolated mitochondria, probably due to its capacity to increase membrane fluidity (Uribe et al., 1985). Herein, two mitochondrial proteins species increased in P. janczewskii mycelial proteome on LA, namely, outer mitochondrial membrane protein porin (spot 3851, FC = 4.2) and succinate dehydrogenase (spot 4041, FC = 4.8) (Table 2). The last enzyme catalyzes the oxidation of succinate to fumarate with the reduction of ubiquinone to ubiquinol, coupling the citric acid cycle and the electron transport chain in the inner mitochondrial membrane. Mitochondrial ubiquinol acts as antioxidant molecule, participating in non-enzymatic anti-ROS mechanisms (Bai et al., 2003; Dinkova-Kostova and Talalay, 2010). LA apparently altered fungal mitochondrial proteins and led to the generation of ROS, notwithstanding its effects were less severe than those induced by more lipophilic terpenoids such as α-pinene and limonene (Abrahim et al., 2000).

NAD(P)H-quinone oxidoreductase is a widely distributed FAD-dependent flavoprotein that promotes obligatory reductions of quinones (as well as of other aromatic molecules) thus dropping their intracellular levels and minimizing the generation of reactive oxygen intermediates by redox cycling (Dinkova-Kostova and Talalay, 2010). It plays major antioxidant roles fighting oxidative stress and is induced as a stress-responsive protein, e.g., in human cells exposed to Angelica sinensis extracts (Dietz et al., 2008) and in A. nidulans (AN0297) expose to either menadione (Pusztahelyi et al., 2011) or ionic liquids (Martins et al., 2013). LA also increased this stress response protein in the mycelia of P. janczewskii (spot 1304, FC = 1.9, Table 2).

A link between oxidative stress and increase of virulence has been suggested (Nikolaou et al., 2009), which might explain the differential accumulation of allergens in the mycelial proteome on LA, namely, the 16 kDa allergen and the allergen Asp f 15 (FC = 3.4 and FC = 4.6, respectively) (Table 2). Other stress-related proteins that increased on LA include the ubiquitous and highly conserved nucleoside diphosphate kinase (spot 4012, FC = 3.9, Table 2). This enzyme is crucial for the cellular homeostasis of nucleosides triphosphate and diphosphate (Lee et al., 2006). As an example, the transcription of its encoding gene is developmentally regulated and induced during stress response in Aspergillus spp. (Malavazi et al., 2006).

In this study, major accumulation of the unknown function protein containing a Bys1 domain (spot 4014, FC = 6.7, Table 2) was found in the mycelial proteome on LA. The antifungal caspofungin has been also shown to increase the levels of a protein containing a Bys 1 domain in A. fumigatus (Cagas et al., 2011). This class of proteins has been associated to morphogenesis, in particular the encoding gene increases during the yeast phase of the dimorphic fungus Blastomyces dermatitidis (Burg and Smith, 1994), notwithstanding the protein could not be identified in a proteomic deep analysis of dimorphism in P. marneffei (Chandler et al., 2008).

The differential protein spot number 3998 (FC = 5.3, Table 2) in the mycelial proteome retrieved two protein identifications, hampering the differential analysis. Both were identified as unknown function protein: one as hypothetical protein with no characterized function; and, the other with high homology to Grg1—a general stress protein involved in lifespan control, which was also identified in the spot 4023 (FC = 8.8). grg1 has been previously identified in dormant conidia of A. fumigatus (Sugui et al., 2008) and levels increased during starvation and asexual development in Neurospora crassa (ccg-1) (Xie et al., 2004). In this study, at the end of the cultivation conidia were visible in the mycelial mat at the liquid air interface (data not shown).

Additional spots found to significantly increase in the mycelial proteome on LA media (3840, 3839, and 4060 with FC = 2,118, 317.8, and 2.8, respectively, Table 2), were however identified as uncharacterized proteins. Future functional and kinetic studies on fungal proteins are essential to elucidate their roles during fungal growth in media supplemented with the plant terpenoid.

LA Hydroxylation

The differential accumulation of a putative P450 (spot 126, FC = 1.17) was found in the fungal secretome on LA (Table 3). This poorly characterized protein shows homology to a predicted P450 in Cordyceps militaris that is thought to contain only 290 amino acid residues. CYPs putative sequence lengths in A. nidulans (CYPome) were manually corrected from (286–750) to (417–607) amino acid residues (e.g., correction of incorrect exons and start/stop sequences) (Kelly et al., 2009). To the best of our knowledge, a similar analysis is missing in the CYPome of C. militaris, hence this predicted P450 possibly contains more than 290 residues. It shows also significant homology to CYP503B1 from Beauveria bassiana (Table 3); a class usually associated with the production of secondary metabolites in A. nidulans (Kelly et al., 2009; Moktali et al., 2012).

Genome sequencing projects on fungi are continuously revealing an increasing number of P450s (Park et al., 2008; Urlacher and Girhard, 2012; Chen et al., 2014), regardless the function of many remain unknown (Urlacher and Eiben, 2006; Chen et al., 2014). In fungi these enzymes play critical roles in the adaptation to specific ecological and/or nutritional niches by modifying potentially harmful chemicals, as well as in the production of an array of secondary metabolites. Interestingly, pine terpenoids increased the expression levels of several P450 genes of the secondary metabolism in G. clavigera, linking the encoded enzymes to their biotransformation (Lah et al., 2013).

P450s catalyze many reactions, such as oxygenation, dealkylation, epoxidation, reduction, dehalogenation and carbon hydroxylation (Urlacher and Girhard, 2012), including high stereo- and regio- selective hydroxylation of non-activated carbons, e.g., in Cochliobolus lunatus (Lah et al., 2011). They participate in the detoxification of a broad range of xenobiotics, particularly by mediating their initial modification (phase I), before transferase activity (conjugation with, e.g., glucoside, glucuronide, and sulfate—phase II) and cellular excretion (phase III) (Schäfer et al., 2004; Harms et al., 2011). None of the P. janczewskii differential accumulated proteins species on LA was identified as a transferase. However, conjugates of either LA or 3β-hydroxy-labdanolic acid, if any, might have been de-conjugated upon their cellular excretion through a reaction mediated by a catalase (Campoy et al., 2009). In fact, catalases were identified in multiple protein spots in the differential secretome (Table 3), and extracellular catalase activity was also detected.

The putative P450 identified herein (Table 3) is devoid of typical signal peptide sequence, notwithstanding its cellular compartmentalization is largely unknown (WoLF-PSORT predicted cellular locations: 12 plasma membrane bound; 6 extracellular; 3 endoplasmatic reticulum; 2 peroxisome; and 2 vacuolar). Though usually P450s are regarded as membrane-associated enzymes, they are increasingly being associated to the secretome (Druzhinina et al., 2012) as a consequence of autolysis or mycelial fragmentation (Kubicek, 2013) as well as secretion, e.g., secretome of Postia placenta growth on wood (5 days) (Vanden Wymelenberg et al., 2010). P450s found in fungal secretomes match multiple classes including some typically regarded as intracellular, emphasizing that current knowledge is still scarce for predicting the function of such uncharacterized enzymes. Several P450s have been linked to lignin degradation upon their detection in secretomes of fungi grown in lignocellulosic substrates (similar or twofold incubation time compared to that used here), e.g., Phanerochaete carnosa (Mahajan and Master, 2010) and Trametes trogii (Ji et al., 2012). In these studies, none of the P450 typical partners (i.e., NADPH-cytochrome P450 oxidoreductase) could be found in the fungal secretomes, similar to that observed here (Table 3).

The stereo-selective hydroxylation of pure LA was substantially inhibited when P450 activity was inhibited in vivo. LA hydroxylation levels were ca. 40 and 18% in the absence and presence of metyrapone (Table 4). The identity of 3β-hydroxy-labdanolic acid as well as LA were verified by mass spectrometry retrieving ion masses for [M + Na]+ of 363.2506 and 347.2570, respectively (full spectra and the isotopic pattern obtained for the extracts of the cultures is depicted in Figure 4). In the presence of metyrapone, the intensity of the ion chromatogram of 3β-hydroxy-labdanolic acid (Figure 4D) is comparable to the noise intensity; accordingly the corresponding isotopic pattern could not be extracted. This observation is consistent with the involvement of P450s in the biotransformation of LA.

Overall, the proteomic data provided novel information on the major cellular responses stimulated by the plant terpenoid during fungal growth, opening new perspectives of their role during plant–fungi interactions. Advanced processes for the biotransformation of terpenoids may require the use of antioxidants. We hope our study may promote the development of novel valorization strategies through biotechnology for the major terpenoids of the widespread resource Cistus ladanifer L.

Author Contributions

IM executed all the proteome analyses and prepared the first draft of the manuscript. IM and AV executed the fungal cultivation and prepared the extracts for analytics. LF and ME produced the labdanolic extracts and provided analytical support. JR and SP performed the MS analyses for protein identification. CSP and CA conceived the study and performed critical revision of the data. CSP elaborated the final manuscript. All the authors contributed to the analysis and interpretation of data and have read and approved the final version of the manuscript.

Conflict of Interest Statement

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Funding

IM and LF are grateful to Fundação para a Ciência e a Tecnologia (FCT), Portugal, for the fellowships SFRH/BPD/110841/2015 and SFRH/BPD/99851/2014. The work was partially supported by FCT: PTDC/QUI/73061/2006, PTDC/QUI-QUI/119823/2010, and Pest-OE/SAU/UI4013/2011. The authors acknowledge funding from FCT through grant UID/Multi/04551/2013 (Research unit GREEN-it “Bioresources for Sustainability”). The authors wish to acknowledge the technical support provided by Mª Cristina Leitão for executing the HPLC analyses and Dr. Ricardo Gomes for the HRMS data obtained (UniMS—Mass Spectrometry Unit, ITQB NOVA/iBET, Oeiras, Portugal).

Supplementary Material

The Supplementary Material for this article can be found online at http://journal.frontiersin.org/article/10.3389/fbioe.2017.00045/full#supplementary-material.

Footnotes

References

Abrahim, D., Braguini, W. L., Kelmer-Bracht, A. M., and Ishii-Iwamoto, E. L. (2000). Effects of four monoterpenes on germination, primary root growth, and mitochondrial respiration of maize. J. Chem. Ecol. 26, 611–624. doi: 10.1023/A:1005467903297

CrossRef Full Text | Google Scholar

Abrahim, D., Francischini, A. C., Pergo, E. M., Kelmer-Bracht, A. M., and Ishii-Iwamoto, E. L. (2003). Effects of α-pinene on the mitochondrial respiration of maize seedlings. Plant Physiol. Biochem. 41, 985–991. doi:10.1016/j.plaphy.2003.07.003

CrossRef Full Text | Google Scholar

Adav, S. S., Chao, L. T., and Sze, S. K. (2012). Quantitative secretomic analysis of Trichoderma reesei strains reveals enzymatic composition for lignocellulosic biomass degradation. Mol. Cell. Proteomics 11, M111.012419. doi:10.1074/mcp.M111.012419

PubMed Abstract | CrossRef Full Text | Google Scholar

Aebi, H. (1984). Catalase in vitro. Method Enzymol. 105, 121–126. doi:10.1016/S0076-6879(84)05016-3

CrossRef Full Text | Google Scholar

Aguirre, J., Hansberg, W., and Navarro, R. (2006). Fungal responses to reactive oxygen species. Med. Mycol. 44, S101–S107. doi:10.1080/13693780600900080

CrossRef Full Text | Google Scholar

Bai, Z. H., Harvey, L. M., and McNeil, B. (2003). Physiological responses of chemostat cultures of Aspergillus niger (B1-D) to simulated and actual oxidative stress. Biotechnol. Bioeng. 82, 691–701. doi:10.1002/bit.10614

PubMed Abstract | CrossRef Full Text | Google Scholar

Belozerskaya, T. A., and Gessler, N. N. (2007). Reactive oxygen species and the strategy of antioxidant defense in fungi: a review. Appl. Biochem. Microbiol. 43, 506–515. doi:10.1134/S0003683807050031

CrossRef Full Text | Google Scholar

Bicas, J. L., Dionisio, A. P., and Pastore, G. M. (2009). Bio-oxidation of terpenes: an approach for the flavor industry. Chem. Rev. 109, 4518–4531. doi:10.1021/cr800190y

CrossRef Full Text | Google Scholar

Bolster, M. G., Jansen, B. J. M., and de Groot, A. (2001). The synthesis of (-)-Ambrox® starting from labdanolic acid. Tetrahedron 57, 5657–5662. doi:10.1016/S0040-4020(01)00493-8

CrossRef Full Text | Google Scholar

Burg, E. F., and Smith, L. H. (1994). Cloning and characterization of bys1, a temperature-dependent cDNA specific to the yeast phase of the pathogenic dimorphic fungus Blastomyces dermatitidis. Infect. Immun. 62, 2521–2528.

PubMed Abstract | Google Scholar

Cagas, S. E., Jain, M. R., Li, H., and Perlin, D. S. (2011). Profiling the Aspergillus fumigatus proteome in response to caspofungin. Antimicrob. Agents Chemother. 55, 146–154. doi:10.1128/AAC.00884-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Campoy, S., Alvarez-Rodríguez, M. L., Recio, E., Rumbero, A., and Coque, J. J. (2009). Biodegradation of 2,4,6-TCA by the white-rot fungus Phlebia radiata is initiated by a phase I (O-demethylation)-phase II (O-conjugation) reactions system: implications for the chlorine cycle. Environ. Microbiol. 11, 99–110. doi:10.1111/j.1462-2920.2008.01744.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, M. B., Martins, I., Leitão, M. C., Garcia, H., Rodrigues, C., San Romão, V., et al. (2009). Screening pentachlorophenol degradation ability by environmental fungal strains belonging to the phyla Ascomycota and Zygomycota. J. Ind. Microbiol. Biotechnol. 36, 1249–1256. doi:10.1007/s10295-009-0603-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Carvalho, M. B., Martins, I., Medeiros, J., Tavares, S., Planchon, S., Renaut, J., et al. (2013). The response of Mucor plumbeus to pentachlorophenol: a toxicoproteomics study. J. Proteomics 78, 159–171. doi:10.1016/j.jprot.2012.11.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Chandler, J. M., Treece, E. R., Trenary, H. R., Brenneman, J. L., Flickner, T. J., Frommelt, J. L., et al. (2008). Protein profiling of the dimorphic pathogenic fungus Penicillium marneffei. Proteome Sci. 6, 17. doi:10.1186/1477-5956-6-17

PubMed Abstract | CrossRef Full Text | Google Scholar

Chen, W., Lee, M. K., Jefcoate, C., Kim, S. C., Chen, F., and Yu, J. H. (2014). Fungal cytochrome P450 monooxygenases: their distribution, structure, functions, family expansion, and evolutionary origin. Genome Biol. Evol. 6, 1620–1634. doi:10.1093/gbe/evu132

PubMed Abstract | CrossRef Full Text | Google Scholar

Dietz, B. M., Liu, D., Hagos, G. K., Yao, P., Schinkovitz, A., Pro, S. M., et al. (2008). Angelica sinensis and its alkylphthalides induce the detoxification enzyme NAD(P)H: quinone oxidoreductase 1 by alkylating keap1. Chem. Res. Toxicol. 21, 1939–1948. doi:10.1021/tx8001274

PubMed Abstract | CrossRef Full Text | Google Scholar

Dinkova-Kostova, A. T., and Talalay, P. (2010). NAD(P)H:quinone acceptor oxidoreductase 1 (NQO1), a multifunctional antioxidant enzyme and exceptionally versatile cytoprotector. Arch. Biochem. Biophys. 501, 116–123. doi:10.1016/j.abb.2010.03.019

PubMed Abstract | CrossRef Full Text | Google Scholar

Druzhinina, I. S., Shelest, E., and Kubicek, C. P. (2012). Novel traits of Trichoderma predicted through the analysis of its secretome. FEMS Microbiol. Lett. 337, 1–9. doi:10.1111/j.1574-6968.2012.02665.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Frija, L. M. T., Frade, R. F. M., and Afonso, C. A. M. (2011). Isolation, chemical and biotransformation routes of labdane-type diterpenes. Chem. Rev. 111, 4418–4452. doi:10.1021/cr100258k

CrossRef Full Text | Google Scholar

Frija, L. M. T., Garcia, H., Rodrigues, C., Martins, I., Candeias, N. R., André, V., et al. (2013). Short synthesis of the natural product 3 β-hydroxy-labd-8(17)-en-15-oic acid via microbial transformation of labdanolic acid. Phytochem. Lett. 6, 165–169. doi:10.1016/j.phytol.2012.12.005

CrossRef Full Text | Google Scholar

Gershenzon, J., and Dudareva, N. (2007). The function of terpene natural products in the natural world. Nat. Chem. Biol. 3, 408–414. doi:10.1038/nchembio.2007.5

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghasemi, S., Mohajeri, M., and Habibi, Z. (2014). Biotransformation of testosterone and testosterone heptanoate by four filamentous fungi. Steroids 92, 7–12. doi:10.1016/j.steroids.2014.09.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Ghosh, S., and Rangan, L. (2014). Unveiling the effect of plant derived labdane diterpenes as promising bactericidal and anticancer agents. Planta Med. 80, 783–783. doi:10.1055/s-0034-1382451

CrossRef Full Text | Google Scholar

Harms, H., Schlosser, D., and Wick, L. Y. (2011). Untapped potential: exploiting fungi in bioremediation of hazardous chemicals. Nat. Rev. Microbiol. 9, 177–192. doi:10.1038/nrmicro2519

PubMed Abstract | CrossRef Full Text | Google Scholar

Hawas, U., Eltomy, S., Nassif, R., El-Hossary, G., and AbouZid, S. (2015). Retraction. A new antifungal labdane diterpene from the leaves of Saraca indica. Nat. Prod. Res. 29, i–vi. doi:10.1080/14786419.2013.855931

CrossRef Full Text | Google Scholar

Hesse-Orce, U., DiGuistini, S., Keeling, C. I., Wang, Y., Li, M., Henderson, H., et al. (2010). Gene discovery for the bark beetle-vectored fungal tree pathogen Grosmannia clavigera. BMC Genomics 11:536. doi:10.1186/1471-2164-11-536

PubMed Abstract | CrossRef Full Text | Google Scholar

Horton, P., Park, K. J., Obayashi, T., Fujita, N., Harada, H., Adams-Collier, C. J., et al. (2007). WoLF PSORT: protein localization predictor. Nucleic Acids Res. 35, 585–587. doi:10.1093/nar/gkm259

PubMed Abstract | CrossRef Full Text | Google Scholar

Huang, Z., Zhu, Z. X., Li, Y. T., Pang, D. R., Zheng, J., Zhang, Q., et al. (2015). Anti-inflammatory labdane diterpenoids from Leonurus macranthus. J. Nat. Prod. 78, 2276–2285. doi:10.1021/acs.jnatprod.5b00635

PubMed Abstract | CrossRef Full Text | Google Scholar

Jayaprakasam, B., Alexander-Lindo, R. L., DeWitt, D. L., and Nair, M. G. (2007). Terpenoids from stinking toe (Hymneae courbaril) fruits with cyclooxygenase and lipid peroxidation inhibitory activities. Food Chem. 105, 485–490. doi:10.1016/j.foodchem.2007.04.004

CrossRef Full Text | Google Scholar

Ji, X.-L., Zhang, W.-T., Gai, Y.-P., Lua, B. Y., Yuana, C. Z., Liu, Q. X., et al. (2012). Patterns of lignocellulose degradation and secretome analysis of Trametes trogii MT. Int. Biodeterior. Biodegradation 75, 55–62. doi:10.1016/j.ibiod.2012.09.001

CrossRef Full Text | Google Scholar

Kapewangolo, P., Omolo, J. J., Bruwer, R., Fonteh, P., and Meyer, D. (2015). Antioxidant and anti-inflammatory activity of Ocimum labiatum extract and isolated labdane diterpenoid. J. Inflamm. (Lond) 12, 4. doi:10.1186/s12950-015-0049-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Kelly, D. E., Krasevec, N., Mullins, J., and Nelson, D. R. (2009). The CYPome (cytochrome P450 complement) of Aspergillus nidulans. Fungal Genet. Biol. 46, S53–S61. doi:10.1016/j.fgb.2008.08.010

PubMed Abstract | CrossRef Full Text | Google Scholar

Kemper, K., Hirte, M., Reinbold, M., Fuchs, M., and Brück, T. (2017). Opportunities and challenges for the sustainable production of structurally complex diterpenoids in recombinant microbial systems. Beilstein J. Org. Chem. 13, 845–854. doi:10.3762/bjoc.13.85

PubMed Abstract | CrossRef Full Text | Google Scholar

Kubicek, C. P. (2013). Systems biological approaches towards understanding cellulase production by Trichoderma reesei. J. Biotechnol. 163, 133–142. doi:10.1016/j.jbiotec.2012.05.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Kwasna, H. (2004). Natural shifts in communities of rhizosphere fungi of common oak after felling. Plant Soil 264, 209–218. doi:10.1023/B:PLSO.0000047752.41575.c7

CrossRef Full Text | Google Scholar

Lah, L., Haridas, S., Bohlmann, J., and Breuil, C. (2013). The cytochromes P450 of Grosmannia clavigera: genome organization, phylogeny, and expression in response to pine host chemicals. Fungal Genet. Biol. 50, 72–81. doi:10.1016/j.fgb.2012.10.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Lah, L., Podobnik, B., Novak, M., Korošec, B., Berne, S., Vogelsang, M., et al. (2011). The versatility of the fungal cytochrome P450 monooxygenase system is instrumental in xenobiotic detoxification. Mol. Microbiol. 81, 1374–1389. doi:10.1111/j.1365-2958.2011.07772.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Lee, B., Yoshida, Y., and Hasunuma, K. (2006). Photomorphogenetic characteristics are severely affected in nucleoside diphosphate kinase-1 (NDK-1)-disrupted mutants in Neurospora crassa. Mol. Genet. Genomics 275, 9–17. doi:10.1007/s00438-005-0044-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Liu, J., Zhu, Y., Du, G., Zhou, J., and Chen, J. (2013). Response of Saccharomyces cerevisiae to d-limonene-induced oxidative stress. Appl. Microbiol. Biotechnol. 97, 6467–6475. doi:10.1007/s00253-013-4931-9

PubMed Abstract | CrossRef Full Text | Google Scholar

Madi, L., and Katan, J. (1998). Penicillium janczewskii and its metabolites, applied to leaves, elicit systemic acquired resistance to stem rot caused by Rhizoctonia solani. Physiol. Mol. Plant Pathol. 53, 163–175. doi:10.1006/pmpp.1998.0174

CrossRef Full Text | Google Scholar

Mahajan, S., and Master, E. R. (2010). Proteomic characterization of lignocellulose-degrading enzymes secreted by Phanerochaete carnosa grown on spruce and microcrystalline cellulose. Appl. Microbiol. Biotechnol. 86, 1903–1914. doi:10.1007/s00253-010-2516-4

PubMed Abstract | CrossRef Full Text | Google Scholar

Malavazi, I., Savoldi, M., Di Mauro, S. M., Menck, C. F., Harris, S. D., Goldman, M. H., et al. (2006). Transcriptome analysis of Aspergillus nidulans exposed to camptothecin-induced DNA damage. Eukaryotic Cell 5, 1688–1704. doi:10.1128/EC.00167-06

PubMed Abstract | CrossRef Full Text | Google Scholar

Mann, M., and Jensen, O. N. (2003). Proteomic analysis of post-translational modifications. Nat. Biotechnol. 21, 255–261. doi:10.1038/nbt0303-255

PubMed Abstract | CrossRef Full Text | Google Scholar

Marschall, R., and Tudzynski, P. (2016). Reactive oxygen species in development and infection processes. Semin. Cell Dev. Biol. 57, 138–146. doi:10.1016/j.semcdb.2016.03.020

PubMed Abstract | CrossRef Full Text | Google Scholar

Martins, A. N. C., Simeonov, S. P., Frija, L. M. T., Viveirosb, R., Lourençob, A., da Silva, M. S., et al. (2014a). Isolation, analytical quantification and seasonal variation of labdanolic acid from the Portuguese-grown Cistus ladaniferus. Indl. Crops Prod. 60, 226–232. doi:10.1016/j.indcrop.2014.06.012

CrossRef Full Text | Google Scholar

Martins, I., Garcia, H., Varela, A., Núñez, O., Planchon, S., Galceran, M. T., et al. (2014b). Investigating Aspergillus nidulans secretome during colonisation of cork cell walls. J. Proteomics 98, 175–188. doi:10.1016/j.jprot.2013.11.023

PubMed Abstract | CrossRef Full Text | Google Scholar

Martins, I., Hartmann, D. O., Alves, P. C., Planchon, S., Renaut, J., Leitão, M. C., et al. (2013). Proteomic alterations induced by ionic liquids in Aspergillus nidulans and Neurospora crassa. J. Proteomics 94, 262–278. doi:10.1016/j.jprot.2013.09.015

PubMed Abstract | CrossRef Full Text | Google Scholar

Matsingou, C., and Demetzos, C. (2007). Effect of the nature of the 3 β-substitution in manoyl oxides on the thermotropic behavior of DPPC lipid bilayer and on DPPC liposomes. J. Liposome Res. 17, 89–105. doi:10.1080/08982100701375076

CrossRef Full Text | Google Scholar

Mendoza, L., Sepulveda, C., Melo, R., and Cotoras, M. (2015). Characterization of the antifungal activity against botrytis cinerea of sclareol and 13-epi-sclareol, two labdane-type diterpenoids. J. Chil. Chem. Soc. 60, 3024–3028. doi:10.4067/S0717-97072015000300010

CrossRef Full Text | Google Scholar

Moktali, V., Park, J., Fedorova-Abrams, N. D., Park, B., Choi, J., Lee, Y. H., et al. (2012). Systematic and searchable classification of cytochrome P450 proteins encoded by fungal and oomycete genomes. BMC Genomics 13:525. doi:10.1186/1471-2164-13-525

PubMed Abstract | CrossRef Full Text | Google Scholar

Mutafova, B., Fernandes, P., Mutafov, S., Berkov, S., and Pavlov, A. (2016). “Microbial transformations of plant secondary metabolites,” in Bioprocessing of Plant In Vitro Systems, eds A. Pavlov and T. Bley (Cham: Springer International Publishing), 1–41.

Google Scholar

Nikolaou, E., Agrafioti, I., Stumpf, M., Quinn, J., Stansfield, I., and Brown, A. J. (2009). Phylogenetic diversity of stress signalling pathways in fungi. BMC Evol. Biol. 9:44. doi:10.1186/1471-2148-9-44

PubMed Abstract | CrossRef Full Text | Google Scholar

Papaefthimiou, D., Papanikolaou, A., Falara, V., Givanoudi, S., Kostas, S., and Kanellis, A. K. (2014). Genus Cistus: a model for exploring labdane-type diterpenes’ biosynthesis and a natural source of high value products with biological, aromatic, and pharmacological properties. Front. Chem. 2:35. doi:10.3389/fchem.2014.00035

PubMed Abstract | CrossRef Full Text | Google Scholar

Park, J., Lee, S., Choi, J., Ahn, K., Park, B., Park, J., et al. (2008). Fungal cytochrome P450 database. BMC Genomics 9:402. doi:10.1186/1471-2164-9-402

CrossRef Full Text | Google Scholar

Pessoni, R. A. B., Braga, M. R., and Cassia, R. L. (2007). Purification and properties of exo-inulinases from Penicillium janczewskii growing on distinct carbon sources. Mycologia 99, 493–503. doi:10.1080/15572536.2007.11832543

PubMed Abstract | CrossRef Full Text | Google Scholar

Pusztahelyi, T., Klement, E., Szajli, E., Klem, J., Miskei, M., Karányi, Z., et al. (2011). Comparison of transcriptional and translational changes caused by long-term menadione exposure in Aspergillus nidulans. Fungal Genet. Biol. 48, 92–103. doi:10.1016/j.fgb.2010.08.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Schäfer, S., Schrader, J., and Sell, D. (2004). Oxygen uptake rate measurements to monitor the activity of terpene transforming fungi. Proc. Biochem. 39, 2221–2228. doi:10.1016/j.procbio.2003.11.023

CrossRef Full Text | Google Scholar

Schrader, J., and Bohlmann, J. (2015). Biotechnology of Isoprenoids. Hannover, Germany: Springer International Publishing.

Google Scholar

Silva, E. O., Furtado, N. A. J. C., Aleu, J., and Collado, I. G. (2013). Terpenoid biotransformations by Mucor species. Phytochem. Rev. 12, 857–876. doi:10.1007/s11101-013-9313-5

CrossRef Full Text | Google Scholar

Singh, S., Kumar, J. K., Saikia, D., Shanker, K., Thakur, J. P., Negi, A. S., et al. (2010). A bioactive labdane diterpenoid from Curcuma amada and its semisynthetic analogues as antitubercular agents. Eur. J. Med. Chem. 45, 4379–4382. doi:10.1016/j.ejmech.2010.06.006

PubMed Abstract | CrossRef Full Text | Google Scholar

Strauber, H., Muller, R. H., and Babel, W. (2003). Evidence of cytochrome P450-catalyzed cleavage of the ether bond of phenoxybutyrate herbicides in Rhodococcus erythropolis K2-3. Biodegradation 14, 41–50. doi:10.1023/A:1023550209155

PubMed Abstract | CrossRef Full Text | Google Scholar

Sugui, J. A., Kim, H. S., Zarember, K. A., Chang, Y. C., Gallin, J. I., Nierman, W. C., et al. (2008). Genes differentially expressed in conidia and hyphae of Aspergillus fumigatus upon exposure to human neutrophils. PLoS ONE 3:e2655. doi:10.1371/journal.pone.0002655

CrossRef Full Text | Google Scholar

Terrasan, C. R., Temer, B., Duarte, M. C., and Carmona, E. C. (2010). Production of xylanolytic enzymes by Penicillium janczewskii. Bioresour. Technol. 101, 4139–4143. doi:10.1016/j.biortech.2010.01.011

PubMed Abstract | CrossRef Full Text | Google Scholar

Terrasan, C. R. F., Guisan, J. M., and Carmona, E. C. (2016). Xylanase and β-xylosidase from Penicillium janczewskii: purification, characterization and hydrolysis of substrates. Electron. J. Biotechnol. 23, 54–62. doi:10.1016/j.ejbt.2016.08.001

CrossRef Full Text | Google Scholar

Terrasan, C. R. F., Romero-Fernández, M., Orrego, A. H., Oliveira, S. M., Pessela, B. C., Carmona, E. C., et al. (2017). Immobilization and stabilization of beta-xylosidases from Penicillium janczewskii. Appl. Biochem. Biotechnol. 182, 349–366. doi:10.1007/s12010-016-2331-1

PubMed Abstract | CrossRef Full Text | Google Scholar

Tholl, D. (2015). “Biosynthesis and biological functions of terpenoids in plants,” in Biotechnology of Isoprenoids, eds J. Schrader and J. Bohlmann (Cham: Springer), 63–106.

Google Scholar

Uribe, S., Ramirez, J., and Pena, A. (1985). Effects of β-pinene on yeast membrane functions. J. Bacteriol. 161, 1195–1200.

Google Scholar

Urlacher, V. B., and Eiben, S. (2006). Cytochrome P450 monooxygenases: perspectives for synthetic application. Trends Biotechnol. 24, 324–330. doi:10.1016/j.tibtech.2006.05.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Urlacher, V. B., and Girhard, M. (2012). Cytochrome P450 monooxygenases: an update on perspectives for synthetic application. Trends Biotechnol. 30, 26–36. doi:10.1016/j.tibtech.2011.06.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Valledor, L., Jorrín, J. V., Rodríguez, J. L., Lenz, C., Meijón, M., Rodríguez, R., et al. (2010). Combined proteomic and transcriptomic analysis identifies differentially expressed pathways associated to Pinus radiata needle maturation. J. Proteome Res. 9, 3954–3979. doi:10.1021/pr1001669

PubMed Abstract | CrossRef Full Text | Google Scholar

Vanden Wymelenberg, A., Gaskell, J., Mozuch, M., Sabat, G., Ralph, J., Skyba, O., et al. (2010). Comparative transcriptome and secretome analysis of wood decay fungi Postia placenta and Phanerochaete chrysosporium. Appl. Environ. Microbiol. 76, 3599–3610. doi:10.1128/AEM.00058-10

PubMed Abstract | CrossRef Full Text | Google Scholar

Vödisch, M., Albrecht, D., Lessing, F., Schmidt, A. D., Winkler, R., Guthke, R., et al. (2009). Two-dimensional proteome reference maps for the human pathogenic filamentous fungus Aspergillus fumigatus. Proteomics 9, 1407–1415. doi:10.1002/pmic.200800394

PubMed Abstract | CrossRef Full Text | Google Scholar

Xie, X., Wilkinson, H. H., Correa, A., Lewis, Z. A., Bell-Pedersen, D., and Ebbole, D. J. (2004). Transcriptional response to glucose starvation and functional analysis of a glucose transporter of Neurospora crassa. Fungal Genet. Biol. 41, 1104–1119. doi:10.1016/j.fgb.2004.08.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Yang, N. Y., Liu, L., Tao, W. W., Duan, J. A., and Tian, L. J. (2010). Diterpenoids from Pinus massoniana resin and their cytotoxicity against A431 and A549 cells. Phytochemistry 71, 1528–1533. doi:10.1016/j.phytochem.2010.06.008

PubMed Abstract | CrossRef Full Text | Google Scholar

Zhang, B., Chen, Y. Q., Li, Z. M., Lu, W., and Cao, Y. (2011). Cloning and functional expression of α-galactosidase cDNA from Penicillium janczewskii Zaleski. Biologia 66, 205–212. doi:10.2478/s11756-011-0014-5

CrossRef Full Text | Google Scholar

Keywords: labdanolic acid, stress response, terpenoids biotransformation, proteomics, Penicillium janczewskii, cytochrome P450 monooxygenase

Citation: Martins I, Varela A, Frija LMT, Estevão MAS, Planchon S, Renaut J, Afonso CAM and Silva Pereira C (2017) Proteomic Insights on the Metabolism of Penicillium janczewskii during the Biotransformation of the Plant Terpenoid Labdanolic Acid. Front. Bioeng. Biotechnol. 5:45. doi: 10.3389/fbioe.2017.00045

Received: 13 October 2016; Accepted: 10 July 2017;
Published: 31 July 2017

Edited by:

Didier Laurent Buisson, Centre National de la Recherche Scientifique (CNRS), France

Reviewed by:

Noha M. Mesbah, Suez Canal University, Egypt
Gustavo Molina, Universidade Federal dos Vales do Jequitinhonha e Mucuri (UFVJM), Brazil

Copyright: © 2017 Martins, Varela, Frija, Estevão, Planchon, Renaut, Afonso and Silva Pereira. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) or licensor are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Cristina Silva Pereira, c3BlcmVpcmEmI3gwMDA0MDtpdHFiLnVubC5wdA==

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.