Myxobacterial Genomics and Post-Genomics: A Review of Genome Biology, Genome Sequences and Related 'Omics Studies

Long-read sequencing sheds light on key bacteria contributing to deadwood decomposition processes

BACKGROUND

Deadwood decomposition is an essential ecological process in forest ecosystems, playing a key role in nutrient cycling and carbon sequestration by enriching soils with organic matter. This process is driven by diverse microbial communities encompassing specialized functions in breaking down organic matter, but the specific roles of individual microorganisms in this process are still not fully understood.

RESULTS

Here, we characterized the deadwood microbiome in a natural mixed temperate forest in Central Europe using PacBio HiFi long-read sequencing and a genome-resolved transcriptomics approach in order to uncover key microbial contributors to wood decomposition. We obtained high quality assemblies, which allowed attribution of complex microbial functions such as nitrogen fixation to individual microbial taxa and enabled the recovery of metagenome-assembled genomes (MAGs) from both abundant and rare deadwood bacteria. We successfully assembled 69 MAGs (including 14 high-quality and 7 single-contig genomes) from 4 samples, representing most of the abundant bacterial phyla in deadwood. The MAGs exhibited a rich diversity of carbohydrate-active enzymes (CAZymes), with Myxococcota encoding the highest number of CAZymes and the full complement of enzymes required for cellulose decomposition. For the first time we observed active nitrogen fixation by Steroidobacteraceae, as well as hemicellulose degradation and chitin recycling by Patescibacteria. Furthermore, PacBio HiFi sequencing identified over 1000 biosynthetic gene clusters, highlighting a vast potential for secondary metabolite production in deadwood, particularly in Pseudomonadota and Myxococcota.

CONCLUSIONS

PacBio HiFi long-read sequencing offers comprehensive insights into deadwood decomposition processes by advancing the identification of functional features involving multiple genes. It represents a robust tool for unraveling novel microbial genomes in complex ecosystems and allows the identification of key microorganisms contributing to deadwood decomposition.

Heterologous biosynthesis of myxobacterial lanthipeptides melittapeptins

The myxobacteria are an attractive bioresource for bioactive compounds since the large size genome contains many biosynthetic gene clusters of secondary metabolites. The genome of the myxobacterium Melittangium boletus contains three biosynthetic gene clusters for lanthipeptide production. One of the gene clusters includes genes coding lanthipeptide precursor (melA), class II lanthipeptide synthetase (melM), and transporter (melT). The amino acid sequence of melA indicated similarity with that of known lanthipeptides mersacidin and lichenicidin A1 by the alignment. To perform heterologous production of new lanthipeptides, the expression vector containing the essential genes (melA and melM) was constructed by utilizing codon-optimized synthetic genes. The co-expression of two genes in the host bacterial cells of Escherichia coli BL21 (DE3) afforded new lanthipeptides named melittapeptins A-C. The structures of melittapeptins A-C including lanthionine/methyllanthionine bridge pattern were proposed based on protease digestion and MS/MS experiments. The native strain of M. boletus did not produce melittapeptins A-C, so heterologous production using the biosynthetic gene cluster was effective in obtaining the lanthipeptides. Melittapeptins A-C showed specific and potent antibacterial activity to the Gram-positive bacterium Micrococcus luteus. To the best of our knowledge, this is the first report of antibacterial lanthipeptides derived from myxobacterial origin. KEY POINTS: • New lanthipeptides melittapeptins were heterologously produced in Escherichia coli. • Melittapeptins showed specific antibacterial activity against Micrococcus luteus. • Melittapeptins were the first antibacterial lanthipeptides of myxobacterial origin.

Combining culture optimization and synthetic biology to improve production and detection of secondary metabolites in Myxococcus xanthus: application to myxoprincomide

Microbial secondary metabolites play crucial ecological roles in governing species interactions and contributing to their defense strategies. Their unique structures and potent bioactivities have been key in discovering antibiotics and anticancer drugs. Genome sequencing has undoubtedly revealed that myxobacteria constitute a huge reservoir of secondary metabolites as the well-known producers, actinomycetes. However, because most secondary metabolites are not produced in the laboratory context, the natural products from myxobacteria characterized to date represent only the tip of the iceberg. By combining the engineering of a dedicated Myxococcus xanthus DZ2 chassis strain with a two-step growth medium protocol, we provide a new approach called two-step Protocol for Resource Integration and Maximization-Biomolecules Overproduction and Optimal Screening Therapeutics (2PRIM-BOOST) for the production of non-ribosomal peptides synthetases (NRPS)/polyketides synthases (PKS) secondary metabolites from myxobacteria. We further show that the 2PRIM-BOOST strategy will facilitate the screening of secondary metabolites for biological activities of medical interest. As proof of concept, using a constitutive strong promoter, the myxoprincomide from M. xanthus DZ2 has been efficiently produced and its biosynthesis has been enhanced using the 2PRIM-BOOST approach, allowing the identification of new features of myxoprincomide. This strategy should allow the chances to produce and discover new NRPS, PKS, and mixed NRPS/PKS hybrid natural metabolites that are currently considered as cryptic and are the most represented in myxobacteria.IMPORTANCEMicrobial secondary metabolites are important in species interactions and are also a prolific source of drugs. Myxobacteria are ubiquitous soil-dwelling bacteria constituting a huge reservoir of secondary metabolites. However, because most of these molecules are not produced in the laboratory context, one can estimate that only one-tenth have been characterized to date. Here, we developed a new strategy called two-step Protocol for Resource Integration and Maximization-Biomolecules Overproduction and Optimal Screening Therapeutics (2PRIM-BOOST) that combines the engineering of a dedicated Myxococcus xanthus chassis strain together with growth medium optimization. By combining these strategies with the insertion of a constitutive promoter upstream the biosynthetic gene cluster (BGC), the production of myxoprincomide, a characterized low-produced secondary metabolite, was successfully and significantly increased. The 2PRIM-BOOST enriches the toolbox used to produce previously cryptic metabolites, unveil their ecological role, and provide new molecules of medical interest.

Genome and pan-genome analysis of a new exopolysaccharide-producing bacterium Pyschrobacillus sp. isolated from iron ores deposit and insights into iron uptake

Bacterial exopolysaccharides (EPS) have emerged as one of the key players in the field of heavy metal-contaminated environmental bioremediation. This study aimed to characterize and evaluate the metal biosorption potential of EPS produced by a novel Psychrobacillus strain, NEAU-3TGS, isolated from an iron ore deposit at Tamra iron mine, northern Tunisia. Genomic and pan-genomic analysis of NEAU-3TGS bacterium with nine validated published Psychrobacillus species was also performed. The results showed that the NEAU-3TGS genome (4.48 Mb) had a mean GC content of 36%, 4,243 coding sequences and 14 RNA genes. Phylogenomic analysis and calculation of nucleotide identity (ANI) values (less than 95% for new species with all strains) confirmed that NEAU-3TGS represents a potential new species. Pangenomic analysis revealed that Psychrobacillus genomic diversity represents an "open" pangenome model with 33,091 homologous genes, including 65 core, 3,738 shell, and 29,288 cloud genes. Structural EPS characterization by attenuated total reflectance-Fourier transform infrared (ATR-FTIR) spectroscopy showed uronic acid and α-1,4-glycosidic bonds as dominant components of the EPS. X-ray diffraction (XRD) analysis revealed the presence of chitin, chitosan, and calcite CaCO3 and confirmed the amorphous nature of the EPS. Heavy metal bioabsorption assessment showed that iron and lead were more adsorbed than copper and cadmium. Notably, the optimum activity was observed at 37°C, pH=7 and after 3 h contact of EPS with each metal. Genomic insights on iron acquisition and metabolism in Psychrobacillus sp. NEAU-3TGS suggested that no genes involved in siderophore biosynthesis were found, and only the gene cluster FeuABCD and trilactone hydrolase genes involved in the uptake of siderophores, iron transporter and exporter are present. Molecular modelling and docking of FeuA (protein peptidoglycan siderophore-binding protein) and siderophores ferrienterobactine [Fe+3 (ENT)]-3 and ferribacillibactine [Fe+3 (BB)]-3 ligand revealed that [Fe+3 (ENT)]-3 binds to Phe122, Lys127, Ile100, Gln314, Arg215, Arg217, and Gln252. Almost the same for [Fe+3 (ENT)]-3 in addition to Cys222 and Tyr229, but not Ile100.To the best of our knowledge, this is the first report on the characterization of EPS and the adsorption of heavy metals by Psychrobacillus species. The heavy metal removal capabilities may be advantageous for using these organisms in metal remediation.

MyxoPortal: a database of myxobacterial genomic features

Myxobacteria are predatory bacteria with antimicrobial activity, utilizing complex mechanisms to kill their prey and assimilate their macromolecules. Having large genomes encoding hundreds of secondary metabolites, hydrolytic enzymes and antimicrobial peptides, these organisms are widely studied for their antibiotic potential. MyxoPortal is a comprehensive genomic database hosting 262 genomes of myxobacterial strains. Datasets included provide genome annotations with gene locations, functions, amino acids and nucleotide sequences, allowing analysis of evolutionary and taxonomical relationships between strains and genes. Biosynthetic gene clusters are identified by AntiSMASH, and dbAMP-generated antimicrobial peptide sequences are included as a resource for novel antimicrobial discoveries, while curated datasets of CRISPR/Cas genes, regulatory protein sequences, and phage associated genes give useful insights into each strain's biological properties. MyxoPortal is an intuitive open-source database that brings together application-oriented genomic features that can be used in taxonomy, evolution, predation and antimicrobial research. MyxoPortal can be accessed at http://dicsoft1.physics.iisc.ac.in/MyxoPortal/. Database URL:  http://dicsoft1.physics.iisc.ac.in/MyxoPortal/. Graphical Abstract.

The genetic basis of predation by myxobacteria

Myxobacteria (phylum Myxococcota) are abundant and virtually ubiquitous microbial predators. Facultatively multicellular organisms, they are able to form multicellular fruiting bodies and swarm across surfaces, cooperatively hunting for prey. Myxobacterial communities are able to kill a wide range of prey microbes, assimilating their biomass to fuel population growth. Their mechanism of predation is exobiotic - hydrolytic enzymes and toxic metabolites are secreted into the extracellular environment, killing and digesting prey cells from without. However, recent observations of single-cell predation and contact-dependent prey killing challenge the dogma of myxobacterial predation being obligately cooperative. Regardless of their predatory mechanisms, myxobacteria have a broad prey range, which includes Gram-negative bacteria, Gram-positive bacteria and fungi. Pangenome analyses have shown that their extremely large genomes are mainly composed of accessory genes, which are not shared by all members of their species. It seems that the diversity of accessory genes in different strains provides the breadth of activity required to prey upon such a smorgasbord of microbes, and also explains the considerable strain-to-strain variation in predatory efficiency against specific prey. After providing a short introduction to general features of myxobacterial biology which are relevant to predation, this review brings together a rapidly growing body of work into the molecular mechanisms and genetic basis of predation, presenting a summary of current knowledge, highlighting trends in research and suggesting strategies by which we can potentially exploit myxobacterial predation in the future.

The RNA cargo of Myxococcus outer membrane vesicles

The outer membrane vesicles (OMVs) secreted by some Gram-negative bacteria contain RNA cargo, which can be introduced into target cells, affecting their cellular processes. To test whether the antimicrobial OMVs secreted by predatory myxobacteria might contain cargo RNA with a role in prey killing, we purified OMVs and cells from four different strains of Myxococcus spp. for RNA-seq transcriptome sequencing. Myxobacterial OMVs contained distinct sets of RNA molecules. The abundance of major cellular transcripts correlated strongly with their abundance in OMVs, suggesting non-specific packaging into OMVs. However, many major cellular transcripts were absent entirely from OMVs and some transcripts were found exclusively in OMVs, suggesting OMV RNA cargo loading is not simply a consequence of sampling the cellular transcriptome. Despite considerable variation in OMV RNA cargo between biological replicates, a small number of transcripts were found consistently in replicate OMV preparations. These 'core' OMV transcripts were often found in the OMVs from multiple strains, and sometimes enriched relative to their abundance in cellular transcriptomes. In addition to providing the first transcriptomes for myxobacterial OMVs, and the first cellular transcriptomes for three strains of Myxococcus spp., we highlight five transcripts for further study. These transcripts are 'core' for at least two of the three strains of M. xanthus studied, and encode two alkyl hydroperoxidase proteins (AhpC and AhpD), two ribosome-associated inhibitors (RaiA-like) and a DO-family protease. It will be interesting to test whether the transcripts serve a biological function within OMVs, potentially being transported into prey cells for translation into toxic proteins.

Bacterial Diversity in Old Hydrocarbon Polluted Sediments of Ecuadorian Amazon River Basins

The Ecuadorian Amazon rainforest stands out as one of the world's most biodiverse regions, yet faces significant threats due to oil extraction activities dating back to the 1970s in the northeastern provinces. This research investigates the environmental and societal consequences of prolonged petroleum exploitation and oil spills in Ecuador's Amazon. Conducted in June 2015, the study involved a comprehensive analysis of freshwater sediment samples from 24 locations in the Rio Aguarico and Napo basins. Parameters such as water and air temperature, conductivity, soil pH, and hydrocarbon concentrations were examined. Total petroleum hydrocarbon (TPH) concentrations ranged from 9.4 to 847.4 mg kg-1, with polycyclic aromatic hydrocarbon (PAH) levels varying from 10.15 to 711.1 mg kg-1. The pristane/phytane ratio indicated historic hydrocarbon pollution in 8 of the 15 chemically analyzed sediments. Using non-culturable techniques (Illumina), bacterial analyses identified over 350 ASV, with prominent families including Comamonadaceae, Chitinophagaceae, Anaeromyxobacteraceae, Sphingomonadaceae, and Xanthobacteraceae. Bacterial diversity, assessed in eight samples, exhibited a positive correlation with PAH concentrations. The study provides insights into how microbial communities respond to varying levels of hydrocarbon pollution, shedding light on the enduring impact of oil exploitation in the Amazonian region. Its objective is to deepen our understanding of the environmental and human well-being in the affected area, underscoring the pressing need for remedial actions in the face of ongoing ecological challenges.

Identification of secondary metabolites containing a diketopiperazine core in extracts from myxobacterial strains with growth inhibition activity against a range of prey species

Myxobacteria produce a variety of bioactive secondary metabolites, and with a wealth of under-researched species they hold vast potential for undiscovered compounds. With the ever-increasing need for new antibiotics, the development of novel therapeutics is vitally important. Therefore, this study aimed to extract and elucidate antimicrobial metabolites from the following myxobacteria: Myxococcus xanthus CA010 and AB022; Corallococcus exiguus DSM 14696T; Myxococcus stipitatus DSM 14675T; and Corallococcus aberystwythensis AB050AT. Metabolite mixtures were extracted in acetone from XAD-16 resin incubated in liquid cultures and analysed using GC-MS. Bioactivity was identified using a growth inhibition assay against a panel of clinically relevant prey species including Gram-positive and Gram-negative bacteria and a fungus. Growth of Klebsiella pneumoniae and Enterococcus faecalis was most affected by the metabolite mixtures and the mixtures from AB022 and AB050AT were effective against the most prey. GC-MS analysis revealed metabolites with roles in the synthesis and degradation of amino acids and fatty acids, but also identified compounds A and B with a diketopiperazine (DKP) core. With previously confirmed bioactivity of compound A, it is suggested that these DKP compounds are contributing to the antimicrobial activity observed. Furthermore, many compounds could not be identified and so these unknowns present further potential for novel bioactive compounds.

Improved electrochemical properties of graphite electrodes incubated with iron powders in rice-paddy fields boost power outputs from microbial fuel cells

Studies have shown that the supplementation of anode-surrounding soil with zero-valent iron (ZVI) boosts power outputs from rice paddy-field microbial fuel cells (RP-MFCs). In order to understand mechanisms by which ZVI boosts outputs from RP-MFCs, the present study operated RP-MFCs with and without ZVI, and compositions of anode-associated bacteria and electrochemical properties of graphite anodes were analyzed after 3-month operation. Metabarcoding using 16S rRNA gene fragments showed that bacterial compositions did not largely differ among these RP-MFCs. Cyclic voltammetry showed improved electrochemical properties of anodes recovered from ZVI-supplemented RP-MFCs, and this was attributed to the adhesion of iron-oxide films onto graphite surfaces. Bioelectrochemical devices equipped with graphite anodes recovered from ZVI-supplemented RP-MFCs generated higher currents than those with fresh graphite anodes. These results suggest that ZVI is oxidized to iron oxides in paddy-field soil and adheres onto graphite anodes, resulting in the boost of power outputs from RP-MFCs.

Phylogenetic Revisit to a Review on Predatory Bacteria

Predatory bacteria, along with the biology of their predatory behavior, have attracted interest in terms of their ecological significance and industrial applications, a trend that has been even more pronounced since the comprehensive review in 2016. This mini-review does not cover research trends, such as the role of outer membrane vesicles in myxobacterial predation, but provides an overview of the classification and newly described taxa of predatory bacteria since 2016, particularly with regard to phylogenetic aspects. Among them, it is noteworthy that in 2020 there was a major phylogenetic reorganization that the taxa hosting Bdellovibrio and Myxococcus, formerly classified as Deltaproteobacteria, were proposed as the new phyla Bdellovibrionota and Myxococcota, respectively. Predatory bacteria have been reported from other phyla, especially from the candidate divisions. Predatory bacteria that prey on cyanobacteria and predatory cyanobacteria that prey on Chlorella have also been found. These are also covered in this mini-review, and trans-phylum phylogenetic trees are presented.

Whole-genome sequencing of a biocontrol Myxococcus xanthus R31 isolate and comparative genomic analysis

Tomato bacterial wilt (TBW) caused by Ralstonia solanacearum is one of the most destructive soil-borne diseases. Myxococcus xanthus R31, isolated from healthy tomato rhizosphere soil using the R. solanacearum baiting method, exhibiting good biocontrol efficacy against TBW. However, the genomic information and evolutionary features of R31 are largely unclear. Here, the high-quality genome assembly of R31 was presented. Using Nanopore sequencing technology, we assembled the 9.25 Mb complete genome of R31 and identified several extracellular enzyme proteins, including carbohydrate-active enzymes (CAZymes) and peptidases. We also performed a comparative genome analysis of R31 and 17 other strains of M. xanthus with genome sequences in the NCBI database to gain insights into myxobacteria predation and genome size expansion. Average nucleotide identity and digital DNA-DNA hybridization calculation and phylogenetic analysis indicated that R31 was closely related to the species M. xanthus. Further comparative genomics analysis suggested that, in addition to characteristics of predatory microorganisms, R31 contains many strain-specific genes, which may provide a genetic basis for its proficient predatory ability. This study provides new insights into R31 and other closely related species and facilitates studies using genetic approaches to further elucidate the predation mechanism of myxobacteria.

Predatory Strategies of Myxococcus xanthus: Prey Susceptibility to OMVs and Moonlighting Enzymes

Predatory outer membrane vesicles (OMVs) secreted by myxobacteria fuse readily with the outer membranes of Gram-negative bacteria, introducing toxic cargo into their prey. Here we used a strain of the myxobacterium Myxococcus xanthus that produces fluorescent OMVs to assay the uptake of OMVs by a panel of Gram-negative bacteria. M. xanthus strains took up significantly less OMV material than the tested prey strains, suggesting that re-fusion of OMVs with producing organisms is somehow inhibited. The OMV killing activity against different prey correlated strongly with the predatory activity of myxobacterial cells, however, there was no correlation between OMV killing activity and their propensity to fuse with different prey. It has previously been proposed that M. xanthus GAPDH stimulates the predatory activity of OMVs by enhancing OMV fusion with prey cells. Therefore, we expressed and purified active fusion proteins of M. xanthus glyceraldehyde-3-phosphate dehydrogenase and phosphoglycerate kinase (GAPDH and PGK; moonlighting enzymes with additional activities beyond their roles in glycolysis/gluconeogenesis) to investigate any involvement in OMV-mediated predation. Neither GAPDH nor PGK caused lysis of prey cells or enhanced OMV-mediated lysis of prey cells. However, both enzymes were found to inhibit the growth of Escherichia coli, even in the absence of OMVs. Our results suggest that fusion efficiency is not a determinant of prey killing, but instead resistance to the cargo of OMVs and co-secreted enzymes dictates whether organisms can be preyed upon by myxobacteria.

Constructing a Myxobacterial Natural Product Database to Facilitate NMR-Based Metabolomics Bioprospecting of Myxobacteria

Myxobacteria are fascinating prokaryotes featuring a potent capacity for producing a wealth of bioactive molecules with intricate chemical topology as well as intriguing enzymology, and thus it is critical to developing an efficient pipeline for bioprospecting. Herein, we construct the database MyxoDB, the first public compendium solely dedicated to myxobacteria, which enabled us to provide an overview of the structural diversity and taxonomic distribution of known myxobacterial natural products. Moreover, we demonstrated that the cutting-edge NMR-based metabolomics was effective to differentiate the biosynthetic priority of myxobacteria, whereby MyxoDB could greatly streamline the dereplication of multifarious known compounds and accordingly speed up the discovery of new compounds. This led to the rapid identification of a class of linear di-lipopeptides (archangimins) and a rare rearranged sterol (corasterol) that were endowed with unique chemical architectures and/or biosynthetic enzymology. We also showcased that NMR-based metabolomics, MyxoDB, and genomics can also work concertedly to accelerate the targeted discovery of a polyketidic compound pyxipyrrolone C. All in all, this study sets the stage for the discovery of many more novel natural products from underexplored myxobacterial resources.

Myxobacterial Predation: A Standardised Lawn Predation Assay Highlights Strains with Unusually Efficient Predatory Activity

Myxobacteria prey upon a broad range of microorganisms. Lawn assays are commonly used to quantify myxobacterial predation-myxobacterial suspensions are spotted onto prey lawns, and monitored via spot expansion. The diversity in motility behaviours of myxobacterial strains and differing assay protocols in myxobacteriology laboratories led us to develop a highly-specified assay, which was applied to 28 myxobacterial strains preying on seven phytopathogenic prey species. Generally, prey organisms showed no qualitative differences in their susceptibility/resistance to myxobacterial predation. For most myxobacteria, prey did not stimulate, and in ~50% of cases actively hindered colony expansion. Only ~25% of predator/prey strain combinations exhibited greater colony expansion than in the absence of nutrients. The activity of predatory strains against different prey correlated, implying effective predators may have relatively non-specific predation mechanisms (e.g., broad specificity proteases/lipases), but no correlation was observed between predatory activity and phylogeny. Predation on dead (but intact) or lysed prey cells gave greater colony expansion than on live prey. Occasional strains grew substantially faster on dead compared to lysed cells, or vice-versa. Such differences in accessing nutrients from live, dead and lysed cells indicates there are strain-specific differences in the efficiencies/machineries of prey killing and nutrient acquisition, which has important implications for the ecology of myxobacterial predators and their prey.

In Silico and In Vitro Analyses Reveal Promising Antimicrobial Peptides from Myxobacteria

Antimicrobial resistance (AMR) is a global concern, and as soon as new antibiotics are introduced, resistance to those agents emerges. Therefore, there is an increased appetite for alternative antimicrobial agents to traditional antibiotics. Here, we used in silico methods to investigate potential antimicrobial peptides (AMPs) from predatory myxobacteria. Six hundred seventy-two potential AMP sequences were extracted from eight complete myxobacterial genomes. Most putative AMPs were predicted to be active against Klebsiella pneumoniae with least activity being predicted against Staphylococcus aureus. One hundred seventeen AMPs (defined here as 'potent putative AMPs') were predicted to have very good activity against more than two bacterial pathogens, and these were characterized further in silico. All potent putative AMPs were predicted to have anti-inflammatory and antifungal properties, but none was predicted to be active against viruses. Twenty six (22%) of them were predicted to be hemolytic to human erythrocytes, five were predicted to have anticancer properties, and 56 (47%) were predicted to be biofilm active. In vitro assays using four synthesized AMPs showed high MIC values (e.g. So_ce_56_913 250 µg/ml and Coral_AMP411 125 µg/ml against E. coli). However, antibiofilm assays showed a substantial reduction in numbers (e.g. Coral_AMP411 and Myxo_mac104 showed a 69% and 73% reduction, respectively, at the lowest concentration against E. coli) compared to traditional antibiotics. Fourteen putative AMPs had high sequence similarity to proteins which were functionally associated with proteins of known function. The myxobacterial genomes also possessed a variety of biosynthetic gene clusters (BGCs) that can encode antimicrobial secondary metabolites, but their numbers did not correlate with those of the AMPs. We suggest that AMPs from myxobacteria are a promising source of novel antimicrobial agents with a plethora of biological properties.