Genome evolution and the emergence of fruiting body development in Myxococcus xanthus

The diversity of physiology and metabolism in chlorophototrophic bacteria

Photosynthesis by (bacterio)chlorophyll-producing organisms ("chlorophototrophy") sustains virtually all life on Earth, providing the biosphere with food and energy. The oxygenic process carried out by plants, algae and cyanobacteria also generates the oxygen we breathe, and ancient cyanobacteria were responsible for oxygenating the atmosphere, creating the conditions that allowed the evolution of complex life. Cyanobacteria were also the endosymbiotic progenitors of chloroplasts, play major roles in biogeochemical cycles and as primary producers in aquatic ecosystems, and act as genetically tractable model organisms for studying oxygenic photosynthesis. In addition to the Cyanobacteriota, eight other bacterial phyla, namely Proteobacteria/Pseudomonadota, Chlorobiota, Chloroflexota, Bacillota, Acidobacteriota, Gemmatimonadota, Vulcanimicrobiota and Myxococcota contain at least one putative chlorophototrophic species, all of which perform a variant of anoxygenic photosynthesis, which does not yield oxygen as a by-product. These chlorophototrophic organisms display incredible diversity in the habitats that they colonise, and in their biochemistry, physiology and metabolism, with variation in the light-harvesting complexes and pigments they produce to utilise solar energy. Whilst some are very well understood, such as the proteobacterial 'purple bacteria', others have only been identified in the last few years and therefore relatively little is known about them - especially those that have not yet been isolated and cultured. In this chapter, we aim to summarise and compare the photosynthetic physiology and central metabolic processes of chlorophototrophic members from the nine phyla in which they are found, giving both a short historical perspective and highlighting gaps in our understanding.

Deletion of an sRNA primes development in a multicellular bacterium

Small non-coding RNAs (sRNAs) regulate gene expression of many biological processes. During growth, some myxobacteria produce an sRNA-Pxr-that blocks fruiting-body development, an aggregative multicellular process typically triggered by starvation. Deleting the pxr gene allows Myxococcus xanthus to develop despite nutrient availability, but Pxr binding targets and the genes regulated by Pxr remain unknown. Here, after showing that Pxr controls the temporal dynamics of development, we compare the transcriptomes of vegetative M. xanthus cells possessing vs. lacking pxr. Over half of the genes impacted by pxr deletion are linked to development, including known and previously undiscovered critical regulators. Pxr also positively regulates genes associated with general metabolic processes. Our study discovers phenotypic effects of Pxr regulation with ecological importance, identifies the suite of genes this sRNA controls during vegetative growth and reveals a previously unknown developmental regulator. These findings provide insights into the molecular mechanism controlling myxobacterial development.

Deciphering the biosynthetic landscape of biofilms in glacier-fed streams

Glacier-fed streams are permanently cold, ultra-oligotrophic, and physically unstable environments, yet microbial life thrives in benthic biofilm communities. Within biofilms, microorganisms rely on secondary metabolites for communication and competition. However, the diversity and genetic potential of secondary metabolites in glacier-fed stream biofilms remain poorly understood. In this study, we present the first large-scale exploration of biosynthetic gene clusters (BGCs) from benthic glacier-fed stream biofilms sampled by the Vanishing Glaciers project from the world's major mountain ranges. We found a remarkable diversity of BGCs, with more than 8,000 of them identified within 2,868 prokaryotic metagenome-assembled genomes, some of them potentially conferring ecological advantages, such as UV protection and quorum sensing. The BGCs were distinct from those sourced from other aquatic microbiomes, with over 40% of them being novel. The glacier-fed stream BGCs exhibited the highest similarity to BGCs from glacier microbiomes. BGC composition displayed geographic patterns and correlated with prokaryotic alpha diversity. We also found that BGC diversity was positively associated with benthic chlorophyll a and prokaryotic diversity, indicative of more biotic interactions in more extensive biofilms. Our study provides new insights into a hitherto poorly explored microbial ecosystem, which is now changing at a rapid pace as glaciers are shrinking due to climate change.

IMPORTANCE

Glacier-fed streams are characterized by low temperatures, high turbidity, and high flow. They host a unique microbiome within biofilms, which form the foundation of the food web and contribute significantly to biogeochemical cycles. Our investigation into secondary metabolites, which likely play an important role in these complex ecosystems, found a unique genetic potential distinct from other aquatic environments. We found the potential to synthesize several secondary metabolites, which may confer ecological advantages, such as UV protection and quorum sensing. This biosynthetic diversity was positively associated with the abundance and complexity of the microbial community, as well as concentrations of chlorophyll a. In the face of climate change, our study offers new insights into a vanishing ecosystem.

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.

Multiple roads lead to Rome: unique morphology and chemistry of endospores, exospores, myxospores, cysts and akinetes in bacteria

The production of specialized resting cells is a remarkable survival strategy developed by many organisms to withstand unfavourable environmental factors such as nutrient depletion or other changes in abiotic and/or biotic conditions. Five bacterial taxa are recognized to form specialized resting cells: Firmicutes, forming endospores; Actinobacteria, forming exospores; Cyanobacteria, forming akinetes; the δ-Proteobacterial order Myxococcales, forming myxospores; and Azotobacteraceae, forming cysts. All these specialized resting cells are characterized by low-to-absent metabolic activity and higher resistance to environmental stress (desiccation, heat, starvation, etc.) when compared to vegetative cells. Given their similarity in function, we tested the potential existence of a universal morpho-chemical marker for identifying these specialized resting cells. After the production of endospores, exospores, akinetes and cysts in model organisms, we performed the first cross-species morphological and chemical comparison of bacterial sporulation. Cryo-electron microscopy of vitreous sections (CEMOVIS) was used to describe near-native morphology of the resting cells in comparison to the morphology of their respective vegetative cells. Resting cells shared a thicker cell envelope as their only common morphological feature. The chemical composition of the different specialized resting cells at the single-cell level was investigated using confocal Raman microspectroscopy. Our results show that the different specialized cells do not share a common chemical signature, but rather each group has a unique signature with a variable conservation of the signature of the vegetative cells. Additionally, we present the validation of Raman signatures associated with calcium dipicolinic acid (CaDPA) and their variation across individual cells to develop specific sorting thresholds for the isolation of endospores. This provides a proof of concept of the feasibility of isolating bacterial spores using a Raman-activated cell-sorting platform. This cross-species comparison and the current knowledge of genetic pathways inducing the formation of the resting cells highlights the complexity of this convergent evolutionary strategy promoting bacterial survival.

Hidden paths to endless forms most wonderful: ecology latently shapes evolution of multicellular development in predatory bacteria

Ecological causes of developmental evolution, for example from predation, remain much investigated, but the potential importance of latent phenotypes in eco-evo-devo has received little attention. Using the predatory bacterium Myxococcus xanthus, which undergoes aggregative fruiting body development upon starvation, we tested whether adaptation to distinct growth environments that do not induce development latently alters developmental phenotypes under starvation conditions that do induce development. In an evolution experiment named MyxoEE-3, growing M. xanthus populations swarmed across agar surfaces while adapting to conditions varying at factors such as surface stiffness or prey identity. Such ecological variation during growth was found to greatly impact the latent evolution of development, including fruiting body morphology, the degree of morphological trait correlation, reaction norms, degrees of developmental plasticity and stochastic diversification. For example, some prey environments promoted retention of developmental proficiency whereas others led to its systematic loss. Our results have implications for understanding evolutionary interactions among predation, development and motility in myxobacterial life cycles, and, more broadly, how ecology can profoundly shape the evolution of developmental systems latently rather than by direct selection on developmental features.

Genome size distributions in bacteria and archaea are strongly linked to evolutionary history at broad phylogenetic scales

The evolutionary forces that determine genome size in bacteria and archaea have been the subject of intense debate over the last few decades. Although the preferential loss of genes observed in prokaryotes is explained through the deletional bias, factors promoting and preventing the fixation of such gene losses often remain unclear. Importantly, statistical analyses on this topic typically do not consider the potential bias introduced by the shared ancestry of many lineages, which is critical when using species as data points because of the potential dependence on residuals. In this study, we investigated the genome size distributions across a broad diversity of bacteria and archaea to evaluate if this trait is phylogenetically conserved at broad phylogenetic scales. After model fit, Pagel's lambda indicated a strong phylogenetic signal in genome size data, suggesting that the diversification of this trait is influenced by shared evolutionary histories. We used a phylogenetic generalized least-squares analysis (PGLS) to test whether phylogeny influences the predictability of genome size from dN/dS ratios and 16S copy number, two variables that have been previously linked to genome size. These results confirm that failure to account for evolutionary history can lead to biased interpretations of genome size predictors. Overall, our results indicate that although bacteria and archaea can rapidly gain and lose genetic material through gene transfers and deletions, respectively, phylogenetic signal for genome size distributions can still be recovered at broad phylogenetic scales that should be taken into account when inferring the drivers of genome size evolution.

A noncanonical cytochrome c stimulates calcium binding by PilY1 for type IVa pili formation

Type IVa pili (T4aP) are versatile bacterial cell surface structures that undergo extension/adhesion/retraction cycles powered by the cell envelope-spanning T4aP machine. In this machine, a complex composed of four minor pilins and PilY1 primes T4aP extension and is also present at the pilus tip mediating adhesion. Similar to many several other bacteria, Myxococcus xanthus contains multiple minor pilins/PilY1 sets that are incompletely understood. Here, we report that minor pilins and PilY1 (PilY1.1) of cluster_1 form priming and tip complexes contingent on calcium and a noncanonical cytochrome c (TfcP) with an unusual His/Cys heme ligation. We provide evidence that TfcP is unlikely to participate in electron transport and instead stimulates calcium binding by PilY1.1 at low-calcium concentrations, thereby stabilizing PilY1.1 and enabling T4aP function in a broader range of calcium concentrations. These results not only identify a previously undescribed function of cytochromes c but also illustrate how incorporation of an accessory factor expands the environmental range under which the T4aP system functions.

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

Complete genome sequence and identification of polyunsaturated fatty acid biosynthesis genes of the myxobacterium Minicystis rosea DSM 24000T

Background

Myxobacteria harbor numerous biosynthetic gene clusters that can produce a diverse range of secondary metabolites. Minicystis rosea DSM 24000^T is a soil-dwelling myxobacterium belonging to the suborderSorangiineae and family Polyangiaceae and is known to produce various secondary metabolites as well as polyunsaturated fatty acids (PUFAs). Here, we use whole-genome sequencing to explore the diversity of biosynthetic gene clusters in M. rosea.

Results

Using PacBio sequencing technology, we assembled the 16.04 Mbp complete genome of M. rosea DSM 24000^T, the largest bacterial genome sequenced to date. About 44% of its coding potential represents paralogous genes predominantly associated with signal transduction, transcriptional regulation, and protein folding. These genes are involved in various essential functions such as cellular organization, diverse niche adaptation, and bacterial cooperation, and enable social behavior like gliding motility, sporulation, and predation, typical of myxobacteria. A profusion of eukaryotic-like kinases (353) and an elevated ratio of phosphatases (8.2/1) in M. rosea as compared to other myxobacteria suggest gene duplication as one of the primary modes of genome expansion. About 7.7% of the genes are involved in the biosynthesis of a diverse array of secondary metabolites such as polyketides, terpenes, and bacteriocins. Phylogeny of the genes involved in PUFA biosynthesis (pfa) together with the conserved synteny of the complete pfa gene cluster suggests acquisition via horizontal gene transfer from Actinobacteria.

Conclusion

Overall, this study describes the complete genome sequence of M. rosea, comparative genomic analysis to explore the putative reasons for its large genome size, and explores the secondary metabolite potential, including the biosynthesis of polyunsaturated fatty acids.

Supplementary Information

The online version contains supplementary material available at 10.1186/s12864-021-07955-x.

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.

Insights into the persistence and phenotypic effects of the endogenous and cryptic plasmid pMF1 in its host strain Myxococcus fulvus 124B02

Many endogenous plasmids carry no noticeable benefits for their bacterial hosts, and the persistence of these 'cryptic plasmids' and their functional impacts are mostly unclear. In this study, we investigated these uncertainties using the social bacterium Myxococcus fulvus 124B02 and its endogenous plasmid pMF1. pMF1 possesses diverse genes that originated from myxobacteria, suggesting a longstanding co-existence of the plasmid with various myxobacterial species. The curing of pMF1 from 124B02 had almost no phenotypic effects on the host. Laboratory evolution experiments showed that the 124B02 strain retained pMF1 when subcultured on dead Escherichia coli cells but lost pMF1 when subcultured on living E. coli cells or on casitone medium; these results indicated that the persistence of pMF1 in 124B02 was environment-dependent. Curing pMF1 caused the mutant to lose the ability to predate and develop fruiting bodies more quickly than the pMF1-containing strain after they were subcultured on dead E. coli cells, which indicated that the presence of pMF1 in M. fulvus 124B02 has some long-term effects on its host. The results provide some new insights into the persistence and impacts of cryptic plasmids in their natural bacterial cells.

Evolutionary dynamics of natural product biosynthesis in bacteria

Covering: 2008 up to 2019The forces of biochemical adaptive evolution operate at the level of genes, manifesting in complex phenotypes and the global biodiversity of proteins and metabolites. While evolutionary histories have been deciphered for some other complex traits, the origins of natural product biosynthesis largely remain a mystery. This fundamental knowledge gap is surprising given the many decades of research probing the genetic, chemical, and biophysical mechanisms of bacterial natural product biosynthesis. Recently, evolutionary thinking has begun to permeate this otherwise mechanistically dominated field. Natural products are now sometimes referred to as 'specialized' rather than 'secondary' metabolites, reinforcing the importance of their biological and ecological functions. Here, we review known evolutionary mechanisms underlying the overwhelming chemical diversity of bacterial secondary metabolism, focusing on enzyme promiscuity and the evolution of enzymatic domains that enable metabolic traits. We discuss the mechanisms that drive the assembly of natural product biosynthetic gene clusters and propose formal definitions for 'specialized' and 'secondary' metabolism. We further explore how biosynthetic gene clusters evolve to synthesize related molecular species, and in turn how the biological and ecological roles that emerge from metabolic diversity are acted on by selection. Finally, we reconcile chemical, functional, and genetic data into an evolutionary model, the dynamic chemical matrix evolutionary hypothesis, in which the relationships between chemical distance, biomolecular activity, and relative fitness shape adaptive landscapes.

Peripheral rods: a specialized developmental cell type in Myxococcus xanthus

In response to nutrient deprivation, the ubiquitous Gram-negative soil bacterium Myxococcus xanthus undergoes a well-characterized developmental response, resulting in the formation of a multicellular fruiting body. The center of the fruiting body consists of myxospores; surrounding this structure are rod-shaped peripheral cells. Unlike spores, the peripheral rods are a metabolically active cell type that inhabits nutrient-deprived environments. The survival characteristics exhibited by peripheral rods, protection from oxidative stress and heat shock, are common survival characteristics exhibited by cells in stationary phase including modifications to morphology and metabolism. Vegetative M. xanthus cells undergo a number of physiological changes during the transition into stationary phase similar to other proteobacteria. In M. xanthus, stationary-phase cells are not considered a component of the developmental response and occur when cells are grown on nutrient-rich plates or in dispersed aqueous media. However, this cell type is not routinely studied and little of its physiology is known. Similarities between these two stress-induced cell types led to the question of whether peripheral rods are actually a distinct developmental cell type or simply cells in stationary phase. In this study, we examine the transcriptome of peripheral rods and its relationship to development. This work demonstrates that peripheral rods are in fact a distinct developmentally differentiated cell type. Although peripheral rods and stationary phase cells display similar characteristics, each transcriptomic pattern is unique and quite different from that of any other M. xanthus cell type.

Orphan Hybrid Histidine Protein Kinase SinK Acts as a Signal Integrator To Fine-Tune Multicellular Behavior in Myxococcus xanthus

His-Asp phosphorelay (also known as two-component signal transduction) proteins are the predominant mechanism used in most bacteria to control behavior in response to changing environmental conditions. In addition to systems consisting of a simple two-component system utilizing an isolated histidine kinase/response regulator pair, some bacteria are enriched in histidine kinases that serve as signal integration proteins; these kinases are usually characterized by noncanonical domain architecture, and the responses that they regulate may be difficult to identify. The environmental bacterium Myxococcus xanthus is highly enriched in these noncanonical histidine kinases. M. xanthus is renowned for a starvation-induced multicellular developmental program in which some cells are induced to aggregate into fruiting bodies and then differentiate into environmentally resistant spores. Here, we characterize the M. xanthus orphan hybrid histidine kinase SinK (Mxan_4465), which consists of a histidine kinase transmitter followed by two receiver domains (REC₁ and REC₂). Nonphosphorylatable sinK mutants were analyzed under two distinct developmental conditions and using a new high-resolution developmental assay. These assays revealed that SinK autophosphorylation and REC₁ impact the onset of aggregation and/or the mobility of aggregates, while REC₂ impacts sporulation efficiency. SinK activity is controlled by a genus-specific hypothetical protein (SinM; Mxan_4466). We propose that SinK serves to fine-tune fruiting body morphology in response to environmental conditions.IMPORTANCE Biofilms are multicellular communities of microorganisms that play important roles in host disease or environmental biofouling. Design of preventative strategies to block biofilms depends on understanding the molecular mechanisms used by microorganisms to build them. The production of biofilms in bacteria often involves two-component signal transduction systems in which one protein component (a kinase) detects an environmental signal and, through phosphotransfer, activates a second protein component (a response regulator) to change the transcription of genes necessary to produce a biofilm. We show that an atypical kinase, SinK, modulates several distinct stages of specialized biofilm produced by the environmental bacterium Myxococcus xanthus SinK likely integrates multiple signals to fine-tune biofilm formation in response to distinct environmental conditions.

The plethora of membrane respiratory chains in the phyla of life

The diversity of microbial cells is reflected in differences in cell size and shape, motility, mechanisms of cell division, pathogenicity or adaptation to different environmental niches. All these variations are achieved by the distinct metabolic strategies adopted by the organisms. The respiratory chains are integral parts of those strategies especially because they perform the most or, at least, most efficient energy conservation in the cell. Respiratory chains are composed of several membrane proteins, which perform a stepwise oxidation of metabolites toward the reduction of terminal electron acceptors. Many of these membrane proteins use the energy released from the oxidoreduction reaction they catalyze to translocate charges across the membrane and thus contribute to the establishment of the membrane potential, i.e. they conserve energy. In this work we illustrate and discuss the composition of the respiratory chains of different taxonomic clades, based on bioinformatic analyses and on biochemical data available in the literature. We explore the diversity of the respiratory chains of Animals, Plants, Fungi and Protists kingdoms as well as of Prokaryotes, including Bacteria and Archaea. The prokaryotic phyla studied in this work are Gammaproteobacteria, Betaproteobacteria, Epsilonproteobacteria, Deltaproteobacteria, Alphaproteobacteria, Firmicutes, Actinobacteria, Chlamydiae, Verrucomicrobia, Acidobacteria, Planctomycetes, Cyanobacteria, Bacteroidetes, Chloroflexi, Deinococcus-Thermus, Aquificae, Thermotogae, Deferribacteres, Nitrospirae, Euryarchaeota, Crenarchaeota and Thaumarchaeota.

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.

An Evo-Devo Perspective on Multicellular Development of Myxobacteria

The transition to multicellularity, recognized as one the major transitions in evolution, has occurred independently several times. While multicellular development has been extensively studied in zygotic organisms including plant and animal groups, just a few aggregative multicellular organisms have been employed as model organisms for the study of multicellularity. Studying different evolutionary origins and modes of multicellularity enables comparative analyses that can help identifying lineage-specific aspects of multicellular evolution and generic factors and mechanisms involved in the transition to multicellularity. Among aggregative multicellular organisms, myxobacteria are a valuable system to explore the particularities that aggregation confers to the evolution of multicellularity and mechanisms shared with clonal organisms. Moreover, myxobacteria species develop fruiting bodies displaying a range of morphological diversity. In this review, we aim to synthesize diverse lines of evidence regarding myxobacteria development and discuss them in the context of Evo-Devo concepts and approaches. First, we briefly describe the developmental processes in myxobacteria, present an updated comparative analysis of the genes involved in their developmental processes and discuss these and other lines of evidence in terms of co-option and developmental system drift, two concepts key to Evo-Devo studies. Next, as has been suggested from Evo-Devo approaches, we discuss how broad comparative studies and integration of diverse genetic, physicochemical, and environmental factors into experimental and theoretical models can further our understanding of myxobacterial development, phenotypic variation, and 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.

Complete Genome Sequence and Comparative Genomics of a Novel Myxobacterium Myxococcus hansupus

Myxobacteria, a group of Gram-negative aerobes, belong to the class δ-proteobacteria and order Myxococcales. Unlike anaerobic δ-proteobacteria, they exhibit several unusual physiogenomic properties like gliding motility, desiccation-resistant myxospores and large genomes with high coding density. Here we report a 9.5 Mbp complete genome of Myxococcus hansupus that encodes 7,753 proteins. Phylogenomic and genome-genome distance based analysis suggest that Myxococcus hansupus is a novel member of the genus Myxococcus. Comparative genome analysis with other members of the genus Myxococcus was performed to explore their genome diversity. The variation in number of unique proteins observed across different species is suggestive of diversity at the genus level while the overrepresentation of several Pfam families indicates the extent and mode of genome expansion as compared to non-Myxococcales δ-proteobacteria.

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

BACKGROUND

Two-component systems (TCSs) are abundant prokaryotic signaling pathways, whose evolution is of particular importance because of their role in bacterial pathogenicity. Comparative genomics can provide important insights into the evolution of these genes, but inferences are dependent on the relatedness of the compared genomes. This study investigated the relationship between evolutionary distance and TCS evolution in myxobacterial genomes, of which there are several sequenced examples, of varying relatedness, and which encode large numbers of TCSs.

METHODS

Myxobacterial TCS gene sets were compared, orthologues defined, and changes in TCS properties such as gene organisation, domain architecture and size identified.

RESULTS

Genome relatedness/evolutionary distance was found to have a large effect on the apparent frequency of evolutionary events affecting TCS genes, but not on the relative dominance of different types of mutations. Large (≥1 gene) indels were the most common changes, often giving rise to gene organisation changes. Smaller indels were also common, sometimes changing domain architecture, and/or leading to pseudogene formation. Individuality of myxobacterial TCS gene sets seems primarily due to lineage specific gene loss. However, there is also evidence of extensive acquisition of genes by lateral transfer, with gene duplication also creating new TCS genes.

CONCLUSIONS

This study provides catalogues of myxobacterial TCS gene sets and their orthology relationships, benchmarked against genome relatedness. It also provides insights into the relationship between evolutionary distance and the inference of TCS estudies of TCS evolution beyond the myxobacteriavolution, which may be important for studies of TCS evolutiThe online version of this articleon beyond the myxobacteria.

Allopatric integrations selectively change host transcriptomes, leading to varied expression efficiencies of exotic genes in Myxococcus xanthus

Background

Exotic genes, especially clustered multiple-genes for a complex pathway, are normally integrated into chromosome for heterologous expression. The influences of insertion sites on heterologous expression and allotropic expressions of exotic genes on host remain mostly unclear.

Results

We compared the integration and expression efficiencies of single and multiple exotic genes that were inserted into Myxococcus xanthus genome by transposition and attB-site-directed recombination. While the site-directed integration had a rather stable chloramphenicol acetyl transferase (CAT) activity, the transposition produced varied CAT enzyme activities. We attempted to integrate the 56-kb gene cluster for the biosynthesis of antitumor polyketides epothilones into M. xanthus genome by site-direction but failed, which was determined to be due to the insertion size limitation at the attB site. The transposition technique produced many recombinants with varied production capabilities of epothilones, which, however, were not paralleled to the transcriptional characteristics of the local sites where the genes were integrated. Comparative transcriptomics analysis demonstrated that the allopatric integrations caused selective changes of host transcriptomes, leading to varied expressions of epothilone genes in different mutants.

Conclusions

With the increase of insertion fragment size, transposition is a more practicable integration method for the expression of exotic genes. Allopatric integrations selectively change host transcriptomes, which lead to varied expression efficiencies of exotic genes.

Electronic supplementary material

The online version of this article (doi:10.1186/s12934-015-0294-5) contains supplementary material, which is available to authorized users.

On the evolution of bacterial multicellularity

Multicellularity is one of the most prevalent evolutionary innovations and nowhere is this more apparent than in the bacterial world, which contains many examples of multicellular organisms in a surprising array of forms. Due to their experimental accessibility and the large and diverse genomic data available, bacteria enable us to probe fundamental aspects of the origins of multicellularity. Here we discuss examples of multicellular behaviors in bacteria, the selective pressures that may have led to their evolution, possible origins and intermediate stages, and whether the ubiquity of apparently convergent multicellular forms argues for its inevitability.

Characterization of four type IV pilin homologues in Stigmatella aurantiaca DSM17044 by heterologous expression in Myxococcus xanthus

As prokaryotic models for multicellular development, Stigmatella aurantiaca and Myxococcus xanthus share many similarities in terms of social behaviors, such as gliding motility. Our current understanding of myxobacterial grouped-cell motilities comes mainly from the research on M. xanthus, which shows that filamentous type IV pili (TFP), composed of type IV pilin (also called PilA protein) subunits, are the key apparatus for social motility (S-motility). However, little is known about the pilin protein in S. aurantiaca. We cloned and sequenced four genes (pilA(Sa1~4)) from S. aurantiaca DSM17044 that are homologous to pilA(Mx) (pilA gene in M. xanthus DK1622). The homology and similarities among pilA(Sa) proteins and other myxobacterial homologues were systematically analyzed. To determine their potential biological functions, the four pilA(Sa) genes were expressed in M. xanthus DK10410 (ΔpilA(Mx)), which did not restore S-motility on soft agar or EPS production to host cells. After further analysis of the motile behaviors in a methylcellulose solution, the M. xanthus strains were categorized into three types. YL6101, carrying pilA(Sa1), and YL6104, carrying pilA(Sa4), produced stable but unretractable surface pili; YL6102, carrying pilA(Sa2), produced stable surface pili and exhibited reduced TFP-dependent motility in methylcellulose; YL6103, carrying pilA(Sa3), produced unstable surface pili. Based on these findings, we propose that pilA(Sa2) might be responsible for the type IV pilin production involved in group motility in S. aurantiaca DSM17044. After examining the developmental processes, it was suggested that the expression of PilA(Sa4) protein might have positive effects on the fruiting body formation of M. xanthus DK10410 cells. Moreover, the formation of fruiting body in M. xanthus cells with stable exogenous TFPSa were compensated by mixing them with S. aurantiaca DSM17044 cells. Our results shed some light on the features and functions of type IV pilin homologues in S. aurantiaca.

The clpB gene is involved in the stress response of Myxococcus xanthus during vegetative growth and development

The Clp/HSP100 family of molecular chaperones is ubiquitous in both prokaryotes and eukaryotes. These proteins play important roles in refolding, disaggregating and degrading proteins damaged by stress. As a subclass of the Clp/HSP100 family, ClpB has been shown to be involved in various stress responses as well as other functions in bacteria. In the present study, we investigated the role of a predicted ClpB-encoding gene, MXAN5092, in the stress response during vegetative growth and development of Myxococcus xanthus. Transcriptional analysis confirmed induction of this clpB homologue under different stress conditions, and further phenotypic analysis revealed that an in-frame deletion mutant of MXAN5092 was more sensitive to various stress treatments than the wild-type strain during vegetative growth. Moreover, the absence of the MXAN5092 gene resulted in decreased heat tolerance of myxospores, indicating the involvement of this clpB homologue in the stress response during the development of myxospores. The M. xanthus recombinant ClpB (MXAN5092) protein also showed a general chaperone activity in vitro. Overall, our genetic and phenotypic analysis of the predicted ATP-dependent chaperone protein ClpB (MXAN5092) demonstrated that it functions as a chaperone protein and plays an important role in cellular stress tolerance during both vegetative growth and development of M. xanthus.

Natural transformation of Myxococcus xanthus

Myxococcus xanthus belongs to the delta class of the proteobacteria and is notable for its complex life-style with social behaviors and relatively large genome. Although previous observations have suggested the existence of horizontal gene transfer in M. xanthus, its ability to take up exogenous DNA via natural transformation has not been experimentally demonstrated. In this study, we achieved natural transformation in M. xanthus using the autonomously replicating myxobacterial plasmid pZJY41 as donor DNA. M. xanthus exopolysaccharide (EPS) was shown to be an extracellular barrier for transformation. Cells deficient in EPS production, e.g., mutant strains carrying ΔdifA or ΔepsA, became naturally transformable. Among the inner barriers to transformation were restriction-modification systems in M. xanthus, which could be partially overcome by methylating DNA in vitro using cell extracts of M. xanthus prior to transformation. In addition, the incubation time of DNA with cells and the presence of divalent magnesium ion affected transformation frequency of M. xanthus. Furthermore, we also observed a potential involvement of the type IV pilus system in the DNA uptake machinery of M. xanthus. The natural transformation was totally eliminated in the ΔpilQ/epsA and Δtgl/epsA mutants, and null mutation of pilB or pilC in an ΔepsA background diminished the transformation rate. Our study, to the best of our knowledge, provides the first example of a naturally transformable species among deltaproteobacteria.

Comparative genomic analysis of fruiting body formation in Myxococcales

Genetic programs underlying multicellular morphogenesis and cellular differentiation are most often associated with eukaryotic organisms, but examples also exist in bacteria such as the formation of multicellular, spore-filled fruiting bodies in the order Myxococcales. Most members of the Myxococcales undergo a multicellular developmental program culminating in the formation of spore-filled fruiting bodies in response to starvation. To gain insight into the evolutionary history of fruiting body formation in Myxococcales, we performed a comparative analysis of the genomes and transcriptomes of five Myxococcales species, four of these undergo fruiting body formation (Myxococcus xanthus, Stigmatella aurantiaca, Sorangium cellulosum, and Haliangium ochraceum) and one does not (Anaeromyxobacter dehalogenans). Our analyses show that a set of 95 known M. xanthus development-specific genes--although suffering from a sampling bias--are overrepresented and occur more frequently than an average M. xanthus gene in S. aurantiaca, whereas they occur at the same frequency as an average M. xanthus gene in S. cellulosum and in H. ochraceum and are underrepresented in A. dehalogenans. Moreover, genes for entire signal transduction pathways important for fruiting body formation in M. xanthus are conserved in S. aurantiaca, whereas only a minority of these genes are conserved in A. dehalogenans, S. cellulosum, and H. ochraceum. Likewise, global gene expression profiling of developmentally regulated genes showed that genes that upregulated during development in M. xanthus are overrepresented in S. aurantiaca and slightly underrepresented in A. dehalogenans, S. cellulosum, and H. ochraceum. These comparative analyses strongly indicate that the genetic programs for fruiting body formation in M. xanthus and S. aurantiaca are highly similar and significantly different from the genetic program directing fruiting body formation in S. cellulosum and H. ochraceum. Thus, our analyses reveal an unexpected level of plasticity in the genetic programs for fruiting body formation in the Myxococcales and strongly suggest that the genetic program underlying fruiting body formation in different Myxococcales is not conserved. The evolutionary implications of this finding are discussed.

The alternative complex III: a different architecture using known building modules

Until recently cytochrome bc(1) complexes were the only enzymes known to be able to transfer electrons from reduced quinones to cytochrome c. However, a complex with the same activity and with a unique subunit composition was purified from the membranes of Rhodothermus marinus. This complex, named alternative complex III (ACIII) was then biochemical, spectroscopic and genetically characterized. Later it was observed that the presence of ACIII was not exclusive of R. marinus being the genes coding for ACIII widespread, at least in the Bacteria domain. In this work, a comprehensive description of the current knowledge on ACIII is presented. The relation of ACIII with members of the complex iron-sulfur molybdoenzyme family is investigated by analyzing all the available completely sequenced genomes. It is concluded that ACIII is a new complex composed by a novel combination of modules already identified in other respiratory complexes.

Computational identification of rare codons of Escherichia coli based on codon pairs preference

Background

Codon bias is believed to play an important role in the control of gene expression. In Escherichia coli, some rare codons, which can limit the expression level of exogenous protein, have been defined by gene engineering operations. Previous studies have confirmed the existence of codon pair's preference in many genomes, but the underlying cause of this bias has not been well established. Here we focus on the patterns of rarely-used synonymous codons. A novel method was introduced to identify the rare codons merely by codon pair bias in Escherichia coli.

Results

In Escherichia coli, we defined the "rare codon pairs" by calculating the frequency of occurrence of all codon pairs in coding sequences. Rare codons which are disliked in genes could make great contributions to forming rare codon pairs. Meanwhile our investigation showed that many of these rare codon pairs contain termination codons and the recognized sites of restriction enzymes. Furthermore, a new index (F_rare) was developed. Through comparison with the classical indices we found a significant negative correlation between F_rare and the indices which depend on reference datasets.

Conclusions

Our approach suggests that we can identify rare codons by studying the context in which a codon lies. Also, the frequency of rare codons (F_rare) could be a useful index of codon bias regardless of the lack of expression abundance information.

Complete genome sequence of Haliangium ochraceum type strain (SMP-2)

Haliangium ochraceum Fudou et al. 2002 is the type species of the genus Haliangium in the myxococcal family 'Haliangiaceae'. Members of the genus Haliangium are the first halophilic myxobacterial taxa described. The cells of the species follow a multicellular lifestyle in highly organized biofilms, called swarms, they decompose bacterial and yeast cells as most myxobacteria do. The fruiting bodies contain particularly small coccoid myxospores. H. ochraceum encodes the first actin homologue identified in a bacterial genome. Here we describe the features of this organism, together with the complete genome sequence, and annotation. This is the first complete genome sequence of a member of the myxococcal suborder Nannocystineae, and the 9,446,314 bp long single replicon genome with its 6,898 protein-coding and 53 RNA genes is part of the Genomic Encyclopedia of Bacteria and Archaea project.

Sociobiology of the myxobacteria

Cooperation is integral to much of biological life but can be threatened by selfish evolutionary strategies. Diverse cooperative traits have evolved among microbes, but particularly sophisticated forms of sociality have arisen in the myxobacteria, including group motility and multicellular fruiting body development. Myxobacterial cooperation has succeeded against socially destructive cheaters and can readily re-evolve from some socially defective genotypes. However, social harmony does not extend far. Spatially structured natural populations of the model species Myxococcus xanthus have fragmented into a large number of socially incompatible genotypes that exclude, exploit, and/or antagonize one another, including genetically similar neighbors. Here, we briefly review basic social evolution concepts as they pertain to microbes, discuss potential benefits of myxobacterial social traits, highlight recent empirical studies of social evolution in M. xanthus, and consider their implications for how myxobacterial cooperation and conflict evolve in the wild.

Eukaryotic-like protein kinases in the prokaryotes and the myxobacterial kinome

Ser/Thr/Tyr kinases, which together comprise a major class of regulatory proteins in eukaryotes, were not believed to play an important role in prokaryotes until recently. However, our analysis of 626 prokaryotic genomes reveals that eukaryotic-like protein kinases (ELKs) are found in nearly two-thirds of the sequenced strains. We have identified 2697 ELKs, most of which are encoded by multicellular strains of the phyla Proteobacteria (Myxococcales), Actinobacteria, Cyanobacteria, and Chloroflexi, and 2 Acidobacteria and 1 Planctomycetes. Astonishingly, 7 myxobacterial strains together encode 892 ELKs, with 4 of the strains exhibiting a genomic ELK density similar to that observed in eukaryotes. Most myxobacterial ELKs show a modular organization in which the kinase domain is located at the N terminus. The C-terminal portion of the ELKs is highly diverse and often contains sequences with similarity to characterized domains, most of them involved in signaling mechanisms or in protein-protein interactions. However, many of these architectures are unique to the myxobacteria, an observation that suggests that this group exploits sophisticated and novel signal transduction systems. Phylogenetic reconstruction using the kinase domains revealed many orthologous sequence pairs and a huge number of gene duplications that probably occurred after speciation. Furthermore, studies of the microsynteny in the ELK-encoding regions reveal only low levels of synteny among Myxococcus xanthus, Plesiocystis pacifica, and Sorangium cellulosum. However, extensive similarities between M. xanthus, Stigmatella aurantiaca, and 3 Anaeromyxobacter strains were observed, indicating that they share regulatory pathways involving various ELKs.