cblaster: a remote search tool for rapid identification and visualization of homologous gene clusters

Sustained bacterial N2O reduction at acidic pH

Nitrous oxide (N2O) is a climate-active gas with emissions predicted to increase due to agricultural intensification. Microbial reduction of N2O to dinitrogen (N2) is the major consumption process but microbial N2O reduction under acidic conditions is considered negligible, albeit strongly acidic soils harbor nosZ genes encoding N2O reductase. Here, we study a co-culture derived from acidic tropical forest soil that reduces N2O at pH 4.5. The co-culture exhibits bimodal growth with a Serratia sp. fermenting pyruvate followed by hydrogenotrophic N2O reduction by a Desulfosporosinus sp. Integrated omics and physiological characterization revealed interspecies nutritional interactions, with the pyruvate fermenting Serratia sp. supplying amino acids as essential growth factors to the N2O-reducing Desulfosporosinus sp. Thus, we demonstrate growth-linked N2O reduction between pH 4.5 and 6, highlighting microbial N2O reduction potential in acidic soils.

BGCFlow: systematic pangenome workflow for the analysis of biosynthetic gene clusters across large genomic datasets

Genome mining is revolutionizing natural products discovery efforts. The rapid increase in available genomes demands comprehensive computational platforms to effectively extract biosynthetic knowledge encoded across bacterial pangenomes. Here, we present BGCFlow, a novel systematic workflow integrating analytics for large-scale genome mining of bacterial pangenomes. BGCFlow incorporates several genome analytics and mining tools grouped into five common stages of analysis such as: (i) data selection, (ii) functional annotation, (iii) phylogenetic analysis, (iv) genome mining, and (v) comparative analysis. Furthermore, BGCFlow provides easy configuration of different projects, parallel distribution, scheduled job monitoring, an interactive database to visualize tables, exploratory Jupyter Notebooks, and customized reports. Here, we demonstrate the application of BGCFlow by investigating the phylogenetic distribution of various biosynthetic gene clusters detected across 42 genomes of the Saccharopolyspora genus, known to produce industrially important secondary/specialized metabolites. The BGCFlow-guided analysis predicted more accurate dereplication of BGCs and guided the targeted comparative analysis of selected RiPPs. The scalable, interoperable, adaptable, re-entrant, and reproducible nature of the BGCFlow will provide an effective novel way to extract the biosynthetic knowledge from the ever-growing genomic datasets of biotechnologically relevant bacterial species.

Chalkophomycin Biosynthesis Revealing Unique Enzyme Architecture for a Hybrid Nonribosomal Peptide Synthetase and Polyketide Synthase

Chalkophomycin is a novel chalkophore with antibiotic activities isolated from Streptomyces sp. CB00271, while its potential in studying cellular copper homeostasis makes it an important probe and drug lead. The constellation of N-hydroxylpyrrole, 2H-oxazoline, diazeniumdiolate, and methoxypyrrolinone functional groups into one compact molecular architecture capable of coordinating cupric ions draws interest to unprecedented enzymology responsible for chalkophomycin biosynthesis. To elucidate the biosynthetic machinery for chalkophomycin production, the chm biosynthetic gene cluster from S. sp. CB00271 was identified, and its involvement in chalkophomycin biosynthesis was confirmed by gene replacement. The chm cluster was localized to a ~31 kb DNA region, consisting of 19 open reading frames that encode five nonribosomal peptide synthetases (ChmHIJLO), one modular polyketide synthase (ChmP), six tailoring enzymes (ChmFGMNQR), two regulatory proteins (ChmAB), and four resistance proteins (ChmA'CDE). A model for chalkophomycin biosynthesis is proposed based on functional assignments from sequence analysis and structure modelling, and is further supported by analogy to over 100 chm-type gene clusters in public databases. Our studies thus set the stage to fully investigate chalkophomycin biosynthesis and to engineer chalkophomycin analogues through a synthetic biology approach.

The shikimate pathway: gateway to metabolic diversity

Covering: 1997 to 2023The shikimate pathway is the metabolic process responsible for the biosynthesis of the aromatic amino acids phenylalanine, tyrosine, and tryptophan. Seven metabolic steps convert phosphoenolpyruvate (PEP) and erythrose 4-phosphate (E4P) into shikimate and ultimately chorismate, which serves as the branch point for dedicated aromatic amino acid biosynthesis. Bacteria, fungi, algae, and plants (yet not animals) biosynthesize chorismate and exploit its intermediates in their specialized metabolism. This review highlights the metabolic diversity derived from intermediates of the shikimate pathway along the seven steps from PEP and E4P to chorismate, as well as additional sections on compounds derived from prephenate, anthranilate and the synonymous aminoshikimate pathway. We discuss the genomic basis and biochemical support leading to shikimate-derived antibiotics, lipids, pigments, cofactors, and other metabolites across the tree of life.

Bioinformatics-Guided Reconstitution of Biosynthetic Machineries of Fungal Eremophilane Sesquiterpenes

Eremophilanes exhibit diverse biological activities and chemical structures. This study reports the bioinformatics-guided reconstitution of the biosynthetic machinery of fungal eremophilanes, eremofortin C and sporogen-AO1, to elucidate their biosynthetic pathways. Their biosyntheses include P450-catalyzed multistep oxidation and enzyme-catalyzed isomerization by the DUF3237 family protein. Successful characterization of six P450s enabled us to discuss the functions of eremophilane P450s in putative eremophilane biosynthetic gene clusters, providing opportunities to understand the oxidative modification pathways of fungal eremophilanes.

A Timeline of Biosynthetic Gene Cluster Discovery in Aspergillus fumigatus: From Characterization to Future Perspectives

In 1999, the first biosynthetic gene cluster (BGC), synthesizing the virulence factor DHN melanin, was characterized in Aspergillus fumigatus. Since then, 19 additional BGCs have been linked to specific secondary metabolites (SMs) in this species. Here, we provide a comprehensive timeline of A. fumigatus BGC discovery and find that initial advances centered around the commonly expressed SMs where chemical structure informed rationale identification of the producing BGC (e.g., gliotoxin, fumigaclavine, fumitremorgin, pseurotin A, helvolic acid, fumiquinazoline). Further advances followed the transcriptional profiling of a ΔlaeA mutant, which aided in the identification of endocrocin, fumagillin, hexadehydroastechrome, trypacidin, and fumisoquin BGCs. These SMs and their precursors are the commonly produced metabolites in most A. fumigatus studies. Characterization of other BGC/SM pairs required additional efforts, such as induction treatments, including co-culture with bacteria (fumicycline/neosartoricin, fumigermin) or growth under copper starvation (fumivaline, fumicicolin). Finally, four BGC/SM pairs were discovered via overexpression technologies, including the use of heterologous hosts (fumicycline/neosartoricin, fumihopaside, sphingofungin, and sartorypyrone). Initial analysis of the two most studied A. fumigatus isolates, Af293 and A1160, suggested that both harbored ca. 34-36 BGCs. However, an examination of 264 available genomes of A. fumigatus shows up to 20 additional BGCs, with some strains showing considerable variations in BGC number and composition. These new BGCs present a new frontier in the future of secondary metabolism characterization in this important species.

Noncontiguous operon atlas for the Staphylococcus aureus genome

Bacteria synchronize the expression of genes with related functions by organizing genes into operons so that they are cotranscribed together in a single polycistronic messenger RNA. However, some cellular processes may benefit if the simultaneous production of the operon proteins coincides with the inhibition of the expression of an antagonist gene. To coordinate such situations, bacteria have evolved noncontiguous operons (NcOs), a subtype of operons that contain one or more genes that are transcribed in the opposite direction to the other operon genes. This structure results in overlapping transcripts whose expression is mutually repressed. The presence of NcOs cannot be predicted computationally and their identification requires a detailed knowledge of the bacterial transcriptome. In this study, we used direct RNA sequencing methodology to determine the NcOs map in the Staphylococcus aureus genome. We detected the presence of 18 NcOs in the genome of S. aureus and four in the genome of the lysogenic prophage 80α. The identified NcOs comprise genes involved in energy metabolism, metal acquisition and transport, toxin-antitoxin systems, and control of the phage life cycle. Using the menaquinone operon as a proof of concept, we show that disarrangement of the NcO architecture results in a reduction of bacterial fitness due to an increase in menaquinone levels and a decrease in the rate of oxygen consumption. Our study demonstrates the significance of NcO structures in bacterial physiology and emphasizes the importance of combining operon maps with transcriptomic data to uncover previously unnoticed functional relationships between neighbouring genes.

Evolution-inspired engineering of nonribosomal peptide synthetases

Many clinically used drugs are derived from or inspired by bacterial natural products that often are produced through nonribosomal peptide synthetases (NRPSs), megasynthetases that activate and join individual amino acids in an assembly line fashion. In this work, we describe a detailed phylogenetic analysis of several bacterial NRPSs that led to the identification of yet undescribed recombination sites within the thiolation (T) domain that can be used for NRPS engineering. We then developed an evolution-inspired "eXchange Unit between T domains" (XUT) approach, which allows the assembly of NRPS fragments over a broad range of GC contents, protein similarities, and extender unit specificities, as demonstrated for the specific production of a proteasome inhibitor designed and assembled from five different NRPS fragments.

Biosynthesis of the Azoxy Compound Azodyrecin from Streptomyces mirabilis P8-A2

Azoxy compounds are a distinctive group of bioactive secondary metabolites characterized by a unique RN═N+(O-)R moiety. The azoxy moiety is present in various classes of metabolites that exhibit various biological activities. The enzymatic mechanisms underlying azoxy bond formation remain enigmatic. Azodyrecins are cytotoxic azoxy metabolites produced by Streptomyces mirabilis P8-A2. Here, we cloned and confirmed the putative azd biosynthetic gene cluster through CATCH cloning followed by expression and production of azodyrecins in two heterologous hosts, S. albidoflavus J1074 and S. coelicolor M1146, respectively. We explored the function of 14 enzymes in azodyrecin biosynthesis through gene knockout using CRISPR-Cas9 base editing in the native producer, S. mirabilis P8-A2. The key intermediates were analyzed in the mutants through MS/MS fragmentation studies, revealing azoxy bond formation via the conversion of hydrazine to an azo compound followed by further oxygenation. Enzymes involved in modifications of the precursor could be postulated based on their predicted function and the intermediates identified in the knockout strains. Moreover, the distribution of the azoxy biosynthetic gene clusters across Streptomyces spp. genomes is explored, highlighting the presence of these clusters in over 20% of the Streptomyces spp. genomes and revealing that azoxymycin and valanimycin are scarce, while azodyrecin and KA57A-like clusters are widely distributed across the phylogenetic tree.

Discovery and heterologous biosynthesis of glycosylated polyketide luteodienoside A reveals unprecedented glucinol-mediated product offloading by a fungal carnitine O-acyltransferase domain

Luteodienoside A is a novel glycosylated polyketide produced by the Australian fungus Aspergillus luteorubrus MST-FP2246, consisting of an unusual 1-O-β-d-glucopyranosyl-myo-inositol (glucinol) ester of 3-hydroxy-2,2,4-trimethylocta-4,6-dienoic acid. Mining the genome of A. luteorubrus identified a putative gene cluster for luteodienoside A biosynthesis (ltb), harbouring a highly reducing polyketide synthase (HR-PKS, LtbA) fused at its C-terminus to a carnitine O-acyltransferase (cAT) domain. Heterologous pathway reconstitution in Aspergillus nidulans, substrate feeding assays and gene truncation confirmed the identity of the ltb cluster and demonstrated that the cAT domain is essential for offloading luteodienoside A from the upstream HR-PKS. Unlike previously characterised cAT domains, the LtbA cAT domain uses glucinol as an offloading substrate to release the product from the HR-PKS. Furthermore, the PKS methyltransferase (MT) domain is capable of catalysing gem-dimethylation of the 3-hydroxy-2,2,4-trimethylocta-4,6-dienoic acid intermediate, without requiring reversible product release and recapture by the cAT domain. This study expands the repertoire of polyketide modifications known to be catalysed by cAT domains and highlights the potential of mining fungal genomes for this subclass of fungal PKSs to discover new structurally diverse secondary metabolites.

Global genome mining-driven discovery of an unusual biosynthetic logic for fungal polyketide-terpenoid hybrids

Genome mining has facilitated the efficient discovery of untapped natural products. We performed global genome mining in fungi and discovered a series of biosynthetic gene clusters (BGCs) that appeared to afford polyketide-terpenoid hybrids via a distinct biosynthetic mechanism from those adopted by known pathways. Characterization of one of the BGCs revealed that it yields the drimane-phthalide hybrid 1. During the biosynthesis of 1, the farnesyl group is unusually introduced by the dimethylallyltryptophan synthase-type prenyltransferase MfmD and is then cyclized by the Pyr4-family terpene cyclase MfmH. The replacement of MfmH with its homologue OcdTC gave another hybrid molecule with a monocyclic terpenoid moiety. Moreover, PsetPT, an MfmD homologue, was found to perform dimethylallylation and was then engineered to install a geranyl group. Our study unraveled an unusual biosynthetic mechanism for fungal phthalide-terpenoid hybrids and provided insights into how their structural diversification could be achieved.

Rare NRPS Gene Cluster for Desferriferrichrome Biosynthesis Controls the Conflict between Trap Formation and Nematicidal Activity in Arthrobotrys oligospora

The formation of the trapping device induced by nematodes has been assumed as an indicator for a switch from saprophytic to predacious lifestyles for nematode-trapping fungi. However, fungal nematocidal activity is not completely synonymous with fungal trap formation. We found that the predominant nematode-trapping fungus Arthrobotrys oligospora harbored a rare NRPS (Ao415) gene cluster that was mainly distributed in nematode-trapping fungi. The gene Ao415 putatively encodes a protein with a unique domain organization, distinct from other NRPSs in other fungi. Mutation of the two key biosynthetic genes Ao415 and Ao414 combined with nontarget metabolic analysis revealed that the Ao415 gene cluster was responsible for the biosynthesis of a hydroxamate siderophore, desferriferrichrome (1). Lack of desferriferrichrome (1) and its hydroxamate precursor (3) could lead to significantly increased Fe3+ content, which induced fungal trap formation without a nematode inducer. Furthermore, the addition of Fe3+ strongly improved fungal trap formation but deleteriously caused broken traps. The addition of 1 significantly attenuated trap formation but enhanced fungal nematicidal activity. Our findings indicate that iron is a key factor for trap formation and provide a new insight into the underlying mechanism of siderophores in nematode-trapping fungi.

Production improvement of an antioxidant in cariogenic Streptococcus mutans UA140

AIMS

This study aimed to improve the production of mutantioxidin, an antioxidant encoded by a biosynthetic gene cluster (mao) in Streptococcus mutans UA140, through a series of optimization methods.

METHOD AND RESULTS

Through the construction of mao knockout strain S. mutans UA140∆mao, we identified mutantioxidin as the antioxidant encoded by mao and verified its antioxidant activity through a reactive oxygen species (ROS) tolerance assay. By optimizing the culture medium and fermentation time, 72 h of fermentation in chemically defined medium (CDM) medium was determined as the optimal fermentation conditions. Based on two promoters commonly used in Streptococcus (ldhp and xylS1p), eight promoter refactoring strains were constructed, nevertheless all showed impaired antioxidant production. In-frame deletion and complementation experiments demonstrated the positive regulatory role of mao1 and mao2, on mao. Afterward, the mao1 and mao2, overexpression strain S. mutans UA140/pDL278:: mao1mao2, were constructed, in which the production of mutantioxidin was improved significantly.

CONCLUSIONS

In this study, through a combination of varied strategies such as optimization of fermentation conditions and overexpression of regulatory genes, production of mutantioxidin was increased by 10.5 times ultimately.

Bacterial polyynes uncovered: a journey through their bioactive properties, biosynthetic mechanisms, and sustainable production strategies

Covering: up to 2023Conjugated polyynes are natural compounds characterized by alternating single and triple carbon-carbon bonds, endowing them with distinct physicochemical traits and a range of biological activities. While traditionally sourced mainly from plants, recent investigations have revealed many compounds originating from bacterial strains. This review synthesizes current research on bacterial-derived conjugated polyynes, delving into their biosynthetic routes, underscoring the variety in their molecular structures, and examining their potential applications in biotechnology. Additionally, we outline future directions for metabolic and protein engineering to establish more robust and stable platforms for their production.

ABC-HuMi: the Atlas of Biosynthetic Gene Clusters in the Human Microbiome

The human microbiome has emerged as a rich source of diverse and bioactive natural products, harboring immense potential for therapeutic applications. To facilitate systematic exploration and analysis of its biosynthetic landscape, we present ABC-HuMi: the Atlas of Biosynthetic Gene Clusters (BGCs) in the Human Microbiome. ABC-HuMi integrates data from major human microbiome sequence databases and provides an expansive repository of BGCs compared to the limited coverage offered by existing resources. Employing state-of-the-art BGC prediction and analysis tools, our database ensures accurate annotation and enhanced prediction capabilities. ABC-HuMi empowers researchers with advanced browsing, filtering, and search functionality, enabling efficient exploration of the resource. At present, ABC-HuMi boasts a catalog of 19 218 representative BGCs derived from the human gut, oral, skin, respiratory and urogenital systems. By capturing the intricate biosynthetic potential across diverse human body sites, our database fosters profound insights into the molecular repertoire encoded within the human microbiome and offers a comprehensive resource for the discovery and characterization of novel bioactive compounds. The database is freely accessible at https://www.ccb.uni-saarland.de/abc_humi/.

Bringing up to date the toolkit for the catabolism of aromatic compounds in fungi: The unexpected 1,2,3,5-tetrahydroxybenzene central pathway

Saprophytic fungi are able to catabolize many plant-derived aromatics, including, for example, gallate. The catabolism of gallate in fungi is assumed to depend on the five main central pathways, i.e., of the central intermediates' catechol, protocatechuate, hydroxyquinol, homogentisate and gentisate, but a definitive demonstration is lacking. To shed light on this process, we analysed the transcriptional reprogramming of the growth of Aspergillus terreus on gallate compared with acetate as the control condition. Surprisingly, the results revealed that the five main central pathways did not exhibit significant positive regulation. Instead, an in-depth analysis identified four highly expressed and upregulated genes that are part of a conserved gene cluster found in numerous species of fungi, though not in Aspergilli. The cluster comprises a monooxygenase gene and a fumarylacetoacetate hydrolase-like gene, which are recognized as key components of catabolic pathways responsible for aromatic compound degradation. The other two genes encode proteins with no reported enzymatic activities. Through functional analyses of gene deletion mutants in Aspergillus nidulans, the conserved short protein with no known domains could be linked to the conversion of the novel metabolite 5-hydroxydienelatone, whereas the DUF3500 gene likely encodes a ring-cleavage enzyme for 1,2,3,5-tetrahydroxybenzene. These significant findings establish the existence of a new 1,2,3,5-tetrahydroxybenzene central pathway for the catabolism of gallate and related compounds (e.g. 2,4,6-trihydroxybenzoate) in numerous fungi where this catabolic gene cluster was observed.

Engineered and total biosynthesis of fungal specialized metabolites

Filamentous fungi produce a very wide range of complex and often bioactive metabolites, demonstrating their inherent ability as hosts of complex biosynthetic pathways. Recent advances in molecular sciences related to fungi have afforded the development of new tools that allow the rational total biosynthesis of highly complex specialized metabolites in a single process. Increasingly, these pathways can also be engineered to produce new metabolites. Engineering can be at the level of gene deletion, gene addition, formation of mixed pathways, engineering of scaffold synthases and engineering of tailoring enzymes. Combination of these approaches with hosts that can metabolize low-value waste streams opens the prospect of one-step syntheses from garbage.

Molecular and evolutionary basis of O-antigenic polysaccharide-driven phage sensitivity in environmental pseudomonads

IMPORTANCE

The application of plant-beneficial microorganisms to protect crop plants is a promising alternative to the usage of chemicals. However, biocontrol research often faces difficulties in implementing this approach due to the inconsistency of the bacterial inoculant to establish itself within the root microbiome. Beneficial bacterial inoculants can be decimated by the presence of their natural predators, notably bacteriophages (also called phages). Thus, it is important to gain knowledge regarding the mechanisms behind phage-bacteria interactions to overcome this challenge. Here, we evidence that the major long O-antigenic polysaccharide (O-PS, O-antigen) of the widely used model plant-beneficial bacterium Pseudomonas protegens CHA0 is the receptor of its natural predator, the phage ΦGP100. We examined the distribution of the gene cluster directing the synthesis of this O-PS and identified signatures of horizontal gene acquisitions. Altogether, our study highlights the importance of bacterial cell surface structure variation in the complex interplay between phages and their Pseudomonas hosts.

Extending the Structural Diversity of Labdane Diterpenoids from Marine-Derived Fungus Talaromyces sp. HDN151403 Using Heterologous Expression

Heterologous biosynthesis has become an effective means to activate fungal silent biosynthetic gene clusters (BGCs) and efficiently utilize fungal genetic resources. Herein, thirteen labdane diterpene derivatives, including five undescribed ones named talarobicins A-E (3-7), were discovered via heterologous expression of a silent BGC (labd) in Aspergillus nidulans. Their structures with absolute configurations were elucidated using extensive MS and NMR spectroscopic methods, as well as electronic circular dichroism (ECD) calculations. These labdanes belong to four skeleton types, and talarobicin B (4) is the first 3,18-dinor-2,3:4,18-diseco-labdane diterpene with the cleavage of the C2-C3 bond in ring A and the decarboxylation at C-3 and C-18. Talarobicin B (4) represents the key intermediate in the biosynthesis of penioxalicin and compound 13. The combinatorial heterologous expression and feeding experiments revealed that the cytochrome P450 enzymes LabdC, LabdE, and LabdF were responsible for catalyzing various chemical reactions, such as oxidation, decarboxylation, and methylation. All of the compounds are noncytotoxic, and compounds 2 and 8 displayed inhibitory effects against methicillin-resistant coagulase-negative staphylococci (MRCNS) and Bacillus cereus.

Distribution and diversity of type VI secretion system clusters in Enterobacter bugandensis and Enterobacter cloacae

Gram-negative bacteria use type VI secretion systems (T6SSs) to antagonize neighbouring cells. Although primarily involved in bacterial competition, the T6SS is also implicated in pathogenesis, biofilm formation and ion scavenging. Enterobacter species belong to the ESKAPE pathogens, and while their antibiotic resistance has been well studied, less is known about their pathogenesis. Here, we investigated the distribution and diversity of T6SS components in isolates of two clinically relevant Enterobacter species, E. cloacae and E. bugandensis. T6SS clusters are grouped into four types (T6SSi-T6SSiv), of which type i can be further divided into six subtypes (i1, i2, i3, i4a, i4b, i5). Analysis of a curated dataset of 31 strains demonstrated that most of them encode T6SS clusters belonging to the T6SSi type. All T6SS-positive strains possessed a conserved i3 cluster, and many harboured one or two additional i2 clusters. These clusters were less conserved, and some strains displayed evidence of deletion. We focused on a pathogenic E. bugandensis clinical isolate for comprehensive in silico effector prediction, with comparative analyses across the 31 isolates. Several new effector candidates were identified, including an evolved VgrG with a metallopeptidase domain and a Tse6-like protein. Additional effectors included an anti-eukaryotic catalase (KatN), M23 peptidase, PAAR and VgrG proteins. Our findings highlight the diversity of Enterobacter T6SSs and reveal new putative effectors that may be important for the interaction of these species with neighbouring cells and their environment.

Tartrolon sensing and detoxification by the Listeria monocytogenes timABR resistance operon

Listeria monocytogenes is a foodborne bacterium that naturally occurs in the soil. Originating from there, it contaminates crops and infects farm animals and their consumption by humans may lead to listeriosis, a systemic life-threatening infectious disease. The adaptation of L. monocytogenes to such contrastive habitats is reflected by the presence of virulence genes for host infection and other genes for survival under environmental conditions. Among the latter are ABC transporters for excretion of antibiotics produced by environmental competitors; however, most of these transporters have not been characterized. Here, we generated a collection of promoter-lacZ fusions for genes encoding ABC-type drug transporters of L. monocytogenes and screened this reporter strain collection for induction using a library of natural compounds produced by various environmental microorganisms. We found that the timABR locus (lmo1964-lmo1962) was induced by the macrodiolide antibiotic tartrolon B, which is synthesized by the soil myxobacterium Sorangium cellulosum. Tartrolon B resistance of L. monocytogenes was dependent on timAB, encoding the ATPase and the permease component of a novel ABC transporter. Moreover, transplantation of timAB was sufficient to confer tartrolon B resistance to Bacillus subtilis. Expression of the timABR locus was found to be auto-repressed by the TimR repressor, whose repressing activity was lost in the presence of tartrolon B. We also demonstrate that tartrolon sensitivity was suppressed by high external potassium concentrations, suggesting that tartrolon acts as potassium ionophore. Our results help to map the ecological interactions of an important human pathogen with its co-residing species within their joint natural reservoir.

Functional diversification despite structural congruence in the HipBST toxin-antitoxin system of Legionella pneumophila

IMPORTANCE

Toxin-antitoxin (TA) systems are parasitic genetic elements found in almost all bacterial genomes. They are exchanged horizontally between cells and are typically poorly conserved across closely related strains and species. Here, we report the characterization of a tripartite TA system in the bacterial pathogen Legionella pneumophila that is highly conserved across Legionella species genomes. This system (denoted HipBSTLp) is a distant homolog of the recently discovered split-HipA system in Escherichia coli (HipBSTEc). We present bioinformatic, molecular, and structural analyses of the divergence between these two systems and the functionality of this newly described TA system family. Furthermore, we provide evidence to refute previous claims that the toxin in this system (HipTLp) possesses bifunctionality as an L. pneumophila virulence protein. Overall, this work expands our understanding of the split-HipA system architecture and illustrates the potential for undiscovered biology in these abundant genetic elements.

Fungal Duel between Penicillium brasilianum and Aspergillus nomius Results in Dual Induction of Miktospiromide A and Kitrinomycin A

Cocultivation of the fungi Penicillium brasilianum MST-FP1927 and Aspergillus nomius MST-FP2004 resulted in the reciprocal induction of two new compounds, miktospiromide A (1) from A. nomius and kitrinomycin A (2) from P. brasilianum. A third new compound, kitrinomycin B (3), was also identified from an axenic culture of P. brasilianum, along with the previously reported compounds austalide K (4), 17S-dihydroaustalide K (5), verruculogen (6), and fumitremorgin B (7). The structures of 1-3 were elucidated by detailed spectroscopic analysis and DFT calculations, while 4-7 were identified by comparison to authentic standards. The genome of A. nomius MST-FP2004 was sequenced, and a putative biosynthetic gene cluster for 1 was identified. Compound 2 showed activity against murine melanoma NS-1 cells (LD99 7.8 μM) and the bovine parasite Tritrichomonas foetus (LD99 4.8 μM).

Exploring the biosynthetic gene clusters in Brevibacterium: a comparative genomic analysis of diversity and distribution

Exploring Brevibacterium strains from various ecosystems may lead to the discovery of new antibiotic-producing strains. Brevibacterium sp. H-BE7, a strain isolated from marine sediments from Northern Patagonia, Chile, had its genome sequenced to study the biosynthetic potential to produce novel natural products within the Brevibacterium genus. The genome sequences of 98 Brevibacterium strains, including strain H-BE7, were selected for a genomic analysis. A phylogenomic cladogram was generated, which divided the Brevibacterium strains into four major clades. A total of 25 strains are potentially unique new species according to Average Nucleotide Identity (ANIb) values. These strains were isolated from various environments, emphasizing the importance of exploring diverse ecosystems to discover the full diversity of Brevibacterium. Pangenome analysis of Brevibacterium strains revealed that only 2.5% of gene clusters are included within the core genome, and most gene clusters occur either as singletons or as cloud genes present in less than ten strains. Brevibacterium strains from various phylogenomic clades exhibit diverse BGCs. Specific groups of BGCs show clade-specific distribution patterns, such as siderophore BGCs and carotenoid-related BGCs. A group of clade IV-A Brevibacterium strains possess a clade-specific Polyketide synthase (PKS) BGCs that connects with phenazine-related BGCs. Within the PKS BGC, five genes, including the biosynthetic PKS gene, participate in the mevalonate pathway and exhibit similarities with the phenazine A BGC. However, additional core biosynthetic phenazine genes were exclusively discovered in nine Brevibacterium strains, primarily isolated from cheese. Evaluating the antibacterial activity of strain H-BE7, it exhibited antimicrobial activity against Salmonella enterica and Listeria monocytogenes. Chemical dereplication identified bioactive compounds, such as 1-methoxyphenazine in the crude extracts of strain H-BE7, which could be responsible of the observed antibacterial activity. While strain H-BE7 lacks the core phenazine biosynthetic genes, it produces 1-methoxyphenazine, indicating the presence of an unknown biosynthetic pathway for this compound. This suggests the existence of alternative biosynthetic pathways or promiscuous enzymes within H-BE7's genome.

Design and validation of a PCR screen for γ-butyrolactone-like regulatory systems in Streptomyces

γ-butyrolactone and related signalling systems are found in Streptomyces and other actinobacteria where they control the production of secondary or specialized metabolites such as antibiotics. Genetic manipulation of these regulatory systems therefore leads to changes in the secondary metabolite profile of a strain and has been used to activate previously silent secondary metabolite gene clusters. However, there is no easy way to assess the presence of γ-butyrolactone-like systems in Streptomyces strains without whole-genome sequencing. We have therefore developed and tested a PCR screen that is able to detect homologues of the commonly co-located butenolide synthase and γ-butyrolactone receptor genes. This PCR screen could be employed for the screening of strain libraries to detect signalling systems without the necessity for whole-genome sequencing.

Marine Bacteroidetes enzymatically digest xylans from terrestrial plants

Marine Bacteroidetes that degrade polysaccharides contribute to carbon cycling in the ocean. Organic matter, including glycans from terrestrial plants, might enter the oceans through rivers. Whether marine bacteria degrade structurally related glycans from diverse sources including terrestrial plants and marine algae was previously unknown. We show that the marine bacterium Flavimarina sp. Hel_I_48 encodes two polysaccharide utilization loci (PULs) which degrade xylans from terrestrial plants and marine algae. Biochemical experiments revealed activity and specificity of the encoded xylanases and associated enzymes of these PULs. Proteomics indicated that these genomic regions respond to glucuronoxylans and arabinoxylans. Substrate specificities of key enzymes suggest dedicated metabolic pathways for xylan utilization. Some of the xylanases were active on different xylans with the conserved β-1,4-linked xylose main chain. Enzyme activity was consistent with growth curves showing Flavimarina sp. Hel_I_48 uses structurally different xylans. The observed abundance of related xylan-degrading enzyme repertoires in genomes of other marine Bacteroidetes indicates similar activities are common in the ocean. The here presented data show that certain marine bacteria are genetically and biochemically variable enough to access parts of structurally diverse xylans from terrestrial plants as well as from marine algal sources.

Genomes and transcriptomes help unravel the complex life cycle of the blastoclad fungus, Coelomomyces lativittatus, an obligate parasite of mosquitoes and microcrustaceans

Species of the phylum Blastocladiomycota, early-diverging zoosporic (flagellated) lineages of fungi, are vastly understudied. This phylum includes the genus Coelomomyces, which consists of more than 80 fungal species that are obligate parasites of arthropods. Known Coelomomyces species lack a complete asexual life cycle, instead surviving through an obligate heterecious alternation of generations life cycle. Despite their global distribution and interesting life cycle, little is known about the genomics of any Coelomomyces species. To address this, we generated three draft-level genomes and annotations for C. lativittatus representing its haploid meiospore, orange gamete, and amber gamete life stages. These draft genome assemblies ranged in size from 5002 to 5799 contigs, with a total length of 19.8-22.8 Mb and a mean of 7416 protein-coding genes. We then demonstrated the utility of these genomes by combining the draft annotations as a reference for analysis of C. lativittatus transcriptomes. We analyzed transcriptomes from across host-associated life stages, including infected larvae and excised mature sporangia from the mosquito Anopheles quadrimaculatus. We identified differentially expressed genes and enriched GO terms both across and within life stages and used these to make hypotheses about C. lativittatus biology. Generally, we found the C. lativittatus transcriptome to be a complex and dynamic expression landscape; GO terms related to metabolism and transport processes were enriched during infection and terms related to dispersal were enriched during sporulation. We further identified five high mobility group (HMG)-box genes in C. lativittatus, three belonging to clades with mating type (MAT) loci from other fungi, as well as four ortholog expansions in C. lativittatus compared with other fungi. The C. lativittatus genomes and transcriptomes reported here are a valuable resource and may be leveraged toward furthering understanding of the biology of these and other early-diverging fungal lineages.

Relation of pest insect-killing and soilborne pathogen-inhibition abilities to species diversification in environmental Pseudomonas protegens

Strains belonging to the Pseudomonas protegens phylogenomic subgroup have long been known for their beneficial association with plant roots, notably antagonising soilborne phytopathogens. Interestingly, they can also infect and kill pest insects, emphasising their interest as biocontrol agents. In the present study, we used all available Pseudomonas genomes to reassess the phylogeny of this subgroup. Clustering analysis revealed the presence of 12 distinct species, many of which were previously unknown. The differences between these species also extend to the phenotypic level. Most of the species were able to antagonise two soilborne phytopathogens, Fusarium graminearum and Pythium ultimum, and to kill the plant pest insect Pieris brassicae in feeding and systemic infection assays. However, four strains failed to do so, likely as a consequence of adaptation to particular niches. The absence of the insecticidal Fit toxin explained the non-pathogenic behaviour of the four strains towards Pieris brassicae. Further analyses of the Fit toxin genomic island evidence that the loss of this toxin is related to non-insecticidal niche specialisation. This work expands the knowledge on the growing Pseudomonas protegens subgroup and suggests that loss of phytopathogen inhibition and pest insect killing abilities in some of these bacteria may be linked to species diversification processes involving adaptation to particular niches. Our work sheds light on the important ecological consequences of gain and loss dynamics for functions involved in pathogenic host interactions of environmental bacteria.

Genomic Based Analysis of the Biocontrol Species Trichoderma harzianum: A Model Resource of Structurally Diverse Pharmaceuticals and Biopesticides

Fungi represents a rich repository of taxonomically restricted, yet chemically diverse, secondary metabolites that are synthesised via specific metabolic pathways. An enzyme's specificity and biosynthetic gene clustering are the bottleneck of secondary metabolite evolution. Trichoderma harzianum M10 v1.0 produces many pharmaceutically important molecules; however, their specific biosynthetic pathways remain uncharacterised. Our genomic-based analysis of this species reveals the biosynthetic diversity of its specialised secondary metabolites, where over 50 BGCs were predicted, most of which were listed as polyketide-like compounds associated clusters. Gene annotation of the biosynthetic candidate genes predicted the production of many medically/industrially important compounds including enterobactin, gramicidin, lovastatin, HC-toxin, tyrocidine, equisetin, erythronolide, strobilurin, asperfuranone, cirtinine, protoilludene, germacrene, and epi-isozizaene. Revealing the biogenetic background of these natural molecules is a step forward towards the expansion of their chemical diversification via engineering their biosynthetic genes heterologously, and the identification of their role in the interaction between this fungus and its biotic/abiotic conditions as well as its role as bio-fungicide.

Anaerobic phloroglucinol degradation by Clostridium scatologenes

Polyphenols are abundant in nature, and their anaerobic biodegradation by gut and soil bacteria is a topic of great interest. The O2 requirement of phenol oxidases is thought to explain the microbial inertness of phenolic compounds in anoxic environments, such as peatlands, termed the enzyme latch hypothesis. A caveat of this model is that certain phenols are known to be degraded by strict anaerobic bacteria, although the biochemical basis for this process is incompletely understood. Here, we report the discovery and characterization of a gene cluster in the environmental bacterium Clostridium scatologenes for the degradation phloroglucinol (1,3,5-trihydroxybenzene), a key intermediate in the anaerobic degradation of flavonoids and tannins, which constitute the most abundant polyphenols in nature. The gene cluster encodes the key C-C cleavage enzyme dihydrophloroglucinol cyclohydrolase, as well as (S)-3-hydroxy-5-oxo-hexanoate dehydrogenase and triacetate acetoacetate-lyase, which enable phloroglucinol to be utilized as a carbon and energy source. Bioinformatics studies revealed the presence of this gene cluster in phylogenetically and metabolically diverse gut and environmental bacteria, with potential impacts on human health and carbon preservation in peat soils and other anaerobic environmental niches. IMPORTANCE This study provides novel insights into the microbiota's anaerobic metabolism of phloroglucinol, a critical intermediate in the degradation of polyphenols in plants. Elucidation of this anaerobic pathway reveals enzymatic mechanisms for the degradation of phloroglucinol into short-chain fatty acids and acetyl-CoA, which are used as a carbon and energy source for bacterium growth. Bioinformatics studies suggested the prevalence of this pathway in phylogenetically and metabolically diverse gut and environmental bacteria, with potential impacts on carbon preservation in peat soils and human gut health.

The Achromobacter type 3 secretion system drives pyroptosis and immunopathology via independent activation of NLRC4 and NLRP3 inflammasomes

How the opportunistic Gram-negative pathogens of the genus Achromobacter interact with the innate immune system is poorly understood. Using three Achromobacter clinical isolates from two species, we show that the type 3 secretion system (T3SS) is required to induce cell death in human macrophages by inflammasome-dependent pyroptosis. Macrophages deficient in the inflammasome sensors NLRC4 or NLRP3 undergo pyroptosis upon bacterial internalization, but those deficient in both NLRC4 and NLRP3 do not, suggesting either sensor mediates pyroptosis in a T3SS-dependent manner. Detailed analysis of the intracellular trafficking of one isolate indicates that the intracellular bacteria reside in a late phagolysosome. Using an intranasal mouse infection model, we observe that Achromobacter damages lung structure and causes severe illness, contingent on a functional T3SS. Together, we demonstrate that Achromobacter species can survive phagocytosis by promoting macrophage cell death and inflammation by redundant mechanisms of pyroptosis induction in a T3SS-dependent manner.

Genomic Analysis of Bacteriocin-Producing Staphylococci: High Prevalence of Lanthipeptides and the Micrococcin P1 Biosynthetic Gene Clusters

Bacteriocins are antimicrobial peptides produced by bacteria. This study aimed to in silico analyze the presence of bacteriocin gene clusters (BGCs) among the genomes of 22 commensal Staphylococcus isolates from different origins (environment/human/food/pet/wild animals) previously identified as bacteriocin producers. The resistome and plasmidome were studied in all isolates. Five types of BGC were detected in 18 genomes of the 22 bacteriocin-producing staphylococci included in this study: class I (Lanthipeptides), class II, circular bacteriocins, the non-ribosomal-peptide lugdunin and the thiopeptide micrococcin P1 (MP1). A high frequency of lanthipeptides was detected in this collection: BGC variants of BSA, bacCH91, and epilancin15X were identified in two Staphylococcus aureus and one Staphylococcus warneri isolates from food and wild animals. Moreover, two potentially new lanthipeptide-like BGCs with no identity to database entries were found in Staphylococcus epidermidis and Staphylococcus simulans from food and wild animal, respectively. Interestingly, four isolates (one S. aureus and one Staphylococcus hominis, environmental origin; two Staphylococcus sciuri, food) carried the MP1 BGC with differences to those previously described. On the other hand, seven of the 22 genomes (~32%) lacked known genes related with antibiotic or disinfectant-acquired resistance mechanisms. Moreover, the potential carriage of plasmids was evaluated, and several Rep-proteins were identified (~73% of strains). In conclusion, a wide variety of BGCs has been observed among the 22 genomes, and an interesting relationship between related Staphylococcus species and the type of bacteriocin has been revealed. Therefore, bacteriocin-producing Staphylococcus and especially coagulase-negative staphylococci (CoNS) can be considered good candidates as a source of novel bacteriocins.

Tolyporphins-Exotic Tetrapyrrole Pigments in a Cyanobacterium-A Review

Tolyporphins were discovered some 30 years ago as part of a global search for antineoplastic compounds from cyanobacteria. To date, the culture HT-58-2, comprised of a cyanobacterium-microbial consortium, is the sole known producer of tolyporphins. Eighteen tolyporphins are now known-each is a free base tetrapyrrole macrocycle with a dioxobacteriochlorin (14), oxochlorin (3), or porphyrin (1) chromophore. Each compound displays two, three, or four open β-pyrrole positions and two, one, or zero appended C-glycoside (or -OH or -OAc) groups, respectively; the appended groups form part of a geminal disubstitution motif flanking the oxo moiety in the pyrroline ring. The distinct structures and repertoire of tolyporphins stand alone in the large pigments-of-life family. Efforts to understand the cyanobacterial origin, biosynthetic pathways, structural diversity, physiological roles, and potential pharmacological properties of tolyporphins have attracted a broad spectrum of researchers from diverse scientific areas. The identification of putative biosynthetic gene clusters in the HT-58-2 cyanobacterial genome and accompanying studies suggest a new biosynthetic paradigm in the tetrapyrrole arena. The present review provides a comprehensive treatment of the rich science concerning tolyporphins.

Paraburkholderia sabiae Uses One Type VI Secretion System (T6SS-1) as a Powerful Weapon against Notorious Plant Pathogens

Paraburkholderia sabiae LMG24235 is a nitrogen-fixing betaproteobacterium originally isolated from a root nodule of Mimosa caesalpiniifolia in Brazil. We show here that this strain effectively kills strains from several bacterial families (Burkholderiaceae, Pseudomonadaceae, Enterobacteriaceae) which include important plant pathogens in a contact-dependent manner. De novo assembly of the first complete genome of P. sabiae using long sequencing reads and subsequent annotation revealed two gene clusters predicted to encode type VI secretion systems (T6SS), which we named T6SS-1 and T6SS-3 according to previous classification methods (G. Shalom, J. G. Shaw, and M. S. Thomas, Microbiology, 153:2689-2699, 2007, https://doi.org/10.1099/mic.0.2007/006585-0). We created P. sabiae with mutations in each of the two T6SS gene clusters that abrogated their function, and the T6SS-1 mutant was no longer able to outcompete other strains in a contact-dependent manner. Notably, our analysis revealed that T6SS-1 is essential for competition against several important plant pathogens in vitro, including Burkholderia plantarii, Ralstonia solanacearum, Pseudomonas syringae, and Pectobacterium carotovorum. The 9-log reduction in P. syringae cells in the presence of P. sabiae was particularly remarkable. Importantly, in an in vivo assay, P. sabiae was able to protect potato tubers from bacterial soft rot disease caused by P. carotovorum, and this protection was partly dependent on T6SS-1. IMPORTANCE Rhizobia often display additional beneficial traits such as the production of plant hormones and the acquisition of limited essential nutrients that improve plant growth and enhance plant yields. Here, we show that the rhizobial strain P. sabiae antagonizes important phytopathogens such as P. carotovorum, P. syringae, and R. solanacearum and that this effect is due to contact-dependent killing mediated by one of two T6SS systems identified in the complete, de novo assembled genome sequence of P. sabiae. Importantly, co-inoculation of Solanum tuberosum tubers with P. sabiae also resulted in a drastic reduction of soft rot caused by P. carotovorum in an in vivo model system. This result highlights the protective potential of P. sabiae against important bacterial plant diseases, which makes it a valuable candidate for application as a biocontrol agent. It also emphasizes the particular potential of rhizobial inoculants that combine several beneficial effects such as plant growth promotion and biocontrol for sustainable agriculture.

Uncovering Genomic Features and Biosynthetic Gene Clusters in Endophytic Bacteria from Roots of the Medicinal Plant Alkanna tinctoria Tausch as a Strategy To Identify Novel Biocontrol Bacteria

The world's population is increasing at a rate not seen in the past. Agriculture, providing food for this increasing population, is reaching its boundaries of space and natural resources. In addition, changing legislation and increased ecological awareness are forcing agriculture to reduce its environmental impact. This entails the replacement of agrochemicals with nature-based solutions. In this regard, the search for effective biocontrol agents that protect crops from pathogens is in the spotlight. In this study, we have investigated the biocontrol activity of endophytic bacteria isolated from the medicinal plant Alkanna tinctoria Tausch. To do so, an extensive collection of bacterial strains was initially genome sequenced and in silico screened for features related to plant stimulation and biocontrol. Based on this information, a selection of bacteria was tested in vitro for antifungal activity using direct antagonism in a plate assay and in planta with a detached-leaf assay. Bacterial strains were tested individually and in combinations to assess the best-performing treatments. The results revealed that many bacteria could produce metabolites that efficiently inhibit the proliferation of several fungi, especially Fusarium graminearum. Among these, Pseudomonas sp. strain R-71838 showed a strong antifungal effect, in both dual-culture and in planta assays, making it the most promising candidate for biocontrol application. Using microbes from medicinal plants, this study highlights the opportunities of using genomic information to speed up the screening of a taxonomically diverse set of bacteria with biocontrol properties. IMPORTANCE Phytopathogenic fungi are a major threat to global food production. The most common management practice to prevent plant infections involves the intensive use of fungicides. However, with the growing awareness of the ecological and human impacts of chemicals, there is a need for alternative strategies, such as the use of bacterial biocontrol agents. Limitations in the design of bacterial biocontrol included the need for labor-intensive and time-consuming experiments to test a wide diversity of strains and the lack of reproducibility of their activity against pathogens. Here, we show that genomic information is an effective tool to select bacteria of interest quickly. Also, we highlight that the strain Pseudomonas sp. R-71838 produced a reproducible antifungal effect both in vitro and in planta. These findings build a foundation for designing a biocontrol strategy based on Pseudomonas sp. R-71838.

Mining for a new class of fungal natural products: the evolution, diversity, and distribution of isocyanide synthase biosynthetic gene clusters

The products of non-canonical isocyanide synthase (ICS) biosynthetic gene clusters (BGCs) mediate pathogenesis, microbial competition, and metal-homeostasis through metal-associated chemistry. We sought to enable research into this class of compounds by characterizing the biosynthetic potential and evolutionary history of these BGCs across the Fungal Kingdom. We amalgamated a pipeline of tools to predict BGCs based on shared promoter motifs and located 3800 ICS BGCs in 3300 genomes, making ICS BGCs the fifth largest class of specialized metabolites compared to canonical classes found by antiSMASH. ICS BGCs are not evenly distributed across fungi, with evidence of gene-family expansions in several Ascomycete families. We show that the ICS dit1/2 gene cluster family (GCF), which was prior only studied in yeast, is present in ∼30% of all Ascomycetes. The dit variety ICS exhibits greater similarity to bacterial ICS than other fungal ICS, suggesting a potential convergence of the ICS backbone domain. The evolutionary origins of the dit GCF in Ascomycota are ancient and these genes are diversifying in some lineages. Our results create a roadmap for future research into ICS BGCs. We developed a website (https://isocyanides.fungi.wisc.edu/) that facilitates the exploration and downloading of all identified fungal ICS BGCs and GCFs.

In host evolution of Exophiala dermatitidis in cystic fibrosis lung micro-environment

Individuals with cystic fibrosis (CF) are susceptible to chronic lung infections that lead to inflammation and irreversible lung damage. While most respiratory infections that occur in CF are caused by bacteria, some are dominated by fungi such as the slow-growing black yeast Exophiala dermatitidis. Here, we analyze isolates of E. dermatitidis cultured from two samples, collected from a single subject 2 years apart. One isolate genome was sequenced using long-read Nanopore technology as an in-population reference to use in comparative single nucleotide polymorphism and insertion-deletion variant analyses of 23 isolates. We then used population genomics and phylo-genomics to compare the isolates to each other as well as the reference genome strain E. dermatitidis NIH/UT8656. Within the CF lung population, three E. dermatitidis clades were detected, each with varying mutation rates. Overall, the isolates were highly similar suggesting that they were recently diverged. All isolates were MAT 1-1, which was consistent with their high relatedness and the absence of evidence for mating or recombination between isolates. Phylogenetic analysis grouped sets of isolates into clades that contained isolates from both early and late time points indicating there are multiple persistent lineages. Functional assessment of variants unique to each clade identified alleles in genes that encode transporters, cytochrome P450 oxidoreductases, iron acquisition, and DNA repair processes. Consistent with the genomic heterogeneity, isolates showed some stable phenotype heterogeneity in melanin production, subtle differences in antifungal minimum inhibitory concentrations, and growth on different substrates. The persistent population heterogeneity identified in lung-derived isolates is an important factor to consider in the study of chronic fungal infections, and the analysis of changes in fungal pathogens over time may provide important insights into the physiology of black yeasts and other slow-growing fungi in vivo.

The conserved iol gene cluster in Pseudomonas is involved in rhizosphere competence

The Pseudomonas genus has shown great potential as a sustainable solution to support agriculture through its plant-growth-promoting and biocontrol activities. However, their efficacy as bioinoculants is limited by unpredictable colonization in natural conditions. Our study identifies the iol locus, a gene cluster in Pseudomonas involved in inositol catabolism, as a feature enriched among superior root colonizers in natural soil. Further characterization revealed that the iol locus increases competitiveness, potentially caused by an observed induction of swimming motility and the production of fluorescent siderophore in response to inositol, a plant-derived compound. Public data analyses indicate that the iol locus is broadly conserved in the Pseudomonas genus and linked to diverse host-microbe interactions. Together, our findings suggest the iol locus as a potential target for developing more effective bioinoculants for sustainable agriculture.

Recently evolved combination of unique sulfatase and amidase genes enables bacterial degradation of the wastewater micropollutant acesulfame worldwide

Xenobiotics often challenge the principle of microbial infallibility. One example is acesulfame introduced in the 1980s as zero-calorie sweetener, which was recalcitrant in wastewater treatment plants until the early 2010s. Then, efficient removal has been reported with increasing frequency. By studying acesulfame metabolism in alphaproteobacterial degraders of the genera Bosea and Chelatococcus, we experimentally confirmed the previously postulated route of two subsequent hydrolysis steps via acetoacetamide-N-sulfonate (ANSA) to acetoacetate and sulfamate. Genome comparison of wildtype Bosea sp. 100-5 and an acesulfame degradation-defective mutant revealed the involvement of two plasmid-borne gene clusters. The acesulfame-hydrolyzing sulfatase is strictly manganese-dependent and belongs to the metallo beta-lactamase family. In all degraders analyzed, it is encoded on a highly conserved gene cluster embedded in a composite transposon. The ANSA amidase, on the other hand, is an amidase signature domain enzyme encoded in another gene cluster showing variable length among degrading strains. Transposition of the sulfatase gene cluster between chromosome and plasmid explains how the two catabolic gene clusters recently combined for the degradation of acesulfame. Searching available genomes and metagenomes for the two hydrolases and associated genes indicates that the acesulfame plasmid evolved and spread worldwide in short time. While the sulfatase is unprecedented and unique for acesulfame degraders, the amidase occurs in different genetic environments and likely evolved for the degradation of other substrates. Evolution of the acesulfame degradation pathway might have been supported by the presence of structurally related natural and anthropogenic compounds, such as aminoacyl sulfamate ribonucleotide or sulfonamide antibiotics.

zol & fai: large-scale targeted detection and evolutionary investigation of gene clusters

Many universally and conditionally important genes are genomically aggregated within clusters. Here, we introduce fai and zol, which together enable large-scale comparative analysis of different types of gene clusters and mobile-genetic elements (MGEs), such as biosynthetic gene clusters (BGCs) or viruses. Fundamentally, they overcome a current bottleneck to reliably perform comprehensive orthology inference at large scale across broad taxonomic contexts and thousands of genomes. First, fai allows the identification of orthologous or homologous instances of a query gene cluster of interest amongst a database of target genomes. Subsequently, zol enables reliable, context-specific inference of protein-encoding ortholog groups for individual genes across gene cluster instances. In addition, zol performs functional annotation and computes a variety of statistics for each inferred ortholog group. These programs are showcased through application to: (i) longitudinal tracking of a virus in metagenomes, (ii) discovering novel population-genetic insights of two common BGCs in a fungal species, and (iii) uncovering large-scale evolutionary trends of a virulence-associated gene cluster across thousands of genomes from a diverse bacterial genus.

Discovery of Antimicrobial Phosphonopeptide Natural Products from Bacillus velezensis by Genome Mining

Phosphonate natural products are renowned for inhibitory activities which underly their development as antibiotics and pesticides. Although most phosphonate natural products have been isolated from Streptomyces, bioinformatic surveys suggest that many other bacterial genera are replete with similar biosynthetic potential. While mining actinobacterial genomes, we encountered a contaminated Mycobacteroides data set which included a biosynthetic gene cluster predicted to produce novel phosphonate compounds. Sequence deconvolution revealed that the contig containing this cluster, as well as many others, belonged to a contaminating Bacillus and is broadly conserved among multiple species, including the epiphyte Bacillus velezensis. Isolation and structure elucidation revealed a new di- and tripeptide composed of l-alanine and a C-terminal l-phosphonoalanine which we name phosphonoalamides E and F. These compounds exhibit broad-spectrum antibacterial activity, including strong inhibition against the agricultural pests responsible for vegetable soft rot (Erwinia rhapontici), onion rot (Pantoea ananatis), and American foulbrood (Paenibacillus larvae). This work expands our knowledge of phosphonate metabolism and underscores the importance of including underexplored microbial taxa in natural product discovery. IMPORTANCE Phosphonate natural products produced by bacteria have been a rich source of clinical antibiotics and commercial pesticides. Here, we describe the discovery of two new phosphonopeptides produced by B. velezensis with antibacterial activity against human and plant pathogens, including those responsible for widespread soft rot in crops and American foulbrood. Our results provide new insight on the natural chemical diversity of phosphonates and suggest that these compounds could be developed as effective antibiotics for use in medicine or agriculture.

Bioinformatics of cyanophycin metabolism genes and characterization of promiscuous isoaspartyl dipeptidases that catalyze the final step of cyanophycin degradation

Cyanophycin is a bacterial biopolymer used for storage of fixed nitrogen. It is composed of a backbone of L-aspartate residues with L-arginines attached to each of their side chains. Cyanophycin is produced by cyanophycin synthetase 1 (CphA1) using Arg, Asp and ATP, and is degraded in two steps. First, cyanophycinase breaks down the backbone peptide bonds, releasing β-Asp-Arg dipeptides. Then, these dipeptides are broken down into free Asp and Arg by enzymes with isoaspartyl dipeptidase activity. Two bacterial enzymes are known to possess promiscuous isoaspartyl dipeptidase activity: isoaspartyl dipeptidase (IadA) and isoaspartyl aminopeptidase (IaaA). We performed a bioinformatic analysis to investigate whether genes for cyanophycin metabolism enzymes cluster together or are spread around the microbial genomes. Many genomes showed incomplete contingents of known cyanophycin metabolizing genes, with different patterns in various bacterial clades. Cyanophycin synthetase and cyanophycinase are usually clustered together when recognizable genes for each are found within a genome. Cyanophycinase and isoaspartyl dipeptidase genes typically cluster within genomes lacking cphA1. About one-third of genomes with genes for CphA1, cyanophycinase and IaaA show these genes clustered together, while the proportion is around one-sixth for CphA1, cyanophycinase and IadA. We used X-ray crystallography and biochemical studies to characterize an IadA and an IaaA from two such clusters, in Leucothrix mucor and Roseivivax halodurans, respectively. The enzymes retained their promiscuous nature, showing that being associated with cyanophycin-related genes did not make them specific for β-Asp-Arg dipeptides derived from cyanophycin degradation.

Ketoreductase Domain-Catalyzed Polyketide Chain Release in Fungal Alkyl Salicylaldehyde Biosynthesis

Alkyl salicylaldehyde derivatives are polyketide natural products, which are widely distributed in fungi and exhibit great structural diversity. Their biosynthetic mechanisms have recently been intensively studied; however, how the polyketide synthases (PKSs) involved in the fungal alkyl salicylaldehyde biosyntheses release their products remained elusive. In this study, we discovered an orphan biosynthetic gene cluster of salicylaldehyde derivatives in the fungus Stachybotrys sp. g12. Intriguingly, the highly reducing PKS StrA, encoded by the gene cluster, performs a reductive polyketide chain release, although it lacks a C-terminal reductase domain, which is typically required for such a reductive release. Our study revealed that the chain release is achieved by the ketoreductase (KR) domain of StrA, which also conducts cannonical β-keto reductions during polyketide chain elongation. Furthermore, we found that the cupin domain-containing protein StrC plays a critical role in the aromatization reaction. Collectively, we have provided an unprecedented example of a KR domain-catalyzed polyketide chain release and a clearer image of how the salicylaldehyde scaffold is generated in fungi.

Human-associated bacteria adopt an unusual route for synthesizing 3-acetylated tetramates for environmental adaptation

BACKGROUND

Tetramates or tetramic acid-containing compounds (TACs) are a group of bioactive natural products featuring a pyrrolidine-2,4-dione ring acknowledged being closed via Dieckmann cyclization. The cariogenic Streptococcus mutans strains bearing a muc biosynthetic gene cluster (BGC) can synthesize mutanocyclin (MUC), a 3-acetylated TAC that can inhibit both leukocyte chemotaxis and filamentous development in Candida albicans. Some strains can also accumulate reutericyclins (RTCs), the intermediates of MUC biosynthesis with antibacterial activities. However, the formation mechanism of the pyrrolidine-2,4-dione ring of MUC and the distribution of muc-like BGCs along with their ecological functions has not been explored extensively.

RESULTS

We demonstrated that a key intermediate of MUC biosynthesis, M-307, is installed by a hybrid nonribosomal peptide synthetase-polyketide synthase assembly line and its pyrrolidine-2,4-dione ring is closed via an unprecedented lactam bond formation style. Subsequent C-3 acetylation will convert M-307 to RTCs, which is then hydrolyzed by a deacylase, MucF, to remove the N-1 fatty acyl appendage to generate MUC. Distribution analysis showed that the muc-like BGCs distribute predominantly in human-associated bacteria. Interestingly, most of the muc-like BGCs possessing a mucF gene were isolated from human or livestock directly, indicating their involvement in alleviating the host's immune attacks by synthesizing MUC; while those BGCs lacking mucF gene distribute mainly in bacteria from fermented products, suggesting that they tend to synthesize RTCs to compete with neighboring bacteria. It is noteworthy that many bacteria in the same habitats (e.g., the oral cavity) lack the muc-like BGC, but possess functional MucF homologues to "detoxify" RTCs to MUC, including several competitive bacteria of S. mutans. We also comparably studied the distribution of TAS1, a fungal enzyme responsible for the production of phytotoxic tenuazonic acids (TeAs), a class of 3-acetylated TACs with similar structure but distinct biosynthetic mechanism to MUC, and found that it mainly exists in plants or crops.

CONCLUSIONS

The in vivo and in vitro experiments revealed that the pyrrolidine-2,4-dione ring of MUC is closed via lactam bond formation, which may be adopted by many TACs without 3-acyl decorations. Besides, we found that muc-like BGCs are widespread in human-associated bacteria and their shapes and main products can be influenced by the habitat environment and vice versa. By comparing with TeAs, we provided thought-provoking insights into how ecological and evolutionary forces drive bacteria and fungi to construct a common 3-acetylated pyrrolidine-2,4-dione core through different routes, and how the biosynthetic processes are delicately controlled to generate diverse 3-acetylated TACs for environmental adaptation. Video Abstract.

Characterization of three succinyl-CoA acyltransferases involved in polyketide chain assembly

Polyketides are a class of natural products with astonishing structural diversities, fascinating biological activities, and a versatile of applications. In polyketides biosynthesis, acyltransferases (ATs) are the 'gatekeeping' enzymes selecting the specific CoA-activated acyl groups as building blocks and transferring them onto the phosphopantetheine arm of acyl carrier proteins (ACPs) to enable the following condensation reactions to assemble the polyketide chain. Herein, the Art2 protein from aurantinins, a group of the antibacterial polyketides, is characterized in vitro as an AT that can load a CoA-activated succinyl unit onto the first ACP domain of Art17 (ACP_Art17-1). In addition, another two proteins, GbnB and EtnB, involved in the biosynthesis of gladiolin and etnangien respectively, were traced by literature mining, homologous searching, and product structure analysis and then identified as functional succinyl-CoA ATs by the ACP_Art17-1 assays. Taken together, by the assay method employing ACP_Art17-1 as an acyl acceptor, we identified three ATs that can introduce a succinyl unit into the polyketide assembly line, which enriches the toolbox of polyketide biosynthetic enzymes and sets a stage for incorporating a succinyl unit into polyketide backbones in synthetic biological manners. KEY POINTS: • Three acyltransferases that are able to load ACP with a succinyl unit were characterized in vitro. • ACP_Art17-1 can be used as an acceptor to assay succinyl-CoA AT from different polyketides. • The succinyl unit can be incorporated into polyketides assembly process.

Evolutionary investigations of the biosynthetic diversity in the skin microbiome using lsaBGC

Bacterial secondary metabolites, synthesized by enzymes encoded in biosynthetic gene clusters (BGCs), can underlie microbiome homeostasis and serve as commercialized products, which have historically been mined from a select group of taxa. While evolutionary approaches have proven beneficial for prioritizing BGCs for experimental characterization efforts to uncover new natural products, dedicated bioinformatics tools designed for comparative and evolutionary analysis of BGCs within focal taxa are limited. We thus developed lineage specific analysis of BGCs (lsaBGC; https://github.com/Kalan-Lab/lsaBGC) to aid exploration of microdiversity and evolutionary trends across homologous groupings of BGCs, gene cluster families (GCFs), in any bacterial taxa of interest. lsaBGC enables rapid and direct identification of GCFs in genomes, calculates evolutionary statistics and conservation for BGC genes, and builds a framework to allow for base resolution mining of novel variants through metagenomic exploration. Through application of the suite to four genera commonly found in skin microbiomes, we uncover new insights into the evolution and diversity of their BGCs. We show that the BGC of the virulence-associated carotenoid staphyloxanthin in Staphylococcus aureus is ubiquitous across the genus Staphylococcus . While one GCF encoding the biosynthesis of staphyloxanthin showcases evidence for plasmid-mediated horizontal gene transfer (HGT) between species, another GCF appears to be transmitted vertically amongst a sub-clade of skin-associated Staphylococcus . Further, the latter GCF, which is well conserved in S. aureus , has been lost in most Staphylococcus epidermidis , which is the most common Staphylococcus species on human skin and is also regarded as a commensal. We also identify thousands of novel single-nucleotide variants (SNVs) within BGCs from the Corynebacterium tuberculostearicum sp. complex, a narrow, multi-species clade that features the most prevalent Corynebacterium in healthy skin microbiomes. Although novel SNVs were approximately 10 times as likely to correspond to synonymous changes when located in the top five percentile of conserved sites, lsaBGC identified SNVs that defied this trend and are predicted to underlie amino acid changes within functionally key enzymatic domains. Ultimately, beyond supporting evolutionary investigations of BGCs, lsaBGC also provides important functionalities to aid efforts for the discovery or directed modification of natural products.

The genome of Bacillus tequilensis EA-CB0015 sheds light into its epiphytic lifestyle and potential as a biocontrol agent

Different Bacillus species have successfully been used as biopesticides against a broad range of plant pathogens. Among these, Bacillus tequilensis EA-CB0015 has shown to efficiently control Black sigatoka disease in banana plants, presumably by mechanisms of adaptation that involve modifying the phyllosphere environment. Here, we report the complete genome of strain EA-CB0015, its precise taxonomic identity, and determined key genetic features that may contribute to its effective biocontrol of plant pathogens. We found that B. tequilensis EA-CB0015 harbors a singular 4 Mb circular chromosome, with 3,951 protein-coding sequences. Multi-locus sequence analysis (MLSA) and average nucleotide identity (ANI) analysis classified strain EA-CB0015 as B. tequilensis. Encoded within its genome are biosynthetic gene clusters (BGCs) for surfactin, iturin, plipastatin, bacillibactin, bacilysin, subtilosin A, sporulation killing factor, and other natural products that may facilitate inter-microbial warfare. Genes for indole-acetic acid (IAA) synthesis, the use of diverse carbon sources, and a multicellular lifestyle involving motility, biofilm formation, quorum sensing, competence, and sporulation suggest EA-CB0015 is adept at colonizing plant surfaces. Defensive mechanisms to survive invading viral infections and preserve genome integrity include putative type I and type II restriction modification (RM) and toxin/antitoxin (TA) systems. The presence of bacteriophage sequences, genomic islands, transposable elements, virulence factors, and antibiotic resistance genes indicate prior occurrences of genetic exchange. Altogether, the genome of EA-CB0015 supports its function as a biocontrol agent against phytopathogens and suggest it has adapted to thrive within phyllosphere environments.

Discovery of diphenyl ether-degrading Streptomyces strains by direct screening based on ether bond-cleaving activity

Diphenyl ethers (DEs), which are widely used in the agricultural and chemical industries, have become hazardous contaminants in the environment. Although several DE-degrading bacteria have been reported, discovering new types of such microorganisms could enhance understanding of the degradation mechanism in the environment. In this study, we used a direct screening method based on detection of ether bond-cleaving activity to screen for microorganisms that degrade 4,4'-dihydroxydiphenyl ether (DHDE) as a model DE. Microorganisms isolated from soil samples were incubated with DHDE, and strains producing hydroquinone via ether bond cleavage were selected using hydroquinone-sensitive Rhodanine reagent. This screening procedure resulted in the isolation of 3 bacteria and 2 fungi that transform DHDE. Interestingly, all of the isolated bacteria belonged to one genus, Streptomyces. To our knowledge, these are the first microorganisms of the genus Streptomyces shown to degrade a DE. Streptomyces sp. TUS-ST3 exhibited high and stable DHDE-degrading activity. HPLC, LC-MS, and GC-MS analyses revealed that strain TUS-ST3 converts DHDE to its hydroxylated analogue and generates hydroquinone as an ether bond-cleavage product. Strain TUS-ST3 also transformed DEs other than DHDE. In addition, glucose-grown TUS-ST3 cells began to transform DHDE after incubation with this compound for 12 h, and produced 75 μM hydroquinone in 72 h. These activities of streptomycetes may play an important role in DE degradation in the environment. We also report the whole genome sequence of strain TUS-ST3.

Stereodivergent Hydroxylation of Berkeleylactones by Penicillium turbatum

Penicillium turbatum has previously been reported to produce A26771B, a 16-membered macrocyclic polyketide with activity against Gram-positive bacteria, mycoplasma, and fungi, as well as the structurally related compounds berkeleylactone E and berkeleylactones I-O. In this work, large-scale cultivation of P. turbatum NRRL 5630 on rice yielded seven new berkeleylactone analogues, berkeleylactone E methyl ester, 14-epi-berkeleylactone F, berkeleylactones P-R, 12-epi-berkeleylactone Q, and 13-epi-berkeleylactone R, and six previously reported analogues, A26771B and berkeleylactones E-G and J-K. The structures of the berkeleylactones were elucidated by detailed analysis of spectroscopic data, molecular modeling, and comparison with literature values. Interestingly, six of the berkeleylactone analogues were isolated as pairs of hydroxy epimers, highlighting how Nature can exploit stereodivergence in biosynthetic pathways to increase chemical diversity. The genome of P. turbatum was sequenced, and a putative gene cluster (bekl) responsible for the biosynthesis of the berkeleylactones was identified. The new berkeleylactone analogues exhibited no significant biological activity against a panel of bacteria, fungi, the parasite Giardia duodenalis, or NS-1 murine myeloma cells, suggesting a hitherto undiscovered biological role.

Bacteroides thetaiotaomicron, a Model Gastrointestinal Tract Species, Prefers Heme as an Iron Source, Yields Protoporphyrin IX as a Product, and Acts as a Heme Reservoir

Members of the phylum Bacteroidetes are abundant in healthy gastrointestinal (GI) tract flora. Bacteroides thetaiotaomicron is a commensal heme auxotroph and representative of this group. Bacteroidetes are sensitive to host dietary iron restriction but proliferate in heme-rich environments that are also associated with colon cancer. We hypothesized that B. thetaiotaomicron may act as a host reservoir for iron and/or heme. In this study, we defined growth-promoting quantities of iron for B. thetaiotaomicron. B. thetaiotaomicron preferentially consumed and hyperaccumulated iron in the form of heme when presented both heme and nonheme iron sources in excess of its growth needs, leading to an estimated 3.6 to 8.4 mg iron in a model GI tract microbiome consisting solely of B. thetaiotaomicron. Protoporphyrin IX was identified as an organic coproduct of heme metabolism, consistent with anaerobic removal of iron from the heme leaving the intact tetrapyrrole as the observed product. Notably, no predicted or discernible pathway for protoporphyrin IX generation exists in B. thetaiotaomicron. Heme metabolism in congeners of B. thetaiotaomicron has previously been associated with the 6-gene hmu operon, based on genetic studies. A bioinformatics survey demonstrated that the intact operon is widespread in but confined to members of the Bacteroidetes phylum and ubiquitous in healthy human GI tract flora. Anaerobic heme metabolism by commensal Bacteroidetes via hmu is likely a major contributor to human host metabolism of the heme from dietary red meat and a driver for the selective growth of these species in the GI tract consortium. IMPORTANCE Research on bacterial iron metabolism has historically focused on the host-pathogen relationship, where the host suppresses pathogen growth by cutting off access to iron. Less is known about how host iron is shared with bacterial species that live commensally in the anaerobic human GI tract, typified by members of phylum Bacteroidetes. While many facultative pathogens avidly produce and consume heme iron, most GI tract anaerobes are heme auxotrophs whose metabolic preferences we aimed to describe. Understanding iron metabolism by model microbiome species like Bacteroides thetaiotaomicron is essential for modeling the ecology of the GI tract, which serves the long-term biomedical goals of manipulating the microbiome to facilitate host metabolism of iron and remediate dysbiosis and associated pathologies (e.g., inflammation and cancer).