Mechanisms driving spatial distribution of residents in colony biofilms: an interdisciplinary perspective

Mechanical properties of Staphylococcus aureus and Pseudomonas aeruginosa dual-species biofilms grown in chronic wound-based models

Wound infections become chronic due to biofilm formation by pathogenic bacteria; two such pathogens are Staphylococcus aureus and Pseudomonas aeruginosa. These bacteria are known to form polymicrobial biofilms in wounds, which exhibit increased colonization rates, enhanced chronicity, and greater resistance to treatment. Previously, the impacts of a wound bed environment on the mechanical properties of P. aeruginosa biofilms have been explored, and in this work the role of a wound bed environment in the viscoelasticity and microstructure of polymicrobial biofilms is characterized. We hypothesize that common wound bed proteins mediate interactions between S. aureus and P. aeruginosa to enable the formation of more elastic and stiff biofilms. Growth media with varying protein content as well as additional collagen, a protein associated with a wound extracellular matrix, were utilized to test our hypothesis. Microrheology indicates that both P. aeruginosa and S. aureus form relatively stiffer single-species biofilms in a wound environment with collagen. S. aureus produced stiffer biofilms in the presence of collagen, regardless of other wound proteins, likely due to its interactions with collagen. When both species were grown together in wound-like media, synergistic effects led to stiffer dual-species biofilms compared to their single-species forms. Under all growth conditions, collagen significantly contributed to stiffening P. aeruginosa/S. aureus dual-species biofilms, suggesting that it mediates complex interspecies interactions. High-resolution imaging and analysis revealed that collagen also influenced the microstructures of P. aeruginosa/S. aureus dual-species biofilms. In media containing wound proteins and collagen, S. aureus clusters were larger and exhibited more complex shapes. These results indicate that the wound bed environment not only provides improved antibacterial resistance due to cooperative interactions, but also improved mechanical protection, which impact common treatment methods like debridement.

Disentangling the feedback loops driving spatial patterning in microbial communities

The properties of multispecies biofilms are determined by how species are arranged in space. How these patterns emerge is a complex and largely unsolved problem. Here, we synthesize the known factors affecting pattern formation, identify the interdependencies and feedback loops coupling them, and discuss approaches to disentangle their effects. Finally, we propose an interdisciplinary research program that could create a predictive understanding of pattern formation in microbial communities.

Spatial population dynamics of bacterial colonies with social antibiotic resistance

Bacteria frequently inhabit surface-attached communities where rich "social" interactions can significantly alter their population-level behavior, including their response to antibiotics. Understanding these collective effects in spatially heterogeneous communities is an ongoing challenge. Here, we investigated the spatial organization that emerges from antibiotic exposure in initially randomly distributed communities containing antibiotic-resistant and -sensitive strains of Enterococcus faecalis, an opportunistic pathogen. We identified that a range of complex spatial structures emerged in the population homeland-the inoculated region that microbes inhabit prior to range expansion-which depended on initial colony composition and antibiotic concentration. We found that these arrangements were explained by cooperative interactions between resistant and sensitive subpopulations with a variable spatial scale, the result of dynamic zones of protection afforded to sensitive cells by growing populations of enzyme-producing resistant neighbors. Using a combination of experiments and mathematical models, we explored the complex spatiotemporal interaction dynamics that create these patterns, and predicted spatial arrangements of sensitive and resistant subpopulations under new conditions. We illustrated how spatial population dynamics in the homeland affect subsequent range expansion, both because they modulate the composition of the initial expanding front, and through long-range cooperation between the homeland and the expanding region. Finally, we showed that these spatial constraints resulted in populations whose size and composition differed markedly from matched populations in well-stirred (planktonic) cultures. These findings underscore the importance of spatial structure and cooperation, long-studied features in theoretical ecology, for determining the fate of bacterial communities under antibiotic exposure.

Intrinsically robust and scalable biofilm segmentation under diverse physical growth conditions

Developmental patterning is a shared feature across biological systems ranging from vertebrates to bacterial biofilms. While vertebrate patterning benefits from well-controlled homeostatic environments, bacterial biofilms can grow in diverse physical contexts. What mechanisms provide developmental robustness under diverse environments remains an open question. We show that a native clock-and-wavefront mechanism robustly segments biofilms in both solid-air and solid-liquid interfaces. Biofilms grown under these distinct physical conditions differ 4-fold in size yet exhibit robust segmentation. The segmentation pattern scaled with biofilm growth rate in a mathematically predictable manner independent of habitat conditions. We show that scaling arises from the coupling between wavefront speed and biofilm growth rate. In contrast to the complexity of scaling mechanisms in vertebrates, our data suggests that the minimal bacterial clock-and-wavefront mechanism is intrinsically robust and scales in real time. Consequently, bacterial biofilms robustly segment under diverse conditions without requiring cell-to-cell signaling to track system size.

Growing bacterial colonies harness emergent genealogical demixing to regulate organizational entropy

Spatiotemporal organization of individuals within growing bacterial colonies is a key determinant of intraspecific interactions and colony-scale heterogeneities. The evolving cellular distribution, in relation to the genealogical lineage, is thus central to our understanding of bacterial fate across scales. Yet, how bacteria self-organize genealogically as a colony expands has remained unknown. Here, by developing a custom-built label-free algorithm, we track and study the genesis and evolution of emergent self-similar genealogical enclaves, whose dynamics are governed by biological activity. Topological defects at enclave boundaries tune finger-like morphologies of the active interfaces. The Shannon entropy of cell arrangements reduce over time; with faster-dividing cells possessing higher spatial affinity to genealogical relatives, at the cost of a well-mixed, entropically favorable state. Our coarse-grained lattice model demonstrates that genealogical enclaves emerge due to an interplay of division-mediated dispersal, stochasticity of division events, and cell-cell interactions. The study reports so-far hidden emergent self-organizing features arising due to entropic suppression, ultimately modulating intraspecific genealogical distances within bacterial colonies.

Spatial population dynamics of bacterial colonies with social antibiotic resistance

Bacteria frequently inhabit surface-attached communities where rich "social" interactions can significantly alter their population-level behavior, including their response to antibiotics. Understanding these collective effects in spatially heterogeneous communities is an ongoing challenge. Here, we investigated the spatial organization that emerges from antibiotic exposure in initially randomly distributed communities containing antibiotic-resistant and -sensitive strains of E. faecalis , an opportunistic pathogen. We identified that a range of complex spatial structures emerged in the population homeland-the inoculated region that microbes inhabit prior to range expansion-, which depended on initial colony composition and antibiotic concentration. We found that these arrangements were explained by cooperative interactions between resistant and sensitive subpopulations with a variable spatial scale, the result of dynamic zones of protection afforded to sensitive cells by growing populations of enzyme-producing resistant neighbors. Using a combination of experiments and mathematical models, we explored the complex spatiotemporal interaction dynamics that create these patterns, and predicted spatial arrangements of sensitive and resistant subpopulations under new conditions. We illustrated how spatial population dynamics in the homeland affect subsequent range expansion, both because they modulate the composition of the initial expanding front, and through long-range cooperation between the homeland and the expanding region. Finally, we showed that these spatial constraints resulted in populations whose size and composition differed markedly from matched populations in well-stirred (planktonic) cultures. These findings underscore the importance of spatial structure and cooperation, long-studied features in theoretical ecology, for determining the fate of bacterial communities under antibiotic exposure.

Significance

Interactions between bacteria are common, particularly in the crowded surface-associated communities that occur anywhere from natural ecosystems to the human body to medical devices. Antibiotic resistance can be influenced by these "social" interactions, making it difficult to predict how spatial communities respond to antibiotic. Here, we show that complex spatial arrangements emerge when initially randomly distributed populations of antibiotic-resistant and -sensitive E. faecalis , a microbial pathogen, are exposed to antibiotic. Using mathematical models and experiments, we show how local competition and dynamic-range cross-protection drive pattern formation. As a result, these spatially structured populations respond differently to antibiotics than well-mixed communities. Our findings elucidate how "social" antibiotic resistance affects spatially structured bacterial communities, a step towards predicting and controlling resistance.

Quantifying the fractal complexity of nutrient transport channels in Escherichia coli biofilms under varying cell shape and growth environment

Recent mesoscopic characterization of nutrient-transporting channels in Escherichia coli has allowed the identification and measurement of individual channels in whole mature colony biofilms. However, their complexity under different physiological and environmental conditions remains unknown. Analysis of confocal micrographs of colony biofilms formed by cell shape mutants of E. coli shows that channels have high fractal complexity, regardless of cell phenotype or growth medium. In particular, colony biofilms formed by the mutant strain ΔompR, which has a wide-cell phenotype, have a higher fractal dimension when grown on rich medium than when grown on minimal medium, with channel complexity affected by glucose and agar concentrations in the medium. Osmotic stress leads to a dramatic reduction in the ΔompR cell size but has a limited effect on channel morphology. This work shows that fractal image analysis is a powerful tool to quantify the effect of phenotypic mutations and growth environment on the morphological complexity of internal E. coli biofilm structures. If applied to a wider range of mutant strains, this approach could help elucidate the genetic determinants of channel formation in E. coli colony biofilms.

Essential genes for Haemophilus parainfluenzae survival and biofilm growth

Haemophilus parainfluenzae (Hp) is a Gram-negative, highly prevalent, and abundant commensal in the human oral cavity, and an infrequent extraoral opportunistic pathogen. Hp occupies multiple niches in the oral cavity, including the supragingival plaque biofilm. Little is known about how Hp interacts with its neighbors in healthy biofilms nor its mechanisms of pathogenesis as an opportunistic pathogen. To address this, we identified the essential genome and conditionally essential genes in in vitro biofilms aerobically and anaerobically. Using transposon insertion sequencing (TnSeq) with a highly saturated mariner transposon library in two strains, the ATCC33392 type-strain (Hp 392) and oral isolate EL1 (Hp EL1), we show that the essential genomes of Hp 392 and Hp EL1 are composed of 395 (20%) and 384 (19%) genes, respectively. The core essential genome, consisting of 341 (17%) essential genes conserved between both strains, was composed of genes associated with genetic information processing, carbohydrate, protein, and energy metabolism. We also identified conditionally essential genes for aerobic and anaerobic biofilm growth, which were associated with carbohydrate and energy metabolism in both strains. RNAseq analysis determined that most genes upregulated during anaerobic growth are not essential for Hp 392 anaerobic survival. The completion of this library and analysis under these conditions gives us a foundational insight into the basic biology of H. parainfluenzae in differing oxygen conditions, similar to its in vivo habitat. This library presents a valuable tool for investigation into conditionally essential genes for an organism that lives in close contact with many microbial species in the human oral habitat.IMPORTANCEHaemophilus parainfluenzae is a highly abundant human commensal microbe, present in most healthy individuals where it colonizes the mouth. H. parainfluenzae correlates with good oral health and may play a role in preservation of healthy host status. Also, H. parainfluenzae can cause opportunistic infections outside of the oral cavity. To date, little is known about how H. parainfluenzae colonizes the human host, despite being such a frequent and abundant part of our human microbiome. Here, we demonstrate the creation and use of a powerful tool, a TnSeq library, used to identify genes necessary for both the outright growth of this organism and also genes conditionally essential for growth in varying oxygen status which it can encounter in the human host. This tool and these data serve as a foundation for further study of this relatively unknown organism that may play a role in preserving human health.

Cyclic olefin copolymer (COC) as a promising biomaterial for affecting bacterial colonization: investigation on Vibrio campbellii

Cyclic olefin copolymer (COC) has emerged as an interesting biocompatible material for Organ-on-a-Chip (OoC) devices monitoring growth, viability, and metabolism of cells. Despite ISO 10993 approval, systematic investigation of bacteria grown onto COC is a still not documented issue. This study discusses biofilm formations of the canonical wild type BB120 Vibrio campbellii strain on a native COC substrate and addresses the impact of the physico-chemical properties of COC compared to conventional hydroxyapatite (HA) and poly(dimethylsiloxane) (PDMS) surfaces. An interdisciplinary approach combining bacterial colony counting, light microscopy imaging and advanced digital image processing remarks interesting results. First, COC can reduce biomass adhesion with respect to common biopolymers, that is suitable for tuning biofilm formations in the biological and medical areas. Second, remarkably different biofilm morphology (dendritic complex patterns only in the case of COC) was observed among the examined substrates. Third, the observed biofilm morphogenesis was related to the interaction of COC with the conditioning layer of the planktonic biological medium. Fourth, Level Co-occurrence Matrix (CGLM)-based analysis enabled quantitative assessment of the biomass textural fractal development under different coverage conditions. All of this is of key practical relevance in searching innovative biocompatible materials for pharmaceutical, implantable and medical products.

A Label-Free Approach for Relative Spatial Quantitation of c-di-GMP in Microbial Biofilms

Microbial biofilms represent an important lifestyle for bacteria and are dynamic three-dimensional structures. Cyclic dimeric guanosine monophosphate (c-di-GMP) is a ubiquitous signaling molecule that is known to be tightly regulated with biofilm processes. While measurements of global levels of c-di-GMP have proven valuable toward understanding the genetic control of c-di-GMP production, there is a need for tools to observe the local changes of c-di-GMP production in biofilm processes. We have developed a label-free method for the direct detection of c-di-GMP in microbial colony biofilms using matrix-assisted laser desorption/ionization mass spectrometry imaging (MALDI-MSI). We applied this method to the enteric pathogen Vibrio cholerae, the marine symbiont V. fischeri, and the opportunistic pathogen Pseudomonas aeruginosa PA14 and detected spatial and temporal changes in c-di-GMP signal that accompanied genetic alterations in factors that synthesize and degrade the compound. We further demonstrated how this method can be simultaneously applied to detect additional metabolites of interest from a single sample.

Resource limitation and population fluctuation drive spatiotemporal order in microbial communities

Microbial communities display complex spatiotemporal behaviors leading to spatially-structured and ordered organization driven by species interactions and environmental factors. Resource availability plays a pivotal role in shaping the dynamics of bacterial colonies. In this study, we delve into the intricate interplay between resource limitation and the emergent properties of a growing colony of two visually distinct bacterial strains having similar growth and mechanical properties. Employing an agent-based modeling and computer simulations, we analyze the resource-driven effect on segregation and sectoring, cell length regulation and nematic ordering within a growing colony. We introduce a dimensionless parameter referred to as the active layer thickness, derived from nutrient diffusion equations, indicating effective population participation due to local resource availability. Our results reveal that lower values of active layer thickness arising from decreased resource abundance lead to rougher colony fronts, fostering heightened population fluctuations within the colony and faster spatial genetic diversity loss. Our temporal analyses unveil the dynamics of mean cell length and fluctuations, showcasing how initial disturbances evolve as colonies are exposed to nutrients and subsequently settle. Furthermore, examining microscopic details, we find that lower resource levels yield diverse cell lengths and enhanced nematic ordering, driven by the increased prevalence of longer rod-shaped cells. Our investigation sheds light on the multifaceted relationship between resource constraints and bacterial colony dynamics, revealing insights into their spatiotemporal organization.

Dissecting cell heterogeneities in bacterial biofilms and their implications for antibiotic tolerance

Bacterial biofilms consist of large, self-formed aggregates where resident bacteria can exhibit very different physiological states and phenotypes. This heterogeneity of cell types is crucial for many structural and functional emergent properties of biofilms. Consequently, it becomes essential to understand what drives cells to differentiate and how they achieve it within the three-dimensional landscape of the biofilms. Here, we discuss recent advances in comprehending two forms of cell heterogeneity that, while recognized to coexist within biofilms, have proven challenging to distinguish. These two forms include cell heterogeneity arising as a consequence of bacteria physiologically responding to resource gradients formed across the biofilms and cell-to-cell phenotypic heterogeneity, which emerges locally within biofilm subzones among neighboring bacteria due to stochastic variations in gene expression. We describe the defining features and concepts related to both forms of cell heterogeneity and discuss their implications, with a particular focus on antibiotic tolerance.

Where bacteria and eukaryotes meet

The international workshop "Interdisciplinary life of microbes: from single cells to multicellular aggregates," following a virtual preassembly in November 2021, was held in person in Dresden, from 9 to 13 November 2022. It attracted not only prominent experts in biofilm research but also researchers from broadly neighboring disciplines, such as medicine, chemistry, and theoretical and experimental biophysics, both eukaryotic and prokaryotic. Focused brainstorming sessions were the special feature of the event and are at the heart of this commentary.

Genetic mixing and demixing on expanding spherical frontiers

Genetic fluctuation during range expansion is a key process driving evolution. When a bacterial population is expanding on a 2D surface, random fluctuations in the growth of the pioneers at the front line cause a strong demixing of genotypes. Even when there is no selective advantage, sectors of low genetic diversity are formed. Experimental studies of range expansions in surface-attached colonies of fluorescently labelled micro-organisms have contributed significantly to our understanding of fundamental evolutionary dynamics. However, experimental studies on genetic fluctuations in 3D range expansions have been sparse, despite their importance for tumour or biofilm development. We encapsulated populations of two fluorescent Escherichia coli strains in inoculation droplets (volumes [Formula: see text] nl). The confined ensemble of cells grew when embedded in a hydrogel-with nutrients-and developed 3D colonies with well-defined, sector-like regions. Using confocal laser scanning microscopy, we imaged the development of 3D colonies and the emergence of sectors. We characterized how cell concentration in the inoculation droplet controls sectors, growth rate, and the transition from branched colonies to quasi-spherical colonies. We further analysed how sectors on the surface change over time. We complement these experimental results with a modified 3D Eden growth model. The model in 3D spherical growth predicts a phase, where sectors are merging, followed by a steady increase (constant rate), and the experimentally analysed sectors were consistent with this prediction. Therefore, our results demonstrate qualitative differences between radial (2D) and spherical (3D) range expansions and their importance in gene fixation processes.

A label-free approach for relative spatial quantitation of c-di-GMP in microbial biofilms

Microbial biofilms represent an important lifestyle for bacteria and are dynamic three dimensional structures. Cyclic dimeric guanosine monophosphate (c-di-GMP) is a ubiquitous signaling molecule that is known to be tightly regulated with biofilm processes. While measurements of global levels of c-di-GMP have proven valuable towards understanding the genetic control of c-di-GMP production, there is a need for tools to observe the local changes of c-di-GMP production in biofilm processes. We have developed a label-free method for the direct detection of c-di-GMP in microbial colony biofilms using matrix-assisted laser desorption ionization mass spectrometry imaging (MALDI-MSI). We applied this method to the enteric pathogen Vibrio cholerae, the marine symbiont V. fischeri, and the opportunistic pathogen Pseudomonas aeruginosa PA14 and detected spatial and temporal changes in c-di-GMP signal that accompanied genetic alterations in factors that synthesize and degrade the compound. We further demonstrated how this method can be simultaneously applied to detect additional metabolites of interest in a single experiment.

The putative role of the epipeptide EpeX in Bacillus subtilis intra-species competition

Bacteria engage in competitive interactions with neighbours that can either be of the same or different species. Multiple mechanisms are deployed to ensure the desired outcome and one tactic commonly implemented is the production of specialised metabolites. The Gram-positive bacterium Bacillus subtilis uses specialized metabolites as part of its intra-species competition determinants to differentiate between kin and non-kin isolates. It is, however, unknown if the collection of specialized metabolites defines competitive fitness when the two isolates start as a close, interwoven community that grows into a densely packed colony biofilm. Moreover, the identity of specialized metabolites that have an active role in defining the outcome of an intra-species interaction has not been revealed. Here, we determine the competition outcomes that manifest when 21 environmental isolates of B. subtilis are individually co-incubated with the model isolate NCIB 3610 in a colony biofilm. We correlated these data with the suite of specialized metabolite biosynthesis clusters encoded by each isolate. We found that the epeXEPAB gene cluster was primarily present in isolates with a strong competitive phenotype. This cluster is responsible for producing the epipeptide EpeX. We demonstrated that EpeX is a competition determinant of B. subtilis in an otherwise isogenic context for NCBI 3610. However, when we competed the NCIB 3610 EpeX-deficient strain against our suite of environmental isolates we found that the impact of EpeX in competition is isolate-specific, as only one of the 21 isolates showed increased survival when EpeX was lacking. Taken together, we have shown that EpeX is a competition determinant used by B. subtilis that impacts intra-species interactions but only in an isolate-specific manner.

Taxonomical composition and functional analysis of biofilms sampled from a nuclear storage pool

Sampling small amounts of biofilm from harsh environments such as the biofilm present on the walls of a radioactive material storage pool offers few analytical options if taxonomic characterization and estimation of the different biomass contributions are the objectives. Although 16S/18S rRNA amplification on extracted DNA and sequencing is the most widely applied method, its reliability in terms of quantitation has been questioned as yields can be species-dependent. Here, we propose a tandem-mass spectrometry proteotyping approach consisting of acquiring peptide data and interpreting then against a generalist database without any a priori. The peptide sequence information is transformed into useful taxonomical information that allows to obtain the different biomass contributions at different taxonomical ranks. This new methodology is applied for the first time to analyze the composition of biofilms from minute quantities of material collected from a pool used to store radioactive sources in a nuclear facility. For these biofilms, we report the identification of three genera, namely Sphingomonas, Caulobacter, and Acidovorax, and their functional characterization by metaproteomics which shows that these organisms are metabolic active. Differential expression of Gene Ontology GOslim terms between the two main microorganisms highlights their metabolic specialization.