Identifying the Gene Regulatory Network of the Starvation-Induced Transcriptional Activator Nla28.

Myxococcus xanthus copes with starvation by producing fruiting bodies filled with dormant and stress-resistant spores. Here, we aimed to better define the gene regulatory network associated with Nla28, a transcriptional activator/enhancer binding protein (EBP) and a key regulator of the early starvation response. Previous work showed that Nla28 directly regulates EBP genes that are important for fruiting body development. However, the Nla28 regulatory network is likely to be much larger because hundreds of starvation-induced genes are downregulated in a nla28 mutant strain. To identify candidates for direct Nla28-mediated transcription, we analyzed the downregulated genes using a bioinformatics approach. Nine potential Nla28 target promoters (29 genes) were discovered. The results of in vitro promoter binding assays, coupled with in vitro and in vivo mutational analyses, suggested that the nine promoters along with three previously identified EBP gene promoters were indeed in vivo targets of Nla28. These results also suggested that Nla28 used tandem, imperfect repeats of an 8-bp sequence for promoter binding. Interestingly, eight of the new Nla28 target promoters were predicted to be intragenic. Based on mutational analyses, the newly identified Nla28 target loci contained at least one gene that was important for starvation-induced development. Most of these loci contained genes predicted to be involved in metabolic or defense-related functions. Using the consensus Nla28 binding sequence, bioinformatics, and expression profiling, 58 additional promoters and 102 genes were tagged as potential Nla28 targets. Among these putative Nla28 targets, functions, such as regulatory, metabolic, and cell envelope biogenesis, were assigned to many genes. IMPORTANCE In bacteria, starvation leads to profound changes in behavior and physiology. Some of these changes have economic and health implications because the starvation response has been linked to the formation of biofilms, virulence, and antibiotic resistance. To better understand how starvation contributes to changes in bacterial physiology and resistance, we identified the putative starvation-induced gene regulatory network associated with Nla28, a transcriptional activator from the bacterium Myxoccocus xanthus. We determined the mechanism by which starvation-responsive genes were activated by Nla28 and showed that several of the genes were important for the formation of a highly resistant cell type.

Profiling Myxococcus xanthus Swarming Phenotypes through Mutation and Environmental Variation.

Myxococcus xanthus is a bacterium that lives on surfaces as a predatory biofilm called a swarm. As a growing swarm feeds on prey and expands, it displays dynamic multicellular patterns such as traveling waves called ripples and branching protrusions called flares. The rate at which a swarm expands across a surface, and the emergence of the coexisting patterns, are all controlled through coordinated cell movement. M. xanthus cells move using two motility systems known as adventurous (A) and social (S). Both are involved in swarm expansion and pattern formation. In this study, we describe a set of M. xanthus swarming genotype-to-phenotype associations that include both genetic and environmental perturbations. We identified new features of the swarming phenotype, recorded and measured swarm expansion using time-lapse microscopy, and compared the impact of mutations on different surfaces. These observations and analyses have increased our ability to discriminate between swarming phenotypes and provided context that allows us to identify some phenotypes as improbable outliers within the M. xanthus swarming phenome. IMPORTANCE Myxococcus xanthus grows on surfaces as a predatory biofilm called a swarm. In nature, a feeding swarm expands by moving over and consuming prey bacteria. In the laboratory, a swarm is created by spotting cell suspension onto nutrient agar in lieu of prey. The suspended cells quickly settle on the surface as the liquid is absorbed into the agar, and the new swarm then expands radially. An assay that measures the expansion rate of a swarm of mutant cells is the first, and sometimes only, measurement used to decide whether a particular mutation impacts swarm motility. We have broadened the scope of this assay by increasing the accuracy of measurements and introducing prey, resulting in new identifiable and quantifiable features that can be used to improve genotype-to-phenotype associations.

Self-Driven Phase Transitions Drive Myxococcus xanthus Fruiting Body Formation.

Combining high-resolution single cell tracking experiments with numerical simulations, we show that starvation-induced fruiting body formation in Myxococcus xanthus is a phase separation driven by cells that tune their motility over time. The phase separation can be understood in terms of cell density and a dimensionless Péclet number that captures cell motility through speed and reversal frequency. Our work suggests that M. xanthus takes advantage of a self-driven nonequilibrium phase transition that can be controlled at the single cell level.

Xanthusbase after five years expands to become Openmods.

Xanthusbase (http://www.xanthusbase.org), a model organism database for the bacterium Myxococcus xanthus, functions as a collaborative information repository based on Wikipedia principles. It was created more than 5 years ago to serve as a cost-effective reference database for M. xanthus researchers, an education tool for undergraduate students to learn about genome annotation, and a means for the community of researchers to collaboratively improve their organism's annotation. We have achieved several goals and are seeking creative solutions to ongoing challenges. Along the way we have made several important improvements to Xanthusbase related to stability, security and usability. Most importantly, we have designed and implemented an installer that enables other microbial model organism communities to use it as a MOD. This version, called Openmods, has already been used to create Xenorhabdusbase (http://xenorhabdusbase.bact.wisc.edu), Caulobacterbase (http://caulobacterbase.bsd.uchicago.edu) and soon Bdellovibriobase.

Reconciliation of sequence data and updated annotation of the genome of Agrobacterium tumefaciens C58, and distribution of a linear chromosome in the genus Agrobacterium.

Two groups independently sequenced the Agrobacterium tumefaciens C58 genome in 2001. We report here consolidation of these sequences, updated annotation, and additional analysis of the evolutionary history of the linear chromosome, which is apparently limited to the biovar I group of Agrobacterium.

Cell behaviors underlying Myxococcus xanthus aggregate dispersal.

IMPORTANCE: Understanding the processes behind bacterial biofilm formation, maintenance, and dispersal is essential for addressing their effects on health and ecology. Within these multicellular communities, various cues can trigger differentiation into distinct cell types, allowing cells to adapt to their specific local environment. The soil bacterium Myxococcus xanthus forms biofilms in response to starvation, marked by cells aggregating into mounds. Some aggregates persist as spore-filled fruiting bodies, while others disperse after initial formation for unknown reasons. Here, we use a combination of cell tracking analysis and computational simulations to identify behaviors at the cellular level that contribute to aggregate dispersal. Our results suggest that cells in aggregates actively determine whether to disperse or persist and undergo a transition to sporulation based on a self-produced cue related to the aggregate size. Identifying these cues is an important step in understanding and potentially manipulating bacterial cell-fate decisions.

Mechanobiology as a tool for addressing the genotype-to-phenotype problem in microbiology.

The central hypothesis of the genotype-phenotype relationship is that the phenotype of a developing organism (i.e., its set of observable attributes) depends on its genome and the environment. However, as we learn more about the genetics and biochemistry of living systems, our understanding does not fully extend to the complex multiscale nature of how cells move, interact, and organize; this gap in understanding is referred to as the genotype-to-phenotype problem. The physics of soft matter sets the background on which living organisms evolved, and the cell environment is a strong determinant of cell phenotype. This inevitably leads to challenges as the full function of many genes, and the diversity of cellular behaviors cannot be assessed without wide screens of environmental conditions. Cellular mechanobiology is an emerging field that provides methodologies to understand how cells integrate chemical and physical environmental stress and signals, and how they are transduced to control cell function. Biofilm forming bacteria represent an attractive model because they are fast growing, genetically malleable and can display sophisticated self-organizing developmental behaviors similar to those found in higher organisms. Here, we propose mechanobiology as a new area of study in prokaryotic systems and describe its potential for unveiling new links between an organism's genome and phenome.

A Clp/Hsp100 chaperone functions in Myxococcus xanthus sporulation and self-organization.

The Clp/Hsp100 proteins are chaperones that play a role in protein degradation and reactivation. In bacteria, they exhibit a high degree of pleiotropy, affecting both individual and multicellular phenotypes. In this article, we present the first characterization of a Clp/Hsp100 homolog in Myxococcus xanthus (MXAN_4832 gene locus). Deletion of MXAN_4832 causes defects in both swarming and aggregation related to cell motility and the production of fibrils, which are an important component of the extracellular matrix of a swarm. The deletion also affects the formation of myxospores during development, causing them to become sensitive to heat. The protein product of MXAN_4832 can act as a chaperone in vitro, providing biochemical evidence in support of our hypothesis that MXAN_4832 is a functional Clp/Hsp100 homolog. There are a total of 12 Clp/Hsp100 homologs in M. xanthus, including MXAN_4832, and, based on its mutational and biochemical characterization, they may well represent an important group.

Spreading rates of bacterial colonies depend on substrate stiffness and permeability.

The ability of bacteria to colonize and grow on different surfaces is an essential process for biofilm development. Here, we report the use of synthetic hydrogels with tunable stiffness and porosity to assess physical effects of the substrate on biofilm development. Using time-lapse microscopy to track the growth of expanding Serratia marcescens colonies, we find that biofilm colony growth can increase with increasing substrate stiffness, unlike what is found on traditional agar substrates. Using traction force microscopy-based techniques, we find that biofilms exert transient stresses correlated over length scales much larger than a single bacterium, and that the magnitude of these forces also increases with increasing substrate stiffness. Our results are consistent with a model of biofilm development in which the interplay between osmotic pressure arising from the biofilm and the poroelastic response of the underlying substrate controls biofilm growth and morphology.

Quantifying aggregation dynamics during Myxococcus xanthus development.

Under starvation conditions, a swarm of Myxococcus xanthus cells will undergo development, a multicellular process culminating in the formation of many aggregates called fruiting bodies, each of which contains up to 100,000 spores. The mechanics of symmetry breaking and the self-organization of cells into fruiting bodies is an active area of research. Here we use microcinematography and automated image processing to quantify several transient features of developmental dynamics. An analysis of experimental data indicates that aggregation reaches its steady state in a highly nonmonotonic fashion. The number of aggregates rapidly peaks at a value 2- to 3-fold higher than the final value and then decreases before reaching a steady state. The time dependence of aggregate size is also nonmonotonic, but to a lesser extent: average aggregate size increases from the onset of aggregation to between 10 and 15 h and then gradually decreases thereafter. During this process, the distribution of aggregates transitions from a nearly random state early in development to a more ordered state later in development. A comparison of experimental results to a mathematical model based on the traffic jam hypothesis indicates that the model fails to reproduce these dynamic features of aggregation, even though it accurately describes its final outcome. The dynamic features of M. xanthus aggregation uncovered in this study impose severe constraints on its underlying mechanisms.

The σ54 system directly regulates bacterial natural product genes.

Bacterial-derived polyketide and non-ribosomal peptide natural products are crucial sources of therapeutics and yet little is known about the conditions that favor activation of natural product genes or the regulatory machinery controlling their transcription. Recent findings suggest that the σ54 system, which includes σ54-loaded RNA polymerase and transcriptional activators called enhancer binding proteins (EBPs), might be a common regulator of natural product genes. Here, we explored this idea by analyzing a selected group of putative σ54 promoters identified in Myxococcus xanthus natural product gene clusters. We show that mutations in putative σ54-RNA polymerase binding regions and in putative Nla28 EBP binding sites dramatically reduce in vivo promoter activities in growing and developing cells. We also show in vivo promoter activities are reduced in a nla28 mutant, that Nla28 binds to wild-type fragments of these promoters in vitro, and that in vitro binding is lost when the Nla28 binding sites are mutated. Together, our results indicate that M. xanthus uses σ54 promoters for transcription of at least some of its natural product genes. Interestingly, the vast majority of experimentally confirmed and putative σ54 promoters in M. xanthus natural product loci are located within genes and not in intergenic sequences.

Quantification of Myxococcus xanthus Aggregation and Rippling Behaviors: Deep-Learning Transformation of Phase-Contrast into Fluorescence Microscopy Images.

Myxococcus xanthus bacteria are a model system for understanding pattern formation and collective cell behaviors. When starving, cells aggregate into fruiting bodies to form metabolically inert spores. During predation, cells self-organize into traveling cell-density waves termed ripples. Both phase-contrast and fluorescence microscopy are used to observe these patterns but each has its limitations. Phase-contrast images have higher contrast, but the resulting image intensities lose their correlation with cell density. The intensities of fluorescence microscopy images, on the other hand, are well-correlated with cell density, enabling better segmentation of aggregates and better visualization of streaming patterns in between aggregates; however, fluorescence microscopy requires the engineering of cells to express fluorescent proteins and can be phototoxic to cells. To combine the advantages of both imaging methodologies, we develop a generative adversarial network that converts phase-contrast into synthesized fluorescent images. By including an additional histogram-equalized output to the state-of-the-art pix2pixHD algorithm, our model generates accurate images of aggregates and streams, enabling the estimation of aggregate positions and sizes, but with small shifts of their boundaries. Further training on ripple patterns enables accurate estimation of the rippling wavelength. Our methods are thus applicable for many other phenotypic behaviors and pattern formation studies.

Genome sequences of three agrobacterium biovars help elucidate the evolution of multichromosome genomes in bacteria.

The family Rhizobiaceae contains plant-associated bacteria with critical roles in ecology and agriculture. Within this family, many Rhizobium and Sinorhizobium strains are nitrogen-fixing plant mutualists, while many strains designated as Agrobacterium are plant pathogens. These contrasting lifestyles are primarily dependent on the transmissible plasmids each strain harbors. Members of the Rhizobiaceae also have diverse genome architectures that include single chromosomes, multiple chromosomes, and plasmids of various sizes. Agrobacterium strains have been divided into three biovars, based on physiological and biochemical properties. The genome of a biovar I strain, A. tumefaciens C58, has been previously sequenced. In this study, the genomes of the biovar II strain A. radiobacter K84, a commercially available biological control strain that inhibits certain pathogenic agrobacteria, and the biovar III strain A. vitis S4, a narrow-host-range strain that infects grapes and invokes a hypersensitive response on nonhost plants, were fully sequenced and annotated. Comparison with other sequenced members of the Alphaproteobacteria provides new data on the evolution of multipartite bacterial genomes. Primary chromosomes show extensive conservation of both gene content and order. In contrast, secondary chromosomes share smaller percentages of genes, and conserved gene order is restricted to short blocks. We propose that secondary chromosomes originated from an ancestral plasmid to which genes have been transferred from a progenitor primary chromosome. Similar patterns are observed in select Beta- and Gammaproteobacteria species. Together, these results define the evolution of chromosome architecture and gene content among the Rhizobiaceae and support a generalized mechanism for second-chromosome formation among bacteria.

Xanthusbase: adapting wikipedia principles to a model organism database.

xanthusBase (http://www.xanthusbase.org) is the official model organism database (MOD) for the social bacterium Myxococcus xanthus. In many respects, M.xanthus represents the pioneer model organism (MO) for studying the genetic, biochemical, and mechanistic basis of prokaryotic multicellularity, a topic that has garnered considerable attention due to the significance of biofilms in both basic and applied microbiology research. To facilitate its utility, the design of xanthusBase incorporates open-source software, leveraging the cumulative experience made available through the Generic Model Organism Database (GMOD) project, MediaWiki (http://www.mediawiki.org), and dictyBase (http://www.dictybase.org), to create a MOD that is both highly useful and easily navigable. In addition, we have incorporated a unique Wikipedia-style curation model which exploits the internet's inherent interactivity, thus enabling M.xanthus and other myxobacterial researchers to contribute directly toward the ongoing genome annotation.

Chemotaxis as an emergent property of a swarm.

We have characterized and quantified a form of bacterial chemotaxis that manifests only as an emergent property by measuring symmetry breaking in a swarm of Myxococcus xanthus exposed to a two-dimensional nutrient gradient from within an agar substrate. M. xanthus chemotaxis requires cell-cell contact and coordinated motility, as individual motile cells exhibit only nonvectorial movement in the presence of a nutrient gradient. Genes that specifically affect M. xanthus chemotaxis include at least 10 of the 53 that express enhancer binding proteins of the NtrC-like class, an indication that this behavior is controlled through transcription, most likely by a complex signal transduction network.

Functional genome annotation through phylogenomic mapping.

Accurate determination of functional interactions among proteins at the genome level remains a challenge for genomic research. Here we introduce a genome-scale approach to functional protein annotation--phylogenomic mapping--that requires only sequence data, can be applied equally well to both finished and unfinished genomes, and can be extended beyond single genomes to annotate multiple genomes simultaneously. We have developed and applied it to more than 200 sequenced bacterial genomes. Proteins with similar evolutionary histories were grouped together, placed on a three dimensional map and visualized as a topographical landscape. The resulting phylogenomic maps display thousands of proteins clustered in mountains on the basis of coinheritance, a strong indicator of shared function. In addition to systematic computational validation, we have experimentally confirmed the ability of phylogenomic maps to predict both mutant phenotype and gene function in the delta proteobacterium Myxococcus xanthus.

Spatial organization of Myxococcus xanthus during fruiting body formation.

Microcinematography was used to examine fruiting body development of Myxococcus xanthus. Wild-type cells progress through three distinct phases: a quiescent phase with some motility but little aggregation (0 to 8 h), a period of vigorous motility leading to raised fruiting bodies (8 to 16 h), and a period of maturation during which sporulation is initiated (16 to 48 h). Fruiting bodies are extended vertically in a series of tiers, each involving the addition of a cell monolayer on top of the uppermost layer. A pilA (MXAN_5783) mutant produced less extracellular matrix material and thus allowed closer examination of tiered aggregate formation. A csgA (MXAN_1294) mutant exhibited no quiescent phase, aberrant aggregation in phase 2, and disintegration of the fruiting bodies in the third phase.

Inter-laboratory evolution of a model organism and its epistatic effects on mutagenesis screens.

In theory, a few naturally occurring evolutionary changes in the genome of a model organism may have little or no observable impact on its wild type phenotype, and yet still substantially impact the phenotypes of mutant strains through epistasis. To see if this is happening in a model organism, we obtained nine different laboratories' wild type Myxococcus xanthus DK1622 "sublines" and sequenced each to determine if they had evolved after their physical separation. Under a common garden experiment, each subline satisfied the phenotypic prerequisites for wild type, but many differed to a significant degree in each of the four quantitative phenotypic traits we measured, with some sublines differing by several-fold. Genome resequencing identified 29 variants between the nine sublines, and eight had at least one unique variant within an Open Reading Frame (ORF). By disrupting the ORF MXAN7041 in two different sublines, we demonstrated substantial epistasis from these naturally occurring variants. The impact of such inter-laboratory wild type evolution is important to any genotype-to-phenotype study; an organism's phenotype may be sensitive to small changes in genetic background, so that results from phenotypic screens and other related experiments might not agree with prior published results or the results from other laboratories.

Phenotypic profiling of ABC transporter coding genes in Myxococcus xanthus.

Information about a gene sometimes can be deduced by examining the impact of its mutation on phenotype. However, the genome-scale utility of the method is limited because, for nearly all model organisms, the majority of mutations result in little or no observable phenotypic impact. The cause of this is often attributed to robustness or redundancy within the genome, but that is only one plausible hypothesis. We examined a standard set of phenotypic traits, and applied statistical methods commonly used in the study of natural variants to an engineered mutant strain collection representing disruptions in 180 of the 192 ABC transporters within the bacterium Myxococcus xanthus. These strains display continuous variation in their phenotypic distributions, with a small number of "outlier" strains at both phenotypic extremes, and the majority within a confidence interval about the mean that always includes wild type. Correlation analysis reveals substantial pleiotropy, indicating that the traits do not represent independent variables. The traits measured in this study co-cluster with expression profiles, thereby demonstrating that these changes in phenotype correspond to changes at the molecular level, and therefore can be indirectly connected to changes in the genome. However, the continuous distributions, the pleiotropy, and the placement of wild type always within the confidence interval all indicate that this standard set of M. xanthus phenotypic assays is measuring a narrow range of partially overlapping traits that do not directly reflect fitness. This is likely a significant cause of the observed small phenotypic impact from mutation, and is unrelated to robustness and redundancy.

Describing Myxococcus xanthus aggregation using Ostwald ripening equations for thin liquid films.

When starved, a swarm of millions of Myxococcus xanthus cells coordinate their movement from outward swarming to inward coalescence. The cells then execute a synchronous program of multicellular development, arranging themselves into dome shaped aggregates. Over the course of development, about half of the initial aggregates disappear, while others persist and mature into fruiting bodies. This work seeks to develop a quantitative model for aggregation that accurately simulates which will disappear and which will persist. We analyzed time-lapse movies of M. xanthus development, modeled aggregation using the equations that describe Ostwald ripening of droplets in thin liquid films, and predicted the disappearance and persistence of aggregates with an average accuracy of 85%. We then experimentally validated a prediction that is fundamental to this model by tracking individual fluorescent cells as they moved between aggregates and demonstrating that cell movement towards and away from aggregates correlates with aggregate disappearance. Describing development through this model may limit the number and type of molecular genetic signals needed to complete M. xanthus development, and it provides numerous additional testable predictions.