Global analysis of biosynthetic gene clusters reveals vast potential of secondary metabolite production in Penicillium species

Filamentous fungi produce a wide range of bioactive compounds with important pharmaceutical applications, such as antibiotic penicillins and cholesterol-lowering statins. However, less attention has been paid to fungal secondary metabolites compared to those from bacteria. In this study, we sequenced the genomes of 9 Penicillium species and, together with 15 published genomes, we investigated the secondary metabolism of Penicillium and identified an immense, unexploited potential for producing secondary metabolites by this genus. A total of 1,317 putative biosynthetic gene clusters (BGCs) were identified, and polyketide synthase and non-ribosomal peptide synthetase based BGCs were grouped into gene cluster families and mapped to known pathways. The grouping of BGCs allowed us to study the evolutionary trajectory of pathways based on 6-methylsalicylic acid (6-MSA) synthases. Finally, we cross-referenced the predicted pathways with published data on the production of secondary metabolites and experimentally validated the production of antibiotic yanuthones in Penicillia and identified a previously undescribed compound from the yanuthone pathway. This study is the first genus-wide analysis of the genomic diversity of Penicillia and highlights the potential of these species as a source of new antibiotics and other pharmaceuticals.

[1]  U. Mortensen,et al.  Genetics of Polyketide Metabolism in Aspergillus nidulans , 2012, Metabolites.

[2]  Cole Trapnell,et al.  TopHat2: accurate alignment of transcriptomes in the presence of insertions, deletions and gene fusions , 2013, Genome Biology.

[3]  Marnix H Medema,et al.  Comprehensive curation and analysis of fungal biosynthetic gene clusters of published natural products. , 2016, Fungal genetics and biology : FG & B.

[4]  P. Bork,et al.  ETE 3: Reconstruction, Analysis, and Visualization of Phylogenomic Data , 2016, Molecular biology and evolution.

[5]  Andreas Klitgaard,et al.  Accurate Dereplication of Bioactive Secondary Metabolites from Marine-Derived Fungi by UHPLC-DAD-QTOFMS and a MS/HRMS Library , 2014, Marine drugs.

[6]  Mark Yandell,et al.  MAKER2: an annotation pipeline and genome-database management tool for second-generation genome projects , 2011, BMC Bioinformatics.

[7]  Xin Chen,et al.  dbCAN: a web resource for automated carbohydrate-active enzyme annotation , 2012, Nucleic Acids Res..

[8]  Axel A. Brakhage,et al.  Regulation of fungal secondary metabolism , 2012, Nature Reviews Microbiology.

[9]  T. Gabaldón,et al.  Genome sequencing and secondary metabolism of the postharvest pathogen Penicillium griseofulvum , 2016, BMC Genomics.

[10]  Rolf Apweiler,et al.  InterProScan - an integration platform for the signature-recognition methods in InterPro , 2001, Bioinform..

[11]  Jens Nielsen,et al.  Development of fungal cell factories for the production of secondary metabolites: Linking genomics and metabolism , 2017, Synthetic and systems biotechnology.

[12]  J. Frisvad,et al.  New sections in Penicillium containing novel species producing patulin, pyripyropens or other bioactive compounds , 2016, Persoonia.

[13]  R. Breitling,et al.  Detecting Sequence Homology at the Gene Cluster Level with MultiGeneBlast , 2013, Molecular biology and evolution.

[14]  K. Bruno,et al.  Molecular Genetic Characterization of Terreic Acid Pathway in Aspergillus terreus , 2014, Organic letters.

[15]  J. Frisvad,et al.  Reclassification of the Penicillium roqueforti group into three species on the basis of molecular genetic and biochemical profiles. , 1996, Microbiology.

[16]  Steven J. M. Jones,et al.  Abyss: a Parallel Assembler for Short Read Sequence Data Material Supplemental Open Access , 2022 .

[17]  U. Mortensen,et al.  Characterization of four new antifungal yanuthones from Aspergillus niger , 2014, The Journal of Antibiotics.

[18]  U. Mortensen,et al.  Molecular and chemical characterization of the biosynthesis of the 6-MSA-derived meroterpenoid yanuthone D in Aspergillus niger. , 2014, Chemistry & biology.

[19]  I. Oswald,et al.  Molecular cloning and functional characterization of two CYP619 cytochrome P450s involved in biosynthesis of patulin in Aspergillus clavatus. , 2009, Microbiology.

[20]  Bartek Wilczynski,et al.  Biopython: freely available Python tools for computational molecular biology and bioinformatics , 2009, Bioinform..

[21]  H. Laatsch,et al.  Sorbicillin analogues and related dimeric compounds from Penicillium notatum. , 2005, Journal of natural products.

[22]  Carla S. Jones,et al.  Minimum Information about a Biosynthetic Gene cluster. , 2015, Nature chemical biology.

[23]  Alexandros Stamatakis,et al.  RAxML-VI-HPC: maximum likelihood-based phylogenetic analyses with thousands of taxa and mixed models , 2006, Bioinform..

[24]  J. Barrios-González,et al.  Biotechnological production and applications of statins , 2009, Applied Microbiology and Biotechnology.

[25]  K. Katoh,et al.  MAFFT Multiple Sequence Alignment Software Version 7: Improvements in Performance and Usability , 2013, Molecular biology and evolution.

[26]  Intawat Nookaew,et al.  The RAVEN Toolbox and Its Use for Generating a Genome-scale Metabolic Model for Penicillium chrysogenum , 2013, PLoS Comput. Biol..

[27]  Y. Lippi,et al.  Sequencing, physical organization and kinetic expression of the patulin biosynthetic gene cluster from Penicillium expansum. , 2014, International journal of food microbiology.

[28]  Michael A Fischbach,et al.  Computational approaches to natural product discovery. , 2015, Nature chemical biology.

[29]  M. Andersen,et al.  Combining Stable Isotope Labeling and Molecular Networking for Biosynthetic Pathway Characterization. , 2015, Analytical chemistry.

[30]  Yu Shen,et al.  Status and prospect of lignocellulosic bioethanol production in China. , 2010, Bioresource technology.

[31]  C. Stoeckert,et al.  OrthoMCL: identification of ortholog groups for eukaryotic genomes. , 2003, Genome research.

[32]  Kai Blin,et al.  antiSMASH 3.0—a comprehensive resource for the genome mining of biosynthetic gene clusters , 2015, Nucleic Acids Res..

[33]  U. Mortensen,et al.  Investigation of a 6‐MSA Synthase Gene Cluster in Aspergillus aculeatus Reveals 6‐MSA‐derived Aculinic Acid, Aculins A–B and Epi‐Aculin A , 2015, Chembiochem : a European journal of chemical biology.

[34]  Andrej Sali,et al.  A Systematic Computational Analysis of Biosynthetic Gene Cluster Evolution: Lessons for Engineering Biosynthesis , 2014, PLoS Comput. Biol..

[35]  Rustam I. Aminov,et al.  A Brief History of the Antibiotic Era: Lessons Learned and Challenges for the Future , 2010, Front. Microbio..

[36]  Tilmann Weber,et al.  Phylogenetic analysis of condensation domains in NRPS sheds light on their functional evolution , 2007, BMC Evolutionary Biology.

[37]  Ralph A. Cacho,et al.  Identification of the viridicatumtoxin and griseofulvin gene clusters from Penicillium aethiopicum. , 2010, Chemistry & biology.

[38]  T. Gabaldón,et al.  Genome, Transcriptome, and Functional Analyses of Penicillium expansum Provide New Insights Into Secondary Metabolism and Pathogenicity. , 2015, Molecular plant-microbe interactions : MPMI.

[39]  Jens Christian Frisvad,et al.  Dereplication of microbial natural products by LC-DAD-TOFMS. , 2011, Journal of natural products.

[40]  M. Borodovsky,et al.  Integration of mapped RNA-Seq reads into automatic training of eukaryotic gene finding algorithm , 2014, Nucleic acids research.

[41]  Robert C. Edgar,et al.  MUSCLE: multiple sequence alignment with high accuracy and high throughput. , 2004, Nucleic acids research.

[42]  J. Frisvad,et al.  Mycotoxins, drugs and other extrolites produced by species in Penicillium subgenus Penicillium , 2004 .

[43]  S. Kroken,et al.  Phylogenomic analysis of type I polyketide synthase genes in pathogenic and saprobic ascomycetes , 2003, Proceedings of the National Academy of Sciences of the United States of America.

[44]  Myung-soo Park,et al.  A New Record of Penicillium antarcticum from Marine Environments in Korea , 2014, Mycobiology.

[45]  Seung Y. Rhee,et al.  Genomic Signatures of Specialized Metabolism in Plants , 2014, Science.

[46]  Thomas Wetter,et al.  Genome Sequence Assembly Using Trace Signals and Additional Sequence Information , 1999, German Conference on Bioinformatics.

[47]  C. Walsh,et al.  Fungal indole alkaloid biosynthesis: genetic and biochemical investigation of the tryptoquialanine pathway in Penicillium aethiopicum. , 2011, Journal of the American Chemical Society.

[48]  J. Bennett,et al.  Fungal secondary metabolism — from biochemistry to genomics , 2005, Nature Reviews Microbiology.

[49]  Wei Qian,et al.  Selection of conserved blocks from multiple alignments for their use in phylogenetic analysis. , 2000, Molecular biology and evolution.

[50]  O. Gascuel,et al.  New algorithms and methods to estimate maximum-likelihood phylogenies: assessing the performance of PhyML 3.0. , 2010, Systematic biology.

[51]  Sergey I. Nikolenko,et al.  SPAdes: A New Genome Assembly Algorithm and Its Applications to Single-Cell Sequencing , 2012, J. Comput. Biol..

[52]  Cole Trapnell,et al.  Transcript assembly and quantification by RNA-Seq reveals unannotated transcripts and isoform switching during cell differentiation. , 2010, Nature biotechnology.

[53]  R. Simpson,et al.  Soil Microorganisms Mediating Phosphorus Availability Update on Microbial Phosphorus , 2011, Plant Physiology.

[54]  Bud Mishra,et al.  Reevaluating Assembly Evaluations with Feature Response Curves: GAGE and Assemblathons , 2012, PloS one.

[55]  A. Dobson,et al.  Cloning and characterization of novel methylsalicylic acid synthase gene involved in the biosynthesis of isoasperlactone and asperlactone in Aspergillus westerdijkiae. , 2009, Fungal genetics and biology : FG & B.

[56]  B. Le Bizec,et al.  Patulin and secondary metabolite production by marine-derived Penicillium strains. , 2012, Fungal biology.

[57]  María Martín,et al.  UniProt: A hub for protein information , 2015 .

[58]  The Uniprot Consortium,et al.  UniProt: a hub for protein information , 2014, Nucleic Acids Res..

[59]  Ramón Doallo,et al.  ProtTest 3: fast selection of best-fit models of protein evolution , 2011, Bioinform..

[60]  Takayuki Itoh,et al.  Reconstitution of a fungal meroterpenoid biosynthesis reveals the involvement of a novel family of terpene cyclases. , 2010, Nature chemistry.

[61]  J. Nielsen,et al.  Penicillium arizonense, a new, genome sequenced fungal species, reveals a high chemical diversity in secreted metabolites , 2016, Scientific Reports.

[62]  Tilmann Weber,et al.  The evolution of genome mining in microbes - a review. , 2016, Natural product reports.

[63]  J. Varga,et al.  Identification and nomenclature of the genus Penicillium , 2014, Studies in mycology.

[64]  Evgeny M. Zdobnov,et al.  BUSCO: assessing genome assembly and annotation completeness with single-copy orthologs , 2015, Bioinform..

[65]  Jennifer H. Wisecaver,et al.  Fungal metabolic gene clusters—caravans traveling across genomes and environments , 2015, Front. Microbiol..

[66]  Krystle L. Chavarria,et al.  Diversity and evolution of secondary metabolism in the marine actinomycete genus Salinispora , 2014, Proceedings of the National Academy of Sciences.

[67]  C. García-Estrada,et al.  Biosynthetic gene clusters for relevant secondary metabolites produced by Penicillium roqueforti in blue cheeses , 2016, Applied Microbiology and Biotechnology.

[68]  J. Kang,et al.  New polyoxygenated farnesylcyclohexenones, deacetoxyyanuthone A and its hydro derivative from the marine-derived fungus Penicillium sp. , 2003, Journal of natural products.

[69]  P. Galtier,et al.  Biosynthesis and Toxicological Effects of Patulin , 2010, Toxins.

[70]  Alexey A. Gurevich,et al.  QUAST: quality assessment tool for genome assemblies , 2013, Bioinform..

[71]  M S Waterman,et al.  Identification of common molecular subsequences. , 1981, Journal of molecular biology.

[72]  J. Badger,et al.  The Natural Product Domain Seeker NaPDoS: A Phylogeny Based Bioinformatic Tool to Classify Secondary Metabolite Gene Diversity , 2012, PloS one.

[73]  Yan Wu,et al.  Isolation and phosphate‐solubilizing ability of a fungus, Penicillium sp. from soil of an alum mine , 2011, Journal of basic microbiology.

[74]  Jian Wang,et al.  SOAPdenovo2: an empirically improved memory-efficient short-read de novo assembler , 2012, GigaScience.