A Genomic Survey of Signalling in the Myxococcaceae

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.

Lysis profile and preference of Myxococcus sp. PT13 for typical soil bacteria

Introduction

Myxococcus sp. PT13 is a wild strain with multiple predatory properties that prey on multiple model microorganisms preserved in the laboratory. However, the lysis spectrum of PT13 on typical soil bacteria and its driving effect on soil microecosystems are still unclear.

Methods

In this study, the lawn predation method was used to determine the predation diameter of 62 typical soil bacteria by myxobacteria PT13 and analyze their lysis spectra.

Results and Discussion

The results showed that PT13 had a predation diameter greater than 15 mm against typical soil microorganisms such as Aeromonas, Bacillus, Brevibacterium, Fictibacillus, Glutamicibacter, Herbaspirillum, and Leifsonia and had an outstanding lysis effect but a significant preference (p < 0.05). Absolute high-throughput sequencing results showed that PT13 predation drove the microcosmic system composed of 16 bacterial genera, with a significant decrease in the Shannon index by 11.8% (CK = 2.04, D = 1.80) and a significant increase in the Simpson index by 45.0% (CK = 0.20, D = 0.29). The results of principal coordinate analysis (PCoA) showed that myxobacterial addition significantly disturbed the microcosmic microbial community structure (ANOSIM, p < 0.05). LEfSe analysis showed that the relative and absolute abundances (copy numbers) of Bacillus, Pedobacter, Staphylococcus, Streptomyces and Fictibacillus decreased significantly very likely due to myxobacterial predation (p < 0.05). However, the predatory effect of PT13 also increased the relative or absolute abundances of some species, such as Sphingobacterium, Paenarthrobacter, Microbacterium, and Leifsonia. It can be concluded that PT13 has a broad-spectrum lysis spectrum but poor cleavage ability for Streptomyces, and the interaction between complex microorganisms limits the predation effect of PT13 on some prey bacteria. This in turn allows some prey to coexist with myxobacteria. This paper will lay a theoretical foundation for the regulation of soil microecology dominated by myxobacteria.

Genomes of Novel Myxococcota Reveal Severely Curtailed Machineries for Predation and Cellular Differentiation

Cultured Myxococcota are predominantly aerobic soil inhabitants, characterized by their highly coordinated predation and cellular differentiation capacities. Little is currently known regarding yet-uncultured Myxococcota from anaerobic, nonsoil habitats. We analyzed genomes representing one novel order (o__JAFGXQ01) and one novel family (f__JAFGIB01) in the Myxococcota from an anoxic freshwater spring (Zodletone Spring) in Oklahoma, USA. Compared to their soil counterparts, anaerobic Myxococcota possess smaller genomes and a smaller number of genes encoding biosynthetic gene clusters (BGCs), peptidases, one- and two-component signal transduction systems, and transcriptional regulators. Detailed analysis of 13 distinct pathways/processes crucial to predation and cellular differentiation revealed severely curtailed machineries, with the notable absence of homologs for key transcription factors (e.g., FruA and MrpC), outer membrane exchange receptor (TraA), and the majority of sporulation-specific and A-motility-specific genes. Further, machine learning approaches based on a set of 634 genes informative of social lifestyle predicted a nonsocial behavior for Zodletone Myxococcota. Metabolically, Zodletone Myxococcota genomes lacked aerobic respiratory capacities but carried genes suggestive of fermentation, dissimilatory nitrite reduction, and dissimilatory sulfate-reduction (in f_JAFGIB01) for energy acquisition. We propose that predation and cellular differentiation represent a niche adaptation strategy that evolved circa 500 million years ago (Mya) in response to the rise of soil as a distinct habitat on Earth. IMPORTANCE The phylum Myxococcota is a phylogenetically coherent bacterial lineage that exhibits unique social traits. Cultured Myxococcota are predominantly aerobic soil-dwelling microorganisms that are capable of predation and fruiting body formation. However, multiple yet-uncultured lineages within the Myxococcota have been encountered in a wide range of nonsoil, predominantly anaerobic habitats, and the metabolic capabilities, physiological preferences, and capacity of social behavior of such lineages remain unclear. Here, we analyzed genomes recovered from a metagenomic analysis of an anoxic freshwater spring in Oklahoma, USA, that represent novel, yet-uncultured, orders and families in the Myxococcota. The genomes appear to lack the characteristic hallmarks for social behavior encountered in Myxococcota genomes and displayed a significantly smaller genome size and a smaller number of genes encoding biosynthetic gene clusters, peptidases, signal transduction systems, and transcriptional regulators. Such perceived lack of social capacity was confirmed through detailed comparative genomic analysis of 13 pathways associated with Myxococcota social behavior, as well as the implementation of machine learning approaches to predict social behavior based on genome composition. Metabolically, these novel Myxococcota are predicted to be strict anaerobes, utilizing fermentation, nitrate reduction, and dissimilarity sulfate reduction for energy acquisition. Our results highlight the broad patterns of metabolic diversity within the yet-uncultured Myxococcota and suggest that the evolution of predation and fruiting body formation in the Myxococcota has occurred in response to soil formation as a distinct habitat on Earth.

Phylogenomic Analysis Substantiates the gyrB Gene as a Powerful Molecular Marker to Efficiently Differentiate the Most Closely Related Genera Myxococcus, Corallococcus, and Pyxidicoccus

Rapid and accurate strain identification of the most closely related genera Myxococcus, Corallococcus, and Pyxidicoccus can enhance the efficiency of the mining of novel secondary metabolites through dereplication. However, the commonly used 16S rRNA gene sequencing cannot accurately differentiate members of the three genera above, and the whole-genome sequencing is unable to rapidly and inexpensively provide species assignation toward a large number of isolates. To overcome the limitations, the gyrB gene was investigated as a candidate genetic marker for exploring the phylogenetic relationships of bacteria within the three genera and for developing the gyrB-based typing method. Here, the bacterial phylogeny and species affiliations of the three genera were determined based on the phylogenomic reconstruction and the analysis of digital DNA–DNA hybridization values among 90 genomes, further confirming nine novel taxa and assigning over one-third of genomes to defined species. The phylogenetic relationships of these strains based on the gyrB gene sequences were congruent with those based on their genome sequences, allowing the use of the gyrB gene as a molecular marker. The gyrB gene-specific primers for the PCR-amplification and sequencing of bacteria within the three genera were designed and validated for 31 isolates from our group collection. The gyrB-based taxonomic tool proved to be able to differentiate closely related isolates at the species level. Based on the newly proposed 98.6% identity threshold for the 966-bp gyrB gene and the phylogenetic inference, these isolates were assigned into two known species and eight additional putative new species. In summary, this report demonstrated that the gyrB gene is a powerful phylogenetic marker for taxonomy and phylogeny of bacteria within the closely related genera Myxococcus, Corallococcus, and Pyxidicoccus, particularly in the case of hundreds or thousands of isolates in environmental studies.

Differential response to prey quorum signals indicates predatory specialization of myxobacteria and ability to predate Pseudomonas aeruginosa

Multiomic analysis of transcriptional and metabolic responses from the predatory myxobacteria Myxococcus xanthus and Cystobacter ferrugineus exposed to prey signalling molecules of the acylhomoserine lactone and quinolone quorum signalling classes provided insight into predatory specialization. Acylhomoserine lactone quorum signals elicited a general response from both myxobacteria. We suggest that this is likely due to the generalist predator lifestyles of myxobacteria and ubiquity of acylhomoserine lactone signals. We also provide data that indicates the core homoserine lactone moiety included in all acylhomoserine lactone scaffolds to be sufficient to induce this general response. Comparing both myxobacteria, unique transcriptional and metabolic responses were observed from Cystobacter ferrugineus exposed to the quinolone signal 2‐heptylquinolin‐4(1H)‐one (HHQ) natively produced by Pseudomonas aeruginosa. We suggest that this unique response and ability to metabolize quinolone signals contribute to the superior predation of P. aeruginosa observed from C. ferrugineus. These results further demonstrate myxobacterial eavesdropping on prey signalling molecules and provide insight into how responses to exogenous signals might correlate with prey range of myxobacteria.

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

Myxobacteria are fascinating and complex microbes. They prey upon other members of the soil microbiome by secreting antimicrobial proteins and metabolites, and will undergo multicellular development if starved. The genome sequence of the model myxobacterium Myxococcus xanthus DK1622 was published in 2006 and 15 years later, 163 myxobacterial genome sequences have now been made public. This explosion in genomic data has enabled comparative genomics analyses to be performed across the taxon, providing important insights into myxobacterial gene conservation and evolution. The availability of myxobacterial genome sequences has allowed system-wide functional genomic investigations into entire classes of genes. It has also enabled post-genomic technologies to be applied to myxobacteria, including transcriptome analyses (microarrays and RNA-seq), proteome studies (gel-based and gel-free), investigations into protein–DNA interactions (ChIP-seq) and metabolism. Here, we review myxobacterial genome sequencing, and summarise the insights into myxobacterial biology that have emerged as a result. We also outline the application of functional genomics and post-genomic approaches in myxobacterial research, highlighting important findings to emerge from seminal studies. The review also provides a comprehensive guide to the genomic datasets available in mid-2021 for myxobacteria (including 24 genomes that we have sequenced and which are described here for the first time).

Whole-Genome Sequencing of Corallococcus sp. Strain EGB Reveals the Genetic Determinants Linking Taxonomy and Predatory Behavior

Corallococcus sp. strain EGB is a Gram-negative myxobacteria isolated from saline soil, and has considerable potential for the biocontrol of phytopathogenic fungi. However, the detailed mechanisms related to development and predatory behavior are unclear. To obtain a comprehensive overview of genetic features, the genome of strain EGB was sequenced, annotated, and compared with 10 other Corallococcus species. The strain EGB genome was assembled as a single circular chromosome of 9.4 Mb with 7916 coding genes. Phylogenomics analysis showed that strain EGB was most closely related to Corallococcus interemptor AB047A, and it was inferred to be a novel species within the Corallococcus genus. Comparative genomic analysis revealed that the pan-genome of Corallococcus genus was large and open. Only a small proportion of genes were specific to strain EGB, and most of them were annotated as hypothetical proteins. Subsequent analyses showed that strain EGB produced abundant extracellular enzymes such as chitinases and β-(1,3)-glucanases, and proteases to degrade the cell-wall components of phytopathogenic fungi. In addition, 35 biosynthetic gene clusters potentially coding for antimicrobial compounds were identified in the strain EGB, and the majority of them were present in the dispensable pan-genome with unexplored metabolites. Other genes related to secretion and regulation were also explored for strain EGB. This study opens new perspectives in the greater understanding of the predatory behavior of strain EGB, and facilitates a potential application in the biocontrol of fungal plant diseases in the future.

Identification of a solo acylhomoserine lactone synthase from the myxobacterium Archangium gephyra

Considered a key taxon in soil and marine microbial communities, myxobacteria exist as coordinated swarms that utilize a combination of lytic enzymes and specialized metabolites to facilitate predation of microbes. This capacity to produce specialized metabolites and the associated abundance of biosynthetic pathways contained within their genomes have motivated continued drug discovery efforts from myxobacteria. Of all myxobacterial biosynthetic gene clusters deposited in the antiSMASH database, only one putative acylhomoserine lactone (AHL) synthase, agpI, was observed, in genome data from Archangium gephyra. Without an AHL receptor also apparent in the genome of A. gephyra, we sought to determine if AgpI was an uncommon example of an orphaned AHL synthase. Herein we report the bioinformatic assessment of AgpI and discovery of a second AHL synthase from Vitiosangium sp. During axenic cultivation conditions, no detectible AHL metabolites were observed in A. gephyra extracts. However, heterologous expression of each synthase in Escherichia coli provided detectible quantities of 3 AHL signals including 2 known AHLs, C8-AHL and C9-AHL. These results suggest that A. gephyra AHL production is dormant during axenic cultivation. The functional, orphaned AHL synthase, AgpI, is unique to A. gephyra, and its utility to the predatory myxobacterium remains unknown.