Genomics-driven discovery of the pneumocandin biosynthetic gene cluster in the fungus Glarea lozoyensis
- Li Chen†1, 2,
- Qun Yue†1,
- Xinyu Zhang1,
- Meichun Xiang1,
- Chengshu Wang3,
- Shaojie Li1,
- Yongsheng Che4,
- Francisco Javier Ortiz-López5,
- Gerald F Bills5, 6,
- Xingzhong Liu1Email author and
- Zhiqiang An6Email author
© Chen et al.; licensee BioMed Central Ltd. 2013
Received: 19 December 2012
Accepted: 9 May 2013
Published: 20 May 2013
The antifungal therapy caspofungin is a semi-synthetic derivative of pneumocandin B0, a lipohexapeptide produced by the fungus Glarea lozoyensis, and was the first member of the echinocandin class approved for human therapy. The nonribosomal peptide synthetase (NRPS)-polyketide synthases (PKS) gene cluster responsible for pneumocandin biosynthesis from G. lozoyensis has not been elucidated to date. In this study, we report the elucidation of the pneumocandin biosynthetic gene cluster by whole genome sequencing of the G. lozoyensis wild-type strain ATCC 20868.
The pneumocandin biosynthetic gene cluster contains a NRPS (GLNRPS4) and a PKS (GLPKS4) arranged in tandem, two cytochrome P450 monooxygenases, seven other modifying enzymes, and genes for L-homotyrosine biosynthesis, a component of the peptide core. Thus, the pneumocandin biosynthetic gene cluster is significantly more autonomous and organized than that of the recently characterized echinocandin B gene cluster. Disruption mutants of GLNRPS4 and GLPKS4 no longer produced the pneumocandins (A0 and B0), and the Δglnrps4 and Δglpks4 mutants lost antifungal activity against the human pathogenic fungus Candida albicans. In addition to pneumocandins, the G. lozoyensis genome encodes a rich repertoire of natural product-encoding genes including 24 PKSs, six NRPSs, five PKS-NRPS hybrids, two dimethylallyl tryptophan synthases, and 14 terpene synthases.
Characterization of the gene cluster provides a blueprint for engineering new pneumocandin derivatives with improved pharmacological properties. Whole genome estimation of the secondary metabolite-encoding genes from G. lozoyensis provides yet another example of the huge potential for drug discovery from natural products from the fungal kingdom.
Several cases of in vivo caspofungin resistance have been reported for Candida and Aspergillus species caused by mutations that reduce the drug sensitivity of the glucan synthase by several thousand-fold [9–12]. A compensatory cell wall remodeling mechanism elevating the chitin content has been found to be associated with caspofungin resistance in C. albicans[13–15]. Generation of pneumocandin derivatives with more desirable pharmacological properties via medicinal chemistry approaches has proven difficult [16, 17]. Elucidation of the biosynthetic pathway to pneumocandins is the first step in applying pathway manipulation and biocombinatorial chemistry approaches to engineer new derivatives with broader spectra of activity and improved physiochemical characteristics to meet the challenges of broader efficacy and clinical resistance.
Based on the structure of pneumocandin, participation of a nonribosomal peptide synthetase (NRPS) and a polyketide synthase (PKS) are predicted for biosynthesis of the cyclic hexapeptide and the 10,12-dimethylmyristoyl side chain [18, 19], respectively. Previous attempts to clone the NRPS and PKS gene cluster responsible for pneumocandin biosynthesis from G. lozoyensis have been unsuccessful [20, 21]. Whole genome sequencing has proven to be an efficient approach in the identification of gene clusters of fungal secondary metabolites, such as PKSs and NRPSs . A recent genomic sequencing project of a pneumocandin B0-overproducing mutant (ATCC 74030) derived from the wild-type (WT) strain of G. lozoyensis was inconclusive in identifying the pneumocandin biosynthetic cluster due to insufficient genome coverage . In this study, we report the elucidation of the pneumocandin biosynthetic gene cluster by genome sequencing of the G. lozoyensis WT strain ATCC 20868. We also compare gene cluster organization with that of the recently published echinocandin B biosynthetic cluster [8, 24]. In addition, analysis of the G. lozoyensis genome revealed a rich repertoire of secondary metabolite-encoding genes that once again illustrates the huge potential for drug discovery from natural products from the fungal kingdom.
The genome characteristics of G. lozoyensis
General features of the G. lozoyensis genome
Assembly size (Mb)
Scaffold N50 (kb)
G + C content (%)
GC exonic (%)
GC intronic (%)
Repeat rate (%)
Gene density (per Mb)
Exons per gene
Chalcone or stilbene synthase gene
Secondary metabolite gene clusters
G. lozoyensisgenome revealed a rich repertoire of secondary metabolite-encoding genes
To identify the pathways involved in the synthesis of secondary metabolites in G. lozoyensis, we searched the genome for genes encoding key enzymes such as NRPS, PKS, terpene synthase (TS), and dimethylallyl tryptophan synthase (DMATS), which are essential for the biosynthesis of peptides, polyketides, terpenes, and alkaloids, respectively. The following secondary metabolite-encoding genes were dispersed among 49 gene clusters: six NRPSs, 24 PKSs, five polyketide synthase-nonribosomal peptide synthase hybrids (PKS-NRPS hybrids), 14 TSs, two DMATSs, 13 NRPS-like, one PKS-like, and one chalcone/stilbene synthase gene (Table 1). In addition to genes encoding the core enzyme(s), the majority of the 49 secondary metabolism gene clusters in G. lozoyensis contained genes encoding other biosynthesis enzymes, transcription regulators, and transporters. For example, about half of the gene clusters contained a gene encoding a Zn2/Cys6 or a C2H2 and C2HC zinc transcriptional factor that could control the expression of genes within of its own cluster. Also, about 60% of the secondary metabolism clusters contained a gene encoding an ABC or a MFS transporter(s) that could export the metabolites produced by the enzymes encoded by the gene cluster (Additional file 1: Figure S1).
Biosynthetic capabilities of G. lozoyensis
To detect the classes of terpene synthases (TSs) in G. lozoyensis, the homologous sequences were analyzed by using BLAST at NCBI (http://www.ncbi.nlm.nih.gov/) (Additional file 2: Table S1). The richness of TSs, compared to related genome-sequenced fungi , revealed a great potential for G. lozoyensis to produce terpenoids. Three TS genes (GLAREA03340, GLAREA04931, GLAREA10578) encoded geranylgeranyl pyrophosphate synthase and geranylgeranyl transferase, and indicated these genes may be responsible for diterpene and carotenoid biosynthesis . Two genes (GLAREA04679, GLAREA02940) encoding farnesyl pyrophosphate synthetase and farnesyl transferase indicated that sesquiterpenes may be formed . Among these TS genes, only three (GLAREA11903, GLAREA03340, GLAREA08044) were located in gene clusters (Additional file 1: Figure S1). Two DMATS genes were found in the G. lozoyensis genome, and one gene (GLAREA04251) was located in a cluster downstream of another core PKS gene GLPKS9, signifying that a polyketide linked with dimethylallyl tryptophan may be the cluster’s end product.
Identification of GLNRPS4 involving in pneumocandin biosynthesis in G. lozoyensis
The lipohexapeptide pneumocandin consists of two key components: a six-amino acid cyclic peptide and a 10,12-dimethylmyristoyl polyketide side chain . Even though no such products are currently known from functionally characterized PKS-PKS hydrids , it is reasonable to consider that pneumocandins might be encoded by one of the PKS-NRPS hybrid proteins. However, in echinocandin B, the lipid side chain was thought to be derived from the cytoplasmic fatty acid pool . Furthermore, domain analysis precluded the five PKS-NRPS hybrid proteins from pneumocandin biosynthesis because the hybrids contained only one A-T-C module, which could only incorporate one amino acid residue in the polyketide chain (Figure 7). Domain analysis of the six NRPS proteins showed that locus GLAREA10035 contained a NRPS with six A-T-C modules (designated as glnrps4 and boxed) (Figure 7). Therefore, locus GLAREA10035 was the only plausible candidate. GLNRPS4, inferred to be responsible for the biosynthesis of the cyclic-hexapeptide core of the pneumocandins, comprised 7,192 amino acids and was encoded by a gene with two introns (Additional file 1: Figure S3). GLNRPS4 encompassed 20 domains grouped into six modules each corresponding to one of the six amino acid incorporated monomers (Figure 7). The first module of GLNRPS4 had a unique T-C-A-T-C domain structure that differed from the other five modules which contained A-T-C domain structures. Two bioinformatics programs were used for substrate prediction, and both predicted that the third module encoded for proline [58, 59]. However, neither program consistently predicted substrate specificities for the other five modules.
Analysis of the PKS-NRPS gene cluster for pneumocandin biosynthesis
Gene analysis of 50 kb of DNA flanking GLNRPS4 revealed a typical gene cluster for fungal secondary metabolite biosynthesis (Figure 8a). Immediately upstream of GLNRPS4 was the glpks4 gene which encodes a PKS of 2,531 amino acids with eight introns (Additional file 1: Figure S3). Moreover, the PKS encoded by glpks4 contained a methyltransferase domain that would be required for the biosynthesis of methyl group-containing fungal polyketides; the pneumocandin polyketide side chain contains two methyl groups (Figure 1a, Additional file 1: Figure S3) [5, 57]. In addition to GLNRPS4 and GLPKS4, two other genes in this cluster stood out, GLAREA10021 encoding an acyltransferase and GLAREA10043 encoding an acyl-CoA ligase (Figure 8a). Labeling experiments at Merck revealed that GLPKS4 assembled a myristate from an acetyl starter, whereas methionines provided two methyl groups to form the 10,12-dimethylmyristoyl side chain . Although functional characterizations will be necessary to define how each gene contributes to the biosynthetic mechanism, based on the above analyses and those of the echinocandin B and emericellamide pathways [24, 60], a hypothetical model of the pneumocandin biosynthetic pathway can be formulated from the four genes, GLNRPS4, GLPKS4, acyltransferase (GLAREA10021), and acyl-CoA ligase (GLAREA10043). The model predicts that 10,12-dimethylmyristoyl side chain is released from GLPKS4 as a carboxylic acid that is converted to a CoA thioester by the acyl-CoA ligase (GLAREA10043), and then loaded onto the acyltransferase (GLAREA10021). The polyketide intermediate could then be shuttled to the first thiolation (T) domain of GLNRPS4, followed by its acylation to 4,5-dihydroxyorinithine to trigger elongation of the cyclic hexapeptide. Like other fungal NRPS and PKS gene clusters, the glpks4 and glnrps4 are positioned within a cluster that contains genes encoding for one or more cytochrome P450s, clavaminate synthase-like proteins (oxygenases), zinc finger transcription factors, and an ABC transporter (Figure 8a). It has been demonstrated that proline 3-hydroxylase and proline 4-hydroxylase, which are members of the 2-oxoglutarate-dependent dioxygenase class, can convert proline to 3-hydroxyproline and 4-hydroxyproline . Two of the four oxygenases (GLAREA10033, GLAREA10041, GLAREA10042, and GLAREA10044) in the gene cluster were presumed to be involved in proline conversion. Two cytochrome P450 monooxygenases (GLAREA10030 and GLAREA10031) were classified in the CYP 512A family by the P450 database (http://www.cyped.uni-stuttgart.de/) which might be responsible for the hydroxylation of the amino acids. These oxygenases were also presumably involved in an oxidative mechanism for the conversion of leucine to methyl proline . The putative zinc finger transcription regulator (GLAREA10050) belongs to the C2H2 and C2HC zinc finger superfamily which are DNA-binding proteins and transcription factors . Some members of this family are pathway-specific transcription regulators of secondary metabolite biosynthesis, e.g., Rua1 that activates the ustilagic acid biosynthesis gene cluster in Ustilago maydis. Therefore, the zinc finger protein GLAREA10050 most likely regulates transcription of the glpks4 and glnrps4 genes. ABC transporters are ubiquitous membrane proteins with the ability to pump a variety of substrate specificities of endogenous and exogenous toxic compounds [64, 65]. The ABC transporter (GLAREA10036) in the cluster possibly secretes antifungal pneumocandins, thus avoiding of intracellular accumulation and ameliorating the toxicity to the producing cells.
Finally, a putative biosynthetic pathway for L-homotyrosine, the non-proteinogenic amino acid in the pneumocandin peptide core’s fourth position, sits downstream of GLNPRS4 (Figure 8a). This set of five contiguous genes showed significant identity to the L-homotyrosine pathway of E. rugulosa (Figure 8b), although the direction of transcription was inverted in two of the five genes, and consisted of GLAREA10037, an aconitase (62% identity to hytD), GLAREA 10038, an isopropyl malate dehydrogenase (71% identity to htyC), GLAREA10039, a 2-isopropyl malate synthase (63.7% identity to hytA), and GLAREA10040, an aminotransferase (64% identity to hytB), and GLAREA10041, a non-heme dioyxgenase (63.8% identity to hytE). However, unlike the L-homotyrosine pathway of E. rugulosa, the cytochrome P450 oxygenase gene corresponding to hytF, was absent (Figure 8b).
Functional analysis of glpks4 and glnrps4in pneumocandin biosynthesis
Sequenced genomes are yielding substantial evidence for a richness of secondary metabolite pathways among the major kinds of fungi, well beyond that imagined to date, and the number of sequenced genomes is growing exponentially [66, 67]. With the advance of next-generation sequencing technology, genome sequencing is evolving as an essential tool to decipher novel genes and gene clusters involved in biosynthesis of different metabolites in fungi [22, 68]. For example, the biosynthetic pathway of the insecticidal cyclodepsipeptide destruxins was recently elucidated in the insect fungal pathogen Metarhizium robertsii by genomic sequencing . Genomic mining of several Aspergillus spp. has led to the elucidation of biosynthetic pathways of multiple bioactive compounds, including terrequinone A , emericellamide , aspyridones ,pyripyropene A  and echinocandin B .
Genomic sequence analysis showed that G. lozoyensis has the potential to produce a diverse array of natural products. The genome was predicted to encode 49 gene clusters that contribute to its secondary metabolome, significantly higher than that of A. sarcoides, also of the Helotiaceae , and in the same order of magnitude as that of B. cinerea, S. sclerotiorum, and other sequenced Leotiomycetes [27, 73]. Most of the ketosynthase domains of the 24 PKSs and five PKS-NRPS hybrids could be clustered with PKSs that were responsible for the biosynthesis of bioactive polyketides and polyketide-nonribosomal peptide hybrids (Figure 6). However, biosynthetic functions for only two of the 49 secondary metabolite-encoding genes in G. lozoyensis were previously validated (GLPKS1 for melanin and GLPKS2 for 6-methylsalicylic acid) [20, 21]. Many secondary metabolites are fusions of nonribosomal peptides and polyketides, in which a PKSs and NRPSs interface and contribute to the same pathway end product [38, 60, 74]. Because the NRPS portion in each of the five PKS-NRPS hybrids in G. lozoyensis genome contains only one A-T-C module, one amino acid is predicted to be added to the polyketide produced by the PKS portion of the cluster, similar to ApdA in A. nidulans and ATEG00325 in A. terreus, which are involved in the biosynthesis of aspyridones and flavipucine, respectively [71, 75].
Comparing the rich genetic potential for secondary metabolites in the G. lozoyensis genome, only pneumocandins were previously identified from the fungus. In an attempt to find additional chemistries, we identified isolecanoric acid and pseudogyrophoric acid as two new fermentation products of G. lozoyensis (Additional file 1: Figure S2). Therefore, majority of the secondary metabolites in G. lozoyensis remain to be characterized. Despite the advances in the field of microbial secondary metabolite biosynthesis, how the basic biology, ecology, and trophic strategies of microorganisms relate to their secondary metabolite production remains poorly understood. Application of efficient strategies to mine the metabolite-encoding gene clusters in G. lozoyensis and other poorly known fungi, while identifying their corresponding metabolites, presents a challenge and opportunity for natural products discovery.
GLNRPS4 and GLPKS4 are centrally located in the pneumocandin biosynthetic gene cluster, and how they cooperate with other genes in the cluster is still speculative. Even though they are independently transcribed and translated, their transcription is likely to be synchronized or co-regulated. The first module in GLNRPS4 has a unique T-C-A-T-C structure, and the first T domain in the T-C-A-T-C module is suggest to accept thiolated intermediates as found in emericellamide biosynthesis  or adenylated substrates similar to yersiniabactin biosynthesis . Thus, the first T domain in the T-C-A-T-C module of GLNRPS4 could be responsible for accepting the incoming 10,12-dimethylmyristoyl side chain intermediate, whereas the second T domain would accept the 4,5-dihydroxyornithine adenylated by the module’s A domain. Threonine, 4-hydroxyproline, 4,5-dihydroxyhomotyrosine, 3-hydroxyglutamine and 3-hydroxyproline/3-hydroxy-4-methylproline would be sequentially added to the growing chain consistent with the in silico prediction that the A3 in GLNRPS4 is specific to proline [8, 58, 59]. Like many other NRPSs [60, 77, 78], the carboxyl terminal of GLNRPS4 lacks a thioesterase (TE) domain, suggesting that a dedicated TE is not required for pneumocandin cyclization. The last C domain of GLNRPS4 is proposed to be responsible for cyclization by condensation to form the peptide bond between 4,5-dihydroxyornithine and 3-hydroxyproline/3-hydroxy-4-methylproline. This proposal is consistent with the fact that the C domain has a HAEYD motif similar to the active site signature in the terminal C domain of cyclosporine synthetase (HSLYD) which is responsible for cyclization of cyclosporine in Tolypocladium inflatum and siderophore synthase SidC (HSLYD) involved in cyclization of the siderophore ferricrocin in A. nidulans[79, 80]. The proposed biosynthetic sequence also parallels that proposed for echinocandin B [8, 24]. Five of the six amino acids in the cyclic hexapeptide were hydroxylated, and hydroxylations of the two proline residues in pneumocandin B0 were catalyzed by a proline-3-hydoxylase and a proline-4-hydoxylase . The enzyme responsible for hydroxylation of 4-methylproline derived from leucine in pneumocandin A0 may also be a proline 3-hydroxylase as 4-methylproline is an analogue of l-proline .
Other genes downstream of the GLNRPS4 that are likely involved the biosynthesis are the putative acyl-CoA ligase GLAREA10043 which shares 43% identity with EasD which converts polyketide carboxylic acid to a CoA thioester in emericellamide biosynthesis in A. nidulans. The putative acyltransferase GLAREA10021 in the cluster shares more than 65% identity with the cholesterol acyltransferases from Cordyceps militaris. Existence of these two genes suggests that the polyketide intermediate was first synthesized by GLPKS4, and then shuttled to the first T domain of GLNRPS4 mediated by the two enzymes, in a fashion similar to the emericellamide biosynthetic pathway . Surprisingly and unlike the echinocandin B pathway , the putative pathway for the homotyrosine residue of the pneumocandin peptide core also sits downstream, and presumably L-homotyrosine biosynthesis is synchronized with the rest of the pathway.
The pneumocandin and echinocandin B pathways have some striking commonalities, yet obviously differ in their organization. The most obvious similarity is the high degree of identity between ecdA and glnrps4 (60.8% identity over 22.7 kb, 55.2% identity over 7218 aa), and both have the same orientation in transcription and functional modules (TCATCATCATCATCATCATCT). Likewise, the genes of the L-homotyrosine pathway are highly similar, although their physical proximities to the core NRPS differ. Both pathways also contain a number of oxygenases that, in the case of echinocandin B, tailor the multiple hydroxyl or diol groups of the amino acid core, but once again their physical location and order are significantly rearranged. However, the inclusion a PKS for side chain biosynthesis and its proximity for immediate loading onto the first thiolation domain, along with close proximity of the L-homotyrosine gene cluster and a possible zinc finger regulatory protein would likely confer greater metabolic autonomy to the pneumocandin pathway. The remarkable similarity between the echinocandin and pneumocandin pathways and especially the high degree of sequence homology between the amp-binding domains of GLNRPS4 and EcdA raises questions about pathway acquisition through horizontal gene transfer among fungi [82, 83]. However, with only two echinocandin type pathways characterized thus far, speculation on why fungi from evolutionary lineages, Eurotiomycete (E. rugulosa) versus Leotiomycete (G. lozoyensis) that diverged 100s of millions of years ago, would share or converge on such similar molecular scaffolds is still premature. Elucidation of additional echinocandin type pathways in the Eurotiomycete, e.g., aculeacin and mulundocandin, and in the Leotiomycetes, e.g. FR901379 (WF11899A) and cryptocandin would yield evidence to determine a possible echinocandin progenitor and the probable directionality in gene recruitment or losses during the evolution of the echinocandin-pneumocandin gene clusters, as well as the significance of these potent cell wall-modifying metabolites to the fungi that produce them.
Elucidation of the pneumocandin biosynthetic pathway in G. lozoyensis paves the way for designing experimental procedures to enhance the production titer of the pneumocandins or engineering analogues with improved oral availability or broader spectrum of antifungal activities. Deletion of other PKS and NRPS genes could potentially reduce metabolic competition for substrates to GLPKS4 and GLNRPS4 and therefore increase the titers of pneumocandin B0, in a manner similar to the disruption of GLPKS1 melanin gene in G. lozoyensis which doubled pneumocandin production titer . Elimination, inactivation, addition or modification of the specificity of domains to GLPKS4 and GLNRPS4 could result in new pneumocandin derivatives via biocombinatorial chemistry approaches for the discovery and development of improved antifungal therapy.
The Glarea lozoyensis genome was sequenced, completely assembled and thoroughly annotated. The menu of secondary metabolites encoding genes was predicted from the genome, thus providing a greater understanding the complexity of primary and secondary metabolism in fungi from the yet poorly studied Leotiomycetes. The biosynthetic gene cluster responsible for pneumocandin was predicted in silico and identified by core gene glpks4 and glnrps4 knockouts and bioassay experiments. The data from this study will form the basis for a more detailed functional analysis of pneumocandin biosynthetic pathways and enable the identification of other antifungal lipohexapeptide pathways in other fungi, of which both will be essential for increasing pneumocandin production and for generating new pneumocandin and echinocandin derivatives via biocombinatorial chemistry approaches.
Fungal and bacterial strains, vectors, and other reagents
The original pneumocandin producing strain of G. lozoyensis ATCC 20868 was obtained from American Type Culture Collection (ATCC) and was used as the wild-type recipient in Agrobacterium-mediated transformation experiments. The Escherichia coli strain DH5α was used in plasmid manipulations. Agrobacterium tumefaciens AGL-1 was described by Lazo et al.. Plasmid pAg1-H3 was described by Zhang et al. , pEASY-T3 vector was from TransGen Biotech (Beijing, China), and pMD18-T vector was from Takara Biotech (Dalian, China). The pneumocandin B0 standard was from Molcan Corporation (Ontario, Canada). LYCP-5 medium, FGY medium and conditions for G. lozoyensis fermentation were described by Connors et al. . M-100 and IMAS mediums were described by Wang et al. . Potato dextrose agar (PDA) and Sabouraud dextrose agar (SDA) were from Becton Dickinson (Franklin Lakes, New Jersey, USA). E. coli and A. tumefaciens AGL-1 were cultured as described by Zhang et al. . Restriction endonucleases and DNA modifying enzymes were from New England Biolabs (Beverly, Massachusetts, USA).
DNA isolation and sequencing
Genomic DNA of G. lozoyensis ATCC 20868 was extracted as previously described by Zhang et al. . Genomic DNA libraries with 500–800 bp inserts were constructed and sequenced with a Roche 454 GS FLX at the Chinese National Human Genome Center in Shanghai. A library with 3 kb inserts was constructed and sequenced with Illumina Genome Analyzer using the protocols as described for genomic sequencing of Cordyceps militaris. The genome sequences were assembled using Newbler software (ver. 2.3) and SSPACE (http://www.baseclear.com/dna-sequencing/data-analysis/bioinformatics-tools/).
G. lozoyensisgenome annotation, orthology and phylogenomic analyses
The G. lozoyensis genome was annotated with Augustus (http://bioinf.uni-greifswald.de/augustus) by referencing annotated genome of Botrytis cinerea. GeneID and GeneMark-ES were additionally used for open reading frames prediction in G. lozoyensis[86, 87]. Repetitive sequences in the genome were identified by BLAST against the RepeatMasker library  and by de novo repetitive sequence search using RepeatModeler (http://www.repeatmasker.org/RepeatModeler.html). Transfer RNAs (tRNAs) were identified with tRNAscan-SE . Ribosomal RNAs (rRNAs) were predicted by a BLAST search with known rRNA modules from other fungal genomes. Whole genome protein families were classified by InterproScan analysis (http://www.ebi.ac.uk/interpro/) and BLAST against Kyoto Encyclopedia of Genes and Genomes (KEGG) database using KEGG Automatic Annotation Server (KAAS: http://www.genome.jp/kegg/kaas/). Carbohydrate-active enzymes from G. lozoyensis and reference fungi (Additional file 2: Table S2) were classified by local Blastp searching against a library of catalytic and carbohydrate-binding module enzymes . PKS, NRPS, DMATS and related gene clusters were predicted by programs SMURF and anti-SMASH  and by manual annotation.
A total of 878 common orthologous genes were identified using the InParanoid pipeline in the selected fungal genomes (Additional file 2: Table S2) , and aligned with Clustal W (ver. 2.0). A maximum likelihood phylogenomic tree was created using the concatenated amino acid sequences in PAUP* 4.0 (beta 10 Win) with heuristic searches . Characters were treated as unordered and gaps were regarded as missing data. Bootstrap support for internal branches was estimated by analysis of 1,000 pseudo replicates. Reference fungi used to construct the phylogenetic tree were described elsewhere (Additional file 2: Table S2 and ref. ). All internal transcribed spacer (ITS) sequences were aligned with Clustal W (ver. 2.0), and a neighbor-joining phylogenetic tree was generated with the program PAUP* 4.0 (beta 10 Win) using 1,000 bootstrap replicates and a Jukes-Cantor substitution model with pairwise deletion for gaps or missing data .
Phylogenetic analysis of PKS and PKS-NRPS genes
KS domains from fungi with PKS genes proven to be responsible for metabolites biosynthesis (Additional file 2: Table S3), PKSs and PKS-NRPS hybrids in G. lozoyensis (Additional file 2: Table S4) were identified by the program anti-SMASH  or visually in multiple alignments. All KS domains from PKS were aligned with Clustal X (ver. 2.0), and analyzed phylogenetically with MEGA 5.0 using a Jones-Taylor-Thornton substitution model, a pair-wise deletion for gaps or missing data, and a 1,000 bootstrap replications test . The tree was rooted with the KS domain of the rat fatty acid synthase (Figure 6).
Gene knockout of glnrps4 and glpks4
To verify function of the predicted pneumocandin gene cluster, glnrps4 and glpks4 deletion were conducted using method reported by Zhang et al. . To verify function of the predicted pneumocandin gene cluster, gene knockout constructs for glnrps4 and glpks4 were created. Briefly, the flanking regions of the target genes were amplified using different primer pairs (Additional file 2: Table S5) and ligated into the binary vector of pAg1-H3 containing the hygromycin resistance gene to form pAg1-H3-nrps4 and pAg1-H3-pks4. The constructs were introduced into G. lozoyensis by Agrobacterium-mediated transformation using method reported by Zhang et al.  with slight modification. Conidia for transformation were harvested into sterile 0.05% Tween-20 followed with 2 times of wash with distilled water and then suspended into 0.5-1.0 mL sterile water (106 spores mL-1). One hundred microliters of G. lozoyensis and 100 μL of A. tumefaciens (OD660 nm = 0.6-0.8) were mixed and spread on the IMAS agar plate and co-incubated at 28°C for 2 d. The co-culture of A.tumefaciens and G. lozoyensis was covered with M-100 supplemented with 300 μg mL-1 cefotaxime and 100 μg mL-1 hygromycin B, and incubated at 25°C for 2–3 weeks before isolating hygB resistant colonies. The transformants were purified by single conidium isolation and the gene knockout transformants were verified by PCR using multiple primers (Additional file 2: Table S5).
HPLC-MS analysis of pneumocandins
Fermentation and pneumocandin extraction protocols were described by Petersen et al. . HPLC separation was performed on an Agilent Zorbax Extend-C18 1.8 μm 2.1 × 50 mm column using an Agilent 1200 Series system (Agilent, USA). The total flow rate was 0.3 mL min-1; mobile phase A was with 0.1% formic acid and mobile phase B was acetonitrile. The total elution program was 25 min. Gradient elution began with 30% B for 0.5 min, changed to 70% B over 3.5 min, changed to 100% B over 8 min, maintained at 100% B for 5 min, to 30% B over 0.5 min, and re-stabilized for 7.5 min prior the next injection. The column temperature was maintained at 40°C. The injection volume was 10 μL.
Mass spectra were acquired with an Agilent Accurate-Mass Quadrupole–Time-of-Flight mass spectrometry (Q-TOF/MS) 6520 system in the positive ionization mode. For Q-TOF/MS conditions, fragmentor and capillary voltages were kept at 130 and 3,500 V, respectively. Nitrogen was supplied as the nebulizing and drying gas. Temperature of the drying gas was set at 30°C. The flow rate of the drying gas and the pressure of the nebulizer were 10 L min-1 and 25 psi, respectively. Full-scan spectra were acquired over a scan range of m/z 80–1,200 at 1.03 spectra s-1.
Candida albicanszone of inhibition (ZOI) assays
Antifungal activity of the WT, glnrps4 and glpks4 gene deletion mutants of G. lozoyensis was measured by a zone of inhibition assay against the human fungal pathogen Candida albicans SC 5314. Ten-mL liquid culture of the wild-type or glnrps4 and glpks4 gene deletion mutants of G. lozoyensis were lyophilized in a vacuum freeze dryer, and 10 mL methanol were added and thoroughly mixed. After 1 h of orbital shaking, the mixtures were first centrifuged at low speed, the supernatant was transferred to glass tubes, and then DMSO (2 mL) was added to solubilize any metabolites precipitated during evaporation. The samples were concentrated to 2 mL under a warm N2 stream during orbital shaking. The final samples were 5× whole broth equivalents including 100% DMSO relative to original culture volume.
Candida albicans SC 5314 cells grown on SDA plates were inoculated into 10 mL of Sabouraud dextrose broth and incubated overnight at 30°C. The C. albicans suspension was adjusted to an optical density of 0.4 at 660 nm and added to SDA in the proportion of 30 mL L-1. Twenty-mL aliquots of the seeded agar media were poured into 9-cm Petri plates. Pneumocandin B0 (5 mg mL-1) and 100% DMSO were used as positive and negative controls. The extracts prepared from liquid culture of G. lozoyensis and the controls (10 μL) were applied to paper discs on the surface of the seeded assay plates. The plates were incubated at 30°C for approximately 20 h and ZOIs were measured and photographed.
Production, purification and identification of isolecanoric and pseudogyrophoric acids
Isolecanoric acid and the new compound pseudogyrophoric acid were isolated from the extract of G. lozoyensis ATCC 20868 grown in MV8 medium (V8 juice 200 mL, maltose 75 g, soy flour 1 g, L-proline 3 g, MES 16.2 g, distilled H2O 800 mL) at 22°C on a rotary shaker at 220 rpm for 14 d. The isolation procedure, mass spectra are summarized in Additional file 1: Figure S2.
Non-ribosomal peptide synthetase
Carbohydrate active enzymes
Carbohydrate binding modules
Dimethylallyl tryptophane synthase
- PKS-NRPS hybrids:
Polyketide synthase-nonribosomal peptide synthetase hybrids
Highly reduced PKS
Non reduced PKS
Internal transcribed spacer
Zone of inhibition.
We thank Drs. Jan Tkacz, Prakash Masurekar, Neal C. Connors and one of the anonymous reviewers for their critical comments on the manuscript; Dr. Wenzhao Wang for HPLC-MS technical support; and Prof. Scott A. Strobel for sharing the Ascocoryne sarcoides genome data. This work was supported by the National Natural Scientific Foundation of China Grant 30625001 to X.L. and the Welch Foundation Grant AU00024 to Z.A.
The Whole Genome Shotgun project has been deposited at DDBJ/EMBL/GenBank under the accession ALVE00000000 Glarea lozoyensis ATCC 20868.
- Denning DW: Echinocandin antifungal drugs. Lancet. 2003, 362: 1142-1151. 10.1016/S0140-6736(03)14472-8.View ArticlePubMedGoogle Scholar
- Kurtz MB, Heath IB, Marrinan J, Dreikorn S, Onishi J, Douglas C: Morphological effects of lipopeptides against Aspergillus fumigatus correlate with activities against (1,3)-beta-D-glucan synthase. Antimicrob Agents Chemother. 1994, 38: 1480-1489. 10.1128/AAC.38.7.1480.PubMed CentralView ArticlePubMedGoogle Scholar
- Cross SA, Scott LJ: Micafungin: a review of its use in adults for the treatment of invasive and oesophageal candidiasis, and as prophylaxis against Candida infections. Drugs. 2008, 68: 2225-2255. 10.2165/00003495-200868150-00010.View ArticlePubMedGoogle Scholar
- Raasch RH: Anidulafungin: review of a new echinocandin antifungal agent. Exp Rev Anti-infec Ther. 2004, 2: 499-508. 10.1586/14787188.8.131.529.View ArticleGoogle Scholar
- Schwartz RE, Sesin DF, Joshua H, Wilson KE, Kempf AJ, Goklen KA, Kuehner D, Gailliot P, Gleason C, White R: Pneumocandins from Zalerion arboricola. I. Discovery and isolation. J Antibiot. 1992, 45: 1853-1866. 10.7164/antibiotics.45.1853.View ArticlePubMedGoogle Scholar
- Nyfeler R, Keller-Schierlein W: Metabolites of microorganisms. 143. Echinocandin B, a novel polypeptide-antibiotic from Aspergillus nidulans var. echinulatus: isolation and structural components. Helv Chim Acta. 1974, 57: 2459-2477. 10.1002/hlca.19740570818.View ArticlePubMedGoogle Scholar
- Iwamoto T, Fujie A, Sakamoto K, Tsurumi Y, Shigematsu N, Yamashita M, Hashimoto S, Okuhara M, Kohsaka M: WF11899A, B and C, novel antifungal lipopeptides. I. Taxonomy, fermentation, isolation and physico-chemical properties. J Antibiot. 1994, 47: 1084-1091. 10.7164/antibiotics.47.1084.View ArticlePubMedGoogle Scholar
- Jiang W, Cacho RA, Chiou G, Garg NK, Tang Y, Walsh CT: EcdGHK are three tailoring iron oxygenases for amino acid building blocks of the echinocandin scaffold. J Am Chem Soc. 2013, 135: 4457-4466. 10.1021/ja312572v.PubMed CentralView ArticlePubMedGoogle Scholar
- Ramos A, Cuervas-Mons V, Noblejas A, Baños I, Duran P, Marcos R, Sánchez-Turrión V, Jiménez M, Arellano B, Corbacho C: Breakthrough rhinocerebral mucormycosis in a liver transplant patient receiving caspofungin. Transplant Proc. 2009, 41: 1972-1975. 10.1016/j.transproceed.2009.01.077.View ArticlePubMedGoogle Scholar
- Madureira A, Bergeron A, Lacroix C, Robin M, Rocha V, de Latour RP, Ferry C, Devergie A, Lapalu J, Gluckman E: Breakthrough invasive aspergillosis in allogeneic haematopoietic stem cell transplant recipients treated with caspofungin. Int J Antimicrob Agents. 2007, 30: 551-554. 10.1016/j.ijantimicag.2007.07.026.View ArticlePubMedGoogle Scholar
- Perlin DS: Resistance to echinocandin-class antifungal drugs. Drug Resist Updat. 2007, 10: 121-130. 10.1016/j.drup.2007.04.002.PubMed CentralView ArticlePubMedGoogle Scholar
- Thompson GR, Wiederhold NP, Vallor AC, Villareal NC, Lewis JS, Patterson TF: Development of caspofungin resistance following prolonged therapy for invasive candidiasis secondary to Candida glabrata infection. Antimicrob Agents Chemother. 2008, 52: 3783-3785. 10.1128/AAC.00473-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Lee KK, MacCallum DM, Jacobsen MD, Walker LA, Odds FC, Gow NAR, Munro CA: Elevated cell wall chitin in Candida albicans confers echinocandin resistance in vivo. Antimicrob Agents Chemother. 2012, 56: 208-217. 10.1128/AAC.00683-11.PubMed CentralView ArticlePubMedGoogle Scholar
- Plaine A, Walker L, Da Costa G, Mora-Montes HM, McKinnon A, Gow NAR, Gaillardin C, Munro CA, Richard ML: Functional analysis of Candida albicans GPI-anchored proteins: roles in cell wall integrity and caspofungin sensitivity. Fungal Gen Biol. 2008, 45: 1404-1414. 10.1016/j.fgb.2008.08.003.View ArticleGoogle Scholar
- Walker LA, Munro CA, De Bruijn I, Lenardon MD, McKinnon A, Gow NAR: Stimulation of chitin synthesis rescues Candida albicans from echinocandins. PloS Path. 2008, 4: e1000040-10.1371/journal.ppat.1000040.View ArticleGoogle Scholar
- Balkovec JM, Black RM, Hammond ML, Heck JV, Zambias RA, Abruzzo G, Bartizal K, Kropp H, Trainor C, Schwartz RE: Synthesis, stability, and biological evaluation of water-soluble prodrugs of a new echinocandin lipopeptide - Discovery of a potential clinical agent for the treatment of systemic candidiasis and Pneumocystis carinii pneumonia (Pcp). J Med Chem. 1992, 35: 194-198. 10.1021/jm00079a027.View ArticlePubMedGoogle Scholar
- Yao J, Liu H, Zhou T, Chen H, Miao Z, Sheng C, Zhang W: Total synthesis and structure–activity relationships of new echinocandin-like antifungal cyclolipohexapeptides. Eur J Med Chem. 2012, 50: 196-208.View ArticlePubMedGoogle Scholar
- Adefarati A, Giacobbe R, Hensens O, Tkacz J: Biosynthesis of L-671,329, an echinocandin-type antibiotic produced by Zalerion arboricola: origins of some of the unusual amino acids and the dimethylmyristic acid side chain. J Am Chem Soc. 1991, 113: 3542-3545. 10.1021/ja00009a048.View ArticleGoogle Scholar
- Adefarati AA, Hensens OD, Jones E, Tkacz J: Pneumocandins from Zalerion arboricola. V. Glutamic acid-and leucine-derived amino acids in pneumocandin A0 (L-671,329) and distinct origins of the substituted proline residues in pneumocandins A0 and B0. J Antibiot. 1992, 45: 1953-1957. 10.7164/antibiotics.45.1953.View ArticlePubMedGoogle Scholar
- Lu P, Zhang A, Dennis LM, Dahl-Roshak AM, Xia YQ, Arison B, An Z, Tkacz JS: A gene (pks2) encoding a putative 6-methylsalicylic acid synthase from Glarea lozoyensis. Mol Gen Genom. 2005, 273: 207-216. 10.1007/s00438-005-1132-y.View ArticleGoogle Scholar
- Zhang A, Lu P, Dahl-Roshak AM, Paress PS, Kennedy S, Tkacz JS, An Z: Efficient disruption of a polyketide synthase gene (pks1) required for melanin synthesis through Agrobacterium-mediated transformation of Glarea lozoyensis. Mol Gen Genom. 2003, 268: 645-655.Google Scholar
- Brakhage AA: Regulation of fungal secondary metabolism. Nat Rev Microbiol. 2013, 11: 21-32.View ArticlePubMedGoogle Scholar
- Youssar L, Gruning BA, Erxleben A, Gunther S, Huttel W: Genome sequence of the fungus Glarea lozoyensis: the first genome sequence of a species from the Helotiaceae family. Euk Cell. 2012, 11: 250-250. 10.1128/EC.05302-11.View ArticleGoogle Scholar
- Cacho RA, Jiang W, Chooi YH, Walsh CT, Tang Y: Identification and characterization of the echinocandin B biosynthetic gene cluster from Emericella rugulosa NRRL 11440. J Am Chem Soc. 2012, 134: 16781-16790. 10.1021/ja307220z.PubMed CentralView ArticlePubMedGoogle Scholar
- Bills GF, Platas G, Peláez F, Masurekar P: Reclassification of a pneumocandin-producing anamorph, Glarea lozoyensis gen. et sp. nov., previously identified as Zalerion arboricola. Mycol Res. 1999, 103: 179-192. 10.1017/S095375629800687X.View ArticleGoogle Scholar
- Peláez F, Collado J, Platas G, Overy DP, Martín J, Vicente F, Basilio A, De la Cruz M, Tormo JR, Gonźlez del Val A: Phylogeny and intercontinental distribution of the pneumocandin-producing anamorphic fungus Glarea lozoyensis. Mycology. 2011, 2: 1-17.View ArticleGoogle Scholar
- Amselem J, Cuomo CA, van Kan JAL, Viaud M, Benito EP, Couloux A, Coutinho PM, de Vries RP, Dyer PS, Fillinger S: Genomic analysis of the necrotrophic fungal pathogens Sclerotinia sclerotiorum and Botrytis cinerea. PloS Gen. 2011, 7: e1002230-10.1371/journal.pgen.1002230.View ArticleGoogle Scholar
- Gianoulis TA, Griffin MA, Spakowicz DJ, Dunican BF, Sboner A, Sismour AM, Kodira C, Egholm M, Church GM, Gerstein MB: Genomic analysis of the hydrocarbon-producing, cellulolytic, endophytic fungus Ascocoryne sarcoides. PloS Gen. 2012, 8: e1002558-10.1371/journal.pgen.1002558.View ArticleGoogle Scholar
- Fisher PJ, Davey RA, Webster J: Degradation of lignin by aquatic and aero-aquatic hyphomycetes. Trans Brit Mycol Soc. 1983, 80: 166-168. 10.1016/S0007-1536(83)80181-8.View ArticleGoogle Scholar
- Fisher PJ, Sharma PD, Webster J: Cellulolytic ability of aero-aquatic hyphomycetes. Trans Brit Mycol Soc. 1977, 69: 495-496. 10.1016/S0007-1536(77)80089-2.View ArticleGoogle Scholar
- Battaglia E, Benoit I, van den Brink J, Wiebenga A, Coutinho PM, Henrissat B, Vries RP: Carbohydrate-active enzymes from the zygomycete fungus Rhizopus oryzae: a highly specialized approach to carbohydrate degradation depicted at genome level. BMC Genom. 2011, 12: 38-10.1186/1471-2164-12-38.View ArticleGoogle Scholar
- van den Brink J, de Vries RP: Fungal enzyme sets for plant polysaccharide degradation. Appl Microbiol Biotechnol. 2011, 91: 1477-1492. 10.1007/s00253-011-3473-2.PubMed CentralView ArticlePubMedGoogle Scholar
- Ohm RA, de Jong JF, Lugones LG, Aerts A, Kothe E, Stajich JE, de Vries RP, Record E, Levasseur A, Baker SE: Genome sequence of the model mushroom Schizophyllum commune. Nat Biotechnol. 2010, 28: 957-U910. 10.1038/nbt.1643.View ArticlePubMedGoogle Scholar
- Cantarel BL, Coutinho PM, Rancurel C, Bernard T, Lombard V, Henrissat B: The Carbohydrate-Active EnZymes database (CAZy): an expert resource for glycogenomics. Nucleic Acids Res. 2009, 37: D233-D238. 10.1093/nar/gkn663.PubMed CentralView ArticlePubMedGoogle Scholar
- Yin YB, Mao XZ, Yang JC, Chen X, Mao FL, Xu Y: dbCAN: a web resource for automated carbohydrate-active enzyme annotation. Nucleic Acids Res. 2012, 40: W445-W451. 10.1093/nar/gks479.PubMed CentralView ArticlePubMedGoogle Scholar
- Martinez D, Berka RM, Henrissat B, Saloheimo M, Arvas M, Baker SE, Chapman J, Chertkov O, Coutinho PM, Cullen D: Genome sequencing and analysis of the biomass-degrading fungus Trichoderma reesei (syn. Hypocrea jecorina). Nat Biotechnol. 2008, 26: 553-560. 10.1038/nbt1403.View ArticlePubMedGoogle Scholar
- Boettger D, Hertweck C: Molecular diversity sculpted by fungal PKS-NRPS hybrids. Chem Bio Chem. 2013, 14: 28-42. 10.1002/cbic.201200624.View ArticlePubMedGoogle Scholar
- Boettger D, Bergmann H, Kuehn B, Shelest E, Hertweck C: Evolutionary imprint of catalytic domains in fungal PKS-NRPS hybrids. Chem Bio Chem. 2012, 13: 2363-2373. 10.1002/cbic.201200449.View ArticlePubMedGoogle Scholar
- Sims JW, Fillmore JP, Warner DD, Schmidt EW: Equisetin biosynthesis in Fusarium heterosporum. Chem Commun. 2005, 14: 186-188.View ArticleGoogle Scholar
- Hendrickson L, Davis CR, Roach C, Nguyen DK, Aldrich T, McAda PC, Reeves CD: Lovastatin biosynthesis in Aspergillus terreus: characterization of blocked mutants, enzyme activities and a multifunctional polyketide synthase gene. Chem Biol. 1999, 6: 429-439. 10.1016/S1074-5521(99)80061-1.View ArticlePubMedGoogle Scholar
- Cox RJ, Glod F, Hurley D, Lazarus CM, Nicholson TP, Rudd BAM, Simpson TJ, Wilkinson B, Zhang Y: Rapid cloning and expression of a fungal polyketide synthase gene involved in squalestatin biosynthesis. Chem Commun. 2004, 10: 2260-2261.View ArticleGoogle Scholar
- Yang G, Rose MS, Turgeon BG, Yoder O: A polyketide synthase is required for fungal virulence and production of the polyketide T-toxin. Plant Cell Online. 1996, 8: 2139-2150.View ArticleGoogle Scholar
- Baker SE, Kroken S, Inderbitzin P, Asvarak T, Li BY, Shi L, Yoder OC, Turgeon BG: Two polyketide synthase-encoding genes are required for biosynthesis of the polyketide virulence factor, T-toxin, by Cochliobolus heterostrophus. Mol Plant-Microbe Interac. 2006, 19: 139-149. 10.1094/MPMI-19-0139.View ArticleGoogle Scholar
- Proctor RH, Desjardins AE, Plattner RD, Hohn TM: A polyketide synthase gene required for biosynthesis of fumonisin mycotoxins in Gibberella fujikuroi mating population A. Fungal Gen Biol. 1999, 27: 100-112. 10.1006/fgbi.1999.1141.View ArticleGoogle Scholar
- Kasahara K, Miyamoto T, Fujimoto T, Oguri H, Tokiwano T, Oikawa H, Ebizuka Y, Fujii I: Solanapyrone synthase, a possible diels-alderase and iterative type I polyketide synthase encoded in a biosynthetic gene cluster from Alternaria solani. Chem Bio Chem. 2010, 11: 1245-1252. 10.1002/cbic.201000173.View ArticlePubMedGoogle Scholar
- Beck J, Ripka S, Siegner A, Schiltz E, Schweizer E: The multifunctional 6-methylsalicylic acid synthase gene of Penicillium patulum - Its gene structure relative to that of other polyketide synthases. Eur J Biochem. 1990, 192: 487-498. 10.1111/j.1432-1033.1990.tb19252.x.View ArticlePubMedGoogle Scholar
- Fujii I, Ono Y, Tada H, Gomi K, Ebizuka Y, Sankawa U: Cloning of the polyketide synthase gene atX from Aspergillus terreus and its identification as the 6-methylsalicylic acid synthase gene by heterologous expression. Mol Gen Genet. 1996, 253: 1-10. 10.1007/s004380050289.View ArticlePubMedGoogle Scholar
- Fulton T, Ibrahim N, Losada M, Grzegorski D, Tkacz J: A melanin polyketide synthase (PKS) gene from Nodulisporium sp. that shows homology to the pks1 gene of Colletotrichum lagenarium. Mol Gen Genet. 1999, 262: 714-720. 10.1007/s004380051133.View ArticlePubMedGoogle Scholar
- Takano Y, Kubo Y, Shimizu K, Mise K, Okuno T, Furusawa I: Structural analysis of PKS1, a polyketide synthase gene involved in melanin biosynthesis in Colletotrichum lagenarium. Mol Gen Genet. 1995, 249: 162-167. 10.1007/BF00290362.View ArticlePubMedGoogle Scholar
- Chooi YH, Cacho R, Tang Y: Identification of the viridicatumtoxin and griseofulvin gene clusters from Penicillium aethiopicum. Chem Biol. 2010, 17: 483-494. 10.1016/j.chembiol.2010.03.015.PubMed CentralView ArticlePubMedGoogle Scholar
- Yu JH, Leonard TJ: Sterigmatocystin biosynthesis in Aspergillus nidulans requires a novel type I polyketide synthase. J Bacteriol. 1995, 177: 4792-4800.PubMed CentralPubMedGoogle Scholar
- Schroeckh V, Scherlach K, Nutzmann HW, Shelest E, Schmidt-Heck W, Schuemann J, Martin K, Hertweck C, Brakhage AA: Intimate bacterial-fungal interaction triggers biosynthesis of archetypal polyketides in Aspergillus nidulans. Proc Nat Acad Sci USA. 2009, 106: 14558-14563. 10.1073/pnas.0901870106.PubMed CentralView ArticlePubMedGoogle Scholar
- Sanchez JF, Chiang YM, Szewczyk E, Davidson AD, Ahuja M, Oakley CE, Bok JW, Keller N, Oakley BR, Wang CCC: Molecular genetic analysis of the orsellinic acid/F9775 gene cluster of Aspergillus nidulans. Mol Biosyst. 2010, 6: 587-593. 10.1039/b904541d.PubMed CentralView ArticlePubMedGoogle Scholar
- Reeves CD, Hu Z, Reid R, Kealey JT: Genes for the biosynthesis of the fungal polyketides hypothemycin from Hypomyces subiculosus and radicicol from Pochonia chlamydosporia. Appl Environ Microbiol. 2008, 74: 5121-5129. 10.1128/AEM.00478-08.PubMed CentralView ArticlePubMedGoogle Scholar
- Zhou H, Qiao KJ, Gao ZZ, Meehan MJ, Li JWH, Zhao XL, Dorrestein PC, Vederas JC, Tang Y: Enzymatic synthesis of resorcylic acid lactones by cooperation of fungal iterative polyketide synthases involved in hypothemycin biosynthesis. J Am Chem Soc. 2010, 132: 4530-4531. 10.1021/ja100060k.PubMed CentralView ArticlePubMedGoogle Scholar
- Keller NP, Turner G, Bennett JW: Fungal secondary metabolism - From biochemistry to genomics. Nat Rev Microbiol. 2005, 3: 937-947. 10.1038/nrmicro1286.View ArticlePubMedGoogle Scholar
- Masurekar PS, Fountoulakis JM, Hallada TC, Sosa MS, Kaplan L: Pneumocandins from Zalerion arboricola. II. Modification of product spectrum by mutation and medium manipulation. J Antibiot. 1992, 45: 1867-1874. 10.7164/antibiotics.45.1867.View ArticlePubMedGoogle Scholar
- Prieto C, Garcia-Estrada C, Lorenzana D, Martin JF: NRPSsp: non-ribosomal peptide synthase substrate predictor. Bioinformatics. 2012, 28: 426-427. 10.1093/bioinformatics/btr659.View ArticlePubMedGoogle Scholar
- Rottig M, Medema MH, Blin K, Weber T, Rausch C, Kohlbacher O: NRPSpredictor2-a web server for predicting NRPS adenylation domain specificity. Nucleic Acids Res. 2011, 39: W362-W367. 10.1093/nar/gkr323.PubMed CentralView ArticlePubMedGoogle Scholar
- Chiang YM, Szewczyk E, Nayak T, Davidson AD, Sanchez JF, Lo HC, Ho WY, Simityan H, Kuo E, Praseuth A: Molecular genetic mining of the Aspergillus secondary metabolome: Discovery of the emericellamide biosynthetic pathway. Chem Biol. 2008, 15: 527-532. 10.1016/j.chembiol.2008.05.010.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen L, Olewinski R, Salmon P, Connors N: Novel proline hydroxylase activities in the pneumocandin-producing fungus Glarea lozoyensis responsible for the formation of trans 3-and trans 4-hydroxyproline. Appl Microbiol Biotechnol. 2003, 62: 263-267. 10.1007/s00253-003-1264-0.View ArticlePubMedGoogle Scholar
- MacPherson S, Larochelle M, Turcotte B: A fungal family of transcriptional regulators: The zinc cluster proteins. Microbiol Mol Biol Rev. 2006, 70: 583-604. 10.1128/MMBR.00015-06.PubMed CentralView ArticlePubMedGoogle Scholar
- Teichmann B, Liu L, Schink KO, Bölker M: The C2H2 zinc finger transcription factor Rua1 activates the ustilagic acid biosynthesis gene cluster in Ustilago maydis. Appl Environ Microbiol. 2010, 76: 2633-2640. 10.1128/AEM.02211-09.PubMed CentralView ArticlePubMedGoogle Scholar
- Kovalchuk A, Driessen A: Phylogenetic analysis of fungal ABC transporters. BMC Genom. 2010, 11: 177-10.1186/1471-2164-11-177.View ArticleGoogle Scholar
- Linton KJ: Structure and function of ABC transporters. Physiology. 2007, 22: 122-130. 10.1152/physiol.00046.2006.View ArticlePubMedGoogle Scholar
- Grigoriev IV, Cullen D, Goodwin SB, Hibbett D, Jeffries TW, Kubicek CP, Kuske C, Magnuson JK, Martin F, Spatafora JW: Fueling the future with fungal genomics. Mycology. 2011, 2: 192-209.Google Scholar
- An Z, Wang C, Liu X, Bennett JW: China’s fungal genomics initiative: a whitepaper. Mycology. 2010, 1: 1-8.View ArticleGoogle Scholar
- Fedorova ND, Moktali V, Medema MH: Bioinformatics approaches and software for detection of secondary metabolic gene clusters. Meth Mol Biol. 2012, 944: 23-45.Google Scholar
- Wang B, Kang QJ, Lu YZ, Bai LQ, Wang CS: Unveiling the biosynthetic puzzle of destruxins in Metarhizium species. Proc Nat Acad Sci USA. 2012, 109: 1287-1292. 10.1073/pnas.1115983109.PubMed CentralView ArticlePubMedGoogle Scholar
- Bok JW, Hoffmeister D, Maggio-Hall LA, Murillo R, Glasner JD, Keller NP: Genomic mining for Aspergillus natural products. Chem Biol. 2006, 13: 31-37. 10.1016/j.chembiol.2005.10.008.View ArticlePubMedGoogle Scholar
- Bergmann S, Schumann J, Scherlach K, Lange C, Brakhage AA, Hertweck C: Genomics-driven discovery of PKS-NRPS hybrid metabolites from Aspergillus nidulans. Nat Chem Biol. 2007, 3: 213-217. 10.1038/nchembio869.View ArticlePubMedGoogle Scholar
- Itoh T, Tokunaga K, Matsuda Y, Fujii I, Abe I, Ebizuka Y, Kushiro T: Reconstitution of a fungal meroterpenoid biosynthesis reveals the involvement of a novel family of terpene cyclases. Nat Chem. 2010, 2: 858-864. 10.1038/nchem.764.View ArticlePubMedGoogle Scholar
- Gao Q, Jin K, Ying SH, Zhang YJ, Xiao GH, Shang YF, Duan ZB, Hu XA, Xie XQ, Zhou G: Genome sequencing and comparative transcriptomics of the model entomopathogenic fungi Metarhizium anisopliae and M. acridum. PloS Gen. 2011, 7: e1001264-10.1371/journal.pgen.1001264.View ArticleGoogle Scholar
- Fischbach MA, Walsh CT: Assembly-line enzymology for polyketide and nonribosomal peptide antibiotics: logic, machinery, and mechanisms. Chem Rev. 2006, 106: 3468-3496. 10.1021/cr0503097.View ArticlePubMedGoogle Scholar
- Gressler M, Zaehle C, Scherlach K, Hertweck C, Brock M: Multifactorial induction of an orphan PKS-NRPS gene cluster in Aspergillus terreus. Chem Biol. 2011, 8: 198-209.View ArticleGoogle Scholar
- Miller DA, Luo LS, Hillson N, Keating TA, Walsh CT: Yersiniabactin synthetase: A four-protein assembly line producing the nonribosomal peptide/polyketide hybrid siderophore of Yersinia pestis. Chem Biol. 2002, 9: 333-344. 10.1016/S1074-5521(02)00115-1.View ArticlePubMedGoogle Scholar
- Bushley KE, Turgeon BG: Phylogenomics reveals subfamilies of fungal nonribosomal peptide synthetases and their evolutionary relationships. BMC Evol Biol. 2010, 10: 26-10.1186/1471-2148-10-26.PubMed CentralView ArticlePubMedGoogle Scholar
- Marahiel MA, Essen LO: Nonribosomal Peptide Synthetases. Mechanistic and Structural Aspects of Essential Domains. Methods Enzymol. 2009, 458: 337-351.View ArticlePubMedGoogle Scholar
- Weber G, Schorgendorfer K, Schneiderscherzer E, Leitner E: The peptide synthetase catalyzing cyclosporine production in Tolypocladium niveum is encoded by a giant 45.8-kilobase open reading frame. Curr Genet. 1994, 26: 120-125. 10.1007/BF00313798.View ArticlePubMedGoogle Scholar
- Eisendle M, Oberegger H, Zadra I, Haas H: The siderophore system is essential for viability of Aspergillus nidulans: functional analysis of two genes encoding l‒ornithine N 5‒monooxygenase (sidA) and a non‒ribosomal peptide synthetase (sidC). Mol Microbiol. 2003, 49: 359-375. 10.1046/j.1365-2958.2003.03586.x.View ArticlePubMedGoogle Scholar
- Zheng P, Xia Y, Xiao G, Xiong C, Hu X, Zhang S, Zheng H, Huang Y, Zhou Y, Wang S: Genome sequence of the insect pathogenic fungus Cordyceps militaris, a valued traditional Chinese medicine. Genom Biol. 2011, 12: R116-10.1186/gb-2011-12-11-r116.View ArticleGoogle Scholar
- Patron NJ, Waller RF, Cozijnsen AJ, Straney DC, Gardiner DM, Nierman WC, Howlett BJ: Origin and distribution of epipolythiodioxopiperazine (ETP) gene clusters in filamentous ascomycetes. BMC Evol Biol. 2007, 7: 174-10.1186/1471-2148-7-174.PubMed CentralView ArticlePubMedGoogle Scholar
- Fitzpatrick DA: Horizontal gene transfer in fungi. FEMS Microbiol Lett. 2012, 329: 1-8. 10.1111/j.1574-6968.2011.02465.x.View ArticlePubMedGoogle Scholar
- Lazo GR, Stein PA, Ludwig RA: A DNA transformation–competent Arabidopsis genomic library in Agrobacterium. Nat Biotechnol. 1991, 9: 963-967. 10.1038/nbt1091-963.View ArticleGoogle Scholar
- Connors N, Petersen L, Hughes R, Saini K, Olewinski R, Salmon P: Residual fructose and osmolality affect the levels of pneumocandins B0 and co-produced by Glarea lozoyensis. Appl Microbiol Biotechnol. 2000, 54: 814-818. 10.1007/s002530000447.View ArticlePubMedGoogle Scholar
- Parra G, Blanco E, Guigó R: Geneid in Drosophila. Genome Res. 2000, 10: 511-515. 10.1101/gr.10.4.511.PubMed CentralView ArticlePubMedGoogle Scholar
- Ter-Hovhannisyan V, Lomsadze A, Chernoff YO, Borodovsky M: Gene prediction in novel fungal genomes using an ab initio algorithm with unsupervised training. Genome Res. 2008, 18: 1979-1990. 10.1101/gr.081612.108.PubMed CentralView ArticlePubMedGoogle Scholar
- Jurka J, Kapitonov VV, Pavlicek A, Klonowski P, Kohany O, Walichiewicz J: Repbase Update, a database of eukaryotic repetitive elements. Cytogenet Genome Res. 2005, 110: 462-467. 10.1159/000084979.View ArticlePubMedGoogle Scholar
- Lowe TM, Eddy SR: tRNAscan-SE: a program for improved detection of transfer RNA genes in genomic sequence. Nucleic Acids Res. 1997, 25: Q955-Q964.View ArticleGoogle Scholar
- Schomburg D, Schomburg I: Enzyme databases. Meth Mol Biol. 2010, 609: 113-128. 10.1007/978-1-60327-241-4_7.View ArticleGoogle Scholar
- Ostlund G, Schmitt T, Forslund K, Kostler T, Messina DN, Roopra S, Frings O, Sonnhammer ELL: InParanoid 7: new algorithms and tools for eukaryotic orthology analysis. Nucleic Acids Res. 2010, 38: D196-D203. 10.1093/nar/gkp931.PubMed CentralView ArticlePubMedGoogle Scholar
- Swofford D: PAUP* version 4.0. Phylogenetic Analysis Using Parsimony (and Other Methods). 2002, Sunderland, Massachusetts, USA: Sinauer AssociatesGoogle Scholar
- Medema MH, Blin K, Cimermancic P, de Jager V, Zakrzewski P, Fischbach MA, Weber T, Takano E, Breitling R: AntiSMASH: rapid identification, annotation and analysis of secondary metabolite biosynthesis gene clusters in bacterial and fungal genome sequences. Nucleic Acids Res. 2011, 39: W339-W346. 10.1093/nar/gkr466.PubMed CentralView ArticlePubMedGoogle Scholar
- Tamura K, Peterson D, Peterson N, Stecher G, Nei M, Kumar S: MEGA5: molecular evolutionary genetics analysis using maximum likelihood, evolutionary distance, and maximum parsimony methods. Mol Biol Evol. 2011, 28: 2731-2739. 10.1093/molbev/msr121.PubMed CentralView ArticlePubMedGoogle Scholar
- Petersen LA, Hughes DL, Hughes R, DiMichele L, Salmon P, Connors N: Effects of amino acid and trace element supplementation on pneumocandin production by Glarea lozoyensis: impact on titer, analogue levels, and the identification of new analogues of pneumocandin B-0. J Indust Microbiol Biotechnol. 2001, 26: 216-221. 10.1038/sj.jim.7000115.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.