Metabolomic Analysis Demonstrates the Impacts of Polyketide Synthases PKS14 and PKS15 on the Production of Beauvericins, Bassianolide, Enniatin A, and Ferricrocin in Entomopathogen Beauveria bassiana.

Beauveria bassiana is a globally distributed entomopathogenic fungus that produces various secondary metabolites to support its pathogenesis in insects. Two polyketide synthase genes, pks14 and pks15, are highly conserved in entomopathogenic fungi and are important for insect virulence. However, understanding of their mechanisms in insect pathogenicity is still limited. Here, we overexpressed these two genes in B. bassiana and compared the metabolite profiles of pks14 and pks15 overexpression strains to those of their respective knockout strains in culture and in vivo using tandem liquid chromatography-mass spectrometry (LC-MS/MS) with Global Natural Products Social Molecular Networking (GNPS). The pks14 and pks15 clusters exhibited crosstalk with biosynthetic clusters encoding insect-virulent metabolites, including beauvericins, bassianolide, enniatin A, and the intracellular siderophore ferricrocin under certain conditions. These secondary metabolites were upregulated in the pks14-overexpressing strain in culture and the pks15-overexpressing strain in vivo. These data suggest that pks14 and pks15, their proteins or their cluster components might be directly or indirectly associated with key pathways in insect pathogenesis of B. bassiana, particularly those related to secondary metabolism. Information about interactions between the polyketide clusters and other biosynthetic clusters improves scientific understanding about crosstalk among biosynthetic pathways and mechanisms of pathogenesis.

Microbial polyketides and their roles in insect virulence: from genomics to biological functions.

Covering: May 1966 up to January 2022Entomopathogenic microorganisms have potential for biological control of insect pests. Their main secondary metabolites include polyketides, nonribosomal peptides, and polyketide-nonribosomal peptide (PK-NRP) hybrids. Among these secondary metabolites, polyketides have mainly been studied for structural identification, pathway engineering, and for their contributions to medicine. However, little is known about the function of polyketides in insect virulence. This review focuses on the role of bacterial and fungal polyketides, as well as PK-NRP hybrids in insect infection and killing. We also discuss gene distribution and evolutional relationships among different microbial species. Further, the role of microbial polyketides and the hybrids in modulating insect-microbial symbiosis is also explored. Understanding the mechanisms of polyketides in insect pathogenesis, how compounds moderate the host-fungus interaction, and the distribution of PKS genes across different fungi and bacteria will facilitate the discovery and development of novel polyketide-derived bio-insecticides.

Iron homeostasis in the absence of ferricrocin and its consequences in fungal development and insect virulence in Beauveria bassiana.

The putative ferricrocin synthetase gene ferS in the fungal entomopathogen Beauveria bassiana BCC 2660 was identified and characterized. The 14,445-bp ferS encodes a multimodular nonribosomal siderophore synthetase tightly clustered with Fusarium graminearum ferricrocin synthetase. Functional analysis of this gene was performed by disruption with the bar cassette. ΔferS mutants were verified by Southern and PCR analyses. HPLC and TLC analyses of crude extracts indicated that biosynthesis of ferricrocin was abolished in ΔferS. Insect bioassays surprisingly indicated that ΔferS killed the Spodoptera exigua larvae faster (LT50 59 h) than wild type (66 h). Growth and developmental assays of the mutant and wild type demonstrated that ΔferS had a significant increase in germination under iron depletion and radial growth and a decrease in conidiation. Mitotracker staining showed that the mitochondrial activity was enriched in ΔferS under both iron excess and iron depletion. Comparative transcriptomes between wild type and ΔferS indicated that the mutant was increased in the expression of eight cytochrome P450 genes and those in iron homeostasis, ferroptosis, oxidative stress response, ergosterol biosynthesis, and TCA cycle, compared to wild type. Our data suggested that ΔferS sensed the iron excess and the oxidative stress and, in turn, was up-regulated in the antioxidant-related genes and those in ergosterol biosynthesis and TCA cycle. These increased biological pathways help ΔferS grow and germinate faster than the wild type and caused higher insect mortality than the wild type in the early phase of infection.

The Fungus Metarhizium sp. BCC 4849 Is an Effective and Safe Mycoinsecticide for the Management of Spider Mites and Other Insect Pests.

Five isolates of Metarhizium sp. were evaluated for their pathogenicity against the spider mite (Tetranychus truncatus Ehara) (Acari: Tetranychidae) and Metarhizium sp. BCC 4849 resulted in the highest mortality (82%) on the 5th day post-inoculation (DPI). Subsequent insect bioassay data indicated similar high virulence against five other insects: African red mites (Eutetranychus africanus Tucker) (Acari: Tetranychidae), bean aphid (Aphis craccivora Koch) (Hemiptera: Aphididae), cassava mealybug (Phenacoccus manihoti Matile-Ferrero) (Hemiptera: Pseudococcidae), sweet potato weevil (Cylas formicarius Fabricius) (Coleoptera: Brentidae), and oriental fruit fly (Bactrocera dorsalis Hendel) (Diptera: Tephritidae), at mortalities of 92-99%, on 3rd-6th DPI, and in laboratory conditions. The pathogenicity assay against E. africanus in hemp plants under greenhouse conditions indicated 85-100% insect mortality on 10th DPI using the fungus alone or in combination with synthetic acaricide. Genome sequencing of Metarhizium sp. BCC 4849 revealed the high abundance of proteins associated with zinc-, heme-, and iron-binding; oxidation-reduction; and transmembrane transport, implicating its versatile mode of interaction with the environment and adaptation to various ion homeostasis. The light and scanning electron microscopy indicated that at 24 h post inoculation (PI), adhesion and appressorial formation occurred, notably near the setae. Most infected mites had stopped moving and started dying by 48-72 h PI. Elongated hyphal bodies and oval blastospores were detected in the legs. At 96-120 h PI or longer, dense mycelia and conidial mass had colonized the interior and exterior of dead mites, primarily at the bottom than the upper part. The shelf-life study also indicated that conidial formulation combined with an oxygen-moisture absorber markedly enhanced the viability and germination after storage at 35 °C for four months. The fungus was tested as safe for humans and animals, according to our toxicological assays.

The polyketide synthase PKS15 has a crucial role in cell wall formation in Beauveria bassiana.

Entomopathogenic fungi utilize specific secondary metabolites to defend against insect immunity, thereby enabling colonization of their specific hosts. We are particularly interested in the polyketide synthesis gene pks15, which is involved in metabolite production, and its role in fungal virulence. Targeted disruption of pks15 followed by genetic complementation with a functional copy of the gene would allow for functional characterization of this secondary metabolite biosynthesis gene. Using a Beauveria bassiana ∆pks15 mutant previously disrupted by a bialophos-resistance (bar) cassette, we report here an in-cis complementation at bar cassette using CRISPR/Cas9 gene editing. A bar-specific short guide RNA was used to target and cause a double-strand break in bar, and a donor DNA carrying a wild-type copy of pks15 was co-transformed with the guide RNA. Isolate G6 of ∆pks15 complemented with pks15 was obtained and verified by PCR, Southern analyses and DNA sequencing. Compared to ∆pks15 which showed a marked reduction in sporulation and insect virulence, the complementation in G6 restored with insect virulence, sporulation and conidial germination to wild-type levels. Atomic force and scanning electron microscopy revealed that G6 and wild-type conidial wall surfaces possessed the characteristic rodlet bundles and rough surface while ∆pks15 walls lacked the bundles and were relatively smoother. Conidia of ∆pks15 were larger and more elongated than that of G6 and the wild type, indicating changes in their cell wall organization. Our data indicate that PKS15 and its metabolite are likely not only important for fungal virulence and asexual reproduction, but also cell wall formation.

Diverse Microbial Community Profiles of Propionate-Degrading Cultures Derived from Different Sludge Sources of Anaerobic Wastewater Treatment Plants.

Anaerobic digestion (AD) has been used for wastewater treatment and production of renewable energy or biogas. Propionate accumulation is one of the important problems leading to an unstable system and low methane production. Revealing propionate-degrading microbiome is necessary to gain a better knowledge for alleviation of the problem. Herein, we systematically investigated the propionate-degrading cultures enriched from various anaerobic sludge sources of agro-industrial wastewater treatment plants using 16S rRNA gene sequencing. Different microbial profiles were shown even though the methanogenic activities of all cultures were similar. Interestingly, non-classical propionate-degrading key players Smithella, Syntrophomonas, and Methanosaeta were observed as common prevalent taxa in our enriched cultures. Moreover, different hydrogenotrophic methanogens were found specifically to the different sludge sources. The enriched culture of high salinity sludge showed a distinct microbial profile compared to the others, containing mainly Thermovirga, Anaerolinaceae, Methanosaeta, Syntrophobactor, and Methanospirillum. Our microbiome analysis revealed different propionate-degrading community profiles via mainly the Smithella pathway and offers inside information for microbiome manipulation in AD systems to increase biogas production corresponding to their specific microbial communities.

Understanding carbon utilization routes between high and low starch-producing cultivars of cassava through Flux Balance Analysis.

Analysis of metabolic flux was used for system level assessment of carbon partitioning in Kasetsart 50 (KU50) and Hanatee (HN) cassava cultivars to understand the metabolic routes for their distinct phenotypes. First, the constraint-based metabolic model of cassava storage roots, rMeCBM, was developed based on the carbon assimilation pathway of cassava. Following the subcellular compartmentalization and curation to ensure full network connectivity and reflect the complexity of eukaryotic cells, cultivar specific data on sucrose uptake and biomass synthesis were input, and rMeCBM model was used to simulate storage root growth in KU50 and HN. Results showed that rMeCBM-KU50 and rMeCBM-HN models well imitated the storage root growth. The flux-sum analysis revealed that both cultivars utilized different metabolic precursors to produce energy in plastid. More carbon flux was invested in the syntheses of carbohydrates and amino acids in KU50 than in HN. Also, KU50 utilized less flux for respiration and less energy to synthesize one gram of dry storage root. These results may disclose metabolic potential of KU50 underlying its higher storage root and starch yield over HN. Moreover, sensitivity analysis indicated the robustness of rMeCBM model. The knowledge gained might be useful for identifying engineering targets for cassava yield improvement.

The reducing clade IIb polyketide synthase PKS14 acts as a virulence determinant of the entomopathogenic fungus Beauveria bassiana.

The reducing clade IIb polyketide synthase gene, pks14, is preserved throughout the evolution of entomopathogenic fungi. We examined the functions of pks14 in Beauveria bassiana using targeted gene disruption, and pks14 disruption was verified by Southern blot and PCR analyses. The radial growth, cell dry weight and conidial germination of Δpks14 were comparable to that of the wild type. Our sequence and gene expression analyses of the pks14 biosynthetic cluster demonstrated: (i) cotranscription and constitutive expression of nearly all the genes of the aforementioned cluster including the C2H2 zinc finger transcription regulator gene, but not pks14 and the cytochrome P450 gene; (ii) expression of the pks14 gene in the insect-containing culture condition only; and (iii) a KAR9-like gene in direct proximity with pks14 is the only gene showing co-regulation. The Δpks14-infected Spodoptera exigua larvae survived significantly longer than those infected by the wild type, indicating a marked reduction in the virulence of Δpks14 against the insect. LT50 of Δpks14 was increased by 1.55 days. Hyphal body formation was decreased in the hemolymph of insects infected by Δpks14 as compared with those inoculated by the wild type. Our results suggest that PKS14-catalyzed polyketide enhances virulence and pathogenicity of B. bassiana on insects.

Culture degeneration in conidia of Beauveria bassiana and virulence determinants by proteomics.

The quality of Beauveria bassiana conidia directly affects the virulence against insects. In this study, continuous subculturing of B. bassiana on both rice grains and potato dextrose agar (PDA) resulted in 55 and 49 % conidial yield reduction after 12 passages and 68 and 60 % virulence reduction after 20 and 12 passages at four d post-inoculation, respectively. The passage through Tenebrio molitor and Spodoptera exigua restored the virulence of rice and PDA subcultures, respectively. To explore the molecular mechanisms underlying the conidial quality and the decline of virulence after multiple subculturing, we investigated the conidial proteomic changes. Successive subculturing markedly increased the protein levels in oxidative stress response, autophagy, amino acid homeostasis, and apoptosis, but decreased the protein levels in DNA repair, ribosome biogenesis, energy metabolism, and virulence. The nitro blue tetrazolium assay verified that the late subculture's colony and conidia had a higher oxidative stress level than the early subculture. A 2A-type protein phosphatase and a Pleckstrin homology domain protein Slm1, effector proteins of the target of rapamycin (TOR) complex 1 and 2, respectively, were dramatically increased in the late subculture. These results suggest that TOR signalling might be associated with ageing in B. bassiana late subculture, in turn affecting its physiological characteristics and virulence.

Targeted disruption of the polyketide synthase gene pks15 affects virulence against insects and phagocytic survival in the fungus Beauveria bassiana.

The reducing clade III polyketide synthase genes, including pks15, are highly conserved among entomopathogenic fungi. To examine the function of pks15, we used targeted disruption to investigate the impact of Beauveria bassiana pks15 on insect pathogenesis. Southern analysis verified that the Δpks15 mutant was disrupted by a single integration of the transformation cassette at the pks15 locus. The Δpks15 mutant had a slight reduction in radial growth, and it produced fewer spores. Our insect bioassays indicated the Δpks15 mutant to be significantly reduced in virulence against beet armyworms compared to wild type (WT), which could be partially accounted for by its markedly decreased ability to survive phagocytosis. Total haemocyte count decreased sharply by 50-fold from days 1-3 post-inoculation in insects infected with WT, compared to a 5-fold decrease in the Δpks15 mutant. The mutant also produced fewer hemolymph hyphal bodies than WT by 3-fold. In co-culture studies with amoebae that have phagocytic ability similar to that of insect haemocytes, at 48 h the mortality rate of amoebae engulfing Δpks15 decreased by 72 %, and Δpks15 CFU decreased by 83 % compared to co-culture with WT. Thus, the Δpks15 mutant had a reduced ability to cope with phagocytosis and highly reduced virulence in an insect host. These data elucidate a mechanism of insect pathogenesis associated with polyketide biosynthesis.

Overcoming Intrinsic Restriction Enzyme Barriers Enhances Transformation Efficiency in Arthrospira platensis C1.

The development of a reliable genetic transformation system for Arthrospira platensis has been a long-term goal, mainly for those trying either to improve its performance in large-scale cultivation systems or to enhance its value as food and feed additives. However, so far, most of the attempts to develop such a transformation system have had limited success. In this study, an efficient and stable transformation system for A. platensis C1 was successfully developed. Based on electroporation and transposon techniques, exogenous DNA could be transferred to and stably maintained in the A. platensis C1 genome. Most strains of Arthrospira possess strong restriction barriers, hampering the development of a gene transfer system for this group of cyanobacteria. By using a type I restriction inhibitor and liposomes to protect the DNA from nuclease digestion, the transformation efficiency was significantly improved. The transformants were able to grow on a selective medium for more than eight passages, and the transformed DNA could be detected from the stable transformants. We propose that the intrinsic endonuclease enzymes, particularly the type I restriction enzyme, in A. platensis C1 play an important role in the transformation efficiency of this industrial important cyanobacterium.

Production and secretion of naphthoquinones is mediated by the MFS transporter MFS1 in the entomopathogenic fungus Ophiocordyceps sp. BCC1869.

Naphthoquinones are deep red polyketide pigments produced by the ant-pathogenic fungus Ophiocordyceps sp. BCC1869. In culture, biosynthesis of these naphthoquinones remains at a low level during the first 20 days and reaches its maximum production level at approximately 50 days. The MFS transporter gene MFS1 was previously identified in Ophiocordyceps sp. BCC1869 from a subtractive EST library between the fungus grown under naphthoquinone-inductive and naphthoquinone-repressive conditions. We cloned and sequenced this transporter gene, which has an open reading frame of 1505 bp and three introns (48, 52, and 58 bp). Phylogenetic analysis showed this MFS transporter was tightly clustered with fungal riboflavin transporters. Functional analysis of this gene was performed by overexpression of MFS1 under the control of a strong, constitutive promoter. We successfully transformed the fungus with this overexpression plasmid using PEG-protoplast transformation, which generated nine transformants per µg of plasmid. RT-PCR indicated that the MFS1 expression level in the overexpressing strains increased 3- to 10-fold compared to the wild type. HPLC analysis of crude extracts of mutants and wild type demonstrated that four naphthoquinone derivatives, erythrostominone, epierythrostominol, deoxyerythrostominone, and deoxyerythrostominol, were the major naphthoquinones produced and excreted in staggering quantities (20- to 2300-fold) in 7-day old liquid cultures by the mutant C7, compared to the wild type. High resolution electrospray ionization mass spectrometry verified mass spectra of these purified metabolites. Three other naphthoquinone derivatives, whose structures have not been identified, were also detected in high amount in the mutant liquid cultures.

Phylogeny of type I polyketide synthases (PKSs) in fungal entomopathogens and expression analysis of PKS genes in Beauveria bassiana BCC 2660.

Entomopathogenic fungi are able to invade and kill insects. Various secondary metabolites can mediate the interaction of a fungal pathogen with an insect host and also help the fungus compete with other microbes. Here we screened 23 isolates of entomopathogenic fungi for polyketide synthase (PKS) genes and amplified 72 PKS gene fragments using degenerate PCR. We performed a phylogenetic analysis of conserved ketosynthase and acyltransferase regions in these 72 sequences and 72 PKSs identified from four insect fungal genome sequences. The resulting genealogy indicated 47 orthologous groups with 99-100 % bootstrap support, suggesting shared biosynthesis of identical or closely related compounds from different fungi. Three insect-specific groups were identified among the PKSs in reducing clades IIa, IIb, and III, which comprised PKSs from 12, 9, and 30 fungal isolates, respectively. A IIa-IIb pair could be found in seven fungi. Expression analyses revealed that eleven out of twelve PKS genes identified in Beauveria bassiana BCC 2660 were expressed in culture. PKS genes from insect-specific clades IIa and IIb were expressed only in insect-containing medium, while others were expressed only in PDB or in CYB, PDB and SDY. The data suggest the potential production of several polyketides in culture.

Tenellin acts as an iron chelator to prevent iron-generated reactive oxygen species toxicity in the entomopathogenic fungus Beauveria bassiana.

Iron is an essential element for life. However, the iron overload can be toxic. Here, we investigated the significant increase of tenellin and iron-tenellin complex production in ferricrocin-deficient mutants of Beauveria bassiana. Our chemical analysis indicated that the ferricrocin-deficient mutants T1, T3 and T5 nearly abolished ferricrocin production. In turn, these mutants had significant accumulation of iron-tenellin complex in their mycelia at 247-289 mg g(-1) cell dry weight under iron-replete condition. Both tenellin and iron-tenellin complex were not detected in the wild-type under such condition. Mass analysis of the mutants' crude extracts demonstrated that tenellin formed a 3:1 complex with iron in the absence of ferricrocin. The unexpected link between ferricrocin and tenellin biosynthesis in ferricrocin-deficient mutants could be a survival strategy during iron-mediated oxidative stress.

Physiological and transcriptional responses to high temperature in Arthrospira (Spirulina) platensis C1.

Arthrospira (Spirulina) platensis is a well-known commercial cyanobacterium that is used as a food and in feed supplements. In this study, we examined the physiological changes and whole-genome expression in A. platensis C1 exposed to high temperature. We found that photosynthetic activity was significantly decreased after the temperature was shifted from 35°C to 42°C for 2 h. A reduction in biomass production and protein content, concomitant with the accumulation of carbohydrate content, was observed after prolonged exposure to high temperatures for 24 h. Moreover, the results of the expression profiling in response to high temperature at the designated time points (8 h) revealed two distinct phases of the responses. The first was the immediate response phase, in which the transcript levels of genes involved in different mechanisms, including genes for heat shock proteins; genes involved in signal transduction and carbon and nitrogen metabolism; and genes encoding inorganic ion transporters for magnesium, nitrite and nitrate, were either transiently induced or repressed by the high temperature. In the second phase, the long-term response phase, both the induction and repression of the expression of genes with important roles in translation and photosynthesis were observed. Taken together, the results of our physiological and transcriptional studies suggest that dynamic changes in the transcriptional profiles of these thermal-responsive genes might play a role in maintaining cell homeostasis under high temperatures, as reflected in the growth and biochemical composition, particularly the protein and carbohydrate content, of A. platensis C1.

Inferring transcriptional gene regulation network of starch metabolism in Arabidopsis thaliana leaves using graphical Gaussian model.

BACKGROUND: Starch serves as a temporal storage of carbohydrates in plant leaves during day/night cycles. To study transcriptional regulatory modules of this dynamic metabolic process, we conducted gene regulation network analysis based on small-sample inference of graphical Gaussian model (GGM). RESULTS: Time-series significant analysis was applied for Arabidopsis leaf transcriptome data to obtain a set of genes that are highly regulated under a diurnal cycle. A total of 1,480 diurnally regulated genes included 21 starch metabolic enzymes, 6 clock-associated genes, and 106 transcription factors (TF). A starch-clock-TF gene regulation network comprising 117 nodes and 266 edges was constructed by GGM from these 133 significant genes that are potentially related to the diurnal control of starch metabolism. From this network, we found that ß-amylase 3 (b-amy3: At4g17090), which participates in starch degradation in chloroplast, is the most frequently connected gene (a hub gene). The robustness of gene-to-gene regulatory network was further analyzed by TF binding site prediction and by evaluating global co-expression of TFs and target starch metabolic enzymes. As a result, two TFs, indeterminate domain 5 (AtIDD5: At2g02070) and constans-like (COL: At2g21320), were identified as positive regulators of starch synthase 4 (SS4: At4g18240). The inference model of AtIDD5-dependent positive regulation of SS4 gene expression was experimentally supported by decreased SS4 mRNA accumulation in Atidd5 mutant plants during the light period of both short and long day conditions. COL was also shown to positively control SS4 mRNA accumulation. Furthermore, the knockout of AtIDD5 and COL led to deformation of chloroplast and its contained starch granules. This deformity also affected the number of starch granules per chloroplast, which increased significantly in both knockout mutant lines. CONCLUSIONS: In this study, we utilized a systematic approach of microarray analysis to discover the transcriptional regulatory network of starch metabolism in Arabidopsis leaves. With this inference method, the starch regulatory network of Arabidopsis was found to be strongly associated with clock genes and TFs, of which AtIDD5 and COL were evidenced to control SS4 gene expression and starch granule formation in chloroplasts.

Draft genome sequence of Arthrospira platensis C1 (PCC9438).

Arthrospira platensis is a cyanobacterium that is extensively cultivated outdoors on a large commercial scale for consumption as a food for humans and animals. It can be grown in monoculture under highly alkaline conditions, making it attractive for industrial production. Here we describe the complete genome sequence of A. platensis C1 strain and its annotation. The A. platensis C1 genome contains 6,089,210 bp including 6,108 protein-coding genes and 45 RNA genes, and no plasmids. The genome information has been used for further comparative analysis, particularly of metabolic pathways, photosynthetic efficiency and barriers to gene transfer.

Biosynthesis of xyrrolin, a new cytotoxic hybrid polyketide/non-ribosomal peptide pyrroline with anticancer potential, in Xylaria sp. BCC 1067.

A gene from Xylaria sp. BCC 1067, pks3, that encodes a putative 3660-residue hybrid polyketide synthase (PKS)/non-ribosomal peptide synthetase (NRPS) was characterised by targeted gene disruption in combination with comprehensive product identification. Studies of the features of a corresponding mutant, YA3, allowed us to demonstrate that pks3 is responsible for the synthesis of a new pyrroline compound, named xyrrolin, in the wild-type Xylaria sp. BCC 1067. The structure of xyrrolin was established by extensive spectroscopic and spectrometric analyses, including low- and high-resolution MS, IR, (1)H NMR, (13)C NMR, (13)C NMR with Dept135, HMQC 2D NMR, HMBC 2D NMR and COSY 2D NMR. On the basis of the Pks3 domain organisation and the chemical structure of xyrrolin, we proposed that biosynthesis of this compound requires the condensation of a tetraketide and an L-serine unit, followed by Dieckmann or reductive cyclisation and enzymatic removal of ketone residue(s). Bioassays of the pure xyrrolin further displayed cytotoxicity against an oral cavity (KB) cancer cell line.