Next Article in Journal
Feasibility of Casein to Record Stable Isotopic Variation of Cow Milk in New Zealand
Next Article in Special Issue
Effects of Ultrasonic-Assisted Extraction on the Yield and the Antioxidative Potential of Bergenia emeiensis Triterpenes
Previous Article in Journal
d-allulose Ameliorates Metabolic Dysfunction in C57BL/KsJ-db/db Mice
Previous Article in Special Issue
Phenolic Composition, Antioxidant and Anticancer Potentials of Extracts from Rosa banksiae Ait. Flowers
 
 
Search for Articles:
Title / Keyword
Author / Affiliation / Email
Journal
Article Type
 
 
Advanced Search
 
Section
Special Issue
Volume
Issue
Number
Page
 
Journals
Molecules
Volume 25
Issue 16
10.3390/molecules25163657
molecules-logo
Submit to this Journal Review for this Journal Propose a Special Issue
► Article Menu

Article Menu

share Share announcement Help format_quote Cite question_answer Discuss in SciProfiles thumb_up
...
Endorse textsms
...
Comment
first_page
settings
Order Article Reprints
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Deletion of the Histone Deacetylase HdaA in Endophytic Fungus Penicillium chrysogenum Fes1701 Induces the Complex Response of Multiple Bioactive Secondary Metabolite Production and Relevant Gene Cluster Expression

by
Zhuang Ding
1,*,
Haibo Zhou
2,
Xiao Wang
1,
Huiming Huang
3,
Haotian Wang
4,
Ruiyan Zhang
1,
Zhengping Wang
1 and
Jun Han
1
1
Institute of BioPharmaceutical Research, Liaocheng University, Liaocheng 252059, China
2
Shandong University-Helmholtz Institute of Biotechnology, State Key Laboratory of Microbial Technology, Shandong University, Qingdao 266237, China
3
School of Life Science, Liaocheng University, Liaocheng 252059, China
4
Faculty of Pharmacy, Bengbu Medical College, Bengbu 233000, China
*
Author to whom correspondence should be addressed.
Molecules 2020, 25(16), 3657; https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25163657
Submission received: 20 July 2020 / Revised: 8 August 2020 / Accepted: 9 August 2020 / Published: 11 August 2020
(This article belongs to the Special Issue Advances on Natural Compounds Research: Chemistry, Biology and Biotechnology)
Browse Figures
Versions Notes

Abstract

:
Epigenetic regulation plays a critical role in controlling fungal secondary metabolism. Here, we report the pleiotropic effects of the epigenetic regulator HdaA (histone deacetylase) on secondary metabolite production and the associated biosynthetic gene clusters (BGCs) expression in the plant endophytic fungus Penicillium chrysogenum Fes1701. Deletion of the hdaA gene in strain Fes1701 induced a significant change of the secondary metabolite profile with the emergence of the bioactive indole alkaloid meleagrin. Simultaneously, more meleagrin/roquefortine-related compounds and less chrysogine were synthesized in the ΔhdaA strain. Transcriptional analysis of relevant gene clusters in ΔhdaA and wild strains indicated that disruption of hdaA had different effects on the expression levels of two BGCs: the meleagrin/roquefortine BGC was upregulated, while the chrysogine BGC was downregulated. Interestingly, transcriptional analysis demonstrated that different functional genes in the same BGC had different responses to the disruption of hdaA. Thereinto, the roqO gene, which encodes a key catalyzing enzyme in meleagrin biosynthesis, showed the highest upregulation in the ΔhdaA strain (84.8-fold). To our knowledge, this is the first report of the upregulation of HdaA inactivation on meleagrin/roquefortine alkaloid production in the endophytic fungus P. chrysogenum. Our results suggest that genetic manipulation based on the epigenetic regulator HdaA is an important strategy for regulating the productions of secondary metabolites and expanding bioactive natural product resources in endophytic fungi.

Graphical Abstract

1. Introduction

Filamentous fungi are well-known producers of diverse secondary metabolites (SMs), which have a wide range of biological activities and can be beneficial or harmful to human beings [1,2]. On the one hand, many beneficial fungi-derived SMs have long been clinically utilized as antibacterials (penicillin and cephalosporin), antifungals (anidulafungin and caspofungin), immunosuppressants (cyclosporin), and antihypercholesterolemic drugs (lovastatin) [3]. On the other hand, some mycotoxins can severely endanger public health, such as aflatoxins, fumonisins, fusarins, and gliotoxin [4]. As an industrial microbial strain used for producing the β-lactam antibiotic penicillin, Penicillium chrysogenum has the capacity to synthesize many SMs with diverse chemical structures and significant bioactivities, such as alkaloids, polyketides, and terpenoids [5,6,7]. In recent decades, genome sequencing and genetic studies of P. chrysogenum have led to the elucidation of many secondary metabolic pathways for SM production [8,9]. Subsequently, the functions of some transcription regulators associated with these biosynthetic pathways have been well studied. For instance, the regulatory factors PcRFX1 and PcFKH1 have been characterized to positively control penicillin biosynthesis in Penicillium chrysogenum [10,11]. Similarly, the deletion of the transcription factor gene laeA resulted in a drastic decrease of penicillin gene expression [12]. In contrast, another transcription factor CreA was reported to negatively regulate penicillin production [13].
Recent studies have shown that the secondary metabolism of filamentous fungi is controlled by a complicated and elaborate regulatory network, which is influenced by not only various transcription factors but also epigenetic regulators [14]. Among various types of epigenetic regulators in eukaryotic cells, histone deacetylases (HDACs) play an important role and profoundly influence DNA replication, transcription, and repair processes [15]. Furthermore, several reports have determined that histone deacetylation induced by HDACs tends to be associated with heterochromatin and gene silencing [16]. Thus, the inactivation of HDACs has been widely considered as an available strategy for the activation of silent SM biosynthetic pathways and increasing the production of bioactive natural products in filamentous fungi. Shwab et al. found that the deletion of the HDAC-encoding gene hdaA resulted in the transcriptional activation of penicillin and sterigmatocystin biosynthetic gene cluster (BGC) in Aspergillus nidulans strain A89 [17]. Deletion of the homologue of hdaA in Aspergillus fumigatus strain AF293 increased the transcription of several NRPS (nonribosomal peptide synthetase) gene clusters [18]. By disrupting hdaA in Calcarisporium arbuscular NRRL3705, Mao et al. activated 75% of the SM biosynthetic genes and found four new compounds [19]. Moreover, the deletion of hdaA in Pestalotiopsis fici CGMCC3.15140 activated the production of a series of macrodiolides [20]. In addition, it was found that HDACs could also regulate the transcription of the genes that were associated with conidiation, sexual reproduction, growth, stress response, and pathogenicity [21].
As a characteristic bioresource, endophytic fungi provide a broad variety of important SMs with diverse bioactivities [22,23]. These SMs are not only relevant to the physiological and ecological peculiarity of fungi, but are also interrelated with human life and health. Recent investigations on endophyte genomes revealed diverse secondary metabolic BGCs [19,24]. However, the considerable biosynthetic potential of endophytic fungi is a reflection of the complex ecological environment, which is difficult to simulate under laboratory conditions, resulting in most BGCs being silent [19]. According to the aforementioned reports, genetic manipulation based on HDAC genes represents a feasible strategy for regulating and activating secondary metabolic BGCs in eukaryotes. However, the application of this regulation technique is limited to endophytic fungi. In this study, we constructed an hdaA homologue deletion mutant of P. chrysogenum Fes1701, an endophytic fungus isolated from rubber tree (Ficus elastica) leaves [25], and investigated the effects of HdaA inactivation on secondary metabolism. Our results indicate that the ΔhdaA strain showed a significant change of secondary metabolic profile with the emergence of some SMs. In addition, we describe the effects of HdaA inactivation on the expression of two SM-associated BGCs of P. chrysogenum Fes1701.

2. Results and Discussion

2.1. Identification and Deletion of the hdaA Gene in P. chrysogenum Fes1701

In the whole genomic sequence of P. chrysogenum Wis54-1255, the hdaA gene is located in chromosome 2 and is designated as Pc21g14570 [8]. The hdaA gene of the strain Fes1701 was located and identified in the genomic sequence via Local-BLAST. Subsequently, we designed special primers to clone the hdaA ORF from the genomic DNA of the strain Fes1701. The obtained PCR fragment was 3039 bp in size, and the predicted coding sequence was 2304 bp, which encoded a 767 amino acid polypeptide. To confirm its identity, bioinformatics analysis of HdaA was performed using an NCBI BLAST search. BLAST analysis indicated that this protein showed 66.8% sequence identity to HdaA (AN8042) of Aspergillus nidulans [26]. The taxonomic relatedness of HdaA and other known homologous protein sequences from other species is shown by a phylogenetic tree in Figure 1. In addition, reverse transcription PCR was performed to determine the transcription level of hdaA. Results indicated that the hdaA transcript was normally expressed when the wild type (WT) of the Fes1701 strain was grown at 25 °C on PDA (potato dextrose agar) medium (Figure S1).
Upstream and downstream flanking fragments of the hdaA ORF were PCR-amplified for the construction of a gene disruption cassette (Table S1). An hdaA deletion mutant (ΔhdaA) was constructed by replacing this gene with a bleomycin resistance gene (bleoR) cassette (Figure 2A). A 0.1% bleomycin-supplemented medium was used for the positive selection of the transformant. The genomic DNA of the selected transformant was extracted and further verified using diagnostic PCR (Figure S2). The 1.9-kb and 2.1-kb fragments could be amplified from the correct ΔhdaA strain using the primers VP1–VP5 and VP2–VP6, respectively, but were absent in the WT. No PCR product was amplified from the ΔhdaA strain using primers VP3–VP4, while a 1.7-kb product could appear using the genomic DNA of the WT (Figure S2). The WT and ΔhdaA strains were cultured on PDA medium at 25 °C for phenotype observation. Comparative results indicate that there was no significant difference in phenotype between the WT and ΔhdaA strains (Figure 2B).

2.2. Effects of hdaA Disruption on SMs Production

To examine the effect of hdaA deletion on secondary metabolism in the strain Fes1701, the SMs extract of the ΔhdaA strain was analyzed using HPLC and compared with that of the WT. After cultivating in PDB (potato dextrose broth) medium for 5 days, the metabolite fingerprints showed obvious differences between the ΔhdaA and WT strains (Figure 3A,B). Four main SMs (1−4) in the HPLC fingerprint of the ΔhdaA strain were further purified for the identification of their chemical structures (Figure 3C). By comparison with the published MS and NMR data in the literature [27,28,29,30], these compounds were identified as chrysogine (1), meleagrin (2), roquefortine F (3), and roquefortine C (4) (Figures S3–S10). The productions of four compounds in the ΔhdaA and WT strains were comparatively analyzed by calculating the HPLC peak area. Results showed that chrysogine (1) in the ΔhdaA strain had a 3.8-fold decrease as compared with the WT. Simultaneously, the production of compounds 2−4 showed 14.4-, 2.1-, and 1.4-fold increases in the ΔhdaA strain (Figure 3B). Similar multiple effects of HDAC on the production of SMs have been found in other fungal species [18,31]. For instance, Lee et al. reported that deletion of the hdaA gene in A. fumigatus increased the production of several SMs but decreased the production of gliotoxin [18].

2.3. Effects of hdaA Disruption on the Transcription of SM Biosynthetic Gene Clusters

Because the deletion of ΔhdaA greatly influenced the production of chrysogine (1) and meleagrin (2), we further investigated the differences in the transcriptional levels of genes from two BGCs encoding the enzymes responsible for the production of these two compounds between the ΔhdaA and WT strains. Chrysogine (1) is a yellow alkaloid produced by several filamentous fungi, such as the genera Penicillium, Aspergillus, and Fusarium [32]. Although this compound was first isolated in 1973, its BGC was not clarified until 2017 [32,33]. Viggiano et al. [33] elucidated the chrysogine biosynthetic pathway in P. chrysogenum and identified the function of each enzyme involved in this pathway. The P. chrysogenum chrysogine biosynthetic pathway contains a NRPS (Pc21g12630) flanked by five associated genes (Pc21g12570, Pc21g12590, Pc21g12600, Pc21g12610, and Pc21g12620), designated as chyA to chyH (Figure 4A). In our study, the transcriptional levels of these genes contained in the chrysogine BGC in ΔhdaA were significantly downregulated as compared with the WT (Figure 4B). Thereinto, the transcriptional levels of Pc21g12630 (chyA), which encodes an NRPS, and Pc21g12620 (chyD), which encodes an amidase, showed a 4.0-fold and 5.0-fold decrease in ΔhdaA, respectively. A similar phenomenon was observed by Guzman-Chavez et al. [34], who also reported that the inactivation of HdaA caused the downregulation of chrysogine BCG expression in the industrial strain P. chrysogenum DS68530.
In P. chrysogenum, the meleagrin/roquefortine BGC is a seven-gene cluster extending from Pc21g15420 to Pc21g15480 (Figure 5A). It encodes proteins from the biosynthetic pathway where, as was confirmed, roquefortine C is a biosynthetic precursor of meleagrin [35]. The biosynthesis of roquefortine depends on the involvement of the upstream genes, including the Pc21g15420 (roqT), Pc21g15430 (roqD), Pc21g15440 (roqN), and Pc21g15480 (roqA) genes. Subsequently, the Pc21g15450 (roqO) and Pc21g15460 (roqM) genes, which respectively encode a P450 scaffold reorganizing oxygenase and a MAK1 monooxygenase, are responsible for the biotransformation from roquefortine to meleagrin [35]. In our study, a comparative analysis of the transcriptional levels showed that these upstream genes responsible for roquefortine were 2.5- to 7.5-fold upregulated in ΔhdaA (Figure 5B). Furthermore, most significant enhancements appeared downstream of this pathway, which were essential for the biosynthesis of meleagrin. In this stage, Pc21g15450 (roqO) and Pc21g15460 (roqM) showed an 84.8- and 43.7-fold increase in ΔhdaA (Figure 5B). This phenomenon was relevant to the production increase of the meleagrin/roquefortine BGC in the ΔhdaA strain. The productions of roquefortine and related compounds exhibited a significant increase, while meleagrin appeared in the metabolite profiling of the ΔhdaA strain as a newly generated compound (Figure 3A). In the previous report, the production of meleagrin/roquefortine-related compounds in P. chrysogenum were unaffected by the inactivation of HdaA [34]. However, the production of these compounds in the ΔhdaA strain were significantly increased in this study. This phenomenon could indicate that the secondary metabolic pathway in the same fungal species derived from different habitats may be affected synergistically by other additional regulatory mechanisms. For example, a heterotrimeric Ga protein Pga1 has been reported to upregulate the biosynthesis of roquefortine in P. chrysogenum Wis54-1255 [36]. Our results demonstrate that it is important to reveal the different effects of epigenetic regulators in the same organism on secondary metabolism.

2.4. Bioactivities of the Metabolites Isolated from the ΔhdaA Strain

The meleagrin and roquefortine groups have been reported to possess various biological properties, such as antibacterial, neurotoxic, cytochrome P450 inhibitory, and tubulin polymerization inhibitory activities [37]. Because of the increasing production of meleagrin and roquefortine compounds (2–4) in the ΔhdaA strain, their antimicrobial and cytotoxic activities were further evaluated in this study (Table 1). Meleagrin (2) exhibited the best growth-inhibitory activity against two different leukemia cells (K562 and HL-60) with IC50 values of 8.9 and 12.7 μM. A similar level of antineoplastic activity was reported by Mady et al. [38], in which meleagrin (2) showed significant inhibitory activities against the proliferation and migration of c-Met-dependent breast malignancy. Roquefortine alkaloids (3,4) showed a medium antibacterial activity against Gram-negative bacterium Escherichia coli and Gram-positive bacterium Staphylococcus aureus. The inhibitory activity of roquefortine alkaloids against Gram-positive bacteria have been reported [39], especially the antibacterial activity of roquefortine C (4) against Bacillus subtilis with an MIC of 6.3 µg/mL [40]. However, all tested compounds did not show obvious antifungal activity in two pathogenic yeast fungi (MIC > 128 μg/mL).

3. Materials and Methods

3.1. Strains and Media

The Penicillium chrysogenum Fes1701 wild type (WT) was originally isolated from fresh leaves of Ficus elastica collected in Liaocheng University Arboretum, Liaocheng, China [24]. The WT and its transformants were grown on PDA medium at 25 °C, 5 days for sporulation, and the fresh spores were harvested with 0.01% Tween 80. The Escherichia coli strain DH5α was cultured in Luria–Bertani liquid medium at 37 °C with 100 µg/mL ampicillin for plasmid propagation.

3.2. Cloning and Identification of the hdaA Gene in Strain Fes1701

Fungal genomic DNA was isolated using the CTAB method [25]. The oligonucleotide sequences of PCR primers are listed in Table S1. The HdaA protein sequence of P. chrysogenum Wisconsin 54-1255 (Pc21g14570) was used as query to perform BLAST searches to find its ortholog in the strain Fes1701. The genome sequence of P. chrysogenum Wisconsin 54-1255 was obtained from the European Molecular Biology Laboratory (EMBL) (Heidelberg, Germany) under accession numbers AM920416–AM920464. The designed PCR primers PchdaA-F and PchdaA-R were used to amplify the hdaA ORF from the genomic DNA of strain Fes1701 using the TransStart® FastPfu DNA Polymerase (TransGen Biotech, Beijing, China). A 3.0-kb fragment was obtained and subcloned using a pEasy® Blunt Cloning Kit (TransGen Biotech, Beijing, China) to yield pEasy-PcHdaA and sequenced by Sangon Biotech (Shanghai, China). The sequence data of the hdaA gene in strain Fes1701 were deposited in the GenBank Database under accession number MT313928. The coding sequence for amino acid polypeptide was predicted by the FGenesh software (Softberry Inc., NY, USA) [41]. Phylogenetic analysis of HdaA orthologs from strain Fes1701 and other fungi was performed using the MEGA 7.0 software [42].

3.3. Creation of the ΔhdaA Strain

Deletion of the hdaA gene in the WT was performed through homologous recombination. A gene disruption cassette was assembled using the fusion PCR method as described previously [43] and inserted into pEasy-Blunt Cloning Vector (TransGen Biotech) to gain the plasmid pEB-HDD. The gene disruption cassette was comprised of an upstream flanking region (2.2 kb, amplified by the primer ΔHdaA-P1 and P2), bleomycin resistance gene (bleoR, 1.5 kb, amplified from the plasmid pZeo [43] by the primer ΔHdaA-P3 and P4), and downstream flanking region (1.9 kb, amplified by the primer ΔHdaA-P5 and P6), consecutively. The disruption cassette was excised from pEB-HDD with the restriction enzymes KpnI and PacI and then transferred into the WT, yielding the ΔhdaA strain by PEG-mediated protoplast transformation. Transformants were selected and single-spore purified in the presence of bleomycin. The genotype of the deletion mutant was confirmed by diagnostic PCR.

3.4. Fermentation Conditions and SM Extraction

For SMs production, the ΔhdaA and WT strains were cultured in PDB medium. The cultures were inoculated with 1 × 105 fresh spores in 250 mL flasks containing 100 mL of PDB medium and grown under shaking condition at 180 rpm, 25 °C for 5 days. Each test contained three replicates for each strain.
The extraction of fermentation products was performed following the protocols established by Ding et al. [25]. After cultivation, the mycelia in the fermentation mixture were broken using a disperser (T18, IKA, Germany), then the fermentation mixture was extracted with an equal volume of ethyl acetate. The extract was evaporated under reduced pressure and redissolved in 2 mL MeOH. Sterile PDB medium as the control was extracted using the same procedure.

3.5. Metabolite Fingerprint Analysis

A 10 µL extract was injected into the HPLC system (Waters Inc., Milford, MA, USA), which contained a model 1525 pump, an ODS column (Pack ODS-A, 250 × 4.6 mm, 5 µm, YMC Co., Ltd., Japan), and a model 2489 UV detector. The gradient increased from 20% to 100% MeOH over 30 min and then was retained at 100% MeOH for 10 min. The fold differences of these compounds between the mutant and WT were calculated by HPLC peak area according to the following formula: [Area (Sample)—Area (Control)]/[Area (WT)–Area (Control)].

3.6. Purification and Identification of Natural Products

A two-liter-scale culture and extract preparation of the ΔhdaA strain were performed using the method mentioned above. The obtained extract (1.3 g) was gained and separated by silica gel vacuum liquid chromatography using MeOH-H2O to give four fractions (Fractions 1–4). Fraction 1 was further separated by Sephadex LH-20 eluted with MeOH and then on a semipreparative HPLC column (Pack ODS-A, 250 × 10 mm, 5 µm, YMC Co., Ltd., Japan) eluted with MeOH-H2O (50:50, 3 mL/min) to provide compound 1 (2.6 mg, tR 8.5 min). Fraction 2 was separated by semipreparative HPLC eluted with MeOH-H2O (60:40, 3 mL/min) to obtain compound 2 (5.7 mg, tR 10.5 min). Fraction 4 was separated on a Sephadex LH-20 column with CH2Cl2-MeOH (1:1) and then on a semipreparative HPLC eluted with MeOH-H2O (75:25, 3 mL/min) to obtain compound 3 (3.5 mg, tR 11.5 min) and compound 4 (3.6 mg, tR 12.5 min). The structures of the compounds were identified using MS and NMR data. MS spectra were recorded on a Q-TOF Ultima Global GAA076 LC mass spectrometer (Wasters Inc., Milford, MA, USA). NMR spectra were collected on a Varian 500 spectrometer (Varian Medical Systems Inc., Palo Alto, CA, USA).

3.7. RNA Extraction and Real-Time PCR Analysis

The relative expression levels of the chrysogine and roquefortine/meleagrin biosynthetic gene clusters in the ΔhdaA and WT strains were analyzed by real-time PCR. The mutant and WT were cultivated in PDB medium as the abovementioned fermentation condition. The mycelia of each strain were collected on the 5th day, and RNA was extracted from the mycelia using the EasyPure® RNA kit for qRT-PCR (ER101, TransGen Biotech, China) following the manufacturer’s protocol. Then, the quality of the RNA was checked by a NanoDrop 2000 nucleotide analyser (Thermo Scientific, Waltham, MA, USA). cDNA was synthesized using the TransScript® kit (AH341, TransGen Biotech, China). Real-time PCR was performed using a CFX96 Real-Time System (Bio-Rad, Hercules, CA, USA) using the TransStart® kit (AQ132, TransGen Biotech, China). The PCR reaction system were 10 μL 2 × qPCR mix, 0.5 μL forward/reverse primer (10 μM), 1 μL template cDNA, and water to 20 μL. The PCR reaction conditions were 95 °C for 2 min, 45 cycles of 95 °C for 5 s, 60 °C for 15 s, 72 °C for 25 s, and 95 °C for 10 s. Three replicates of each cDNA sample were carried out, and the average threshold cycle was calculated. Relative expression levels were calculated using the 2−ΔΔCt method with the expression level of the actin gene as a reference. The primers used for real-time PCR are listed in Table S1.

3.8. Bioactivity Assay

Cytotoxicity against two human tumor cell lines, K562 (human chronic myeloblastic leukemia cells) and HL-60 (human promyelocytic leukemia cells), was evaluated using the methyl-thiazolyl-tetrazolium (MTT) method with Adriamycin as positive control [44].
Antimicrobial activities against four microorganisms, including the Gram-negative bacterium Escherichia coli CMCC 44102, the Gram-positive bacterium Staphylococcus aureus CMCC 26003, the fungus Candida albicans ATCC 10231, and C. glabrata ATCC 15126, were performed using the sequential 2-fold dilutions method as previously reported [45]. Chloramphenicol and fluconazole were used as the positive controls for antibacterial and antifungal assays, respectively.

3.9. Statistical Analysis

Comparisons of results were analyzed using the GraphPad software (version 7.01) (GraphPad software, San Diego, CA, USA) followed by the Tukey–Kramer test at p < 0.01. Values with asterisks are statistically significant.

4. Conclusions

As a characteristic bioresource, endophytic microorganisms have been increasingly recognized as a significant reservoir for exploiting bioactive SMs [46]. With the availability of a huge number of fungal genome sequences, SM BGCs and various metabolic regulators associated with these BGCs have become one of the research hot spots in fungal chemical biology. However, the identification of epigenetic regulators and their effects on secondary metabolism is still limited, especially in endophytic fungi. In the present study, we constructed an hdaA deletion mutant of the endophytic fungus P. chrysogenum Fes1701 and verified its effects on secondary metabolism. SMs profile analysis showed that the ΔhdaA strain produced less chrysogine- and more roquefortine-related compounds. Transcriptional analysis demonstrated that the deletion of hdaA activated the expression of the meleagrin/roquefortine BGC but inhibited the expression of the chrysogine BGC. Results from the current study suggest that the genetic manipulation of the epigenetic regulator HdaA represents a promising approach for activating and promoting yields of bioactive compounds in endophytic fungi.

Supplementary Materials

Table S1: Primers used in this study; Figure S1: Transcriptional analysis of a potential hdaA regulator gene by reverse transcription PCR; Figure S2: PCR verification of the ΔhdaA strain; Figures S3–S6: 1H-NMR spectrum of compounds 1–4 (CH3Cl-d3); Figures S7–S10: MS spectrum of compounds 1–4.

Author Contributions

Z.D. conceived and designed the experiments, prepared the manuscript, and was involved in the isolation of compounds. H.Z. contributed to the determination of the structures of isolated compounds. X.W. and H.H. performed the genetic manipulation and bioactivity assay. H.W. contributed to the preparation of fermentation and extraction. R.Z., Z.W., and J.H. supervised the work and revised the manuscript. All authors have read and approved the final version of the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was financially supported by the Natural Science Foundation of Shandong Province (ZR2017BB077, ZR2018BH043), the Bengbu Medical College Natural Science Fund Project (BYKY1768), and the Tai-Shan Scholar Foundation of Shandong Province (319190201).

Acknowledgments

The authors are grateful to Yi Tang (University of California, Los Angeles) for providing plasmids, and Dehai Li (Ocean University of China, Qingdao) for his technical assistance.

Conflicts of Interest

The authors declare that they have no conflicts of interest.

References

  1. Newman, D.J.; Cragg, G.M. Natural products as sources of new drugs from 1981 to 2014. J. Nat. Prod. 2016, 79, 629–661. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Strobel, G.; Daisy, B.; Castillo, U.; Harper, J. Natural products from endophytic microorganisms. J. Nat. Prod. 2004, 67, 257–268. [Google Scholar] [CrossRef] [PubMed]
  3. Hoffmeister, D.; Keller, N.P. Natural products of filamentous fungi: Enzymes, genes, and their regulation. Nat. Prod. Rep. 2007, 24, 393–416. [Google Scholar] [CrossRef] [PubMed]
  4. Wiemann, P.; Keller, N.P. Strategies for mining fungal natural products. J. Ind. Microbiol. Biotechnol. 2014, 41, 301–313. [Google Scholar] [CrossRef] [PubMed]
  5. Xu, W.F.; Mao, N.; Xue, X.J.; Qi, Y.X.; Wei, M.Y.; Wang, C.Y.; Shao, C.L. Structures and absolute configurations of diketopiperazine alkaloids chrysopiperazines A–C from the gorgonian-derived Penicillium chrysogenum Fungus. Mar. Drugs 2019, 17, 250. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Zhao, D.L.; Yuan, X.L.; Du, Y.M.; Zhang, Z.F.; Zhang, P. Benzophenone derivatives from an algal-endophytic isolate of Penicillium chrysogenum and their cytotoxicity. Molecules 2018, 23, 3378. [Google Scholar] [CrossRef] [Green Version]
  7. Qi, B.; Liu, X.; Mo, T.; Zhu, Z.; Li, J.; Wang, J.; Shi, X.; Zeng, K.; Wang, X.; Tu, P.; et al. 3,5-dimethylorsellinic acid derived meroterpenoids from Penicillium chrysogenum MT-12, an endophytic fungus isolated from Huperzia serrata. J. Nat. Prod. 2017, 80, 2699–2707. [Google Scholar] [CrossRef]
  8. Van den Berg, M.A.; Albang, R.; Albermann, K.; Badger, J.H.; Daran, J.M.; Driessen, A.J.; Garcia-Estrada, C.; Fedorova, N.D.; Harris, D.M.; Heijne, W.H.; et al. Genome sequencing and analysis of the filamentous fungus Penicillium chrysogenum. Nat. Biotechnol. 2008, 26, 1161–1168. [Google Scholar] [CrossRef] [Green Version]
  9. Kozlovskii, A.G.; Antipova, T.V.; Zhelifonova, V.P. Biosynthesis of biologically active low-molecular weight compound by fungi of the genus Penicillium (Review). Appl. Biochem. Microbiol. 2015, 51, 236–242. [Google Scholar] [CrossRef]
  10. Domínguez-Santos, R.; Martín, J.F.; Kosalková, K.; Prieto, C.; Ullán, R.V.; García-Estrada, C. The regulatory factor PcRFX1 controls the expression of the three genes of b-lactam biosynthesis in Penicillium chrysogenum. Fungal Genet. Biol. 2012, 49, 866–881. [Google Scholar] [CrossRef]
  11. Domínguez-Santos, R.; García-Estrada, C.; Kosalková, K.; Prieto, C.; Santamarta, I.; Martín, J.F. PcFKH1, a novel regulatory factor from the forkhead family, controls the biosynthesis of penicillin in Penicillium chrysogenum. Biochimie 2015, 115, 162–176. [Google Scholar] [CrossRef] [PubMed]
  12. Kosalková, K.; García-Estrada, C.; Ullán, R.V.; Godio, R.P.; Feltrer, R.; Teijeira, F.; Mauriz, E.; Martín, J.F. The global regulator LaeA controls penicillin biosynthesis, pigmentation and sporulation, but not roquefortine C synthesis in Penicillium chrysogenum. Biochimie 2009, 91, 214–225. [Google Scholar] [CrossRef]
  13. Cepeda-García, C.; Domínguez-Santos, R.; García-Rico, R.O.; García-Estrada, C.; Cajiao, A.; Fierro, F.; Martín, J.F. Direct involvement of the CreA transcription factor in penicillin biosynthesis and expression of the pcbAB gene in Penicillium chrysogenum. Appl. Microbiol. Biotechnol. 2014, 98, 7113–7124. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  14. Cichewicz, R.H. Epigenome manipulation as a pathway to new natural product scaffolds and their congeners. Nat. Prod. Rep. 2010, 27, 11–22. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  15. Robyr, D.; Suka, Y.; Xenarios, I.; Kurdistani, S.K.; Wang, A.; Suka, N.; Grunstein, M. Microarray deacetylation maps determine genome-wide functions for yeast histone deacetylases. Cell 2002, 109, 437–446. [Google Scholar] [CrossRef] [Green Version]
  16. Bulger, M. Hyperacetylated chromatin domains: Lessons from heterochromatin. J. Biol. Chem. 2005, 280, 21689–21692. [Google Scholar] [CrossRef] [Green Version]
  17. Shwab, E.K.; Bok, J.W.; Tribus, M.; Galehr, J.; Graessle, S.; Keller, N.P. Histone deacetylase activity regulates chemical diversity in Aspergillus. Eukaryot. Cell 2007, 6, 1656–1664. [Google Scholar] [CrossRef] [Green Version]
  18. Lee, I.; Oh, J.H.; Shwab, E.K.; Dagenais, T.R.; Andes, D.; Keller, N.P. HdaA, a class 2 histone deacetylase of Aspergillus fumigatus, affects germination and secondary metabolite production. Fungal Genet. Biol. 2009, 46, 782–790. [Google Scholar] [CrossRef] [Green Version]
  19. Mao, X.M.; Xu, W.; Li, D.; Yin, W.B.; Chooi, Y.H.; Li, Y.Q.; Tang, Y.; Hu, Y. Epigenetic genome mining of an endophytic fungus leads to the pleiotropic biosynthesis of natural products. Angew. Chem. Int. Ed. 2015, 54, 7592–7596. [Google Scholar] [CrossRef] [Green Version]
  20. Wu, G.W.; Zhou, H.C.; Zheng, P.; Wang, X.; Li, W.; Zhang, W.; Liu, X.; Liu, H.W.; Keller, N.P.; An, Z.; et al. Polyketide production of pestaloficiols and macrodiolide ficiolides revealed by manipulation of epigenetic regulators in an endophytic fungus. Org. Lett. 2016, 18, 1832–1835. [Google Scholar] [CrossRef]
  21. Niu, X.; Hao, X.; Hong, Z.; Chen, L.; Yu, X.; Zhu, X. A putative histone deacetylase modulates the biosynthesis of pestalotiollide B and conidiation in Pestalotiopsis microspora. J. Microbiol. Biotechnol. 2015, 25, 579–588. [Google Scholar] [CrossRef] [PubMed]
  22. Stierle, A.; Strobel, G.A.; Stierle, D.B. Taxol and taxane production by Taxomyces andreanae, an endophytic fungus of Pacific yew. Science 1993, 260, 214–216. [Google Scholar] [CrossRef]
  23. Strobel, G.A.; Miller, R.V.; Miller, C.; Condron, M.; Teplow, D.B.; Hess, W.M. Cryptocandin, a potent antimycotic from the endophytic fungus Cryptosporiopsis cf. quercina. Microbiology 1999, 145, 1919–1926. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  24. Xu, W.; Wang, F.; Zhang, M.; Ou, T.; Wang, R.; Strobel, G.; Xiang, Z.; Zhou, Z.; Xie, J. Diversity of cultivable endophytic bacteria in mulberry and their potential for antimicrobial and plant growth-promoting activities. Microbiol. Res. 2019, 229, 126328. [Google Scholar] [CrossRef] [PubMed]
  25. Ding, Z.; Tao, T.; Wang, L.; Zhao, Y.; Huang, H.; Zhang, D.; Liu, M.; Wang, Z.; Han, J. Bioprospecting of novel and bioactive metabolites from endophytic fungi isolated from rubber tree Ficus elastica leaves. J. Microbiol. Biotechnol. 2019, 29, 731–738. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Tribus, M.; Galehr, J.; Trojer, P.; Brosch, G.; Loidl, P.; Marx, F.; Haas, H.; Graessle, S. HdaA, a major class 2 histone deacetylase of Aspergillus nidulans, affects growth under conditions of oxidative stress. Eukaryot Cell 2005, 4, 1736–1745. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Niederer, D.; Tamm, C.; Zurcher, W. Nitrogen containing metabolites of Fusarium sambucinum. Tetrahedron Lett. 1992, 33, 3997–4000. [Google Scholar] [CrossRef]
  28. Kawai, K.; Nozawa, K.; Nakajima, S.; Iitaka, Y. Studies on fungal products. VII. The structures of meleagrin and 9-O-p-bromobenzoylmeleagrin. Chem. Pharm. Bull. 1984, 32, 94–98. [Google Scholar] [CrossRef] [Green Version]
  29. Musuku, A.; Selala, M.I.; de Bruyne, T.; Claeys, M.; Schepens, P.J.C.; Tsatsakis, A.; Shtilman, M.I. Isolation and structure determination of a new roquefortine-related mycotoxin from penicillium verrucosum var. cyclopium isolated from cassava. J. Nat. Prod. 1994, 57, 983–987. [Google Scholar] [CrossRef]
  30. Du, L.; Li, D.; Zhu, T.; Cai, S.; Wang, F.; Xiao, X.; Gu, Q. New alkaloids and diterpenes from a deep ocean sediment derived fungus Penicillium sp. Tetrahedron 2009, 65, 1033–1039. [Google Scholar] [CrossRef]
  31. Studt, L.; Schmidt, F.J.; Jahn, L.; Sieber, C.M.; Connolly, L.R.; Niehaus, E.M.; Freitag, M.; Humpf, H.U.; Tudzynski, B. Two histone deacetylases, FfHda1 and FfHda2, are important for Fusarium fujikuroi secondary metabolism and virulence. Appl. Environ. Microb. 2013, 79, 7719–7734. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Wollenberg, R.D.; Saei, W.; Westphal, K.R.; Klitgaard, C.S.; Nielsen, K.L.; Lysøe, E.; Gardiner, D.M.; Wimmer, R.; Sondergaard, T.E.; Sørensen, J.L. Chrysogine biosynthesis is mediated by a two-module nonribosomal peptide synthetase. J. Nat. Prod. 2017, 80, 2131–2135. [Google Scholar] [CrossRef] [PubMed]
  33. Viggiano, A.; Salo, O.; Ali, H.; Szymanski, W.; Lankhorst, P.P.; Nygård, Y.; Bovenberg, R.A.L.; Driessen, A.J.M. Pathway for the biosynthesis of the pigment chrysogine by penicillium chrysogenum. Appl. Environ. Microbiol. 2018, 84, e02246-17. [Google Scholar] [CrossRef] [Green Version]
  34. Guzman-Chavez, F.; Salo, O.; Samol, M.; Ries, M.; Kuipers, J.; Bovenberg, R.A.L.; Vreeken, R.J.; Driessen, A.J.M. Deregulation of secondary metabolism in a histone deacetylase mutant of Penicillium chrysogenum. MicrobiologyOpen 2018, 7, e00598. [Google Scholar] [CrossRef] [Green Version]
  35. García-Estrada, C.; Ullán, R.V.; Albillos, S.M.; Fernández-Bodega, M.Á.; Durek, P.; von Döhren, H.; Martín, J.F. A single cluster of coregulated genes encodes the biosynthesis of the mycotoxins roquefortine C and meleagrin in Penicillium chrysogenum. Chem. Biol. 2011, 18, 1499–1512. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  36. García-Rico, R.O.; Fierro, F.; Mauriz, E.; Gómez, A.; Fernández-Bodega, M.A.; Martín, J.F. The heterotrimeric Galpha protein pga1 regulates biosynthesis of penicillin, chrysogenin and roquefortine in Penicillium chrysogenum. Microbiology 2008, 154, 3567–3578. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  37. Du, L.; Feng, T.; Zhao, B.; Li, D.; Cai, S.; Zhu, T.; Wang, F.; Xiao, X.; Gu, Q. Alkaloids from a deep ocean sediment-derived fungus Penicillium sp. and their antitumor activities. J. Antibiot. 2010, 63, 165–170. [Google Scholar] [CrossRef] [Green Version]
  38. Mady, M.S.; Mohyeldin, M.M.; Ebrahim, H.Y.; Elsayed, H.E.; Houssen, W.E.; Haggag, E.G.; Soliman, R.F.; EI Sayed, K.A. The indole alkaloid meleagrin, from the olive tree endophytic fungus Penicillium chrysogenum, as a novel lead for the control of c-Met-dependent breast cancer proliferation, migration and invasion. Bioorgan. Med. Chem. 2016, 24, 113–122. [Google Scholar] [CrossRef] [Green Version]
  39. Kopp-Holtwiesche, B.; Rehm, H.J. Antimicrobial action of roquefortine. J. Environ. Pathol. Toxicol. Oncol. 1990, 10, 41–44. [Google Scholar]
  40. Clark, B.; Capon, R.J.; Lacey, E.; Tennant, S.; Gill, J.H. Roquefortine E, a diketopiperazine from an Australian isolate of Gymnoascus reessii. J. Nat. Prod. 2005, 68, 1661–1664. [Google Scholar] [CrossRef]
  41. Solovyev, V.; Kosarev, P.; Seledsov, I.; Vorobyev, D. Automatic annotation of eukaryotic genes, pseudogenes and promoters. Genome Biol. 2006, 7, S10. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Kumar, S.; Stecher, G.; Tamura, K. Molecular Evolutionary Genetics Analysis Version 7.0 for bigger datasets. Mol. Biol. Evol. 2016, 33, 1870–1874. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Tang, M.C.; Cui, X.; He, X.; Ding, Z.; Zhu, T.; Tang, Y.; Li, D. Late-Stage terpene cyclization by an integral membrane cyclase in the biosynthesis of isoprenoid epoxycyclohexenone natural products. Org. Lett. 2017, 19, 5376–5379. [Google Scholar] [CrossRef] [PubMed]
  44. Zilla, M.K.; Qadri, M.; Pathania, A.S.; Strobel, G.A.; Nalli, Y.; Kumar, S.; Guru, S.K.; Bhushan, S.; Singh, S.K.; Vishwakarma, R.A.; et al. Bioactive metabolites from an endophytic Cryptosporiopsis sp. inhabiting Clidemia hirta. Phytochemistry 2013, 95, 291–297. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  45. Zhang, F.Z.; Li, X.M.; Li, X.; Yang, S.Q.; Meng, L.H.; Wang, B.G. Polyketides from the mangrove-derived endophytic fungus Cladosporium cladosporioides. Mar. Drugs 2019, 17, 296. [Google Scholar] [CrossRef] [Green Version]
  46. Strobel, G.A. Endophytes as a source of bioactive products. Microbes Infect. 2003, 6, 535–544. [Google Scholar] [CrossRef]
Sample Availability: Samples of the compounds 1–4 are available from the authors.
Molecules 25 03657 g001 550
Figure 1. Phylogenetic tree analyses of HdaA in the strain Fes170 and its homologs from different species. Branch lengths are in proportion to distance.
Figure 1. Phylogenetic tree analyses of HdaA in the strain Fes170 and its homologs from different species. Branch lengths are in proportion to distance.
Molecules 25 03657 g001
Molecules 25 03657 g002 550
Figure 2. Generation and phenotype of the ΔhdaA strain. (A) Schematic illustration for hdaA disruption. The bleoR gene is amplified from the plasmid pZeo, and the bleomycin is used for the selection of transformants bearing the bleoR gene. Transformation was performed by homologous recombination using the protoplast transformation method. (B) The phenotype of the ΔhdaA and WT strains grown on PDA plates (25 °C for 5 days).
Figure 2. Generation and phenotype of the ΔhdaA strain. (A) Schematic illustration for hdaA disruption. The bleoR gene is amplified from the plasmid pZeo, and the bleomycin is used for the selection of transformants bearing the bleoR gene. Transformation was performed by homologous recombination using the protoplast transformation method. (B) The phenotype of the ΔhdaA and WT strains grown on PDA plates (25 °C for 5 days).
Molecules 25 03657 g002
Molecules 25 03657 g003 550
Figure 3. Secondary metabolite profiles of ΔhdaA and WT strains. (A) HPLC analysis of secondary metabolite profiles. (B) Relative amounts of four main products 1–4 in ΔhdaA compared with the WT. Mean values with asterisks are significant. (C) The chemical structure of the four main products detected in this study: 1, chrysogine; 2, meleagrin; 3, roquefortine C; 4, roquefortine F. The analysis for each strain was performed in triplicate. Mean values with asterisk are significant at p < 0.01.
Figure 3. Secondary metabolite profiles of ΔhdaA and WT strains. (A) HPLC analysis of secondary metabolite profiles. (B) Relative amounts of four main products 1–4 in ΔhdaA compared with the WT. Mean values with asterisks are significant. (C) The chemical structure of the four main products detected in this study: 1, chrysogine; 2, meleagrin; 3, roquefortine C; 4, roquefortine F. The analysis for each strain was performed in triplicate. Mean values with asterisk are significant at p < 0.01.
Molecules 25 03657 g003
Molecules 25 03657 g004 550
Figure 4. Transcriptional analysis of the chrysogine biosynthetic gene cluster. (A) Organization of the chrysogine biosynthetic gene cluster (BGC). (B) Quantitative RT-PCR analysis of the chrysogine BGC. The analysis for each strain was performed in triplicate. Data are shown as fold change relative to the first trial of the WT. Mean values with asterisk are significant at p < 0.01.
Figure 4. Transcriptional analysis of the chrysogine biosynthetic gene cluster. (A) Organization of the chrysogine biosynthetic gene cluster (BGC). (B) Quantitative RT-PCR analysis of the chrysogine BGC. The analysis for each strain was performed in triplicate. Data are shown as fold change relative to the first trial of the WT. Mean values with asterisk are significant at p < 0.01.
Molecules 25 03657 g004
Molecules 25 03657 g005 550
Figure 5. Transcriptional analysis of the meleagrin/roquefortine biosynthetic gene cluster. (A) Organization of the meleagrin/roquefortine BGC. (B) Quantitative RT-PCR analysis of the meleagrin/roquefortine BGC. The analysis for each strain was performed in triplicate. Data are shown as fold change relative to the first trial of the WT. Mean values with asterisk are significant at p < 0.01.
Figure 5. Transcriptional analysis of the meleagrin/roquefortine biosynthetic gene cluster. (A) Organization of the meleagrin/roquefortine BGC. (B) Quantitative RT-PCR analysis of the meleagrin/roquefortine BGC. The analysis for each strain was performed in triplicate. Data are shown as fold change relative to the first trial of the WT. Mean values with asterisk are significant at p < 0.01.
Molecules 25 03657 g005
Table
Table 1. Bioactivities of the meleagrin and roquefortine compounds from the ΔhdaA strain.
Table 1. Bioactivities of the meleagrin and roquefortine compounds from the ΔhdaA strain.
CompoundAntimicrobial Activity (MIC, μg/mL)Cytotoxicity (IC50, μM)
E. coliS. aureusC. albicansC. glabrataK562HL-60
Meleagrin12832>128>1288.912.7
Roquefortine C6416>128>12827.428.1
Roquefortine F6416>128>12822.725.1
Chloramphenicol21n.t.n.t.n.t.n.t.
Fluconazolen.t.n.t.11n.t.n.t.
Adriamycinn.t.n.t.n.t.n.t.0.30.1
Antibacterial chloramphenicol, antifungal fluconazole, and anticancer Adriamycin were used as positive control. n.t.–not tested.

Share and Cite

MDPI and ACS Style

Ding, Z.; Zhou, H.; Wang, X.; Huang, H.; Wang, H.; Zhang, R.; Wang, Z.; Han, J. Deletion of the Histone Deacetylase HdaA in Endophytic Fungus Penicillium chrysogenum Fes1701 Induces the Complex Response of Multiple Bioactive Secondary Metabolite Production and Relevant Gene Cluster Expression. Molecules 2020, 25, 3657. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25163657

AMA Style

Ding Z, Zhou H, Wang X, Huang H, Wang H, Zhang R, Wang Z, Han J. Deletion of the Histone Deacetylase HdaA in Endophytic Fungus Penicillium chrysogenum Fes1701 Induces the Complex Response of Multiple Bioactive Secondary Metabolite Production and Relevant Gene Cluster Expression. Molecules. 2020; 25(16):3657. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25163657

Chicago/Turabian Style

Ding, Zhuang, Haibo Zhou, Xiao Wang, Huiming Huang, Haotian Wang, Ruiyan Zhang, Zhengping Wang, and Jun Han. 2020. "Deletion of the Histone Deacetylase HdaA in Endophytic Fungus Penicillium chrysogenum Fes1701 Induces the Complex Response of Multiple Bioactive Secondary Metabolite Production and Relevant Gene Cluster Expression" Molecules 25, no. 16: 3657. https://0-doi-org.brum.beds.ac.uk/10.3390/molecules25163657

Article Metrics

Back to TopTop