1. Introduction
Most of the biotechnological processes are based on planktonic cells in suspension in the cultivation medium [
1]. Bioreactor operations are often limited to batch and fed-batch processes, although continuous processing would be more cost-efficient due to reduced downtime for the reactor cleaning, preparation and cell growth [
1]. Cell retention and a long-term cell viability represent the main challenges in a continuous reactor [
1].
Natural cell immobilization through biofilm formation presents an interesting alternative technique to design new continuous bioprocesses. In nature, biofilms are the predominant lifestyle of bacteria. A biofilm is a multicellular community of one or several bacterial species that is protected through a self-produced polymer matrix. Thereby, biofilms possess an enhanced tolerance to toxic substrates or products compared to the cells in planktonic state and thus, remain viable under unfavorable conditions as well as are able to regenerate themselves [
2,
3]. Due to the high biomass density in biofilms and their stability, biofilm reactors have a high potential for long-term fermentation processes [
4,
5]. However, the biofilm community is highly heterogeneous due to cell differentiation as a result of adaption to nutrients and oxygen gradients inside the biofilm. This heterogeneity makes it challenging to control the growth of the biofilm in the bioreactor.
Many microorganisms are able to grow naturally on diverse surfaces [
1,
3]. In the medical sector, harmful biofilms are a heavy burden since they provoke severe infections and have detrimental effects on human health [
5]. In industrial installation, biofilms can be responsible for biofouling and contaminations and thus present high hygienic risks [
6,
7]. Yet, many industrial applications exist that are taking advantage of biofilms by using them as workhorses. These beneficial biofilms are for example used in the waste-water treatment, bioremediation or the production of bioenergy [
8,
9,
10].
Bacillus spp. are well known for their ability to produce different families of biosurfactant lipopeptides with high application potential such as surfactins, fengycins and iturins [
11]. Previous works have shown that cell immobilization in biofilm bioreactors is particularly favorable for the production of the above-mentioned compounds and allows the design of bioprocesses avoiding excessive foam formation [
12,
13,
14,
15,
16], although biofilm development is a highly dynamic process with instabilities depending on the environmental conditions, such as the release of cells back into the liquid phase upon biofilm disruption. Biofilm development is difficult to assess during the cultivation due to restricted access to the support where the biofilm is growing. It is thus important to develop new measurement and control strategies for monitoring biofilm development and for designing robust processes.
In environmental biotechnology, mathematical modelling of biofilms is used to plan, design, optimize and evaluate processes in wastewater treatment plants [
17]. The implementation of biofilm models permits to calculate the development over time of microbial species and substrates [
18] and to get insights into relevant parameters that control the performance of the biofilm process [
19]. It is important to select only the most relevant parameter to describe the physiological state of the organism and the behavior of the system to reduce the complexity of the model [
20]. These models are developed through the set-up of mass balance equations for the relevant components involved in the bioprocess and the description of the corresponding kinetics expressions [
20]. The components can generally be divided into two categories: the microorganisms and the consumed or produced materials of the microorganisms [
19]. However, the mathematical modeling of biofilm reactors is not always straight forward due to the complexity of biological reactions involved in substrate conversion and the lack of accurate kinetic parameters for the biofilm development [
21]. The approach of inverse modeling has been shown to be an attractive method for the numerical evaluation of kinetic parameters in biofilm processes. Through the validation of the biofilm model with the measured data, the parameters are determined in the way that the observed process behavior is approximately represented through the model [
21].
In this work, a lab-scale two-compartment microbial system composed of a trickle-bed biofilm bioreactor and a stirred tank reactor was designed for the production of surfactin. Through a continuous operation mode, a strong selective pressure was induced on the cell populations. In the actual bioreactor design, biofilm development is promoted to achieve a high cell density on the packing element to increase the production yield. The planktonic cells, in contrast, are not favored and eliminated through a high dilution rate in order to simplify the downstream process of the secreted product. Experimental data are collected with a strong- and poor-biofilm-forming strain derived from B. subtilis 168 for establishing a growth model in order to get a deeper insight into the populations’ behavior. The model is especially useful for predicting the kinetics of the biofilm development on the packing elements, a parameter difficult to assess during cultivation. Moreover, additional information on the system behavior can be obtained through the processing of the model. This provides important information for further process improvement through strain engineering.
2. Materials and Methods
2.1. Strains
The two
B. subtilis strains used in this study and their corresponding genotype are listed in
Table 1. Both strains were derived from the laboratory strain
B. subtilis 168 (
trpC2,
sfp0,
epsC0).
2.2. Biofilm Growth Visualization on Drip-Flow Reactor Coupons
The two
B. subtilis strains were cultivated in a drip-flow reactor device during 48 h on silicone coupons, exactly as described in [
24]. The biofilm is developing on the surface of the coupons which permits to observe easily different biofilm phenotypes. The biofilm images were taken with a Samsung Dual Pixel 12 MP camera at the end of cultivation.
2.3. Design of the Lab-Scale Trickle-Bed Biofilm Reactor and Culture Conditions
A lab-scale (2 L) trickle-bed biofilm bioreactor has been designed on the basis of previous works carried out on a 20 L bioreactor containing a structured stainless steel packing element [
16,
25]. The experimental set-up of the designed reactor is presented in
Figure 1.
For the lab-scale trickle-bed biofilm bioreactor, the system was separated into two main reactors: one for medium mixing and another that contained a tower of five structured metal packing elements for biofilm formation. The packing elements are composed of assembled corrugated gauze stainless steel sheets, a hydrophobic material with good wettability capacities (Laboratory packings, 83 × 55 mm, Sulzer Chemtech, Winterthur, Switzerland). Moreover, the metal structured packing provides an increased gas/liquid mass transfer.
The medium was recirculated continuously between these two devices with a flow rate of 85 mL min−1. The medium was mixed at 300 rpm in the reactor. The mixing reactor was a conventional 2 L bioreactor (BIOSTAT B Plus, Sartorius Stedim, Schaerbeek, Belgium) whereas the reactor containing the packing elements was composed of a previous 2 L chemical reactor with a double jacket for temperature regulation (Reactor-Ready, Radleys, Shire Hill, Saffron Walden (Essex), UK). Since this type of reactor does not possess a condenser which is crucial to avoid filter clogging and pressure problems due to medium evaporation, the gas outlet was refrigerated by an additional cooling system to reduce evaporation in the packing reactor. The temperature of both reactors was regulated to 37 °C. For security, the gas outlet was connected to a reservoir bottle with filters in case of too strong evaporation to collect the condensate. During the cultivation, there is no aeration in the mixing reactor. Air (1 L min−1) is injected only on the downside of the packing reactor to prevent foam formation. The medium is injected on the upper side of the packing reactor and then flows down by gravity on the packing elements. Oxygen mass transfer is promoted through the counter-current flow of the injected air and the liquid. For the continuous process mode, an entry to and exit from the mixing reactor was added. Samples were taken from the mixing reactor by means of a sterile syringe.
For the reactor inoculation, a series of pre-cultures was prepared. First, 2 mL of lysogeny broth (LB) medium (10 g L−1 tryptone, 5 g L−1 yeast extract, 10 g L−1 NaCl) was inoculated with a colony. The first pre-culture was incubated for about 6 h at 37 °C and 160 rpm. Then, a second pre-culture was prepared by a 10 times dilution of pre-culture I in LB medium. The second pre-culture was incubated overnight at 37 °C and 160 rpm and then 10 times diluted with Landy MOPS medium (20 g L−1 glucose, 5 g L−1 glutamic acid, 1 g L−1 yeast extract, 0.5 g L−1 MgSO4, 1 g L−1 K2HPO4, 0.5 g L−1 KCl, 1.6 mg L−1 CuSO4, 1.2 mg L−1 MnSO4, 0.4 mg L−1 FeSO4, 21 g L−1 MOPS, 1.6 mg L−1 tryptophan) to prepare the main pre-culture. The main pre-culture was grown to an OD600 nm between 2 and 3 and then used to inoculate the reactor (1 L working volume) with an OD600 nm of 0.2 (corresponds to ~0.08 g L−1 cell dry weight). Before inoculation, the cells were washed once in a 0.9% NaCl solution to synchronize the cells and eliminate the produced primary and secondary metabolites. For this purpose, the cell culture was centrifuged (10 min at 2700× g) and the supernatant was discarded. The remaining cell pellet was resuspended in a 0.9% NaCl solution and then used to inoculate the reactor. The reactor contained Landy medium without MOPS buffer. The reactor pH regulation was executed using 1 M H3PO4 as acid and 3 M NaOH as base. The pH in the reactor was set at 7.0. To the reactor medium 50 µL L−1 of a silicone-free organic antifoaming agent (TEGO® Antifoam KS911, Evonik, Essen, Germany) was added. The culture was started with a batch fermentation during 16 h to increase the cell number in the reactor and to promote cell adhesion and biofilm development on the support. Then, the continuous phase was launched during ~28 h with a dilution rate of D = 0.5 h−1 which corresponds to a feeding rate of 500 mL h−1. Two replicates of the biofilm cultivation experiments were performed per strain.
2.4. Determination of the Mean Residence Time in the Packing Tower
For the mean residence time determination in the packing tower, tracer particles (1 µm) were injected on the top of the packing tower with a flow rate of 85 mL min
−1 and collected at the packing tower exit at time intervals of 5 s. The collected particles were counted by flow cytometry (BF Accuri
TM C6, BD Biosciences, Erembodegem-Dorp, Belgium). The mean residence time was then calculated with the measured tracer concentration over the time by Equation (1):
where c represents the measured tracer concentration at time point t.
2.5. Biomass Dry Weight Determination
The cell culture samples were centrifuged (10 min at 2400×
g) and the supernatant was collected to determine the surfactin concentration as described in
Section 2.7. The remaining cell pellets were washed by resuspending them in distilled water followed by centrifugation. After the centrifugation, the supernatant was discarded, and the remaining cell pellet was re-dissolved in water and filled into a pre-weighted aluminum cup. The biomass sample in the aluminum cup has been dried in the oven at 105 °C and weighted after 48 h to determine the corresponding cell dry weight.
After the cultivation, the packing elements were left for two hours in the reactor to let drain the residual liquid before they were weighted with the wet biofilm. The corresponding biofilm dry weight was then calculated using a previously determined biofilm dry weight percentage (see
Data S1).
2.6. Glucose Analysis
Glucose concentration was analyzed in the supernatant using high-performance liquid chromatography (HPLC). A Waters Acquity UPLC® H-Class System (Waters, Zellik, Belgium) with an ion-exchange Aminex HPX-87H column 7.8 × 300 mm (Bio-Rad Laboratories N.V., Temse, Belgium) heated up to 50 °C was used for analysis. A metabolite analysis was carried out with an isocratic flow rate of 0.6 mL min−1 for 25 min. The mobile phase was composed of water containing 5 mM H2SO4. Elution profiles were monitored through a Waters Acquity® Refractive Index Detector (RID) (Waters, Zellik, Belgium). A glucose standard solution (Sigma-Aldrich, Overijse, Belgium) was used to determine the retention time and to establish a calibration curve.
2.7. Surfactin Analysis
The supernatants from the centrifuged cell culture samples were filtered (0.2 µm) and the surfactin concentration was determined by reversed-phase HPLC (Agilent 1100 Series HPLC Value System, Agilent Technologies, Diegem, Belgium) with an Eclipse XDB C−18 column (3.5 µm, 2.1 × 150 mm) (Agilent Technologies, Diegem, Belgium). The HPLC analysis method was based on an isocratic elution profile with a mobile phase composition of 80% acetonitrile and 20% water containing 0.1% trifluoroacetic acid (TFA). The flow rate was set at 0.4 mL min−1 with an analysis time of 22 min per sample. The surfactin molecules were detected by UV at 214 nm. Purified surfactin samples (>98%) (Lipofabrik, Villeneuve d’Ascq, France) were injected to identify the retention time of the surfactin molecules and to determine a calibration curve.
2.8. Total Organic Carbon Analysis and Establishment of the Mass Balance
Total organic carbon (TOC) measurements of the culture medium were performed in order to estimate the TOC consumption of the cells. The planktonic cells were separated from the bulk medium by centrifugation. Subsequently, the TOC content of the culture medium was measured using a Lotix Combustion TOC Analyzer (TELEDYNE TEKMAR, Mason, OH, USA). The diluted culture medium samples were injected into the combustion tube where the samples were completely oxidized to CO2 through catalytic combustion at 720 °C. Subsequently, the produced CO2 was detected by flow-through non-dispersive infrared spectroscopy. The instrument was calibrated with a standard solution of potassium hydrogen phthalate for a calibration range of 0 to 20 ppm.
A TOC mass balance for the batch and continuous cultivation phase has been established to estimate the TOC consumed by the cells present in the system. The TOC consumption in the batch phase was determined by Equation (2):
and for the continuous phase, with Equation (3):
where TOC represents the amount of total organic carbon in g at a certain time point t in h, D is the dilution rate in h
−1 and Δt the difference between time point t
0 and t
1 in h. For the mass balance of the continuous phase, it was assumed that the TOC consumption rate remains constant during the measured time interval.
2.9. Biofilm Reactor Compartment Model
The trickle-bed bioreactor system can be simplified into two main compartments representing the main places of residence in the system for sessile and planktonic cells. The first compartment comprises the sessile cells which form a biofilm on the packing tower where cells have unlimited access to dissolved oxygen. In the second compartment, the planktonic cells are growing under limited dissolved oxygen conditions. The overall growth conditions for the cells are better in the packing tower since there is more dissolved oxygen available as well as enough nutrients since the medium is continuously recirculated. In the present system, the contact between the injected gas and the adhered bacteria on the packing elements is strongly enhanced which favors additionally an interfacial oxygen transfer through a direct bacteria-air contact contributing to an increased total oxygen transfer [
26].
The objective of this process is to increase the adhered biomass on the packing tower and reduce or eliminate the presence of planktonic cells in the stirred tank reactor through a high dilution rate (D > µ). Increased cell density on the packing elements means increased production yield. The elimination of planktonic cells would strongly facilitate the downstream process since the secreted product in the bulk medium could be easily recovered. This means that the aim for this system is not to reach a steady state as in a normal chemostat reactor (µ = D), but a steady state with a planktonic cells number close to zero and a continuously and stable growing biofilm.
In this work, the development of the two populations (sessile and planktonic cells) was investigated in order to acquire more information about their behavior for further process optimization. By means of experimental data, a simple ordinary differential equations (ODEs) model was established on the basis of bacterial growth equations. A schematic description of the model is presented in
Figure 2.
2.10. Determination of the Volumetric Oxygen Mass-Transfer Coefficient KLa in the Stirred Tank Reactor by Dynamic Gassing-In/Gassing-Out Method
The oxygen transfer rate from a gas to a liquid phase is given by Equation (4):
where K
La is the volumetric oxygen mass-transfer coefficient, C
sat the oxygen concentration at saturation in the liquid medium in equilibrium to the gas phase and C
L the dissolved oxygen concentration in the liquid medium [
27]. In a stirred tank reactor where the liquid phase is well mixed, the accumulation of oxygen in the liquid phase can be described through Equation (5):
where OTR is the oxygen transfer rate from the gas to the liquid and the OUR represents the oxygen uptake rate by the biomass [
27]. Since the volumetric oxygen mass-transfer coefficient K
La has been measured in the absence of biomass, OUR = 0, Equation (5) can be simplified to Equation (6):
And thus can be transformed into Equation (7):
The dynamical absorption method [
27] was applied in order to determine the K
La value. This method consists of the elimination of oxygen in the liquid phase to obtain an oxygen concentration close to zero through the injection of nitrogen. This permits to simplify further Equation (7) with t
1 = 0 and C
L1 = 0%. Then, the liquid is again put into contact with air and the increase of the dissolved oxygen concentration is measured over the time. The K
La can then be deviated through the slope of the
vs. time plot. Measurements were performed in triplicates.
2.11. Mathematical Development of a Growth Model to Describe the Microbial Population Dynamics
The following assumptions are made for the model construction: (i) no oxygen limitation in the biofilm compartment with the packing tower, (ii) the oxygen concentration in the planktonic cell compartment is limited, (iii) the dilution rate is affecting directly the planktonic cell compartment but not the biofilm compartment. The used model parameters are listed in
Table 2.
2.11.1. Batch Fermentation
The growth rate of bacteria can be described through the well known Monod equation of growth represented by Equation (8):
The total biomass development in the trickle-bed biofilm reactor can be divided into the growth of sessile and planktonic cells. For the planktonic cells, two limiting factors have to be taken into account: the substrate and dissolved oxygen concentration. If oxygen is a limiting factor, the specific growth rate varies with the dissolved oxygen concentration according to the Monod equation like for any other substrate limitation. In our case, oxygen and substrate are complementary substrates and thus, the product rule is applied [
28]. The growth speed for the planktonic cells is thus given through Equation (9):
For the model, Equation (9) was adapted according to the approach used by Roels [
29], as shown in Equation (10):
Here, the growth speed of the planktonic cells is assumed to be influenced by the more pronounced limiting factor which means the minimum value of the term representing either the substrate limitation or the limited dissolved oxygen availability.
The biomass development for the planktonic cells can be described by the differential Equation (11):
where k
a represents the switching rate from the planktonic state to the sessile state of the cells (adsorption) and k
d the releasing rate of the sessile cells to the planktonic state (detachment). Thus, the term k
aX
p correspond to the number of planktonic cells that adhere to the support whereas k
dX
b describes the sessile cells detaching from the support.
The growth speed for the sessile cells can be described through Equation (12) by taking into account the substrate limitation due to the randomly distributed medium on the packing elements:
In this case, dissolved oxygen limitations are not considered for the sessile cells in the model. It can be assumed that the aeration is very efficient in the packing tower and the biofilm thickness is sufficiently low to neglect oxygen gradients.
The development of the biofilm on the packing elements can be described by Equation (13):
The terms kaXp and kdXb represent the corresponding biomass that is adhering or detaching as described above.
The substrate consumption of the sessile and planktonic cells is given by Equation (14):
The availability of dissolved oxygen can be described by Equation (15):
where the terms
and
represent the oxygen uptake rate of the biofilm and planktonic cells, respectively. For the model, it was assumed that the dissolved oxygen concentrations were equivalent for both compartments due to the continuous recirculation of the medium between the stirred tank reactor and the packing tower.
2.11.2. Continuous Fermentation
For the continuous fermentation, the dilution rate affects only the planktonic phase. The supply and removal of dissolved oxygen through the alimentation and elimination is neglected. This means that Equation (11) describing the development of the planktonic biomass is extended with the term in bold in Equation (16):
And Equation (14) describing the substrate consumption is extended as shown by Equation (17):
The ODEs were coded and solved with Python 3.7 via the Anaconda–Spyder interface using the odeint function (see
Data S2 for the code).
4. Discussion
The objective of this work was to develop a model able to describe the growth dynamics of the biofilm and planktonic population present in the designed trickle-bed biofilm reactor in order to understand better the behavior of the system for further process intensification. In particular, biofilm development on the packing elements gives important information about the process, but is difficult to monitor during cultivation. The growth model was used in order to confirm the two hypotheses that the significant difference in biofilm development of BBG111 and RL5260 is linked to the production or not of EPS, and that the high dilution rate washes out the non-adherent or detaching cells in the designed system.
The experimental data are in good accordance with those obtained with the developed growth model by using a combination of the first hypothesis (different ka/kd ratio during batch cultivation) and the second one (ka = 0 during continuous cultivation) for both strains. This was confirmed by a Chi-square goodness to fit test with a confidence level of α = 0.05. The two hypotheses made initially for the present system have thus been validated. The model also confirmed the biofilm development dynamics determined via experimental TOC measurements and the established TOC mass balance.
The increased k
a/k
d ratio for RL5260 during the batch cultivation was linked to the capacity of EPS secretion which has shown to improve the colonization capacity and reduce cell detachment. The presence of EPS permitted RL5260 to build up a functional biofilm structure and to protect the adhered cells from external influences. Once adhered, the cells produced EPS and proliferated on the packing elements to construct their own environment. Several works on
B. subtilis biofilm formation have shown that EPS production facilitates cell spreading and promotes the colonization of a solid support [
24,
39,
40]. Since BBG111 is a poor biofilm former and does not produce EPS, the biofilm formation capacities were reduced (lower k
a/k
d ratio) and cell detachment occurred more frequently after the cell adhesion step than in the case of RL5260. Moreover, the additional high dilution rate carried out a strong selective pressure on the planktonic cells and limited the re-adherence during the continuous cultivation due to the washing out of the planktonic cells (k
a = 0).
Globally, BBG111 and RL2560 produced comparable amounts of planktonic cells. RL5260 produced more planktonic cells during the batch phase. However, when the continuous cultivation phase was launched, the planktonic cells were mostly washed out for RL5260 whereas the number of planktonic cells of BBG111 increased during the continuous phase. This was probably a result of the limited adhesion capacities of this strain due to the absence of EPS. The maximum specific growth rates of 0.39 h
−1 and 0.38 h
−1 for BBG111 and RL5260 were comparable in the two-compartment system. They were close to the values of 0.35 h
−1 and 0.38 h
−1 determined by Guez et al. [
31] and Martínez et al. [
41] as growth rates for
B. subtilis in glucose-limited fed-batch cultures. However, the EPS
+ strain RL5260 produced about 1.6 more adhered biomass than BBG111 (EPS
-) which resulted in an important difference regarding the biofilm versus planktonic cell ratio. This ratio was nearly two times higher for RL5260.
Although EPS production is advantageous for cell adhesion and leads to enhanced biofilm formation, it is metabolically expensive [
42]. Thus, EPS production may reduce the cell growth and affects negatively the surfactin production. Nevertheless, the results have shown that the mean surfactin productivity of the strain RL5260 with increased biofilm formation capacity through EPS production was about 37% improved during the continuous phase compared to BBG111. This indicates clearly the improved performance of the EPS
+ strain in this system compared to the EPS
- strain.
Surprisingly, both strains showed a similar glucose consumption profile when the concentration was measured in the bulk medium. For the same amount of consumed glucose, RL5260 produced significantly more adhered biomass as well as higher amounts of surfactin than BBG111. This indicated that both strains had a completely different cell physiology in the system due to the differences in EPS production.
Regarding the performed TOC measurements, the TOC consumption profile for BBG111 was significantly increased compared to RL5260 during the continuous phase. This was most likely linked to the different biofilm development capacities due to the production of EPS or not of RL5260 and BBG111. Hence, RL5260 and BBG111 used the available carbon source in the medium in a different way. Given that BBG111 is not able to synthesize a biofilm matrix, the adhered biomass consisted mainly of cells whereas the adhered biomass of RL5260 contained a mixture of cells and biofilm matrix. Biofilm composition measurements of RL5260 that were performed in our laboratory using biofilms developed on drip-flow reactor coupons revealed a relative EPS amount of 81% and a cell content of 19%. Both strains show comparable glucose-to-biomass conversion yields for the cellular production in suspended cell cultures (Y
X/S,cells of 0.22 g g
−1 for BBG111 and Y
x/s,
cells = 0.26 g g
−1 for RL5260). The yields were similar or close to the yield of 0.22 g g
−1 previously reported by Guez et al. [
31] for
B. subtilis ATCC6633 grown in Landy medium in shaking flasks. Assuming a substrate-to-EPS conversion yield that is significantly higher than the conversion yield for cellular production, e.g., Y
x/s,EPS ~ 0.57 g g
−1 as obtained by Huang et al. for the production of poly-γ-glutamic acid (PGA), a major extracellular compound of
B. subtilis CGMCC1250 [
43], RL5260 used in total lesser amounts of carbon sources than BBG111 for the biofilm development. This assumption is further confirmed through the determined yield of the produced biofilm per consumed TOC Y
biofilm/TOC for both strains. RL5260 reached a yield of Y
biofilm/TOC = 0.61 g g
−1 whereas BBG111 reached only 0.25 g g
−1. This shows the lower energy consumption of RL5260 for the biomass production due to the increased biosynthesis of EPS instead of cells.
The reduced energy consumption of RL5260 for the biofilm development resulted in a more efficient surfactin production. It can be considered that surfactin was mainly produced by the cells present in the biofilm since a sufficient aeration is necessary for the production which was not guaranteed for the planktonic cells in the stirred tank reactor. Consequently, the specific surfactin production was significantly increased for RL5260. Hence, RL5260 reached a mean specific surfactin production of 90 mg L−1 h−1 per g of adhered cell dry weight whereas BBG111 produced only 20 mg L−1 h−1 per g of adhered cell dry weight.
In conclusion, the two-compartment biofilm reactor designed in this study has shown to be suitable for continuous surfactin production. The EPS+ strain exhibited significantly improved performances in terms of cell adhesion and surfactin production in this system by comparison with the EPS− strain. The surfactin yield and population stability inside the reactor could be further improved by engineering the biofilm formation capacity of the cells. For a good process performance, a trade-off between enhanced cell adhesion and increased productivity has to be chosen. EPS production could be modulated in favor of surfactin production by guaranteeing a sufficient cell adhesion through a controlled EPS production while increasing the numbers of potential cell factories. Moreover, cell morphology engineering could improve cell adhesion and further reduce cell detachment.