Growth dependence of conjugation explains limited plasmid invasion in biofilms: an individual-based modelling study

Biodegradation of atrazine with biochar-mediated functional bacterial biofilm: Construction, characterization and mechanisms

The abuse and residue of herbicides in the black soil area had seriously affected the soil structure, function and crop growth, posing severe threats to agricultural soil environment and public health. Given the limitation of routine microbial remediation, innovative and eco-friendly functional bacterial biofilm which could adapt under adverse conditions was developed on the biochar to investigate its enhanced bioremediation and metabolic characteristics of typical herbicide atrazine. Results revealed that the atrazine degrading strain Acinetobacter lwoffii had competitive advantage in soil indigenous microorganisms and formed dense biofilms on the biochar which was beneficial to cell viability maintenance and aggregations. Metatranscriptomics and RT-qPCR analysis demonstrated that the biochar-mediated biofilm improved the frequency of intercellular communications through quorum sensing and two-component signal regulation systems, and enhanced the atrazine biodegradation efficiency through horizontal gene transfer in co-metabolism mode, providing important scientific basis for the biological remediation of farmland soil non-point source pollution.

Flagellar interference with plasmid uptake in biofilms: a joint experimental and modeling study

Plasmid conjugation is a key facilitator of horizontal gene transfer (HGT), and plasmids encoding antibiotic resistance drive the increasing prevalence of antibiotic resistance. In natural, engineered, and clinical environments, bacteria often grow in protective biofilms. Therefore, a better understanding of plasmid transfer in biofilms is needed. Our aim was to investigate plasmid transfer in a biofilm-adapted wrinkly colony mutant of Xanthomonas retroflexus (XRw) with enhanced matrix production and reduced motility. We found that XRw biofilms had an increased uptake of the broad host-range IncP-1ϵ plasmid pKJK5 compared to the wild type (WT). Proteomics revealed fewer flagellar-associated proteins in XRw, suggesting that flagella were responsible for reducing plasmid uptake. This was confirmed by the higher plasmid uptake of non-flagellated fliM mutants of the X. retroflexus wrinkly mutant as well as the wild type. Moreover, testing several flagellar mutants of Pseudomonas putida suggested that the flagellar effect was more general. We identified seven mechanisms with the potential to explain the flagellar effect and simulated them in an individual-based model. Two mechanisms could thus be eliminated (increased distances between cells and increased lag times due to flagella). Another mechanism identified as viable in the modeling was eliminated by further experiments. The possibility of steric hindrance of pilus movement and binding by flagella, reducing the frequency of contact and thus plasmid uptake, proved viable, and the three other viable mechanisms had a reduced probability of plasmid transfer in common. Our findings highlight the important yet complex effects of flagella during bacterial conjugation in biofilms.IMPORTANCEBiofilms are the dominant form of microbial life and bacteria living in biofilms are markedly different from their planktonic counterparts, yet the impact of the biofilm lifestyle on horizontal gene transfer (HGT) is still poorly understood. Horizontal gene transfer by conjugative plasmids is a major driver in bacterial evolution and adaptation, as exemplified by the troubling spread of antibiotic resistance. To either limit or promote plasmid prevalence and dissemination, we need a better understanding of plasmid transfer between bacterial cells, especially in biofilms. Here, we identified a new factor impacting the transfer of plasmids, flagella, which are required for many types of bacterial motility. We show that their absence or altered activity can lead to enhanced plasmid uptake in two bacterial species, Xanthomonas retroflexus and Pseudomonas putida. Moreover, we demonstrate the utility of mathematical modeling to eliminate hypothetical mechanisms.

In Situ and Individual-Based Analysis of the Influence of Polystyrene Microplastics on Escherichia coli Conjugative Gene Transfer at the Single-Cell Level

The impact of microplastic particles of micro- and nanometer sizes on microbial horizontal gene transfer (HGT) remains a controversial topic. Existing studies rely on traditional approaches, which analyze population behavior, leading to conflicting conclusions and a limited understanding. The present study addressed these limitations by employing a novel microfluidic chamber system for in situ visualization and precise quantification of the effects of different concentrations of polystyrene (PS) microbeads on microbial HGT at the single-cell level. The statistical analysis indicated no significant difference in the division times of both the donor and recipient bacteria across different PS microbead concentrations. However, as the concentration of PS microbeads increased from 0 to 2000 mg L-1, the average conjugation frequency of Escherichia coli decreased from 0.028 ± 0.015 to 0.004 ± 0.003. Our observations from the microfluidic experiments revealed that 500 nm PS microbeads created a barrier effect on bacterial conjugative transfer. The presence of microbeads resulted in reduced contact and interaction between the donor and recipient strains, thereby causing a decrease in the conjugation transfer frequency. These findings were validated by an individual-based modeling framework parameterized by the data from the individual-level microfluidic experiments. Overall, this study offers a fresh perspective and strategy for investigating the risks associated with the dissemination of antibiotic resistance genes related to microplastics.

Agent-Based Modeling of Microbial Communities

Microbial communities are complex living systems that populate the planet with diverse functions and are increasingly harnessed for practical human needs. To deepen the fundamental understanding of their organization and functioning as well as to facilitate their engineering for applications, mathematical modeling has played an increasingly important role. Agent-based models represent a class of powerful quantitative frameworks for investigating microbial communities because of their individualistic nature in describing cells, mechanistic characterization of molecular and cellular processes, and intrinsic ability to produce emergent system properties. This review presents a comprehensive overview of recent advances in agent-based modeling of microbial communities. It surveys the state-of-the-art algorithms employed to simulate intracellular biomolecular events, single-cell behaviors, intercellular interactions, and interactions between cells and their environments that collectively serve as the driving forces of community behaviors. It also highlights three lines of applications of agent-based modeling, namely, the elucidation of microbial range expansion and colony ecology, the design of synthetic gene circuits and microbial populations for desired behaviors, and the characterization of biofilm formation and dispersal. The review concludes with a discussion of existing challenges, including the computational cost of the modeling, and potential mitigation strategies.

Droplet size and surface hydrophobicity enhance bacterial plasmid transfer rates in microscopic surface wetness

Conjugal plasmids constitute a major engine for horizontal gene transfer in bacteria, and are key drivers of the spread of antibiotic resistance, virulence, and metabolic functions. Bacteria in terrestrial habitats often inhabit surfaces that are not constantly water-saturated, where microscopic surface wetness (MSW), comprised of thin liquid films and microdroplets, permanently or intermittently occurs. How physical properties of microdroplets, and of the surfaces they reside on, affect plasmid transfer rates is not well understood. Here, building on microscopy-based microdroplet experiments, we examined the relation between droplet properties (size and spread) and plasmid transfer rates at single-cell and individual droplet resolution, using Pseudomonas putida as a model species. We show that transfer rates increase with droplet size, due to higher densities of cells on the surface in larger droplets, resulting from lower ratio between the area of the liquid-solid interface and droplet volumes. We further show that surface hydrophobicity promotes transfer rates via the same mechanism. Our results provide new insights into how physical properties of surfaces and MSW affect plasmid transfer rates, and more generally, microbial interactions mediated by cell-to-cell contact, with important implications for our understanding of the ecology and evolution of bacteria in unsaturated environments.

Models for Gut-Mediated Horizontal Gene Transfer by Bacterial Plasmid Conjugation

The emergence of new antimicrobial resistant and virulent bacterial strains may pose a threat to human and animal health. Bacterial plasmid conjugation is a significant contributor to rapid microbial evolutions that results in the emergence and spread of antimicrobial resistance (AR). The gut of animals is believed to be a potent reservoir for the spread of AR and virulence genes through the horizontal exchange of mobile genetic elements such as plasmids. The study of the plasmid transfer process in the complex gut environment is limited due to the confounding factors that affect colonization, persistence, and plasmid conjugation. Furthermore, study of plasmid transfer in the gut of humans is limited to observational studies, leading to the need to identify alternate models that provide insight into the factors regulating conjugation in the gut. This review discusses key studies on the current models for in silico, in vitro, and in vivo modeling of bacterial conjugation, and their ability to reflect the gut of animals. We particularly emphasize the use of computational and in vitro models that may approximate aspects of the gut, as well as animal models that represent in vivo conditions to a greater extent. Directions on future research studies in the field are provided.

Simultaneous estimation of parameters governing the vertical and horizontal transfer of antibiotic resistance genes

The accelerated spread of antibiotic resistance genes (ARG) in the environment occurs mainly through plasmid transfer facilitated via bacterial conjugation. To predict and efficiently counteract the problems associated with ARG transmission, it is important to estimate conjugation rates under different experimental conditions. The classical models typically used to estimate parameters for mating experiments, while pragmatic in calculating growth and plasmid transfer, often ignore processes such as the reduction in growth due to plasmid bearing costs and are non-inclusive of environmental influences like temperature effects. Here, we present a process-based numerical model taking into account the fitness cost associated with plasmid carriage and temperature dependencies in vertical and horizontal gene transfer processes. Observations from liquid culture conjugation experiments using Escherichia coli and the plasmid pB10 were used to validate our proposed model. We present a comparison between the parameters estimated using the existing and the proposed model. Uncertainties in the estimated parameters were quantified using classical and advanced Bayesian methods. For our mating experiments, we found that at temperatures between 20 and 37 °C, the plasmid bearing costs reduced the growth rates by > 35%. The temperature dependency model of conjugation showed a good fit (mean absolute percentage error < 10%) independent of the bacteria and the plasmid under study. The proposed model simultaneously estimates growth and plasmid transfer rate constants for all three strains (donor, recipient, and transconjugant). Simultaneous estimation of growth and conjugation parameters is particularly useful to estimate the spread of ARG when one of the mating partners inhibits the growth of the other, which is common in multi-species mating or when the incurred plasmid costs are situation dependent (e.g., increased plasmid cost in a mating environment) as observed in this study.

Mathematical Models of Plasmid Population Dynamics

With plasmid-mediated antibiotic resistance thriving and threatening to become a serious public health problem, it is paramount to increase our understanding of the forces that enable the spread and maintenance of drug resistance genes encoded in mobile genetic elements. The relevance of plasmids as vehicles for the dissemination of antibiotic resistance genes, in addition to the extensive use of plasmid-derived vectors for biotechnological and industrial purposes, has promoted the in-depth study of the molecular mechanisms controlling multiple aspects of a plasmids' life cycle. This body of experimental work has been paralleled by the development of a wealth of mathematical models aimed at understanding the interplay between transmission, replication, and segregation, as well as their consequences in the ecological and evolutionary dynamics of plasmid-bearing bacterial populations. In this review, we discuss theoretical models of plasmid dynamics that span from the molecular mechanisms of plasmid partition and copy-number control occurring at a cellular level, to their consequences in the population dynamics of complex microbial communities. We conclude by discussing future directions for this exciting research topic.

A Review on Occurrence and Spread of Antibiotic Resistance in Wastewaters and in Wastewater Treatment Plants: Mechanisms and Perspectives

This paper reviews current knowledge on sources, spread and removal mechanisms of antibiotic resistance genes (ARGs) in microbial communities of wastewaters, treatment plants and downstream recipients. Antibiotic is the most important tool to cure bacterial infections in humans and animals. The over- and misuse of antibiotics have played a major role in the development, spread, and prevalence of antibiotic resistance (AR) in the microbiomes of humans and animals, and microbial ecosystems worldwide. AR can be transferred and spread amongst bacteria via intra- and interspecies horizontal gene transfer (HGT). Wastewater treatment plants (WWTPs) receive wastewater containing an enormous variety of pollutants, including antibiotics, and chemicals from different sources. They contain large and diverse communities of microorganisms and provide a favorable environment for the spread and reproduction of AR. Existing WWTPs are not designed to remove micropollutants, antibiotic resistant bacteria (ARB) and ARGs, which therefore remain present in the effluent. Studies have shown that raw and treated wastewaters carry a higher amount of ARB in comparison to surface water, and such reports have led to further studies on more advanced treatment processes. This review summarizes what is known about AR removal efficiencies of different wastewater treatment methods, and it shows the variations among different methods. Results vary, but the trend is that conventional activated sludge treatment, with aerobic and/or anaerobic reactors alone or in series, followed by advanced post treatment methods like UV, ozonation, and oxidation removes considerably more ARGs and ARB than activated sludge treatment alone. In addition to AR levels in treated wastewater, it examines AR levels in biosolids, settled by-product from wastewater treatment, and discusses AR removal efficiency of different biosolids treatment procedures. Finally, it puts forward key-points and suggestions for dealing with and preventing further increase of AR in WWTPs and other aquatic environments, together with a discussion on the use of mathematical models to quantify and simulate the spread of ARGs in WWTPs. Mathematical models already play a role in the analysis and development of WWTPs, but they do not consider AR and challenges remain before models can be used to reliably study the dynamics and reduction of AR in such systems.

The Evolution of Plasmid Transfer Rate in Bacteria and Its Effect on Plasmid Persistence

AbstractPlasmids are extrachromosomal segments of DNA that can transfer genes between bacterial cells. Many plasmid genes benefit bacteria but cause harm to human health by granting antibiotic resistance to pathogens. Transfer rate is a key parameter for predicting plasmid dynamics, but observed rates are highly variable, and the effects of selective forces on their evolution are unclear. We apply evolutionary analysis to plasmid conjugation models to investigate selective pressures affecting plasmid transfer rate, emphasizing host versus plasmid control, the costs of plasmid transfer, and the role of recipient cells. Our analyses show that plasmid-determined transfer rates can be predicted with three parameters (host growth rate, plasmid loss rate, and the cost of plasmid transfer on growth) under some conditions. We also show that low-frequency genetic variation in transfer rate can accumulate, facilitating rapid adaptation to changing conditions. Furthermore, reduced transfer rates due to host control have limited effects on plasmid prevalence until low enough to prevent plasmid persistence. These results provide a framework to predict plasmid transfer rate evolution in different environments and demonstrate the limited impact of host mechanisms to control the costs incurred when plasmids are present.

The ecology of plasmid-coded antibiotic resistance: a basic framework for experimental research and modeling

Many antibiotic resistance genes are associated with plasmids. The ecological success of these mobile genetic elements within microbial communities depends on varying mechanisms to secure their own propagation, not only on environmental selection. Among the most important are the cost of plasmids and their ability to be transferred to new hosts through mechanisms such as conjugation. These are regulated by dynamic control systems of the conjugation machinery and genetic adaptations that plasmid-host pairs can acquire in coevolution. However, in complex communities, these processes and mechanisms are subject to a variety of interactions with other bacterial species and other plasmid types. This article summarizes basic plasmid properties and ecological principles particularly important for understanding the persistence of plasmid-coded antibiotic resistance in aquatic environments. Through selected examples, it further introduces to the features of different types of simulation models such as systems of ordinary differential equations and individual-based models, which are considered to be important tools to understand these complex systems. This ecological perspective aims to improve the way we study and understand the dynamics, diversity and persistence of plasmids and associated antibiotic resistance genes.

From molecules to multispecies ecosystems: the roles of structure in bacterial biofilms

Biofilms are communities of sessile microbes that are bound to each other by a matrix made of biopolymers and proteins. Spatial structure is present in biofilms on many lengthscales. These range from the nanometer scale of molecular motifs to the hundred-micron scale of multicellular aggregates. Spatial structure is a physical property that impacts the biology of biofilms in many ways. The molecular structure of matrix components controls their interaction with each other (thereby impacting biofilm mechanics) and with diffusing molecules such as antibiotics and immune factors (thereby impacting antibiotic tolerance and evasion of the immune system). The size and structure of multicellular aggregates, combined with microbial consumption of growth substrate, give rise to differentiated microenvironments with different patterns of metabolism and gene expression. Spatial association of more than one species can benefit one or both species, while distances between species can both determine and result from the transport of diffusible factors between species. Thus, a widespread theme in the biological importance of spatial structure in biofilms is the effect of structure on transport. We survey what is known about this and other effects of spatial structure in biofilms, from molecules up to multispecies ecosystems. We conclude with an overview of what experimental approaches have been developed to control spatial structure in biofilms and how these and other experiments can be complemented with computational work.

Recent Advances in Anti-virulence Therapeutic Strategies With a Focus on Dismantling Bacterial Membrane Microdomains, Toxin Neutralization, Quorum-Sensing Interference and Biofilm Inhibition

Antimicrobial resistance constitutes one of the major challenges facing humanity in the Twenty-First century. The spread of resistant pathogens has been such that the possibility of returning to a pre-antibiotic era is real. In this scenario, innovative therapeutic strategies must be employed to restrict resistance. Among the innovative proposed strategies, anti-virulence therapy has been envisioned as a promising alternative for effective control of the emergence and spread of resistant pathogens. This review presents some of the anti-virulence strategies that are currently being developed, it will cover strategies focused on quench pathogen quorum sensing (QS) systems, disassemble of bacterial functional membrane microdomains (FMMs), disruption of biofilm formation and bacterial toxin neutralization.

Cell-to-cell bacterial interactions promoted by drier conditions on soil surfaces

Bacterial cell-to-cell interactions are in the core of evolutionary and ecological processes in soil and other environments. Under most conditions, natural soils are unsaturated where the fragmented aqueous habitats and thin liquid films confine bacterial cells within small volumes and close proximity for prolonged periods. We report effects of a range of hydration conditions on bacterial cell-level interactions that are marked by plasmid transfer between donor and recipient cells within populations of the soil bacterium Pseudomonas putida Using hydration-controlled sand microcosms, we demonstrate that the frequency of cell-to-cell contacts under prescribed hydration increases with lowering water potential values (i.e., under drier conditions where the aqueous phase shrinks and fragments). These observations were supported using a mechanistic individual-based model for linking macroscopic soil water potential to microscopic distribution of liquid phase and explicit bacterial cell interactions in a simplified porous medium. Model results are in good agreement with observations and inspire confidence in the underlying mechanisms. The study highlights important physical factors that control short-range bacterial cell interactions in soil and on surfaces, specifically, the central role of the aqueous phase in mediating bacterial interactions and conditions that promote genetic information transfer in support of soil microbial diversity.

Agent-Based Modeling in Molecular Systems Biology

Molecular systems orchestrating the biology of the cell typically involve a complex web of interactions among various components and span a vast range of spatial and temporal scales. Computational methods have advanced our understanding of the behavior of molecular systems by enabling us to test assumptions and hypotheses, explore the effect of different parameters on the outcome, and eventually guide experiments. While several different mathematical and computational methods are developed to study molecular systems at different spatiotemporal scales, there is still a need for methods that bridge the gap between spatially-detailed and computationally-efficient approaches. In this review, we summarize the capabilities of agent-based modeling (ABM) as an emerging molecular systems biology technique that provides researchers with a new tool in exploring the dynamics of molecular systems/pathways in health and disease.

A Bayesian approach to modelling the impact of hydrodynamic shear stress on biofilm deformation

We investigate the feasibility of using a surrogate-based method to emulate the deformation and detachment behaviour of a biofilm in response to hydrodynamic shear stress. The influence of shear force, growth rate and viscoelastic parameters on the patterns of growth, structure and resulting shape of microbial biofilms was examined. We develop a statistical modelling approach to this problem, using combination of Bayesian Poisson regression and dynamic linear models for the emulation. We observe that the hydrodynamic shear force affects biofilm deformation in line with some literature. Sensitivity results also showed that the expected number of shear events, shear flow, yield coefficient for heterotrophic bacteria and extracellular polymeric substance (EPS) stiffness per unit EPS mass are the four principal mechanisms governing the bacteria detachment in this study. The sensitivity of the model parameters is temporally dynamic, emphasising the significance of conducting the sensitivity analysis across multiple time points. The surrogate models are shown to perform well, and produced ≈ 480 fold increase in computational efficiency. We conclude that a surrogate-based approach is effective, and resulting biofilm structure is determined primarily by a balance between bacteria growth, viscoelastic parameters and applied shear stress.

Towards an integrated understanding of how micro scale processes shape groundwater ecosystem functions

Micro scale processes are expected to have a fundamental role in shaping groundwater ecosystems and yet they remain poorly understood and under-researched. In part, this is due to the fact that sampling is rarely carried out at the scale at which microorganisms, and their grazers and predators, function and thus we lack essential information. While set within a larger scale framework in terms of geochemical features, supply with energy and nutrients, and exchange intensity and dynamics, the micro scale adds variability, by providing heterogeneous zones at the micro scale which enable a wider range of redox reactions. Here we outline how understanding micro scale processes better may lead to improved appreciation of the range of ecosystems functions taking place at all scales. Such processes are relied upon in bioremediation and we demonstrate that ecosystem modelling as well as engineering measures have to take into account, and use, understanding at the micro scale. We discuss the importance of integrating faunal processes and computational appraisals in research, in order to continue to secure sustainable water resources from groundwater.

Conjugative plasmids enable the maintenance of low cost non-transmissible plasmids

Some plasmids can be transferred by conjugation to other bacterial hosts. But almost half of the plasmids are non-transmissible. These plasmid types can only be transmitted to the daughter cells of their host after bacterial fission. Previous studies suggest that non-transmissible plasmids become extinct in the absence of selection of their encoded traits, as plasmid-free bacteria are more competitive. Here, we aim to identify mechanisms that enable non-transmissible plasmids to persist, even if they are not beneficial. For this purpose, an individual-based model for plasmid population dynamics was set up and carefully tested for structural consistency and plausibility. Our results demonstrate that non-transmissible plasmids can be stably maintained in a population, even if they impose a substantial burden on their host cells growth. A prerequisite is the co-occurrence of an incompatible and costly conjugative plasmid type, which indirectly facilitates the preservation of the non-transmissible type. We suggest that this constellation might be considered as a potential mechanism maintaining plasmids and associated antibiotic resistances. It should be investigated in upcoming laboratory experiments.

Predictive Modeling of a Batch Filter Mating Process

Quantitative characterizations of horizontal gene transfer are needed to accurately describe gene transfer processes in natural and engineered systems. A number of approaches to the quantitative description of plasmid conjugation have appeared in the literature. In this study, we seek to extend that work, motivated by the question of whether a mathematical model can accurately predict growth and conjugation dynamics in a batch process. We used flow cytometry to make time-point observations of a filter-associated mating between two E. coli strains, and fit ordinary differential equation models to the data. A model comparison analysis identified the model formulation that is best supported by the data. Identifiability analysis revealed that the parameters were estimated with acceptable accuracy. The predictive power of the model was assessed by comparison with test data that demanded extrapolation from the training experiments. This study represents the first attempt to assess the quality of model predictions for plasmid conjugation. Our successful application of this approach lays a foundation for predictive modeling that can be used both in the study of natural plasmid transmission and in model-based design of engineering approaches that employ conjugation, such as plasmid-mediated bioaugmentation.

Growth phase-dependent control of R27 conjugation is mediated by the interplay between the plasmid-encoded regulatory circuit TrhR/TrhY-HtdA and the cAMP regulon

Plasmids of the incompatibility group HI1 (IncHI1) have been isolated from several Gram-negative pathogens and are associated with the spread of multidrug resistance. Their conjugation is tightly regulated and it is inhibited at temperatures higher than 30°C, indicating that conjugation occurs outside warm-blooded hosts. Using R27, the prototype of IncHI1 plasmids, we report that plasmid transfer efficiency in E. coli strongly depends on the physiological state of the donor cells. Conjugation frequency is high when cells are actively growing, dropping sharply when cells enter the stationary phase of growth. Accordingly, our transcriptomic assays show significant downregulation of numerous R27 genes during the stationary phase, including several tra (transfer) genes. Growth phase-dependent regulation of tra genes transcription is independent of H-NS, a silencer of horizontal gene transfer, and ppGpp and RpoS, regulators of the stationary phase, but highly dependent on the plasmid-encoded regulatory circuit TrhR/TrhY-HtdA. The metabolic sensor cAMP, whose synthesis is chromosomally encoded, is also involved in the growth phase regulation of R27 conjugation by modulating htdA expression. Our data suggest that the involvement of regulators encoded by both chromosome and plasmid are required for efficient physiological control of IncHI1 plasmid conjugation.

Mycelia as a focal point for horizontal gene transfer among soil bacteria

Horizontal gene transfer (HGT) is a main mechanism of bacterial evolution endowing bacteria with new genetic traits. The transfer of mobile genetic elements such as plasmids (conjugation) requires the close proximity of cells. HGT between genetically distinct bacteria largely depends on cell movement in water films, which are typically discontinuous in natural systems like soil. Using laboratory microcosms, a bacterial reporter system and flow cytometry, we here investigated if and to which degree mycelial networks facilitate contact of and HGT between spatially separated bacteria. Our study shows that the network structures of mycelia promote bacterial HGT by providing continuous liquid films in which bacterial migration and contacts are favoured. This finding was confirmed by individual-based simulations, revealing that the tendency of migrating bacteria to concentrate in the liquid film around hyphae is a key factor for improved HGT along mycelial networks. Given their ubiquity, we propose that hyphae can act as focal point for HGT and genetic adaptation in soil.

Plasmid transfer in biofilms: a perspective on limitations and opportunities

Biofilms dominate microbial life, and their importance for human health and the environment can no longer be dismissed. Nevertheless many of the processes governing this form of microbial growth are still poorly understood. This includes the horizontal exchange of genetic information, which is a major driver in bacterial evolution and rapid adaptation, exemplified by the alarming spread of multi-drug resistance among pathogens mediated by plasmids. Biofilms are often considered hot spot for horizontal gene transfer, yet several studies have shown that plasmid transfer is limited to the outer layers. On the basis of results from decades of research we analyse this paradox and discuss the mechanisms by which biofilm growth can promote the initial transfer of some plasmids, but also limit further plasmid invasion into the population or community. If we want to adequately promote or combat horizontal gene spread in biofilms, we need to gain better insight into the physicochemical and biological mechanisms that control this process.

Advancing microbial sciences by individual-based modelling

Remarkable technological advances have revealed ever more properties and behaviours of individual microorganisms, but the novel data generated by these techniques have not yet been fully exploited. In this Opinion article, we explain how individual-based models (IBMs) can be constructed based on the findings of such techniques and how they help to explore competitive and cooperative microbial interactions. Furthermore, we describe how IBMs have provided insights into self-organized spatial patterns from biofilms to the oceans of the world, phage-CRISPR dynamics and other emergent phenomena. Finally, we discuss how combining individual-based observations with IBMs can advance our understanding at both the individual and population levels, leading to the new approach of microbial individual-based ecology (μIBE).

Conjugation in Gram-Positive Bacteria

Conjugative transfer is the most important means of spreading antibiotic resistance and virulence factors among bacteria. The key vehicles of this horizontal gene transfer are a group of mobile genetic elements, termed conjugative plasmids. Conjugative plasmids contain as minimum instrumentation an origin of transfer (oriT), DNA-processing factors (a relaxase and accessory proteins), as well as proteins that constitute the trans-envelope transport channel, the so-called mating pair formation (Mpf) proteins. All these protein factors are encoded by one or more transfer (tra) operons that together form the DNA transport machinery, the Gram-positive type IV secretion system. However, multicellular Gram-positive bacteria belonging to the streptomycetes appear to have evolved another mechanism for conjugative plasmid spread reminiscent of the machinery involved in bacterial cell division and sporulation, which transports double-stranded DNA from donor to recipient cells. Here, we focus on the protein key players involved in the plasmid spread through the two different modes and present a new secondary structure homology-based classification system for type IV secretion protein families. Moreover, we discuss the relevance of conjugative plasmid transfer in the environment and summarize novel techniques to visualize and quantify conjugative transfer in situ.

Mathematical modelling of antimicrobial resistance in agricultural waste highlights importance of gene transfer rate

Antimicrobial resistance is of global concern. Most antimicrobial use is in agriculture; manures and slurry are especially important because they contain a mix of bacteria, including potential pathogens, antimicrobial resistance genes and antimicrobials. In many countries, manures and slurry are stored, especially over winter, before spreading onto fields as organic fertilizer. Thus, these are a potential location for gene exchange and selection for resistance. We develop and analyse a mathematical model to quantify the spread of antimicrobial resistance in stored agricultural waste. We use parameters from a slurry tank on a UK dairy farm as an exemplar. We show that the spread of resistance depends in a subtle way on the rates of gene transfer and antibiotic inflow. If the gene transfer rate is high, then its reduction controls resistance, while cutting antibiotic inflow has little impact. If the gene transfer rate is low, then reducing antibiotic inflow controls resistance. Reducing length of storage can also control spread of resistance. Bacterial growth rate, fitness costs of carrying antimicrobial resistance and proportion of resistant bacteria in animal faeces have little impact on spread of resistance. Therefore, effective treatment strategies depend critically on knowledge of gene transfer rates.

Forging patterns and making waves from biology to geology: a commentary on Turing (1952) 'The chemical basis of morphogenesis'

Alan Turing was neither a biologist nor a chemist, and yet the paper he published in 1952, 'The chemical basis of morphogenesis', on the spontaneous formation of patterns in systems undergoing reaction and diffusion of their ingredients has had a substantial impact on both fields, as well as in other areas as disparate as geomorphology and criminology. Motivated by the question of how a spherical embryo becomes a decidedly non-spherical organism such as a human being, Turing devised a mathematical model that explained how random fluctuations can drive the emergence of pattern and structure from initial uniformity. The spontaneous appearance of pattern and form in a system far away from its equilibrium state occurs in many types of natural process, and in some artificial ones too. It is often driven by very general mechanisms, of which Turing's model supplies one of the most versatile. For that reason, these patterns show striking similarities in systems that seem superficially to share nothing in common, such as the stripes of sand ripples and of pigmentation on a zebra skin. New examples of 'Turing patterns' in biology and beyond are still being discovered today. This commentary was written to celebrate the 350th anniversary of the journal Philosophical Transactions of the Royal Society.

Patterning of individual Staphylococcus aureus bacteria onto photogenerated polymeric surface structures

This manuscript describes the fabrication of bacterial surface arrays by using photolithographic techniques having in addition some particularly interesting features.

Novel assay to measure the plasmid mobilizing potential of mixed microbial communities

Mobilizable plasmids lack necessary genes for complete conjugation and are therefore non-self-transmissible. Instead, they rely on the conjugation system of conjugal plasmids to be horizontally transferred to new recipients. While community permissiveness, the fraction of a mixed microbial community that can receive self-transmissible conjugal plasmids, has been studied, the intrinsic ability of a community to mobilize plasmids that lack conjugation systems is unexplored. Here, we present a novel framework and experimental method to estimate the mobilization potential of mixed communities. We compare the transfer frequency of a mobilizable plasmid to that of a mobilizing and conjugal plasmid measured for a model strain and for the assayed community. With Pseudomonas putida carrying the gfp-tagged mobilizable IncQ plasmid RSF1010 as donor strain, we conducted solid surface mating experiments with either a P. putida strain carrying the mobilizing IncP-1α plasmid RP4 or a model bacterial community that was extracted from the inner walls of a domestic shower conduit. Additionally, we estimated the permissiveness of the same community for RP4 using P. putida as donor strain. The permissiveness of the model community for RP4 [at 1.16 × 10^-4 transconjugants per recipient (T/R)] was similar to that previously measured for soil microbial communities. RSF1010 was mobilized by the model community at a frequency of 1.16 × 10^-5 T/R, only one order of magnitude lower than its permissiveness to RP4. This mobilization frequency is unexpectedly high considering that (i) mobilization requires the presence of mobilizing conjugal plasmids within the permissive fraction of the recipients; (ii) in pure culture experiments with P. putida retromobilization of RSF1010 through RP4 only took place in approximately half of the donors receiving the conjugal plasmid in the first step. Further work is needed to establish how plasmid mobilization potential varies within and across microbial communities. This method has the potential to provide such insights; in addition it allows for the direct isolation of in situ mobilizing plasmids together with their endogenous hosts.

Repair rather than segregation of damage is the optimal unicellular aging strategy

BACKGROUND

How aging, being unfavourable for the individual, can evolve is one of the fundamental problems of biology. Evidence for aging in unicellular organisms is far from conclusive. Some studies found aging even in symmetrically dividing unicellular species; others did not find aging in the same, or in different, unicellular species, or only under stress. Mathematical models suggested that segregation of non-genetic damage, as an aging strategy, would increase fitness. However, these models failed to consider repair as an alternative strategy or did not properly account for the benefits of repair. We used a new and improved individual-based model to examine rigorously the effect of a range of aging strategies on fitness in various environments.

RESULTS

Repair of damage emerges as the best strategy despite its fitness costs, since it immediately increases growth rate. There is an optimal investment in repair that outperforms damage segregation in well-mixed, lasting and benign environments over a wide range of parameter values. Damage segregation becomes beneficial, and only in combination with repair, when three factors are combined: (i) the rate of damage accumulation is high, (ii) damage is toxic and (iii) efficiency of repair is low. In contrast to previous models, our model predicts that unicellular organisms should have active mechanisms to repair damage rather than age by segregating damage. Indeed, as predicted, all organisms have evolved active mechanisms of repair whilst aging in unicellular organisms is absent or minimal under benign conditions, apart from microorganisms with a different ecology, inhabiting short-lived environments strongly favouring early reproduction rather than longevity.

CONCLUSIONS

Aging confers no fitness advantage for unicellular organisms in lasting environments under benign conditions, since repair of non-genetic damage is better than damage segregation.

Towards predictive models of the human gut microbiome

The intestinal microbiota is an ecosystem susceptible to external perturbations such as dietary changes and antibiotic therapies. Mathematical models of microbial communities could be of great value in the rational design of microbiota-tailoring diets and therapies. Here, we discuss how advances in another field, engineering of microbial communities for wastewater treatment bioreactors, could inspire development of mechanistic mathematical models of the gut microbiota. We review the state of the art in bioreactor modeling and current efforts in modeling the intestinal microbiota. Mathematical modeling could benefit greatly from the deluge of data emerging from metagenomic studies, but data-driven approaches such as network inference that aim to predict microbiome dynamics without explicit mechanistic knowledge seem better suited to model these data. Finally, we discuss how the integration of microbiome shotgun sequencing and metabolic modeling approaches such as flux balance analysis may fulfill the promise of a mechanistic model.

Metabolic and demographic feedbacks shape the emergent spatial structure and function of microbial communities

Microbes are predominantly found in surface-attached and spatially structured polymicrobial communities. Within these communities, microbial cells excrete a wide range of metabolites, setting the stage for interspecific metabolic interactions. The links, however, between metabolic and ecological interactions (functional relationships), and species spatial organization (structural relationships) are still poorly understood. Here, we use an individual-based modelling framework to simulate the growth of a two-species surface-attached community where food (resource) is traded for detoxification (service) and investigate how metabolic constraints of individual species shape the emergent structural and functional relationships of the community. We show that strong metabolic interdependence drives the emergence of mutualism, robust interspecific mixing, and increased community productivity. Specifically, we observed a striking and highly stable emergent lineage branching pattern, generating a persistent lineage mixing that was absent when the metabolic exchange was removed. These emergent community properties are driven by demographic feedbacks, such that aid from neighbouring cells directly enhances focal cell growth, which in turn feeds back to neighbour fecundity. In contrast, weak metabolic interdependence drives conflict (exploitation or competition), and in turn greater interspecific segregation. Together, these results support the idea that species structural and functional relationships represent the net balance of metabolic interdependencies.

Understanding the structure and functioning of polymicrobial communities is a major challenge in biology, as witnessed by the dramatic yet mysterious roles played by the human microbiome in human health. Microbial multispecies communities often show complex spatial structures and patterns of metabolic exchange, yet our understanding of how species spatial and ecological relationships emerge from the metabolic rules of species interactions is still limited. What mechanisms underlie multispecies community self-organization? In this study, we simulate the growth of a minimal—two species— community and show how the emergent properties of community spatial structure and function depend on the nature of metabolic interactions between the two species. We found that strong mutual need for help (strong metabolic interdependence) favours the emergence of mutualism, increased productivity and lineage mixing via striking and highly stable branching patterns. In contrast, when the mutual need for help is low, conflict dominates and the two species tend to segregate. Finally, we show how the emergent species mixing follows from a positive feedback of providing aid to neighbouring helpers.

Speed of invasion of an expanding population by a horizontally transmitted trait

Range expansions are a ubiquitous phenomenon, leading to the spatial spread of genetic, ecological, and cultural traits. While some of these traits are advantageous (and hence selected), other, nonselected traits can also spread by hitchhiking on the wave of population expansion. This requires us to understand how the spread of a hitchhiking trait is coupled to the wave of advance of its host population. Here, we use a system of coupled Fisher-Kolmogorov-Petrovsky-Piskunov (F-KPP) equations to describe the spread of a horizontally transmitted hitchhiking trait within a population as it expands. We extend F-KPP wave theory to the system of coupled equations to predict how the hitchhiking trait spreads as a wave within the expanding population. We show that the speed of this trait wave is controlled by an intricate coupling between the tip of the population and trait waves. Our analysis yields a new speed selection mechanism for coupled waves of advance and reveals the existence of previously unexpected speed transitions.

Novel multiscale modeling tool applied to Pseudomonas aeruginosa biofilm formation

Multiscale modeling is used to represent biological systems with increasing frequency and success. Multiscale models are often hybrids of different modeling frameworks and programming languages. We present the MATLAB-NetLogo extension (MatNet) as a novel tool for multiscale modeling. We demonstrate the utility of the tool with a multiscale model of Pseudomonas aeruginosa biofilm formation that incorporates both an agent-based model (ABM) and constraint-based metabolic modeling. The hybrid model correctly recapitulates oxygen-limited biofilm metabolic activity and predicts increased growth rate via anaerobic respiration with the addition of nitrate to the growth media. In addition, a genome-wide survey of metabolic mutants and biofilm formation exemplifies the powerful analyses that are enabled by this computational modeling tool.

Multicellular computing using conjugation for wiring

Recent efforts in synthetic biology have focussed on the implementation of logical functions within living cells. One aim is to facilitate both internal "re-programming" and external control of cells, with potential applications in a wide range of domains. However, fundamental limitations on the degree to which single cells may be re-engineered have led to a growth of interest in multicellular systems, in which a "computation" is distributed over a number of different cell types, in a manner analogous to modern computer networks. Within this model, individual cell type perform specific sub-tasks, the results of which are then communicated to other cell types for further processing. The manner in which outputs are communicated is therefore of great significance to the overall success of such a scheme. Previous experiments in distributed cellular computation have used global communication schemes, such as quorum sensing (QS), to implement the "wiring" between cell types. While useful, this method lacks specificity, and limits the amount of information that may be transferred at any one time. We propose an alternative scheme, based on specific cell-cell conjugation. This mechanism allows for the direct transfer of genetic information between bacteria, via circular DNA strands known as plasmids. We design a multi-cellular population that is able to compute, in a distributed fashion, a Boolean XOR function. Through this, we describe a general scheme for distributed logic that works by mixing different strains in a single population; this constitutes an important advantage of our novel approach. Importantly, the amount of genetic information exchanged through conjugation is significantly higher than the amount possible through QS-based communication. We provide full computational modelling and simulation results, using deterministic, stochastic and spatially-explicit methods. These simulations explore the behaviour of one possible conjugation-wired cellular computing system under different conditions, and provide baseline information for future laboratory implementations.

The interconnection between biofilm formation and horizontal gene transfer

Recent research has revealed that horizontal gene transfer and biofilm formation are connected processes. Although published research investigating this interconnectedness is still limited, we will review this subject in order to highlight the potential of these observations because of their believed importance in the understanding of the adaptation and subsequent evolution of social traits in bacteria. Here, we discuss current evidence for such interconnectedness centred on plasmids. Horizontal transfer rates are typically higher in biofilm communities compared with those in planktonic states. Biofilms, furthermore, promote plasmid stability and may enhance the host range of mobile genetic elements that are transferred horizontally. Plasmids, on the other hand, are very well suited to promote the evolution of social traits such as biofilm formation. This, essentially, transpires because plasmids are independent replicons that enhance their own success by promoting inter-bacterial interactions. They typically also carry genes that heighten their hosts' direct fitness. Furthermore, current research shows that the so-called mafia traits encoded on mobile genetic elements can enforce bacteria to maintain stable social interactions. It also indicates that horizontal gene transfer ultimately enhances the relatedness of bacteria carrying the mobile genetic elements of the same origin. The perspective of this review extends to an overall interconnectedness between horizontal gene transfer, mobile genetic elements and social evolution of bacteria.

iDynoMiCS: next-generation individual-based modelling of biofilms

Individual-based modelling of biofilms accounts for the fact that individual organisms of the same species may well be in a different physiological state as a result of environmental gradients, lag times in responding to change, or noise in gene expression, which we have become increasingly aware of with the advent of single-cell microbiology. But progress in developing and using individual-based modelling has been hampered by different groups writing their own code and the lack of an available standard model. We therefore set out to merge most features of previous models and incorporate various improvements in order to provide a common basis for further developments. Four improvements stand out: the biofilm pressure field allows for shrinking or consolidating biofilms; the continuous-in-time extracellular polymeric substances excretion leads to more realistic fluid behaviour of the extracellular matrix, avoiding artefacts; the stochastic chemostat mode allows comparison of spatially uniform and heterogeneous systems; and the separation of growth kinetics from the individual cell allows condition-dependent switching of metabolism. As an illustration of the model's use, we used the latter feature to study how environmentally fluctuating oxygen availability affects the diversity and composition of a community of denitrifying bacteria that induce the denitrification pathway under anoxic or low oxygen conditions. We tested the hypothesis that the existence of these diverse strategies of denitrification can be explained solely by assuming that faster response incurs higher costs. We found that if the ability to switch metabolic pathways quickly incurs no costs the fastest responder is always the best. However, if there is a trade-off where faster switching incurs higher costs, then there is a strategy with optimal response time for any frequency of environmental fluctuations, suggesting that different types of denitrifying strategies win in different environments. In a single environment, biodiversity of denitrifiers is higher in biofilms than chemostats, higher with than without costs and higher at intermediate frequency of change. The highly modular nature of the new computational model made this case study straightforward to implement, and reflects the sort of novel studies that can easily be executed with the new model.