Further Biochemical Profiling of Hypholoma fasciculare Metabolome Reveals Its Chemogenetic Diversity

Introduction

Mushroom-forming fungi are recognized for producing a plethora of chemical compounds to help defeat competitive organisms that coexist in their ecosystem (Schmidt-Dannert, 2016). These are specialized secondary metabolites (SMs) or natural products, which are known for their wide-ranging useful biological activities including antimicrobial, antitumor, and insecticide properties (Wilkins and Harris, 1944; Wilkins, 1946; de Mattos-Shipley et al., 2016). Over the recent decades, antibiotic drug discovery has mostly focussed on bacteria and ascomycetes fungi, overlooking the potential of basidiomycetes (Robbins et al., 1947). Reports in the literature state that only 20% of the current antibiotics are of fungal origin (Newman and Cragg, 2016). This is most likely because fungi primarily produce their SMs to outcompete their competitors in their environment and that many of these triggers or cues are not reproducible under laboratory conditions (Jakubczyk and Dussart, 2020). Compared to other groups of fungi, mushroom-forming basidiomycetes grow at a slow rate and reproduce using dikaryotic cell types, making genetic investigations more complex and time-consuming. Therefore, it was predicted that these basidiomycetes may produce exploitable novel compounds (Rahi and Malik, 2016). Although basidiomycetes are known to produce mycotoxins with significant biological activities with both medicinal and agricultural applications, for example psilocybin from Psilocybe spp., strobilurin from Strobilurus tenacellus, and pleuromutilin from Clitopilus passeckerianus (Bailey et al., 2016; Nofiani et al., 2018; Fricke et al., 2019), the main groups of natural products that have been isolated from basidiomycetes fungi are typically halogenated compounds and terpenoids. These include tetrachlorinated phenols, illudanes, sterpurenes, and illudalanes, all of which are largely produced exclusively by fungi, with one exception being illudalanes, which are also produced by some plants (De Jong and Field, 1997; Quin et al., 2014).

Given the number of unstudied Agaricales species, there is a high probability of finding SMs with useful biological activities, such as antibiotics, within this group of organisms. H. fasciculare (Naematoloma fasciculare) is an inedible mushroom-forming basidiomycete found on decaying wood. It is commonly named “sulfur tuft” due to its growth pattern (tight clusters or tufts of mushrooms) and the bright sulfur yellow color of its cap. Several field studies have shown the capability of this fungus to control the colonization of other wood decay organisms (De Jong et al., 1994). H. fasciculare has previously been reported in the literature as a rich source of terpenoid and organohalogen natural products including fascicularones and anisaldehyde metabolites; however, the pharmaceutical properties and the biological synthesis of these SMs remain unstudied (Supplementary Figure 1 lists all chemicals characterized from H. fasciculare from 1967 to 2019). To gain a greater understanding of the biogenetics and biochemistry of the H. fasciculare metabolome, a series of genetic manipulations, bioactivity assays, and chemical analyses of crude and pure extracts of H. fasciculare were performed to detect any potential antimicrobial activities and to further investigate the metabolic potential of this fungus. The overall aim was to attempt to bridge the gap between natural bioactive molecules and combating antibacterial resistance.

Materials and Methods

H. fasciculare Genome Mining

Our previous genomic investigation of Hypholoma suggested that only terpenoid compounds were produced, with a range of cyclization patterns (Al-Salihi et al., 2019). However, a subsequent in-depth BLAST search of functionally characterized core enzymes selected from different fungi resulted in the identification of additional biosynthetic gene clusters (BGCs) in both Hypholoma species (see Supplementary Tables 1, 2). The introns and exons of selected scaffolds were predicted using a combination of Softberry and Local BLAST searches, allowing the subsequent functional analysis of the predicted biosynthetic gene clusters. Annotation of the predicted open reading frames of each of the submitted contigs was then carried out using Artemis. Protein BLAST search on NCBI was undertaken for each gene for possible function prediction. At least 10 genes either side of the predicted SM core enzymes were annotated. The Hypholoma BGCs were then manually curated by BLAST searches against the Hypholoma sublateritium genome on JGI.

Antisense Plasmid Construction and Agrobacterium Transformation

Construction of Antisense Vector Targeting Argininosuccinate Synthetase Gene

Argininosuccinate synthetase is an essential protein for fungal growth, and successful silencing would present as a starved growth pattern, allowing a quick and simple assessment of any potential silencing. The published sequences of H. sublateritium argininosuccinate synthetase were BLAST searched against the H. fasciculare genome. A gene with 93% identity was identified in H. fasciculare contig 63. The argininosuccinate antisense plasmid consisted of a pCAMBIA0380YA backbone, 500 bp of the H. fasciculare argininosuccinate gene, which was inserted (in the antisense orientation) between the H. sublateritium gpd promoter and TrpC terminator, and the hygromycin cassette (hph gene under the Agaricus bisporus gpdII promoter and CaMV35S terminator). The verified argininosuccinate synthetase-silencing construct was moved into the Agrobacterium tumefaciens strain LBA4404 and used to transform H. fasciculare (Supplementary Table 3 lists the primers used for the construction and confirmation of the antisense plasmids).

Silencing of H. fasciculare Terpene Synthases

During the in silico analysis, it appeared that terpene synthases are the most common H. fasciculare SM gene clusters. RNA interference (RNAi)-mediated gene silencing of the core terpene synthases was performed in an attempt to link each predicted terpene synthase gene to at least one of the previously reported natural molecules from H. fasciculare. Genes were selected according to the predicted enzymatic carbon cyclization pattern, including five representatives predicted to exhibit 1,11 carbon cyclization (HfasTerp-255, HfasTerp-94A, HfasTerp-94B, and HfasTerp-105). The remaining genes were as follows: HfasTerp-147 for 1,10 3RNNP, HfasTerp-804 for 1,6 3R/S-NPP, and HfasTerp-342 and HfasTerp-179 for 1,10, E,E,farnesyl diphosphate (E,E-FPP). The atypical HfasTerp-85b was also included in this investigation.

Prior to antisense plasmid construction, reverse transcription PCR (RT-PCR) was deployed for the genes selected to confirm their predicted splicing patterns. The introns of all nine selected terpene synthases plus two housekeeping genes (gpd and β-tubulin) were predicted using a combination of SoftBerry and Artemis. RNA was extracted from growing mycelial cultures, in which their antimicrobial activity had already been confirmed. Complementary DNAs (cDNAs) were then synthesized with Oligo(dT)18 primers, and amplification of 150- to 250-bp fragments spanning at least one intron was carried out for each gene. These cDNA-derived segments of terpene synthase genes were cloned into silencing vectors as described above.

Chemical Profiling of H. fasciculare Silenced Lines

Mycelial plugs of the silenced transformants were individually inoculated into 100 ml of MEB (15 g/L malt extract broth) in a 250-ml flask and incubated at 25°C and 200 rpm for 21 days. The previously described ethyl acetate metabolite extraction protocol was used (Bailey et al., 2016). The chemical compositions of the wild type and the silenced lines (20 μl, final concentration of 5 mg/ml) of each crude extract were then compared by high-performance liquid chromatography (HPLC) as described in (Al-Salihi et al., 2019).

Expression of Selected Terpene Synthase Enzymes in Aspergillus oryzae

To avoid the potential issue associated with intron miss-splicing, full-length cDNA templates for the selected genes (HfasTerp-94A, HfasTerp94B, HfasTerp179, and HfasTerp344) were synthesized by RT-PCR. The cDNA versions of the sesquiterpene synthases (Cop-1, Cop-2, Cop-3, Cop-4, Omph-6, and Omph-7) were kindly provided by Schmidt’s group (Agger et al., 2009; Wawrzyn et al., 2012). Expression vectors were generated by yeast-based recombination as described in Al-Salihi et al. (2019). A. oryzae transformants were generated for the 10 selected enzymes and chemically analyzed using the protocol described in Al-Salihi et al. (2019).

Results

Bioassay

We assayed nine basidiomycetes to determine their ability to produce bioactive SMs on a range of solid media (see Supplementary Material for details of the method), from which the two Strophariaceae species (H. fasciculare and H. sublateritium) displayed noticeable antimicrobial activity against the three challenged microbes (see Figure 1). In contrast, Paxillus involutus showed no activity against any of the microbes tested. Variable inhibition zones were produced by the remaining basidiomycetes (see Supplementary Figures 2–7 for media, test microbes, and clearing zones diameter description).

FIGURE 1

Figure 1. (A,D) Examples of the zone inhibition plates of Hypholoma fasciculare and Hypholoma sublateritium showing the clearing zone around the fungal colony, indicating the antimicrobial activity of these fungi against Bacillus subtilis (1), Saccharomyces cerevisiae (2), and Escherichia coli (3), respectively. (B) Zone inhibition assay to evaluate the antimicrobial activity of H. fasciculare growing on different media against B. subtilis, E. coli, and S. cerevisiae. Error bars indicate the standard deviations of three technical replicate measurements for both fungal colony diameter (column in blue) and inhibition zone diameter (column in red). (E) Zone inhibition assay of H. sublateritium growing on different media against B. subtilis, E. coli, and S. cerevisiae. Error bars indicate the standard deviations of three technical replicate measurements. (C,F) Thin-layer chromatography (TLC) plates developed in a polar (H. fasciculare) and a semi-polar (H. sublateritium) system and visualized under UV 320 nm showing the separated compounds of H. fasciculare. Twenty-five microliters of 10 mg/ml crude extract of H. fasciculare from five different media was spotted. 1 = CSO supernatant extract; 2 = CSO mycelial extract; 3 = PDB supernatant extract; 4 = PDB mycelial extract; 5 = YMG supernatant extract; 6 = YMG mycelial extract; 7 = CGC supernatant extract; 8 = CGC mycelial extract; 9 = MEB supernatant extract; 10 = MEB mycelial extract. YMG, yeast extract malt; PDB, potato dextrose broth; PDA, potato dextrose agar; MEA, malt extract agar; Pep, peptone agar; YM, yeast malt extract; CDA, Czapek Dox agar; MPGA, malt peptone glucose agar.

Bioautography

Due to the broad similarity between the chemical profiles of the crude extracts of both the mycelia and supernatant of the five culture media tested for both Hypholoma strains (data not shown), it was decided to use homogenized cultures for the extraction of the metabolites. The chemical extracts were then separately spotted on the thin-layer chromatography (TLC) plates and developed in different solvent systems—non-polar, polar, and semi-polar—which were then tested against Bacillus subtilis (see Supplementary Figures 9–12 for TLC plate examples).

The antibacterial activity of the crude extracts separated using a polar system showed one large zone of inhibition (sizes ranged from 20 to 26 mm), and the highest antibacterial activity (26 mm) was observed for the crude extracts of yeast extract malt (YMG) among the other extracts of H. fasciculare and in potato dextrose broth (PDB, 25 mm) crude extract of H. sublateritium. Supplementary Figures 13–15 show some examples of the bioautography plates of H. fasciculare.

Chemical Profiling and Structure Elucidation of 3,5-Dichloro-4-Methoxybenzoic Acid

Analysis of the bioautography assays suggested similar chemical profiles for both Hypholoma species. As a result, the metabolites H. fasciculare species were purified. The variations in the retention factors (RFs) of the observed active compounds in the bioautography plates suggested that they were different metabolites. CSO-1A, YMG, CGC, and MEB were then selected for purification. Initial HPLC analysis of the crude extracts was performed to select the most potent samples for metabolite purification. This led to the detection of five main products: A, B, C, D, and E (Supplementary Figure 16). Subsequent isolation experiments were attempted to purify the five major peaks.

Despite attempting different fractionation programs, only two fractions (A and D) could be purified from the crude extracts. A comparison was then carried out between the mass spectra of the purified metabolites (A and D) and the compounds previously reported from the Hypholoma genus, leading to the characterization of metabolite A as fascicularone G (a known sesquiterpene from Hypholoma). However, for fraction D, no match could be found with any of the previously characterized metabolites from Hypholoma. We therefore subjected fraction D to high-resolution mass spectrometry (HRMS) to confirm its exact mass and predict its molecular formula (see Figure 2). This analysis indicated an aromatic molecule of chemical formula C₈H₆O₃Cl₂ to better than 2 ppm mass accuracy. This is putatively assigned as 3,5-dichloro-4-methoxybenzoic acid (3,5-D). The isotope distribution of the [M-H]^– ion also matches the predicted distribution. The only fragment ion in the spectrum is due to decarboxylation (loss of CO₂), which is typical of an aromatic carboxylic acid in negative ion mode, which provides further evidence for the assignment of metabolite D.

FIGURE 2

Figure 2. Accurate-mass negative ion electrospray ionization (ESI) mass spectrum of fraction C—putative 3,5-dichloro-4-methoxy benzoic acid. The spectrum shows the [M-H]^– at m/z 218.9625 (218.9621 for C₈H₅O₃Cl₂^–) and fragmentation resulting from decarboxylation at m/z 174.9726 (174.9723 for C₇H₅OCl₂^–).

The structure of compound D was further investigated using ¹H and ¹³C NMR, where the presence of a symmetrical benzoic acid molecule, two chlorines, and one methoxy group was confirmed. The correlations between the aromatic protons at 7.9 and 3.8 ppm and the carbon signals were also investigated by heteronuclear multiple bond correlation (HMBC), and the structure of 3,5-dichloro-4-methoxybenzoic acid was confirmed (see Supplementary Figures 17–21 for the NMR characterization).

Antimicrobial Activity of 3,5-D

The antibacterial activity of 3,5-D was assessed by disk diffusion assay using gradient concentrations of 100, 300, and 500 μg. Disks containing methanol and kanamycin were used as negative and positive control, respectively. As shown in Table 1, the antibacterial activity of 3,5-D increased with an increase of its concentration against Staphylococcus aureus and B. subtilis.

TABLE 1

Table 1. Antimicrobial activity of 3,5-dichloro-4-methoxybenzoic acid (3,5-D) tested against Bacillus subtilis and Staphylococcus aureus.

Genomic Investigation

Analysis of the H. fasciculare genome using FungiSmash identified 15 pairs of orthologous BGCs for H. fasciculare and H. sublateritium, as outlined below (detailed description of all predicted BGCs can be found in Supplementary Figures 22–30).

(1) HfasTerp-94: This biosynthetic gene cluster is located in contig 94 of the H. fasciculare genome. Its terpene synthase was predicted as a member of the 1,11 E,E-FPP carbon cyclization enzyme group. This group of enzymes is responsible for the production of trans-humulyl cation sesquiterpenes, such as humulene, caryophyllene, and protoilludane (Quin et al., 2014). When the sequences of these genes were BLAST searched against the H. sublateritium genome, an analogous biosynthetic gene cluster was found, with one exception—SDR3, which could not be found (Figures 3A–C).

FIGURE 3

Figure 3. Selected sesquiterpene biosynthetic gene clusters of 1,11 E,E-FPP carbon cyclization clade identified in the Hypholoma fasciculare genome and their likely initial products. (A) HfasTerp-94a and HfasTerp-94b biosynthetic gene clusters and examples of likely products—caryophyllene, humulene, and protoilludane—illustrating the 1,11 carbon cyclization pattern predicted for this synthase. 1 = Serine/threonine kinase; 2 = Zinc carboxypeptidase; 3 = Hypothetical protein; 4 = Short-chain dehydrogenase 1; 5 = Short-chain dehydrogenase 2; 6 = Short-chain dehydrogenase 3; 7 = Tyrosinase; 8 = F-box domain; 9 = Anchor signaling protein; 10 = Terpene synthase B; 11 = Splicing co-activator; 12 = F-box domain; 13 = Glucose transporter; 14 = Peptide transporter; 15 = Terpene synthase A; 16 = Aromatic ring hydroxylase; 17 = Glycoside hydrolase. (B) HfasTerp-179 biosynthetic gene cluster (BGC) and examples of likely products arising from the predicted 1,10 carbon cyclization for this synthase. 1 = Alpha/beta-hydrolases; 2 = Aspartyl protease; 3 = E3 ubiquitin–protein ligase; 4 = NAD-dependent histone deacetylase; 5 = Hypothetical protein; 6 = Transcription factor; 7 = GMC oxidoreductase; 8 = Hypothetical protein; 9 = Cytochrome P450; 10 = Terpene synthase; 11 = GMC oxidoreductase; 12 = Fructose 2,6-bisphosphatase; 13 = Kinase like protein; 14 = Hypothetical protein; 15 = DNA repair endonuclease; 16 = Telomerase transcriptase; 17 = Monooxygenase FAD; 18 = Glycoside hydrolase; 19 = Alpha ketogluturate. (C) HfasTerp-804 BGC predicted to be responsible for the production of the antitumor clavaric acid or similar derivatives. 1 = Hypothetical protein; 2 = Flavin-containing amine oxidase; 3 = Pyruvate decarboxylase; 4 = Cytoplasmic protein; 5 = Phosphomannomutase; 6 = Delta DNA polymerase; 7 = Putative transcription factor; 8 = 60S ribosomal protein translation; 9 = Actin-related protein; 10 = Oxidosqualene clavarinone cyclase; 11 = Nuclear condensin complex.

(2) HfasTerp-179: This biosynthetic gene cluster was identified in contig 179 of H. fasciculare. This sesquiterpene synthase clade is identified as those based on germacradienyl derivatives. A BLAST search of the core enzyme and its adjacent biosynthetic genes against H. sublateritium genome indicated a matching cluster in H. sublateritium scaffold 11, but this appeared to be a larger cluster with two additional genes, likely beyond the end of the HfasTerp179 contig. Searching the Hfas genome with the two additional genes from the H. sublateritium scaffold 11 cluster identified these genes at the start of contig 615 of H. fasciculare, indicative of an incomplete genome assembly.

(3) HfasTerp-804: This cluster consisted of the terpene synthase that demonstrated high sequence similarity with the experimentally characterized oxidosqualene synthase located in scaffold 133 of H. sublateritium, a key enzyme for the production of the antitumor compound clavaric acid in H. sublateritium (Godio and Martín, 2009).

(4) HfasTerp-255: Manual sequence annotation of H. sublateritium scaffold 30 predicted three terpene synthases: HsubTerp-30a, HsubTerp-30b, and HsubTerp-30c. Phylogenetic analysis suggested HsubTerp-30a to function as squalene synthase and HsubTerp-30c as potential protoilludane synthase. BLAST searches of these genes against the H. fasciculare genome revealed high sequence similarity with several terpene synthases; the highest similarity (88%) was observed between HfasTerp-255 and HsubTerp-30c (Figure 4).

FIGURE 4

Figure 4. Maximum likelihood tree of putative terpene synthases from Hypholoma fasciculare and Hypholoma sublateritium. Contigs or scaffold numbers are shown adjacent to species abbreviations. The evolutionary history was inferred by using the maximum Likelihood method and the Whelan and Goldman model conducted in MEGA X. The percentages of trees in which the associated taxa clustered together are shown next to the branches.

(5) HfasTerp-147: This BGC includes the only predicted terpene enzyme that follows the 1,10 3R neryl diphosphate (NPP) cyclization pattern. A homologous BGC was located in scaffold 11 of the H. sublateritium genome. Subsequent analysis of the tailoring genes of the H. sublateritium biosynthetic gene cluster suggested that the HfasTerp-147 gene cluster was assembled into two different contigs: HfasTerp-458 and HfasTerp-147 (Figure 5A).

FIGURE 5

Figure 5. Schematic illustration of the organization of the putative biosynthetic gene clusters for both HfasTerp-147 and HfasTerp-104. (A) The extended HfasTerp-147 biosynthetic gene cluster. 1 = Complex I intermediate-associated protein; 2 = Protein binding; 3 = F-box; 4 = Histone methyltransferase; 5 = RNA polymerase; 6 = Terpene synthase; 7 = Cystathionine beta-lyase; 8 = Pyruvate carboxyltransferase; 9 = Dolichol kinase; 10 = Zinc finger; 11 = Splicing factor; 12 = Proliferation-associated protein; 13 = Chaperone regulator; 14 = Enoyl-CoA hydratase; 15 = Aspartic peptidase; 16 = Aspartic peptidase. (B) Organization of the putative biosynthetic gene cluster likely responsible in producing 3,5-dichloro-4-methoxybenzoic acid (3,5-D) in both Hypholoma fasciculare and the homologous cluster of Hypholoma sublateritium. 1 = Glycoside; 2 = Glutathione S-transferase; 3 = Anchor protein; 4 = Oxoprolinase; 5 = Monocarboxylate; 6 = Glycoside hydrolase; 7 = F-box; 8 = Anchor protein; 9 = Hypothetical protein; 10 = SDR; 11 = Calmodulin-related protein; 12 = Hypothetical protein; 13 = Terpene synthase; 14 = Regulator of G protein signaling; 15 = Benzoic acid reductase-PKS; 16 = 3-Phosphoshikimate-1-carboxyvinyltransferase.

(6) The 3,5-D putative BGC: Among the shared hybrid biosynthetic gene clusters of Hypholoma spp., HfasTerp-104 and HsubTerp-38 appeared to have the necessary enzymes for the synthesis of the compound 3,5-D, such as 3-phosphoshikimate-1-carboxyvinyltransferase, benzoic acid reductase-polyketide synthase (PKS), short-chain dehydrogenase reductase (SDR), and glycoside hydrolase, highlighting their likely role in halogenated natural product synthesis (Figure 5B).

Gene Clusters for Other Classes of Secondary Metabolites

(7) HfasPKS-221: This biosynthetic gene cluster was positioned in contig 221 of H. fasciculare. Subsequent comparison with the H. sublateritium genome revealed an identical BGC situated in scaffold 53 (Figure 6A).

FIGURE 6

Figure 6. (A) Schematic of putative HfasPKS-221 (74 Kb) biosynthetic gene cluster. 1 = Triphosphate hydrolase; 2 = Multidrug resistance; 3 = Hypothetical protein; 4 = Alpha-beta hydrolase; 5 = PKS1; 6 = Zinc finger H2C2; 7 = Kinesin-like protein; 8 = Hypothetical protein; 9 = Actin binding protein; 10 = DNA binding protein; 11 = Acetylcholinesterase; 12 = Actin depolymerizing factor; 13 = Metalloproteases; 14 = Metalloproteases; 15 = Metalloproteases; 16 = N-acetylglucosaminyltransferase; 17 = Transcription factor. (B) Schematic of putative HfasNRPS-29 (76 Kb) biosynthetic gene cluster. 1 = mRNA-decapping enzyme; 2 = Peptidase aspartic; 3 = Glycoside hydrolase; 4 = Polyadenylate binding protein; 5 = Mitochondrial carrier protein; 6 = Ubiquitin protease; 7 = Acid phosphatase; 8 = AMP-dependent synthetase and ligase; 9 = Exocyst complex; 10 = DUF domain protein; 11 = Carbohydrate transporter; 12 = Adenine DNA glycosylase; 13 = Cargo transport protein; 14 = WD repeat-containing nucleolar rRNA; 15 = Oxidoreductase; 16 = NRPS; 17 = Acyl-CoA N-acyltransferase. A homologous biosynthetic gene cluster was found in scaffold 100 of the Hypholoma sublateritium genome. (C) Schematic of putative HfasSidA-14 (98 Kb) biosynthetic gene cluster. 1 = Transcription factor; 2 = Sterol desaturase; 3 = Cytochrome P450; 4 = Triphosphate hydrolase; 5 = Calcium transporting; 6 = Glycosyltransferase; 7 = GMC oxidoreductase; 8 = Amino acid transporter; 9 = Zinc binding protein; 10 = ABC transporter; 11 = Sid A; 12 = Glycoside hydrolase; 13 = Transporter; 14 = RNA-associated protein; 15 = F-box.

(8) HfasNRPS-29: This biosynthetic gene cluster was identified in contig 29 of H. fasciculare, and its ortholog cluster was found in scaffold 100 of the H. sublateritium genome (Figure 6B).

(9) HfasSid-14: This biosynthetic gene cluster was predicted in contig 14 of H. fasciculare and its ortholog located in scaffold 11 of H. sublateritium (Figure 6C).

Silencing Experiments

One method to validating gene function would be to silence the expression of each target gene and see whether there was a resulting depletion of a corresponding metabolite. Silencing was first assessed against the gene argininosuccinate synthetase. The silencing cassettes (Figure 7A) were transferred into H. fasciculare using Agrobacterium-mediated transformation following the protocol described in Al-Salihi et al. (2017). Silencing efficiency was assessed by comparing the hyphal growth pattern of the mutants with the wild type, with and without the external addition of arginine.

FIGURE 7

Figure 7. (A) Schematic representing the antisense vector pCAMHsgpdHfas used for targeting argininosuccinate synthetase in Hypholoma fasciculare. (B) H. fasciculare wild type and antisense transformant 14 showing differences in the colony growth rate on potato dextrose agar (PDA) with and without arginine supplementation. 1 = H. fasciculare wild type on PDA medium; 2 = H. fasciculare wild type on PDA medium supplemented with 5 mM of arginine; 3 = H. fasciculare antisense transformant 49 on PDA medium; 4 = H. fasciculare antisense transformant 49 on PDA medium supplemented with 5 mM of arginine.

Several transformants showed a reduction in their colony size, among which transformant Hfas-asTR14 displayed the lowest growth rate (20 mm) compared to the wild type (30 mm), suggesting potential silencing in those lines. To further analyze this phenotype, replicates of transformant Hfas-asTR14 and the wild type (WT) were subcultured on potato dextrose agar (PDA) plates with and without 5 mM of the arginine supplement. Following 2 weeks of incubation at 25°C, the colony diameters were measured, and the results indicated potential successful gene silencing in transformant Hfas-asTR14, whereby silenced colonies demonstrated a slower rate of growth and different form of mycelia compared to the wild type (Figure 7B).

Following the argininosuccinate silencing experiments, terpene synthase silencing transformation was carried out in a similar manner. Classical selection on supplemented PDA plates was performed, and two selected silenced transformants for each line including two argininosuccinate synthetase transformants were further investigated for their biological activity using plate-based bioassay. Despite standardizing the size of the inoculum and culture conditions for all transformants, the silenced lines displayed variations in terms of colony diameter and inhibition zone. The observed morphogenesis of the silenced colonies was rigid, fluffy, wrinkled, condensed, and with loss of aerial mycelium (see Supplementary Figure 31).

Of nine silenced lines, HfasTerp-94bTR6 and HfasTerp-85bTR9 displayed the highest reduction in their zones of inhibition, 38 and 45%, respectively, while transformant HfasTerp-804TR8 and the argininosuccinate antisense Hfas-asTR49 produced larger inhibition zones (25%) compared to the wild type. Reductions from 12 to 32% were demonstrated by the remaining transformants (Table 2).

TABLE 2

Table 2. Average colony and clearing zone diameters of two selected putative antisense transformants alongside the wild type.

Gene Expression Analysis of H. fasciculare Silenced Lines

Genes HfasTerp-94a, HfasTerp-94b, and HfasTerp-105, gpd, and β-tubulin were used to detect their expression levels in the selected silenced transformants alongside the wild type. All transformants displayed reductions in their expression levels, attributed to the successful downregulation of their corresponding genes. Table 3 shows the fold differences of the selected genes over the conserved ones.

TABLE 3

Table 3. RT-qPCR outcomes of the silenced lines.

Chemical Analysis of Silenced Lines

Different levels of SM production were observed among the silenced lines, particularly in transformants Hfas-asTR49, HfasTerp85bTR2, and HfasTerp85bTR9, where the production of a number of the molecules was reduced, including fascicularone G and naematoline. However, the production of the newly characterized (in H. fasciculare) 3,5-D showed no reduction in all transformants, indicating the involvement of a different type of key enzyme in its biological synthesis (Supplementary Figure 32 and Table 4).

TABLE 4

Table 4. Growth pattern and mass detector response count per second of two selected putative antisense transformants for each silenced line alongside the wild type.

Heterologous Expression of Selected Sesquiterpene Synthases

Although silencing constructs has been proven successful for functional studies in H. fasciculare, its role in linking sesquiterpene metabolites to their specific biogenetic genes was inconclusive. We therefore adapted the vector pTYAGS-arg to express the selected sesquiterpene synthases in A. oryzae in order to further assess whether using A. oryzae as the expression host, as well as whether using different isolation techniques, would affect the measurement of the expression outcomes of some selected genes. A. oryzae transformants from a previous work were used (Al-Salihi et al., 2019), as well as six previously characterized fungal sesquiterpene synthases from two different basidiomycetes, Omphalotus olearius and Coprinopsis cinerea, as controls: cop-1, cop-2, cop-3, cop-4, omp-6, and omp-7 (Agger et al., 2009; Wawrzyn et al., 2012). Seven verified constructs—pTYGS-argHfasTerp179, pTYG-argCop-1, pTYG-argCop-2, pTYG-argCop-3, pTYG-argCop-4, pTYG-argOmp6, and pTYG-argOmp7—were independently transferred to A. oryzae strain NSAR1 protoplasts via polyethylene glycol (PEG)-mediated transformation reactions. Subsequent gas chromatography–mass spectrometry (GC-MS) analysis on previous (pTYGS-argHfasTerp94a, pTYGS-argHfasTerp94b, and pTYGS-argHfasTerp344) and current transformants confirmed our previous results of HfasTerp94a, HfasTerp94b, and HfasTerp344 (Al-Salihi et al., 2019) through the replication of humulene and caryophyllene production, although at lower levels. The production of most of the major products of the two controls, pTYG-argCop-3 and pTYG-argCop-4, were confirmed as well as the previously uncharacterized humulene in the Cop-3 extract (Figure 8 and Supplementary Figure 33). However, no product could be observed for the HfasTerp179 A. oryzae transformant.

FIGURE 8

Figure 8. Gas chromatography–mass spectrometry (GC-MS) comparison of the sesquiterpene synthase Cop4. (A) GC-MS analysis of the NSAR1-Cop4 transformant. Eight compounds [muurola-4(14)diene, humulene, naphthalene, gleenol, cubenene, bicylosesquiphellandrene, Cubenol, and alpha-copaene] were characterized using HP-5 MS quartz capillary column (30 m × 0.25 mm, 0.25-μm film thickness). (B) GC-MS analysis of the NSAR1-Cop4 transformant. Three compounds (alpha-copaene, Cubenol, and sigma-muurolene) were characterized using HP-1 (50 m × 0.32 mm, 0.17-μm phase thickness).

Co-expression and Chemical Analysis of Selected Biosynthetic Genes of the HfasTerp94 Gene Cluster

Adjacent genes (SDR1, SDR2, SDR3, and tyrosinase) of HfasTerp94 were selected for co-expression with NSAR1-humulene synthase. Due to the unsuccessful attempts of full-length cDNA amplification of the selected genes, an alternative approach of fragment amplification was chosen. In silico analysis predicted two exons for each SDR gene. Accordingly, four pairs of primers with 60 bp were used to amplify the two exons of each SDR from H. fasciculare genomic DNA (gDNA). A pTYGS-arg backbone was used for fragment recombination (Supplementary Figure 33). However, due to the prediction of several introns within the tyrosinase gene, a synthetic version was used.

Following successful transformation of the selected genes into A. oryzae, mass spectrum comparison between the five generated transformants (NSAR1-humulene synthase-SDR1, NSAR1-humulene synthase-SDR2, NSAR1-humulene synthase-SDR3, NSAR1-humulene synthase-SDR1-SDR2, and NSAR1-humulene synthase-SDR1-SDR2-Tyrosinase) and NSAR1-humulene synthase, crude extracts were evaluated; in total, seven new peaks were produced. Since the analysis was performed in electrospray ionization (ESI) negative mode, it was assumed that the observed m/z values would be 46 mass units higher due to the formation of a formic acid adduct. For the humulene synthase-SDR1 transgenic, one new peak (metabolite 1) was observed at 14.47 min, with an m/z of 489. In addition to metabolite 1, the chromatograms of both NSAR1-humulene synthase-SDR2 and NSAR1-humulene synthase-SDR3 produced three new peaks, eluting at 12.90 min (metabolite 2), 13.30 min (metabolite 3), and 14.90 min (metabolite 4), with m/z 487, 489, and 473, respectively. Although the chemical profile of the transgenic NSAR1 with humulene synthase-SDR1-SDR2 showed no accumulation for both metabolites 3 and 4, new peaks at 15.32 and 15.98 min with m/z 487 and 295, respectively, were observed. The fragmentation patterns of metabolites 1 and 2 were identical to metabolites 3 and 5, respectively, indicative that they may be isomeric compounds (Figures 9, 10).

FIGURE 9

Figure 9. Negative ion LC-MS spectrum (-ESI) comparison of NSAR1-WT with transgenics of single, double, and triple biosynthetic genes of the humulene biosynthetic gene cluster. The genes investigated in this experiment were: one terpene cyclase (humulene synthase), three genes of short-chain dehydrogenase reductase (SDR), and one tyrosinase, all selected from the terpene synthase biosynthetic gene cluster located in contig 94 of the Hypholoma fasciculare genome. (A) NSAR1-WT. (B) NSAR1-humulene synthase. (C) NSAR1-humulene synthase-SDR1. (D) NSAR1-humulene synthase-SDR2. (E) NSAR1-humulene synthase-SDR3. (F) NSAR1-humulene synthase-SDR1-SDR2. (G) NSAR1-humulene synthase-SDR1-SDR2-Tyrosinase. The mass spectrum was selected from min 11 to min 17.50 to avoid peaks overlapping. In total, seven new peaks were detected at retention times (RTs) of 12.87, 13.30, 14.47, 14.90, 15.32, 15.98, and 11.90 within this comparison. The transgenic NSAR1-humulene synthase-SDR1 showed the lowest number of new peaks compared to the other transgenics: NSAR1-humulene synthase-SDR2, NSAR1-humulene synthase-SDR3, and NSAR1-humulene synthase-SDR1-SDR2.

FIGURE 10

Figure 10. Negative ESI mass spectra of NSAR1 transgenics. Seven new metabolites were detected as formate adduct [M++HCO₂] at: metabolite 1, m/z 490; metabolite 2, m/z 311; metabolite 3, m/z 489; metabolite 4, m/z 473; metabolite 5, m/z 487; metabolite 6, m/z 295; and metabolite 7, m/z 279.

Genome Accession Number

The genome sequences of H. fasciculare can be freely accessed on the NCBI platform under the accession number JAEKIV010000000.

Discussion

The continued rise in antibiotic resistance has prompted scientists to look for new ways to discover potential new drugs. While mushroom-forming fungi have been reported as a rich repository of structurally diverse natural products, the limited access to their metabolomic and genomic data poses a bottleneck to the development of their potential. In this study, the prospective of selected fungi as a source of novel antibiotics is investigated. Eight of the nine species selected displayed antimicrobial properties, of which Hypholoma showed a notable range of antimicrobial activity. This is reflected in its metabolic variety and the flexibility in their response to various competitive ecological niches (Boddy and Watkinson, 1995), which could be potentially important in defeating resistant microbes.

Direct bioautography and subsequent chromatographic, mass spectral, and NMR spectroscopic analysis enabled us to unveil the antibacterial potency of an aromatic compound, 3,5-D. This molecule was first characterized from Stropharia sp. through the observation by Thines et al. (1995) of two aromatic chemical shifts (aromatic proton and the methoxy group) and the carbon they are attached to at 130.4 and 155.8 ppm, respectively. Here, we fruitfully used both 1D and 2D spectroscopic analyses to provide solid proof of the presence of the other substituent groups (Cl, CH₃, and COOH) and to allow the first elucidation of 3,5-D in Hypholoma genus and second within fungi kingdom. Previous studies have reported that chlorinated compounds are synthesized through the shikimate pathway (Jensen et al., 1994; Mester et al., 1997; Hage et al., 1999). This pathway involves the conversion of the phenylalanine precursor to benzoic acid derivatives via sequential condensation, hydroxylation, and chlorination. Interestingly, the predicted HfasTerp104 gene cluster includes enzymes that are likely responsible for the synthesis of 3,5-D, including the benzoic acid reductase-PKS, SDR, glycoside hydrolase, and the multifunctional 3-phosphoshikimate-1-carboxyvinyltransferase. 3-Phosphoshikimate-1-carboxyvinyltransferase is a multidomain enzyme with a main role in catalyzing the conversion of phenylalanine-like compounds to cinnamic acid derivatives. Subsequent hydroxylation and reduction on the resulting acids lead to the production of anisaldehyde isomers such as 3,5-D (Field et al., 1996). The biological activity of chlorinated natural products is well documented (Hautzel and Anke, 1990; Pfefferle et al., 1990; Becker et al., 1994). Identifying their BGCs is a step forward in the biochemical manipulation of such understudied potential potent antimicrobial agents.

Unlike H. sublateritium, the H. fasciculare genome was assembled into greater number of short contigs. It is likely that our assembly method dispersed the resulted gene clusters of H. fasciculare into several short contigs rather than one long scaffold, as the case of its related H. sublateritium genome assembled by JGI. Nevertheless, phylogenetic comparisons and in-depth manual curation with the H. sublateritium genome sequences guided the annotation of a set of hidden biosynthetic gene clusters in both species. Most of the 15 putative gene clusters contain key genes encoding terpene, PKS, non-ribosomal peptide synthase (NRPS), Sid, and their related modifying genes, which represent the main features of fungal SM biosynthetic pathways. The in silico analysis suggests that the predicted terpene clusters are likely responsible for the production of terpene-like compounds, particularly sesquiterpenes and triterpenes (previously reported compounds in Hypholoma spp.). A putative siderophore biosynthetic pathway and a multimodule NRPS core enzyme were also found in the H. fasciculare genome, and subsequent manual curation led to the identification of a putative biosynthetic cluster in both Hypholoma species. Most often, multidomain NRPSs are involved in the synthesis of SM with medicinal properties (Sayari et al., 2019), indicative that the Hypholoma species may produce NRPSs and siderophore-based metabolites that require further investigation. Several hybrid clusters containing two different cyclase enzymes (e.g., a predicted PKS and a terpene synthase) were also distinguished. Such biosynthetic clusters often coordinate the production of meroterpenoid molecules such as melleolide (Engels et al., 2011), as well as the newly characterized 3,5-D. Further annotation of the clusters suggested that the PKS of Hfasterp104 is actually a benzoic acid reductase that produces benzoate precursors. This enzyme, along with the multifunctional gene 3-phosphoshikimate-1-carboxyvinyltransferase, catalyzes five sequential steps in the shikimate pathway and is the key biosynthetic producer of phenolic compounds in plants (Santos-Sánchez et al., 2019).

Gene silencing was performed to downregulate the expressions of some selected terpene synthetase genes (alongside argininosuccinate synthetase), demonstrating that RNA silencing is an effective tool in knocking down transcripts of argininosuccinate. The silenced transformants showed variability in their phenotypes, represented by the slow growth rates or the production of unusual forms of mycelia. These phenotypic variations are congruent with other published silencing works in basidiomycetes (Eastwood et al., 2008; Costa et al., 2009). However, extended silencing work on sesquiterpene synthases reveals a crossed role of such enzymes in the synthesis of the H. fasciculare SMs, where the production of several SMs was prevented better as a result of silencing one enzyme. For example, silencing HfasTerp84b resulted in the absence of most metabolites that were observed in the chemical profile of the WT. This promiscuous feature of basidiomycetes sesquiterpene synthase has been reported previously (Agger et al., 2009; Wawrzyn et al., 2012), where independent heterologous expressions of a group of sesquiterpene synthases from two different basidiomycetes, O. olearius and C. cinerea, in E. coli resulted in several chemically distinguishable sesquiterpene molecules for each enzyme. Interestingly, the GC-MS analysis of A. oryzae transgenics of some of the enzymes used in Schmidt’s work (Agger et al., 2009; Wawrzyn et al., 2012) confirmed the fact that A. oryzae is a feasible tool in linking SM to their specific enzymes and that its metabolomic system has no effect on the final product of sesquiterpene core enzymes, as most major products of Schmidt’s group were replicated here. The absence of the remaining compounds from the observation in the other Aspergillus transgenics is likely due to the use of a different isolation method, as the detection of compounds is influenced by the isolation method and the chromatographic system used during the characterization process (Jeleń, 2003).

The detection of humulene in the GC-MS analysis of HfasTerp94 motivated the further expression of its related biosynthetic genes in A. oryzae. In this cluster, three SDR genes were predicted alongside other distinctive genes including zinc carboxypeptidase and tyrosinase. SDR genes belong to a large family of oxidoreductases that have a key role in oxidoreductase reactions through synthesizing the oxidation of metabolites (Kallberg et al., 2010). Co-expression of SDR genes with the humulene synthase was therefore hypothesized to result in the production of oxidized humulene derivatives. Different combinations of genes (dual and triple gene vectors) were therefore transformed into NSAR1 to potentially reveal the function of each enzyme. The crude chemical extract of the positive Aspergillus transformants were then analyzed via GC-MS, leading to the detection of several new peaks. This suggested that the SDR and tyrosinase protein sequences were correctly annotated and that they can further modify the humulene backbone. Based on previously reported results, humulene can be directly deprotonated to produce sesquiterpene by-products (Agger et al., 2009). Although the observation of two new metabolites (5 and 6) in the chemical profile of NSAR1-SDR1-SDR2-Tyrosinase may indicate the potential conversion of metabolites 3 and 4 of NSAR1-SDR1-SDR2, further evidence of the compound’s identity (such as a full NMR characterization) is needed to precisely link the new SMs to the corresponding genes. It is also worth noting that the expressions of some modifying enzymes (e.g., cytochrome P450 and SDR) may impact the metabolic pathway of the utilized host, A. oryzae (Shinohara et al., 2016), and that the compounds observed are the result of combinatorial efforts between the inserted genes and the native genes of the host. This hypothesis, however, needs further laboratory-based evidence.

Conclusion

Mushroom-forming fungi offer the potential for the production of new molecules with novel biological activities that have the potential to be medicinally useful, especially in the development of new antibiotics. However, only a small number of these fungi has been fully explored at both the genetic and chemical levels. To address problems associated with SM production in H. fasciculare (low yield and difficulties in growth), an in-depth sequence analysis of the H. fasciculare genome was performed. This led to the characterization of several putative biosynthetic gene clusters that are likely to be involved in the production of the previously reported compounds in this fungus. When we experimentally tested our in silico genome prediction via RNAi-mediated and heterologous expression in A. oryzae, we further demonstrated that RNAi silencing is still an efficient tool in triggering the downregulation of gene expression and could be deployed in other genetic manipulations as in H. fasciculare and that sesquiterpene synthases are promiscuous enzymes. The study also permitted the chemical characterization of a biologically active meroterpenoid (3,5-D) as well as the in silico prediction of the biosynthetic gene cluster associated with it. Prominently, this cumulative knowledge of the bioactive molecules’ structures and their biosynthesis and pharmacokinetic properties will help in developing a molecular framework to prolifically employ this fungus in the battle against antimicrobial resistance.

Data Availability Statement

The raw data supporting the conclusions of this article will be made available by the authors, without undue reservation.

Author Contributions

SA-S designed and performed the experiments, analyzed data, and wrote the manuscript. IB provided some GC-MS samples analysis. RA-S and KS participated in solving the chemical structure of 3,5-D compound. PG participated in 3,5-D HPLC run and provided feedback on its fragmentation. AB and GF devised and supervised the project. All authors contributed to the article and approved the submitted version.

Funding

SA-S was supported financially by the Ministry of Higher Education and Scientific Research, Iraq. All research was conducted at the University of Bristol, United Kingdom.

Conflict of Interest

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.

Acknowledgments

We thank Fabrizio Alberti, Claudio Greco, and Colin Lazarus for their assistance throughout the study. This manuscript has been released as a preprint at biorxiv, doi: https://doi.org/10.1101/2020.05.28.122176 (Al-Salihi et al., 2020). Further biochemical profiling of the Hypholoma fasciculare metabolome reveals its chemogenetic diversity (Al-Salihi et al., 2020).

Supplementary Material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fbioe.2021.567384/full#supplementary-material

References

Agger, S., Lopez-Gallego, F., and Schmidt-Dannert, C. (2009). Diversity of sesquiterpene synthases in the basidiomycete Coprinus cinereus. Mol. Microbiol. 72, 1181–1195. doi: 10.1111/j.1365-2958.2009.06717.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Salihi, S. A., Bull, I., Al-Salhi, R. A., Gates, P. J., Salih, K., Bailey, A. M., et al. (2020). Further biochemical profiling of Hypholoma fasciculare metabolome reveals its chemogenetic diversity. biorxiv doi: 10.1101/2020.05.28.122176

CrossRef Full Text | Google Scholar

Al-Salihi, S. A., Dao, T. T., Williams, K., Bailey, A. M., and Foster, G. D. (2019). The Biogenetic Origin of the Biologically Active Naematolin of Hypholoma Species Involves an Unusual Sesquiterpene Synthase. Mol. Biotechnol. 61, 754–762. doi: 10.1007/s12033-019-00199-x

PubMed Abstract | CrossRef Full Text | Google Scholar

Al-Salihi, S. A., Scott, T. A., Bailey, A. M., and Foster, G. D. (2017). Improved vectors for Agrobacterium mediated genetic manipulation of Hypholoma spp. and other homobasidiomycetes. J. Microbiol. Methods 142, 4–9. doi: 10.1016/j.mimet.2017.08.014

PubMed Abstract | CrossRef Full Text | Google Scholar

Bailey, A. M., Alberti, F., Kilaru, S., Collins, C. M., de Mattos-Shipley, K., Hartley, A. J., et al. (2016). Identification and manipulation of the pleuromutilin gene cluster from Clitopilus passeckerianus for increased rapid antibiotic production. Sci. Rep. 6:25202.

Google Scholar

Becker, U., Anke, T., and Sterner, O. (1994). A novel bioactive illudalane sesquiterpene from the fungus Pholiota destruens. Nat. Prod. Lett. 5, 171–174. doi: 10.1080/10575639408044054

CrossRef Full Text | Google Scholar

Boddy, L., and Watkinson, S. C. (1995). Wood decomposition, higher fungi, and their role in nutrient redistribution. Can. J. Bot. 73, 1377–1383. doi: 10.1139/b95-400

CrossRef Full Text | Google Scholar

Costa, A. S., Thomas, D. J., Eastwood, D., Cutler, S. B., Bailey, A. M., Foster, G. D., et al. (2009). Quantifiable downregulation of endogenous genes in Agaricus bisporus mediated by expression of RNA hairpins. J. Microbiol. Biotechnol. 19, 271–276.

Google Scholar

De Jong, E., and Field, J. A. (1997). Sulfur tuft and turkey tail: biosynthesis and biodegradation of organohalogens by basidiomycetes. Annu. Rev. Microbiol. 51, 375–414. doi: 10.1146/annurev.micro.51.1.375

PubMed Abstract | CrossRef Full Text | Google Scholar

De Jong, E. D., Field, J. A., Spinnler, H. E., Wijnberg, J. B., and de Bont, J. A. (1994). Significant biogenesis of chlorinated aromatics by fungi in natural environments. Appl. Environ. Microbiol. 60, 264–270. doi: 10.1128/aem.60.1.264-270.1994

PubMed Abstract | CrossRef Full Text | Google Scholar

de Mattos-Shipley, K. M. J., Ford, K. L., Alberti, F., Banks, A. M., Bailey, A. M., and Foster, G. D. (2016). The good, the bad and the tasty: the many roles of mushrooms. Stud. Mycol. 85, 125–157. doi: 10.1016/j.simyco.2016.11.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Eastwood, D. C., Challen, M. P., Zhang, C., Jenkins, H., Henderson, J., and Burton, K. S. (2008). Hairpin-mediated down-regulation of the urea cycle enzyme argininosuccinate lyase in Agaricus bisporus. Mycol. Res. 112, 708–716. doi: 10.1016/j.mycres.2008.01.009

PubMed Abstract | CrossRef Full Text | Google Scholar

Engels, B., Heinig, U., Grothe, T., Stadler, M., and Jennewein, S. (2011). Cloning and characterization of an Armillaria gallica cDNA encoding protoilludene synthase, which catalyzes the first committed step in the synthesis of antimicrobial melleolides. J. Biol. Chem. 286, 6871–6878. doi: 10.1074/jbc.m110.165845

PubMed Abstract | CrossRef Full Text | Google Scholar

Field, J. A., Verhagen, F. J. M., Mester, T., Swarts, H. J., Wijnberg, J. B. P. A., and de Jong, E. (1996). “Organohalogen Metabolites of Basidiomycetes,” in Mechanisms of Biohalogenation and Dehalogenation, ed. R. Wever (North-Holland: Amsterdam), 33–43.

Google Scholar

Fricke, J., Lenz, C., Wick, J., Blei, F., and Hoffmeister, D. (2019). Production Options for Psilocybin: Making of the Magic. Chemistry 25, 897–903. doi: 10.1002/chem.201802758

PubMed Abstract | CrossRef Full Text | Google Scholar

Godio, R. P., and Martín, J. F. (2009). Modified oxidosqualene cyclases in the formation of bioactive secondary metabolites: biosynthesis of the antitumor clavaric acid. Fungal Genet. Biol. 46, 232–242. doi: 10.1016/j.fgb.2008.12.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Hage, A., Schoemaker, H. E., and Field, J. A. (1999). Reduction of aryl acids by white-rot fungi for the biocatalytic production of aryl aldehydes and alcohols. Appl. Microbiol. Biotechnol. 52, 834–838. doi: 10.1007/s002530051600

PubMed Abstract | CrossRef Full Text | Google Scholar

Hautzel, R., and Anke, H. (1990). Screening of basidiomycetes and ascomycetes for plant growth regulating substances. Introduction of the gibberellic acid induced de-novo synthesis of hydrolytic enzymes in embryoless seeds of Triticum aestivum as test system. Z. Naturforsch. C 45, 1093–1098. doi: 10.1515/znc-1990-11-1204

CrossRef Full Text | Google Scholar

Jakubczyk, D., and Dussart, F. (2020). Selected Fungal Natural Products with Antimicrobial Properties. Molecules 25:911. doi: 10.3390/molecules25040911

PubMed Abstract | CrossRef Full Text | Google Scholar

Jeleń, H. H. (2003). Use of solid phase microextraction (SPME) for profiling fungal volatile metabolites. Lett. Appl. Microbiol. 36, 263–267. doi: 10.1046/j.1472-765x.2003.01305.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Jensen, K. A., Evans, K. M., Kirk, T. K., and Hammel, K. E. (1994). Biosynthetic pathway for veratryl alcohol in the ligninolytic fungus Phanerochaete chrysosporium. Appl. Environ. Microbiol. 60, 709–714. doi: 10.1128/aem.60.2.709-714.1994

PubMed Abstract | CrossRef Full Text | Google Scholar

Kallberg, Y., Oppermann, U., and Persson, B. (2010). Classification of the short-chain dehydrogenase/reductase superfamily using hidden Markov models. FEBS J. 277, 2375–2386. doi: 10.1111/j.1742-4658.2010.07656.x

PubMed Abstract | CrossRef Full Text | Google Scholar

Mester, T., Swarts, H. J., i Sole, S. R., De Bont, J. A., and Field, J. A. (1997). Stimulation of aryl metabolite production in the basidiomycete Bjerkandera sp. strain BOS55 with biosynthetic precursors and lignin degradation products. Appl. Environ. Microbiol. 63, 1987–1994. doi: 10.1128/aem.63.5.1987-1994.1997

PubMed Abstract | CrossRef Full Text | Google Scholar

Newman, D. J., and Cragg, G. M. (2016). Natural products as sources of new drugs from 1981 to 2014. J. Nat. Prod. 79, 629–661. doi: 10.1021/acs.jnatprod.5b01055

PubMed Abstract | CrossRef Full Text | Google Scholar

Nofiani, R., de Mattos-Shipley, K., Lebe, K. E., Han, L. C., Iqbal, Z., Bailey, A. M., et al. (2018). Strobilurin biosynthesis in Basidiomycete fungi. Nat. Commun. 9:3940.

Google Scholar

Pfefferle, W., Anke, H., Bross, M., and Steglich, W. (1990). Inhibition of solubilized chitin synthase by chlorinated aromatic compounds isolated from mushroom cultures. Agric. Biol. Chem. 54, 1381–1384. doi: 10.1271/bbb1961.54.1381

CrossRef Full Text | Google Scholar

Quin, M. B., Flynn, C. M., and Schmidt-Dannert, C. (2014). Traversing the fungal terpenome. Nat. Prod. Rep. 31, 1449–1473. doi: 10.1039/c4np00075g

PubMed Abstract | CrossRef Full Text | Google Scholar

Rahi, D. K., and Malik, D. (2016). Diversity of mushrooms and their metabolites of nutraceutical and therapeutic significance. J. Mycol. 2016, 1–18. doi: 10.1155/2016/7654123

CrossRef Full Text | Google Scholar

Robbins, W. J., Kavanagh, F., and Hervey, A. (1947). Production of antibiotic substances by Basidiomycetes. Proc.Natl. Acad. Sci. U. S. A. 33, 171–176. doi: 10.1073/pnas.33.6.171

PubMed Abstract | CrossRef Full Text | Google Scholar

Santos-Sánchez, N. F., Salas-Coronado, R., Hernández-Carlos, B., and Villanueva-Cańongo, C. (2019). “Shikimic Acid Pathway in Biosynthesis of Phenolic Compounds,” in Plant Physiological Aspects of Phenolic Compounds, 2nd Edn, ed. M. Soto-Hernández (London: IntechOpen), 1–28. doi: 10.1016/b978-1-4831-9843-9.50005-0

CrossRef Full Text | Google Scholar

Sayari, M., van der Nest, M. A., Steenkamp, E. T., Soal, N. C., Wilken, P. M., and Wingfield, B. D. (2019). Distribution and Evolution of Nonribosomal Peptide Synthetase Gene Clusters in the Ceratocystidaceae. Genes 10:328. doi: 10.3390/genes10050328

PubMed Abstract | CrossRef Full Text | Google Scholar

Schmidt-Dannert, C. (2016). Biocatalytic portfolio of Basidiomycota. Curr. Opin. Chem. Biol. 31, 40–49. doi: 10.1016/j.cbpa.2016.01.002

PubMed Abstract | CrossRef Full Text | Google Scholar

Shinohara, Y., Kawatani, M., Futamura, Y., Osada, H., and Koyama, Y. (2016). An overproduction of astellolides induced by genetic disruption of chromatin-remodeling factors in Aspergillus oryzae. J. Antibiot. 69, 4–8. doi: 10.1038/ja.2015.73

PubMed Abstract | CrossRef Full Text | Google Scholar

Thines, E., Daußmann, T., Sterner, O., Semar, M., and Anke, H. (1995). Fungal melanin biosynthesis inhibitors: introduction of a test system based on the production of dihydroxynaphthalene (DHN) melanin in agar cultures. Z. Naturforsch. C 50, 813–819. doi: 10.1515/znc-1995-11-1212

CrossRef Full Text | Google Scholar

Wawrzyn, G. T., Quin, M. B., Choudhary, S., López-Gallego, F., and Schmidt-Dannert, C. (2012). Draft genome of Omphalotus olearius provides a predictive framework for sesquiterpenoid natural product biosynthesis in Basidiomycota. Chem. Biol. 19, 772–783. doi: 10.1016/j.chembiol.2012.05.012

PubMed Abstract | CrossRef Full Text | Google Scholar

Wilkins, W. H. (1946). Investigation into the production of bacteriostatic substances by fungi. Preliminary examination of the fifth 100 species, all Basidiomycetes, mostly of the wood-destroying type. Br. J. Exp. Pathol. 27:140.

Google Scholar

Wilkins, W. H., and Harris, G. C. M. (1944). Investigation into the production of bacteriostatic substances by fungi VI. Examination of the larger Basidiomycetes. Ann. Appl. Biol. 31, 261–270. doi: 10.1111/j.1744-7348.1944.tb06737.x

CrossRef Full Text | Google Scholar