A barrier to homologous recombination between sympatric strains of the cooperative soil bacterium Myxococcus xanthus

Mutation Rate and Effective Population Size of the Model Cooperative Bacterium Myxococcus xanthus

Intrinsic rates of genetic mutation have diverged greatly across taxa and exhibit statistical associations with several other parameters and features. These include effective population size (Ne), genome size, and gametic multicellularity, with the latter being associated with both increased mutation rates and decreased effective population sizes. However, data sufficient to test for possible relationships between microbial multicellularity and mutation rate (µ) are lacking. Here, we report estimates of two key population-genetic parameters, Ne and µ, for Myxococcus xanthus, a bacterial model organism for the study of aggregative multicellular development, predation, and social swarming. To estimate µ, we conducted an ∼400-day mutation accumulation experiment with 46 lineages subjected to regular single colony bottlenecks prior to clonal regrowth. Upon conclusion, we sequenced one clonal-isolate genome per lineage. Given collective evolution for 85,323 generations across all lines, we calculate a per base-pair mutation rate of ∼5.5 × 10-10 per site per generation, one of the highest mutation rates among free-living eubacteria. Given our estimate of µ, we derived Ne at ∼107 from neutral diversity at four-fold degenerate sites across two dozen M. xanthus natural isolates. This estimate is below average for eubacteria and strengthens an already clear negative correlation between µ and Ne in prokaryotes. The higher and lower than average mutation rate and Ne for M. xanthus, respectively, amplify the question of whether any features of its multicellular life cycle-such as group-size reduction during fruiting-body development-or its highly structured spatial distribution have significantly influenced how these parameters have evolved.

Dynamics of bacterial recombination in the human gut microbiome

Horizontal gene transfer (HGT) is a ubiquitous force in microbial evolution. Previous work has shown that the human gut is a hotspot for gene transfer between species, but the more subtle exchange of variation within species-also known as recombination-remains poorly characterized in this ecosystem. Here, we show that the genetic structure of the human gut microbiome provides an opportunity to measure recent recombination events from sequenced fecal samples, enabling quantitative comparisons across diverse commensal species that inhabit a common environment. By analyzing recent recombination events in the core genomes of 29 human gut bacteria, we observed widespread heterogeneities in the rates and lengths of transferred fragments, which are difficult to explain by existing models of ecological isolation or homology-dependent recombination rates. We also show that natural selection helps facilitate the spread of genetic variants across strain backgrounds, both within individual hosts and across the broader population. These results shed light on the dynamics of in situ recombination, which can strongly constrain the adaptability of gut microbial communities.

Ongoing diversification of the global fish pathogen Piscirickettsia salmonis through genetic isolation and transposition bursts

The management of bacterial pathogens remains a key challenge of aquaculture. The marine gammaproteobacterium Piscirickettsia salmonis is the etiological agent of piscirickettsiosis and causes multi-systemic infections in different salmon species, resulting in considerable mortality and substantial commercial losses. Here, we elucidate its global diversity, evolution, and selection during human interventions. Our comprehensive analysis of 73 closed, high quality genome sequences covered strains from major outbreaks and was supplemented by an analysis of all P. salmonis 16S rRNA gene sequences and metagenomic reads available in public databases. Genome comparison showed that Piscirickettsia comprises at least three distinct, genetically isolated species of which two showed evidence for continuing speciation. However, at least twice the number of species exist in marine fish or seawater. A hallmark of Piscirickettsia diversification is the unprecedented amount and diversity of transposases which are particularly active in subgroups undergoing rapid speciation and are key to the acquisition of novel genes and to pseudogenization. Several group-specific genes are involved in surface antigen synthesis and may explain the differences in virulence between strains. However, the frequent failure of antibiotic treatment of piscirickettsiosis outbreaks cannot be explained by horizontal acquisition of resistance genes which so far occurred only very rarely. Besides revealing a dynamic diversification of an important pathogen, our study also provides the data for improving its surveillance, predicting the emergence of novel lineages, and adapting aquaculture management, and thereby contributes towards the sustainability of salmon farming.

Two reasons to kill: predation and kin discrimination in myxobacteria

Myxobacteria are social microbial predators that use cell-cell contacts to identify bacterial or fungal prey and to differentiate kin relatives to initiate cellular responses. For prey killing, they assemble Tad-like and type III-like secretion systems at contact sites. For kin discrimination (KD), they assemble outer membrane exchange complexes composed of the TraA and TraB receptors at contacts sites. A type VI secretion system and Rhs proteins also mediate KD. Following cellular recognition, these systems deliver appropriate effectors into target cells. For prey, this leads to cell death and lysis for nutrient consumption by myxobacteria. In KD, a panel of effectors are delivered, and if adjacent cells are clonal cells, resistance ensues because they express a cognate panel of immunity factors; while nonkin lack complete immunity and are intoxicated. This review compares and contrasts recent findings from these systems in myxobacteria.

Social Diversification Driven by Mobile Genetic Elements

Social diversification in microbes is an evolutionary process where lineages bifurcate into distinct populations that cooperate with themselves but not with other groups. In bacteria, this is frequently driven by horizontal transfer of mobile genetic elements (MGEs). Here, the resulting acquisition of new genes changes the recipient's social traits and consequently how they interact with kin. These changes include discriminating behaviors mediated by newly acquired effectors. Since the producing cell is protected by cognate immunity factors, these selfish elements benefit from selective discrimination against recent ancestors, thus facilitating their proliferation and benefiting the host. Whether social diversification benefits the population at large is less obvious. The widespread use of next-generation sequencing has recently provided new insights into population dynamics in natural habitats and the roles MGEs play. MGEs belong to accessory genomes, which often constitute the majority of the pangenome of a taxon, and contain most of the kin-discriminating loci that fuel rapid social diversification. We further discuss mechanisms of diversification and its consequences to populations and conclude with a case study involving myxobacteria.

Allopatric divergence of cooperators confers cheating resistance and limits effects of a defector mutation

BACKGROUND

Social defectors may meet diverse cooperators. Genotype-by-genotype interactions may constrain the ranges of cooperators upon which particular defectors can cheat, limiting cheater spread. Upon starvation, the soil bacterium Myxococcus xanthus cooperatively develops into spore-bearing fruiting bodies, using a complex regulatory network and several intercellular signals. Some strains (cheaters) are unable to sporulate effectively in pure culture due to mutations that reduce signal production but can exploit and outcompete cooperators within mixed groups.

RESULTS

In this study, interactions between a cheater disrupted at the signaling gene csgA and allopatrically diversified cooperators reveal a very small cheating range. Expectedly, the cheater failed to cheat on all natural-isolate cooperators owing to non-cheater-specific antagonisms. Surprisingly, some lab-evolved cooperators had already exited the csgA mutant's cheating range after accumulating fewer than 20 mutations and without experiencing cheating during evolution. Cooperators might also diversify in the potential for a mutation to reduce expression of a cooperative trait or generate a cheating phenotype. A new csgA mutation constructed in several highly diverged cooperators generated diverse sporulation phenotypes, ranging from a complete defect to no defect, indicating that genetic backgrounds can limit the set of genomes in which a mutation creates a defector.

CONCLUSIONS

Our results demonstrate that natural populations may feature geographic mosaics of cooperators that have diversified in their susceptibility to particular cheaters, limiting defectors' cheating ranges and preventing them from spreading. This diversification may also lead to variation in the phenotypes generated by any given cooperation-gene mutation, further decreasing the chance of a cheater emerging which threatens the persistence of cooperation in the system.

Long-term excessive application of K2SO4 fertilizer alters bacterial community and functional pathway of tobacco-planting soil

To improve tobacco leaf quality, excessive K₂SO₄ fertilizers were applied to soils in major tobacco-planting areas in China. However, the effects of K₂SO₄ application on soil microbial community and functions are still unclear. An eight-year field experiment with three kinds of K₂SO₄ amounts (low amount, K₂O 82.57 kg hm^-2, LK; moderate amount, K₂O 165.07 kg hm^-2, MK; high amount, K₂O 247.58 kg hm^-2, HK) was established to assess the effects of K₂SO₄ application on the chemical and bacterial characteristics of tobacco-planting soil using 16S rRNA gene and metagenomic sequencing approaches. Results showed that HK led to lower pH and higher nitrogen (N), potassium (K), sulfur(S) and organic matter contents of the soil than LK. The bacterial community composition of HK was significantly different from those of MK and LK, while these of MK and LK were similar. Compared to LK, HK increased the relative abundance of predicted copiotrophic groups (e.g. Burkholderiaceae, Rhodospirillaceae families and Ellin6067 genus) and potentially beneficial bacteria (e.g. Gemmatimonadetes phylum and Bacillus genus) associated with pathogens and heavy metal resistance, N fixation, dissolution of phosphorus and K. While some oligotrophic taxa (e.g. Acidobacteria phylum) related to carbon, N metabolism exhibited adverse responses to HK. Metagenomic analysis suggested that the improvement of pathways related to carbohydrate metabolism and genetic information processing by HK might be the self-protection mechanism of microorganisms against environmental stress. Besides, the redundancy analysis and variation partitioning analysis showed that soil pH, available K and S were the primary soil factors in shifting the bacterial community and KEGG pathways. This study provides a clear understanding of the responses of soil microbial communities and potential functions to excessive application of K₂SO₄ in tobacco-planting soil.

High genomic differentiation and limited gene flow indicate recent cryptic speciation within the genus Laspinema (cyanobacteria)

The sympatric occurrence of closely related lineages displaying conserved morphological and ecological traits is often characteristic of free-living microbes. Gene flow, recombination, selection, and mutations govern the genetic variability between these cryptic lineages and drive their differentiation. However, sequencing conservative molecular markers (e.g., 16S rRNA) coupled with insufficient population-level sampling hindered the study of intra-species genetic diversity and speciation in cyanobacteria. We used phylogenomics and a population genomic approach to investigate the extent of local genomic diversity and the mechanisms underlying sympatric speciation of Laspinema thermale. We found two cryptic lineages of Laspinema. The lineages were highly genetically diverse, with recombination occurring more frequently within than between them. That suggests the existence of a barrier to gene flow, which further maintains divergence. Genomic regions of high population differentiation harbored genes associated with possible adaptations to high/low light conditions and stress stimuli, although with a weak diversifying selection. Overall, the diversification of Laspinema species might have been affected by both genomic and ecological processes.

Coastal Transient Niches Shape the Microdiversity Pattern of a Bacterioplankton Population with Reduced Genomes

Globally dominant marine bacterioplankton lineages are often limited in metabolic versatility, owing to their extensive genome reductions, and thus cannot take advantage of transient nutrient patches. It is therefore perplexing how the nutrient-poor bulk seawater sustains the pelagic streamlined lineages, each containing numerous populations. Here, we sequenced the genomes of 33 isolates of the recently discovered CHUG lineage (~2.6 Mbp), which have some of the smallest genomes in the globally abundant Roseobacter group (commonly over 4 Mbp). These genome-reduced bacteria were isolated from a transient habitat: seawater surrounding the brown alga, Sargassum hemiphyllum. Population genomic analyses showed that: (i) these isolates, despite sharing identical 16S rRNA genes, were differentiated into several genetically isolated populations through successive speciation events; (ii) only the first speciation event led to the genetic separation of both core and accessory genomes; and (iii) populations resulting from this event are differentiated at many loci involved in carbon utilization and oxygen respiration, corroborated by BiOLOG phenotype microarray assays and oxygen uptake kinetics experiments, respectively. These differentiated traits match well with the dynamic nature of the macroalgal seawater, in which the quantity and quality of carbon sources and the concentration of oxygen likely vary spatially and temporally, though other habitats, like fresh organic aggregates, cannot be ruled out. Our study implies that transient habitats in the overall nutrient-poor ocean can shape the microdiversity and population structure of genome-reduced bacterioplankton lineages.

Signatures of selection in core and accessory genomes indicate different ecological drivers of diversification among Bacillus cereus clades

Bacterial clades are often ecologically distinct, despite extensive horizontal gene transfer (HGT). How selection works on different parts of bacterial pan‐genomes to drive and maintain the emergence of clades is unclear. Focusing on the three largest clades in the diverse and well‐studied Bacillus cereus sensu lato group, we identified clade‐specific core genes (present in all clade members) and then used clade‐specific allelic diversity to identify genes under purifying and diversifying selection. Clade‐specific accessory genes (present in a subset of strains within a clade) were characterized as being under selection using presence/absence in specific clades. Gene ontology analyses of genes under selection revealed that different gene functions were enriched in different clades. Furthermore, some gene functions were enriched only amongst clade‐specific core or accessory genomes. Genes under purifying selection were often clade‐specific, while genes under diversifying selection showed signs of frequent HGT. These patterns are consistent with different selection pressures acting on both the core and the accessory genomes of different clades and can lead to ecological divergence in both cases. Examining variation in allelic diversity allows us to uncover genes under clade‐specific selection, allowing ready identification of strains and their ecological niche.

Population Genomics of Microbial Biostalactites: Non-recombinogenic Genome Islands and Microdiversification by Transposons

Intrapopulation genetic variability in prokaryotes is receiving increasing attention thanks to improving sequencing methods; however, the ability to distinguish intrapopulation variability from species clusters or initial stages of gene flow barrier development remains insufficient. To overcome this limitation, we took advantage of the lifestyle of Ferrovum myxofaciens, a species that may represent 99% of prokaryotic microbiome of biostalactites growing at acid mine drainage springs. We gained four complete and one draft metagenome-assembled F. myxofaciens genomes using Oxford Nanopore and Illumina sequencing and mapped the reads from each sample on the reference genomes to assess the intrapopulation variability. We observed two phenomena associated with intrapopulation variability: hypervariable regions affected by mobilome expansion called “scrapyards,” and variability in gene disruptions caused by transposons within each population. Both phenomena were previously described in prokaryotes. However, we present here for the first time scrapyard regression and the development of a new one. Nearly complete loss of intrapopulation short sequence variability in the old scrapyard and high variability in the new one suggest that localized gene flow suppression is necessary for scrapyard formation. Concerning the variable gene disruptions, up to 9 out of 41 occurrences per sample were located in highly conserved diguanylate cyclases/phosphodiesterases. We propose that microdiversification of life strategies may be an adaptive outcome of random diguanylate cyclase elimination. The mine biostalactites thus proved as a unique model system for describing genomic intrapopulation processes, as they offer easily sampleable units enriched in a single microbial species.

Population differentiation of Rhodobacteraceae along with coral compartments

Coral mucus, tissue, and skeleton harbor compositionally different microbiota, but how these coral compartments shape the microbial evolution remains unexplored. Here, we sampled bacteria inhabiting a prevalent coral species Platygyra acuta and sequenced genomes of 234 isolates comprising two populations in Rhodobacteraceae, an alphaproteobacterial lineage representing a significant but variable proportion (5-50%) of the coral microbiota. The Ruegeria population (20 genomes) contains three clades represented by eight, six, and six isolates predominantly sampled from the skeleton (outgroup), mucus (clade-M), and skeleton (clade-S), respectively. The clade-M possesses functions involved in the utilization of coral osmolytes abundant in the mucus (e.g., methylamines, DMSP, taurine, and L-proline), whereas the clade-S uniquely harbors traits that may promote adaptation to the low-energy and diurnally anoxic skeleton (e.g., sulfur oxidation and swimming motility). These between-clade genetic differences were largely supported by physiological assays. Expanded analyses by including genomes of 24 related isolates (including seven new genomes) from other marine environments suggest that clade-M and clade-S may have diversified in non-coral habitats, but they also consolidated a key role of distinct coral compartments in diversifying many of the above-mentioned traits. The unassigned Rhodobacteraceae population (214 genomes) varies only at a few dozen nucleotide sites across the whole genomes, but the number of between-compartment migration events predicted by the Slatkin-Maddison test supported that dispersal limitation between coral compartments is another key mechanism diversifying microbial populations. Collectively, our results suggest that different coral compartments represent ecologically distinct and microgeographically separate habitats that drive the evolution of the coral microbiota.

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

Kin recognition and outer membrane exchange (OME) in myxobacteria

Myxobacteria conduct complex social traits that requires populations to be highly related and devoid of exploiters. To enrich for clonal cells in populations, they employ kin discrimination mechanisms. One key system involves a polymorphic cell surface receptor, TraA, which recognizes self by homotypic interactions with neighboring myxobacterial cells. Recent studies revealed that TraA and its partner TraB are fluid outer membrane proteins that coalesce into foci upon recognition of kin. The formation of foci leads to transient membrane fusion junctions and the bidirectional exchange of outer membrane components that facilitates cooperative behaviors. Additionally, expansive suites of polymorphic lipoprotein toxins are exchanged, which act as self-identity barcodes that exquisitely discriminate against nonself to assemble homogenous populations.

Cryptic speciation of a pelagic Roseobacter population varying at a few thousand nucleotide sites

A drop of seawater contains numerous microspatial niches at the scale relevant to microbial activities. Examples are abiotic niches such as detrital particles that show different sizes and organic contents, and biotic niches resulting from bacteria-phage and bacteria-phytoplankton interactions. A common practice to investigate the impact of microenvironments on bacterial evolution is to separate the microenvironments physically and compare the bacterial inhabitants from each. It remains poorly understood, however, which microenvironment primarily drives bacterioplankton evolution in the pelagic ocean. By applying a dilution cultivation approach to an undisturbed coastal water sample, we isolate a bacterial population affiliated with the globally dominant Roseobacter group. Although varying at just a few thousand nucleotide sites across the whole genomes, members of this clonal population are diverging into two genetically separated subspecies. Genes underlying speciation are not unique to subspecies but instead clustered at the shared regions that represent ~6% of the genomic DNA. They are primarily involved in vitamin synthesis, motility, oxidative defense, carbohydrate, and amino acid utilization, consistent with the known strategies that roseobacters take to interact with phytoplankton and particles. Physiological assays corroborate that one subspecies outcompetes the other in these traits. Our results indicate that the microenvironments in the pelagic ocean represented by phytoplankton and organic particles are likely important niches that drive the cryptic speciation of the Roseobacter population, though microhabitats contributed by other less abundant pelagic hosts cannot be ruled out.

Soil bacterial populations are shaped by recombination and gene-specific selection across a grassland meadow

Soil microbial diversity is often studied from the perspective of community composition, but less is known about genetic heterogeneity within species. The relative impacts of clonal interference, gene-specific selection, and recombination in many abundant but rarely cultivated soil microbes remain unknown. Here we track genome-wide population genetic variation for 19 highly abundant bacterial species sampled from across a grassland meadow. Genomic inferences about population structure are made using the millions of sequencing reads that are assembled de novo into consensus genomes from metagenomes, as each read pair describes a short genomic sequence from a cell in each population. Genomic nucleotide identity of assembled genomes was significantly associated with local geography for over half of the populations studied, and for a majority of populations within-sample nucleotide diversity could often be as high as meadow-wide nucleotide diversity. Genes involved in metabolite biosynthesis and extracellular transport were characterized by elevated nucleotide diversity in multiple species. Microbial populations displayed varying degrees of homologous recombination and recombinant variants were often detected at 7-36% of loci genome-wide. Within multiple populations we identified genes with unusually high spatial differentiation of alleles, fewer recombinant events, elevated ratios of nonsynonymous to synonymous variants, and lower nucleotide diversity, suggesting recent selective sweeps for gene variants. Taken together, these results indicate that recombination and gene-specific selection commonly shape genetic variation in several understudied soil bacterial lineages.

Rapid diversification of wild social groups driven by toxin-immunity loci on mobile genetic elements

Many species form distinct social groups that provide fitness advantages to individuals. However, the evolutionary processes that generate new social groups are not well understood. Here we examined recently diverged natural isolates of the model social bacterium, Myxococcus xanthus, to probe the genetic mechanisms and evolutionary processes of kin discrimination that occurred naturally in soil. We show that social incompatibilities were formed from horizontal gene transfer of effectors belonging to three distinct polymorphic toxin systems; outer membrane exchange, type VI secretion and rearrangement hotspot systems. Strikingly, the unique toxin effectors and their respective immunity genes that are responsible for social incompatibilities reside on mobile genetic elements, which make up nearly all of the genotypic variation between isolates within clades. By disrupting these three toxin systems, we engineered social harmony between strains that were originally incompatible. In addition, a horizontal allele swap of a single kin recognition receptor changed social interactions and competition outcomes. Our results provide a case study for how horizontal gene transfer led to social diversification in a natural context. Finally, we show how genomic information of kin discriminatory loci can be used to predict social interactions.

Mining the Biosynthetic Potential for Specialized Metabolism of a Streptomyces Soil Community

The diversity and distribution of specialized metabolite gene clusters within a community of bacteria living in the same soil habitat are poorly documented. Here we analyzed the genomes of 8 Streptomyces isolated at micro-scale from a forest soil that belong to the same species or to different species. The results reveal high levels of diversity, with a total of 261 biosynthesis gene clusters (BGCs) encoding metabolites such as terpenes, polyketides (PKs), non-ribosomal peptides (NRPs) and ribosomally synthesized and post-translationally modified peptides (RiPPs) with potential bioactivities. A significant part of these BGCs (n = 53) were unique to only one strain when only 5 were common to all strains. The metabolites belong to very diverse chemical families and revealed that a large diversity of metabolites can potentially be produced in the community. Although that analysis of the global metabolome using GC-MS revealed that most of the metabolites were shared between the strains, they exhibited a specific metabolic pattern. We also observed that the presence of these accessory pathways might result from frequent loss and gain of genes (horizontal transfer), showing that the potential of metabolite production is a dynamic phenomenon in the community. Sampling Streptomyces at the community level constitutes a good frame to discover new biosynthetic pathways and it appears as a promising reservoir for the discovery of new bioactive compounds.

Kin discrimination and outer membrane exchange in Myxococcus xanthus: Experimental analysis of a natural population

In some species of myxobacteria, adjacent cells sufficiently similar at the adhesin protein TraA can exchange components of their outer membranes. The primary benefits of such outer membrane exchange (OME) in natural populations are unclear, but in some OME interactions, transferred OM content can include SitA toxins that kill OME participants lacking an appropriate immunity gene. Such OME-dependent toxin transfer across Myxococcus xanthus strains that differ only in their sitBAI toxin/antitoxin cassette can mediate inter-strain killing and generate colony-merger incompatibilities (CMIs)-inter-colony border phenotypes between distinct genotypes that differ from respective self-self colony interfaces. Here we ask whether OME-dependent toxin transfer is a common cause of prevalent CMIs and antagonisms between M. xanthus natural isolates identical at TraA. We disrupted traA in eleven isolates from a cm-scale soil population and assayed whether traA disruption eliminated or reduced CMIs between swarming colonies or antagonisms between strains in mixed cultures. Among 33 isolate pairs identical at traA that form clear CMIs, in no case did functional disruption of traA in one partner detectably alter CMI phenotypes. Further, traA disruption did not alleviate strong antagonisms observed during starvation-induced fruiting-body development in seven pairs of strains identical at traA. Collectively, our results suggest that most mechanisms of interference competition and inter-colony kin discrimination in natural populations of myxobacteria do not require OME. Finally, our experiments also indicate that several closely related laboratory reference strains kill some natural isolates by toxins delivered by a shared, OME-independent type VI secretion system (T6SS), suggesting that some antagonisms between sympatric natural isolates may also involve T6SS toxins.

Population Genetics in the Human Microbiome

While the human microbiome's structure and function have been extensively studied, its within-species genetic diversity is less well understood. However, genetic mutations in the microbiome can confer biomedically relevant traits, such as the ability to extract nutrients from food, metabolize drugs, evade antibiotics, and communicate with the host immune system. The population genetic processes by which these traits evolve are complex, in part due to interacting ecological and evolutionary forces in the microbiome. Advances in metagenomic sequencing, coupled with bioinformatics tools and population genetic models, facilitate quantification of microbiome genetic variation and inferences about how this diversity arises, evolves, and correlates with traits of both microbes and hosts. In this review, we explore the population genetic forces (mutation, recombination, drift, and selection) that shape microbiome genetic diversity within and between hosts, as well as efforts towards predictive models that leverage microbiome genetics.

Self-identity barcodes encoded by six expansive polymorphic toxin families discriminate kin in myxobacteria

Myxobacteria are an example of how single-cell individuals can transition into multicellular life by an aggregation strategy. For these and all organisms that consist of social groups of cells, discrimination against, and exclusion of, nonself is critical. In myxobacteria, TraA is a polymorphic cell surface receptor that identifies kin by homotypic binding, and in so doing exchanges outer membrane (OM) proteins and lipids between cells with compatible receptors. However, TraA variability alone is not sufficient to discriminate against all cells, as traA allele diversity is not necessarily high among local strains. To increase discrimination ability, myxobacteria include polymorphic OM lipoprotein toxins called SitA in their delivered cargo, which poison recipient cells that lack the cognate, allele-specific SitI immunity protein. We previously characterized 3 SitAI toxin/immunity pairs that belong to 2 families. Here, we discover 4 additional SitA families. Each family is unique in sequence, but share the characteristic features of SitA: OM-associated toxins delivered by TraA. We demonstrate that, within a SitA family, C-terminal nuclease domains are polymorphic and often modular. Remarkably, sitA loci are strikingly numerous and diverse, with most genomes possessing >30 and up to 83 distinct sitAI loci. Interestingly, all SitA protein families are serially transferred between cells, allowing a SitA inhibitor cell to poison multiple targets, including cells that never made direct contact. The expansive suites of sitAI loci thus serve as identify barcodes to exquisitely discriminate against nonself to ensure populations are genetically homogenous to conduct cooperative behaviors.

Maintenance of Sympatric and Allopatric Populations in Free-Living Terrestrial Bacteria

Due to the promiscuous exchange of genetic material and asexual reproduction, delineating microbial species (and, by extension, populations) remains challenging. Because of this, the vast majority of microbial studies assessing population structure often compare divergent strains from disparate environments under varied selective pressures. Here, we investigated the population structure within a single bacterial ecotype, a unit equivalent to a eukaryotic species, defined as highly clustered genotypic and phenotypic strains with the same ecological niche. Using a combination of genomic and computational analyses, we assessed the phylogenetic structure, extent of recombination, and flexible gene content of this genomic diversity to infer patterns of gene flow. To our knowledge, this study is the first to do so for a dominant soil bacterium. Our results indicate that bacterial soil populations, similarly to those in other environments, are structured by gene flow discontinuities and exhibit distributional patterns consistent with both isolation by distance and isolation by environment. Thus, both dispersal limitation and local environments contribute to the divergence among closely related soil bacteria as observed in macroorganisms.

For free-living bacteria and archaea, the equivalent of the biological species concept does not exist, creating several obstacles to the study of the processes contributing to microbial diversification. These obstacles are particularly high in soil, where high bacterial diversity inhibits the study of closely related genotypes and therefore the factors structuring microbial populations. Here, we isolated strains within a single Curtobacterium ecotype from surface soil (leaf litter) across a regional climate gradient and investigated the phylogenetic structure, recombination, and flexible gene content of this genomic diversity to infer patterns of gene flow. Our results indicate that microbial populations are delineated by gene flow discontinuities, with distinct populations cooccurring at multiple sites. Bacterial population structure was further delineated by genomic features allowing for the identification of candidate genes possibly contributing to local adaptation. These results suggest that the genetic structure within this bacterium is maintained both by ecological specialization in localized microenvironments (isolation by environment) and by dispersal limitation between geographic locations (isolation by distance).

Sexual Selection in Bacteria?

A main mechanism of lateral gene transfer in bacteria is transformation, where cells take up free DNA from the environment which subsequently can be recombined into the genome. Bacteria are also known to actively release DNA into the environment through secretion or lysis, which could aid uptake via transformation. Various evolutionary benefits of DNA uptake and DNA release have been proposed but these have all been framed in the context of natural selection. Here, we interpret bacterial DNA uptake and release in the context of sexual selection theory, which has been central to our understanding of the bewildering diversity of traits associated with sexual reproduction in the eukaryote world but has never been applied to prokaryotes. Specifically, we explore potential scenarios where bacteria releasing DNA into the environment could compete for successful uptake by other cells, or where bacteria could selectively take up DNA to enhance their fitness. We conclude that there is potential for sexual selection to act in bacteria, and that this might in part explain the considerable diversity in transformation-related behaviours.

Massive Gene Flux Drives Genome Diversity between Sympatric Streptomyces Conspecifics

Horizontal gene transfer is a rapid and efficient way to diversify bacterial gene pools. Currently, little is known about this gene flux within natural soil populations. Using comparative genomics of Streptomyces strains belonging to the same species and isolated at microscale, we reveal frequent transfer of a significant fraction of the pangenome. We show that it occurs at a time scale enabling the population to diversify and to cope with its changing environment, notably, through the production of public goods.

In this work, by comparing genomes of closely related individuals of Streptomyces isolated at a spatial microscale (millimeters or centimeters), we investigated the extent and impact of horizontal gene transfer in the diversification of a natural Streptomyces population. We show that despite these conspecific strains sharing a recent common ancestor, all harbored significantly different gene contents, implying massive and rapid gene flux. The accessory genome of the strains was distributed across insertion/deletion events (indels) ranging from one to several hundreds of genes. Indels were preferentially located in the arms of the linear chromosomes (ca. 12 Mb) and appeared to form recombination hot spots. Some of them harbored biosynthetic gene clusters (BGCs) whose products confer an inhibitory capacity and may constitute public goods that can favor the cohesiveness of the bacterial population. Moreover, a significant proportion of these variable genes were either plasmid borne or harbored signatures of actinomycete integrative and conjugative elements (AICEs). We propose that conjugation is the main driver for the indel flux and diversity in Streptomyces populations.

Multifactorial control of the expression of a CRISPR-Cas system by an extracytoplasmic function σ/anti-σ pair and a global regulatory complex

Expression of CRISPR-Cas systems is a prerequisite for their defensive role against invading genetic elements. Yet, much remains unknown about how this crucial step is regulated. We describe a new mechanism controlling CRISPR-cas expression, which requires an extracytoplasmic function (ECF) σ factor (DdvS), its membrane-bound anti-σ (DdvA) and a global regulatory complex (CarD–CarG). Transcriptomic analyses revealed that the DdvS/CarD/CarG-dependent regulon comprises a type III-B CRISPR-Cas system in Myxococcus xanthus. We mapped four DdvS-driven CarD/CarG-dependent promoters, with one lying immediately upstream of the cas cluster. Consistent with direct action, DdvS and CarD–CarG localize at these promoters in vivo. The cas genes are transcribed as a polycistronic mRNA that reads through the leader into the CRISPR array, a putative σ^A-dependent promoter in the leader having negligible activity in vivo. Consequently, expression of the entire CRISPR-Cas system and mature CRISPR-RNA (crRNA) production is DdvS/CarD/CarG-dependent. DdvA likely uses its large C-terminal domain to sense and transduce the extracytoplasmic signal triggering CRISPR-cas expression, which we show is not starvation-induced multicellular development. An ECF-σ/anti-σ pair and a global regulatory complex provide an effective mechanism to coordinate signal-sensing with production of precursor crRNA, its processing Cas6 endoribonuclease and other Cas proteins for mature crRNA biogenesis and interference.

Diversity of Myxobacteria-We Only See the Tip of the Iceberg

The discovery of new antibiotics is mandatory with regard to the increasing number of resistant pathogens. One approach is the search for new antibiotic producers in nature. Among actinomycetes, Bacillus species, and fungi, myxobacteria have been a rich source for bioactive secondary metabolites for decades. To date, about 600 substances could be described, many of them with antibacterial, antifungal, or cytostatic activity. But, recent cultivation-independent studies on marine, terrestrial, or uncommon habitats unequivocally demonstrate that the number of uncultured myxobacteria is much higher than would be expected from the number of cultivated strains. Although several highly promising myxobacterial taxa have been identified recently, this so-called Great Plate Count Anomaly must be overcome to get broader access to new secondary metabolite producers. In the last years it turned out that especially new species, genera, and families of myxobacteria are promising sources for new bioactive metabolites. Therefore, the cultivation of the hitherto uncultivable ones is our biggest challenge.

Multifactorial control of the expression of a CRISPR-Cas system by an extracytoplasmic function σ/anti-σ pair and a global regulatory complex

Expression of CRISPR-Cas systems is a prerequisite for their defensive role against invading genetic elements. Yet, much remains unknown about how this crucial step is regulated. We describe a new mechanism controlling CRISPR-cas expression, which requires an extracytoplasmic function (ECF) σ factor (DdvS), its membrane-bound anti-σ (DdvA) and a global regulatory complex (CarD-CarG). Transcriptomic analyses revealed that the DdvS/CarD/CarG-dependent regulon comprises a type III-B CRISPR-Cas system in Myxococcus xanthus. We mapped four DdvS-driven CarD/CarG-dependent promoters, with one lying immediately upstream of the cas cluster. Consistent with direct action, DdvS and CarD-CarG localize at these promoters in vivo. The cas genes are transcribed as a polycistronic mRNA that reads through the leader into the CRISPR array, a putative σA-dependent promoter in the leader having negligible activity in vivo. Consequently, expression of the entire CRISPR-Cas system and mature CRISPR-RNA (crRNA) production is DdvS/CarD/CarG-dependent. DdvA likely uses its large C-terminal domain to sense and transduce the extracytoplasmic signal triggering CRISPR-cas expression, which we show is not starvation-induced multicellular development. An ECF-σ/anti-σ pair and a global regulatory complex provide an effective mechanism to coordinate signal-sensing with production of precursor crRNA, its processing Cas6 endoribonuclease and other Cas proteins for mature crRNA biogenesis and interference.

Bacteriophages of Myxococcus xanthus, a Social Bacterium

Bacteriophages have been used as molecular tools in fundamental biology investigations for decades. Beyond this, however, they play a crucial role in the eco-evolutionary dynamics of bacterial communities through their demographic impact and the source of genetic information they represent. The increasing interest in describing ecological and evolutionary aspects of bacteria–phage interactions has led to major insights into their fundamental characteristics, including arms race dynamics and acquired bacterial immunity. Here, we review knowledge on the phages of the myxobacteria with a major focus on phages infecting Myxococcus xanthus, a bacterial model system widely used to study developmental biology and social evolution. In particular, we focus upon the isolation of myxophages from natural sources and describe the morphology and life cycle parameters, as well as the molecular genetics and genomics of the major groups of myxophages. Finally, we propose several interesting research directions which focus on the interplay between myxobacterial host sociality and bacteria–phage interactions.

Kin discrimination and outer membrane exchange in Myxococcus xanthus: A comparative analysis among natural isolates

Genetically similar cells of the soil bacterium Myxococcus xanthus cooperate at multiple social behaviours, including motility and multicellular development. Another social interaction in this species is outer membrane exchange (OME), a behaviour of unknown primary benefit in which cells displaying closely related variants of the outer membrane protein TraA transiently fuse and exchange membrane contents. Functionally incompatible TraA variants do not mediate OME, which led to the proposal that TraA incompatibilities determine patterns of intercellular cooperation in nature, but how this might occur remains unclear. Using natural isolates from a centimetre-scale patch of soil, we analyse patterns of TraA diversity and ask whether relatedness at TraA is causally related to patterns of kin discrimination in the form of both colony-merger incompatibilities (CMIs) and interstrain antagonisms. A large proportion of TraA functional diversity documented among global isolates is predicted to be contained within this cm-scale population. We find evidence of balancing selection on the highly variable PA14-portion of TraA and extensive transfer of traA alleles across genomic backgrounds. CMIs are shown to be common among strains identical at TraA, suggesting that CMIs are not generally caused by TraA dissimilarity. Finally, it has been proposed that interstrain antagonisms might be caused by OME-mediated toxin transfer. However, we predict that most strain pairs previously shown to exhibit strong antagonisms are incapable of OME due to TraA dissimilarity. Overall, our results suggest that most documented patterns of kin discrimination in a natural population of M. xanthus are not causally related to the TraA sequences of interactants.

Competitive Interactions Between Incompatible Mutants of the Social Bacterium Myxococcus xanthus DK1622

Due to the high similarity in their requirements for space and food, close bacterial relatives may be each other's strongest competitors. Close bacterial relatives often form visible boundaries to separate their swarming colonies, a phenomenon termed colony-merger incompatibility. While bacterial species are known to have many incompatible strains, it is largely unclear which traits lead to multiple incompatibilities and the interactions between multiple incompatible siblings. To investigate the competitive interactions of closely related incompatible strains, we mutated Myxococcus xanthus DK1622, a predatory bacterium with complex social behavior. From 3392 random transposon mutations, we obtained 11 self-identification (SI) deficient mutants that formed unmerged colony boundaries with the ancestral strain. The mutations were at nine loci with unknown functions and formed nine independent SI mutants. Compared with their ancestral strain, most of the SI mutants showed reduced growth, swarming and development abilities, but some remained unchanged from their monocultures. When pairwise mixed with their ancestral strain for co-cultivation, these mutants exhibited improved, reduced or unchanged competitive abilities compared with the ancestral strain. The sporulation efficiencies were affected by the DK1622 partner, ranging from almost complete inhibition to 360% stimulation. The differences in competitive growth between the SI mutants and DK1622 were highly correlated with the differences in their sporulation efficiencies. However, the competitive efficiencies of the mutants in mixture were inconsistent with their growth or sporulation abilities in monocultures. We propose that the colony-merger incompatibility in M. xanthus is associated with multiple independent genetic loci, and the incompatible strains hold competitive interaction abilities, which probably determine the complex relationships between multiple incompatible M. xanthus strains and their co-existence strategies.

Mechanism of Kin-Discriminatory Demarcation Line Formation between Colonies of Swarming Bacteria

Swarming bacteria use kin discrimination to preferentially associate with their clonemates for certain cooperative behaviors. Kin discrimination can manifest as an apparent demarcation line (a region lacking cells or with much lower cell density) between antagonist strains swarming toward each other. In contrast, two identical strains merge with no demarcation. Experimental studies suggest contact-dependent killing between different strains as a mechanism of kin discrimination, but it is not clear whether this killing is sufficient to explain the observed patterns. Here, we investigate the formation of demarcation line with a mathematical model. First, using data from competition experiments between kin discriminating strains of Myxococcus xanthus and Proteus mirabilis, we found the rates of killing between the strains to be highly asymmetric, i.e., one strain kills another at a much higher rate. Then, to investigate how such asymmetric interactions can lead to a stable demarcation line, we construct reaction-diffusion models for colony expansion of kin-discriminatory strains. Our results demonstrate that a stable demarcation line can form when both cell movement and cell growth cease at low nutrient levels. Further, our study suggests that, depending on the initial separation between the inoculated colonies, the demarcation line may move transiently before stabilizing. We validated these model predictions by observing dynamics of merger between two M. xanthus strains, where one strain expresses a toxin protein that kills a second strain lacking the corresponding antitoxin. Our study therefore provides a theoretical understanding of demarcation line formation between kin-discriminatory populations, and can be used for analyzing and designing future experiments.

Multilocus sequence analysis of homologous recombination and diversity in Arthrobacter sensu lato named species and glacier-inhabiting strains

Members of the bacterial genus Arthrobacter sensu lato are Gram-positive actinomycetes distributed worldwide and found in numerous environments including soil, water, glacier ice, and sewage. Homologous recombination is an important driving force in bacterial evolution, but its impact on Arthrobacter sensu lato evolution is poorly understood. We evaluated homologous recombination among 41 Arthrobacter sensu lato named species, using multilocus sequence analysis (MLSA). A high level of recombination was found, associated with strong diversification and a reticulate evolutionary pattern of Arthrobacter sensu lato. We also collected a total of 31 cold-adapted Arthrobacter sensu lato strains from two cold glaciers located in northwest China and two temperate glaciers in southwest China, and evaluated their diversity and population structure by MLSA. The glacier strains displayed high diversity, but rates of recombination among the four glacier groups were quite low, indicating that barriers to homologous recombination formed in the past among the populations on different glaciers. Our findings indicate that historical glaciation events shaped the contemporary distributions, taxonomic relationships, and phylogeographic patterns of Arthrobacter sensu lato species on glaciers.

Prokaryote genome fluidity is dependent on effective population size

Many prokaryote species are known to have fluid genomes, with different strains varying markedly in accessory gene content through the combined action of gene loss, gene gain via lateral transfer, as well as gene duplication. However, the evolutionary forces determining genome fluidity are not yet well understood. We here for the first time systematically analyse the degree to which this distinctive genomic feature differs between bacterial species. We find that genome fluidity is positively correlated with synonymous nucleotide diversity of the core genome, a measure of effective population size N_e. No effects of genome size, phylogeny or homologous recombination rate on genome fluidity were found. Our findings are consistent with a scenario where accessory gene content turnover is for a large part dictated by neutral evolution.

Self-identity reprogrammed by a single residue switch in a cell surface receptor of a social bacterium

The ability to recognize close kin confers survival benefits on single-celled microbes that live in complex and changing environments. Microbial kinship detection relies on perceptible cues that reflect relatedness between individuals, although the mechanisms underlying recognition in natural populations remain poorly understood. In myxobacteria, cells identify related individuals through a polymorphic cell surface receptor, TraA. Recognition of compatible receptors leads to outer membrane exchange among clonemates and fitness consequences. Here, we investigated how a single receptor creates a diversity in recognition across myxobacterial populations. We first show that TraA requires its partner protein TraB to function in cell-cell adhesion. Recognition is shown to be traA allele-specific, where polymorphisms within TraA dictate binding selectivity. We reveal the malleability of TraA recognition, and seemingly minor changes to its variable region reprogram recognition outcomes. Strikingly, we identify a single residue (A/P205) as a molecular switch for TraA recognition. Substitutions at this position change the specificity of a diverse panel of environmental TraA receptors. In addition, we engineered a receptor with unique specificity by simply creating an A205P substitution, suggesting that modest changes in TraA can lead to diversification of new recognition groups in nature. We hypothesize that the malleable property of TraA has allowed it to evolve and create social barriers between myxobacterial populations and in turn avoid adverse interactions with relatives.

Infectious polymorphic toxins delivered by outer membrane exchange discriminate kin in myxobacteria

Myxobacteria are known for complex social behaviors including outer membrane exchange (OME), in which cells exchange large amounts of outer membrane lipids and proteins upon contact. The TraA cell surface receptor selects OME partners based on a variable domain. However, traA polymorphism alone is not sufficient to precisely discriminate kin. Here, we report a novel family of OME-delivered toxins that promote kin discrimination of OME partners. These SitA lipoprotein toxins are polymorphic and widespread in myxobacteria. Each sitA is associated with a cognate sitI immunity gene, and in some cases a sitB accessory gene. Remarkably, we show that SitA is transferred serially between target cells, allowing the toxins to move cell-to-cell like an infectious agent. Consequently, SitA toxins define strong identity barriers between strains and likely contribute to population structure, maintenance of cooperation, and strain diversification. Moreover, these results highlight the diversity of systems evolved to deliver toxins between bacteria.

Evolution by flight and fight: diverse mechanisms of adaptation by actively motile microbes

Evolutionary adaptation can be achieved by mechanisms accessible to all organisms, including faster growth and interference competition, but self-generated motility offers additional possibilities. We tested whether 55 populations of the bacterium Myxococcus xanthus that underwent selection for increased fitness at the leading edge of swarming colonies adapted by swarming faster toward unused resources or by other means. Populations adapted greatly but diversified markedly in both swarming phenotypes and apparent mechanisms of adaptation. Intriguingly, although many adapted populations swarm intrinsically faster than their ancestors, numerous others do not. Some populations evolved interference competition toward their ancestors, whereas others gained the ability to facultatively increase swarming rate specifically upon direct interaction with ancestral competitors. Our results both highlight the diverse range of mechanisms by which actively motile organisms can adapt evolutionarily and help to explain the high levels of swarming-phenotype diversity found in local soil populations of M. xanthus.

Kin Recognition in Bacteria

The ability of bacteria to recognize kin provides a means to form social groups. In turn these groups can lead to cooperative behaviors that surpass the ability of the individual. Kin recognition involves specific biochemical interactions between a receptor(s) and an identification molecule(s). Recognition specificity, ensuring that nonkin are excluded and kin are included, is critical and depends on the number of loci and polymorphisms involved. After recognition and biochemical perception, the common ensuing cooperative behaviors include biofilm formation, quorum responses, development, and swarming motility. Although kin recognition is a fundamental mechanism through which cells might interact, microbiologists are only beginning to explore the topic. This review considers both molecular and theoretical aspects of bacterial kin recognition. Consideration is also given to bacterial diversity, genetic relatedness, kin selection theory, and mechanisms of recognition.