# Differential homotypic and heterotypic interactions of antigen 43 (Ag43) variants in autotransporter-mediated bacterial autoaggregation

Valentin Ageorges, Marion Schiavone, Grégory Jubelin, Nelly Caccia, Philippe

Ruiz, Ingrid Chafsey, Xavier Bailly, Etienne Dague, Sabine Leroy, Jason

Paxman, et al.

► To cite this version:

Valentin Ageorges, Marion Schiavone, Grégory Jubelin, Nelly Caccia, Philippe Ruiz, et al. Differential homotypic and heterotypic interactions of antigen 43 (Ag43) variants in autotransporter-mediated bacterial autoaggregation. Scientific Reports, 2019, 9 (1), pp.11100. 10.1038/s41598-019-47608-4. hal-02335796

HAL Id: hal-02335796

https://hal.science/hal-02335796

Submitted on 28 Oct 2019

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
OPEN

# Differential homotypic and heterotypic interactions of antigen 43 (Ag43) variants in autotransporter-mediated bacterial autoaggregation

Received: 5 March 2019

Accepted: 18 June 2019

Published online: 31 July 2019

Valentin Ageorges¹, Marion Schiavone²,³, Grégory Jubelin¹, Nelly Caccia¹, Philippe Ruiz¹, Ingrid Chafsey¹, Xavier Bailly¹,⁵, Etienne Dague⁶, Sabine Leroy¹, Jason Paxman¹,⁶, Begoña Heras⁶, Frédérique Chaucheyras-Durand¹,², Amanda E. Rossiter⁷, Ian R. Henderson⁷,⁸ & Mickaël Desvaux¹

Antigen 43 (Ag43) is a cell-surface exposed protein of *Escherichia coli* secreted by the Type V, subtype a, secretion system (T5aSS) and belonging to the family of self-associating autotransporters (SAATs). These modular proteins, comprising a cleavable N-terminal signal peptide, a surface-exposed central passenger and an outer membrane C-terminal translocator, self-recognise in a Velcro-like handshake mechanism. A phylogenetic network analysis focusing on the passenger revealed for the first time that they actually distribute into four distinct classes, namely C1, C2, C3 and C4. Structural alignment and modelling analyses demonstrated these classes arose from shuffling of two different subdomains within the Ag43 passengers. Functional analyses revealed that homotypic interactions occur for all Ag43 classes but significant differences in the sedimentation kinetics and aggregation state were present when Ag43C3 was expressed. In contrast, heterotypic interaction occurred in a very limited number of cases. Single cell-force spectroscopy demonstrated the importance of specific as well as nonspecific interactions in mediating Ag43-Ag43 recognition. We propose that structural differences in the subdomains of the Ag43 classes account for different autoaggregation dynamics and propensities to co-interact.

In lipopolysaccharidic (LPS) diderm bacteria (archetypical Gram-negative bacteria), nine protein secretion systems (T1SS-T9SS) are currently recognised as enabling protein transport across the outer membrane (OM), whereas the Sec and Tat export pathways enable protein translocation across the cytoplasmic membrane, also called the inner membrane (IM)¹⁻³. Among those, the T5SS represents the simplest and most widely used mechanism for protein secretion in LPS-diderm bacteria⁴⁻⁷, where proteins are translocated across the OM via a β-barrel with the assistance of some periplasmic chaperones (e.g. Skp, SurA or DegP), the BAM (β-barrel assembly machinery) and TAM (translocation and assembly module) complexes along the process⁸⁻¹¹. Based on phylogenetic analyses and structural features, the T5SS is subdivided into the (i) classical autotransporters (T5aSS), (ii) two-partner passenger-translocators (T5bSS), (iii) trimeric autotransporters (T5cSS), (iv) hybrid autotransporters (T5dSS), and (v) inverted autotransporters (T5eSS)¹¹⁻¹⁴.

¹Université Clermont Auvergne, INRA, UMR454 MEDI5, 63000, Clermont-Ferrand, France. ²Lallemand Animal Nutrition, 31702, Blagnac, Cedex, France. ³LISBP, Université de Toulouse, CNRS UMR5504, INRA UMR792, INSA, 31077, Toulouse, Cedex, France. ⁴LAAS-CNRS, Université de Toulouse, 31077, Toulouse, Cedex, France. ⁵INRA UR346 Epidémiologie Animale, 63122, Saint Genès Champanelle, France. ⁶Department of Biochemistry and Genetics, La Trobe Institute for Molecular Science, La Trobe University, Bundoora, Victoria, 3086, Australia. ⁷Institute of Microbiology and Infection, University of Birmingham, B152TT, Birmingham, United Kingdom. ⁸Institute for Molecular Biosciences, University of Queensland, St. Lucia, Queensland, 4067, Australia. Correspondence and requests for materials should be addressed to M.D. (email: mickael.desvaux@inra.fr)
A classical autotransporter (AT) is a modular protein constituted of three main regions, i.e. (i) a N-terminal signal peptide (SP), (ii) a central passenger, and (iii) a translocator at the C-terminus<sup>15,16</sup>. As such, AT export via the Sec translocase is mediated by the SP before the passenger is secreted through the β-barrel pore formed by the translocator in the OM. Once secreted, the passenger can have three different fates as it may be (i) displayed at the bacterial cell surface, (ii) cleaved from the translocator and released as an exoprotein into the extracellular milieu, or (iii) processed but remain associated to the cognate translocator<sup>10</sup>. The functions of ATs are quite diverse and can be broadly categorised as (i) self-associating ATs (SAATs), (ii) protease ATs, (iii) phosphatase/hydrolase ATs, (iv) lipase/esterase ATs, (v) vacuolating ATs, (vi) adhesin ATs, and (vii) ATs of unknown function<sup>3,10</sup>. SAATs promote bacterial aggregation by cell self-association, a phenomenon called autoaggregation<sup>17,18</sup>. Besides floculation and sedimentation of bacterial cells, autoaggregation is an important mechanism with consequences for bacterial virulence and infection due to biofilm formation, protection from environmental stresses, antibiotic resistance and defence against predators<sup>18</sup>. The antigen 43 (Ag43) is a prototypical member of the SAATs; this category of ATs also include AIDA (adhesin involved in diffuse adherence) and TibA (toxigenic invasion locus b)<sup>17</sup>. Besides homotypic interactions, heterotypic interactions were demonstrated between AIDA-Ag43, AIDA-TibA and Ag43-TibA<sup>17,19</sup>. By promoting interspecies cell-to-cell contact, these interactions would have important ecophysiological roles in inducing co-operative multispecies biofilms that allows for more successful microbial colonisation of various niches<sup>20</sup>.

Besides the basic and common modular architecture of an AT, Ag43 exhibits some additional features. At the N-terminus, Ag43 displays a highly conserved extended SP region (ESPR) that slows the rate of export across the IM through the Sec-YidC pathway to prevent protein misfolding in the periplasm and consequently sustain the secretion-competent state of the passenger across the OM through the translocator<sup>21,22</sup>. Of note, the ESPR is not exclusive to Ag43 as it occurs in many different, but not all, ATs. N-terminal to the translocator, Ag43 exhibits a type 1 autochaperone domain (AC1), which is considered as essential for folding of β-helical passengers<sup>3,23</sup>. Like most ATs, the Ag43 passenger folds into a single-stranded right-handed parallel β-helix with triangular-shaped coil cross-sections<sup>3,24</sup>. Most recently, the crystal structure of the Ag43 from uropathogenic *E. coli* (UPEC) CFT073 (Ag43a) was solved and revealed the passenger displayed a twisted L-shape β-helical architecture<sup>25</sup>. The L-shaped conformation of the passenger further appeared as a requirement for driving the formation of cell aggregates and self-association in a Velcro-like handshake mechanism. While the function of Ag43 was primarily investigated in *Escherichia coli* K12<sup>24</sup>, its implication in autoaggregation, adhesion and biofilm formation was further confirmed in some pathogenic species, namely enterohaemorrhagic *E. coli* (EHEC) and UPEC<sup>26,27</sup>. The Ag43 is encoded by two identical alleles in EHEC O157:H7 EDL933 but by two different alleles in UPEC CFT073; the UPEC Ag43a variant mediates a strong autoaggregation and biofilm formation phenotypes, whereas the Ag43b variant is a less efficient aggregating and biofilm forming factor<sup>27</sup>. Phylogenetic analyses had revealed the existence of two sub-families (SF) of Ag43, namely SF-I and SF-II<sup>24</sup>.

While homotypic Ag43 interactions have previously been noted, heterotypic interactions of Ag43 variants had not been extensively characterised as yet, it prompted us to investigate these aspects in a systematic study. In a cooperative evolution hypothesis that “birds of a feather flock together”, different combinations of Ag43 variants might select bacterial interactions with siblings to the exclusion of others. From phylogenetic analyses of Ag43 passengers, it first appeared that the grouping was not dichotomous but clustered into four different classes associated to structural differences. Following the expression of representative Ag43 members, different propensities of the variants to induce autoaggregation were evidenced. By monitoring the formation of cell aggregates by epifluorescence microscopy and flow cytometry, homotypic and heterotypic interactions between different Ag43 variants were for the first time identified and characterised. Moreover, biophysical studies by atomic force microscopy (AFM) revealed different dynamics of Ag43 interactions suggesting that structural differences in the N-terminal and C-terminal subdomains of the Ag43 passengers can modulate the association of the different variants.

## Results

**Ag43 passengers cluster into four distinct phylogenetic classes with combination of two different subdomains.** The functional properties of the modular protein Ag43 are determined by the passenger, which protrudes from the bacterial OM upon translocation. To apprehend the diversity of Ag43, a phylogenetic network analysis of the passenger was applied to a dataset of 106 distinct passengers generated from previously reported Ag43 sequences in the literature<sup>24,28</sup> (Supplementary Table S1). The network obtained defined for four main classes of passengers (C1 to C4) supported by robust splits (Fig. 1A).

C1 harbours the Ag43 passenger from *E. coli* K12 MG1655 (UPI00003B25EE), which is among the first identified and the most investigated Ag43<sup>24</sup>. C2 includes the Ag43 passenger from EHEC O157:H7 EDL933 (UPI00000D07E5), which is the first Ag43 investigated in a pathogenic *E. coli*<sup>26</sup>. C3 harbours the Ag43a (UPI00000E4A66) and Ag43b (UPI00000E49D2) from UPEC CFT073, the former being the first and only structurally characterised Ag43 passenger (PDB: 4KH3)<sup>25</sup>. To date, no Ag43 homologues from C4 have been functionally or structurally characterised. The original SF-I grouping corresponds to the gathering of C1 and C3, whereas SF-II grouping reflects C2 and C4<sup>24</sup>. The initial clustering correlated with the most obvious difference among the Ag43 passengers, their sizes; while Ag43 passengers from C1 and C3 contain around 500 amino acids, C2 and C4 homologues are smaller with about 430 amino acids (Supplementary Table S1).

Our updated phylogenetic network further revealed a strong ambiguity in the relatedness patterns among Ag43 passengers from C1, C2, C3, and C4, illustrated by the diamond shape of the network with almost no specific divergence. Therefore, a more detailed analysis of subdomains of the passengers was performed to identify the source of this ambiguity. Homologies to pectin lyase fold (IPR011050) overlapping with P22 tailspike protein (IPR012332) and AT adhesin (IPR030930) domains were systematically predicted for all Ag43 passengers, highlighting the β-helix topology encountered in SAATs (Supplementary Table S1). Modelling of the 3D structure
<figure><img src="image_2.png" /><figcaption>**Figure 1.** Phylogenetic analysis of Ag43 passengers and their subdomains. The amino acid sequences of the passengers were aligned using T-Coffee and structural information of Ag43<sup>C3-UP100000E4A66</sup> (PDB: 4KH3) as template. Based on the L-shaped Ag43 structure (Fig. 2), the amino acid sequences of the SL, EJ and BL subdomains were recovered for alignment. Trees were generated using SplitsTree applying bootstrapping and 95% confidence for the observed nodes. Inserts represent the trees before applying bootstrapping and confidence threshold as provided in Supplementary Fig. S1. (A) Tree of the Ag43 passengers highlighting the four classes, C1, C2, C3 and C4. (B) Tree of the SL subdomain highlighting the two subtypes, SL1 and SL2. (C) Tree of the EJ subdomain highlighting the two subtypes, EJ1 and EJ2. (D) Tree of the BL subdomain highlighting the two subtypes, BL1 and BL2. For clarity, some sequences were removed from the trees but are readily available in Supplementary Table S1 and Fig. S1. The scale bars represent the number of substitutions per site.</figcaption></figure>

of the Ag43 passengers confirmed that they are all predicted to exhibit a right-handed parallel β-helix fold with an L-shaped structure as originally observed in Ag43<sup>C3-UP100000E4A625</sup> (Fig. 2A). Following the analysis of the coil cross-sections along the Ag43 passengers<sup>29</sup>, it appeared the N-terminal part, namely the SL (for stem of the letter L) region, systematically exhibited triangular (T) β-solenoids of subtype 4 (T4) (Fig. 2 and Supplementary Table S1). A phylogenetic network analysis of the SL subdomain identified two robust and unambiguously differentiated clusters of sequences, SL of type 1 (SL1) and SL2 (Fig. 1B). The two subtypes were related by a single branch, and remaining splits were located within SL2. Sequences from SL1 and SL2 subdomains had mean percentages of similarity of 93 and 88%, respectively, but these percentages fell at 75% when comparing SL1 versus SL2 (Supplementary Table S1). While SL1 clustered subdomains corresponding to Ag43<sup>C1</sup> and Ag43<sup>C4</sup> passengers, SL2 clustered subdomains from Ag43<sup>C2</sup> and Ag43<sup>C3</sup> passengers. Nevertheless, the SL region appeared slightly longer by 15–17 amino acids in Ag43<sup>C1</sup> and Ag43<sup>C3</sup> passengers.

The elbow region joining (EJ) the upper and lower parts of the L-shape was always bordered by β-hairpin motifs with three-rung parallel T5 β-solenoids in-between, which exhibited a highly conserved sequence in the first rung, namely VVLE[N/S]GGRLDVL. The EJ region reveals a remarkably conserved size of 85 amino acids. The phylogenetic network analysis of the EJ subdomain showed a different relatedness pattern compared to the SL subdomain. Two robust sequence subtypes were clustered, namely (i) EJ of type 1 (EJ1), which included Ag43 passenger subdomains from C1 and C3, and (ii) EJ2, which clustered subdomains from the C2 and C4 (Fig. 1C). In average, sequences from EJ1 and EJ2 subdomains had percentages of similarity of 97 and 94%, respectively, but these percentages fell at 53% when comparing EJ1 versus EJ2 (Supplementary Table S1).
<figure><img src="image_3.png" /><figcaption>**Figure 2.** Modular structure of the Ag43 passengers belonging to classes 1, 2, 3 and 4. (A) Based on the experimentally determined L-shaped structure of Ag43<sup>C3-UPI00000E4A66</sup> (PDB: 4KH3), three subdomains could be defined SL (for stem of the letter L), EJ (for elbow joint), and BL (for bottom of the letter L)<sup>25</sup>. Both SL and BL subdomains exhibit a right-handed parallel β-helix, the cross-sections of each rung along SL form triangular (T) β-solenoids of subtype 4 (T4) and of T5-type β-solenoids for BL (Supplementary Table S1). The EJ subdomain is flanked by β-hairpin motifs with three-rung parallel T5 β-solenoids. (B) Molecular modelling of the structures of Ag43 passengers from C1 to C4 using Phyre<sup>54,55</sup>. The SL domain is shown in shades of blue, with electric blue corresponding to SL1 and cyan to SL2. The EJ domain is depicted in shades of orange, with bright orange corresponding to EJ1 and light orange to EJ2. The BL domain is represented in shades of red, with bright red corresponding to BL1 and magenta to BL2. PT: prototypical variant.</figcaption></figure>

The BL (for bottom of the letter L) subdomain in the C-terminal part of the passengers is composed of T5-type β-solenoids (Fig. 2 and Supplementary Table S1). The phylogenetic network analysis of the BL was coherent with the analysis of the EJ subdomain. Two robust sequence subtypes were clustered: (i) BL of type 1 (BL1), which included Ag43 passengers subdomains from C1 and C3, and (ii) BL2, which clustered subdomains from C2 and C4 (Fig. 1D). Sequences from BL1 and BL2 subdomains had percentages of similarity from and above 95 and 93% (mean values of 99 and 98%), respectively, but these percentages fell below 30% (28% in average) when comparing SL1 versus SL2 (Supplementary Table S1). Furthermore, a difference of 53 amino-acids among the length of BL subdomains was observed. It accounted for most of the observed variation in the length of Ag43 passengers between C1/C3 over C2/C4 (Supplementary Table S1). This difference indeed resulted from a truncation at the C-terminus of BL2 compared to BL1.

These data suggest that the two SL phylogenetic clusters and the two coherent phylogenetic clusters of EJ and BL were shuffled during the evolution of Ag43 proteins explaining the observed ambiguity in relatedness pattern of the passengers (Fig. 1A). It resulted in four combinations of Ag43 passengers subdomains: (i) SL1-EJ1-BL1 corresponding to Ag43<sup>C1</sup>, (ii) SL2-EJ2-BL2 corresponding to Ag43<sup>C2</sup>, (iii) SL2-EJ1-BL1 corresponding to Ag43<sup>C3</sup>, and (iv) SL1-EJ2-BL2 corresponding to Ag43<sup>C4</sup>.

**Autoaggregation occurs through homotypic interactions for all Ag43 variants but the size of the aggregates is significantly much smaller with Ag43<sup>C3</sup>.** Our new classification of Ag43 delineated the proteins into four classes based on sequence similarity and consequently structural differences. We hypothesised that these variations would result in functional differences. Therefore, we investigated the autoaggregation phenotype of three selected Ag43 homologues from each of the four classes (Fig. 2 and Table 1) when expressed in the non-aggregative *E. coli* SRUD<sup>30</sup>. When bacterial cell cultures were left standing statically, all bacterial cells expressing the Ag43 flocculated and sedimented at the bottom of tubes contrary to the control SRUD pBAD (with empty vector pBAB/Myc-HisA, i.e. Ag43<sup>-</sup>) (Fig. 3A,B) (Table 1). No significant differences could be observed for the sedimentation kinetics of bacterial cells expressing Ag43 from C1, C2 and C4 with all three having compact sedimentation pellets. Microscopic observations indicated the flocculation and sedimentation of the bacterial cultures resulted from bacterial cell autoaggregation as no cell aggregates could be observed with SRUD pBAD, control (Fig. 3C).

In contrast, strains expressing the three C3 homologues (Fig. 2), namely Ag43<sup>C3-UPI00000E4A66</sup>, Ag43<sup>C3-UPI00000E49D2</sup> and Ag43<sup>C3-UPI0001EB6ED2</sup>, showed a significant delay in the autoaggregation kinetics compared to other Ag3 expressing strains (Fig. 3A). In addition, the sedimentation pellets were looser and cell aggregates appeared smaller, with the presence of individual single cells (Fig. 3B,C). Within C3, some significant differences were observed as bacterial cultures of SRUD expressing Ag43<sup>C3-UPI00000E4A66</sup> sedimented significantly quicker than
<table><caption>**Table 1.** Genotype and phenotype of most relevant *E. coli* strains used in the study. <sup>a</sup>The complete list of *E. coli* strains and plasmids used in this study is available in Supplementary Table S2. PT: prototypical variant.</caption><thead><tr><th>Name of E. coli strain<sup>a</sup></th><th>Relevant genotype and/or phenotype</th></tr></thead><tbody><tr><td>SRUD</td><td>E. coli K12 MG1655 Δ fim Δ agn43<sup>30</sup></td></tr><tr><td>SRUD-Green pBAD</td><td>Empty pBAD/Myc-HisA in E. coli SRUD-Green</td></tr><tr><td>SRUD-Red pBAD</td><td>Empty pBAD/Myc-HisA in E. coli SRUD-Red</td></tr><tr><td>SRUD-Blue pBAD</td><td>Empty pBAD/Myc-HisA in E. coli SRUD-Blue</td></tr><tr><td>SRUD Ag43<sup>C1-PT</sup></td><td>pBAD-Ag43<sup>C1-UP100003B25EE</sup> in E. coli SRUD</td></tr><tr><td>SRUD-Green Ag43<sup>C1-PT</sup></td><td>pBAD-Ag43<sup>C1-UP100003B25EE</sup> in E. coli SRUD-Green</td></tr><tr><td>SRUD-Red Ag43<sup>C1-PT</sup></td><td>pBAD-Ag43<sup>C1-UP100003B25EE</sup> in E. coli SRUD-Red</td></tr><tr><td>SRUD-Blue Ag43<sup>C1-PT</sup></td><td>pBAD-Ag43<sup>C1-UP100003B25EE</sup> in E. coli SRUD-Blue</td></tr><tr><td>SRUD Ag43<sup>C2-PT</sup></td><td>pBAD-Ag43<sup>C2-UP100000D07E5</sup> in E. coli SRUD</td></tr><tr><td>SRUD-Green Ag43<sup>C2-PT</sup></td><td>pBAD-Ag43<sup>C2-UP100000D07E5</sup> in E. coli SRUD-Green</td></tr><tr><td>SRUD-Red Ag43<sup>C2-PT</sup></td><td>pBAD-Ag43<sup>C2-UP100000D07E5</sup> in E. coli SRUD-Red</td></tr><tr><td>SRUD-Blue Ag43<sup>C2-PT</sup></td><td>pBAD-Ag43<sup>C2-UP100000D07E5</sup> in E. coli SRUD-Blue</td></tr><tr><td>SRUD Ag43<sup>C3-PT</sup></td><td>pBAD-Ag43<sup>C3-UP100000E49D2</sup> in E. coli SRUD</td></tr><tr><td>SRUD-Green Ag43<sup>C3-PT</sup></td><td>pBAD-Ag43<sup>C3-UP100000E49D2</sup> in E. coli SRUD-Green</td></tr><tr><td>SRUD-Red Ag43<sup>C3-PT</sup></td><td>pBAD-Ag43<sup>C3-UP100000E49D2</sup> in E. coli SRUD-Red</td></tr><tr><td>SRUD-Blue Ag43<sup>C3-PT</sup></td><td>pBAD-Ag43<sup>C3-UP100000E49D2</sup> in E. coli SRUD-Blue</td></tr><tr><td>SRUD Ag43<sup>C4-PT</sup></td><td>pBAD-Ag43<sup>C4-UP10001E6D993</sup> in E. coli SRUD</td></tr><tr><td>SRUD-Green Ag43<sup>C4-PT</sup></td><td>pBAD-Ag43<sup>C4-UP10001E6D993</sup> in E. coli SRUD-Green</td></tr><tr><td>SRUD-Red Ag43<sup>C4-PT</sup></td><td>pBAD-Ag43<sup>C4-UP10001E6D993</sup> in E. coli SRUD-Red</td></tr><tr><td>SRUD-Blue Ag43<sup>C4-PT</sup></td><td>pBAD-Ag43<sup>C4-UP10001E6D993</sup> in E. coli SRUD-Blue</td></tr></tbody></table>

those expressing Ag43<sup>C3-UP100000E49D2</sup> or Ag43<sup>C3-UP10001E86ED2</sup>. To confirm this slower autoaggregation phenotype for bacterial cells expressing Ag43<sup>C3</sup>, three additional homologues from this class were tested; in each case, autoaggregation was slower (Fig. 3A).

To further characterise the homotypic interactions that are responsible for the autoaggregation phenotype, one prototypical (PT) representative was selected from each class (Table 1), namely Ag43<sup>C1-PT</sup> (UPI00003B25EE), Ag43<sup>C2-PT</sup> (UPI00000D07E5), Ag43<sup>C3-PT</sup> (UPI00000E49D2) and Ag43<sup>C4-PT</sup> (UPI0001E6D993). Fluorescence microscopy revealed that in the absence of Ag43 (Ag43<sup>−</sup>, i.e. SRUD-Green pBAD and SRUD-Red pBAD) no cell aggregates could be observed (Fig. 4A). In contrast aggregates were detected for cells expressing Ag43 variants from each group (Fig. 4B–E). The combination of Ag43<sup>−</sup> red-fluorescent-tagged cells mixed with Ag43<sup>+</sup> green-fluorescent-tagged cells (i.e. SRUD-Green pBAD-Ag43) systematically exhibited only one type of bacterial cell aggregate, consisting uniquely of green-fluorescent-tagged cells where only a few Ag43<sup>−</sup> cells could be observed as trapped within the aggregates (Fig. 4B–E). When green- and red-fluorescent cells expressing the same Ag43 were mixed equally, both fluorescent cells were found in the aggregates (Fig. 4F–I). These results confirmed that all tested Ag43 proteins, Ag43<sup>C1-PT</sup>, Ag43<sup>C2-PT</sup>, Ag43<sup>C3-PT</sup> or Ag43<sup>C4-PT</sup>, enable homotypic interactions inducing bacterial autoaggregation. However, examination of the obtained micrographs seemed to indicate that more individual cells were present when Ag43<sup>C3-PT</sup> was expressed (Fig. 4H).

To gain quantitative insight, fluorescent nonaggregative and aggregative *E. coli* cells were analysed using flow cytometry. SSC-A (side scatter area) and FSC-A (forward scatter area) parameters are commonly used to discriminate different types of cells within a heterogeneous population with respect to granularity and cell size, respectively. Here, SSC-A over FSC-A dot plots allowed differentiating single cells from aggregated cells. As shown in Supplementary Fig. S3, flow cytometry analysis of the non-aggregative *E. coli* SRUD pBAD revealed a FCS-A/SSC-A pattern that is clearly distinct from the one obtained using the aggregative *E. coli* SRUD expressing Ag43<sup>C1-PT</sup>. Comparison of these 2 samples allowed discriminating single cells from cell aggregates, the latter forming a comet tail in the upper part of the SSC-A over FSC-A representation (Supplementary Fig. S3). In agreement with microscopic examination, no cell aggregates were observed in the absence of Ag43, whereas cells expressing Ag43<sup>C1-PT</sup>, Ag43<sup>C2-PT</sup>, Ag43<sup>C3-PT</sup> and Ag43<sup>C4-PT</sup> formed aggregates (Supplementary Fig. S4). Furthermore, when green- with blue-fluorescent-tagged SRUD cells expressing the same Ag43 allele were combined, mixed aggregates formed (Fig. 5A); these represented around 90% of the total aggregates (Fig. 5B). As determined by flow cytometry analyses, quantification of the concentration of single individual cells showed a significant difference for cells expressing Ag43<sup>C3-PT</sup> compared to other Ag43 variants (Fig. 5C). Further, the shorter comet tail for *E. coli* SRUD Ag43<sup>C3-PT</sup> was indicative of smaller aggregates (Fig. 5A) and the mean size of the aggregates was confirmed to be significantly smaller for cells expressing Ag43<sup>C3-PT</sup> compared to any of the other Ag43 variants (Fig. 5D). To ensure the differences reported for Ag43<sup>C3-PT</sup> were not influenced by the protein expression levels, Western blot analyses were performed and indicated the Ag43 variants were expressed at similar levels in *E. coli* SRUD and displayed at the OM (Fig. S5).

**Autoaggregation occurs only for some combinations of heterotypic interactions.** To study heterotypic interactions, all six possible combinations of interactions between the prototypical Ag43 proteins from C1 to C4 were tested (Fig. 6). Fluorescent microscopic examinations revealed three major trends (i) large clumps
<figure><img src="image_4.png" /><figcaption>**Figure 3.** Autoaggregation of *E. coli* SRUD expressing different Ag43 variants. (A) The sedimentation assay was performed in LB at 37°C over 6 h with SRUD expressing variants of Ag43 from C1 (shades of green), Ag43<sup>C2</sup> (shades of grey), (iii) Ag43<sup>C3</sup> (shades of red), Ag43<sup>C4</sup> (shades of blue). The nonaggregative SRUD pBAD was used as negative control. (B) Sedimentation was visualised after 6 h of static incubation. (C) Phase-contrast microscopy observation from sampling at 6 h of static incubation. PT: prototypical variant (Table 1).</figcaption></figure>

of bacteria containing only red- or green-fluorescent-tagged cells, suggesting only homotypic interactions occur (Fig. 6A,B), (ii) large aggregates containing only one type of fluorescent-tagged cell with smaller aggregates of the other tagged cells associated with the periphery, suggesting a weak propensity for heterotypic interactions (Fig. 6C,D), and (iii) clumps containing both green and red fluorescent cells suggesting the occurrence of heterotypic interactions (Fig. 6E,F).

Flow cytometry analyses confirmed the absence of heterotypic interaction between cells expressing Ag43<sup>C1-PT</sup> and Ag43<sup>C2-PT</sup> or cells expressing Ag43<sup>C2-PT</sup> and Ag43<sup>C4-PT</sup> (Fig. 7A) as the percentage of mixed aggregates were systematically estimated at less than 3% (Fig. 7B). Similarly, flow cytometry analyses confirmed the absence of significant heterotypic interactions as observed by epifluorescent microscopy for combinations of cells expressing Ag43<sup>C1-PT</sup>/Ag43<sup>C3-PT</sup> or Ag43<sup>C2-PT</sup>/Ag43<sup>C3-PT</sup>, with a percentage of mixed aggregates estimated at less than 5% (Fig. 7A,B). However, combinations of cells expressing Ag43<sup>C1-PT</sup>/Ag43<sup>C4-PT</sup> or Ag43<sup>C3-PT</sup>/Ag43<sup>C4-PT</sup> showed significant heterotypic interactions with mixed cell aggregates estimated at just above 50% and 12%, respectively (Fig. 7A,B). In both cases, the percentage of mixed cell aggregates were significantly lower than the corresponding homotypic combinations (Fig. 5B). For all heterotypic combinations, no significant differences could be observed in the concentration of single cells (Fig. 7C). Heterotypic combinations resulted in aggregates with different sizes, some showing a mean size significantly larger than the corresponding homotypic combinations (i.e. Ag43<sup>C1-PT</sup>/Ag43<sup>C2-PT</sup>, Ag43<sup>C1-PT</sup>/Ag43<sup>C3-PT</sup>, and Ag43<sup>C1-PT</sup>/Ag43<sup>C4-PT</sup>) (Figs 5D and 7D). Interestingly, the association of cells expressing Ag43<sup>C3-PT</sup> with cells expressing any other class of Ag43 resulted in significantly lower concentration of single individual cells (Figs 5C and 7C) as well as significantly larger cell aggregates (Figs 5D and 7D) when compared to homotypic interactions.

Overall, heterotypic interactions between Ag43 from different classes is rather compromised, since no heterotypic interactions were observed in four out of the six possible combinations of prototypical Ag43 from C1 to C4. When heterotypic interactions occurred, the percentages of mixed aggregates were significantly lower than for homotypic interactions. Nonetheless, the highest percentage of mixed aggregates resulted from cells expressing Ag43<sup>C1-PT</sup>/Ag43<sup>C4-PT</sup>, both passengers consisting of the SL1 subdomain (Fig. 2B). Conversely, and despite exhibiting the same SL2 subdomain, no heterotypic interactions could be observed for cells expressing Ag43<sup>C2-PT</sup>/Ag43<sup>C3-PT</sup>. Surprisingly, some heterotypic interactions were observed between cells expressing Ag43<sup>C3-PT</sup>/Ag43<sup>C4-PT</sup>, which would involve interactions between SL1 and SL2 subdomains.
<figure><img src="image_5.png" /><figcaption>**Figure 4.** Fluorescent microscopy analysis of homotypic interactions of bacterial cells expressing different variants of Ag43. Epifluorescence microscopic observations were performed from bacterial cells grown in LB and after 6 h of static incubation at 37 °C. Co-culture of nonaggregative (NA) SRUD-Red pBAD with (A) SRUD-Green pBAD, (B) SRUD-Green Ag43C1-PT, (C) SRUD-Green Ag43C2-PT, (D) SRUD-Green Ag43C3-PT, or (E) SRUD-Green Ag43C4-PT. Co-culture of SRUD-Red and SRUD-Green both expressing (F) Ag43C1-PT, (G) Ag43C2-PT, (H) Ag43C3-PT, or (I) Ag43C4-PT. The NA, C1, C2, C3 and C4 labels match the colour of the fluorescent bacterial strains.</figcaption></figure>

**Ag43-mediated bacterial autoaggregation involves both specific and nonspecific interactions.**

Our data suggest the Ag43-Ag43 interactions could be stronger or weaker for some pairings. To quantify this and reinforce the hypothesis, single cell-force spectroscopy experiments were carried out to gain further insight
<figure><img src="image_6.png" /><figcaption>**Figure 5.** Characterisation of homotypic interactions of Ag43 variants by flow cytometry. Flow cytometry analyses were performed by mixing SRUD-Green and SRUD-Blue both expressing the same Ag43, namely Ag43<sup>C1-PT</sup> (C1), Ag43<sup>C2-PT</sup> (C2), Ag43<sup>C3-PT</sup> (C3) or Ag43<sup>C4-PT</sup> (C4). (A) Based on the SSC-A over FSC-A representation (pictured in the insert), the fluorescent FITC-H over PB450-H representation focused on the gate corresponding to the cell aggregates as previously defined (Supplementary Fig. S3). The green cloud corresponds to cell aggregates of SRUD-Green, the blue cloud to cell aggregates of SRUD-Blue, and the orange cloud to mixed cell aggregates of SRUD-Green and SRUD-Blue. Graphs correspond to one representative experiment. (B) The percentage of mixed cell aggregates was calculated from the events detected in orange cloud over the total number of events corresponding to aggregates. (C) The concentration of cells was calculated from the gates corresponding to single individual cells as previously defined (Supplementary Fig. S3). (D) The size of aggregates was recovered from the gate corresponding to cell aggregates as previously defined (Supplementary Fig. S3).</figcaption></figure>

in the adhesive forces at play in the autoaggregation process. First, the self-interaction between cells expressing different Ag43 variants was investigated by recording multiple force curves to measure the frequencies of adhesive events and adhesion forces between two *E. coli* cells (Table 2). For homotypic interactions, Ag43<sup>C2-PT</sup> expressing cells exhibited the highest percentage of adhesive events, whereas bacterial cells expressing Ag43<sup>C3-PT</sup> exhibited the lowest (approximately 72 and 28%, respectively). For homotypic associations of Ag43<sup>C1-PT</sup> cells, the distribution of the adhesion forces show the highest most probable rupture force with 113 pN (Table 2 and Supplementary Fig. S6). Interestingly, adhesion forces for Ag43<sup>C2-PT</sup> homotypic interactions showed a bimodal distribution with two maxima at 91 and 223 pN, respectively, which suggest that two different interaction conformations can exist for this protein (Table 2 and Supplementary Fig. S6). Force curves show that multiple adhesion events took place during the self-interaction of Ag43<sup>C1-PT</sup> or Ag43<sup>C2-PT</sup> expressing cells, whereas single interaction events were mostly recorded for Ag43<sup>C3-PT</sup> or Ag43<sup>C4-PT</sup> self-interaction (Fig. 8A). For cells expressing Ag43<sup>C1-PT</sup> or Ag43<sup>C2-PT</sup> but not Ag43<sup>C3-PT</sup> or Ag43<sup>C4-PT</sup>, statistical analysis of the number of rupture events detected per force curves evidenced a maximum of 10 rupture events, indicating sequential unfolding of these homotypic interactions (Fig. 8B). Using the WLC (worm-like chain) model (allowing to measure the contour length of the
<figure><img src="image_7.png" /><figcaption>**Figure 6.** Fluorescent microscopy analysis of heterotypic interactions of *E. coli* cells expressing different variants of Ag43. (A) Co-culture of SRUD-Green Ag43<sup>C1-PT</sup> (C1) and SRUD-Red Ag43<sup>C2-PT</sup> (C2). (B) Co-culture of SRUD-Green Ag43<sup>C2-PT</sup> and SRUD-Red Ag43<sup>C4-PT</sup> (C4). (C) Co-culture of SRUD-Green Ag43<sup>C1-PT</sup> and SRUD-Red Ag43<sup>C3-PT</sup> (C3). (D) Co-culture of SRUD-Green Ag43<sup>C2-PT</sup> and SRUD-Red Ag43<sup>C3-PT</sup>. (E) Co-culture of SRUD-Green Ag43<sup>C1-PT</sup> and SRUD-Red Ag43<sup>C4-PT</sup>. (F) Co-culture of SRUD-Green Ag43<sup>C3-PT</sup> and SRUD-Red Ag43<sup>C4-PT</sup>. The NA, C1, C2, C3 and C4 labels match the colour of the fluorescent bacterial strains.</figcaption></figure>

unfolded molecule, which is defined as the length of the fully extended protein), the distribution of adhesion forces as a function of the contour length of unfolded Ag43 displayed different type of interactions (Fig. 8C). For Ag43<sup>C2-PT</sup>, Ag43<sup>C3-PT</sup> and Ag43<sup>C4-PT</sup>, the presence of adhesive events associated with unfolding at contour length of >400 nm and rupture forces (50–200 pN) suggests the presence of specific interactions between two cells via multiple electrostatic and/or hydrogen bonds. At the same time, adhesive events with contour length of <100 nm were observed for all Ag43 variants and suggests non-specific hydrophobic interactions. Taking together, these observations could explain the highest frequency of adhesive events measured for Ag43<sup>C2-PT</sup>.

As mixed cell aggregates were observed between cells expressing Ag43<sup>C4-PT</sup> and (i) Ag43<sup>C1-PT</sup>, as well as at lesser extent with (ii) Ag43<sup>C3-PT</sup>, but not with (iii) Ag43<sup>C2-PT</sup>, these heterotypic interactions were further investigated. A high frequency of adhesive events reaching approximately 74% was recorded between cells expressing Ag43<sup>C4-PT</sup> and Ag43<sup>C1-PT</sup> (Table 2). Conversely, this frequency remained moderate between cells expressing Ag43<sup>C4-PT</sup> and Ag43<sup>C3-PT</sup> (approximately 32%) and as expected extremely poor between Ag43<sup>C4-PT</sup> and Ag43<sup>C2-PT</sup> (approximately 4%), the latter indicative of an absence of significant interactions. Interestingly, the distribution of adhesion forces in heterotypic interactions showed two maxima corresponding to a bimodal distribution indicating the presence of two Ag43 binding conformations (Supplementary Fig. S7). For heterotypic interactions
<figure><img src="image_8.png" /><figcaption>**Figure 7.** Characterisation of heterotypic interactions of Ag43 variants by flow cytometry. Flow cytometry analyses were performed by mixing SRUD-Green and SRUD-Blue expressing respectively Ag43<sup>C1-PT</sup> (C1) and Ag43<sup>C2-PT</sup> (C2), Ag43<sup>C2-PT</sup> and Ag43<sup>C4-PT</sup> (C4), Ag43<sup>C1-PT</sup> and Ag43<sup>C3-PT</sup> (C3), Ag43<sup>C2-PT</sup> and Ag43<sup>C3-PT</sup>, Ag43<sup>C1-PT</sup> and Ag43<sup>C4-PT</sup>, or Ag43<sup>C3-PT</sup> and Ag43<sup>C4-PT</sup>. (A) Based on the SSC-A over FSC-A representation (pictured in the insert), the fluorescent FITC-H over PB450-H representation focused on the gate corresponding to the cell aggregates as previously defined (Supplementary Fig. S3). Graphs correspond to one representative experiment. (B) The percentage of mixed cell aggregates was calculated from the events detected in orange cloud over the total number of events corresponding to bacterial cell aggregates. (C) The concentration of single individual cells was calculated from the corresponding gate as previously defined (Supplementary Fig. S3). (D) The size of the aggregates was recovered from the gate for cell aggregates as previously defined (Supplementary Fig. S3).</figcaption></figure>

between cells expressing Ag43<sup>C4-PT</sup> and Ag43<sup>C1-PT</sup>, multiple rupture events were observed with a maximum of 21 rupture events (Fig. 9A,B). Additionally, Ag43<sup>C1-PT</sup> and Ag43<sup>C4-PT</sup> heterotypic association seemed to also exhibit both non-specific hydrophobic and several specific interactions via multiple electrostatic and/or hydrogen bonds (Fig. 9C). In summary, whenever the occurrence of homotypic or heterotypic interactions, AFM analyses emphasised the important role of both specific and nonspecific interactions in the Ag43-mediated bacterial autoaggregation.

Discussion

Previous phylogenetic analysis had simply divided Ag43 homologues into two distinct clusters, namely SF-I and SF-II<sup>24</sup>. Our further phylogenetic network analysis on the functional Ag43 passenger domains responsible for bacterial autoaggregation, has now revealed for the first time the existence of four distinct Ag43 phylogenetic classes, C1 to C4. Despite these divisions, structural modelling predicted that all Ag43 passengers form a similar L-shaped turn, the passenger structure themselves being subdivided into SL, EJ and BL subdomains. Ag43 passengers consistently displayed overall right-handed parallel β-helix fold with the exception of the two hairpins
<table><caption>**Table 2.** Percentage of adhesive events and values of the rupture force measured by single-cell force spectroscopy from bacteria expressing different Ag43 variants. <sup>A</sup>E. Coli SRUD was used to express the prototypical Ag43 variants from the different classes, namely Ag43<sup>C1-PT</sup> (C1), Ag43<sup>C2-PT</sup> (C2), Ag43<sup>C3-PT</sup> (C3) and Ag43<sup>C4-PT</sup> (C4). <sup>B</sup>Data correspond to mean value ± standard deviation and were extracted from AFM force-distance curves analysed with WLC model as described in the Methods section (Figs S4 and S5). <sup>C</sup>Data correspond to the most probable rupture force (most frequent and maximum value) ± the range of the distribution of the adhesion forces and were extracted from AFM force-distance curves analysed with WLC model as described in the Methods section (Figs S4 and S5). From bimodal distributions, two maxima could be calculated.</caption><thead><tr><th>Binding Couple<sup>A</sup></th><th>Adhesive Events (%)<sup>B</sup></th><th>Rupture Force (PN)<sup>C</sup></th></tr></thead><tbody><tr><td colspan="3">Homotypic Interactions</td></tr><tr><td>C1/C1</td><td>44.0 ± 9.0</td><td>113 ± 13</td></tr><tr><td>C2/C2</td><td>71.9 ± 14.2</td><td>91 ± 30; 223 ± 175</td></tr><tr><td>cC3/C3</td><td>27.8 ± 1.5</td><td>78 ± 16</td></tr><tr><td>C4/C4</td><td>37.3 ± 10.3</td><td>84 ± 108</td></tr><tr><td colspan="3">Heterotypic Interactions</td></tr><tr><td>C4/C2</td><td>4.2 ± 2.6</td><td>149 ± 90</td></tr><tr><td>C4/c1</td><td>73.8 ± 6.1</td><td>67 ± 30; 350 ± 200</td></tr><tr><td>C4/c3</td><td>31.9 ± 6.2</td><td>68 ± 40; 315 ± 100</td></tr></tbody></table>

(each formed of two antiparallel β-sheets) delimiting the EJ subdomain. While there was evidence of structural variation reported in the literature<sup>30,31</sup>, the present investigation evidenced that each SL, EJ and BL subdomain belong to two distinct phylogenetic types, whereby combinations of SL with EJ-BL subdomains divide the four different classes of Ag43 passengers. This new finding represents a major breakthrough in the field as it provides a rational framework for directed functional investigations of Ag43.

With this new understanding of the diversity of Ag43, the autoaggregation phenotype was reinvestigated based on this updated Ag43<sup>C1</sup>, Ag43<sup>C2</sup>, Ag43<sup>C3</sup> and Ag43<sup>C4</sup> classification. Following bacterial cell expression of several Ag43 from each of these classes, the most striking difference was based in the significant delay of the sedimentation kinetics for Ag43 proteins belonging to the class 3, compared to the Ag43<sup>C1</sup>, Ag43<sup>C2</sup> and Ag43<sup>C4</sup>. While phenotypical analyses using different Ag43 have been reported in the scientific literature<sup>30,31</sup>, no ration-ale had been previously described. Further, some contradictory results have also been conveyed<sup>24</sup> that could be attributed to differences in critical parameters in the experimental set-up, such as the culture media (e.g. complex over chemically defined) or the bacterial genetic background (e.g. the use of E. coli laboratory strain for Ag43 expression still containing *agn43* and/or *fim*). Here, using the nonaggregative E. coli SRUD strain and the same growth medium, the Ag43a from UPEC CFT073 (UPI00000E4A66) appeared to exhibit a slow sed-imentation as originally reported<sup>31</sup>. Together with Ag43b from UPEC CFT073 (UPI00000E49D2) and other Ag43<sup>C3</sup>, their sedimentation kinetics were significantly slower than Ag43<sup>C1</sup>, Ag43<sup>C2</sup> or Ag43<sup>C4</sup>. Further, these aggregation experiments also revealed differences within C3 with Ag43a showing the strongest sedimentations properties in this class compared to other proteins like Ag43b, which showed weaker cell settling properties. Therefore, the latter Ag43 homologue better represents the phenotypical properties of this class and could be considered a prototypical Ag43<sup>C3</sup>. Most of the current understanding on Ag43 stems from the work done on the prototypical Ag43<sup>C1-UPI00003B25EE</sup> from E. coli K12 MG1655, which has a SL1-EJ1-BL1 passenger. To date, only one article has been reported characterising Ag43<sup>C2-UPI00000D07E5</sup> from EHEC O157:H7 EDL933, which presents a SL2-EJ2-BL2 passenger<sup>26</sup>. Our work is the first to report the characterisation of Ag43 homologues from C4, such as Ag43<sup>C4-UPI0001E6D993</sup> from enterotoxigenic E. coli O78:H11 H10407, which exhibits a SL1-EJ2-BL2 pas-senger<sup>32</sup>. We show that bacterial cells expressing different Ag43 variants were able to induce autoaggregation through homotypic interactions as supported by microscopic examination and flow cytometry analyses (Fig. 10). However, differences were observed among the four Ag43 classes, with Ag43<sup>C3-PT</sup> showing slower sedimentation kinetics, looser autoaggregation pellets and significantly smaller aggregates associated with lower rupture force revealed by AFM analyses. With regard to single-molecule force spectroscopy analysis, it is worth noting that for Ag43<sup>C1-PT</sup> similar force curves and rupture forces of 100–200 pN were previously recorded<sup>33</sup>.

Considering that heterologous associations have been reported between very different members of the SAATs family, namely Ag43, AIDA-I and TibA<sup>17,19</sup>, the paucity of heterotypic interactions between Ag43 variants evi-denced in this work was unexpected. While the recognition and involvement of binding interfaces in both homo-typic and heterotypic interactions were suggested in the case of co-interacting SAATs, the restricted number of binding pairs respective to the different Ag43 variants suggest the involvement of different binding interfaces. Structural studies of homotypic interactions in Ag43<sup>C3-UPI00000E4A66</sup> revealed a Velcro-like handshake mechanism via trans-association of the SL subdomains stabilised by multiple hydrogen bonds and electrostatic interactions<sup>25</sup>. Ag43<sup>C1-PT</sup> and Ag43<sup>C4-PT</sup> share the SL1 subdomain and we could observe some heterotypic interactions between cells expressing those Ag43 homologues. Nonetheless, the level and force of this heterotypic association are signif-icantly weaker than the corresponding homotypic interactions. Differences at amino acid level between Ag43<sup>C1-PT</sup> and Ag43<sup>C4-PT</sup> proteins probably result in a reduced number of interactions stabilising the hetero-dimers com-pared to the homodimers. Additionally, the EJ-BL subdomains could also play an indirect role in the Ag43-Ag43 interaction. Indeed, previous work showed that modification of the EJ subdomain altered the angle between the SL and BL subdomains and the overall L-shaped structure, completely abolishing the ability of cells expressing
<figure><img src="image_9.png" /><figcaption>**Figure 8.** Retraction force curve and analysis of the rupture events during homotypic interactions of *E. coli* cells expressing different Ag43 variants. (A) Retraction force curves. (B) Rupture events measured per force curve. (C) Distribution of adhesion force measured between two *E. coli* cells expressing different Ag43 as function of the contour length for each force-distance curves. Analyses were performed using SRUD expressing the Ag43 from the different classes, namely Ag43<sup>C1-PT</sup> (C1), Ag43<sup>C2-PT</sup> (C2), Ag43<sup>C3-PT</sup> (C3) or Ag43<sup>C4-PT</sup> (C4).</figcaption></figure>

Ag43 to form aggregates²⁵. The absence of aggregates between cells expressing Ag43<sup>C2-PT</sup> and Ag43<sup>C3-PT</sup>, despite both having a SL2, suggests this subdomain is less tolerant towards heterotypic associations and may be due to the binding interface regions being more selective in forming interactions. Similarly, no heterotypic interactions could be detected for combinations of Ag43 exhibiting different SL subdomain, namely SL1/SL2, except in the case of cells expressing Ag43<sup>C3-PT</sup> and Ag43<sup>C4-PT</sup>. Considering the only two observed heterotypic interactions involved Ag43<sup>C4-PT</sup>, it is tempting to hypothesise that the amino acid composition of SL1 and the EJ2-BL2 architecture in the passenger are more permissive to heterotypic interactions. Diverse EJ-BL subdomains may induce different bending angles to the L-shape Ag43 structure and modulate the ability to form heterologous associations. However, the molecular modelling approach used here to model the structure of different Ag43 homologues was not able to detect any differences as the angle between SL and BL biased the bending towards the Ag43<sup>C3-UPI100000E49D2</sup> template. Solving the experimental structures of Ag43 belonging to classes 1, 2 and 4, as well as other SAATs is critical to further inform the role of the different passenger subdomains in inducing autoaggregation and potentially predict heterologous associations among different SAATs.

Autoaggregation is a widespread phenomenon in bacteria but its ecophysiological importance is not fully elucidated. Several direct and side effects of autoaggregation have been reported in the literature, e.g. mechanism of adaptation against different stresses such as host defences as well as promoting cell invasion and colonisation¹⁸. The differential autoaggregation as well as the paucity of heterotypic interactions suggests that the Ag43 variants could play a role in structuring the *E. coli* population in an ecosystem by rallying or excluding some bacterial strains and thus stratifying the bacterial community. Considering that up to three *agn43* alleles can be present in one bacterial genome²⁴, further genomic investigations would be required to characterise the prevalence, distribution and diversity of Ag43, such as the combinations of classes and domain shuffling, among *E. coli* and within the Bacteria kingdom, as well the correlations with some features of the bacterial strains, such as the origin, the level of virulence or serotype. The reason why *E. coli* can carry different *agn43* alleles in one genome and express so many different variants of Ag43 is still unclear²⁴ but their implication in social behaviour especially through selective altruism following a greenbeard effect is an attractive concept³⁴⁻³⁶. Considering *agn43* can identify the presence of copies of themselves in other cells through Ag43-Ag43 interactions, this might induce some kind of nepotism behaviour toward those individual cells. The diversity of Ag43 variants and multiplicity of alleles could then favour differential and beneficial interactions between individual cells in a bacterial population.

These different aspects would undoubtedly require further investigations to understand how it could contribute to colonisation process and specific adaptation of *E. coli* within an ecosystem, both considering commensal
<figure><img src="image_10.png" /><figcaption>**Figure 10.** Summary of the different combinations of homotypic and heterotypic interactions of *E. coli* cells expressing Ag43 from the different classes in relation with the model structure of the passengers. Size of the black arrows schematically represent the force of the interactions. Grey lines indicate the tested combinations where no interaction could be experimentally evidenced. All structures of the Ag43 passengers were modelled using Phyre<sup>54,55</sup>.</figcaption></figure>

To obtain fluorescent bacterial strains, genes encoding for mRuby2, mTagBFP2 and mEmerald fluorescent proteins (GeneBank: AFR60232.1, AIQ82697.1 and AIA24532.1, respectively) were synthesised by GeneArt including a constitutive promoter (iGEM: BBa-J23119) and ribosome binding site (RBS) from pBAD/Myc-HisA, prior to cloning into the pACYC184 *Ava*I and *Bsu*36I/*Tth*111I restriction sites to generate pVaAgRed, pVaAgBlue and pVaAgGreen respectively, whose p15a replication origin belongs to a different incompatibility group from the pBR322 replication origin in pBAD-vector series. Plasmid constructs were transformed by electroporation in *E. coli* SRUD with or without a pVaAg-vector series and selection was carried on LB agar supplemented with chloramphenicol and when required with the additional supplementation of ampicillin.

**Immunodetection and subcellular localisation of Ag43 variants.** Bacterial cells were grown with arabinose till late exponential phase as described above. Cells were harvested by centrifugation (7 000 g, 5 min, 4°C), washed and resuspended in PBS (0.01 M, 1 mM PMSF, pH8). Cells were disrupted with glass microbeads using FastPrep (20 s twice at 6 m.s⁻¹). Part of lysed whole cells was reserved and the other part subjected to cell fractionation using a modified sarkosyl (sodium lauroyl sarcosinate) method<sup>39,40</sup>. Briefly, unbroken cells and cellular debris were removed by centrifugation (12 000 g, 10 min, 4°C) before the supernatant was ultracentrifuged (100 000 g, 40 min, 4°C). The supernatant was retained as the soluble fraction (cytoplasm and periplasm), and the pellet was resuspended in 2% sarkosyl (20 min, room temperature) followed by another ultracentrifugation. The pellet corresponded to the sarcosyl-insoluble fraction gathering the OM. Cold acetone was added to the supernatant (overnight at -20°C), centrifuged and the pellet corresponding to the sarcosyl-soluble fraction gathering the IM. Protein concentrations were determined using the Bradford method<sup>41</sup>.

Proteins from the different fractions were separated by SDS-PAGE. Cell fractions were diluted in Laemmli buffer and heated up at 95 °C for 5 min. From prior testing with increasing amounts of whole cell fractions, 1 µg of total proteins was determined as a non-saturating concentration and equal amounts of each fraction based on the starting material were loaded per lane using Mini-ProTEAN Precast TGX Gels 4–15% (Bio-Rad). Precision Plus
Protein Dual Color Standards (Bio-Rad) was used to enable molecular mass estimation as well as transfer control. For Western blotting, proteins were dry transferred from SDS-PAGE on iBlot Transfer Stacks using iBlot 2 Dry Blotting System (Invitrogen). After saturating unoccupied protein-binding sites with blocking agents, membranes were probed with specific anti-Ag43 primary antibodies (dilution 1/5000 in PBS with 0.1% tween and 0.1% milk, 1 h, room temperature). Rabbit Ag43 antibodies were obtained as previously described<sup>42,43</sup> and are crossed reactive between Ag43 variants<sup>44</sup>. Goat anti-Rabbit IgG antibodies conjugated with horseradish peroxidase was used as secondary antibodies (dilution of 1/10 000, 1 h, room temperature) and TMB (3,3',5,5'-tetramethylbenzidine) was used to detect the enzymatic activity. Enzymatic reaction was stopped by rinsing the membranes with milliQ water. Western blots were imaged under white light conversion screen using E-BOX (VILBER).

**Phylogenetic analysis of protein sequences.** The Ag43 protein sequences were retrieved from the literature<sup>24,28</sup> and the corresponding passengers recovered by removing the SP, AC1 and translocator regions. For identification of the SPs in Ag43 proteins, the sequences submitted to SignalP v4.1<sup>45</sup> were deleted from the ESPR (InterPro: IPR024973), which was first identified using InterProScan v5.22<sup>46</sup>. The AC1 domains were identified based on a structural alignment as previously described<sup>3</sup>. The translocators were identified using the twin HMM (Hidden Markov Model) strategy for AT detection<sup>47</sup>. Finally, truncated and redundant Ag43 passenger sequences were removed, resulting in a dataset of 106 distinct Ag43 passengers (Supplementary Table S1). Then, the amino acid sequences were aligned with T-Coffee in the Expresso mode for automatic incorporation of structural information<sup>48,49</sup> using the solved structure of Ag43 (PDB: 4KH3) as template. The relatedness among protein sequences was studied using a Neighbor-Net phylogenetic network approach with SplitsTree v4.14.8 using the Hamming uncorrected-P distance<sup>50</sup>. The most robust splits of the network were identified by bootstrap using a 95% threshold over 1,000 pseudo-replicates. A similar phylogenetic analysis was repeated for the different passenger subdomains. Pairwise sequence similarities and identities were obtained using Needleman-Wunsch alignment of two sequences with BLAST<sup>51,52</sup> (Supplementary Table S1).

**Identification of conserved motifs and structural modelling.** Conserved motifs were searched using InterProScan v5.22<sup>46</sup> and interrogating InterPro (IPR) v61.0 database<sup>53</sup>. Three-dimensional structures were modelled with Phyre v2.0<sup>54,55</sup> as well as I-TASSER v5.1<sup>56,57</sup>. The types of β-solenoid motif present in the β-helices were predicted using PfScan<sup>29</sup>. In β-helices, rungs were defined as adjacent parallel β-strands. For selected variants as presented in Fig. 2, Ag43 homolog models were checked and modified in the program COOT<sup>58</sup> using its Ramachandran, geometry and rotamer analysis validation software. The quality of the models were then assessed using ProCheck<sup>39</sup> with its Ramachandran errors reported as percentages of generously allowed and disallowed amino acids (Supplementary Table S3). Protein structures were visualised using PyMOL v1.7.4<sup>60,61</sup>.

**Autoaggregation assay.** Autoaggregation assays were performed as previously described<sup>62</sup>. Briefly, overnight precultures were diluted 1:100 and grown in LB at 37°C under orbital shaking. Arabinose was added in early exponential phase (OD<sub>600 nm</sub> = 0.25) to induce the Ag43 expression from pBAD-vector series (Supplementary Table S2). Cell suspensions were placed in conical tubes at the same OD<sub>600 nm</sub>, incubated statically and sampled every hour from the top of the tube to measure the OD<sub>600 nm</sub>. Observations in phase-contrast microscopy were made at the end of the experiments from sediment at the bottom of the tube to confirm the bacterial cell autoaggregation.

**Image acquisition and analysis of phase contrast and epifluorescence microscopic observations.** Microscopic observations of bacterial cells were performed in phase contrast and phase contrast epifluorescence reflected light. Image acquisition was performed using an Axioplan 2 microscope coupled to an AxioCam MRm camera driven by the AxioVision 4.8 software (Carl Zeiss). For phase contrast microscopy, a 40 × /0.75 objective (Plan-NEOFLUAR) was used and the fluorescence light source was a mercury lamp (HBO100W/2, AttoArc, Carl Zeiss). For fluorescence image acquisition a SYTO9-AO-FITC cube was fitted with a 450–490 nm exciter, 515–565 nm emitter, and 510 nm dichroic mirror as well as a Propidium iodide-Rhodamine (Con A) cube containing 546/12 nm exciter, 590 nm emitter, and 580 nm dichroic mirror (Carl Zeiss).

**Quantification of bacterial cell interactions by flow cytometry.** To investigate interactions of bacterial cells expressing Ag43 variants, *E. coli* SRUD-Green (carrying pVaAgGreen) or *E. coli* SRUD-Blue (carrying pVagAgBlue) were transformed with pBAD/Myc-HisA vectors encoding the different Ag43 variants or empty vector (negative control with non-aggregative *E. coli* SRUD bacterial cells) (Supplementary Table S2). Overnight *E. coli* SRUD-Green and *E. coli* SRUD-Blue precultures were diluted 1:100 and grown as described above in LB until early exponential phase. At this point, samples from these two different bacterial cultures were adjusted at an OD<sub>600 nm</sub> of 0.25, mixed together at a 50/50 (v/v) ratio with arabinose and grown at 37 °C under shaking until an OD<sub>600 nm</sub> of 0.5. Conical tubes were then incubated statically for 6 h and bacterial cells fixed with paraformaldehyde (2% final concentration) for 15 min. Cells were collected in the bottom of each tube and diluted in PBS (phosphate-buffered saline).

Samples were filtered on a 40-µm-nylon cell strainer prior to flow cytometry analysis (CytoFLEX, Beckman Coulter). FSC (forward scatter), SSC (side scatter), FITC (fluorescein isothiocyanate; 488 nm exciter, 525/540 nm emitter) and PB450 (Pacific Blue 450 nm; 405 nm exciter, 450/445 nm emitter) axes were set to log display. Areas for FSC and SSC (FSC-A and SSC-A, respectively) were used to discriminate aggregates from single individual cells, whereas heights for FITC and PB450 (FITC-H and PB450-H, respectively) were used to quantify co-interactions within aggregates (Supplementary Fig. S3). At least 3000 events/sec were recorded for each sample and data were analysed with CytExpert software v2.2 (Beckman Coulter). Results were expressed as percentages of mixed aggregates over the total aggregates.
**Atomic force microscopy analyses.** For cell probe preparation, bacterial cells were handled as for the autoaggregation assay prior to dilution in PBS buffer to an OD$_{600}$$_{m}$ of 0.2. One mL of bacterial suspension was directly deposited onto glass slide pre-coated with polyethylenimin (PEI) as previously described$^{63}$. Briefly, cleaned glass slides are immersed in 0.2% PEI solution for 1 h, then extensively rinsed with Milli-Q water and finally dried under nitrogen) for 30 min and then the surface was rinsed twice with PBS buffer to get rid of un-adhered free bacteria and remaining residues from growth medium.

A Nanowizard III AFM (JPK-Bruker) was used to attach to atipless silicon nitride cantilevers (NP-010, Bruker) 6 µm borosilicate microspheres using UV curing resin (Norland optical adhesive 63 Glue, Norland Edmund Optics). Sphere attached cantilevers were cured under UV light (λ = 400 nm) for 3 min. The cantilever was then immersed for 1 h in a 0.2% PEI solution and washed. Functionalised cantilevers were calibrated before analysis by the thermal noise method$^{64}$ and spring constants were found to be 0.048 ± 0.011 N/m. Colloidal probe$^{65}$ was brought into contact with an isolated cell for 3 min, and the obtained cell probe was then moved to another cell to perform single cell-force spectroscopy measurements. Suitable parameters for cell attachment are a contact force of 1 nN, a contact time of 2 s and an approach and retract rate of 1 µms$^{-1}$.

Single cell-force spectroscopy experiments were carried out with a Nanowizard III atomic force microscope (JPK-Bruker). Measurements were performed at room temperature in PBS buffer at pH 7.2. Force-distance curves were recorded on different cells, using a maximum applied force of 1 nN and a constant approach and retraction speed of 1 µm.s$^{-1}$. The time of contact was set at 2 s to enhance specific interactions as described previously$^{66,67}$. The type of interactions, i.e. specific or non-specific, could be deduced from the patterns of the force-distance curves$^{68–70}$.

The obtained force versus distance curves were analysed with the JPK Data Processing software on the basis of the worm-like chain (WLC) model, which is suitable to describe the extension of protein. The extension x of the molecule versus the pulling force (or adhesion force) is given by the following equation:

$$F(x) = \frac{k_B T}{l_p} \left[ \frac{1}{4} \left( 1 - \frac{x}{L_c} \right)^{-2} + \frac{x}{L_c} - \frac{1}{4} \right]$$

where k$_{B}$ is the Boltzmann constant, T the absolute temperature, l$_{p}$ is the persistent length and L$_{c}$ is the contour length. The persistent length gives information about the degree of structural rigidity of the polymer chain and is defined as the longest segment below which the chain can be considered as a rigid rod, the extension of the polymer until the point at which the force necessary to extend further and in our case to break the interaction with the functionalised AFM tip corresponds to the contour length.

Adhesion peaks were generated using Origin 8 software. Statistical differences of AFM data were assessed by analysis of the variance (one-way ANOVA test) using Origin 8. To check the specificity of the measured adhesion forces and avoid artifacts that can be associated with cell probe preparation, control experiments were performed. Use of PEI-coated probes without bacterial cell interacting on PEI-coated glass slide led to a dramatic decrease of the adhesion frequency with only 3% of the force curves presenting an adhesion peak.

**Statistical analysis.** Statistical analyses were performed using GraphPad Prism v6.01. At least five biological replicates were carried out for all statistical analyses. Error bars on the figures represent the 95% confidence interval for the mean of the independent experiments. Data were statistically analysed following ANOVA one way or two way test, followed by the post-hoc Tukey test, with differences considered significant (p < 0.05, * or lowercase letter), very significant (p < 0.01, ** or bold lowercase letter), highly significant (p < 0.001, *** or capital letter) or very highly significant (p < 0.0001, **** or bold capital letter).

## Data Availability

All data generated or analysed during this study are included in this published article and its Supplementary Material files. Other data are available from the corresponding author upon reasonable request.

## References

1. Desvaux, M., Parham, N. J., Scott-Tucker, A. & Henderson, I. R. The general secretory pathway: a general misnomer? *Trends in microbiology* **12**, 306–309, https://doi.org/10.1016/j.tim.2004.05.002 (2004).

2. Desvaux, M., Hébraud, M., Talon, R. & Henderson, I. R. Secretion and subcellular localizations of bacterial proteins: a semantic awareness issue. *Trends in microbiology* **17**, 139–145, https://doi.org/10.1016/j.tim.2009.01.004 (2009).

3. Rojas-Lopez, M. *et al*. Identification of the Autochaperone Domain in the Type Va Secretion System (T5aSS): Prevalent Feature of Autotransporters with a β-Helical Passenger. *Frontiers in microbiology* **8**, 2607, https://doi.org/10.3389/fmicb.2017.02607 (2017).

4. Yen, M. R. *et al*. Protein-translocating outer membrane porins of Gram-negative bacteria. *Biochimica et biophysica acta* **1562**, 6–31, https://doi.org/10.1016/S0005-2736(02)00359-0 (2002).

5. Abby, S. S. *et al*. Identification of protein secretion systems in bacterial genomes. *Scientific reports* **6**, 23080, https://doi.org/10.1038/srep23080 (2016).

6. Henderson, I. R. & Lam, A. C. Polymorphic proteins of *Chlamydia* spp.-autotransporters beyond the *Proteobacteria*. *Trends in microbiology* **9**, 573–578, S0966-842X(01)02234-X (2001).

7. Desvaux, M., Khan, A., Beatson, S. A., Scott-Tucker, A. & Henderson, I. R. Protein secretion systems in *Fusobacterium nucleatum*: genomic identification of Type 4 piliation and complete Type V pathways brings new insight into mechanisms of pathogenesis. *Biochim Biophys Acta* **1713**, 92–112, https://doi.org/10.1016/j.bbamem.2005.05.002 (2005).

8. Leyton, D. L., Rossiter, A. E. & Henderson, I. R. From self sufficiency to dependence: mechanisms and factors important for autotransporter biogenesis. *Nat Rev Microbiol* **10**, 213–225, https://doi.org/10.1038/nrmicro2733 (2012).

9. Selkrig, J. *et al*. Discovery of an archetypal protein transport system in bacterial outer membranes. *Nat Struct Mol Biol* **19**(506–510), S501, https://doi.org/10.1038/nsmb.2261 (2012).
10. Henderson, I. R., Navarro-Garcia, F., Desvaux, M., Fernandez, R. C. & Ala'Aldeen, D. Type V protein secretion pathway: the autotransporter story. *Microbiology and molecular biology reviews* **MMBR 68**, 692–744, https://doi.org/10.1128/MMBR.68.4.692-744.2004 (2004).

11. Leo, J. C., Grin, I. & Linke, D. Type V secretion: mechanism(s) of autotransport through the bacterial outer membrane. *Philosophical transactions of the Royal Society of London. Series B, Biological sciences* **367**, 1088–1101, https://doi.org/10.1098/rstb.2011.0208 (2012).

12. Desvaux, M., Parham, N. J. & Henderson, I. R. Le système de sécrétion de type V chez les bactéries Gram-négatives (Type V secretion system in Gram-negative bacteria). *Biofutur* **237**, 34–37 (2003).

13. Desvaux, M., Parham, N. J. & Henderson, I. R. Type V protein secretion: simplicity gone awry? *Curr Issues Mol Biol* **6**, 111–124 (2004).

14. Tsai, J. C. et al. The bacterial intimins and invasins: a large and novel family of secreted proteins. *PloS one* **5**, e14403, https://doi.org/10.1371/journal.pone.0014403 (2010).

15. Desvaux, M., Parham, N. J. & Henderson, I. R. The autotransporter secretion system. *Research in microbiology* **155**, 53–60, https://doi.org/10.1016/j.resmic.2003.10.002 (2004).

16. Drobnak, I. et al. Of linkers and autochaperones: an unambiguous nomenclature to identify common and uncommon themes for autotransporter secretion. *Molecular microbiology* **95**, 1–16, https://doi.org/10.1111/mmi.12838 (2015).

17. Klemm, P., Vejborg, R. M. & Sherlock, O. Self-associating autotransporters, SAATs: functional and structural similarities. *Int J Med Microbiol* **296**, 187–195, https://doi.org/10.1016/j.ijmm.2005.10.002 (2006).

18. Trunk, T., Khalil, H. S. & Leo, J. C. Bacterial autoaggregation. *Aims Microbiol* **4**, 140–164, https://doi.org/10.3934/microbiol.2018.1.140 (2018).

19. Sherlock, O., Schembri, M. A., Reisner, A. & Klemm, P. Novel roles for the AIDA adhesin from diarrheagenic *Escherichia coli*: cell aggregation and biofilm formation. *Journal of bacteriology* **186**, 8058–8065, https://doi.org/10.1128/JB.186.23.8058-8065.2004 (2004).

20. Kjaergaard, K., Schembri, M. A., Ramos, C., Molin, S. & Klemm, P. Antigen 43 facilitates formation of multispecies biofilms. *Environ Microbiol* **2**, 695–702, https://doi.org/10.1046/j.1462-2920.2000.00152.x (2000).

21. Desvaux, M. et al. A conserved extended signal peptide region directs posttranslational protein translocation via a novel mechanism. *Microbiology* **153**, 59–70, https://doi.org/10.1099/mic.0.29091-0 (2007).

22. Jong, W. S. et al. YidC is involved in the biogenesis of the secreted autotransporter hemoglobin protease. *The Journal of biological chemistry* **285**, 39682–39690, https://doi.org/10.1074/jbc.M110.167650 (2010).

23. Oliver, D. C., Huang, G., Nodel, E., Pleasance, S. & Fernandez, R. C. A conserved region within the *Bordetella pertussis* autotransporter BrkA is necessary for folding of its passenger domain. *Molecular microbiology* **47**, 1367–1383, https://doi.org/10.1046/j.1365-2958.2003.03377.x (2003).

24. van der Woude, M. W. & Henderson, I. R. Regulation and function of Ag43 (flu). *Annual review of microbiology* **62**, 153–169, https://doi.org/10.1146/annurev.micro.62.081307.162938 (2008).

25. Heras, B. et al. The antigen 43 structure reveals a molecular Velcro-like mechanism of autotransporter-mediated bacterial clumping. *Proceedings of the National Academy of Sciences of the United States of America* **111**, 457–462, https://doi.org/10.1073/pnas.1311592111 (2014).

26. Torres, A. G. et al. Characterization of Cah, a calcium-binding and heat-extractable autotransporter protein of enterohaemorrhagic *Escherichia coli*. *Molecular microbiology* **45**, 951–966, https://doi.org/10.1046/j.1365-2958.2002.03094.x (2002).

27. Ulett, G. C. et al. Functional analysis of antigen 43 in uropathogenic *Escherichia coli* reveals a role in long-term persistence in the urinary tract. *Infection and immunity* **75**, 3233–3244, https://doi.org/10.1128/IAI.01952-06 (2007).

28. Vo, J. L. et al. Autotransporter Adhesins in *Escherichia coli* Pathogenesis. *Proteomics* **17**, https://doi.org/10.1002/pmic.201600431 (2017).

29. Kajava, A. V. & Steven, A. C. The turn of the screw: variations of the abundant beta-solenoid motif in passenger domains of Type V secretory proteins. *J Struct Biol* **155**, 306–315, https://doi.org/10.1016/j.jsb.2006.01.015 (2006).

30. Reidl, S., Lehmann, A., Schiller, R., Salam Khan, A. & Dobrindt, U. Impact of O-glycosylation on the molecular and cellular adhesion properties of the *Escherichia coli* autotransporter protein Ag43. *Int J Med Microbiol* **299**, 389–401, https://doi.org/10.1016/j.ijmm.2009.01.001 (2009).

31. Klemm, P., Hjerrild, L., Gjermansen, M. & Schembri, M. A. Structure-function analysis of the self-recognizing Antigen 43 autotransporter protein from *Escherichia coli*. *Molecular microbiology* **51**, 283–296, https://doi.org/10.1046/j.1365-2958.2003.03833.x (2004).

32. Crossman, L. C. et al. A commensal gone bad: complete genome sequence of the prototypical enterotoxigenic Escherichia coli strain H10407. *Journal of bacteriology* **192**, 5822–5831, https://doi.org/10.1128/JB.00710-10 (2010).

33. Jacquot, A. et al. The dynamics and pH-dependence of Ag43 adhesins' self-association probed by atomic force spectroscopy. *Nanoscale* **6**, 12665–12681, https://doi.org/10.1039/c4nr03312d (2014).

34. Gardner, A. & West, S. A. Greenbeards. *Evolution* **64**, 25–38, https://doi.org/10.1111/j.1558-5646.2009.00842.x (2010).

35. West, S. A. & Gardner, A. Altruism, spite, and greenbeards. *Science* **327**, 1341–1344, https://doi.org/10.1126/science.1178332 (2010).

36. West, S. A., Griffin, A. S. & Gardner, A. Social semantics: altruism, cooperation, mutualism, strong reciprocity and group selection. *J Evol Biol* **20**, 415–432, https://doi.org/10.1111/j.1420-9101.2006.01258.x (2007).

37. Bertani, G. Studies on lysogenesis. I. The mode of phage liberation by lysogenic *Escherichia coli*. *Journal of bacteriology* **62**, 293–300 (1951).

38. Bertani, G. Lysogeny at mid-twentieth century: P1, P2, and other experimental systems. *Journal of bacteriology* **186**, 595–600, https://doi.org/10.1128/jb.186.3.595-600.2004 (2004).

39. Filip, C., Fletcher, G., Wulff, J. L. & Earhart, C. F. Solubilization of the cytoplasmic membrane of *Escherichia coli* by the ionic detergent sodium-lauryl sarcosinate. *Journal of bacteriology* **115**, 717–722 (1973).

40. Huang, C. Z. et al. Systematic identification of the subproteome of *Escherichia coli* cell envelope reveals the interaction network of membrane proteins and membrane-associated peripheral proteins. *J Proteome Res* **5**, 3268–3276, https://doi.org/10.1021/pr060257h (2006).

41. Bradford, M. M. A rapid and sensitive method for the quantitation of microgram quantities of protein utilizing the principle of protein-dye binding. *Anal Biochem* **72**, 248–254, https://doi.org/10.1016/0003-2697(76)90527-3 (1976).

42. Owen, P., Caffrey, P. & Josefsson, L. G. Identification and partial characterization of a novel bipartite protein antigen associated with the outer membrane of *Escherichia coli*. *Journal of bacteriology* **169**, 3770–3777, https://doi.org/10.1128/jb.169.8.3770-3777.1987 (1987).

43. Henderson, I. R., Meehan, M. & Owen, P. Antigen 43, a phase-variable bipartite outer membrane protein, determines colony morphology and autoaggregation in *Escherichia coli* K-12. *FEMS microbiology letters* **149**, 115–120, S0378-1097(97)00065-7 (1997).

44. Roche, A., McFadden, J. & Owen, P. Antigen 43, the major phase-variable protein of the *Escherichia coli* outer membrane, can exist as a family of proteins encoded by multiple alleles. *Microbiology* **147**, 161–169, https://doi.org/10.1099/00221287-147-1-161 (2001).

45. Nielsen, H. Predicting Secretory Proteins with SignalP. *Methods Mol Biol* **1611**, 59–73, https://doi.org/10.1007/978-1-4939-7015-5_6 (2017).

46. Jones, P. et al. InterProScan 5: genome-scale protein function classification. *Bioinformatics* **30**, 1236–1240, https://doi.org/10.1093/bioinformatics/btu031 (2014).
47. Celik, N. et al. A bioinformatic strategy for the detection, classification and analysis of bacterial autotransporters. *PloS one* 7, e43245, https://doi.org/10.1371/journal.pone.0043245 (2012).

48. Armougom, F. et al. Expresso: automatic incorporation of structural information in multiple sequence alignments using 3D-Coffee. *Nucleic acids research* 34, W604–608, https://doi.org/10.1093/nar/gkl092 (2006).

49. Di Tommaso, P. et al. T-Coffee: a web server for the multiple sequence alignment of protein and RNA sequences using structural information and homology extension. *Nucleic acids research* 39, W13–17, https://doi.org/10.1093/nar/gkr245 (2011).

50. Huson, D. H. & Bryant, D. Application of phylogenetic networks in evolutionary studies. *Molecular biology and evolution* 23, 254–267, https://doi.org/10.1093/molbev/msj030 (2006).

51. Needleman, S. B. & Wunsch, C. D. A general method applicable to the search for similarities in the amino acid sequence of two proteins. *Journal of molecular biology* 48, 443–453, https://doi.org/10.1016/0022-2836(70)90057-4 (1970).

52. Tatusova, T. A. & Madden, T. L. Blast 2 sequences - a new tool for comparing protein and nucleotide sequences. *FEMS microbiology letters* 174, 247–250, https://doi.org/10.1111/j.1574-6968.1999.tb13575.x (1999).

53. Finn, R. D. et al. InterPro in 2017-beyond protein family and domain annotations. *Nucleic acids research* 45, D190–D199, https://doi.org/10.1093/nar/gkw1107 (2017).

54. Kelley, L. A., Mezulis, S., Yates, C. M., Wass, M. N. & Sternberg, M. J. The Phyre2 web portal for protein modeling, prediction and analysis. *Nat Protoc* 10, 845–858, https://doi.org/10.1038/nprot.2015.053 (2015).

55. Kelley, L. A. & Sternberg, M. J. Protein structure prediction on the Web: a case study using the Phyre server. *Nat Protoc* 4, 363–371, https://doi.org/10.1038/nprot.2009.2 (2009).

56. Zhang, Y. I-TASSER server for protein 3D structure prediction. *BMC Bioinformatics* 9, 40, https://doi.org/10.1186/1471-2105-9-40 (2008).

57. Yang, J. et al. The I-TASSER Suite: protein structure and function prediction. *Nat Methods* 12, 7–8, https://doi.org/10.1038/nnmeth.3213 (2015).

58. Emsley, P. & Cowtan, K. Coot: model-building tools for molecular graphics. *Acta Crystallogr D Biol Crystallogr* 60, 2126–2132, https://doi.org/10.1107/S0907444904019158 (2004).

59. Laskowski, R. A., MacArthur, M. W., Moss, D. S. & Thornton, J. M. PROCHECK - a program to check the stereochemical quality of protein structures. *J. Appl. Cryst.* 26, 283–291, https://doi.org/10.1107/S0021889892009944 (1993).

60. Yuan, S., Chan, H. C. S., Filipek, S. & Vogel, H. PyMOL and Inkscape Bridge the Data and the Data Visualization. *Structure* 24, 2041–2042, https://doi.org/10.1016/j.str.2016.11.012 (2016).

61. Rigsby, R. E. & Parker, A. B. Using the PyMOL application to reinforce visual understanding of protein structure. *Biochem Mol Biol Educ* 44, 433–437, https://doi.org/10.1002/bmb.20966 (2016).

62. Chagnot, C. et al. *In vitro* colonization of the muscle extracellular matrix components by *Escherichia coli* O157:H7: the influence of growth medium, temperature and pH on initial adhesion and induction of biofilm formation by collagens I and III. *PloS one* 8, e59386, https://doi.org/10.1371/journal.pone.0059386 (2013).

63. Vadillo-Rodriguez, V. et al. Comparison of atomic force microscopy interaction forces between bacteria and silicon nitride substrata for three commonly used immobilization methods. *Applied and environmental microbiology* 70, 5441–5446, https://doi.org/10.1128/AEM.70.9.5441-5446.2004 (2004).

64. Hutter, J. L. & Bechhoefer, J. Calibration of Atomic-Force Microscope Tips. *Rev Sci Instrum* 64, 3342–3342, https://doi.org/10.1063/1.1144449 (1993).

65. Dicker, W. A., Senden, T. J. & Pashley, R. M. Direct Measurement of Colloidal Forces Using an Atomic Force Microscope. *Nature* 353, 239–241, https://doi.org/10.1038/353239a0 (1991).

66. Beaussart, A. et al. Quantifying the forces guiding microbial cell adhesion using single-cell force spectroscopy. *Nature Protocols* 9, 1049–1055, https://doi.org/10.1038/nprot.2014.066 (2014).

67. Zeng, G. H., Muller, T. & Meyer, R. L. Single-Cell Force Spectroscopy of Bacteria Enabled by Naturally Derived Proteins. *Langmuir: the ACS journal of surfaces and colloids* 30, 4019–4025, https://doi.org/10.1021/la404673q (2014).

68. Gaboriaud, F. & Dufrene, Y. F. Atomic force microscopy of microbial cells: application to nanomechanical properties, surface forces and molecular recognition forces. *Colloids Surf B Biointerfaces* 54, 10–19, https://doi.org/10.1016/j.col surfb.2006.09.014 (2007).

69. Ahimou, F., Denis, F. A., Touhami, A. & Dufrêne, Y. F. Probing microbial cell surface charges by atomic force microscopy. *Langmuir*, https://doi.org/10.1021/la026273k (2002).

70. Alsteens, D., Dague, E., Rouxhet, P. G., Baulard, A. R. & Dufrene, Y. F. Direct measurement of hydrophobic forces on cell surfaces using AFM. *Langmuir* 23, 11977–11979, https://doi.org/10.1021/la702765c (2007).

## Acknowledgements

This work was supported in part by INRA (Institut National de la Recherche Agronomique) and by the CoMBa ("Colonisation de la Matrice extracellulaire Bacterienne") project funded by the "Région Auvergne" FRI ("Fond Régional Innovation") cluster IRP ("Institut de Recherche Pharmabiotique"). Support from the Australian Research Council (ARC) (DP150102287, DP180102987, FT130100580) is also acknowledged. Valentin AGEORGES is a PhD Research Fellow granted by the "Conseil Régional d'Auvergne - FEDER (Fonds Européen de Développement Régional)". The authors thank Ulrich Dobrindt for providing with *E. coli* K12 MG1655 Δ*fimΔagn43*. The authors are grateful to the INRA MIGALE bioinformatics platform (http://migale.jouy.inra.fr) for providing computational resources as well as to the INRA MEDI*S* microscopy platform for access to the epifluorescence microscope.

## Author Contributions

V.A. and M.D. conceptualise the overarching aims of the research study. V.A., M.S., E.D., G.J., B.H., I.R.H. and M.D. conceived and designed the experiments. V.A., M.S., N.C., P.R., I.C., X.B. and A.E.R. performed the experiments and data acquisition. V.A., M.S., G.J., P.R., I.C., X.B., S.L., J.P., B.H., F.C.D., I.R.H. and M.D. analysed and interpreted the data. M.D. had management as well as coordination responsibility for the execution of the research work. M.D. and B.H. contributed to the acquisition of the financial supports and resources leading to this publication. V.A., M.S., G.J., X.B., E.D., S.L., J.P., B.H., F.C.D., I.R.H. and M.D. wrote the article, including drafting and revising critically the manuscript for important intellectual content.

## Additional Information

Supplementary information accompanies this paper at https://doi.org/10.1038/s41598-019-47608-4.

**Competing Interests:** F.C.D. and M.S. are employees of the Lallemand SAS Company, which provided support in the form of salaries and partly financed their work. All other authors declare no competing interests.
**Publisher's note:** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_11.png" /></figure>

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International

License, which permits use, sharing, adaptation, distribution and reproduction in any medium or
format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Cre-
ative Commons license, and indicate if changes were made. The images or other third party material in this
article are included in the article's Creative Commons license, unless indicated otherwise in a credit line to the
material. If material is not included in the article's Creative Commons license and your intended use is not per-
mitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the
copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

© The Author(s) 2019