Investigation of the mesoscale structure and volumetric features of biofilms using optical coherence tomography

Optical coherence tomography in soft matter

Optical coherence tomography (OCT) has become an indispensable tool for investigating mesoscopic features in soft matter and fluid mechanics. Its ability to provide high-resolution, non-invasive measurements in both spatial and temporal domains bridges critical gaps in experimental instrumentation, enabling the study of complex, confined, and dynamic systems. This review serves as both an introduction to OCT and a practical guide for researchers seeking to adopt this technology. A set of tutorials, complemented by Python scripts, is provided for both intensity- and Doppler-based techniques. The versatility of OCT is illustrated through case studies, including time-resolved velocimetry, particle-based velocity measurements, slip velocity characterization, detection of shear-induced structures, and analysis of fluid-fluid and fluid-structure interactions. Drawing on our experiences, we also present a set of practical guidelines for avoiding common pitfalls.

The role of fluid friction in streamer formation and biofilm growth

Biofilms constitute one of the most common forms of living matter, playing an increasingly important role in technology, health, and ecology. While it is well established that biofilm growth and morphology are highly dependent on the external flow environment, the precise role of fluid friction has remained elusive. We grew Bacillus subtilis biofilms on flat surfaces of a channel in a laminar flow at wall shear stresses spanning one order of magnitude (τw = 0.068 Pa to τw = 0.67 Pa). By monitoring the three-dimensional distribution of biofilm over seven days, we found that the biofilms consist of smaller microcolonies, shaped like leaning pillars, many of which feature a streamer in the form of a thin filament that originates near the tip of the pillar. While the shape, size, and distribution of these microcolonies depend on the imposed shear stress, the same structural features appear consistently for all shear stress values. The formation of streamers occurs after the development of a base structure, suggesting that the latter induces a secondary flow that triggers streamer formation. Moreover, we observed that the biofilm volume grows approximately linearly over seven days for all shear stress values, with a growth rate inversely proportional to the wall shear stress. We develop a scaling model, providing insight into the mechanisms by which friction limits biofilm growth.

Automated image processing algorithm for 3D OCT images of fouling in spacer-filled membrane filtration channels

OCT imaging is an important technique to study fouling in spacer-filled channels of reverse osmosis systems for seawater desalination. However, OCT imaging of membrane filtration channels with feed spacers is challenging because the spacer material can be (partly) transparent, making it difficult to detect and possibly mistaken for fouling, and the longer optical pathway through the spacer material distorts the image below the spacer. This study presents an automated 3D OCT image processing method in MATLAB for visualization and quantification of biofouling in spacer-filled channels. First, a spacer template of arbitrary size and rotation was generated from a CT scan of the feed spacer. Second, background noise and file size were reduced by representing the OCT image with a list of discrete reflectors. Finally, the spacer template was overlayed with the feed spacer in the 3D OCT image, enabling automated visualization of the feed spacer and correction of the distortions. Moreover, the method allows the selection of datasets with the same location relative to the position of the spacer, enabling systematic comparison between datasets and quantitative analysis.•A spacer template of arbitrary size and rotation was generated from a CT scan.•The background noise was removed, and the file size was reduced by representing the OCT dataset with a list of discrete reflectors.•The spacer template was overlayed with the feed spacer in the 3D OCT image.

Biofilm growth dynamics in a micro-irrigation with reclaimed wastewater in the field scale

The dripper clogging due to the development of biofilm can reduce the benefits of micro-irrigation technology implementation using reclaimed wastewater. The narrow cross-section and labyrinth geometry of the dripper channel enhance the fouling mechanisms. The aim of this study was to evaluate the water distribution and biofouling of drip irrigation systems at the field scale during irrigation with treated wastewater. Six 100 m lines of commercial pipes with two pressure-compensating dripper types (flow rate, Q, of 0.65 L h-1 and 1.5 L h-1, respectively) were monitored for four months. Different zones along the pipes were selected to evaluate the influence of hydrodynamical conditions (Reynolds number = 5400 to 0) on biofouling. Destructive methods involving the biofilm extraction by mechanical means, showed little biofilm development without significant differences in dry and organic matter content in function of the sampling location along the pipe or dripper flow rate (Q0.65 and Q1.5). These results were confirmed by non-destructive methods, such as optical coherence tomography, that nevertheless showed that biofouling concerned 15-20% of the total dripper labyrinth volume. Total organic carbon monitoring and its composition (by three-dimensional excitation and emission matrix fluorescence microscopy) showed that the biofilm did not significantly influence the organic matter nature. Our results indicated that the biological activity and biofilm development in irrigation systems were more affected by the environmental conditions, particularly water temperature, rather than flow conditions. This confirmed that treated wastewater with low organic content can be used in micro-irrigation systems without significant loss of efficiency, even in conditions requiring intensive irrigation, such as the Mediterranean climate.

Mass transfer in heterogeneous biofilms: Key issues in biofilm reactors and AI-driven performance prediction

Biofilm reactors, known for utilizing biofilm formation for cell immobilization, offer enhanced biomass concentration and operational stability over traditional planktonic systems. However, the dense nature of biofilms poses challenges for substrate accessibility to cells and the efficient release of products, making mass transfer efficiency a critical issue in these systems. Recent advancements have unveiled the intricate, heterogeneous architecture of biofilms, contradicting the earlier view of them as uniform, porous structures with consistent mass transfer properties. In this review, we explore six biofilm reactor configurations and their potential combinations, emphasizing how the spatial arrangement of biofilms within reactors influences mass transfer efficiency and overall reactor performance. Furthermore, we discuss how to apply artificial intelligence in processing biofilm measurement data and predicting reactor performance. This review highlights the role of biofilm reactors in environmental and energy sectors, paving the way for future innovations in biofilm-based technologies and their broader applications.

A review of mechanism analysis methods in multi-species biofilm of foodborne pathogens

Biofilms are an aggregation of microorganisms that have high resistance to antimicrobial agents. In the food industry, it has been widely studied that foodborne pathogens on both food surfaces and food-contact surfaces can form biofilms thereby threatening the safety of the food. In the natural environment, multi-species biofilms formed by more than two different microorganisms are abundant. In addition, the resistance of multi-species biofilms to antimicrobial agents is higher than that of mono-species biofilms. Therefore, studies to elucidate the mechanisms of multi-species biofilms formed by foodborne pathogens are still required in the food industry. In this review paper, we summarized the novel analytical methods studied to evaluate the mechanisms of multi-species biofilms formed by foodborne pathogens by dividing them into four categories: spatial distribution, bacterial interaction, extracellular polymeric substance production and quorum sensing analytical methods.

Non-invasive investigation of the morphology and optical properties of the upside-down jellyfish Cassiopea with optical coherence tomography

The jellyfish Cassiopea largely cover their carbon demand via photosynthates produced by microalgal endosymbionts, but how holobiont morphology and tissue optical properties affect the light microclimate and symbiont photosynthesis in Cassiopea remain unexplored. Here, we use optical coherence tomography (OCT) to study the morphology of Cassiopea medusae at high spatial resolution. We include detailed 3D reconstructions of external micromorphology, and show the spatial distribution of endosymbionts and white granules in the bell tissue. Furthermore, we use OCT data to extract inherent optical properties from light-scattering white granules in Cassiopea, and show that granules enhance local light-availability for symbionts in close proximity. Individual granules had a scattering coefficient of µs = 200-300 cm-1, and scattering anisotropy factor of g = 0.7, while large tissue-regions filled with white granules had a lower µs = 40-100 cm-1, and g = 0.8-0.9. We combined OCT information with isotopic labelling experiments to investigate the effect of enhanced light-availability in whitish tissue regions. Endosymbionts located in whitish tissue exhibited significantly higher carbon fixation compared to symbionts in anastomosing tissue (i.e. tissue without light-scattering white granules). Our findings support previous suggestions that white granules in Cassiopea play an important role in the host modulation of the light-microenvironment.

Role of Biofilms in Waste Water Treatment

Biofilm cells have a different physiology than planktonic cells, which has been the focus of most research. Biofilms are complex biostructures that form on any surface that comes into contact with water on a regular basis. They are dynamic, structurally complex systems having characteristics of multicellular animals and multiple ecosystems. The three themes covered in this review are biofilm ecology, biofilm reactor technology and design, and biofilm modeling. Membrane-supported biofilm reactors, moving bed biofilm reactors, granular sludge, and integrated fixed-film activated sludge processes are all examples of biofilm reactors used for water treatment. Biofilm control and/or beneficial application in membrane processes are improving. Biofilm models have become critical tools for biofilm foundational research as well as biofilm reactor architecture and design. At the same time, the differences between biofilm modeling and biofilm reactor modeling methods are acknowledged.

Selective Optical Imaging for Detection of Bacterial Biofilms in Tissues

SIGNIFICANCE

The development of an imaging technique to accurately identify biofilm regions on tissues and in wounds is crucial for the implementation of precise surface-based treatments, leading to better patient outcomes and reduced chances of infection.

AIM

The goal of this study was to develop an imaging technique that relies on selective trypan blue (TB) staining of dead cells, necrotic tissues, and bacterial biofilms, to identify biofilm regions on tissues and wounds.

APPROACH

The study explored combinations of ambient multi-colored LED lights to obtain maximum differentiation between stained biofilm regions and the underlying chicken tissue or glass substrate during image acquisition. The TB imaging results were then visually and statistically compared to fluorescence images using a shape similarity measure.

RESULTS

The comparisons between the proposed TB staining method and the fluorescence standard used to detect biofilms on tissues and glass substrates showed up to 97 percent similarity, suggesting that the TB staining method is a promising technique for identifying biofilm regions.

CONCLUSIONS

The TB staining method demonstrates significant potential as an effective imaging technique for the identification of fluorescing and non-fluorescing biofilms on tissues and in wounds. This approach could lead to improved precision in surface-based treatments and better patient outcomes.

Formation and emergent dynamics of spatially organized microbial systems

Spatial organization is the norm rather than the exception in the microbial world. While the study of microbial physiology has been dominated by studies in well-mixed cultures, there is now increasing interest in understanding the role of spatial organization in microbial physiology, coexistence and evolution. Where studied, spatial organization has been shown to influence all three of these aspects. In this mini review and perspective article, we emphasize that the dynamics within spatially organized microbial systems (SOMS) are governed by feedbacks between local physico-chemical conditions, cell physiology and movement, and evolution. These feedbacks can give rise to emergent dynamics, which need to be studied through a combination of spatio-temporal measurements and mathematical models. We highlight the initial formation of SOMS and their emergent dynamics as two open areas of investigation for future studies. These studies will benefit from the development of model systems that can mimic natural ones in terms of species composition and spatial structure.

Separate effect of turbulent pulsation on internal mass transfer in porous biofilms

Turbulence is considered to be the result of coupled time-averaged and pulsating velocities, making it difficult to distinguish the respective effects, and the quantitative effect of turbulent pulsation alone on mass transfer within biofilms has not been discussed in the literature. In this study, we constructed a special oscillating grid biofilm reactor combining Particle Image Velocimetry (PIV) measurements and Computational Fluid Dynamics (CFD) simulations to achieve nearly isotropic turbulence in a designed ambient without time-averaged velocity and shear stress. Subsequently, velocity and contaminant concentration distributions were obtained by solving a mass transfer model with a k-ε turbulence model, combined with measurements of biofilm structure parameters. The results showed that the increase in turbulent pulsation intensity led to a significant stratification of the percolation velocity gradient in biofilms, which enhanced convective mass transfer. The changes of biofilm density and porosity under turbulent pulsation were more strongly correlated with convective mass transfer. When the turbulent intensity (q) increased to 2.50 cm s^-1, the removal rate reached the highest value of 96.93%, accelerating the migration of contaminant concentration and the diffusive mass transfer effect was obvious. In addition, the trend of biofilm thickness under turbulent pulsation was consistent with the change of contaminant concentration distribution, and the correlation between them was greater. In summary, at q of 2.50 cm s^-1, there was a positive effect on both convection and diffusion mechanisms in biofilms, and the contaminant removal rate and biofilm thickness reached the maximum, which was the recommended turbulent pulsation conditions.

Peaks, pores, and dragon eggs: Uncovering and quantifying the heterogeneity of treatment wetland biofilm matrices

Biofilms serve to house diverse microbial communities, which are responsible for the majority of wastewater constituent degradation and transformation in treatment wetlands (TWs). TW biofilm has been generally conceptualized as a relatively uniform film covering available surfaces. However, no studies attaining direct visual 3D representations of biofilm morphology have been conducted. This study focuses on imaging the morphology of detached, gravel-associated, and rhizospheric (Phalaris arundinacea) biofilms from subsurface TW mesocosms. Images obtained through both traditional light microscopy, environmental scanning electron microscopy (E-SEM) and Wet-SEM revealed that TW biofilms are structurally heterogeneous ranging from corrugated films to clusters of aggregates. Features such as water channels and pores were observed suggesting that pollutant transport inside biofilms is complex, and that the interfacial surface area between water and biofilm is much larger than previously understood. Biofilm thickness generally ranged between 170 and 240 μm, with internal biofilm porosities estimated as 34 ± 10 %, reaching a maximum of 50 %. Internal biofilm matrix pore diameters ranged from 1 to 205.2 μm, with a distribution that favored pores and channels smaller than 10 μm, and a mean equivalent spherical diameter of 8.6 μm. Based on the large variation in pore and channel sizes it is expected that a variety of flow regimes and therefore pollutant dynamics are likely to occur inside TW biofilm matrices. Based on the visual evidence and analysis, a new conceptual model was created to reflect the microscale TW biofilm dynamics and morphology. This new conceptual model will serve to inform future biokinetic modelling, microscale hydrology, microbial community assessment, and pollutant treatment studies.

Endotracheal tube biofilm in critically ill patients during the COVID-19 pandemic : description of an underestimated microbiological compartment

Biofilm (BF) growth is believed to play a major role in the development of ventilator-associated pneumonia (VAP) in the intensive care unit. Despite concerted efforts to understand the potential implication of endotracheal tube (ETT)-BF dispersal, clinically relevant data are lacking to better characterize the impact of its mesostructure and microbiological singularity on the occurrence of VAP. We conducted a multicenter, retrospective observational study during the third wave of the COVID-19 pandemic, between March and May 2021. In total, 64 ETTs collected from 61 patients were included in the present BIOPAVIR study. Confocal microscopy acquisitions revealed two main morphological aspects of ETT-deposited BF: (1) a thin, continuous ribbon-shaped aspect, less likely monobacterial and predominantly associated with Enterobacter spp., Streptococcus pneumoniae or Viridans streptococci, and (2) a thicker, discontinuous, mushroom-shaped appearance, more likely characterized by the association of bacterial and fungal species in respiratory samples. The microbiological characterization of ETT-deposited BF found higher acquired resistance in more than 80% of analyzed BF phenotypes, compared to other colonization sites from the patient's environment. These findings reveal BF as a singular microbiological compartment, and are of added clinical value, with a view to future ETT-deposited BF-based antimicrobial stewardship in critically ill patients. Trial registration NCT04926493. Retrospectively registered 15 June 2021.

Synthetic periphyton as a model system to understand species dynamics in complex microbial freshwater communities

Phototrophic biofilms, also known as periphyton, are microbial freshwater communities that drive crucial ecological processes in streams and lakes. Gaining a deep mechanistic understanding of the biological processes occurring in natural periphyton remains challenging due to the high complexity and variability of such communities. To address this challenge, we rationally developed a workflow to construct a synthetic community by co-culturing 26 phototrophic species (i.e., diatoms, green algae, and cyanobacteria) that were inoculated in a successional sequence to create a periphytic biofilm on glass slides. We show that this community is diverse, stable, and highly reproducible in terms of microbial composition, function, and 3D spatial structure of the biofilm. We also demonstrate the ability to monitor microbial dynamics at the single species level during periphyton development and how their abundances are impacted by stressors such as increased temperature and a herbicide, singly and in combination. Overall, such a synthetic periphyton, grown under controlled conditions, can be used as a model system for theory testing through targeted manipulation.

In vivo characterization of electroactive biofilms inside porous electrodes with MR Imaging

Identifying the limiting processes of electroactive biofilms is key to improve the performance of bioelectrochemical systems (BES). For modelling and developing BES, spatial information of transport phenomena and biofilm distribution are required and can be determined by Magnetic Resonance Imaging (MRI) in vivo, in situ and in operando even inside opaque porous electrodes. A custom bioelectrochemical cell was designed that allows MRI measurements with a spatial resolution of 50 μm inside a 500 μm thick porous carbon electrode. The MRI data showed that only a fraction of the electrode pore space is colonized by the Shewanella oneidensis MR-1 biofilm. The maximum biofilm density was observed inside the porous electrode close to the electrode-medium interface. Inside the biofilm, mass transport by diffusion is lowered down to 45% compared to the bulk growth medium. The presented data and the methods can be used for detailed models of bioelectrochemical systems and for the design of improved electrode structures.

The use of magnetic resonance imaging can contribute to a better understanding of limiting processes occurring in electroactive biofilms especially inside opaque porous electrodes.

Characterization of biofilm structure and properties via processing of 2D optical coherence tomography images in BISCAP

MOTIVATION

Processing of Optical Coherence Tomography (OCT) biofilm images is currently restricted to a set of custom-made MATLAB scripts. None of the tools currently available for biofilm image processing (including those developed for Confocal Laser Scanning Microscopy-CLSM) enable a fully automatic processing of 2D OCT images.

RESULTS

A novel software tool entitled Biofilm Imaging and Structure Classification Automatic Processor (BISCAP) is presented. It was developed specifically for the automatic processing of 2D OCT biofilm images. The proposed approach makes use of some of the key principles used in CLSM image processing, and introduces a novel thresholding algorithm and substratum detection strategy. Two complementary pixel continuity checks are executed, enabling very detailed pixel characterizations. BISCAP delivers common structural biofilm parameters and a set of processed images for biofilm analysis. A novel biofilm 'compaction parameter' is suggested. The proposed strategy was tested on a set of 300 images with highly satisfactory results obtained. BISCAP is a Python-based standalone application, not requiring any programming knowledge or property licenses, and where all operations are managed via an intuitive Graphical User Interface. The automatic nature of this image processing strategy decreases biasing problems associated to human-perception and allows a reliable comparison of outputs.

AVAILABILITY AND IMPLEMENTATION

BISCAP and a collection of biofilm images obtained from OCT scans can be found at: https://github.com/diogonarciso/BISCAP.

SUPPLEMENTARY INFORMATION

Supplementary data are available at Bioinformatics online.

Controlling the hydraulic resistance of membrane biofilms by engineering biofilm physical structure

The application of membrane technology for water treatment and reuse is hampered by the development of a microbial biofilm. Biofilm growth in micro-and ultrafiltration (MF/UF) membrane modules, on both the membrane surface and feed spacer, can form a secondary membrane and exert resistance to permeation and crossflow, increasing energy demand and decreasing permeate quantity and quality. In recent years, exhaustive efforts were made to understand the chemical, structural and hydraulic characteristics of membrane biofilms. In this review, we critically assess which specific structural features of membrane biofilms exert resistance to forced water passage in MF/UF membranes systems applied to water and wastewater treatment, and how biofilm physical structure can be engineered by process operation to impose less hydraulic resistance ("below-the-pain threshold"). Counter-intuitively, biofilms with greater thickness do not always cause a higher hydraulic resistance than thinner biofilms. Dense biofilms, however, had consistently higher hydraulic resistances compared to less dense biofilms. The mechanism by which density exerts hydraulic resistance is reported in the literature to be dependant on the biofilms' internal packing structure and EPS chemical composition (e.g., porosity, polymer concentration). Current reports of internal porosity in membrane biofilms are not supported by adequate experimental evidence or by a reliable methodology, limiting a unified understanding of biofilm internal structure. Identifying the dependency of hydraulic resistance on biofilm density invites efforts to control the hydraulic resistance of membrane biofilms by engineering internal biofilm structure. Regulation of biofilm internal structure is possible by alteration of key determinants such as feed water nutrient composition/concentration, hydraulic shear stress and resistance and can engineer biofilm structural development to decrease density and therein hydraulic resistance. Future efforts should seek to determine the extent to which the concept of "biofilm engineering" can be extended to other biofilm parameters such as mechanical stability and the implication for biofilm control/removal in engineered water systems (e.g., pipelines and/or, cooling towers) susceptible to biofouling.

Data-driven modeling of heterogeneous viscoelastic biofilms

Biofilms are typically heterogeneous in morphology, structure, and composition, resulting in non-uniform mechanical properties. The distribution of mechanical properties, in turn, determines the biofilm behavior, such as deformation and detachment. Most biofilm models neglect biofilm heterogeneity, especially at the microscale. In this study, an image-based modeling approach was developed to transform two-dimensional optical coherence tomography (OCT) biofilm images to a pixel-scale non-Newtonian viscosity map of the biofilm. The map was calibrated using the bulk viscosity data from rheometer tests, based on assumed maximum and minimum viscosities and a relationship between OCT image intensity signals and non-Newtonian viscosity. While not quantitatively measuring biofilm viscosity for each pixel, it allows a rational spatial allocation of viscosities within the biofilm: areas with lower cell density, e.g., voids, are assigned lower viscosities, and areas with high cell densities are assigned higher viscosities. The spatial distribution of non-Newtonian viscosity was applied in an established Oldroyd-B constitutive model and implemented using the phase-field continuum approach for the deformation and stress analysis. The heterogeneous model was able to predict deformations more accurately than a homogenous one. Stress distribution in the heterogeneous biofilm displayed better characteristics than that in the homogeneous one, because it is highly dependent on the viscosity distribution. This work, using a pixel-scale, image-based approach to map the mechanical heterogeneity of biofilms for computational deformation and stress analysis, provides a novel modeling approach that allows the consideration of biofilm structural and mechanical heterogeneity. Future research should better characterize the relationship between OCT signal and viscosity, and consider other constitutive models for biofilm mechanical behavior. This article is protected by copyright. All rights reserved.

Low-Field Nuclear Magnetic Resonance Characteristics of Biofilm Development Process

To in situ and noninvasively monitor the biofilm development process by low-field nuclear magnetic resonance (NMR), experiments should be made to determine the mechanisms responsible for the T₂ signals of biofilm growth. In this paper, biofilms were cultivated in both fluid media and saturated porous media. T₂ relaxation for each sample was measured to investigate the contribution of the related processes to T₂ relaxation signals. In addition, OD values of bacterial cell suspensions were measured to provide the relative number of bacterial cells. We also obtained SEM photos of the biofilms after vacuum freeze-drying the pure sand and the sand with biofilm formation to confirm the space within the biofilm matrix and identify the existence of biofilm formation. The T₂ relaxation distribution is strongly dependent on the density of the bacterial cells suspended in the fluid and the stage of biofilm development. The peak time and the peak percentage can be used as indicators of the biofilm growth states.

Insights into the Development of Phototrophic Biofilms in a Bioreactor by a Combination of X-ray Microtomography and Optical Coherence Tomography

As productive biofilms are increasingly gaining interest in research, the quantitative monitoring of biofilm formation on- or offline for the process remains a challenge. Optical coherence tomography (OCT) is a fast and often used method for scanning biofilms, but it has difficulty scanning through more dense optical materials. X-ray microtomography (μCT) can measure biofilms in most geometries but is very time-consuming. By combining both methods for the first time, the weaknesses of both methods could be compensated. The phototrophic cyanobacterium Tolypothrix distorta was cultured in a moving bed photobioreactor inside a biocarrier with a semi-enclosed geometry. An automated workflow was developed to process µCT scans of the biocarriers. This allowed quantification of biomass volume and biofilm-coverage on the biocarrier, both globally and spatially resolved. At the beginning of the cultivation, a growth limitation was detected in the outer region of the carrier, presumably due to shear stress. In the later phase, light limitations could be found inside the biocarrier. µCT data and biofilm thicknesses measured by OCT displayed good correlation. The latter could therefore be used to rapidly measure the biofilm formation in a process. The methods presented here can help gain a deeper understanding of biofilms inside a process and detect any limitations.

Visualizing, quantifying, and controlling local hydrodynamic effects on biofilm accumulation in complex flow paths

Complex flow paths (CFPs) are commonly applied in precision equipment to accurately supply controllable fluids with designed structures. However, the presence of biofilms in CFPs causes quite a few unwanted issues, such as bio-erosion, clogging, or even health risks. To date, visualizing and quantifying the interaction between biofilm distribution and local hydrodynamics remains difficult, and the mechanism during the process is unclear. In this paper, the remodeling simulation method (3D industrial computed tomography scanning-inverse modeling-numerical simulation) and 16S rRNA high-throughput sequencing were integrated. The results indicated that local hydrodynamic characteristics significantly affected biofilm thicknesses on CFP surfaces (relative differences of 41.3-71.2%), which inversely influenced the local turbulence intensity. The average biofilm thicknesses exhibited a significant quadratic correlation with the near-wall hydraulic shear forces (r > 0.72, p < 0.05), and the biofilm reached a maximum thickness at 0.36-0.45 Pa. On the other hand, the near-wall hydraulic shear forces not only affected microbial community characteristics of biofilms, but they also influenced the number of microorganisms involved, which determined the biofilm accumulation thereafter. The PHYLUM Firmicutes and Proteobacteria were the dominant bacteria during the process. The results obtained in this paper could provide practical conceptions for the targeted control of biofilms and put forward more efficient controlling methods in commonly applied CFP systems.

Assessment of the Impact of Temperature on Biofilm Composition with a Laboratory Heat Exchanger Module

Temperature change over the length of heat exchangers might be an important factor affecting biofouling. This research aimed at assessing the impact of temperature on biofilm accumulation and composition with respect to bacterial community and extracellular polymeric substances. Two identical laboratory-scale plate heat exchanger modules were developed and tested. Tap water supplemented with nutrients was fed to the two modules to enhance biofilm formation. One “reference” module was kept at 20.0 ± 1.4 °C and one “heated” module was operated with a counter-flow hot water stream resulting in a bulk water gradient from 20 to 27 °C. Biofilms were grown during 40 days, sampled, and characterized using 16S rRNA gene amplicon sequencing, EPS extraction, FTIR, protein and polysaccharide quantifications. The experiments were performed in consecutive triplicate. Monitoring of heat transfer resistance in the heated module displayed a replicable biofilm growth profile. The module was shown suitable to study the impact of temperature on biofouling formation. Biofilm analyses revealed: (i) comparable amounts of biofilms and EPS yield in the reference and heated modules, (ii) a significantly different protein to polysaccharide ratio in the EPS of the reference (5.4 ± 1.0%) and heated modules (7.8 ± 2.1%), caused by a relatively lower extracellular sugar production at elevated temperatures, and (iii) a strong shift in bacterial community composition with increasing temperature. The outcomes of the study, therefore, suggest that heat induces a change in biofilm bacterial community members and EPS composition, which should be taken into consideration when investigating heat exchanger biofouling and cleaning strategies. Research potential and optimization of the heat exchanger modules are discussed.

Morphogenesis and oxygen dynamics in phototrophic biofilms growing across a gradient of hydraulic conditions

Biofilms are surface-attached and matrix-enclosed microbial communities that dominate microbial life in numerous ecosystems. Using flumes and automated optical coherence tomography, we studied the morphogenesis of phototrophic biofilms along a gradient of hydraulic conditions. Compact and coalescent biofilms formed under elevated bed shear stress, whereas protruding clusters separated by troughs formed under reduced shear stress. This morphological differentiation did not linearly follow the hydraulic gradient, but a break point emerged around a shear stress of ~0.08 Pa. While community composition did not differ between high and low shear environments, our results suggest that the morphological differentiation was linked to biomass displacement and reciprocal interactions between the biofilm structure and hydraulics. Mapping oxygen concentrations within and around biofilm structures, we provide empirical evidence for biofilm-induced alterations of oxygen mass transfer. Our findings suggest that architectural plasticity, efficient mass transfer, and resistance to shear stress contribute to the success of phototrophic biofilms.

Non-Invasive Measurement, Mathematical Simulation and In Situ Detection of Biofilm Evolution in Porous Media: A Review

The development of biofilms and the related changes in porous media in the subsurface cannot be directly observed and evaluated. The primary reason that the mechanism of biofilm clogging in porous media cannot be clearly demonstrated is due to the opacity and structural complexity of three-dimensional pore space. Interest in exploring methods to overcome this limitation has been increasing. In the first part of this review, we introduce the underlying characteristics of biofilm in porous media. Then, we summarize two approaches, non-invasive measurement methods and mathematical simulation strategies, for studying fluid–biofilm–porous medium interaction with spatiotemporal resolution. We also discuss the advantages and limitations of these approaches. Lastly, we provide a perspective on opportunities for in situ monitoring at the field site.

Hydrodynamics and surface properties influence biofilm proliferation

A biofilm is an interface-associated colloidal dispersion of bacterial cells and excreted polymers in which microorganisms find protection from their environment. Successful colonization of a surface by a bacterial community is typically a detriment to human health and property. Insight into the biofilm life-cycle provides clues on how their proliferation can be suppressed. In this review, we follow a cell through the cycle of attachment, growth, and departure from a colony. Among the abundance of factors that guide the three phases, we focus on hydrodynamics and stratum properties due to the synergistic effect such properties have on bacteria rejection and removal. Cell motion, whether facilitated by the environment via medium flow or self-actuated by use of an appendage, drastically improves the survivability of a bacterium. Once in the vicinity of a stratum, a single cell is exposed to near-surface interactions, such as van der Waals, electrostatic and specific interactions, similarly to any other colloidal particle. The success of the attachment and the potential for detachment is heavily influenced by surface properties such as material type and topography. The growth of the colony is similarly guided by mainstream flow and the convective transport throughout the biofilm. Beyond the growth phase, hydrodynamic traction forces on a biofilm can elicit strongly non-linear viscoelastic responses from the biofilm soft matter. As the colony exhausts the means of survival at a particular location, a set of trigger signals activates mechanisms of bacterial release, a life-cycle phase also facilitated by fluid flow. A review of biofilm-relevant hydrodynamics and startum properties provides insight into future research avenues.

Exploring flow-biofilm-sediment interactions: Assessment of current status and future challenges

Biofilm activities and their interactions with physical, chemical and biological processes are of great importance for a variety of ecosystem functions, impacting hydrogeomorphology, water quality and aquatic ecosystem health. Effective management of water bodies requires advancing our understanding of how flow influences biofilm-bound sediment and ecosystem processes and vice-versa. However, research on this triangle of flow-biofilm-sediment is still at its infancy. In this Review, we summarize the current state of the art and methodological approaches in the flow-biofilm-sediment research with an emphasis on biostabilization and fine sediment dynamics mainly in the benthic zone of lotic and lentic environments. Example studies of this three-way interaction across a range of spatial scales from cell (nm - µm) to patch scale (mm - dm) are highlighted in view of the urgent need for interdisciplinary approaches. As a contribution to the review, we combine a literature survey with results of a pilot experiment that was conducted in the framework of a joint workshop to explore the feasibility of asking interdisciplinary questions. Further, within this workshop various observation and measuring approaches were tested and the quality of the achieved results was evaluated individually and in combination. Accordingly, the paper concludes by highlighting the following research challenges to be considered within the forthcoming years in the triangle of flow-biofilm-sediment: i) Establish a collaborative work among hydraulic and sedimentation engineers as well as ecologists to study mutual goals with appropriate methods. Perform realistic experimental studies to test hypotheses on flow-biofilm-sediment interactions as well as structural and mechanical characteristics of the bed. ii) Consider spatially varying characteristics of flow at the sediment-water interface. Utilize combinations of microsensors and non-intrusive optical methods, such as particle image velocimetry and laser scanner to elucidate the mechanism behind biofilm growth as well as mass and momentum flux exchanges between biofilm and water. Use molecular approaches (DNA, pigments, staining, microscopy) for sophisticated community analyses. Link varying flow regimes to microbial communities (and processes) and fine sediment properties to explore the role of key microbial players and functions in enhancing sediment stability (biostabilization). iii) Link laboratory-scale observations to larger scales relevant for management of water bodies. Conduct field experiments to better understand the complex effects of variable flow and sediment regimes on biostabilization. Employ scalable and informative observation techniques (e.g., hyperspectral imaging, particle tracking) that can support predictions on the functional aspects, such as metabolic activity, bed stability, nutrient fluxes under variable regimes of flow-biofilm-sediment.

Hydrodynamic effect on biofouling of milli-labyrinth channel and bacterial communities in drip irrigation systems fed with reclaimed wastewater

The clogging of drippers due to the development of biofilms reduces the benefits and is an obstacle to the implementation of drip irrigation technology in a reclaimed water context. The narrow section and labyrinth geometry of the dripper channel results the development of a heterogeneous flow behaviours with the vortex zones which it enhance the fouling mechanisms. The objective of this study was to analyse the influence of the three dripper types, defined by their geometric and hydraulic parameters, fed with reclaimed wastewater, on the biofouling kinetics and the bacterial communities. Using optical coherence tomography, we demonstrated that the inlet of the drippers (mainly the first baffle) and vortex zones are the most sensitive area for biofouling. Drippers with the lowest Reynolds number and average cross-section velocity v (1 l·h^-1) were the most sensible to biofouling, even if detachment events seemed more frequent in this dripper type. Therefore, dripper flow path with larger v should be consider to improve the anti-clogging performance. In addition, the dripper type and the geometry of the flow path influenced the structure of the bacterial communities from dripper biofilms. Relative abundancy of filamentous bacteria belonging to Chloroflexi phylum was higher in 1 l·h^-1 drippers, which presented a higher level of biofouling. However, further research on the role of this phylum in dripper biofouling is required.

Intra-colony channels in E. coli function as a nutrient uptake system

The ability of microorganisms to grow as aggregated assemblages has been known for many years, however their structure has remained largely unexplored across multiple spatial scales. The development of the Mesolens, an optical system which uniquely allows simultaneous imaging of individual bacteria over a 36 mm2 field of view, has enabled the study of mature Escherichia coli macro-colony biofilm architecture like never before. The Mesolens enabled the discovery of intra-colony channels on the order of 10 μm in diameter, that are integral to E. coli macro-colony biofilms and form as an emergent property of biofilm growth. These channels have a characteristic structure and re-form after total mechanical disaggregation of the colony. We demonstrate that the channels are able to transport particles and play a role in the acquisition of and distribution of nutrients through the biofilm. These channels potentially offer a new route for the delivery of dispersal agents for antimicrobial drugs to biofilms, ultimately lowering their impact on public health and industry.

Effect of slow biodegradable substrate addition on biofilm structure and reactor performance in two MBBRs filled with different support media

In this study, two moving-bed biofilm reactors (MBBR₁ and MBBR₂) filled with different size of carrier media (Kaldnes K₁ and Kaldnes K₁ micro, respectively) were subjected to soluble (sugar and sodium acetate (Ac)) substrate and mixture of soluble and particulate (particulate potato starch (PS)) substrate in a very high organic loading rate (12 kgCOD/m³·d) at different temperatures (26 and 15°C, in MBBR₁ and MBBR₂, respectively). The effects of carrier type and substrate on biofilm structure and reactor performance have been studied. Starch was removed by adsorption at the biofilm surface and hydrolyzed which caused substrate gradient in MBBR₁, however, hydrolyzed uniformly within biofilm in MBBR₂. The biofilm of MBBR₁ was irregular due to filamentous structure growth due to the substrate gradient, while, it was regular in MBBR₂ due to uniform distribution of substrate. The performance of both MBBRs in ammonium, COD and TN removal decreased significantly when the amount of small particles in the reactor increased owing to feeding by starch, which led to biomass density decline. The type of media affected the quantity and distribution of attached biomass, which in turn influenced the activity of specific microbial functional groups in the biofilm. The biofilm in MBBR₂ was thicker and consequently nitrogen removal by denitrification was much higher. The lower temperature did not affect negatively the reactor performance in MBBR₂.

Visualizing the macroscale spatial distributions of biofilms in complex flow channels using industrial computed tomography

The presence of biofilms in enclosed pipelines can lead to numerous deleterious issues. To date, it has been difficult to use optical imaging techniques to monitor the macroscale spatial distributions of biofilms. To address this concern, a combination of industrial computed tomography (ICT) and a contrast agent was explored to noninvasively visualize biofilms in three types of drip irrigation emitters. The results showed that ICT successfully observed and quantified the macroscale spatial distributions of biofilms. The complex hydrodynamic characteristics in the emitter channels affected the local distributions of biofilms. Biofilms were mainly attached to the lateral and medial faces and biomass decreased along the flow directions. Based on the distributions of biofilms, some emitter structural design defects were further diagnosed. Applying ICT in combination with the contrast agent could potentially provide a visual and effective way to reveal the formation mechanisms of biofilms and to optimize flow channel structures to avoid biofilm accumulations.

Computational and Experimental Investigation of Biofilm Disruption Dynamics Induced by High-Velocity Gas Jet Impingement

Knowledge of mechanisms promoting disruption though mechanical forces is essential in optimizing biofilm control strategies which rely on fluid shear. Our results provide insight into how biofilm disruption dynamics is governed by applied forces and fluid properties, revealing a mechanism for ripple formation and fluid-biofilm mixing. These findings have important implications for the rational design of new biofilm cleaning strategies with fluid jets, such as determining optimal parameters (e.g., jet velocity and position) to remove the biofilm from a certain zone (e.g., in dental hygiene or debridement of surgical site infections) or using antimicrobial agents which could increase the interfacial area available for exchange, as well as causing internal mixing within the biofilm matrix, thus disrupting the localized microenvironment which is associated with antimicrobial tolerance. The developed model also has potential application in predicting drag and pressure drop caused by biofilms on bioreactor, pipeline, and ship hull surfaces.

Experimental data showed that high-speed microsprays can effectively disrupt biofilms on their support substratum, producing a variety of dynamic reactions such as elongation, displacement, ripple formation, and fluidization. However, the mechanics underlying the impact of high-speed turbulent flows on biofilm structure is complex under such extreme conditions, since direct measurements of viscosity at these high shear rates are not possible using dynamic testing instruments. Here, we used computational fluid dynamics simulations to assess the complex fluid interactions of ripple patterning produced by high-speed turbulent air jets impacting perpendicular to the surface of Streptococcus mutans biofilms, a dental pathogen causing caries, captured by high-speed imaging. The numerical model involved a two-phase flow of air over a non-Newtonian biofilm, whose viscosity as a function of shear rate was estimated using the Herschel-Bulkley model. The simulation suggested that inertial, shear, and interfacial tension forces governed biofilm disruption by the air jet. Additionally, the high shear rates generated by the jet impacts coupled with shear-thinning biofilm property resulted in rapid liquefaction (within milliseconds) of the biofilm, followed by surface instability and traveling waves from the impact site. Our findings suggest that rapid shear thinning under very high shear flows causes the biofilm to behave like a fluid and elasticity can be neglected. A parametric sensitivity study confirmed that both applied force intensity (i.e., high jet nozzle air velocity) and biofilm properties (i.e., low viscosity and low air-biofilm surface tension and thickness) intensify biofilm disruption by generating large interfacial instabilities.

3D biofilm visualization and quantification on granular bioanodes with magnetic resonance imaging

The use of microbial fuel cells (MFCs) for wastewater treatment fits in a circular economy context, as they can produce electricity by the removal of organic matter in the wastewater. Activated carbon (AC) granules are an attractive electrode material for bioanodes in MFCs, as they are cheap and provide electroactive bacteria with a large surface area for attachment. The characterization of biofilm growth on AC granules, however, is challenging due to their high roughness and three-dimensional structure. In this research, we show that 3D magnetic resonance imaging (MRI) can be used to visualize biofilm distribution and determine its volume on irregular-shaped single AC granules in a non-destructive way, while being combined with electrochemical and biomass analyses. Ten AC granules with electroactive biofilm (i.e. granular bioanodes) were collected at different growth stages (3 to 21 days after microbial inoculation) from a multi-anode MFC and T₁-weighted 3D-MRI experiments were performed for three-dimensional biofilm visualization. With time, a more homogeneous biofilm distribution and an increased biofilm thickness could be observed in the 3D-MRI images. Biofilm volumes varied from 0.4 μL (day 4) to 2 μL (day 21) and were linearly correlated (R² = 0.9) to the total produced electric charge and total nitrogen content of the granular bioanodes, with values of 66.4 C μL^-1 and 17 μg N μL^-1, respectively. In future, in situ MRI measurements could be used to monitor biofilm growth and distribution on AC granules.

Long term performance and dynamics of microbial biofilm communities performing sulfur-oxidizing autotrophic denitrification in a moving-bed biofilm reactor

Sulfide-oxidizing autotrophic denitrification (SOAD) implemented in a moving-bed biofilm reactor (MBBR) is a promising alternative to conventional heterotrophic denitrification in mainstream biological nitrogen removal. The sulfide-oxidation intermediate - elemental sulfur - is crucial for the kinetic and microbial properties of the sulfur-oxidizing bacterial communities, but its role is yet to be studied in depth. Hence, to investigate the performance and microbial communities of the aforementioned new biosystem, we operated for a long term a laboratory-scale (700 d) SOAD MBBR to treat synthetic saline domestic sewage, with an increase of the surface loading rate from 8 to 50 mg N/(m²·h) achieved by shortening the hydraulic retention time from 12 h to 2 h. The specific reaction rates of the reactor were eventually increased up to 0.37 kg N/(m³·d) and 0.73 kg S/(m³·d) for nitrate reduction and sulfide oxidation with no significant sulfur elemental accumulation. Two sulfur-oxidizing bacterial (SOB) clades, Sox-independent SOB (SOB_I) and Sox-dependent SOB (SOB_II), were responsible for indirect two-step sulfur oxidation (S^2-→S⁰→SO₄^2-) and direct one-step sulfur oxidation (S^2-→SO₄^2-), respectively. The SOB_II biomass-specific electron transfer capacity could be around 2.5 times greater than that of SOB_I (38 mmol e^-/(gSOB_II·d) versus 15 mmol e^-/(gSOB_I·d)), possibly resulting in the selection of SOB_II over SOB_I under stress conditions (such as a shorter HRT). Further studies on the methods and mechanism of selecting of SOB_II over SOB_I in biofilm reactors are recommended. Overall, the findings shed light on the design and operation of MBBR-based SOAD processes for mainstream biological denitrification.

Viscoelastic response of Escherichia coli biofilms to genetically altered expression of extracellular matrix components

How the viscoelastic properties of the extracellular matrix affect the various biological functions conferred by biofilms is an important question in microbiology. In this study, the viscoelastic response of Escherichia coli biofilms to the genetically altered expression of extracellular matrix components was studied. Biofilms of the wild type E. coli MG1655 and its mutant strains producing different amounts of extracellular matrix components (curli, colanic acid, and poly-β-1,6-N-acetyl-d-glucosamine) were used to examine the viscoelastic behavior of biofilms grown at the solid-atmosphere interface. The results suggest that the presence of curli proteins dominates biofilm mechanical behavior. The rheological data indicate that the cohesive energy of the biofilm was the highest in the wild type strain. The results demonstrate the importance of extracellular matrix composition for biofilm mechanical properties. We propose that by genetically altering the expression of extracellular matrix polymers, bacteria are able to modulate the mechanical properties of their local environment in accordance with bulk environmental conditions.

Biofilm compressibility in ultrafiltration: A relation between biofilm morphology, mechanics and hydraulic resistance

Poroelastic fluid-structure interaction models were coupled to experimental data to determine the effects of biofilm spatial distribution of mechanical and hydraulic properties on the biofilm hydraulic resistance and compressibility in membrane filtration processes. Biofilms were cultivated on ultrafiltration membranes for 20 and 30 days under high (0.28 bar) and low (0.06 bar) transmembrane pressure (TMP), in dead-end filtration mode. Subsequently, biofilms were subjected to a compression/relaxation cycles by step-wise TMP changes. Structural deformation of biofilms during compression was observed in-situ using optical coherence tomography. Experimental results show that the observed increase in the biofilm hydraulic resistance during compression is not necessarily accompanied by a detectable biofilm thickness reduction. A dual-layer biofilm model with a dense base and porous top layer could explain these observed results. Because porosity controls indirectly the mechanical response of biofilms under compression, results could be described without assuming a gradient in mechanical properties within the biofilm. The biofilm surface roughness did not significantly influence the water flux in this study. However, the fraction of biofilm base layer directly exposed to bulk liquid could be a good indicator in the determination of water flux. The main implications of this study for the design and operation of low-pressure membrane systems (e.g., MF and UF with fouling layer being the main filtration resistance) lays in the selection of favorable operational TMP and biofilm morphology.

Predation creates unique void layer in membrane-aerated biofilms

The membrane-aerated biofilm reactor (MABR) is a novel wastewater treatment technology based on oxygen-supplying membranes. The counter diffusion of oxygen and electron donors in MABRs leads to unique behavior, and we hypothesized it also could impact predation. We used optical coherence tomography (OCT), microsensor analyses, and mathematical modeling to investigate predation in membrane-aerated biofilms (MABs). When protozoa were excluded from the inoculum, the MAB's OCT-observable void fraction was around 5%. When protozoa were included, the void fraction grew to nearly 50%, with large, continuous voids at the base of the biofilm. Real-time OCT imaging showed highly motile protozoa in the voids. MABs with protozoa and a high bulk COD (270 mg/L) only had 4% void fraction. DNA sequencing revealed a high relative abundance of amoeba in both high and low-COD MABs. Flagellates were only abundant in the low-COD MAB. Modeling also suggested a relationship between substrate concentrations, diffusion mode (co- or counter-diffusion), and biofilm void fraction. Results suggest that amoeba proliferate in the biofilm interior, especially in the aerobic zones. Voids form once COD limitation at the base of MABs allows predation rates to exceed microbial growth rates. Once formed, the voids provide a niche for motile protozoa, which expand the voids into a large, continuous gap. This increases the potential for biofilm sloughing, and may have detrimental effects on slow-growing, aerobic microorganisms such as nitrifying bacteria.

Gravity-driven membrane filtration for water and wastewater treatment: A review

Gravity-driven membrane (GDM) filtration has been investigated for almost 10 years. The technology is characterized not only by relatively lower transmembrane pressures which can be achieved by gravity (extremely low energy consumption), but also by the phenomenon of flux stabilization: A biofilm is allowed to form on the membrane and a stabilization of flux occurs which is related to biological processes within the biofilm layer on the membrane. This enables stable operation during a year or longer without any cleaning or flushing. Initially, the technology was developed mainly for household drinking water treatment, but in the meantime, the research and application has expanded to the treatment of greywater, rainwater, and wastewater as well as the pretreatment of seawater for desalination. This review covers the field from the rather fundamental research on biofilm morphology and microbial community analysis to the impact of feedwater composition, process parameters and organic removal performance. Not only household applications, but also for community-scale treatment and full-scale applications are discussed. In addition, the application potential is highlighted in comparison to conventional ultrafiltration. Finally, an overall assessment is illustrated and the research and development needs are identified.

Determination of mechanical properties of biofilms by modelling the deformation measured using optical coherence tomography

The advantage of using non-invasive imaging such as optical coherence tomography (OCT) to asses material properties from deformed biofilm geometries can be compromised by the assumptions made on fluid forces acting on the biofilm. This study developed a method for the determination of elastic properties of biofilms by modelling the biofilm deformation recorded by OCT imaging with poroelastic fluid-structure interaction computations. Two-dimensional biofilm geometries were extracted from OCT scans of non-deformed and deformed structures as a result of hydrodynamic loading. The biofilm geometries were implemented in a model coupling fluid dynamics with elastic solid mechanics and Darcy flow in the biofilm. The simulation results were compared with real deformed geometries and a fitting procedure allowed estimation of the Young's modulus in given flow conditions. The present method considerably improves the estimation of elastic moduli of biofilms grown in mini-fluidic rectangular channels. This superior prediction is based on the relaxation of several simplifying assumptions made in past studies: shear stress is not anymore taken constant over the biofilm surface, total stress including also pressure is accounted for, any biofilm shape can be used in the determinations, and non-linear behavior of mechanical properties can be estimated. Biofilm elastic moduli between 70 and 700 Pa were obtained and biofilm hardening at large applied stress due to increasing flow velocity was quantified. The work performed here opens the way for in-situ determination of other mechanical properties (e.g., viscoelastic properties, relaxation times, plastic yields) and provides data for modelling biofilm deformation and detachment with eventual applications in biofilm control and removal strategies.

Physical structure determines compression of membrane biofilms during Gravity Driven Membrane (GDM) ultrafiltration

Increasing transmembrane pressure (TMP) can compress and increase the hydraulic resistance of membrane biofilms. The purpose of the present study is to evaluate how compression of membrane biofilms occurs and how structural rearrangement can affect hydraulic resistance. Biofilms with heterogeneous and homogeneous physical structures were grown in membrane fouling simulators (MFS) in dead-end mode for 20 days with either (i) a nutrient enriched condition with a nutrient ratio of 100:30:10 (C: N: P), (ii) a phosphorus limitation (C: N: P ratio: 100:30:0), or (iii) river water (C: N: P ratio: ca. 100:10:1). The structural and hydraulic response of membrane biofilms to (a) changes in transmembrane pressures (0.06-0.1-0.5-0.1-0.06 bar) and (b) changes in permeate flux (10-15-20-15-10 L/m²/h) were investigated. Optical coherence tomography (OCT) was used to monitor biofilm structural response, and OCT images were processed to quantify changes in the mean biofilm thickness and relative roughness. Nutrient enriched and river water biofilms had heterogeneous physical structures with greater surface roughness (Ra' > 0.2) than homogeneous P limiting biofilms (Ra' < 0.2). Compression of biofilms with rough heterogeneous structures (Ra' > 0.2) was irreversible, indicated by irreversible decrease in surface roughness, partial relaxation in mean biofilm thickness and irreversible increase in hydraulic resistance. Compression of homogeneous biofilm (Ra' < 0.2) was on the other hand reversible, indicated by full relaxation of the biofilms structure and restoration of initial hydraulic resistance. Hydraulic response (i.e., change in the specific biofilm resistance) did not correspond with the change in physical structure of heterogeneous biofilms. The presented study provides a fundamental understanding of how biofilm physical structure can affect the biofilm's response to a change in TMP, with practical relevance for the operation of GDM filtration systems.