Early biofilm and streamer formation is mediated by wall shear stress and surface wettability: A multifactorial microfluidic study

Abstract Biofilms are intricate communities of microorganisms encapsulated within a self‐produced matrix of extra‐polymeric substances (EPS), creating complex three‐dimensional structures allowing for liquid and nutrient transport through them. These aggregations offer constituent microorganisms enhanced protection from environmental stimuli—like fluid flow—and are also associated with higher resistance to antimicrobial compounds, providing a persistent cause of concern in numerous sectors like the marine (biofouling and aquaculture), medical (infections and antimicrobial resistance), dentistry (plaque on teeth), food safety, as well as causing energy loss and corrosion. Recent studies have demonstrated that biofilms interact with microplastics, often influencing their pathway to higher trophic levels. Previous research has shown that initial bacterial attachment is affected by surface properties. Using a microfluidic flow cell, we have investigated the relationship between both wall shear stress (τw) and surface properties (surface wettability) upon biofilm formation of two species (Cobetia marina and Pseudomonas aeruginosa). We investigated biofilm development on low‐density polyethylene (LDPE) membranes, Permanox® slides, and glass slides, using nucleic acid staining and end‐point confocal laser scanning microscopy. The results show that flow conditions affect biomass, maximum thickness, and surface area of biofilms, with higher τw (5.6 Pa) resulting in thinner biofilms than lower τw (0.2 Pa). In addition, we observed differences in biofilm development across the surfaces tested, with LDPE typically demonstrating more overall biofilm in comparison to Permanox® and glass. Moreover, we demonstrate the formation of biofilm streamers under laminar flow conditions within straight micro‐channels.


| INTRODUCTION
Biofilms are collections of microorganisms adhered to a surface (living or inanimate), or as flocs, and encapsulated by a network of extrapolymeric substances (EPS), providing the constituent organisms with enhanced protection from both environmental stressors and antimicrobial substances (Costerton, 1999;H. Flemming & Wingender, 2010). Biofilms are the prominent growth form of bacteria and can be characterized by intricate, three-dimensional microstructures which allow for liquid through flow and the generation of nutrient gradients (Hans Curt Flemming & Wuertz, 2019;Kolter & Greenberg, 2006).
The detrimental effects of biofilm development are prevalent in a vast range of fields, including industrial, ecological, and medicinal settings (Donlan, 2001;Michael P. Schultz, 2007;Salta, Wharton, Blache, et al., 2013;Vertes et al., 2012). These profound negative consequences can be demonstrated clearly in a range of scenarios and sectors, for example in biofilm-related infections (medical), food processing, to the maritime sector (H. C. Flemming, 2002;Guaglianone et al., 2010;Høiby et al., 2011). Biofouling is the progressive accumulation of organisms upon submerged surfaces and biofilms, aka microfouling, act as the precursor stage (M. P. Schultz et al., 2011). The incremental build-up of fouling leads to rapid system clogging, biocorrosion, reduction of operational sensitivity in environmental sensors, and increased rates of hydrodynamic drag on ships (Delauney et al., 2010;Neria-González et al., 2006;M. P. Schultz et al., 2011). Within medical settings, the presence and persistence of biofilms create a serious concern (Monteiro et al., 2009), for example, implants, such as urinary catheters, arterial stents, artificial joints, and dental implants are highly susceptible to biofilm formation, providing an artificial surface for microorganisms to colonize (Guaglianone et al., 2010;Percival et al., 2017). With the established increased rates of antimicrobial resistance associated with biofilms, indwelling infections are much harder to treat and remove (Vertes et al., 2012).
Research in this field has previously focused on improving and modifying the materials used in implants to increase their antibiofilm capability, with the incorporation of antimicrobial compounds (Jordan et al., 2015;Von Borowski et al., 2019).
Initial bacterial attachment and subsequent biofilm formation are intrinsically linked with the dynamic conditions of the surrounding environment (Dunsmore et al., 2002;Shumi et al., 2010). Factors such as pH, temperature, hydrodynamic forces, and the properties of the surface the biofilm colonizes dictate both the rate and eventual extent of biofilm formation (Jeong et al., 2014;Karimi et al., 2015;Stewart, 2012). Biofilm development is strongly influenced by the surrounding flow and the associated wall shear stress (τ w ) levels (Tsagkari & Sloan, 2018). It has been shown that the flow field can impact features such as community composition, physical structures, and growth (Purevdorj et al., 2002;Rupp et al., 2004;Wang et al., 2018). It has also been established that τ w can modulate the growth stages in biofilm development, with an extended immature stage demonstrating a clear adaptation to flow conditions (Rickard et al., 2004;Rochex et al., 2008). In the current research, we include a series of τ w levels that are simultaneously applied to investigate their effect on biofilm early establishment, allowing for detailed and direct comparisons between flow conditions. The adhesion of bacteria to a surface is the fundamental and primary process in biofilm development and depends upon the surface properties of any substrate being colonized, with attachment relying upon attraction forces between bacterial cells and the surface, described by the extended DLVO theory (named after Boris Derjaguin and Lev Landau, Evert Verwey and Theodoor Overbeek) (Katsikogianni et al., 2004;Tuson & Weibel, 2013). This close relationship between surface properties and bacterial attachment has been targeted within anti-biofilm research to identify strategies aimed at limiting the extent of biofilm formation (Bohinc et al., 2014).
Factors such as surface topography, wettability, and surface energy have been modified, alongside the addition of embedded antimicrobial substances, in successful attempts to reduce rates of bacterial attachment while also producing a surface from which biofilms can be more readily removed (Arpa-Sancet et al., 2012;Pasmore et al., 2002;Sanchis et al., 2006). Surfaces selected in the current work are characterized by different surface energies, allowing for a comparative quantification of the effect of substrate properties on biofilm growth.
Several biofilm morphologies associated with continued biofilm development have been shown to increase the detrimental impact of biofilms within flow cell environments (Drescher et al., 2013;Marty et al., 2012). Biofilm streamers, which are extensions of the bulk biomass that are suspended in the surrounding flow have gained notoriety, with increasing research aiming to characterize early formation and development due to their established links to higher rates of clogging, increased maintenance costs, decreased flow rates, limits to the efficiency of self-cleaning systems, and a reduction in an operational lifetime (Drescher et al., 2013;Stoodley et al., 1998).
Previous research into the fundamental aspects of streamer formation has concluded that the convergence of flow over biofilm can extend the biomass and cause an elongation of the EPS which trails behind, suspended in the flow (Rusconi et al., 2010). Experimental and mathematical simulations have demonstrated the formation of streamers creates a more streamlined biofilm profile, which reduces the hydrodynamic force the bulk biomass is subject to, resulting in a more resilient biofilm (Taherzadeh et al., 2010).
While previous research has evaluated the role of either fluid shear stress or surface properties on biofilm development, we have utilized a microfluidic-based platform to investigate these effects simultaneously (B. Li & Logan, 2004;Shumi et al., 2013). Microfluidic devices can be purpose-built and are entirely customizable to suit experimental designs across a vast range of fields including biotechnology, microbiology, and pharmaceutics (K. P. Kim et al., 2010;Sackmann et al., 2014;Zheng et al., 2016). Flow cells are a class of devices that involve the manipulation of liquids through a defined experimental area, to replicate flow metrics that are relevant to specific applications, i.e. τ w or fluid velocity in physiological systems (X. Li et al., 2012;Nance et al., 2013;Runyon et al., 2008).
These small-scale devices have allowed for detailed examination within biofilm research, on aspects such as bacterial attachment and biofilm development (K. P. Kim et al., 2010;Samarian et al., 2014).
In the current work, three different surfaces were used to test initial bacterial attachment and biofilm development; these substrates were chosen to provide a range of surface properties representative of the wide settings where biofilms are a cause for concern. Low-Density Polyethylene (LDPE) is a thermoplastic ubiquitously used in numerous applications and products (plastic wraps and bags, squeeze bottles, toys, and gas and water pipes).
Plastic and microplastic pollution is a major issue in the marine environment, therefore shedding light on bacterial colonization and biofilm formation on this material is crucial (Kooi et al., 2017;Michels et al., 2018;Rummel et al., 2017). The second surface used was Permanox ® , which is an inert surface frequently used for cell attachment and growth. Finally, the third surface used was glass, which is a standardized material used in every laboratory worldwide. Using a range of surfaces with different surface energies, we aimed to identify any difference in biofilm development and, in turn, demonstrate the close relationship with substrate properties.
Moreover, the occurrence of streamers in straight channels remains a largely unexplored area, which we further investigate in the current work. Existing research on biofilm development routinely utilizes micro-channels with complex internal geometries to induce flow disruption over and around biofilms, therefore increasing the prevalence of biofilm streamers (Rusconi et al., 2010(Rusconi et al., , 2011. With a significant increase in clogging associated with biofilm streamers, it is crucial to examine the processes that govern their formation to direct future research aimed at their prevention (Drescher et al., 2013;Stoodley et al., 1998). Furthermore, results from this multi-parametric experimental investigation could potentially inform the development of predictive computational models that could be employed in the design of effective anti-biofilm surfaces for a range of applications.

| Flow cell design and fabrication
The flow cell used within these experiments is a second-generation device, designed and fabricated in the same way as the precursor model, which is detailed in Salta, Wharton, Blache, et al. (2013).
Briefly, the channel architecture was micro-milled in a layer of polymethyl methacrylate (PMMA), using a Datron CAT3D-M6 milling machine (Datron Dynamic, Inc.). A recess was also milled within this layer, and a custom-built silicone gasket was positioned in the recess to allow for effective sealing between the PMMA layer and the substrate surface. The gasket was fabricated from polydimethylsiloxane (PDMS, Sylgard 184, Dow Corning), using a weight ratio of 10:1 between PDMS monomer and curing agent.
The device was designed to investigate the performance of foulingcontrol surfaces/coatings under different wall shear stress levels and to create a standardized test for the development of new coatings for the industries that eventually deploy them. While the number of parallel channels has decreased from six to four from the initial design, the original premise of channels with decreasing heights generating a range of wall shear stress levels remains. The internal dimensions of these channels have been altered and optimized to extend the range of τ w generated. The τ w values generated within these chambers are relative to those present in tidal, blood flow, or cooling systems generating highly representative results for each of these settings (Cowle et al., 2020;Ku et al., 1985;Manuel et al., 2007). Moreover, while in the original design the channel height decreased in a step-like fashion, in the present device a more gradual (e.g., tapered) transition between chambers of different heights was established. This was aimed at preventing the onset of vortical flow in these regions of the device, which may potentially act as entrapment sites for flowing bacteria.

| Contact angle measurements
Contact angle (CA) measurements for each of the three test surfaces used in this study were performed to relate surface characteristics to biofilm development (KSV Instruments LTD, CAM101). The static measurements were conducted with ultra-filtered deionized water and were repeated 10 times for each surface, and the average and standard deviation of the contact angle values were calculated.
Images were captured the moment the water droplet touched the surface and were used to calculate the contact angle.

| Bacterial attachment and biofilm development assays
The species used in these experiments included Cobetia marina ATCC25374, as it has been previously employed as a model species in multiple attachment assays (Mieszkin et al., 2012;Salta, Wharton, Dennington et al., 2013). C. marina aliquots were taken from cryopreserved stocks stored at -80°C, plated onto marine agar (BS DifcoTM Marine Agar 2216), and incubated at 25°C. After an initial growth period of 48h, a single colony was used to inoculate Sea Salt Peptone (SSP), made using 35 g/L of Sea salts (S9883, Sigma Aldrich) and 18 g/L of Peptone (LP0037, Oxoid). The liquid culture was incubated at 25°C under agitation at 80 rpm, for over 12 h and bacterial growth was measured using a Synergy H1 microplate reader (BioTek ® ) with optical density (OD) at λ = 600 nm (OD 600 ).
The second species used was Pseudomonas aeruginosa ATCC25668; this species has been used in medical biofilm research, with certain strains causing infections (Pasmore et al., 2002;Purevdorj et al., 2002). P. aeruginosa was taken from cryopreserved stocks, stored at -80°C, plated onto nutrient agar (Oxoid LP0013), and incubated at 37°C. After 48h, a single colony was transferred to nutrient broth (NB) (Oxoid CM0067) and incubated at 37°C for 12 h.
With both bacteria, experiments were initiated once the liquid cultures reached an optical density (OD) at λ = 600 nm (OD 600 ) of 0.2.
A media reservoir containing 500 mL of SSP for C. marina or NB for P. aeruginosa was connected to a peristaltic pump (Watson-Marlow series 323S) to generate a continuous flow through the closed set-up, as shown in Figure 2. The addition of a 0.22 μm sterile filter prevented contamination while allowing aeration of the liquid media. A dampener was used in this setup to create a steady flow, attenuating the pressure wave generated by the pump. All components were connected using silicone rubber tubing (Masterflex, with an internal diameter of 3.1 mm). The flow cell and test surfaces were disinfected using 70% ethanol while all tubing was autoclaved before the start of each experiment. Low-density polyethylene (LDPE) membrane material was obtained from Fisher Scientific, Nunc™ Permanox ® microscope slides from Thermo Scientific, and glass microscope slides from Jaytec Glass Ltd).
An initial 1-h period was used to condition both the test surfaces and the channels with the appropriate media alone (SSP or NB). A sterile syringe was then used to introduce 5 ml of the desired F I G U R E 2 Schematic illustrating the experimental set-up, showing the constituent parts including an image of the flow cell. Arrows illustrate the direction of flow. Dotted lines indicate the option of adding a CCD camera for real-time measurements. bacterial culture (OD 600 = 0.2), after which the bacteria were left to attach for 30 min under static conditions. The flow was then restarted and continued for three and a half hours. The average flow rate achieved through these experiments (14.85 ml/h) was later used in the numerical simulations to accurately quantify the wall shear stress levels on the bottom surface of the device. The flow rate was calculated as described in Table A1. All experiments were conducted at room temperature. Upon completion of the experiment, the nucleic acid stain SYTO™9 (green fluorescent nucleic acid stain that stains all bacteria within the sample; Molecular Probes) was introduced into the flow cell using a sterile syringe, and the set-up was stored in a dark environment for 25 min at room temperature under static conditions. After this incubation period, any excess stain was washed away by conveying 5 mL of PBS through the flow cell. This procedure did not remove any established biofilm as the applied flow rate for the washing step for both species was the same as the experimental. The experiments were repeated four times per surface, while the multiple channel design shown in Figure 1, allowed for additional four replications within each experiment, enabling robust statistical comparisons between experimental groups (see next section for detailed experimental replication).

| Microscopy and image processing
The samples were analyzed using a confocal laser scanning microscope, Zeiss LSM 5 Pascal (Carl Zeiss). For each surface and in each of the four chambers, 20 Z-stacks were taken (N = 20), totaling N = 160 stacks per experiment, resulting in a total of N = 320 stacks per surface (experiments were repeated four times). All Z-stack images were recorded at Z-intervals of 0.5 μm, using a 63X magnification water immersion lens. Images were taken at the center of the channels, at a distance of 0.2 mm away from the sidewalls and the tapered region connecting different chambers in a channel, which can be seen in Figure 1. These measures were taken to avoid recording biomass that may have been attached to the lateral or top surfaces of the channel and not the surfaces being tested. The light exposure for SYTO9 was at λ EX = 488 nm, and the emission was collected between λ EX = 500-600 nm. Volocity software (PerkinElmer, Inc.) was used to improve the quality and resolution of 3D image data sets.
The Z-stacks were processed using COMSTAT2, a plugin within the image processing software ImageJ (MacBiophotonics ImageJ) (Heydorn et al., 2000). A fixed threshold value and connected volume filtration were used throughout all image processing and analysis. Using COMSTAT2, the biofilm biomass (μm 3 /μm 2 ), maximum thickness (μm), area occupied in layers, and surface area (μm 2 ) were determined. Biomass is calculated as the volume of all voxels, above a threshold calculated as per Otsu's method (Otsu, 1979), that contains biomass divided by the total image area. Maximum thickness is defined by the highest point of the biofilm relevant to the substratum. The area occupied in layers is the biomass recorded within each slice of a Z-stack, and the surface area accounts for the area occupied by the biofilm.

| Statistical analysis
Data processed with COMSTAT2 (biomass, surface area, and maximum thickness) were analyzed for statistical differences using IBM SPSS statistics 24. To determine the homogeneity of variances, Levene's test was used. In cases where the data failed to meet the homogeneity of variances, a Kruskal-Wallis test was used; otherwise, a one-way analysis of variance (ANOVA) was used. In cases with multiple parameters, a multivariate analysis of variance (MANOVA) was applied. All conclusions were based on a 95% confidence level Regression analysis (linear fit) was performed to determine significant differences among the different shears and scatterplots were generated using OriginPro 2020b.

| Flow characterization
The computer-aided design elements of the flow cell, shown in Figure 1d,e, were used in CFD simulations to predict the wall shear stress field acting over the biofilm surface in the experiments. The flow regime generated in our experiments was characterized as steady and laminar throughout the whole device. As shown in Table 1, the τ w generated from a volumetric flow rate of 14.85 ml/min ranged from 0.1 Pa to 5.6 Pa along the flow channel. Given that the fluid velocity has a parabolic profile, the wall shear stress in close proximity to the side walls of the channel is lower when compared to the central region. The reported wall shear stress values have thus been determined at a distance >0.2 mm away from the lateral walls, where the wall shear stress is substantially uniform. Notably, the wall shear stress values achieved are within the range biofilms experience in dental environments and physiological or pathological arterial blood flow (Guaglianone et al., 2010;Ku et al., 1985;Nance et al., 2013;Samarian et al., 2014). It should also be noted that the overall residence time of flowing cells within the device is <1 s. In multiple cases for biomass, maximum thickness, and surface area, significant differences can be found when τ w increases by an order of magnitude of two (i.e., 0.2 Pa is significantly different from 4.1 Pa, but not 2.2 Pa), as shown in Figure 3 and Tables A2-A3. These observations suggest that the end-point biofilm formation metrics are inversely proportional to F I G U R E 3 Influence of shear stress on Cobetia marina on the left and Pseudomonas aeruginosa on the right-hand side when exposed to different surfaces. Results collected from each test surface are compiled above including the R² value for the linear fit, in the order of LDPE, Permanox ® and Glass, respectively. Error bars ± SE. Note: the associated statistical analysis can be found in Tables A2-A5. SE, standard error wall shear stress. Such trends support previous findings suggesting that changing shear stress alters formation patterns of bacterial biofilms, with higher τ w restricting formation more than lower τ w (Liu & Tay, 2002;Paul et al., 2012). For C. marina biomass on LDPE, a significant difference was observed between the lowest and highest τ w , with lower τ w supporting greater biomass (Table A2-

| Biofilm development is dependent upon the surface type
Biofilm formation was impacted by surface type, as shown in Figure 3; these findings are consistent with the current understanding that surface wettability plays an intrinsic role in biofilm formation and development (Pasmore et al., 2002;Michael P. Schultz, 2007;Tuson & Weibel, 2013;Zheng et al., 2021). This variation between surfaces is associated with surface properties and, by extension, the way cells interact with the surfaces themselves (Katsikogianni et al., 2004;Tuson & Weibel, 2013 In terms of maximum thickness, however, the type of surface had a more profound effect. C. marina biofilms showed thicker structures at the lower τ w on LDPE followed by Permanox ® , while the thinnest biofilms formed on glass (Figures 3 and 4 and Table A4). Once again, as τ w increased, the surface properties had a reduced effect as wall shear stress became the dominant governing factor.
The results for P. aeruginosa displayed a less varied biofilm colony distribution than the marine species, where differences between surfaces became less evident. Under all four wall shear stress values, significant differences were found between the two hydrophobic surfaces for biomass, with LDPE showing higher biomass in all cases (LDPE vs. Permanox ® , p < 0.001). As the τ w increased, the underlying effect of surface properties upon resulting biomass became evident in the case of the hydrophilic glass. Biofilm development on glass resulted in the sharpest decline with increasing wall shear stress (Figures 3 and 4), illustrating that both shear forces and surface properties affect biofilm development. For maximum thickness, the only significant difference between surfaces is found between LDPE and glass, at τ w of 0.2 Pa, 2.2 Pa, and 5.6 Pa (Table A5). There were some clear differences between the two species selected, since P. aeruginosa formed thinner biofilms on LDPE than on glass, while the opposite was found for C. marina (LDPE showing thicker biofilms than on glass). Surface roughness and stiffness appear to influence P. aeruginosa adhesion and c-di-GMP production (Zheng et al., 2021).
For instance, P. aeruginosa illustrated a reduced adhesion on a stainless steel rough surface (Ra: 172.5 nm) when compared with a polished one (Ra: 84.4-45.2 nm); while it has been illustrated that softer material promotes c-di-GMP expression and bacterial adhesion (Song et al., 2015). Material stiffness, as an influencing factor towards bacterial adhesion, is still being explored as thoroughly reviewed by CHUN ET AL.  Figure 4 (and Figures   A1-A6), at the lower τ w (0.2, 2.2 Pa) both species generated comparatively uniform, thicker biofilms covering a relatively larger area of the channel. As the τ w increased, the biomass decreased, promoting more dispersed and overall smaller clusters, as previously observed (Salta, Capretto, Carugo, et al., 2013). At the higher τ w levels (4.1 Pa, 5.6 Pa), the biofilm showed a smaller vertical profile, F I G U R E 4 Confocal laser scanning microscopy (CLSM) images of Cobetia marina and Pseudomonas aeruginosa biofilm on low-density polyethylene (LDPE), Permanox ® , and glass showing both the XY and XZ planes. In all images, the flow was oriented from left to right, with scale bars of 50 μm. These images demonstrate that overall biomass and biofilm thickness decrease as the wall shear stress levels increase.
indicating that the increased wall shear stress limited the vertical development of biofilms ( Figures A4-A6). A recent study illustrated that in biofilm formation under temperatures lower than 25°C, the level of intracellular 3',5'-cyclic diguanylate (c-di-GMP), which controls biofilm formation and total exopolysaccharide production, rapidly increases, resulting in more and better-structured biofilms (Kim et al., 2020). At this point, it should be noted that P. aeruginosa biofilm morphology may have also been impacted by the temperature of the experimental setup, which has not been cardinal for this species. And although a temperature effect on biofilm formation cannot be excluded, our experimental duration has been significantly shorter (a total of 4 h) when compared to the study by Kim et al. (2020)  where some exemptions include modified glass and metallic surfaces (Bohinc et al., 2014;Oder et al., 2018). Some studies have also looked at biofilm growth on plastic surfaces targeting water distribution systems, therefore exploring different flow and shear stress regimes (not within microfluidic devices) (Cowle et al., 2020;Manuel et al., 2007). Previous research has compared the attachment of a marine Pseudomonas sp. upon a range of surfaces, including polyethylene and glass, and recorded higher attachment to hydrophobic surfaces than hydrophilic (Fletcher & Loeb, 1979). The majority of bacteria genera have a negative net charge, typically defined by zeta potential measurements (Katsikogianni & Missirlis, 2010;Renner & Weibel, 2011); it is well established that electrostatic forces are key in determining bacterial cell attachment to a surface which is also expected to be charged. This interaction can be influenced by the medium's ionic strength, while small molecules, proteins, and ions can alter the surface chemistry and charge via diffusion and mass transport (Renner & Weibel, 2011). In the current study, the surface charge of our test material was not assessed, however, it has been shown that high shear stress can impact the expected surface/ bacterial interactions that are based on colloid theories (e.g., DLVO and extended DLVO) and macromolecule binding considerations. Katsikogianni and Missirlis (2010)

| Observation of streamers
Biofilm streamers of both species were observed upon all surfaces, but exclusively under the highest wall shear stress value (5.6 Pa).
These features were observed as small extensions of biofilm clusters, oriented in a similar direction as the flow, as can be seen in Figure 5 (and Figures A1 and A2). Interestingly, a recent study by Zhang et al. (2017) found that at shear stresses exceeding 200 Pa, up to 25% of P. aeruginosa cells adhered tenaciously on a range of tested surfaces even at shear stresses as high as 2000 Pa. It was shown that this subpopulation of P. aeruginosa In our study, streamer formation was also present on all surfaces at the highest τ w regardless of species and surface type. Therefore, it would be very interesting to explore the possibility of a shearselective presence of SSP cells with increasing shear stress (as supported by our flow cell design) and a potential for streamer formation by SSP cells. In addition, a study by Rodesney et al. (2017) revealed an increase of c-di-GMP with shear for P. aeruginosa, therefore it would be interesting to explore the role of this intracellular secondary messenger in streamer formation under shear stress conditions.

| SUMMARY AND CONSIDERATIONS FOR FUTURE WORK
In this study, we have demonstrated that higher wall shear stress levels produce overall thinner biofilms than lower wall shear stress; these findings are consistent with previous research showing that increased shear reduced biofilm thickness (Liu & Tay, 2002;Paul et al., 2012). We have created flow dynamic conditions characterized by a range of wall shear stress levels (0.2-5.6 Pa); varying the inlet flow rate will allow for the recreation of an even wider range of flow conditions which will ultimately relate to a broader field of applications (Samarian et al., 2014). An increased range of wall shear stress values will also allow for a comprehensive analysis of the relationship between shear stress and biofilm development, creating a unique database that could be applicable to an extensive range of settings.
We have shown that the developed flow cell enables high throughput evaluation of different surface types and that the surface properties play a significant role in biofilm development. Therefore, future experiments may investigate a wider range of surface types, relevant to a spectrum of different applications (e.g., other polymers, coatings, and metals). Moreover, alongside contact angle calculations, surface roughness should be evaluated as it can affect bacterial cell attachment (Gharechahi et al., 2012;Song et al., 2015). This is also evident from the current study; for instance, Permanox ® and LDPE short-term experiments (Jeong et al., 2014). The inclusion of modified bacterial strains may also highlight some of the key factors governing formation responses exhibited through experimental investigations (Drescher et al., 2013). While experiments in this study concern mono-species biofilms, the dominant form of wild-type biofilms is often multi-species (Elias & Banin, 2012;Rendueles & Ghigo, 2012;Rickard et al., 2003). The complex composition of such biofilms critically influences their form, with constituent species playing specific roles within Yang et al., 2011). The flow cell used in the current work could be used to investigate the effect of flow upon multi-species biofilms, later drawing direct comparisons against experiments using individual species. The incorporation of molecular analysis techniques could also address whether gene expression differs from multi-to mono-species, aiming to define any potential mechanism that increases resistance to shear stress. Existing research has shown that wall shear stress also influences biofilm community composition, typically reducing overall diversity (Rickard et al., 2004;Rochex et al., 2008). It has been shown that shear stress maintains the biofilm in a young state which is characterized by lower diversity (Rochex et al., 2008

CONFLICTS OF INTEREST
None declared.

DATA AVAILABILITY STATEMENT
All data generated or analyzed during this study are included in this published article.

ETHICS STATEMENT
None required.

APPENDIX A
See Tables A1-A5 and Figures A1-A7. T A B L E A1 The flow rate was calculated at the end of each experiment with the average time (N=5) taken to fill a 10 ml vessel recorded and used to determine the overall millilitre per minute