Genome-wide analysis of myxobacterial two-component systems: genome relatedness and evolutionary changes

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.

Myxobacteria: biology and bioactive secondary metabolites

Myxobacteria are Gram-negative eubacteria and they thrive in a variety of habitats including soil rich in organic matter, rotting wood, animal dung and marine environment. Myxobacteria are a promising source of new compounds associated with diverse bioactive spectrum and unique mode of action. The genome information of myxobacteria has revealed many orphan biosynthetic pathways indicating that these bacteria can be the source of several novel natural products. In this review, we highlight the biology of myxobacteria with emphasis on their habitat, life cycle, isolation methods and enlist all the bioactive secondary metabolites purified till date and their mode of action.

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).

Assessment of Evolutionary Relationships for Prioritization of Myxobacteria for Natural Product Discovery

Discoveries of novel myxobacteria have started to unveil the potentially vast phylogenetic diversity within the family Myxococcaceae and have brought about an updated approach to myxobacterial classification. While traditional approaches focused on morphology, 16S gene sequences, and biochemistry, modern methods including comparative genomics have provided a more thorough assessment of myxobacterial taxonomy. Herein, we utilize long-read genome sequencing for two myxobacteria previously classified as Archangium primigenium and Chondrococcus macrosporus, as well as four environmental myxobacteria newly isolated for this study. Average nucleotide identity and digital DNA-DNA hybridization scores from comparative genomics suggest previously classified as A. primigenium to instead be a novel member of the genus Melittangium, C. macrosporus to be a potentially novel member of the genus Corallococcus with high similarity to Corallococcus exercitus, and the four isolated myxobacteria to include another novel Corallococcus species, a novel Pyxidicoccus species, a strain of Corallococcus exiguus, and a potentially novel Myxococcus species with high similarity to Myxococcus stipitatus. We assess the biosynthetic potential of each sequenced myxobacterium and suggest that genus-level conservation of biosynthetic pathways support our preliminary taxonomic assignment. Altogether, we suggest that long-read genome sequencing benefits the classification of myxobacteria and improves determination of biosynthetic potential for prioritization of natural product discovery.

Methods for Identifying Microbial Natural Product Compounds that Target Kinetoplastid RNA Structural Motifs by Homology and De Novo Modeled 18S rRNA

The development of novel anti-infectives against Kinetoplastids pathogens targeting proteins is a big problem occasioned by the antigenic variation in these parasites. This is also a global concern due to the zoonosis of these parasites, as they infect both humans and animals. Therefore, we need not only to create novel antibiotics, but also to speed up the development pipeline for these antibiotics. This may be achieved by using novel drug targets for Kinetoplastids drug discovery. In this study, we focused our attention on motifs of rRNA molecules that have been created using homology modeling. The RNA is the most ambiguous biopolymer in the kinetoplatid, which carries many different functions. For instance, tRNAs, rRNAs, and mRNAs are essential for gene expression both in the pro-and eukaryotes. However, all these types of RNAs have sequences with unique 3D structures that are specific for kinetoplastids only and can be used to shut down essential biochemical processes in kinetoplastids only. All these features make RNA very potent targets for antibacterial drug development. Here, we combine in silico methods combined with both computational biology and structure prediction tools to address our hypothesis. In this study, we outline a systematic approach for identifying kinetoplastid rRNA-ligand interactions and, more specifically, techniques that can be used to identify small molecules that target particular RNA. The high-resolution optimized model structures of these kineoplastids were generated using RNA 123, where all the stereochemical conflicts were solved and energies minimized to attain the best biological qualities. The high-resolution optimized model's structures of these kinetoplastids were generated using RNA 123 where all the stereochemical conflicts were solved and energies minimized to attain the best biological qualities. These models were further analyzed to give their docking assessment reliability. Docking strategies, virtual screening, and fishing approaches successfully recognized novel and myriad macromolecular targets for the myxobacterial natural products with high binding affinities to exploit the unmet therapeutic needs. We demonstrate a sensible exploitation of virtual screening strategies to 18S rRNA using natural products interfaced with classical maximization of their efficacy in phamacognosy strategies that are well established. Integration of these virtual screening strategies in natural products chemistry and biochemistry research will spur the development of potential interventions to these tropical neglected diseases.

Global gene expression analysis of the Myxococcus xanthus developmental time course

To accurately identify the genes and pathways involved in the initiation of the Myxococcus xanthus multicellular developmental program, we have previously reported a method of growing vegetative populations as biofilms within a controllable environment. Using a modified approach to remove up to ~90% rRNAs, we report a comprehensive transcriptional analysis of the M. xanthus developmental cycle while comparing it with the vegetative biofilms grown in rich and poor nutrients. This study identified 1522 differentially regulated genes distributed within eight clusters during development. It also provided a comprehensive overview of genes expressed during a nutrient-stress response, specific development time points, and during development initiation and regulation. We identified several differentially expressed genes involved in key central metabolic pathways suggesting their role in regulating myxobacterial development. Overall, this study will prove an important resource for myxobacterial researchers to delineate the regulatory and functional pathways responsible for development from those of the general nutrient stress response.

A Genomic Survey of Signalling in the Myxococcaceae

As prokaryotes diverge by evolution, essential 'core' genes required for conserved phenotypes are preferentially retained, while inessential 'accessory' genes are lost or diversify. We used the recently expanded number of myxobacterial genome sequences to investigate the conservation of their signalling proteins, focusing on two sister genera (Myxococcus and Corallococcus), and on a species within each genus (Myxococcus xanthus and Corallococcus exiguus). Four new C. exiguus genome sequences are also described here. Despite accessory genes accounting for substantial proportions of each myxobacterial genome, signalling proteins were found to be enriched in the core genome, with two-component system genes almost exclusively so. We also investigated the conservation of signalling proteins in three myxobacterial behaviours. The linear carotenogenesis pathway was entirely conserved, with no gene gain/loss observed. However, the modular fruiting body formation network was found to be evolutionarily plastic, with dispensable components in all modules (including components required for fruiting in the model myxobacterium M. xanthus DK1622). Quorum signalling (QS) is thought to be absent from most myxobacteria, however, they generally appear to be able to produce CAI-I (cholerae autoinducer-1), to sense other QS molecules, and to disrupt the QS of other organisms, potentially important abilities during predation of other prokaryotes.

An ambruticin-sensing complex modulates Myxococcus xanthus development and mediates myxobacterial interspecies communication

Starvation induces cell aggregation in the soil bacterium Myxococcus xanthus, followed by formation of fruiting bodies packed with myxospores. Sporulation in the absence of fruiting bodies can be artificially induced by high concentrations of glycerol through unclear mechanisms. Here, we show that a compound (ambruticin VS-3) produced by a different myxobacterium, Sorangium cellulosum, affects the development of M. xanthus in a similar manner. Both glycerol (at millimolar levels) and ambruticin VS-3 (at nanomolar concentrations) inhibit M. xanthus fruiting body formation under starvation, and induce sporulation in the presence of nutrients. The response is mediated in M. xanthus by three hybrid histidine kinases (AskA, AskB, AskC) that form complexes interacting with two major developmental regulators (MrpC, FruA). In addition, AskB binds directly to the mrpC promoter in vitro. Thus, our work indicates that the AskABC-dependent regulatory pathway mediates the responses to ambruticin VS-3 and glycerol. We hypothesize that production of ambruticin VS-3 may allow S. sorangium to outcompete M. xanthus under both starvation and growth conditions in soil.

Starvation induces cell aggregation and formation of spore-containing fruiting bodies in the bacterium Myxococcus xanthus. Here, the authors show that a different myxobacterial species produces a compound that inhibits the development of fruiting bodies in M. xanthus, by affecting the function of histidine kinases and major regulators.

A survey of non-coding RNAs in the social and predatory myxobacterium Myxococcus xanthus DK1622

Prokaryotic ncRNAs are important regulators of gene expression, and can be involved in complex signalling networks. The myxobacteria are model organisms for studies into multicellular development and microbial predation, being particularly renowned for their large genomes and exceptionally sophisticated signalling networks. However, apart from two specific examples, little is known about their regulatory ncRNAs. Here, we integrate bioinformatic predictions and transcriptome sequence data to provide a comprehensive survey of the ncRNAs made by the exemplar myxobacterium M. xanthus DK1622. M. xanthus RNA-seq data from four experimental conditions was interrogated to identify transcripts mapping outside coding sequences and to known ncRNAs. The resulting 37 ncRNAs were clustered on the genome and most (30/37) were conserved across the myxobacteria. A majority of ncRNAs (22/37) were intergenic, while 13 were at least partially antisense to protein-coding genes. Predicted promoter and terminator sequences explained the start/stop sites of 18 ncRNAs. mRNA targets for the ncRNAs were predicted, including plausible candidates for a known regulatory ncRNA. 22 ncRNAs were differentially expressed by nutrient availability and expression of 25 predicted targets was found to correlate strongly with that of their regulatory ncRNAs. Sharing of predicted mRNA targets by multiple ncRNAs suggests that some ncRNAs might regulate each other within signalling networks. This genomic survey of M. xanthus ncRNA biology provides a starting point for further studies of myxobacterial ncRNAs, which are likely to have important functions in these industrially important and sophisticated organisms.

Interspecies conflict affects RNA expression

Predation is an extreme form of competition between bacteria, involving the secretion of antimicrobial substances by predators, often packaged within outer membrane vesicles (OMVs). Recent studies into the Myxococcus xanthus/Escherichia coli predator/prey relationship have illuminated transcriptional changes during predation, identifying likely targets of predatory attack in the prey and nutrient assimilation strategies of the predator. Abundant non-coding RNAs can be observed in the predator and prey transcriptomes, with evidence of predation-dependent regulation of RNA levels. Given the observed secretion of regulatory RNAs within OMVs by bacteria, it will next be exciting to test whether the intercellular trafficking of regulatory RNAs is employed by predator and/or prey in their survival struggles.

Is "Wolf-Pack" Predation by Antimicrobial Bacteria Cooperative? Cell Behaviour and Predatory Mechanisms Indicate Profound Selfishness, Even when Working Alongside Kin

For decades, myxobacteria have been spotlighted as exemplars of social "wolf-pack" predation, communally secreting antimicrobial substances into the shared public milieu. This behavior has been described as cooperative, becoming more efficient if performed by more cells. However, laboratory evidence for cooperativity is limited and of little relevance to predation in a natural setting. In contrast, there is accumulating evidence for predatory mechanisms promoting "selfish" behavior during predation, which together with conflicting definitions of cooperativity, casts doubt on whether microbial "wolf-pack" predation really is cooperative. Here, it is hypothesized that public-goods-mediated predation is not cooperative, and it is argued that a holistic model of microbial predation is needed, accounting for predator and prey relatedness, social phenotypes, spatial organization, activity/specificity/transport of secreted toxins, and prey resistance mechanisms. Filling such gaps in our knowledge is vital if the evolutionary benefits of potentially costly microbial behaviors mediated by public goods are to be properly understood.

Genome Analysis, Metabolic Potential, and Predatory Capabilities of Herpetosiphon llansteffanense sp. nov

Predatory bacteria are able to kill and consume other microbes and are therefore of interest as potential sources of new antimicrobial substances for applications in the clinic. “Wolf pack” predators kill prey by secreting antimicrobial substances into their surroundings, and those substances can kill prey organisms independently of the predatory cells. The genus Herpetosiphon exhibits wolf pack predation, yet its members are poorly described compared to other wolf pack predators, such as the myxobacteria. By providing a thorough characterization of a novel Herpetosiphon species, including its predatory, biochemical, and genomic features, this study increases our understanding of genomic variation within the Herpetosiphon genus and how that variation affects predatory activity. This will facilitate future rational exploitation of genus members (and other wolf pack predators) as sources of novel antimicrobials.

Herpetosiphon spp. are ubiquitous, chemoheterotrophic, filamentous gliding bacteria with the ability to prey on other microbes through a “wolf pack” mechanism. The genus currently comprises four known species (H. aurantiacus, H. geysericola, H. giganteus, and H. gulosus), which produce antimicrobial secondary metabolites such as siphonazole. As part of a study isolating myxobacterial wolf pack predators, we serendipitously isolated a novel environmental strain (CA052B) from the edge of a stream at Llansteffan, United Kingdom, which was identified as a member of the Herpetosiphon genus. A lawn culture method was utilized to analyze the predatory activity of CA052B against 10 prey organisms of clinical relevance. CA052B was found to prey on Escherichia coli, Klebsiella pneumoniae, Proteus mirabilis, Staphylococcus aureus, Staphylococcus epidermidis, Staphylococcus saprophyticus, Enterococcus faecalis, Bacillus subtilis, and Candida albicans. Purified CA052B outer membrane vesicles also exhibited killing activity against the prey organisms when tested by flow cytometry. 16S rRNA sequencing of CA052B showed 98 to 99% identity with other Herpetosiphon species members. Comparing the genome of CA052B with the publicly available genomes of H. aurantiacus and H. geysericola revealed average nucleotide identities of only 84% and 91%, respectively, whereas the genome-to-genome distance calculation showed sequence identities of 28.2% and 46.6%, respectively. Biochemical characterization also revealed distinctions between CA052B and both H. gulosus and H. giganteus. Thus, strain CA052B^T (= DSM 107618^T = NBRC 113495^T) is proposed to be the type strain of a novel species, Herpetosiphon llansteffanense sp. nov. The genome sequence of CA052B also revealed diverse secondary metabolite biosynthetic clusters, encouraging further exploration of its antibiotic production potential.

IMPORTANCE Predatory bacteria are able to kill and consume other microbes and are therefore of interest as potential sources of new antimicrobial substances for applications in the clinic. “Wolf pack” predators kill prey by secreting antimicrobial substances into their surroundings, and those substances can kill prey organisms independently of the predatory cells. The genus Herpetosiphon exhibits wolf pack predation, yet its members are poorly described compared to other wolf pack predators, such as the myxobacteria. By providing a thorough characterization of a novel Herpetosiphon species, including its predatory, biochemical, and genomic features, this study increases our understanding of genomic variation within the Herpetosiphon genus and how that variation affects predatory activity. This will facilitate future rational exploitation of genus members (and other wolf pack predators) as sources of novel antimicrobials.

Diversity and Evolution of Myxobacterial Type IV Pilus Systems

Type IV pili (T4P) are surface-exposed protein fibers that play key roles in the bacterial life cycle via surface attachment/adhesion, biofilm formation, motility, and development. The order Myxococcales (myxobacteria) are members of the class Deltaproteobacteria and known for their large genome size and complex social behaviors, including gliding motility, fruiting body formation, biofilm production, and prey hunting. Myxococcus xanthus, the best-characterized member of the order, relies on the appropriate expression of 17 type IVa (T4aP) genes organized in a single cluster plus additional genes (distributed throughout the genome) for social motility and development. Here, we compared T4aP genes organization within the myxobacteria to understand their evolutionary origins and diversity. We found that T4aP genes are organized as large clusters in suborder Cystobacterineae, whereas in other two suborders Sorangiineae and Nannocystineae, these genes are dispersed throughout the genome. Based on the genomic organization, the phylogeny of conserved proteins, and synteny studies among 28 myxobacterial and 66 Proteobacterial genomes, we propose an evolutionary model for the origin of myxobacterial T4aP genes independently from other orders in class Deltaproteobacteria. Considering a major role for T4P, this study further proposes the origins and evolution of social motility in myxobacteria and provides a foundation for understanding how complex-behavioral traits, such as gliding motility, multicellular development, etc., might have evolved in this diverse group of complex organisms.

Transcriptional changes when Myxococcus xanthus preys on Escherichia coli suggest myxobacterial predators are constitutively toxic but regulate their feeding

Predation is a fundamental ecological process, but within most microbial ecosystems the molecular mechanisms of predation remain poorly understood. We investigated transcriptome changes associated with the predation of Escherichia coli by the myxobacterium Myxococcus xanthus using mRNA sequencing. Exposure to pre-killed prey significantly altered expression of 1319 predator genes. However, the transcriptional response to living prey was minimal, with only 12 genes being significantly up-regulated. The genes most induced by prey presence (kdpA and kdpB, members of the kdp regulon) were confirmed by reverse transcriptase quantitative PCR to be regulated by osmotic shock in M. xanthus, suggesting indirect sensing of prey. However, the prey showed extensive transcriptome changes when co-cultured with predator, with 40 % of its genes (1534) showing significant changes in expression. Bacteriolytic M. xanthus culture supernatant and secreted outer membrane vesicles (OMVs) also induced changes in expression of large numbers of prey genes (598 and 461, respectively). Five metabolic pathways were significantly enriched in prey genes up-regulated on exposure to OMVs, supernatant and/or predatory cells, including those for ribosome and lipopolysaccharide production, suggesting that the prey cell wall and protein production are primary targets of the predator's attack. Our data suggest a model of the myxobacterial predatome (genes and proteins associated with predation) in which the predator constitutively produces secretions which disable its prey whilst simultaneously generating a signal that prey is present. That signal then triggers a regulated feeding response in the predator.

Comparative Genomics of Myxobacterial Chemosensory Systems

Chemosensory systems (CSS) are among the most complex organizations of proteins functioning cooperatively to regulate bacterial motility and other cellular activities. These systems have been studied extensively in bacteria, and usually, they are present as a single system. Eight CSS, the highest number in bacteria, have been reported in Myxococcus xanthus DK1622 and are involved in coordinating diverse functions. Here, we have explored and compared the CSS in all available genomes of order Myxococcales. Myxococcales members contain 97 to 476 two-component system (TCS) proteins, which assist the bacteria in surviving and adapting to varying environmental conditions. The number of myxobacterial CSS ranges between 1 and 12, with the largest number in family Cystobacteraceae and the smallest in Nannocystaceae CheA protein was used as a phylogenetic marker to infer evolutionary relatedness between different CSS, and six novel CSS ("extra CSS" [ECSS]) were thus identified in the myxobacteria besides the previously reported Che1 to Che8 systems from M. xanthus Che1 to Che8 systems are monophyletic to deltaproteobacteria, whereas the newly identified ECSS form separate clades with different bacterial classes. The comparative modular organization was concordant with the phylogeny. Four clusters lacking CheA proteins were also identified via CheB-based phylogenetic analysis and were categorized as accessory CSS (ACSS). In Archangium, an orphan CSS was identified, in which both CheA and CheB were absent. The novel, accessory, and orphan multimodular CSS identified here suggest the emergence of myxobacterial CSS and could assist in further characterizing their roles.IMPORTANCE This study is focused on chemosensory systems (CSS), which help the bacterium in directing its movement toward or away from chemical gradients. CSS are present as a single system in most of the bacteria except in some groups, including Myxococcus xanthus, which has 8 CSS, the highest number reported to date. This is the first comprehensive study carrying out a comparative analysis of the 22 available myxobacterial genomes, which suggests the evolutionary diversity of these systems. We are interested in understanding the distribution of CSS within all known myxobacteria and their probable evolution.

Myxobacteria Are Able to Prey Broadly upon Clinically-Relevant Pathogens, Exhibiting a Prey Range Which Cannot Be Explained by Phylogeny

Myxobacteria are natural predators of microorganisms and the subjects of concerted efforts to identify novel antimicrobial compounds. Myxobacterial predatory activity seems to require more than just the possession of specific antimicrobial metabolites. Thus a holistic approach to studying predation promises novel insights into antimicrobial action. Here, we report the isolation of 113 myxobacteria from samples of soil taken from a range of habitats in mid Wales. Predatory activity of each isolate was quantified against a panel of clinically important prey organisms, including Klebsiella pneumoniae, Proteus mirabilis, Candida albicans, Enterococcus faecalis, and three species of Staphylococcus. Myxobacterial isolates exhibited a wide range of predation activity profiles against the panel of prey. Efficient predation of all prey by isolates within the collection was observed, with K. pneumoniae and C. albicans proving particularly susceptible to myxobacterial predation. Notably efficient predators tended to be proficient at predating multiple prey organisms, suggesting they possess gene(s) encoding a broad range killing activity. However, predatory activity was not congruent with phylogeny, suggesting prey range is subject to relatively rapid specialization, potentially involving lateral gene transfer. The broad but patchy prey ranges observed for natural myxobacterial isolates also implies multiple (potentially overlapping) genetic determinants are responsible for dictating predatory activity.

MetaPred2CS: a sequence-based meta-predictor for protein-protein interactions of prokaryotic two-component system proteins

MOTIVATION

Two-component systems (TCS) are the main signalling pathways of prokaryotes, and control a wide range of biological phenomena. Their functioning depends on interactions between TCS proteins, the specificity of which is poorly understood.

RESULTS

The MetaPred2CS web-server interfaces a sequence-based meta-predictor specifically designed to predict pairing of the histidine kinase and response-regulator proteins forming TCSs. MetaPred2CS integrates six sequence-based methods using a support vector machine classifier and has been intensively tested under different benchmarking conditions: (i) species specific gene sets; (ii) neighbouring versus orphan pairs; and (iii) k-fold cross validation on experimentally validated datasets.

AVAILABILITY AND IMPLEMENTATION

Web server at: http://metapred2cs.ibers.aber.ac.uk/, Source code: https://github.com/martinjvickers/MetaPred2CS or implemented as Virtual Machine at: http://metapred2cs.ibers.aber.ac.uk/download CONTACT: naf4@aber.ac.ukSupplementary information: Supplementary data are available at Bioinformatics online.

Myxobacteria: Moving, Killing, Feeding, and Surviving Together

Myxococcus xanthus, like other myxobacteria, is a social bacterium that moves and feeds cooperatively in predatory groups. On surfaces, rod-shaped vegetative cells move in search of the prey in a coordinated manner, forming dynamic multicellular groups referred to as swarms. Within the swarms, cells interact with one another and use two separate locomotion systems. Adventurous motility, which drives the movement of individual cells, is associated with the secretion of slime that forms trails at the leading edge of the swarms. It has been proposed that cellular traffic along these trails contributes to M. xanthus social behavior via stigmergic regulation. However, most of the cells travel in groups by using social motility, which is cell contact-dependent and requires a large number of individuals. Exopolysaccharides and the retraction of type IV pili at alternate poles of the cells are the engines associated with social motility. When the swarms encounter prey, the population of M. xanthus lyses and takes up nutrients from nearby cells. This cooperative and highly density-dependent feeding behavior has the advantage that the pool of hydrolytic enzymes and other secondary metabolites secreted by the entire group is shared by the community to optimize the use of the degradation products. This multicellular behavior is especially observed in the absence of nutrients. In this condition, M. xanthus swarms have the ability to organize the gliding movements of 1000s of rods, synchronizing rippling waves of oscillating cells, to form macroscopic fruiting bodies, with three subpopulations of cells showing division of labor. A small fraction of cells either develop into resistant myxospores or remain as peripheral rods, while the majority of cells die, probably to provide nutrients to allow aggregation and spore differentiation. Sporulation within multicellular fruiting bodies has the benefit of enabling survival in hostile environments, and increases germination and growth rates when cells encounter favorable conditions. Herein, we review how these social bacteria cooperate and review the main cell–cell signaling systems used for communication to maintain multicellularity.