Multiscale modeling of bacterial colonies: how pili mediate the dynamics of single cells and cellular aggregates

Diclofenac sensitizes multi-drug resistant Acinetobacter baumannii to colistin

Acinetobacter baumannii causes life-threatening infections that are becoming difficult to treat due to increasing rates of multi-drug resistance (MDR) among clinical isolates. This has led the World Health Organization and the CDC to categorize MDR A. baumannii as a top priority for the research and development of new antibiotics. Colistin is the last-resort antibiotic to treat carbapenem-resistant A. baumannii . Not surprisingly, reintroduction of colistin has resulted in the emergence of colistin-resistant strains. Diclofenac is a nonsteroidal anti-inflammatory drug used to treat pain and inflammation associated with arthritis. In this work, we show that diclofenac sensitizes colistin-resistant A. baumannii clinical strains to colistin, in vitro and in a murine model of pneumonia. Diclofenac also reduced the colistin MIC of Klebsiella pneumoniae and Pseudomonas aeruginosa isolates. Transcriptomic and proteomic analyses revealed an upregulation of oxidative stress-related genes and downregulation of type IV pili induced by the combination treatment. Notably, the concentrations of colistin and diclofenac effective in the murine model were substantially lower than those determined in vitro , implying a stronger synergistic effect in vivo compared to in vitro . A pilA mutant strain, lacking the primary component of the type IV pili, became sensitive to colistin in the absence of diclofenac. This suggest that the downregulation of type IV pili is key for the synergistic activity of these drugs in vivo and indicates that colistin and diclofenac exert an anti-virulence effect. Together, these results suggest that the diclofenac can be repurposed with colistin to treat MDR A. baumannii .

Data-driven modelling makes quantitative predictions regarding bacteria surface motility

In this work, we quantitatively compare computer simulations and existing cell tracking data of P. aeruginosa surface motility in order to analyse the underlying motility mechanism. We present a three dimensional twitching motility model, that simulates the extension, retraction and surface association of individual Type IV Pili (TFP), and is informed by recent experimental observations of TFP. Sensitivity analysis is implemented to minimise the number of model parameters, and quantitative estimates for the remaining parameters are inferred from tracking data by approximate Bayesian computation. We argue that the motility mechanism is highly sensitive to experimental conditions. We predict a TFP retraction speed for the tracking data we study that is in a good agreement with experimental results obtained under very similar conditions. Furthermore, we examine whether estimates for biologically important parameters, whose direct experimental determination is challenging, can be inferred directly from tracking data. One example is the width of the distribution of TFP on the bacteria body. We predict that the TFP are broadly distributed over the bacteria pole in both walking and crawling motility types. Moreover, we identified specific configurations of TFP that lead to transitions between walking and crawling states.

Interfacial morphodynamics of proliferating microbial communities

In microbial communities, various cell types often coexist by occupying distinct spatial domains. What determines the shape of the interface between such domains-which in turn influences the interactions between cells and overall community function? Here, we address this question by developing a continuum model of a 2D spatially-structured microbial community with two distinct cell types. We find that, depending on the balance of the different cell proliferation rates and substrate friction coefficients, the interface between domains is either stable and smooth, or unstable and develops finger-like protrusions. We establish quantitative principles describing when these different interfacial behaviors arise, and find good agreement both with the results of previous experimental reports as well as new experiments performed here. Our work thus helps to provide a biophysical basis for understanding the interfacial morphodynamics of proliferating microbial communities, as well as a broader range of proliferating active systems.

Immobilized Stenotrophomonas maltophilia KB2 in Naproxen Degradation

Immobilization is a commonly used method in response to the need to increase the resistance of microorganisms to the toxic effects of xenobiotics. In this study, a plant sponge from Luffa cylindrica was used as a carrier for the immobilization of the Stenotrophomonas maltophilia KB2 strain since such a carrier meets the criteria for high-quality carriers, i.e., low price and biodegradability. The optimal immobilization conditions were established as a temperature of 30 °C, pH 7.2, incubation time of 72 h, and an optical density of the culture of 1.4. The strain immobilized in such conditions was used for the biodegradation of naproxen, and an average rate of degradation of 3.8 µg/hour was obtained under cometabolic conditions with glucose. The obtained results indicate that a microbiological preparation based on immobilized cells on a luffa sponge can be used in bioremediation processes where it is necessary to remove the introduced carrier.

Arrested coalescence of multicellular aggregates

Multicellular aggregates are known to exhibit liquid-like properties. The fusion process of two cell aggregates is commonly studied as the coalescence of two viscous drops. However, tissues are complex materials and can exhibit viscoelastic behaviour. It is known that elastic effects can prevent the complete fusion of two drops, a phenomenon known as arrested coalescence. Here we study this phenomenon in stem cell aggregates and provide a theoretical framework which agrees with the experiments. In addition, agent-based simulations show that active cell fluctuations can control a solid-to-fluid phase transition, revealing that arrested coalescence can be found in the vicinity of an unjamming transition. By analysing the dynamics of the fusion process and combining it with nanoindentation measurements, we obtain the effective viscosity, shear modulus and surface tension of the aggregates. More generally, our work provides a simple, fast and inexpensive method to characterize the mechanical properties of viscoelastic materials.

Mobility of extracellular DNA within gonococcal colonies

Transformation enables bacteria to acquire genetic information from extracellular DNA (eDNA). Close proximity between bacteria in colonies and biofilms may inhibit escape of eDNA from the colony but it also hinders its diffusion between donor and recipient. In this study, we investigate the mobility of DNA within colonies formed by Neisseria gonorrhoeae, and relate it to transformation efficiency. We characterize the penetration dynamics of fluorescent DNA into the colony at a time scale of hours and find that 300 bp fragments diffuse through the colony without hindrance. For DNA length exceeding 3 kbp, a concentration gradient between the edge and the center of the colony develops, indicating hindered diffusion. Accumulation of DNA within the colony increases with increasing DNA length. The presence of the gonococcal DNA uptake sequence (DUS), which mediates specific binding to type 4 pili (T4P) and uptake into the cell, steepens the radial concentration gradient within the colony, suggesting that the DUS reduces DNA mobility. In particular, DNA of N. gonorrhoeae containing multiple DUS is trapped at the periphery. Under conditions, where DUS containing DNA fragments readily enter the colony center, we investigate the efficiency of transformation. We show that despite rapid diffusion of DNA, the transformation is limited to the edge of young colonies. We conclude that DNA mobility depends on DNA length and specific binding mediated by the DUS, resulting in restricted mobility of gonococcal DNA. Yet gonococcal colonies accumulate DNA, and may therefore act as a reservoir for eDNA.

•

DNA fragments encompassing the length of a typical operon efficiently penetrate bacterial colonies.

•

Bacterial colonies accumulate eDNA with an efficiency that depends on the length and the DNA uptake sequence.

•

Genomic DNA from a distinct species spreads efficiently through gonococcal colonies, while gonococcal DNA and DNA from a closely related species are trapped.

•

Transformation is most efficient at the periphery of freshly assembled gonococcal colonies.

Spatiotemporal dynamics of growth and death within spherical bacterial colonies

Bacterial growth within colonies and biofilms is heterogeneous. Local reduction of growth rates has been associated with tolerance against various antibiotics. However, spatial gradients of growth rates are poorly characterized in three-dimensional bacterial colonies. Here, we report two spatially resolved methods for measuring growth rates in bacterial colonies. As bacteria grow and divide, they generate a velocity field that is directly related to the growth rates. We derive profiles of growth rates from the velocity field and show that they are consistent with the profiles obtained by single-cell-counting. Using these methods, we reveal that even small colonies initiated with a few thousand cells of the human pathogen Neisseria gonorrhoeae develop a steep gradient of growth rates within two generations. Furthermore, we show that stringent response decelerates growth inhibition at the colony center. Based on our results, we suggest that aggregation-related growth inhibition can protect gonococci from external stresses even at early biofilm stages.

Competence pili in Streptococcus pneumoniae are highly dynamic structures that retract to promote DNA uptake

The competence pili of transformable Gram-positive species are phylogenetically related to the diverse and widespread class of extracellular filamentous organelles known as type IV pili. In Gram-negative bacteria, type IV pili act through dynamic cycles of extension and retraction to carry out diverse activities including attachment, motility, protein secretion, and DNA uptake. It remains unclear whether competence pili in Gram-positive species exhibit similar dynamic activity, and their mechanism of action for DNA uptake remains unclear. They are hypothesized to either (1) leave transient cavities in the cell wall that facilitate DNA passage, (2) form static adhesins to enrich DNA near the cell surface for subsequent uptake by membrane-embedded transporters, or (3) play an active role in translocating bound DNA via dynamic activity. Here, we use a recently described pilus labeling approach to demonstrate that competence pili in Streptococcus pneumoniae are highly dynamic structures that rapidly extend and retract from the cell surface. By labeling the principal pilus monomer, ComGC, with bulky adducts, we further demonstrate that pilus retraction is essential for natural transformation. Together, our results suggest that Gram-positive competence pili in other species may also be dynamic and retractile structures that play an active role in DNA uptake.

How Physical Interactions Shape Bacterial Biofilms

Biofilms are structured communities formed by a single or multiple microbial species. Within biofilms, bacteria are embedded into extracellular matrix, allowing them to build macroscopic objects. Biofilm structure can respond to environmental changes such as the presence of antibiotics or predators. By adjusting expression levels of surface and extracellular matrix components, bacteria tune cell-to-cell interactions. One major challenge in the field is the fact that these components are very diverse among different species. Deciphering how physical interactions within biofilms are affected by changes in gene expression is a promising approach to obtaining a more unified picture of how bacteria modulate biofilms. This review focuses on recent advances in characterizing attractive and repulsive forces between bacteria in correlation with biofilm structure, dynamics, and spreading. How bacteria control physical interactions to maximize their fitness is an emerging theme.

Antibiotics modulate attractive interactions in bacterial colonies affecting survivability under combined treatment

Biofilm formation protects bacteria from antibiotics. Very little is known about the response of biofilm-dwelling bacteria to antibiotics at the single cell level. Here, we developed a cell-tracking approach to investigate how antibiotics affect structure and dynamics of colonies formed by the human pathogen Neisseria gonorrhoeae. Antibiotics targeting different cellular functions enlarge the cell volumes and modulate within-colony motility. Focusing on azithromycin and ceftriaxone, we identify changes in type 4 pilus (T4P) mediated cell-to-cell attraction as the molecular mechanism for different effects on motility. By using strongly attractive mutant strains, we reveal that the survivability under ceftriaxone treatment depends on motility. Combining our results, we find that sequential treatment with azithromycin and ceftriaxone is synergistic. Taken together, we demonstrate that antibiotics modulate T4P-mediated attractions and hence cell motility and colony fluidity.

On the strong connection between nanoscale adhesion of Yad fimbriae and macroscale attachment of Yad-decorated bacteria to glycosylated, hydrophobic and hydrophilic surfaces

Yad fimbriae are currently viewed as versatile bacterial adhesins able to significantly mediate host or plant-pathogen recognition and contribute to the persistence of Escherichia coli in both the environment and within hosts. To date, however, the underlying adhesion process of Yad fimbriae on surfaces defined by controlled coating chemistries has not been evaluated on the relevant molecular scale. In this work, the interaction forces operational between Yad fimbriae expressed by genetically modified E. coli and self-assembled monolayers (SAM) differing in terms of charge, hydrophobicity or the nature of decorating sugar units are quantified by Single Molecule Force Spectroscopy (SMFS) on the nanoscale. It is found that the adhesion of Yad fimbriae onto probes functionalized with xylose is as strong as that measured with probes decorated with anti-Yad antibodies (ca. 80 to 300 pN). In contrast, the interactions of Yad with galactose, lactose, mannose, -OH, -NH₂, -COOH and -CH₃ terminated SAMs are clearly non-specific. Interpretation of SMFS measurements on the basis of worm-like-chain modeling for polypeptide nanomechanics further leads to the estimates of the maximal extension of Yad fimbriae upon stretching, of their persistence length and of their polydispersity. Finally, we show for the first time a strong correlation between the adhesion properties of Yad-decorated bacteria determined from conventional macroscopic counting methods and the molecular adhesion capacity of Yad fimbriae. This demonstration advocates the effort that should be made to understand on the nanoscale level the interactions between fimbriae and their cognate ligands. The results could further help the design of potential anti-adhesive molecules or surfaces to better fight against the virulence of bacterial pathogens.

Continuum Theory of Active Phase Separation in Cellular Aggregates

Dense cellular aggregates are common in biology, ranging from bacterial biofilms to organoids, cell spheroids, and tumors. Their dynamics, driven by intercellular forces, is intrinsically out of equilibrium. Motivated by bacterial colonies as a model system, we present a continuum theory to study dense, active, cellular aggregates. We describe the process of aggregate formation as an active phase separation phenomenon, while the merging of aggregates is rationalized as a coalescence of viscoelastic droplets where the key timescales are linked to the turnover of the active force. Our theory provides a general framework for studying the rheology and nonequilibrium dynamics of dense cellular aggregates.

Role of Bacterial Cytoskeleton and Other Apparatuses in Cell Communication

The bacterial cytoskeleton is crucial for sensing the external environment and plays a major role in cell to cell communication. There are several other apparatuses such as conjugation tubes, membrane vesicles, and nanotubes used by bacterial cells for communication. The present review article describes the various bacterial cytoskeletal proteins and other apparatuses, the physical structures they form and their role in sensing environmental stress. The implications of this cellular communication in pathogenicity are discussed.

Type IV pili: dynamics, biophysics and functional consequences

The surfaces of many bacteria are decorated with long, exquisitely thin appendages called type IV pili (T4P), dynamic filaments that are rapidly polymerized and depolymerized from a pool of pilin subunits. Cycles of pilus extension, binding and retraction enable T4P to perform a phenomenally diverse array of functions, including twitching motility, DNA uptake and microcolony formation. On the basis of recent developments, a comprehensive understanding is emerging of the molecular architecture of the T4P machinery and the filament it builds, providing mechanistic insights into the assembly and retraction processes. Combined microbiological and biophysical approaches have revealed how T4P dynamics influence self-organization of bacteria, how bacteria respond to external stimuli to regulate T4P activity for directed movement, and the role of T4P retraction in surface sensing. In this Review, we discuss the T4P machine architecture and filament structure and present current molecular models for T4P dynamics, with a particular focus on recent insights into T4P retraction. We also discuss the functional consequences of T4P dynamics, which have important implications for bacterial lifestyle and pathogenesis.

Modelling bacterial twitching in fluid flows: a CFD-DEM approach

Bacterial habitats are often associated with fluid flow environments. Bacterial twitching is important for initial bacterial colonization and biofilm formation. The existing research about bacteria twitching is largely experimental orientated. There is a lack of models of twitching motility of bacteria in shear flows, which could provide fundamental understanding about how bacterial twitching would be affected by bacteria associated properties such as number of pili and their distribution on the cell body and environmental factors such as flow and surface patterns. In this work, a three-dimensional modelling approach of Computational Fluid Dynamics (CFD) coupled with the Discrete Element Method (DEM) proposed to study bacterial twitching on flat and groove surfaces under shear flow conditions. Rod-shaped bacteria are modelled as groups of spherical particles and Type IV pili attached to bacteria are modelled as dynamic springs which can elongate, retract, attach and detach. The CFD-DEM model of rod-shape bacteria is validated against orbiting of immotile bacteria in shear flows. The effects of fluid flow rate and surface topography on twitching motility are studied. The model can successfully predict upstream twitching motility of rod-shaped bacteria in shear flows. Our model can predict that there would be an optimal range of wall shear stress in which bacterial upstream twitching is most efficient. The results also indicate that when bacteria twitch on groove surfaces, they are likely to accumulate around the downstream side of the groove walls.

Motor Properties of PilT-Independent Type 4 Pilus Retraction in Gonococci

Bacterial type 4 pili (T4P) belong to the strongest molecular machines. The gonococcal T4P retraction ATPase PilT supports forces exceeding 100 pN during T4P retraction. Here, we address the question of whether gonococcal T4P retract in the absence of PilT. We show that pilT deletion strains indeed retract their T4P, but the maximum force is reduced to 5 pN. Similarly, the speed of T4P retraction is lower by orders of magnitude compared to that of T4P retraction driven by PilT. Deleting the pilT paralogue pilT2 further reduces the speed of T4P retraction, yet T4P retraction is detectable in the absence of all three pilT paralogues. Furthermore, we show that depletion of proton motive force (PMF) slows but does not inhibit pilT-independent T4P retraction. We conclude that the retraction ATPase is not essential for gonococcal T4P retraction. However, the force generated in the absence of PilT is too low to support important functions of T4P, including twitching motility, fluidization of colonies, and induction of host cell response.IMPORTANCE Bacterial type 4 pili (T4P) have been termed the "Swiss Army knives" of bacteria because they perform numerous functions, including host cell interaction, twitching motility, colony formation, DNA uptake, protein secretion, and surface sensing. The pilus fiber continuously elongates or retracts, and these dynamics are functionally important. Curiously, only a subset of T4P systems employ T4P retraction ATPases to power T4P retraction. Here, we show that one of the strongest T4P machines, the gonococcal T4P, retracts without a retraction ATPase. Biophysical characterization reveals strongly reduced force and speed compared to retraction with ATPase. We propose that bacteria encode retraction ATPases when T4P have to generate high-force-supporting functions like twitching motility, triggering host cell response, or fluidizing colonies.

Substrate-rigidity dependent migration of an idealized twitching bacterium

Mechanical properties of the extracellular matrix are important determinants of cellular migration in diverse processes, such as immune response, wound healing, and cancer metastasis. Moreover, recent studies indicate that even bacterial surface colonization can depend on the mechanics of the substrate. Here, we focus on physical mechanisms that can give rise to substrate-rigidity dependent migration. We study a "twitcher", a cell driven by extension-retraction cycles, to idealize bacteria and perhaps eukaryotic cells that employ a slip-stick mode of motion. The twitcher is asymmetric and always pulls itself forward at its front. Analytical calculations show that the migration speed of a twitcher depends non-linearly on substrate rigidity. For soft substrates, deformations do not lead to build-up of significant force and the migration speed is therefore determined by stochastic adhesion unbinding. For rigid substrates, forced adhesion rupture determines the migration speed. Depending on the force-sensitivity of front and rear adhesions, forced bond rupture implies an increase or a decrease of the migration speed. A requirement for the occurrence of rigidity-dependent stick-slip migration is a "sticky" substrate, with binding rates being an order of magnitude larger than unbinding rates in absence of force. Computer simulations show that small stall forces of the driving machinery lead to a reduced movement on high rigidities, regardless of force-sensitivities of bonds. The simulations also confirm the occurrence of rigidity-dependent migration speed in a generic model for slip-stick migration of cells on a sticky substrate.

How bacterial cells and colonies move on solid substrates

Many bacteria rely on active cell appendages, such as type IV pili, to move over substrates and interact with neighboring cells. Here, we study the motion of individual cells and bacterial colonies, mediated by the collective interactions of multiple pili. It was shown experimentally that the substrate motility of Neisseria gonorrhoeae cells can be described as a persistent random walk with a persistence length that exceeds the mean pili length. Moreover, the persistence length increases for a higher number of pili per cell. With the help of a simple, tractable stochastic model, we test whether a tug of war without directional memory can explain the persistent motion of single Neisseria gonorrhoeae cells. While persistent motion of single cells indeed emerges naturally in the model, a tug of war alone is not capable of explaining the motility of microcolonies, which becomes weaker with increasing colony size. We suggest sliding friction between the microcolonies and the substrate as the missing ingredient. While such friction almost does not affect the general mechanism of single cell motility, it has a strong effect on colony motility. We validate the theoretical predictions by using a three-dimensional computational model that includes explicit details of the pili dynamics, force generation, and geometry of cells.

Type IV Pilin Post-Translational Modifications Modulate Material Properties of Bacterial Colonies

Bacterial type 4 pili (T4P) are extracellular polymers that initiate the formation of microcolonies and biofilms. T4P continuously elongate and retract. These pilus dynamics crucially affect the local order, shape, and fluidity of microcolonies. The major pilin subunit of the T4P bears multiple post-translational modifications. By interfering with different steps of the pilin glycosylation and phosphoform modification pathways, we investigated the effect of pilin post-translational modification on the shape and dynamics of microcolonies formed by Neisseria gonorrhoeae. Deleting the phosphotransferase responsible for phosphoethanolamine modification at residue serine 68 inhibits shape relaxations of microcolonies after perturbation and causes bacteria carrying the phosphoform modification to segregate to the surface of mixed colonies. We relate these mesoscopic phenotypes to increased attractive forces generated by T4P between cells. Moreover, by deleting genes responsible for the pilin glycan structure, we show that the number of saccharides attached at residue serine 63 affects the ratio between surface tension and viscosity and cause sorting between bacteria carrying different pilin glycoforms. We conclude that different pilin post-translational modifications moderately affect the attractive forces between bacteria but have severe effects on the material properties of microcolonies.

Pili mediated intercellular forces shape heterogeneous bacterial microcolonies prior to multicellular differentiation

Microcolonies are aggregates of a few dozen to a few thousand cells exhibited by many bacteria. The formation of microcolonies is a crucial step towards the formation of more mature bacterial communities known as biofilms, but also marks a significant change in bacterial physiology. Within a microcolony, bacteria forgo a single cell lifestyle for a communal lifestyle hallmarked by high cell density and physical interactions between cells potentially altering their behaviour. It is thus crucial to understand how initially identical single cells start to behave differently while assembling in these tight communities. Here we show that cells in the microcolonies formed by the human pathogen Neisseria gonorrhoeae (Ng) present differential motility behaviors within an hour upon colony formation. Observation of merging microcolonies and tracking of single cells within microcolonies reveal a heterogeneous motility behavior: cells close to the surface of the microcolony exhibit a much higher motility compared to cells towards the center. Numerical simulations of a biophysical model for the microcolonies at the single cell level suggest that the emergence of differential behavior within a multicellular microcolony of otherwise identical cells is of mechanical origin. It could suggest a route toward further bacterial differentiation and ultimately mature biofilms.

Molecular Motors Govern Liquidlike Ordering and Fusion Dynamics of Bacterial Colonies

Bacteria can adjust the structure of colonies and biofilms to enhance their survival rate under external stress. Here, we explore the link between bacterial interaction forces and colony structure. We show that the activity of extracellular pilus motors enhances local ordering and accelerates fusion dynamics of bacterial colonies. The radial distribution function of mature colonies shows local fluidlike order. The degree and dynamics of ordering are dependent on motor activity. At a larger scale, the fusion dynamics of two colonies shows liquidlike behavior whereby motor activity strongly affects surface tension and viscosity.

Simulation of colony pattern formation under differential adhesion and cell proliferation

Proliferation of individual cells is one of the hallmarks of living systems, and collectively the cells within a colony or tissue form highly structured patterns, influencing the properties at the population level. We investigate the joint effect of proliferation in the form of cell division and cell sorting due to differential adhesion using a cellular automaton model. Through simulations and theoretical analysis akin to interface growth, we show that this model gives rise to slower than exponential growth in the case of a single cell type as well as novel colony patterns in the case of two cell types. In particular, engulfment of one cell type by the other is strongly enhanced compared to the prediction from the differential adhesion hypothesis in the absence of proliferation. These observations provide new insights in predicting and characterizing colony morphology using experimentally accessible information such as single cell growth rate and cell adhesion strength.

Twitching motility of bacteria with type-IV pili: Fractal walks, first passage time, and their consequences on microcolonies

A human pathogen, Neisseria gonorrhoeae (NG), moves on surfaces by attaching and retracting polymeric structures called Type IV pili. The tug-of-war between the pili results in a two-dimensional stochastic motion called twitching motility. In this paper, with the help of real-time NG trajectories, we develop coarse-grained models for their description. The fractal properties of these trajectories are determined and their influence on first passage time and formation of bacterial microcolonies is studied. Our main observations are as follows: (i) NG performs a fast ballistic walk on small time scales and a slow diffusive walk over long time scales with a long crossover region; (ii) there exists a characteristic persistent length l_{p}^{*}, which yields the fastest growth of bacterial aggregates or biofilms. Our simulations reveal that l_{p}^{*}∼L^{0.6}, where L×L is the surface on which the bacteria move; (iii) the morphologies have distinct fractal characteristics as a consequence of the ballistic and diffusive motion of the constituting bacteria.

Extracellular Polymeric Substance Production and Aggregated Bacteria Colonization Influence the Competition of Microbes in Biofilms

The production of extracellular polymeric substance (EPS) is important for the survival of biofilms. However, EPS production is costly for bacteria and the bacterial strains that produce EPS (EPS+) grow in the same environment as non-producers (EPS-) leading to competition between these strains for nutrients and space. The outcome of this competition is likely to be dependent on factors such as initial attachment, EPS production rate, ambient nutrient levels and quorum sensing. We use an Individual-based Model (IbM) to study the competition between EPS+ and EPS- strains by varying the nature of initial colonizers which can either be in the form of single cells or multicellular aggregates. The microbes with EPS+ characteristics obtain a competitive advantage if they initially colonize the surface as smaller aggregates and are widely spread-out between the cells of EPS-, when both are deposited on the substratum. Furthermore, the results show that quorum sensing-regulated EPS production may significantly reduce the fitness of EPS producers when they initially deposit as aggregates. The results provide insights into how the distribution of bacterial aggregates during initial colonization could be a deciding factor in the competition among different strains in biofilms.