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.


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.
It is now broadly accepted that bacteria principally exist as surface-associated communities called biofilms 1,2 . Formation and survival of biofilms are a major concern, both in a medical and industrial context [3][4][5] . On the other side, biofilms can also provide useful applications for wastewater treatment 6 and are important for the proper functioning of many ecosystems 7 . The early stages of biofilm development are usually characterized by the formation of tethered small aggregates, so-called microcolonies, either by successive recruitments of new bacteria from the surrounding bulk fluid, multiplication of adhered bacteria or aggregation of bacteria actively moving on a surface 2 . Early microcolonies are comprised of dozens to thousands of cells, are often assembled in matter of hours and have been observed in many different bacteria species 8,9 . Microcolonies represent the first stage of a usually complex development into mature differentiated multicellular biofilms 1 . However, microcolonies are also commonly found by themselves in vivo 8,[10][11][12] . For instance, a recent study has mapped out the advantage that microcolonies could represent for infection of the human pathogen Neisseria meningitidis 13 . Microcolonies thus also play an important role beyond that of an intermediate towards biofilm maturation. The first measurable step towards obtaining cells with different biological functions, the hallmark of differentiation, will be heterogeneous gene expression among the cell population. Understanding the nature of this transition from unicellular to multicellular lifestyle holds the promise for unraveling some of the mechanisms leading to differentiation and controlling biofilm development.
While eukaryotes and prokaryotes are obviously dissimilar, a parallel can be drawn between the development of bacterial biofilm and the development of complex multicellular organisms 2 . The cues that dictate the cascade of signaling that governs multicellular development have long been thought to be of chemical origin but recent studies have brought to light the importance of the mechanical interactions between eukaryotic cells in the differentiation process [14][15][16] . With the knowledge accumulated in the dynamics of eukaryotic cell aggregates in mind, we will explore here the internal dynamics of bacterial microcolonies. In many bacterial species, the events leading to microcolonies are highly regulated with ultimately many different interactions governing the adhesion of bacterial cells to the substrate and between bacterial cells [17][18][19] . Often biofilms are held together by an extracellular matrix composed of DNA, excreted polysaccharides and various bacterial appendages 20,21 . In contrast, the formation of microcolonies of Neisseria gonorrhoeae (Ng) is solely relying on the interactions mediated by a ubiquitous appendage, the Type IV pilus (Tfp) 22,23 . Mutants lacking Tfp are not able to form microcolonies 24 . The unique reliance on Tfp makes Ng an ideal model system to fully understand the dynamics of formation of bacterial microcolonies. In this study, we look experimentally at the dynamics of formation of Ng microcolonies and highlight the crucial role of the mechanical forces generated by retractile Tfp in this process. Our central result is the discovery of emerging heterogeneous behavior within bacterial microcolonies within the first hours of formation. We observe a sharp gradient of bacterial motility from mobile surface layer towards nearly immobile bulk of the microcolony. These results are corroborated by experiments with bacteria incapable of Tfp retraction and comparison with the predictions of the in silico model we recently developed 25 . Ultimately, we see that heterogeneous gene expression follows the heterogeneous motile behavior.

Results and Discussion
Ng microcolonies merge with dynamics consistent with a heterogeneous composition. Tfp are retractile bacterial appendages whose cycles of elongation and retraction enable bacteria to exert forces on their surroundings 23,26 . These polymers have a diameter of molecular size (below 10 nm) and length exceeding the size of the bacteria body (several microns) 23 . An average Ng cell has 10-20 Tfps. Tfp can generate forces up to the nanonewton range when in bundles 27 . In the case of Ng, Tfp are the only motility appendage that the bacteria possess. This leaves the cycles of elongation and retraction of Tfp and the forces that Tfp can exert on their surroundings as the principal agents of microcolony formation. Ng bacteria can form nearly spherical microcolonies of upward to thousands of cells within a few hours, which greatly facilitates their study (See Fig. 1a, Supplementary Movie S1). The active merging of smaller microcolonies into a larger one is the central mechanism responsible for microcolony growth 24 (See Fig. 1a-c, Supplementary Movie S1). We took advantage of the fact that the merger of microcolonies necessitates a complex rearrangement of cells and thus will inform us on the internal dynamics of bacterial microcolonies. To this end, we studied in detail the dynamics of two merging microcolonies. Microcolonies were self-assembled by letting bacteria interact with each other on a surface. Microcolonies of the desired size could be subsequently retrieved and brought into close vicinity and let to interact under a microscope. To quantify the transition of two interacting colonies towards a spherical shape we used the images at the midplane cross section. By fitting an ellipse to the shape of the cross section we measured the aggregate's short and long symmetry axis. Their ratio approaches 1 from below as the colony rounds up (See Figs 2a and S1a,c). Additionally, we measured the height of the "bridge" -a contact area forming between the two touching microcolonies (See Figs 2a and S1a,b). Our analysis shows that the merger occurs with three different dynamical regimes. In the first rapid regime (a few seconds) the two microcolonies are pulled together by retracting pili (retracting  Fig. 2c and Supplementary Movie S2). These findings would be consistent with a heterogeneous microcolony where an envelope of mostly motile cells drives the relaxation in the area with highest surface curvature (at the contact line between microcolonies) and is resisted by a core of much less motile cells responsible for the longer relaxation time. To probe directly for the existence of such "liquid-like" envelope and a more "solid-like" core, we decided to follow the motility of single cells as a function of their position in a microcolony.

Ng microcolonies have an outer layer of motile cells and a core of far less motile cells thus presenting a heterogeneous motility behavior. To look at the motility of single cells within microcolonies
we used a mixture of wild-type (WT) non-labeled cells mixed with a small percentage (5 to 10%) of fluorescently labeled bacteria. We used fluorescent and bright field channels to simultaneously record bacterial tracks and determine their location within the microcolony. To quantify the motility of cells in the microcolonies without being hindered by the global motion of the microcolony, we evaluated the changes in relative distance for pairs of cells at a given depth from the colony surface. The mean squared relative distance increments (MSRD) as a function of time provide information about the motility of cells (See Fig. 3a,b, Supplementary Movie S3 and Methods section) 28 . For short time scales (<1 minute) and for all depths the MSRD grows linearly with time, consistent 2 to the bridge height we were able to estimate the first time scale corresponding to the initial approach of the two colonies and the time scale characterizing the closure of the bridge. The third time scale resulted from a fit to the aspect ratio of the short axis and the long axis of the ellipse and corresponds to the relaxation of the ellipsoidal colony to a spherical shape. with normal diffusive behavior of cells. However, the intensity of the random motion (quantified by the value of the diffusion constant) declines sharply for bacteria at different depth away from the microcolony surface (deeper into the microcolony) ( Fig. 3c,d). These results indicate that Ng bacteria as they assemble into microcolonies create a microenvironment with strikingly different motility behaviors. In order to investigate the nature of the emergence of differential motility within a microcolony we turned to computer modeling. We use a numerical simulation of a bacterial microcolony at the single cell level of detail with explicit pili-pili interactions 25 . This model accounts for known Tfp dynamics and force generation, along with the physical shape of cells and reaches the limits of current computational time constraints (See Methods section and Supplementary Information). The simulations recapitulate the motility behavior observed experimentally: the mean square displacements of bacteria pairs show diffusive motility throughout a microcolony at short time scales, but the value of the associated diffusion constant diminishes greatly for the bacteria tracked deeper towards the core of the colony (See Fig. 3e,f). Importantly, since all cells in our computational study have identical features, the emergence of the heterogeneous behavior highlights the importance to take a dynamical approach when analyzing any aggregation of cells. An energetic approach that would only take into account the difference in adhesion properties between the cells will most likely miss the complex reality that a more precise and dynamical approach might uncover 29 . In the simulations we can correlate the heterogeneous motility of cells with an average number of Tfp bound to pili of other cells and the times of presence and attachment of each pilus. The number of attached pili slowly decreases going from the core of the microcolony to the outer layer as the fluctuations in the number of pili (attached or not) increases (See Fig. S3a,b). Additionally, in the simulations the attachment times rapidly decrease when going from the core to the outer layer (See Fig. S3c). For the cells deep in the core of the colony, even if one of their pili detaches, the multitude of remaining contacts holds the cell in place long enough for the detached pilus to bind again, thus leading to an effectively very slow dynamics in the core. Due to a higher cell density compared to the colony surface, within the bulk the pili density is increased, thus increasing also the re-binding rate of a pilus.
In other words, cells in the core are more likely to be bound to each other than cells on the outside. The qualitative agreement between experiments and simulations indicates that the occurrence of different motility behaviors within a microcolony can be explained solely by the intercellular forces exerted by Tfp. If we use the same model to determine the dynamic of microcolonies merger we obtain similarly good qualitative agreement with the experiments (See Fig. S4). Our model contains just a few free parameters with the other parameters provided by literature (Supplementary Data Table 1), however, while reproducing the data qualitatively it still deviates in numerical values. One observation is that our in silico colonies are less heterogeneous as compared to experiments with a sharper transition. Most likely the discrepancy is due to the simplification of the model allowing for one pilus to have at maximum one contact point with any other pili of other cells. While dictated by computer feasibility it leads to more "dynamic" microcolonies as compared to experiments. Multiple pili contact points in the real setting lead to the formation of active pili network with higher gradient in motility and less dynamic core Recently a different modeling approach averaging the contribution of all pili to an intermittent attractive force 13 similarly led to the existence of a smooth diffusion gradient emphasizing the importance to be able to take into account multiple pili interactions in the future. Importantly, we see that the heterogeneous behavior in the microcolonies emerges during the assembly process of physically identical bacterial cells and this phenomenon is very robust for a wide range of model parameters (see 25 and Supplementary Information).

Non-retracting cells are excluded on the outside of microcolonies. If the forces due to retracting
Tfp play the central role in shaping the microcolony, it would be predicted that cells deprived of the ability to generate those forces would alter the process of colony formation. To test this in the model, we performed simulations of a 1:1 mixture of pulling wild-type cell and cells that possess pili but are unable to retract them. We observed that the non-pulling cells were excluded to the outside of the microcolonies (See Fig. 4c). To find the actual experimental confirmation of this prediction, we used the fact that retraction of Tfp in Ng is under the control of a AAA ATPase pilT 30 . Bacteria deprived of this molecular motor (NgΔpilT) still have Tfp and thus can interact through Tfp-Tfp interactions with other bacteria but they do not retract their pili and thus cannot exert forces on their surroundings 22 . A 1:1 mixture of NgWT cells and NgΔpilT (fluorescently labeled) led, in agreement with the prediction of the computational model, to the sorting of cells that cannot exert forces to the outside of the microcolonies during the self-assembly process (See Fig. 4a,b). Similar sorting of two types of bacteria were observed previously in the case of bacteria with different types of pili and thus different interaction forces between them 31 . Here, the pili are made of the same major pilin in both, pulling and non-pulling bacteria. Thus, the origin of the sorting in our experiments is the absence of retraction forces, not the difference in the interaction forces between different types of pili. A similar discrepancy exists in eukaryotic systems where the difference of adhesion forces between cells has been first postulated to be a driving force in early embryogenesis 32 . More recent theories take also into account the ability of eukaryotic cells to exert forces on their surrounding through contraction of their skeleton. In this case the dynamics of the cells seem crucial to the biological outcome 29,33 . In a control experiment, when the fluorescent markers were used to label two NgWT populations with different fluorophores in a 1:1 mixture, the corresponding forming microcolonies showed a homogeneous repartition of both types of bacteria (See Fig. 4d,e). The homogeneous repartition is also obtained in the case of the simulations (See Fig. 4f). These results exemplify the crucial role that Tfp contractile forces have in shaping Ng microcolonies and demonstrate the predictive power of the computational model.

Heterogeneous gene expression follows heterogeneous motility behavior in microcolonies.
We have demonstrated that Tfp retraction forces in a set of initially identical cells can lead to the formation of a heterogeneous motility behavior as they form a microcolony. To probe whether this motility behavior is associated with differential gene expression across a microcolony we have generated a gene reporter for the pilin gene: the promoter of the pilE gene was incorporated heterologously in the genome of Ng driving a mCherry fluorescent protein (Ng P pilE -mCherry) to get an idea of the expression of pilE across a microcolony. The genomic site of incorporation has been shown to not modify the behavior of the bacteria and the fluorescence will be a proxy for the expression of the pilE gene. When these microcolonies were treated in similar conditions as the ones we used to study the dynamics of microcolonies, the fluorescence reaches the edge of the microcolony and the fluorescence is consistent with a homogeneous expression across the microcolony (See Fig. 5a). Our measurements of metabolic levels among cells within microcolonies during the first few hours of formation show also homogeneous profiles (See Fig. S5a,b). As we have seen previously in the case of Tfp numbers, life and attachment times, the computer model enables us to access quantities that are difficult to measure experimentally. Importantly, we can also compute the mechanical forces within the microcolony. According to the model the dynamical formation of bacterial microcolonies via Tfp interactions will lead not only to a gradient of numbers of Tfp and times of attachments (See Fig. S3a-c) but also to the appearance of a force gradient across a microcolony (See Fig. 5c). So the heterogeneous motility within microcolonies is linked to a force gradient. After a few hours of interactions, the force gradient is present, but the gene expression and metabolism are homogeneous (See Figs 5a and S5a,b). All small diffusible chemical cues are most likely also spatially homogeneous due to the small size of the microcolonies. For bacterial aggregates below 40 microns in diameter diffusion is a very efficient process and both simulations and experiments indicates that oxygen, nutrients and waste products will be homogeneously distributed 34,35 . When microcolonies are left to develop for 7 hours in liquid the pattern of gene expression of pilE becomes heterogeneous as shown by the fluorescence across the microcolony (See Fig. 5b). The pattern of gene expression parallels the pattern of differential motility that we observed. But importantly the pattern of gene expression appears multiple hours after the observation of the motility pattern and the metabolic activity is still homogeneous at that later 7 hour time point (See Fig. S5a,b). These results are indicative that an initial change of motility behavior precedes a differential gene expression within the microcolony. The pursuit of the mechanisms relating these gradients to differential gene expression in a microcolony is beyond the scope of the present study. Nevertheless, keeping in mind the idea from eukaryotic mechanobiology that a gradient of mechanical forces can trigger difference in gene expression, it is tantalizing that the gradient of forces within a bacterial microcolony triggered by the dynamics of Tfp can be the first step towards differentiation as it precedes heterogeneity in gene expression.

Conclusion
The development of multicellular eukaryotic organisms has been a long standing biological question with an overwhelming focus on the presence of gradients of different molecules until recently, when the role of physical forces in development was acknowledged 14 . Bacterial biofilms have been recognized as being similar to multicellular entities and being able to develop differentiated states usually over multiple days 1,2 . The mechanisms at play in this development are the subject of intense scientific inquiry due to the health related importance and ubiquity of biofilms. The role of mechanics in bacteria lifestyle, in particular within biofilm, is being recently appreciated, whether it is the role of motility, hydrodynamical flow or internal forces building within biofilms 36 . Besides, there is now mounting evidence that bacteria have the ability to sense mechanical forces and those forces can in turn change genetic expression 37,38 .
In the case of Neisseria gonorrhoeae, we have shown that a group of identical cells powered by cycles of extensions and retractions of Tfp can self-organize in microcolonies with heterogeneous motility pattern. Such heterogeneous motility behavior correlates with a gradient of mechanical forces within a microcolony. If these mechanical forces can be sensed by bacteria, a feedback mechanism could provide the ability to these forces to take over the control of the spatial development of the microcolony by triggering heterogeneous gene expressions. Ng microcolonies are akin to early developing embryos and eukaryotic cell spheroids, where the forces between cells, both adhesive and contractile, play a crucial role 15,33,[39][40][41] . What might unify these biologically so different systems is the common biophysical mechanism where the mechanical forces generated on the surface of aggregates drive the shape transformation which is resisted by the viscoelastic response of the bulk. This mechanism can apply generally despite the difference in origin of these aggregates whether it is motile aggregation of different cells or successive divisions of a starting cell or a combination of those two modalities. The simple system of Ng bacteria, accompanied by the in silico model, not only enables us to understand better the physiology of an important human disease but it could also give a new insight into the earliest steps of genetic differentiation within a group of identical bacterial cells and ultimately the evolution of multicellularity 42 .

Methods
Bacteria strains and growth conditions. The wild-type (WT) strain used in this study is MS11. The ΔpilT mutant was obtained by an in-frame allelic replacement of the pilT gene by a Kanamycin resistance cassette. Fluorescent proteins (YFP, mCherry or tdTomato) driven by a consensus promoter were incorporated by allelic replacement together with an antibiotic marker (either Kanamycin or Chloramphenicol). Similarly, mCherry driven by the reporter of the pilin gene (370 bp before the beginning of the starting ATG of the pilin ORF) was incorporated by allelic replacement together with a Chloramphenicol marker. Bacteria were grown on GCB-medium base agar plates supplemented with Kellog's supplements at 37 °C and 5% CO 2 . 80 μg/ml of Kanamycin or 7 μg/ml of Chloramphenicol were added when growing mutants with the corresponding antibiotic resistance cassette. Cells were streaked from frozen stock allowed to grow for 24 hours and then lawned onto identical agar plates and used after a 16 to 20 hour growth period. The list of the strains used in this study can be found in the table below: Microcolony formation. Bacteria from lawns on agar plates were resuspended in 1 ml of GCB medium at an optical density of O.D. = 0.7. 100 μl of the suspension was added in the well of the 6 well plate containing 2 ml of GCB medium with a BSA coated coverglass (round 25 mm diameter coverglasses (CS-25R) Warner Instruments) at the bottom or without. The 6 well plate was centrifuged at 1600 × g in a swinging bucket in an 5810 R centrifuge (Eppendorf) for 5 minutes resulting in single bacteria uniformly coating the bottom of the well. For direct imaging the coverglass were transferred to an observation chamber (Attofluor cell chamber, Thermo Fisher Scientific). In the case of mixture the suspension at O.D. = 0.7 of both components of the mixture were prepared and a new 1 ml was prepared by the proper ratio of the two suspension. 100 μl of that new suspension was used similarly to what was described previously.
Microscopy and microcolony merger. All movies were obtained on a Nikon Ti Eclipse inverted microscope equipped for epifluorescence and DIC microscopy and with an optical tweezers setup all under an environmental chamber maintaining temperature, humidity and CO 2 concentration. The objective used is a 60X plan Apo objective. The camera used were either a sCMOS camera (Neo, Andor) or a CMOS USB camera (DCC1240M, Thorlabs). 1 Hz fluorescent movies and 0.1 Hz DIC movies of either microcolony merger or follow up of single cell motility were taken for further analysis. In the case of microcolony merger, microcolonies were preformed and were brought into contact either by optical tweezers or hydrodynamical flow.  Edge detection. In order to detect the edges of single colonies and the merger from DIC data the same algorithm was used. We computed the first derivatives of intensities in x-and y-direction of a Gaussian filtered image and thresholded its absolute value. Afterwards we dilated and eroded the binary image, filled all remaining holes and removed small objects. For all steps internal functions of Matlab were used.

Scanning Electron
Single Cell Tracking inside of Colonies. To track single cells from the fluorescence images, we first computed the center of mass (COM) of the binary shape resulting from the edges and corrected the fluorescence data in order to reduce the effects of the translations of the colony on the tracking algorithm. We interpolated the position of the COM from the DIC data recorded with frequency of 0.1 Hz to match that of the fluorescence data with 1 Hz resolution. A cubic spline data interpolation was applied on the x-and y-component on the COM. We next used the detection and tracking algorithm developed by Blair et al. 45 . For original images we computed the background by applying a large-scale Gaussian filter and substracted its values from the images. Additionally, we used a smaller Gaussian filter for smoothing.
Estimating the MSRD and the diffusion coefficient of cells inside of colonies. The measurements of displacement of individual cells is strongly affected by translation and rotation of the microcolony as a whole. This effect can be circumvented by measuring the absolute value of the distance vector of two individual cells and computing the time-averaged second moment of the relative distance increments, as a function of the time lag. One can show that this quantity is equal to the sum of the mean squared displacements of the two individual cells in the colony [Note: The ensemble averaged mean squared relative displacement coincides with the time averaged quantity as long as the individual cells displacements are smaller than the initial distance between the cells, which is the case for most of pairs of cells in our measurements.]. It thus grows linearly with the lag time and the prefactor of this linear growth quantifies the diffusion of the cells. To identify the location of cells in the microcolony, we computed the distance from the surface by finding the edge shape from DIC images taken in parallel to the fluorescent images and picking the cell pairs that move in a region. By assuming that cells belonging to the same region have the same diffusion coefficient we can read out of individual cells from the fit. The constant offset b accounts for the measurement noise. This method allows us to find the diffusivity of cells as a function of the distance from the surface.
Ellipse fitting. The binary images of the merger movies allowed us to fit an ellipse by computing the central moments, as explained in 46 . This algorithm also allowed us to compute the orientation of a non-spherical colony, the length of the short and long axis of the ellipse and their axis ratio (Fig. 2b).
Bridge height measurement. To compute the height of the bridge forming between the two regions of the binary images we rotated the image so that the colonies were oriented along the x-axis and calculated the COM of the combined regions. Then we moved a line of length L centered around the COM and perpendicular to the axis connecting the two colonies. The bridge was defined as the range for which the whole line could be found inside of the colony region. L was chosen to be small enough to be not affected by the elliptical shape of the colonies. For the late coalescence the results were compared to the short axis of the ellipse.

Simulations. General information.
We implemented the simulation of the model (see Supplementary   Information for a summary of the model) in C++ and used the package OpenMP for parallelization on up to eight CPUs. The simulations were performed on the local computing cluster of the Max Planck Institute for the Physics of Complex Systems consisting of x86-64 GNU/Linux systems.
Single colony modeling. For the simulation of the dynamics of individual cells inside of a colony we randomly initialized 1700 cells inside of a sphere without cell-substrate interactions. First we allowed the cells to repel each other until there was no overlap between neighboring cells. Next, we introduced pili and their dynamics, causing the motion of cells. In order to reproduce the experimental results, we only monitored the cells being positioned <1 µm above and below the midplane and computed the second moment of the relative distance increments (as for the experimental data) of cell pairs, projected to the midplane.
Merger modeling. For the merging simulations we again neglected cell-substrate interactions and initialized two separate colonies each consisting of 1000 randomly distributed cells. The spherical shape of microcolonies indicates cell-substrate interactions are negligible as, if not, they would deform the microcolony towards the substrate 25 .For the first 100 s the pili of the colonies were not allowed to interact with pili of the other colony. This way we allowed for the formation of stable individual microcolonies with a radius of ~7 µm and initial separation of 16 µm between their centers. To compute the bridge height and the eccentricity we created the binary image from the projection of the cell positions and their shapes onto a plane tangential to the axis between the two colonies.
Assembly modeling. For the assembly experiments we distributed 1500 cells on a substrate of the size 98.56 × 98.56 µm 2 with periodic boundary conditions. The pili of a fraction of cells were not able to retract, modeling the ΔpilT mutant.