Expression of Microcystis Biosynthetic Gene Clusters in Natural Populations Suggests Temporally Dynamic Synthesis of Novel and Known Secondary Metabolites in Western Lake Erie.

Microcystis spp. produce diverse secondary metabolites within freshwater cyanobacterial harmful algal blooms (cyanoHABs) around the world. In addition to the biosynthetic gene clusters (BGCs) encoding known compounds, Microcystis genomes harbor numerous BGCs of unknown function, indicating a poorly understood chemical repertoire. While recent studies show that Microcystis produces several metabolites in the lab and field, little work has focused on analyzing the abundance and expression of its broader suite of BGCs during cyanoHAB events. Here, we use metagenomic and metatranscriptomic approaches to track the relative abundance of Microcystis BGCs and their transcripts throughout the 2014 western Lake Erie cyanoHAB. The results indicate the presence of several transcriptionally active BGCs that are predicted to synthesize both known and novel secondary metabolites. The abundance and expression of these BGCs shifted throughout the bloom, with transcript abundance levels correlating with temperature, nitrate, and phosphorus concentrations and the abundance of co-occurring predatory and competitive eukaryotic microorganisms, suggesting the importance of both abiotic and biotic controls in regulating expression. This work highlights the need for understanding the chemical ecology and potential risks to human and environmental health posed by secondary metabolites that are produced but often unmonitored. It also indicates the prospects for identifying pharmaceutical-like molecules from cyanoHAB-derived BGCs. IMPORTANCE Microcystis spp. dominate cyanobacterial harmful algal blooms (cyanoHABs) worldwide and pose significant threats to water quality through the production of secondary metabolites, many of which are toxic. While the toxicity and biochemistry of microcystins and several other compounds have been studied, the broader suite of secondary metabolites produced by Microcystis remains poorly understood, leaving gaps in our understanding of their impacts on human and ecosystem health. We used community DNA and RNA sequences to track the diversity of genes encoding synthesis of secondary metabolites in natural Microcystis populations and assess patterns of transcription in western Lake Erie cyanoHABs. Our results reveal the presence of both known gene clusters that encode toxic secondary metabolites as well as novel ones that may encode cryptic compounds. This research highlights the need for targeted studies of the secondary metabolite diversity in western Lake Erie, a vital freshwater source to the United States and Canada.

Metagenomic and Metatranscriptomic Insights into Population Diversity of Microcystis Blooms: Spatial and Temporal Dynamics of mcy Genotypes, Including a Partial Operon That Can Be Abundant and Expressed.

Cyanobacterial harmful algal blooms (cyanoHABs) degrade freshwater ecosystems globally. Microcystis aeruginosa often dominates cyanoHABs and produces microcystin (MC), a class of hepatotoxins that poses threats to human and animal health. Microcystin toxicity is influenced by distinct structural elements across a diversity of related molecules encoded by variant mcy operons. However, the composition and distribution of mcy operon variants in natural blooms remain poorly understood. Here, we characterized the variant composition of mcy genes in western Lake Erie Microcystis blooms from 2014 and 2018. Sampling was conducted across several spatial and temporal scales, including different bloom phases within 2014, extensive spatial coverage on the same day (2018), and frequent, autonomous sampling over a 2-week period (2018). Mapping of metagenomic and metatranscriptomic sequences to reference sequences revealed three Microcystis mcy genotypes: complete (all genes present [mcyA-J]), partial (truncated mcyA, complete mcyBC, and missing mcyD-J), and absent (no mcy genes). We also detected two different variants of mcyB that may influence the production of microcystin congeners. The relative abundance of these genotypes was correlated with pH and nitrate concentrations. Metatranscriptomic analysis revealed that partial operons were, at times, the most abundant genotype and expressed in situ, suggesting the potential biosynthesis of truncated products. Quantification of genetic divergence between genotypes suggests that the observed strains are the result of preexisting heterogeneity rather than de novo mutation during the sampling period. Overall, our results show that natural Microcystis populations contain several cooccurring mcy genotypes that dynamically shift in abundance spatiotemporally via strain succession and likely influence the observed diversity of the produced congeners. IMPORTANCE Cyanobacteria are responsible for producing microcystins (MCs), a class of potent and structurally diverse toxins, in freshwater systems around the world. While microcystins have been studied for over 50 years, the diversity of their chemical forms and how this variation is encoded at the genetic level remain poorly understood, especially within natural populations of cyanobacterial harmful algal blooms (cyanoHABs). Here, we leverage community DNA and RNA sequences to track shifts in mcy genes responsible for producing microcystin, uncovering the relative abundance, expression, and variation of these genes. We studied this phenomenon in western Lake Erie, which suffers annually from cyanoHAB events, with impacts on drinking water, recreation, tourism, and commercial fishing.

DDAP: docking domain affinity and biosynthetic pathway prediction tool for type I polyketide synthases.

SUMMARY: DDAP is a tool for predicting the biosynthetic pathways of the products of type I modular polyketide synthase (PKS) with the focus on providing a more accurate prediction of the ordering of proteins and substrates in the pathway. In this study, the module docking domain (DD) affinity prediction performance on a hold-out testing dataset reached 0.88 as measured by the area under the receiver operating characteristic (ROC) curve (AUC); the Mean Reciprocal Ranking (MRR) of pathway prediction reached 0.67. DDAP has advantages compared to previous informatics tools in several aspects: (i) it does not rely on large databases, making it a high efficiency tool, (ii) the predicted DD affinity is represented by a probability (0-1), which is more intuitive than raw scores, (iii) its performance is competitive compared to the current popular rule-based algorithm. DDAP is so far the first machine learning based algorithm for type I PKS DD affinity and pathway prediction. We also established the first database of type I modular PKSs, featuring a comprehensive annotation of available docking domains information in bacterial biosynthetic pathways. AVAILABILITY AND IMPLEMENTATION: The DDAP database is available at https://tylii.github.io/ddap. The prediction algorithm DDAP is freely available on GitHub (https://github.com/tylii/ddap) and released under the MIT license. SUPPLEMENTARY INFORMATION: Supplementary data are available at Bioinformatics online.

Structure and Function of NzeB, a Versatile C-C and C-N Bond-Forming Diketopiperazine Dimerase.

The dimeric diketopiperazine (DKPs) alkaloids are a diverse family of natural products (NPs) whose unique structural architectures and biological activities have inspired the development of new synthetic methodologies to access these molecules. However, catalyst-controlled methods that enable the selective formation of constitutional and stereoisomeric dimers from a single monomer are lacking. To resolve this long-standing synthetic challenge, we sought to characterize the biosynthetic enzymes that assemble these NPs for application in biocatalytic syntheses. Genome mining enabled identification of the cytochrome P450, NzeB (Streptomyces sp. NRRL F-5053), which catalyzes both intermolecular carbon-carbon (C-C) and carbon-nitrogen (C-N) bond formation. To identify the molecular basis for the flexible site-selectivity, stereoselectivity, and chemoselectivity of NzeB, we obtained high-resolution crystal structures (1.5 Å) of the protein in complex with native and non-native substrates. This, to our knowledge, represents the first crystal structure of an oxidase catalyzing direct, intermolecular C-H amination. Site-directed mutagenesis was utilized to assess the role individual active-site residues play in guiding selective DKP dimerization. Finally, computational approaches were employed to evaluate plausible mechanisms regarding NzeB function and its ability to catalyze both C-C and C-N bond formation. These results provide a structural and computational rationale for the catalytic versatility of NzeB, as well as new insights into variables that control selectivity of CYP450 diketopiperazine dimerases.

Molecular Basis for Spirocycle Formation in the Paraherquamide Biosynthetic Pathway.

The paraherquamides are potent anthelmintic natural products with complex heptacyclic scaffolds. One key feature of these molecules is the spiro-oxindole moiety that lends a strained three-dimensional architecture to these structures. The flavin monooxygenase PhqK was found to catalyze spirocycle formation through two parallel pathways in the biosynthesis of paraherquamides A and G. Two new paraherquamides (K and L) were isolated from a ΔphqK strain of Penicillium simplicissimum, and subsequent enzymatic reactions with these compounds generated two additional metabolites, paraherquamides M and N. Crystal structures of PhqK in complex with various substrates provided a foundation for mechanistic analyses and computational studies. While it is evident that PhqK can react with various substrates, reaction kinetics and molecular dynamics simulations indicated that the dioxepin-containing paraherquamide L is the favored substrate. Through this effort, we have elucidated a key step in the biosynthesis of the paraherquamides and provided a rationale for the selective spirocyclization of these powerful anthelmintic agents.

Structural basis of the Cope rearrangement and cyclization in hapalindole biogenesis.

Hapalindole alkaloids are a structurally diverse class of cyanobacterial natural products defined by their varied polycyclic ring systems and diverse biological activities. These complex metabolites are generated from a common biosynthetic intermediate by the Stig cyclases in three mechanistic steps: a rare Cope rearrangement, 6-exo-trig cyclization, and electrophilic aromatic substitution. Here we report the structure of HpiC1, a Stig cyclase that catalyzes the formation of 12-epi-hapalindole U in vitro. The 1.5-Å structure revealed a dimeric assembly with two calcium ions per monomer and with the active sites located at the distal ends of the protein dimer. Mutational analysis and computational methods uncovered key residues for an acid-catalyzed [3,3]-sigmatropic rearrangement, as well as specific determinants that control the position of terminal electrophilic aromatic substitution, leading to a switch from hapalindole to fischerindole alkaloids.

Control of Stereoselectivity in Diverse Hapalindole Metabolites is Mediated by Cofactor-Induced Combinatorial Pairing of Stig Cyclases.

Stereospecific polycyclic core formation of hapalindoles and fischerindoles is controlled by Stig cyclases through a three-step cascade involving Cope rearrangement, 6-exo-trig cyclization, and a final electrophilic aromatic substitution. Reported here is a comprehensive study of all currently annotated Stig cyclases, revealing that these proteins can assemble into heteromeric complexes, induced by Ca2+ , to cooperatively control the stereochemistry of hapalindole natural products.

Decoding cyclase-dependent assembly of hapalindole and fischerindole alkaloids.

The formation of C-C bonds in an enantioselective fashion to create complex polycyclic scaffolds in the hapalindole- and fischerindole- type alkaloids from Stigonematales cyanobacteria represents a compelling and urgent challenge in adapting microbial biosynthesis as a catalytic platform in drug development. Here we determine the biochemical basis for tri- and tetracyclic core formation in these secondary metabolites, involving a new class of cyclases that catalyze a complex cyclization cascade.

Structural and stereochemical diversity in prenylated indole alkaloids containing the bicyclo[2.2.2]diazaoctane ring system from marine and terrestrial fungi.

Covering: up to February 2017 Various fungi of the genera Aspergillus, Penicillium, and Malbranchea produce prenylated indole alkaloids possessing a bicyclo[2.2.2]diazaoctane ring system. After the discovery of distinct enantiomers of the natural alkaloids stephacidin A and notoamide B, from A. protuberus MF297-2 and A. amoenus NRRL 35660, another fungi, A. taichungensis, was found to produce their diastereomers, 6-epi-stephacidin A and versicolamide B, as major metabolites. Distinct enantiomers of stephacidin A and 6-epi-stephacidin A may be derived from a common precursor, notoamide S, by enzymes that form a bicyclo[2.2.2]diazaoctane core via a putative intramolecular hetero-Diels-Alder cycloaddition. This review provides our current understanding of the structural and stereochemical homologies and disparities of these alkaloids. Through the deployment of biomimetic syntheses, whole-genome sequencing, and biochemical studies, a unified biogenesis of both the dioxopiperazine and the monooxopiperazine families of prenylated indole alkaloids constituted of bicyclo[2.2.2]diazaoctane ring systems is presented.

Chemoenzymatic Total Synthesis and Structural Diversification of Tylactone-Based Macrolide Antibiotics through Late-Stage Polyketide Assembly, Tailoring, and C-H Functionalization.

Polyketide synthases (PKSs) represent a powerful catalytic platform capable of effecting multiple carbon-carbon bond forming reactions and oxidation state adjustments. We explored the functionality of two terminal PKS modules that produce the 16-membered tylosin macrocycle, using them as biocatalysts in the chemoenzymatic synthesis of tylactone and its subsequent elaboration to complete the first total synthesis of the juvenimicin, M-4365, and rosamicin classes of macrolide antibiotics via late-stage diversification. Synthetic chemistry was employed to generate the tylactone hexaketide chain elongation intermediate that was accepted by the juvenimicin (Juv) ketosynthase of the penultimate JuvEIV PKS module. The hexaketide is processed through two complete modules (JuvEIV and JuvEV) in vitro, which catalyze elongation and functionalization of two ketide units followed by cyclization of the resulting octaketide into tylactone. After macrolactonization, a combination of in vivo glycosylation, selective in vitro cytochrome P450-mediated oxidation, and chemical oxidation was used to complete the scalable construction of a series of macrolide natural products in as few as 15 linear steps (21 total) with an overall yield of 4.6%.

Environmental Pseudomonads Inhibit Cystic Fibrosis Patient-Derived Pseudomonas aeruginosa.

Pseudomonas aeruginosa is an opportunistic pathogen which is evolving resistance to many currently used antibiotics. While much research has been devoted to the roles of pathogenic P. aeruginosa in cystic fibrosis (CF) patients, less is known of its ecological properties. P. aeruginosa dominates the lungs during chronic infection in CF patients, yet its abundance in some environments is less than that of other diverse groups of pseudomonads. Here, we sought to determine if clinical isolates of P. aeruginosa are vulnerable to environmental pseudomonads that dominate soil and water habitats in one-to-one competitions which may provide a source of inhibitory factors. We isolated a total of 330 pseudomonads from diverse habitats of soil and freshwater ecosystems and competed these strains against one another to determine their capacity for antagonistic activity. Over 900 individual inhibitory events were observed. Extending the analysis to P. aeruginosa isolates revealed that clinical isolates, including ones with increased alginate production, were susceptible to competition by multiple environmental strains. We performed transposon mutagenesis on one isolate and identified an ∼14.8-kb locus involved in antagonistic activity. Only two other environmental isolates were observed to carry the locus, suggesting the presence of additional unique compounds or interactions among other isolates involved in outcompeting P. aeruginosa This collection of strains represents a source of compounds that are active against multiple pathogenic strains. With the evolution of resistance of P. aeruginosa to currently used antibiotics, these environmental strains provide opportunities for novel compound discovery against drug-resistant clinical strains. IMPORTANCE: We demonstrate that clinical CF-derived isolates of P. aeruginosa are susceptible to competition in the presence of environmental pseudomonads. We observed that many diverse environmental strains exhibited varied antagonistic profiles against a panel of clinical P. aeruginosa isolates, suggesting the presence of distinct mechanisms of inhibition among these ecological strains. Understanding the properties of these antagonistic events offers the potential for discoveries of antimicrobial compounds or metabolic pathways important to the development of novel treatments for P. aeruginosa infections.

OxaD: A Versatile Indolic Nitrone Synthase from the Marine-Derived Fungus Penicillium oxalicum F30.

Indole alkaloids are a diverse class of natural products known for their wide range of biological activities and complex chemical structures. Rarely observed in this class are indolic nitrones, such as avrainvillamide and waikialoid, which possess potent bioactivities. Herein the oxa gene cluster from the marine-derived fungus Penicillium oxalicum F30 is described along with the characterization of OxaD, a flavin-dependent oxidase that generates roquefortine L, a nitrone-bearing intermediate in the biosynthesis of oxaline. Nitrone functionality in roquefortine L was confirmed by spectroscopic methods and 1,3-dipolar cycloaddition with methyl acrylate. OxaD is a versatile biocatalyst that converts an array of semisynthetic roquefortine C derivatives bearing indoline systems to their respective nitrones. This work describes the first implementation of a nitrone synthase as a biocatalyst and establishes a novel platform for late-stage diversification of a range of complex natural products.

Hapalindole/Ambiguine Biogenesis Is Mediated by a Cope Rearrangement, C-C Bond-Forming Cascade.

Hapalindoles are bioactive indole alkaloids with fascinating polycyclic ring systems whose biosynthetic assembly mechanism has remained unknown since their initial discovery in the 1980s. In this study, we describe the fam gene cluster from the cyanobacterium Fischerella ambigua UTEX 1903 encoding hapalindole and ambiguine biosynthesis along with the characterization of two aromatic prenyltransferases, FamD1 and FamD2, and a previously undescribed cyclase, FamC1. These studies demonstrate that FamD2 and FamC1 act in concert to form the tetracyclic core ring system of the hapalindoles from cis-indole isonitrile and geranyl pyrophosphate through a presumed biosynthetic Cope rearrangement and subsequent 6-exo-trig cyclization/electrophilic aromatic substitution reaction.

Fungal indole alkaloid biogenesis through evolution of a bifunctional reductase/Diels-Alderase.

Prenylated indole alkaloids such as the calmodulin-inhibitory malbrancheamides and anthelmintic paraherquamides possess great structural diversity and pharmaceutical utility. Here, we report complete elucidation of the malbrancheamide biosynthetic pathway accomplished through complementary approaches. These include a biomimetic total synthesis to access the natural alkaloid and biosynthetic intermediates in racemic form and in vitro enzymatic reconstitution to provide access to the natural antipode (+)-malbrancheamide. Reductive cleavage of an L-Pro-L-Trp dipeptide from the MalG non-ribosomal peptide synthetase (NRPS) followed by reverse prenylation and a cascade of post-NRPS reactions culminates in an intramolecular [4+2] hetero-Diels-Alder (IMDA) cyclization to furnish the bicyclo[2.2.2]diazaoctane scaffold. Enzymatic assembly of optically pure (+)-premalbrancheamide involves an unexpected zwitterionic intermediate where MalC catalyses enantioselective cycloaddition as a bifunctional NADPH-dependent reductase/Diels-Alderase. The crystal structures of substrate and product complexes together with site-directed mutagenesis and molecular dynamics simulations demonstrate how MalC and PhqE (its homologue from the paraherquamide pathway) catalyse diastereo- and enantioselective cyclization in the construction of this important class of secondary metabolites.

High-quality draft genome sequence of a biofilm forming lignocellulolytic Aspergillus niger strain ATCC 10864.

Filamentous fungus Aspergillus niger has high industrial value due to their lignocellulolytic enzyme activities and ATCC 10864 is one of the few type strains of A. niger which has a unique biofilm forming capability. Here we report the first draft genome sequence of A. niger ATCC 10864 strain. The genome of A. niger ATCC 10864 is 36,172,237 bp long and comprise of 310 scaffolds with 49.5% average GC content. A total of 10,804 protein-coding genes were predicted among which 10,761 genes were with putative functions. A. niger ATCC 10864 genome coded for 709 putative carbohydrate active enzyme families distributed in six functional categories and among them glycoside hydrolases (GHs) represent the most number of families (279). Genes that include pepA, brlA, exgA, LaeA, rodA, GCN have also been identified in this study, which may play a role in biofilm formation. This high-quality draft genome sequence will facilitate our understanding of the mechanisms behind fungal biofilm formation and higher lignocellulolytic enzyme production.

Insights from the genome of a high alkaline cellulase producing Aspergillus fumigatus strain obtained from Peruvian Amazon rainforest.

Here, we report the complete genome sequence of a high alkaline cellulase producing Aspergillus fumigatus strain LMB-35Aa isolated from soil of Peruvian Amazon rainforest. The genome is ∼27.5mb in size, comprises of 228 scaffolds with an average GC content of 50%, and is predicted to contain a total of 8660 protein-coding genes. Of which, 6156 are with known function; it codes for 607 putative CAZymes families potentially involved in carbohydrate metabolism. Several important cellulose degrading genes, such as endoglucanase A, endoglucanase B, endoglucanase D and beta-glucosidase, are also identified. The genome of A. fumigatus strain LMB-35Aa represents the first whole sequenced genome of non-clinical, high cellulase producing A. fumigatus strain isolated from forest soil.

Discovery of cahuitamycins as biofilm inhibitors derived from a convergent biosynthetic pathway.

Pathogenic microorganisms often have the ability to attach to a surface, building a complex matrix where they colonize to form a biofilm. This cellular superstructure can display increased resistance to antibiotics and cause serious, persistent health problems in humans. Here we describe a high-throughput in vitro screen to identify inhibitors of Acinetobacter baumannii biofilms using a library of natural product extracts derived from marine microbes. Analysis of extracts derived from Streptomyces gandocaensis results in the discovery of three peptidic metabolites (cahuitamycins A-C), with cahuitamycin C being the most effective inhibitor (IC50=14.5 µM). Biosynthesis of cahuitamycin C proceeds via a convergent biosynthetic pathway, with one of the steps apparently being catalysed by an unlinked gene encoding a 6-methylsalicylate synthase. Efforts to assess starter unit diversification through selective mutasynthesis lead to production of unnatural analogues cahuitamycins D and E of increased potency (IC50=8.4 and 10.5 µM).

Developing a genetic system for functional manipulations of FUM1, a polyketide synthase gene for the biosynthesis of fumonisins in Fusarium verticillioides.

Fumonisins are mycotoxins produced by Fusarium verticillioides, a filamentous fungus that is a widespread pathogen of corn. The biosynthesis of fumonisins is catalyzed by an iterative modular polyketide synthase (PKS). The study of the biosynthetic mechanism for these reduced fungal polyketides has been challenging due to the difficulties in detecting the intermediates with a linear carbon chain and manipulating the 7-domain PKS gene from the filamentous fungus. Here, we described the development of a genetic system for functionally manipulating the methyltransferase domain of FUM1 that is responsible for the assembly of a dimethylated 18-carbon chain. Using a two-stage screening strategy, including both positive and negative screenings, we were able to generate mutant strains with a specifically changed active-site in FUM1. LC-MS analyses indicated that biosynthetic intermediates were detectable in the early stage of culture. The results represent the first functional manipulation of the PKS involved in the biosynthesis of fumonisins.

Comparative analysis of the biosynthetic systems for fungal bicyclo[2.2.2]diazaoctane indole alkaloids: the (+)/(-)-notoamide, paraherquamide and malbrancheamide pathways.

The biosynthesis of fungal bicyclo[2.2.2]diazaoctane indole alkaloids with a wide spectrum of biological activities have attracted increasing interest. Their intriguing mode of assembly has long been proposed to feature a non-ribosomal peptide synthetase, a presumed intramolecular Diels-Alderase, a variant number of prenyltransferases, and a series of oxidases responsible for the diverse tailoring modifications of their cyclodipeptide-based structural core. Until recently, the details of these biosynthetic pathways have remained largely unknown due to lack of information on the fungal derived biosynthetic gene clusters. Herein, we report a comparative analysis of four natural product metabolic systems of a select group of bicyclo[2.2.2]diazaoctane indole alkaloids including (+)/(-)-notoamide, paraherquamide and malbrancheamide, in which we propose an enzyme for each step in the biosynthetic pathway based on deep annotation and on-going biochemical studies.

Meta-omic characterization of the marine invertebrate microbial consortium that produces the chemotherapeutic natural product ET-743.

In many macroorganisms, the ultimate source of potent biologically active natural products has remained elusive due to an inability to identify and culture the producing symbiotic microorganisms. As a model system for developing a meta-omic approach to identify and characterize natural product pathways from invertebrate-derived microbial consortia, we chose to investigate the ET-743 (Yondelis) biosynthetic pathway. This molecule is an approved anticancer agent obtained in low abundance (10(-4)-10(-5) % w/w) from the tunicate Ecteinascidia turbinata and is generated in suitable quantities for clinical use by a lengthy semisynthetic process. On the basis of structural similarities to three bacterial secondary metabolites, we hypothesized that ET-743 is the product of a marine bacterial symbiont. Using metagenomic sequencing of total DNA from the tunicate/microbial consortium, we targeted and assembled a 35 kb contig containing 25 genes that comprise the core of the NRPS biosynthetic pathway for this valuable anticancer agent. Rigorous sequence analysis based on codon usage of two large unlinked contigs suggests that Candidatus Endoecteinascidia frumentensis produces the ET-743 metabolite. Subsequent metaproteomic analysis confirmed expression of three key biosynthetic proteins. Moreover, the predicted activity of an enzyme for assembly of the tetrahydroisoquinoline core of ET-743 was verified in vitro. This work provides a foundation for direct production of the drug and new analogues through metabolic engineering. We expect that the interdisciplinary approach described is applicable to diverse host-symbiont systems that generate valuable natural products for drug discovery and development.